微分方程错误..Error using ==> feval

[复制链接]

字体大小: 正常放大

慢跑20

60 主题	8 听众	3684 积分

TA的每日心情

	开心 2017-2-22 14:21

签到天数: 271 天

[LV.8]以坛为家I

群组: 2014年美赛冲刺培训

群组: 物联网工程师考试

群组: 2013年电工杯B题讨论群

群组: 物联网工程师培训

群组: 2013电工杯A题讨论群组

电梯直达

1^#

发表于 2013-5-11 17:18 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

课本传染病模型SIR.
先建立M文件
function y=ill(t,x)
a=1;b=0.3;
y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]';
然后再MATLAB中输入:
ts=0:50;
x0=[0.02,0.98];
[t,x]=ode45('ill',ts,x0);[t,x]
Plot(t,x(:1,),t,x(:,2)),grid,pause
Plot(x(:,2),x(:,1)),grid
function y=ill(t,x)
a=1;b=0.3;
y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]';

就会:
显示错误:
如下:
??? Error using ==> feval
Undefined command/function 'ill'.

Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...

>>
啥问题呢?

zan

Error