ode45関数で実行したら、警告が出て途中で終了してしまいます。

Question

0 个投票

以下のコマンドを実行したら、警告が出ます。 [t,y] = ode45(@odefcn,[0 10],[0; 10]);

警告: t=8.956588e+00 で失敗。時間 t で最小の値で許可された (1.421085e-14) より小さいステップサイズに減らさずに積分の許容誤差を満たすことができません。

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

Answer 1

Takafumi 2017-5-11

2 个投票

まずソルバーを替えます。 ode45 から、ode15s に変更する(スティッフソルバーに)

>> [t,y] = ode15s(@odefcn,[0 10],[0; 10])

それで、解決しない場合は、ソルバーの許容誤差の設定を緩めることで、実行する事ができます。

>> opts = odeset('AbsTol',1e-4,'RelTol',1e-4);
>> [t,y] = ode45(@odefcn,[0 10],[0; 10],opts);

ちなみに、ソルバー選択はマニュアルに記載があります。

https://jp.mathworks.com/help/matlab/math/choose-an-ode-solver.html