传染病模型

附录

附件1：matlab模拟与实际累计病人程序 >> x=[1:1:65];

>> y=[402 610 666 782 863 954 1093 1255 1275 1358 1408 ... 1553 1636 1741 1803 1897 1960 2049 2136 2177 2227 ... 2265 2304 2347 2370 2388 2405 2420 2434 2437 2444 ... 2444 2456 2465 2490 2499 2504 2512 2514 2517 2520 ...

2521 2522 2522 2522 2522 2522 2522 2523 2522 2522 2522 ...

2523 2523 2522 2522 2522 2521 2521 2521 2521 2521 2521 2521 2521]; xi=1:1:64;

>> yi=polyval(p,xi); >> plot(x,y,'*r',xi,yi) >> p=polyfit(x,y,3); xi=1:1:64;

yi=polyval(p,xi); >> grid on

2.q值确定的模拟

y=[0.150 0.141 0.125 0.130 0.133 0.131 0.126 0.120 0.114 0.110 0.115 0.117 0.122 0.121 ... 0.1209 0.1180 0.1230 0.1225 0.1236 0.1979 0.1307 0.1351 0.1462 0.1611 0.1650 0.1662 ...

0.1721 0.1868 0.1968 0.2048 0.2246 0.2455 0.5692 0.3010 0.333 0.3485 0.3654 0.3980 0.4140 0.4386 ... 0.4694 0.5030 0.5174 0.5305 0.5432 0.5726 0.5956 0.6289 0.6457 0.7664 0.7283 0.7665 0.7954 ... 0.8176 0.8450 0.8656 0.8739 0.8897 0.9163 0.9302 0.9369 0.9440 0.9607 0.9710 0.9790]; xi=1:1:65;

?yi=polyval(p,xi); ?plot(x,y,'*r',xi,yi) ?p=polyfit(x,y,3) p =

-0.0000 0.0008 -0.0183 0.2024 ?xi=1:1:65;

?yi=polyval(p,xi); ?plot(x,y,'*r',xi,yi) ?grid on

p(x)=0.00008x^2-0.01813x+0.2024 为一个2次函数

17