- 上一篇：Black-Scholes方程的求解方法分析及应用
- 下一篇：利用插值多项式求解函数问题的综述
同SIR模型一样，假设网路中的节点总数为1000000，最初被感染的计算机台数 ， =0。感染率以及排除率不变。在该模型中，由于采取了了反病毒措施，所以假设采用反病毒措施的概率 。
（2）模型编程
function y=ill(t,x)
N=10^6;
a=0.8/N;b=0.05;k=0.005;
y=[-a*x(1)*x(2)-k*x(1),a*x(1)*x(2)-(b+k)*x(2),k*x(1),(b+k)*x(2)]';
>> ts=0:40;
>> x0=[10^6-1,1,0,0];
>> [t,x]=ode45('ill',ts,x0);[t,x]
>> plot(t,x(:,2)+x(:,4)),grid,pause
（3）程序结果
ans =
1.0e+005 *
0 10.0000 0.0000 0 0
0.0000 9.9501 0.0000 0.0499 0.0000
0.0000 9.9005 0.0000 0.0995 0.0000
0.0000 9.8510 0.0001 0.1489 0.0000
0.0000 9.8018 0.0002 0.1980 0.0000
0.0001 9.7527 0.0004 0.2469 0.0000
0.0001 9.7036 0.0008 0.2955 0.0001
0.0001 9.6543 0.0017 0.3439 0.0001
0.0001 9.6042 0.0034 0.3921 0.0003
0.0001 9.5525 0.0070 0.4400 0.0005
0.0001 9.4972 0.0141 0.4876 0.0011
0.0001 9.4344 0.0285 0.5349 0.0022
0.0001 9.3565 0.0571 0.5819 0.0045
0.0001 9.2487 0.1139 0.6284 0.0090
。。。
0.0002 1.3907 6.4180 0.9113 1.2800
0.0002 0.8194 6.6239 0.9167 1.6400
0.0002 0.4797 6.5960 0.9199 2.0043
0.0003 0.2833 6.4319 0.9218 2.3631
0.0003 0.1700 6.1966 0.9229 2.7106
0.0003 0.1041 5.9282 0.9235 3.0441
0.0003 0.0652 5.6483 0.9239 3.3625
0.0003 0.0418 5.3686 0.9242 3.6654
0.0003 0.0273 5.0951 0.9244 3.9531
0.0003 0.0183 4.8311 0.9245 4.2261
0.0003 0.0125 4.5782 0.9246 4.4848