[MATLAB数学相关] 请问为什么这个常微分方程组运行老是BUSY..... - MATLAB爱好者论坛-LabFans.com

zlwr454 · 2009-11-08, 21:31

function ydot=vdpo1(t,y)
Q11=2500;
Q21=1800;
Q41=700;
Qe1=500;
mu=2;
ydot=[y(2)+y(3);mu*(1-y(1)^2)*y(2)-y(1);(Q11+Q41-Q21-Qe1)/24];

然后再运用ODE45求解：tspan=[0 20];yo=[2;0;30];
[t,y]=ode45(@vdpo1,tspan,yo)
size(t)
plot(t,y(:,1),t,y(:,2),'--')
请问为什么一个简单的微分方程组无法求解？谢谢~！

zlwr454 · 2009-11-08, 22:04

有没有哪位可以解答一下呢~！谢谢~！

zlwr454 · 2009-11-08, 22:06

有没有哪位可以解答一下呢~！谢谢啊！

zlwr454 · 2009-11-10, 22:51

没有人可以解答一下嘛.....

laosam280 · 2009-11-22, 15:46

时间范围取得太大了
因为龙格库塔法的步长很小，0到20计算量是有点大

2009-11-08, 21:31	#1
zlwr454 初级会员注册日期: 2009-11-08 帖子: 5 声望力: 0	请问为什么这个常微分方程组运行老是BUSY..... function ydot=vdpo1(t,y) Q11=2500; Q21=1800; Q41=700; Qe1=500; mu=2; ydot=[y(2)+y(3);mu(1-y(1)^2)y(2)-y(1);(Q11+Q41-Q21-Qe1)/24]; 然后再运用ODE45求解：tspan=[0 20];yo=[2;0;30]; [t,y]=ode45(@vdpo1,tspan,yo) size(t) plot(t,y(:,1),t,y(:,2),'--') 请问为什么一个简单的微分方程组无法求解？谢谢~！

2009-11-08, 22:04	#2
zlwr454 初级会员注册日期: 2009-11-08 帖子: 5 声望力: 0	回复: 请问为什么这个常微分方程组运行老是BUSY..... 有没有哪位可以解答一下呢~！谢谢~！

2009-11-08, 22:06	#3
zlwr454 初级会员注册日期: 2009-11-08 帖子: 5 声望力: 0	回复: 请问为什么这个常微分方程组运行老是BUSY..... 有没有哪位可以解答一下呢~！谢谢啊！

2009-11-10, 22:51	#4
zlwr454 初级会员注册日期: 2009-11-08 帖子: 5 声望力: 0	回复: 请问为什么这个常微分方程组运行老是BUSY..... 没有人可以解答一下嘛.....

2009-11-22, 15:46	#5
laosam280 高级会员注册日期: 2008-11-07 住址: 湖南长沙帖子: 233 声望力: 21	回复: 请问为什么这个常微分方程组运行老是BUSY..... 时间范围取得太大了因为龙格库塔法的步长很小，0到20计算量是有点大 __________________ 坚持就是胜利，努力就有奇迹。