2Â¥: Originally posted by ÔÂÖ»À¶ at 2013-11-27 18:34:09
ûÓнâÎö½â¡£¿ÉÒÔÇóÊýÖµ½â¡£

2Â¥: Originally posted by ÔÂÖ»À¶ at 2013-11-27 18:34:09
ûÓнâÎö½â¡£¿ÉÒÔÇóÊýÖµ½â¡£

4Â¥: Originally posted by watertxf at 2013-11-27 19:07:38
tµÄÈ¡Öµ·¶Î§Îª1~5...

function solve_ode
clear all;clc
y0=0;
tspan=1:0.01:5;
[ts ys]=ode15s(@odefun,tspan,y0);
figure(1),
plot(ts,ys,'-ob')


function dydt=odefun(t,y)
dydt=1-(1+16*(1-exp(-t)-y))*y;

5Â¥: Originally posted by ÔÂÖ»À¶ at 2013-11-27 19:09:07
×Ô±äÁ¿ÊÇt£¬Èç¹ûÊÇʱ¼äµÄ»°£¬¾ÍÊdz£Î¢·Ö·½³Ì³õÖµÎÊÌ⣬ÐèÒª¸ø³öt=0ʱµÄyÖµ£¬ÒÔ¼°Òª»­Í¼ÇúÏßµÄtµÄÈ¡Öµ·¶Î§£¬È»ºóÓÃode45º¯Êý¼´¿É½â¡£

6Â¥: Originally posted by ÔÂÖ»À¶ at 2013-11-27 19:21:53
Èç¹ût=0ʱ£¬y=0µÄ»°£¬³ÌÐòºÍ½á¹ûÈçÏ£º

function solve_ode
clear all;clc
y0=0;
tspan=1:0.01:5;
=ode15s(@odefun,tspan,y0);
figure(1),
plot(ts,ys,'-ob')


function dydt=odefun(t,y)
dydt=1-(1 ...

6Â¥: Originally posted by ÔÂÖ»À¶ at 2013-11-27 19:21:53
Èç¹ût=0ʱ£¬y=0µÄ»°£¬³ÌÐòºÍ½á¹ûÈçÏ£º

function solve_ode
clear all;clc
y0=0;
tspan=1:0.01:5;
=ode15s(@odefun,tspan,y0);
figure(1),
plot(ts,ys,'-ob')


function dydt=odefun(t,y)
dydt=1-(1 ...

9Â¥: Originally posted by watertxf at 2013-11-27 19:30:05
tµÄ·¶Î§²»ÊÇ1~5£¬ÊÇ0~5
ΪʲôÎҸijÉtspan=0:0.01:5;ºóµÃ³öµÄ½á¹û³ÉϱߵÄͼ1ÁË£¿ºÜÃ÷ÏÔ²»¶ÔµÄ£¬tµÄ·¶Î§¾Í²»¶Ô.ÕâÊÇʲôԭÒòÄØ£¿Ð»Ð»£¡

D_V9`YP`W~TCPD[5EA(@AZT.jpg
...

function solve_ode
clear all;clc
y0=0;
tspan=0:0.01:5;
[ts ys]=ode45(@odefun,tspan,y0);
figure(1),
plot(ts,ys,'-ob')


function dydt=odefun(t,y)
dydt=1-(1+16*(1-exp(-t)-y))*y;