matlab病毒扩散,[MATLAB数学相关] 现有Allen的病毒扩散模型SIRS的马尔科夫过程程序，求指导...

现有Allen的病毒扩散模型SIRS的马尔科夫过程程序，求指导

由Allen编写的SIRS模型的马尔可夫过程如下：

for k=1:4 % 解马尔科夫链的样本轨迹

clear t s i r

t(1)=0; i(1)=i0; s(1)=N-i0;r(1)=0;

j=1;

while i(j)>0 && t(j)

nn=s(j)+i(j)+r(j);

u1=rand;u2=rand;

den=((beta/nn)*i(j)*s(j)+gam*i(j)+alpha*i(j)+delta*r(j));%所有patch的概率之和

t(j+1)=-log(u1)/den+t(j); % Time to next event

e1=(beta/nn)*s(j)*i(j)/den;%以下是每个patch概率与总概率之比

e2=e1+gam*i(j)/den;

e3=e2+alpha*i(j)/den;

e4=e3+delta*r(j)/den;

if (u2<=e1) %每个patch情况下的个体变化情况

s(j+1)=s(j)-1;

i(j+1)=i(j)+1;

r(j+1)=r(j);

elseif (u2>e1 && u2<=e2)

s(j+1)=s(j);

i(j+1)=i(j)-1;

r(j+1)=r(j)+1;

elseif (u2>e2 && u2<=e3)

s(j+1)=s(j);

i(j+1)=i(j)-1;

r(j+1)=r(j);

else

s(j+1)=s(j)+1;

i(j+1)=i(j);

r(j+1)=r(j)-1;

end

j=j+1;

end

(省略出图部分)

end

本人菜鸟，对上述标黄部分的分支概率不太理解，变量“den”是各个转换概率之和，e1-e4分别是四种转换概率与“den”的比值，请问各位大神，e1-e4的顺序有区别么？另附本程序的转换图，见“SIR模型”图，若想在本程序的基础上，改写成“SIRP模型”图所应用的程序，应该怎么写呢？

目前存在两个问题：1.S→I之间有两种转换概率的情况该怎么解决呢？

2. P作为独立变量，对S→I有影响，但其增加或减少的变动情况与S,I,R相对独立，P的每一步数量应该单独统计么？

感谢各位！论文卡在这里很久了，请大家帮忙。

SIR模型.jpg

SIRP模型.jpg