ode45(Matlab)のしくみ

解決済

質問者：kiyotamakiyota
質問日時：2012/11/18 21:31
回答数：1件

MATLABの関数ode45は
　dy/dt=F(t,y)
の型の微分方程式の解yを得るのにもちいりますが，基本的なプログラムの表記は
[T,Y]=ode45('F',tspan,yo)
と書きますよね．ここで得られるT,Yは
T:to　t1 t2 t3.....tend
Y:yo y1 y2 y3.....yend
であるわけですが，実際このYはどのようにして得られた値なのでしょうか？

ｙ２なら例えば
　　y2＝dy/dt*(t2-t1)+y1
などのようにして得られているのでしょうか？

通報する

この質問への回答は締め切られました。

質問の本文を隠す

回答 (1件)

No.1ベストアンサー

回答者： masudaya
回答日時：2012/11/20 15:51

ヘルプを見ると

ode45 は、Dormand-Princean の陽的 Runge-Kutta (4,5) 公式に基づきます。これは、単段階ソルバです。すなわち、y(tn) の計算では、直前のy(tn-1) における解のみ必要とします。一般に、 ode45 は、ほとんどの問題に対して最初の試みとして利用すべき最良の関数です。 [3]

[3] Dormand, J. R. and P. J. Prince, “A family of embedded Runge-Kutta formulae,” J. Comp. Appl. Math., Vol. 6, 1980, pp 19-26.

と記載があります.この文献を当たれば分かるのでは無いでしょうか.
添付に関連する情報があります.

http://odeity.googlecode.com/hg/doc/papers/an-ef …
http://www.sciencedirect.com/science/article/pii …