mon890 发表于 2006-11-13 19:21

求助,帮忙解这个微分方程组,万分感谢

dx(t)/dt=(1/2)*Gn*(z(t)-nth)*x(t)+(k/tin)*x(t-tt)*cos(y(t)-y(t-tt));
dy(t)/dt=(1/2)*a*Gn*(z(t)-n0)-(k/tin)*(x(t-tt)/x(t))*sin(y(t)-y(t-tt));
dz(t)/dt=(1.3*Ith)/(e*V)-(1/ts)*z(t)-Gn*(z(t)-n0)*x(t)^2;
Gn=8.4e-13;nth=2.0e+24;k=0.015;tin=8e-12;
n0=1.4e+24;k=0.015;a=6;
e=1.6e-19;V=1.0e-16;ts=2e-9;Ith=54e-3;tt=1.2e-9;

xjzuo 发表于 2006-11-13 19:24

回复

时滞方程,建议用dde23求解.

captainmeng 发表于 2006-11-13 20:11

自己边程序啊,用龙格库塔法,很快的
页: [1]
查看完整版本: 求助,帮忙解这个微分方程组,万分感谢