SI 模型
利用 MATLAB 求解传染病模型中的 SI 模型的解析解:
程序中 a 即 λ,y 即 i
>> y=dsolve('Dy=a*(y-y^2)','y(0)=y0')
y =
1/(1-exp(-a*t)*(-1+y0)/y0)
画图:SI 模型的 i~t 曲线
设 λ=1, i(0)=0.1
>> y=dsolve('Dy=y-y^2','y(0)=0.1')
y =
1/(1+9*exp(-t))
>> x=0:0.01:13;
y=1./(1+9.*exp(-x));
>> plot(x,y)
title('SI 模型的 i~t 曲线');
xlabel('t');
ylabel('i');
axis([0 13 0 1.1]);