shw1220 发表于 2012-4-18 15:56:48

求助:ode45解微分方程组提前中断

m文件为:
function dydt=vdt1(t,y)
dydt=[(7.7794e-14)*360000*sinh(y(2)/21.7455);
y(2)*(7.7794e-14)*360000*sinh(y(2)/21.7455)];
计算命令为:
=ode45(@vdt1,,);
plot(t,y(:,1))
图像显示计算到t为9000多的时候就中断了,改成ode23s和ode15s都显示Unable to meet integration tolerances without reducing the step size below the smallest value allowed at time t.
本人新手,最近需要计算这个方程组才翻开matlab的,诚请各位帮助,谢谢!
页: [1]
查看完整版本: 求助:ode45解微分方程组提前中断