2
y'''' + (a1 + a2*v(t)^2)*y'' + (a3 + a4*v(t)^2)*y = 0 
y(0) = 2, Dy(0) = 0, D2y(0) = 0, D3y(0) = 0

上記の4次常微分方程式をMatlabで解こうとしましdsolveたが、得られた結果は非常に大きかったです。私が今理解しているように、dsolveここでは機能しませんode45。この方程式を解くために使用する必要があります。

ode45でMatlabのドキュメントをフォローしようとしましたが、完全には理解できませんでした。さらに、私の場合、係数も時間に依存します。

この方程式を理解できる形式に変換するにはどうすればよいode45ですか?

4

1 に答える 1

2

これは4次常微分方程式です。直接解くことはできません。これを一次常微分方程式のシステムとして書き直す必要があります。

y_3' = -(a1+a2*v(t))*y_2 - (a3+a4*v(t)^2)*y
y_2' = y_3
y_1' = y_2
y' = y_1

y(0) = 2
y_1(0) = 0
y_2(0) = 0
y_3(0) = 0

これで、ODE45にフィードできます...

于 2013-03-07T09:09:52.553 に答える