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，如果是时间的话，就是常微分方程初值问题，需要给出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
为什么我改成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;