求助龙格库塔解两元一阶微分方程组

pablo_001

求助：请问哪位兄台用过龙格库塔解两元一阶微分方程组，请提供算法程序，我做了好多次解都是１，谢谢

zhchtju

ode45的算法就是45阶Runge－kutta法
直接调用函数解就好了

pablo_001

<P>可我试过好多次都不能得出解，请详细说明，谢谢</P>

AaronSpark

把程序贴出来大家帮你看看，重新写一个太麻烦了

linqus

（1）确认你的方程不是刚性方程，否则没法用explicit Runge-Kutta法求解
可以用implicit Runge-Kutta，gear等等方法求解。
（2）请找一本matlab的书，将用ode45求解常微分的例子自己作一遍，
再做你自己的方程求解。

应该没有什么问题的

pablo_001

function xdot=kaihuan(t,x,flag,A,B,C,D,E,F,w)

xdot=[A*x(1)+B*x(2)+E*w;C*x(1)+D*x(2)+F*w];

t_final=100;

x0=[0;0];

A=-1.2521;

B=9.1980;

C=-0.9801;

D=-1.2425;

E=16.0199;

F=0.0828;

w=30;

[t,x]=ode45('kaihuan',[0,t_final],x0,[],A,B,C,D,E,F);

plot(t,x),

figure;%打开新图形窗口

plot2(x(:,1),x(:,2));

??? Strings passed to EVAL cannot contain function declarations.



请各位高手提供修改建议，谢谢!

pablo_001

刚性方程如何严验证?谢谢

happy

function xdot=kaihuan(t,x,flag,A,B,C,D,E,F,w)

xdot=[A*x(1)+B*x(2)+E*w;C*x(1)+D*x(2)+F*w];

t_final=100;

x0=[0;0];

A=-1.2521;

B=9.1980;

C=-0.9801;

D=-1.2425;

E=16.0199;

F=0.0828;

w=30;

[t,x]=ode45('kaihuan',[0,t_final],x0,[],A,B,C,D,E,F,w);



plot(t,x),

figure;%打开新图形窗口

plot(x(:,1),x(:,2));

pablo_001

谢谢指教!