問題タブ [odeint]

では、curve_fit を使用して、dy/dt = a*y^2 + b のように、微分方程式の定数 a と b を最適化するコードを作成するにはどうすればよいでしょうか。ODE を解くために odeint を使用してから、a と b を最適化するために curve_fit を使用します。この状況についてご意見をお聞かせいただければ幸いです。

ブースト 1.53 アウトと OdeIntがあります。これには、さまざまな数値 ode ソルバーが含まれています。それでも、そこにlsodeアナログがあるのだろうか - あなたが分析ソリューションを保証した場合のために?

odeint C++ ライブラリを使用して、離散化された偏微分方程式から結合された一連の ODE の数値積分を実装しました。うまく機能し、非常に高速ですが、1 つの問題があります。

ODE の私のシステムには、いわゆる吸収境界条件があります。N double のベクトル (人口密度) である状態変数 n の時間導関数は、システム関数で計算されますが、それが起こる前 (または後)時間積分) 設定したい:

ただし、システム関数の状態変数がconstとして宣言されているため、もちろんこれは機能しません。これは、ライブラリをいじる以外に変更できないように見えます...これを回避する方法はありますか?

dndt[N]とを 0に設定dndt[N-1]することは解決策のように見えるかもしれませんが、境界条件を吸収するという概念に反するため、実際には役に立ちません (n[N]そして、任意の時点での の値ではなく、 でn[N-1]の値を常に持つことになります)。間に合うように)、だから私は本当に n を変更したいと思います。t=0n[N-2]

助けてくれてありがとう！

よろしく、

マイケル

odeintNumpyでのシミュレーション中に中間変数を保存する最も簡単な方法は何ですか?

例えば：

xシミュレーション中に保存されたデータを保存する最も簡単な方法は何ですか? tに渡される引数で指定されたポイントが理想的odeintです。

odeint で 2 点境界問題を解こうとしています。私の方程式は次の形をしています

の境界条件がある場合はかなり些細なことですy(x_1) = y_1 , y'(x_2) = y_2が、境界条件がの場合y(x_1) = y_1 , y(x_2) = y_2は失われます。このような問題を odeint や他の科学ライブラリで処理する方法を知っている人はいますか?

ODE ソルバーを使用しboost::numeric::odeintていますが、それらについて質問があります。説明書にはbulirsch_stoer、ステップサイズと順序制御を備えたステッパーであると書かれています。また、順序はユーザーが変更できる変数です。ヘッダーファイルを見ましたが、わかりませんでしたか? 誰にも解決策はありますか？