zlwr454 发表于 2009-11-10 23:17

请问为什么这个常微分方程组运行老是BUSY.....

function ydot=vdpo1(t,y)
Q11=2500;
Q21=1800;
Q41=700;
Qe1=500;
mu=2;
ydot=;

然后再运用ODE45求解：tspan=;yo=;
=ode45(@vdpo1,tspan,yo)；
plot(t,y(:,1),t,y(:,2),'--')
请问为什么一个简单的微分方程组运行一直BUSY？谢谢~！

ChaChing 发表于 2009-11-11 10:48

回复 楼主 zlwr454 的帖子

ydot第三项(Q11+Q41-Q21-Qe1)/24=37.5, 但yo(3)=30;
如此会不会不收敛?
页: [1]
查看完整版本: 请问为什么这个常微分方程组运行老是BUSY.....