pablo_001 发表于 2006-3-17 13:30

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

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

zhchtju 发表于 2006-3-17 15:31

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

pablo_001 发表于 2006-3-18 11:12

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

AaronSpark 发表于 2006-3-18 11:19

回复：（pablo_001）求助龙格库塔解两元一阶微分方程组...

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

linqus 发表于 2006-3-18 11:22

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

应该没有什么问题的

pablo_001 发表于 2006-3-19 13:19

程序

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

xdot=;

t_final=100;

x0=;

A=-1.2521;

B=9.1980;

C=-0.9801;

D=-1.2425;

E=16.0199;

F=0.0828;

w=30;

=ode45('kaihuan',,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 发表于 2006-3-19 13:21

请问如何验证是刚性方程?谢谢

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

happy 发表于 2006-3-19 14:28

回复：（pablo_001）程序

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

xdot=;

t_final=100;

x0=;

A=-1.2521;

B=9.1980;

C=-0.9801;

D=-1.2425;

E=16.0199;

F=0.0828;

w=30;

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



plot(t,x),

figure;%打开新图形窗口

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

pablo_001 发表于 2006-3-21 12:50

谢谢

谢谢指教!
页: [1]
查看完整版本: 求助龙格库塔解两元一阶微分方程组