2012-01-05 16 views
2

私はPythonの初心者です。私は私がmatlabで持っているいくつかのプログラムを切り替えるしようとしています。 私は、入力がすべての行列である堅いode方程式を解く必要があります。私が使用するmatlabでPythonで堅いodeを解決するには?

[ttT,uT] = ode23s('SST',t,fT); 
+3

はあなたにも数学の記法でそれを記述することができますか?誰もがmatlabに精通しているわけではありません。 –

+0

これは有用であり得る:http://stackoverflow.com/questions/2088473/integrate-stiff-odes-with-python –

+0

関数ode23sは、この例ではMatlabの に硬い方程式ソルバーで、SSTは、AUX関数であります状態空間関係を確立する.z = A * z + B、AおよびBは行列である。 Aは正方形でBは列です。 tは時間的増分を有する配列であり、fTは初期条件を有する配列である。 – marco

答えて

2

Matlabで行うことは、PythonのNumPyモジュールで行うことができます。それはhereであることがわかります。

また、同様に関連モジュールSciPyが役立つかもしれません。

PyDSToolもあなたに関連性があるかもしれません。これはRadauソルバーのラッパーです。

その後、プロットするためmatplotlibを試してみたいことがあります。それはMatlabのプロットのようなものです。

次のリンクは、あまりにも、役立つかもしれない:

+0

ありがとうございました。私はすでにnumpyとscipyをインストールしました。 Googleでは、mrgatesページが見つかりました。しかし、彼の "モジュール"のインストール方法はわかりません。私はゲイツのアプローチがもっと簡単になると思う。 私はEPD加入者ではありません。 – marco

+0

最後のリンク、申し訳ありませんが、あまりにもあいまいです。私はどちらが有用かどうか分からなかった!私は「ただ」簡単なことをしたいが、私はPythonで理解している、そうではない! – marco

+0

Gatesのモジュールをインストールするには、以下の手順に従います。1)アーカイブを解凍します。 2)C:\ ode-0.1a4.tarに解凍したとします。このフォルダの中に 'setup.py'というファイルがあることを確認してください。 3)このフォルダでコマンドラインを開き、次のコマンドを実行します: 'python setup.py install'。これはあなたのマシンにモジュールをインストールするはずです。ゲイツ氏が自分のページで説明しているように、このモジュールを読み込んで使用することができます。 –

関連する問題