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ÊDz»ÊÇ´ú±í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£©£¬²»ÖªÊDz»ÊǼÆËãµÄʱºòÓÐʲô¼ÙÉèÌõ¼þ£¿...

6Â¥: Originally posted by ÔÂÖ»À¶ at 2014-11-18 08:52:10
ÎÊÌâ1
´ð£ºÊǵģ¬ÕýÊdzõÖµ¡£
ÎÊÌâ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µÄÊýÖµ¡£...