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')


...

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')


...

4楼: Originally posted by zhuwx1999 at 2014-11-18 08:15:01
谢谢你！
我还有两个小问题：
1、=ode45(@(t,y)a*y^0.266-b*y^3.532+c*y^2.266, ,0.148   )公式中的0.148是不是代表dy/dt的初值？
2、另外，想把公式中的0.266改为-0.266，用y^（-0.266）还是直接y^-0.266？...

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

6楼: Originally posted by 月只蓝 at 2014-11-18 08:52:10
问题1
答：是的，正是初值。
问题2
答：y^(-0.266)比较保险；此外注意要在英文状态下输入。...

8楼: Originally posted by zhuwx1999 at 2014-11-18 13:26:05
可能我没有表达清楚。
问题1：公式中0.148表示的是y的初值还是dy/dt的初值？...

9楼: Originally posted by 月只蓝 at 2014-11-18 14:03:05
对于一阶常微分方程，初值专指t=0时刻，y的数值。...