ode

3熱

1答えて

私はPythonを初めて使用しましたが、C++とMATLABの経験があります。私は現在、dy/dtとdx/dtを含む非線形システムの位相空間内の軌跡をプロットするプログラムを作成しています。しかし、より複雑な関数形式の場合、私はオーバーフローエラーを受け取りました。この問題を回避する方法はありますか？前もって感謝します！これらは私のコードです： fig = plt.figure(figsize

3熱

1答えて

Sympy ODEソルバーの初期条件は無効です

dsolveを使用して初期値の問題を解決したいが、初期条件icsは効果がないようです。この例のdsolveは両方とも同じ結果を示します。私は x f(x) = C1 ℯ C1同じ解を得るいずれの場合も from sympy import * x = symbols('x') f = Function('f') eq = Eq(Derivative(f(x), x), f(x))

0熱

1答えて

GSL ODEソルバへのC++クラスメンバ関数

私は明らかによく知られている問題に苦しんでいます。クラスメンバー関数で定義されたODEシステムがあり、GSLソルバの1つで解決/統合したいと考えています。私の統合ルーチンに続いて class my_model{ ... public: int NEQ = 4; double y[4], dydt[4]; double params[25]; int od

1熱

1答えて

このRK4 ODEソルバでは、内側のラムダが外側のラムダfuctionをどのように使用しますか？ Pythonの

私は数値解析のための私のコードを改善するために、私の検索では、次のコードを見つけました：今 from math import sqrt def RK4(f): return lambda t, y, dt: ( lambda dy1: ( lambda dy2: ( lambda dy3: ( lambda dy4: (dy1 +

0熱

1答えて

C＃のforeachループの出力から配列を作成する方法は？

px [0]が時間とともに変化し、配列が必要なので、foreachの出力（pX [0]など）を配列に入れて計算プロセスを完了し、通常の微分方程式を解く必要がありますそれらすべてを保つ。いくつかのデータを配列に格納する最も良い方法は "forループ"ですが、Microsoft Osloライブラリを使用した場合は、Osloユーザーガイドに従ってこのように記述する必要があります。出力を助けることがで

0熱

1答えて

scipyのsolve_bvpを開始できません

私はscipyのsolve_bvpを使って解決しようとしている8つの境界条件を持つ8つの微分方程式のシステムを持っています。連結軸を除くすべての入力配列の次元：とValueError：それと間違って何かがありますように、私は次のエラーメッセージを取得していますので、それは、見えますしかしあなたは、私が import numpy as np from scipy.integrate import

2熱

1答えて

Pythonの大きな複合配列の指数関数[exp（）]を計算する最も速い方法

私は、scipyのcomplex_odeを使ってODEを統合するコードを開発しています。ここで、integrandにはフーリエ変換と指数演算子があり、値。パフォーマンスを最適化するために、私はこれをプロファイリングし、主なボトルネックはラインに（PyFFTWなどを使用してFFTを最適化した後）であることがわかりました。 val = np.exp(float_value * arr) 私は現在

1熱

1答えて

ode45を使用中にインデックスが行列の次元を超えています

ode45を使用してodeを解決しようとしていますが、表示されています。 function ydot=vdpol(t,y) t=0; y=0; mu=2; ydot=[y(2);mu*(1-y(1)^2)*y(2)-y(1)]; tspan=[0 20]; y0=[2;0]; [t,y]=ode45(vdpol,tspan,y0); size(t) size(y) plot(t,