huzhaolin 发表于 2006-12-20 16:45

各位大哥,如何解此微分方程?

x''+a*y''+c*x'+d*x+e=0;
y''+f*y'+g*y+h=0;
如果用ode45 的话,方程应该这么表示可以吗?
令u1=x,u2=x',u3=y,u4=y';
du1=u2;
du2=-a*du4-c*u2-d*u1-e;
du3=u4;
du4=-f*u4-g*u3-h;

shenyongjun 发表于 2006-12-20 21:05

很简单的问题啊
du1=u2;
du2=-a*u4-c*u2-d*u1-e;
du3=u4;
du4=-f*u4-g*u3-h;

huzhaolin 发表于 2006-12-20 21:15

du2=-a*u4-c*u2-d*u1-e;
好象不是这样的,因为x''=-a*y''-c*x'-d*x;
右边是y'',不是y';
这真是我想问的问题.
如果令u1=x,u2=x',u3=x'',u4=y,u5=y',u6=y''是可以用ode解,但是方程成六元了,解起来慢.不知大狭们怎么解决啊!
页: [1]
查看完整版本: 各位大哥,如何解此微分方程?