5楼: Originally posted by zhuwx1999 at 2014-11-18 08:35:17
不好意思，再问一个问题：
根据实验，y-t曲线应该出现一个峰值（dy/dt小于等于0），不知是不是计算的时候有什么假设条件？...

clear all;clc
a=0.01162;
b=0.00135;
c=0.0000051246;
[t y]=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, [0:100],0.148   );
[t y]
figure(1)
plot(t,y,'b-o')

clear all;clc
a=0.01162;
b=0.00135;
c=0.0000051246;
[t y]=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, [0:1000],0.148   );
[t y]
figure(1)
plot(t,y,'b-o')

3楼: Originally posted by 月只蓝 at 2014-11-17 20:40:45
随着时间的推移，即增大计算时间，y会稳定某一数值：
clear all;clc
a=0.01162;
b=0.00135;
c=0.0000051246;
=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, ,0.148   );

figure(1)
plot(t,y,'b-o')


...