SIR模型matlab程序ill参数的数目不足

ts=0:50

x0=[0.02,0.98]

[t,x]=ode45('ill',ts,x0)

plot(t,x(:,1),t,x(:,2)),grid,pause

plot(x(:,2),x(:,1)),grid,

你的灶颤程序后半部分谈伏，按上含辩携面语句试试

function y=ill(t,x)

a=1b=0.3

y=[a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)]'高闷裂

ts=0:50

x0=[0.02,0.98]

[t,x]=ode45('ill',ts,x0)

plot(t,x(:,1),t,x(:,2)),grid,pause

plot(x(:,2),x(:,1)),grid on

end

以上是sir模罩纯型的方程，参数自己设置！戚闭

根据ode45()函数对微亮虚备分方程的形式要求, y=ill(t,x)实际上是 y'誉明=ill(t,x), y必须是列向量

y(1) = dy1/dt, y(2)=dy2/dt

故这敬毁一行就是 dy1/dt = a* y1*y2-b*y1, dy2/dt= -a*y1*y2

正是sir模型方程(y1=i, y2=s)

---