言能够通过传播者−免疫者恢复或者通过外部干预恢复,前者称为点对点的恢复过程,后
者称为干预恢复过程。除了层内传播和恢复之外,假设谣言还可以在层间传播和恢复。为
了方便,把子网络中的 3 类人群标记为 3 种状态:易感态(SS)、感染态(II)以及恢复态(RR)。
因此,谣言在层内和层间的传播动力可用随机切换事件来描述:
Si+Ii−→βi2IiSi+Ij−→γjIi+IjSi+Ii→βi2IiSi+Ij→γjIi+Ij
Ii+Ri−→αi2RiIi+Rj−→νjRi+RjIi+Ri→αi2RiIi+Rj→νjRi+Rj
Ii−→−μi(t)Rit⩾t0Ii→μi(t)Rit⩾t0
式中,i=1i=1 时,j=2j=2;i=2i=2 时,j=1j=1;式(1)表示层内和层间的传播过程;式
(2)表示层内和层间恢复过程;式(3)表示来自于外部控制措施下的谣言恢复过程;t0t0 表示
开始采取控制措施的时刻。控制意味着人力、物力、财力等资源的投入,从节约资源的角
度,假设采取的控制措施能够影响层内的谣言传播,即层内干预恢复。相应地,
βi(αi)βi(αi)表示层内传播(恢复)速率,γj(νj)γj(νj)表示层间传播(恢复)速率,μi(t)μi(t)表示来
自外部干预下的恢复速率。
假设 PS1,i(t)(PS2,i(t))P1,iS(t)(P2,iS(t)),PI1,i(t)(PI2,i(t))P1,iI(t)(P2,iI(t)),
PR1,i(t)(PR2,i(t))P1,iR(t)(P2,iR(t))分别表示 tt 时刻子网络 1 和 2 上个体 ii 处于 3 个状态的概
率,则概率随时间的演化可用离散动力方程表示为:
PSk,i(t+1)=PSk,i(t)(1−qS,Ik,i(t))PIk,i(t+1)=PIk,i(t)(1−qI,Rk,i(t))+PSk,i(t)qS,Ik,i(t)PRk,i(t+1)=PRk,i(t)+PIk,i(t)qI,Rk,i(t)Pk,iS(t+1)=Pk,iS(t)(1−qk,iS,I(t))Pk,iI(t+1)=Pk,iI(t)(1−qk,iI,R(t))+Pk,iS(t)qk,iS,I(t)Pk,iR(t+1)=Pk,iR(t)+Pk,iI(t)qk,iI,R(t)
式中,k=1,2k=1,2;qS,Ik,i(t)qk,iS,I(t)和 qI,Rk,i(t)qk,iI,R(t)分别是子网络 kk 上从 SS 到
I,以及从 I 到 R 的转移概率。Ω1,iΩ1,i 和 Ω2,iΩ2,i 分别表示个体 ii 在子网络 1 和 2 上的邻
居集。根据微观马氏链方法
[20-22]
,对子网络 1,有:
qS,I1,i(t)=1−∏j∈Ω1,i(1−β1PI1,j(t))∏k∈Ω2,i(1−γ2PI2,k(t))qI,R1,i(t)=1−∏j∈Ω1,i(1−α1PR1,j(t))∏k∈Ω2,i(1−ν2PR2,k(t))×(1−aμ1(t))q1,iS,I(t)=1−∏j∈Ω1,i(1−β1P1,jI(t))∏k∈Ω2,i(1−γ2P2,kI(t))q1,iI,R(t)=1−∏j∈Ω1,i(1−α1P1,jR(t))∏k∈Ω2,i(1−ν2P2,kR(t))×(1−aμ1(t))
式中,∏j∈Ω1,i(1−β1PI1,j(t))∏j∈Ω1,i(1−β1P1,jI(t))和∏k∈Ω2,i(1−γ2PI2,k(t))∏k∈
Ω2,i(1−γ2P2,kI(t))分别表示个体 ii 没有受到子网络 1 和子网络 2 中谣言传播者影响的概
率;∏j∈Ω1,i(1−α1PR1,j(t))∏j∈Ω1,i(1−α1P1,jR(t))和∏k∈Ω2,i(1−ν2PR2,k(t))∏k∈
Ω2,i(1−ν2P2,kR(t))分别表示个体 ii 没有受到子网络 1 和子网络 2 中恢复者影响的概率;
(1−aμ1(t))(1−aμ1(t))表示在外部控制措施出现时,仍没有恢复的概率。a∈{0,1}a∈{0,1},
其中 a=0a=0 表示没有控制措施,a=1a=1 表示采取了控制措施。类似的,对网络层 2,
有:
qS,I2,i(t)=1−∏j∈Ω2,i(1−β2PI2,j(t))∏k∈Ω1,i(1−γ1PI1,k(t))qI,R2,i(t)=1−∏j∈Ω2,i(1−α2PR2,j(t))∏k∈Ω1,i(1−ν1PR1,k(t))×(1−aμ2(t))q2,iS,I(t)=1−∏j∈Ω2,i(1−β2P2,jI(t))∏k∈Ω1,i(1−γ1P1,kI(t))q2,iI,R(t)=1−∏j∈Ω2,i(1−α2P2,jR(t))∏k∈Ω1,i(1−ν1P1,kR(t))×(1−aμ2(t))
分析可知,式(4)的第一行的右边表示个体 ii 仍是易感态的概率;第二行右边的第一项
表示个体 ii 仍是感染态的概率,第二项是从易感态变为感染态的概率;第三行右边第一项
是个体 ii 仍然是恢复态的概率,第二项是从感染态变为恢复态的概率。此外,网络层上 3
类群体的密度可表示为:
ρkS(t)=1N∑i=1NPSk,i(t)ρkI(t)=1N∑i=1NPIk,i(t)ρkR(t)=1N∑i=1NPRk,i(t)ρSk(t)=1N∑i=1NPk,iS(t)ρIk(t)=1N∑i=1NPk,iI(t)ρRk(t)=1N∑i=1NPk,iR(t)
需要指出的是,由于采取的控制策略是动态的,因而式(5)和式(6)中的干预恢复速率
μ1(t)μ1(t)和 μ2(t)μ2(t)是关于时间的函数。通常,当一部分群体受到谣言影响之后,谣言
才会引起有关部门的注意,进而采取措施控制谣言传播。换句话说,从谣言开始出现到采