传染病模型在计算机病毒传播中的应用(10)
同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
下一篇：利用插值多项式求解函数问题的综述