官方媒体影响下的随机谣言传播模型
Stochastic Rumor Propagation Model under the Influence of Official Media
摘要: 在添加官方媒体对谣言传播易感染者影响下的确定性方程中,考虑影响因子γ受环境噪声的扰动,建立官方媒体影响下的随机微分方程,并研究了在无谣言平衡点和正平衡点附近的渐近性质。证明了该随机微分方程的正解存在唯一性。然后,定义了新的基本再生数ℜ'0,并证明了满足一定条件时,系统的解在无谣言平衡点附近是渐近稳定的;而当满足另外四点要求时,系统的解在正谣言平衡点附近是渐近稳定的。
Abstract: In addition to the deterministic equation under the influence of the official media on the rumor propagation of susceptible infected people, the stochastic differential equation under the influence of the official media is established by considering the disturbance of the influence factor by environmental noise, and the asymptotic properties near the rumor free equilibrium point and positive equilibrium point are studied. First, the existence and uniqueness of positive solutions of the stochastic differential equation are proved. Then, a new basic reproduction number '0 is defined. And it is proved that when certain conditions are satisfied, the solution of the system is asymptotically stable near the rumor free equilibrium point. When the other four requirements are met, the solution of the system is asymptotically stable near the positive rumor equilibrium point.
文章引用:胡悦, 胡良剑. 官方媒体影响下的随机谣言传播模型[J]. 应用数学进展, 2020, 9(12): 2317-2328. https://doi.org/10.12677/AAM.2020.912271

1. 引言

2019年10月开始报道武汉不明原因的肺炎,再到后面陆续大量的病人被确诊,网络上各种关于新冠的谣言顿时成为微博、微信和抖音等热门话题。2020年1月31日,世界卫生组织在新闻发布会上正式宣布:新冠疫情被列入“国际关注的突发公共卫生事件”,瞬间网络上就有营销博主称这意味着中国将被列入“疫区国”,这引起民众极大的恐慌;更有盐水漱口能预防感染、电风扇对手和脸部吹30秒能消毒、吃维生素C能预防新冠病毒等谣言,这些迅速在网络平台,特别在是长辈群里传播,造成维生素C等无明显作用药物的哄抢。由于网络的便捷、对用户的评论和发表言论的管理较为松懈,就会大大降低网络消息的真实性,极大地影响了社会的稳定。

1965年Daley和Kendall [1] 在其论文中利用新的方法——“任意常数扩散原理”,首次研究得到:从表面上看,谣言与流行病两者之间有着相似之处,提出了谣言传播模型,这个模型就是“D-K”模型。1973年Maki和Thomson [2] 在D-K模型的基础上,认为第一个人知道并传播谣言,才会有可能知道谣言并不传播谣言,得到了M-T模型,该模型认为这类人不会受到谣言的干扰。2001年,Zanette [3] 在《在小世界网络上传播的关键行为》文章中首次用实验数据证明了一个类似流行病的模型,这个模型可以被解释为谣言传播模型,而且得到了该模型传播的临界指数,研究了临界点的随机性与网络的连通性之间的关系。

Zhao [4] 在《社交网络中的谣言传播模式》指出:当考虑到遗忘与记忆机制的相互作用时,谣言传播模型与流行病传播模型就会存在显著差异,于是提出新的谣言传播模型SIHR模型。该模型考虑到人群中的无知者到顽抗者的过程与一种新型冬眠者的相关联系,然后对方程的稳态进行分析,并研究了在传播率、窒息率、遗忘率和网络平均程度不同状态下的最终传播规模,得到了谣言传播的阈值。最后通过龙格–库塔方法进行数值模拟。Zhu等 [5] 在最基础的SIR模型上,改进了新的谣言传播模型,用于研究复杂社会网络中的谣言传播问题。该论文中考虑连通性随环境的扰动,建立了均匀网络和非均匀网络中谣言传播的随机微分方程,分析新建模型解的性质以及该系统的稳定性。最后通过数据模拟给出了均匀网络下谣言扩散的阈值,并发现在异质网络中无谣言传播的阈值。

Liu等 [6] 提出了一种新的异质网络谣言传播模型,分析了系统的全局动力学行为,给出了该系统的谣言传播的阈值。陈安等 [7] 人以新冠肺炎为例,建立了以SEIR为基础的伪科学网络谣言传播模型,通过具体的数值分析,分析了伪科学网络谣言传播在本次新冠疫情期间的传播特点,给出了理论性的建议。

从谣言传播的研究结果来看,这些模型主要分为三种:第一种是基于传染病模型的改进,第二种是在复杂网络或者小世界网络上的谣言传播模型,第三种是以伪科学这类谣言传播的相关模型。实际上,谣言的传播还与媒体的报道效果以及媒体的官方性有关,其中媒体的相关报道对谣言传播起到了至关重要的作用。宋雪等 [8] [9] 对受媒体正面报道和负面报道的谣言传播模型进行了研究,定义了系统的阈值,讨论了系统平衡点附近的动力学行为。在赵敏 [10] 的确定性模型中,考虑到官方媒体宣传的影响因素,因此原本的易感人群可分为两类,一是未受到官方媒体宣传影响的易感人群,二是受到官方媒体宣传的易感人群,建立了如下的微分方程:

{ d S u ( t ) d t = q Λ β S u ( t ) I ( t ) γ S u ( t ) d S e ( t ) d t = ( 1 q ) Λ β σ S e ( t ) I ( t ) γ S u ( t ) μ S e ( t ) d I ( t ) d t = β S u ( t ) I ( t ) + β σ S e ( t ) I ( t ) γ I ( t ) μ I ( t ) d R ( t ) d t = μ S e ( t ) + μ I ( t ) γ R ( t ) (1)

其中: S u ( t ) 表示易感人群中未受官方媒体影响下的人群数量; S e ( t ) 表示易感人群中受官方媒体影响下的人群数量; I ( t ) 表示谣言传播者的数量; R ( t ) 表示谣言免疫者的数量; Λ 表示该系统原有的总人数; σ 表示媒体宣传程度的有效性; γ 表示影响因子:当到达一定时间段的时候,系统里面的人对社会状况的认知会发生改变,从而影响系统中的四类人离开此系统。作者研究了模型(1)的无谣言平衡点和正平衡点附近的动力学行为。

结合上述分析,本文在官方媒体影响下的谣言传播模型的基础上,考虑影响因子 γ 随环境变化而发生改变,不再是个固定的数值。所以,本文考虑把影响因子 γ 添加一个随机扰动项,即 γ γ + σ 1 B ( t ) ,那么就能得到如下随机微分方程为:

{ d S u ( t ) = ( q Λ β S u ( t ) I ( t ) γ S u ( t ) ) d t σ 1 S u ( t ) d B ( t ) d S e ( t ) = ( ( 1 q ) Λ β σ S e ( t ) I ( t ) γ S e ( t ) μ S e ( t ) ) d t σ 1 S e ( t ) d B ( t ) d I ( t ) = ( β S u ( t ) I ( t ) + β σ S e ( t ) I ( t ) γ I ( t ) μ I ( t ) ) d t σ 1 I ( t ) d B ( t ) d R ( t ) = ( μ S e ( t ) + μ I ( t ) γ R ( t ) ) d t σ 1 R ( t ) d B ( t ) (2)

根据确定性模型的假设 [10],这里也有如下规定:

1) 在这个系统中的所有人都将被认为是易感人群;

2) 这里 Λ 表示的是整个系统的总人数,那么在这个系统中, q Λ 表示最初该系统中受到官方媒体影响的易感人群的数量,而 ( 1 q ) Λ 就表示该系统中未受官方媒体影响的易感人群最初数量;

3) 受官方媒体影响的易感者以 σ β 的概率变成为谣言传播者,未受官方媒体影响的易感人群以概率 β 转化为谣言传播者;

4) 一段时间后,考虑系统中有一定比例 γ 的人会由于社会认知等原因而离开系统;

5) 在谣言传播的整个过程中,会因为其他因素或环境等变化等影响,受到官方媒体影响的易感者和谣言传播者这两类人的数量,将会以概率 μ 转变成谣言免疫者。

上述过程就可以用如下文献 [10] 流程图表示(图1):

Figure 1. Schematic diagram of rumor propagation process

图1. 谣言传播过程示意图

2. 正解的存在唯一性

在研究系统(2)的动力学行为之前,首先需要研究其正解的存在唯一性,这样才能保证方程的现实意义。显然系统(2)的系数是局部李普希兹连续的,但不满足线性增长条件,因此一般的存在唯一性将不再适用于此过程的证明,于是这里采用构造李雅普诺夫函数的方法进行正解存在唯一性的证明。

定理1对于任意给定的初值 ( S u ( 0 ) , S e ( 0 ) , I ( 0 ) , R ( 0 ) ) 系统存在唯一正解 ( S u ( t ) , S e ( t ) , I ( t ) , R ( t ) ) 对所有 t 0 ,并且解以概率1存在。

证明 显然随机微分方程组(2)的系统满足局部 李普希兹条件。对任意给定的初值,有最大的唯一解 ( S u ( t ) , S e ( t ) , I ( t ) , R ( t ) ) [ 0, τ ε ) ,其中 τ ε 是爆破时间,令 k 0 充分大,且 k 0 > 0 ,对 ( S u ( 0 ) , S e ( 0 ) , I ( 0 ) , R ( 0 ) ) [ 1 k 0 , k 0 ) 中。对每个正数 k > k 0 ,定义停时 τ k = inf { t [ 0 , τ e ) : S u ( t ) S e ( t ) I ( t ) R ( t ) ( 1 k , k ) } 。当 k 时, τ k 单调增长,令 τ lim k τ k ,显然 τ τ ε a . s . 。因此只需证 τ = a . s . 。否则,存在常数 T > 0 k 1 > k 0 ε ( 0 , 1 ) 使得 P ( Ω k ) ε ,对所有的 k > k 1 成立,其中 Ω k = { τ k T }

构造李雅普诺夫函数

V ( S u ( t ) , S e ( t ) , I ( t ) , R ( t ) ) = S u + a a ln S u a + S e + a a ln S e a + I + 1 ln I + R + 1 ln R

根据伊藤公式得到

L V = ( 1 a S u ) ( q Λ β S u I γ S u ) + a 2 σ 1 2 + ( 1 a S e ) ( ( 1 q ) Λ β σ S e I γ S e μ S e ) + a 2 σ 1 2 + ( 1 1 I ) ( β S u I + β σ S e I γ I μ I ) + 1 2 σ 1 2 + ( 1 1 R ) ( μ S e + μ I γ R ) + 1 2 σ 1 2 Λ γ S u γ S e γ I γ R + a β I + 2 a γ + a σ 1 2 + a β σ I + σ 1 2 + a μ + μ + 2 γ Λ + 2 a γ + a σ 1 2 + σ 1 2 + a μ + μ + 2 γ + ( 2 a β γ ) I

2 a β γ = 0

a = γ 2 β

因此 L V Q ,其中

Q = Λ + 2 a γ + a σ 1 2 + σ 1 2 + a μ + μ + 2 γ

E V ( S u ( T τ k ) , S e ( T τ k ) , I ( T τ k ) , R ( T τ k ) ) V 0 + E 0 T τ k Q d s ,

E V ( S u ( T τ k ) , S e ( T τ k ) , I ( T τ k ) , R ( T τ k ) ) V 0 + Q T

在Ωk内, τ k T S u ( τ k ) S e ( τ k ) I ( τ k ) R ( τ k ) 至少有一个等于k 1 k ,有

V ( S u ( τ k ) , S e ( τ k ) , I ( τ k ) , R ( τ k ) ) ( k + 1 ln k ) ( 1 k + 1 + ln k ) ( k + a a ln k a ) ( 1 k + a + a ln k a )

因此可以得到:

V 0 + Q T E [ 1 Ω k V ( S u ( τ k ) , S e ( τ k ) , I ( τ k ) , R ( τ k ) ) ] P ( Ω k ) ( k + 1 ln k ) ( 1 k + 1 + ln k ) ( k + a a ln k a ) ( 1 k + a + a ln k a )

k 时, V 0 + Q T = ,这与 V 0 + Q T < 产生矛盾,因此 τ = a . s

系统(1)的边界平衡点为 E 0 = ( q Λ γ , ( 1 q ) Λ γ + μ , 0 , μ Λ ( 1 q ) γ ( γ + μ ) ) ,其含义是指网络中没有谣言进行传播,也叫做无谣言传播平衡点。在文献 [10] 中,研究了此无谣言传播平衡点的存在性,并定义了基本再生数 0 ,当 0 < 1 时,系统(1)的谣言传播就会逐渐消失。那么,在随机微分方程中,考虑了噪声的影响,边界平衡点仍然存在,那么是否可以定义新的基本再生数,使得满足相应条件时,边界平衡点仍具有某种性质?下面本文研究系统(2)在边界平衡点附近解的性态。

3. 在边界平衡点附近解的渐近性态

定理2若 0 = β q Λ γ ( γ + μ ) + β σ Λ ( 1 q ) ( γ + μ ) 2 < 1 且满足条件:

1) ( 1 + c 1 ) σ 1 2 < c 1 ( γ + μ ) + γ

2) 2 σ 1 2 + 1 2 < 2 γ

则对于任何给定的初值 ( S u ( 0 ) , S e ( 0 ) , I ( 0 ) , R ( 0 ) ) R + 4 ,随机微分方程满足如下性质:

lim t sup 1 t E 0 t [ ( S u ( s ) q Λ γ ) 2 + ( S e ( s ) ( 1 q ) Λ γ + μ ) 2 + I ( s ) 2 + ( R ( s ) ( 1 q ) Λ μ ( γ + μ ) γ ) 2 ] d s A B ,

其中

c 1 = 2 γ 2 β q Λ σ > 0 c 2 = 4 γ β σ > 0 A = 2 ( 1 + c 1 ) σ 1 2 q 2 Λ 2 γ 2 + c 2 μ 2 + 2 ( 1 + c 1 ) σ 1 2 ( 1 q ) 2 Λ 2 ( γ + μ ) 2 + 2 σ 1 2 ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2

B = min { 2 γ 2 ( 1 + c 1 ) σ 1 2 + 2 c 1 γ , 2 γ 2 ( 1 + c 1 ) σ 1 2 + 2 c 1 ( γ + μ ) 2 c 2 , 2 γ 2 σ 1 2 2 c 2 , 2 γ 2 σ 1 2 }

证明:令

G ( t ) = S u ( t ) q Λ γ H ( t ) = S e ( t ) ( 1 q ) Λ γ + μ Z ( t ) = I ( t ) M ( t ) = R ( t ) ( 1 q ) μ Λ ( μ + γ ) γ

则原系统(2)转化为:

{ d G ( t ) = ( β G Z ( t ) β q Λ Z γ γ G ) d t σ 1 ( G + q Λ γ ) d B 1 ( t ) d H ( t ) = ( β σ H Z β σ ( 1 q ) Λ Z γ + μ γ H μ H ) d t σ 1 ( H + ( 1 q ) Λ γ + μ ) d B 1 ( t ) d Z ( t ) = ( β G Z + β q Λ Z γ + β σ H Z + β σ ( 1 q ) Λ Z γ + μ γ Z μ Z ) d t σ 1 Z d B 1 ( t ) d M ( t ) = ( μ H + μ Z γ M ) d t σ 1 ( M + ( 1 q ) μ Λ ( μ + γ ) γ ) d B 1 ( t )

其中 G , H , M R , Z > 0 。对此我们考虑如下李雅普诺夫函数:

V ( G , H , Z , M ) = ( G + H + Z + M ) 2 + c 1 ( G + H ) 2 + c 2 ( Z + M )

其中 c 1 , c 2 是大于0的常数。由伊藤公式可得

L V = [ 2 ( G + H + Z + M ) + 2 c 1 G ] [ β G Z β q Λ Z γ γ G ] + ( 1 + c 1 ) σ 1 2 ( G + q Λ γ ) 2 + [ 2 ( G + H + Z + M ) + 2 c 1 H ] [ β σ H Z β σ ( 1 q ) Λ Z γ + μ γ H μ H ] + ( 1 + c 1 ) σ 1 2 ( H + ( 1 q ) Λ γ + μ ) 2 + σ 1 2 Z 2 + [ 2 ( G + H + Z + M ) + c 2 ] [ μ H + μ Z γ M ] + [ 2 ( G + H + Z + M ) + c 2 ] [ β G Z + β q Λ Z γ + β σ H Z + β σ ( 1 q ) Λ Z γ + μ γ Z μ Z ] + σ 1 2 ( M + ( 1 q ) Λ μ ( γ + μ ) γ ) 2

= ( 2 γ + ( 1 + c 1 ) σ 1 2 2 c 1 γ ) G 2 + ( 2 γ + ( 1 + c 1 ) σ 1 2 2 c 1 γ 2 c 1 μ ) H 2 + ( 2 γ + σ 1 2 ) Z 2 + ( 2 γ + σ 1 2 ) M 2 4 G H + ( 4 γ 2 c β 1 q Λ γ + c 2 β ) G Z 4 γ G M + ( 4 γ 2 c 1 β σ ( 1 q ) Λ γ + μ + c 2 β σ ) H Z 4 γ H M 4 γ Z M + 2 ( 1 + c 1 ) σ 1 2 q Λ γ G + ( c 2 μ + 2 ( 1 + c 1 ) σ 1 2 ( 1 q ) Λ γ + μ ) H + ( c 2 q Λ β γ + c 2 ( 1 q ) Λ β σ γ + μ c 2 γ + c 2 μ c 2 μ ) Z

+ 2 σ 1 2 ( 1 q ) μ Λ ( γ + μ ) γ M 2 c 1 β G 2 Z 2 c 1 β σ H 2 Z + ( 1 + c 1 ) σ 1 2 q 2 Λ 2 γ 2 + ( 1 + c 1 ) σ 1 2 ( 1 q ) 2 Λ 2 ( γ + μ ) 2 + σ 1 2 ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2

根据定理给出的条件可知

0 = β q Λ γ ( γ + μ ) + β σ Λ ( 1 q ) ( γ + μ ) 2 < 1

那么

β q Λ γ + β σ Λ ( 1 q ) γ + μ < γ + μ

因此

β q Λ γ + β σ Λ ( 1 q ) γ + μ ( γ + μ ) < 0

利用公式 a b 1 2 ( a 2 + b 2 ) 得到

L V [ 2 γ + ( 1 + c 1 ) σ 1 2 2 c 1 γ ] G 2 + [ 2 γ + ( 1 + c 1 ) σ 1 2 2 c 1 γ 2 c 1 μ ] H 2 + ( 2 γ + σ 1 2 ) Z 2 + ( 2 γ + σ 1 2 ) M 2 + ( 2 c 1 β q Λ γ + c 2 β ) G Z + ( 4 γ + c 2 β σ ) H Z + ( 1 + c 1 ) σ 1 2 ( q 2 Λ 2 γ + G 2 ) + c 2 μ H + 2 ( 1 + c 1 ) σ 1 2 ( 1 q ) Λ γ + μ G + 1 2 c 2 ( μ 2 + Z 2 ) + σ 1 2 [ ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2 + M 2 ] + ( 1 + c 1 ) σ 1 2 q 2 Λ 2 γ 2 + ( 1 + c 1 ) σ 1 2 ( 1 q ) 2 Λ 2 ( γ + μ ) 2 + σ 1 2 ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2

这里令

2 c 1 β q Λ γ + c 2 β = 0 4 γ + c 2 β σ = 0

那么有

c 1 = 2 γ 2 β q Λ σ > 0 c 2 = 4 γ β σ > 0

所以得到

L V [ 2 γ + ( 1 + c 1 ) σ 1 2 2 c 1 γ ] G 2 + [ 2 γ + ( 1 + c 1 ) σ 1 2 2 c 1 γ 2 c 1 μ ] H 2 + ( 2 γ + σ 1 2 ) Z 2 + ( 2 γ + σ 1 2 ) M 2 + ( 1 + c 1 ) σ 1 2 q 2 Λ 2 γ + ( 1 + c 1 ) σ 1 2 G 2 + 1 2 c 2 ( μ 2 + H 2 ) + ( 1 + c 1 ) σ 1 2 ( 1 q ) 2 Λ 2 ( γ + μ ) 2 + ( 1 + c 1 ) σ 1 2 H 2 + 1 2 c 2 ( μ 2 + Z 2 ) + σ 1 2 [ ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2 + M 2 ] + ( 1 + c 1 ) σ 1 2 q 2 Λ 2 γ 2 + ( 1 + c 1 ) σ 1 2 ( 1 q ) 2 Λ 2 ( γ + μ ) 2 + σ 1 2 ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2

2 ( 1 + c 1 ) σ 1 2 q 2 Λ 2 γ 2 + c 2 μ 2 + 2 ( 1 + c 1 ) σ 1 2 ( 1 q ) 2 Λ 2 ( γ + μ ) 2 + 2 σ 1 2 ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2 + [ 2 γ + 2 ( 1 + c 1 ) σ 1 2 2 c 1 γ ] G 2 + [ 2 γ + 2 ( 1 + c 1 ) σ 1 2 2 c 1 ( γ + μ ) + 2 c 2 ] H 2 + ( 2 γ + 2 σ 1 2 + 2 c 2 ) Z 2 + ( 2 γ + 2 σ 1 2 ) M 2

A = 2 ( 1 + c 1 ) σ 1 2 q 2 Λ 2 γ 2 + c 2 μ 2 + 2 ( 1 + c 1 ) σ 1 2 ( 1 q ) 2 Λ 2 ( γ + μ ) 2 + 2 σ 1 2 ( 1 q ) 2 μ 2 Λ 2 ( γ + μ ) 2 γ 2

0 E V ( G ( 0 ) , H ( 0 ) , Z ( 0 ) , M ( 0 ) ) + E 0 t { A + [ 2 γ + 2 ( 1 c 1 ) 2 c 1 γ ] G 2 + [ 2 γ + 2 ( 1 c 1 ) 2 c 1 ( γ + μ ) + 2 c 2 ] H 2 + [ 2 γ + 2 σ 1 2 + 2 c 2 ] Z 2 + [ 2 γ + 2 σ 1 2 ] M 2 } d s

这里取

B = min { 2 γ 2 ( 1 + c 1 ) σ 1 2 + 2 c 1 γ , 2 γ 2 ( 1 + c 1 ) σ 1 2 + 2 c 1 ( γ + μ ) 2 c 2 , 2 γ 2 σ 1 2 2 c 2 , 2 γ 2 σ 1 2 }

因此

0 E V V ( G ( 0 ) , H ( 0 ) , Z ( 0 ) , M ( 0 ) ) E B 0 t ( G 2 + H 2 + Z 2 + M 2 ) d s + A t

进而

E B 0 t ( G 2 + H 2 + Z 2 + M 2 ) d s V [ G ( 0 ) , H ( 0 ) , Z ( 0 ) , M ( 0 ) ] + A t

所以

1 t E 0 t ( G 2 + H 2 + Z 2 + M 2 ) d s V [ G ( 0 ) , H ( 0 ) , Z ( 0 ) , M ( 0 ) ] B t + A B

那么

lim t sup 1 t E 0 t ( G 2 + H 2 + Z 2 + M 2 ) d s A B

定理得证。

4. 在正平衡点附近的渐近性态

文献 [10] 已证明了在确定性方程中的正平衡点 E 1 = ( S u * , S e * , I * , R * ) 的存在性,并且给出了阈值,证明了高于该阈值的时候,系统(1)的正平衡点是全局渐近稳定的,那么在随机微分方程中,此正平衡点将不满足随机微分方程,那在这点附近是否仍然满足条件使之是渐近稳定的呢?下面证明在一定条件下,系统(2)的解在点 E 1 = ( S u * , S e * , I * , R * ) 附近的渐近性。

定理3 若 0 = β q Λ γ ( γ + μ ) + β σ Λ ( 1 q ) ( γ + μ ) 2 > 1 且满足条件:

1) 4 γ b β 2 > 0

2) 2 γ σ 1 2 > 0

3) 2 γ + d γ σ 1 2 d σ 1 2 2 > 0

则对于任何给定的初值 ( S u ( 0 ) , S e ( 0 ) , I ( 0 ) , R ( 0 ) ) R + 4 ,随机微分方程满足如下性质:

lim t sup 1 t E 0 t [ ( S u k 1 2 S u ) 2 + ( S e k 2 S e ) 2 + ( I k 3 I ) 2 + ( R k 4 R ) 2 ] d s N k 5

其中:

k 1 = 8 γ S u + b β S u 4 γ S e * 4 γ R ( 4 γ + b β 2 + σ 1 2 ) S u k 2 = 2 γ 2 γ S u S e 2 γ + σ 1 2 k 3 = 4 γ + b β 2 4 γ + b β 2 + σ 1 2

k 4 = 2 γ d γ 2 γ S u R 2 γ d γ + σ 1 2 + d σ 1 2 2 k 5 = min { 4 γ b β 2 , 2 γ σ 1 2 , 2 γ + d γ σ 1 2 d σ 1 2 2 }

N = k 1 2 4 ( 4 γ b β 2 σ 1 2 ) S u 2 + ( 2 γ σ 1 2 ) k 2 2 S e 2 + ( 4 γ b β 2 σ 1 2 ) k 3 2 I 2 + 1 2 σ 1 2 b I + ( 2 γ + d γ σ 1 2 d σ 1 2 2 ) k 4 2 R 2

证明 设

υ 1 = ( S u S u + S e S e + I I + R R ) 2 , υ 2 = I I I ln I I , υ 3 = 1 2 ( R R ) 2

构造李雅普诺夫函数

V ( S u ( t ) , S e ( t ) , I ( t ) , R ( t ) ) = υ 1 + b υ 2 + d υ 3

由伊藤公式可知:

L υ 1 = 2 ( S u S u + S e S e + I I + R R ) ( Λ γ S u γ S e γ I γ R ) + σ 1 2 S u 2 + σ 1 2 S e 2 + σ 1 2 I 2 + σ 1 2 R 2 = 2 γ ( S u S u ) 2 2 γ ( S e S e ) 2 2 γ ( I I ) 2 2 γ ( R R ) 2 4 γ ( S u S u ) ( S e S e ) 4 γ ( S u S u ) ( I I ) 4 γ ( S u S u ) ( R R ) 4 γ ( S e S e ) ( I I ) 4 γ ( S e S e ) ( R R ) 4 γ ( I I ) ( R R ) + σ 1 2 S u 2 + σ 1 2 S e 2 + σ 1 2 I 2 + σ 1 2 R 2 L υ 2 = ( 1 I I ) ( β S u I + β σ S e I γ I μ I ) + 1 2 σ 1 2 I = ( I I ) ( β S u + β σ S e γ μ ) + 1 2 σ 1 2 I = β ( S u S u ) ( I I ) + β σ ( S e S e ) ( I I ) + 1 2 σ 1 2 I

L υ 3 = ( R R ) ( μ S e + μ I γ R ) + 1 2 σ 1 2 R 2 = μ ( S e S e ) ( R R ) + μ ( I I ) ( R R ) γ ( R R ) 2 + 1 2 σ 1 2 R 2

那么:

L V = L υ 1 + b L υ 2 + d L υ 3 = 2 γ ( S u S u ) 2 2 γ ( S e S e ) 2 2 γ ( I I ) 2 2 γ ( R R ) 2 4 γ ( S u S u ) ( S e S e ) 4 γ ( S u S u ) ( I I ) 4 γ ( S u S u ) ( R R ) 4 γ ( S e S e ) ( I I ) 4 γ ( S e S e ) ( R R ) 4 γ ( I I ) ( R R ) + σ 1 2 S u 2 + σ 1 2 S e 2 + σ 1 2 I 2 + σ 1 2 R 2 + b β ( S u S u ) ( I I ) + b β σ ( S e S e ) ( I I ) + 1 2 σ 1 2 b I + d μ ( S e S e ) ( R R ) + d μ ( I I ) ( R R ) d γ ( R R ) 2 + 1 2 σ 1 2 d R 2

由于 σ 表示媒体宣传程度的有效性,我们有 0 σ 1 ,那么可以得到

4 γ + b β σ 4 γ + b β

4 γ + b β σ = 0 , 4 γ + d μ = 0

得到

d = 4 γ μ , b = 4 γ β σ

那么有

4 γ + b β > 0 , ( 4 γ + b β ) ( S u S u ) ( I I ) 1 2 ( 4 γ + b β ) [ ( S u S u ) 2 + ( I I ) 2 ]

L V 2 γ ( S u S u ) 2 2 γ ( S e S e ) 2 2 γ ( I I ) 2 ( 2 γ + d ) γ ( R R ) 2 + 4 γ S u S e + 4 γ S u S e + 4 γ + b β 2 ( S u S u ) 2 + 4 γ + b β 2 ( I I ) 2 + 4 γ S u R + 4 γ S u R + σ 1 2 S u 2 + σ 1 2 S e 2 + σ 1 2 I 2 + σ 1 2 R 2 + 1 2 σ 1 2 b I + 1 2 σ 1 2 d R 2

整理后得到

L V ( 4 γ b β 2 σ 1 2 ) ( S u k 1 2 S u ) 2 + ( 4 γ + b β 2 ) S u 2 + k 1 2 4 ( 4 γ b β 2 σ 1 2 ) S u 2 ( 2 γ σ 1 2 ) ( S e k 2 S e ) 2 + ( 2 γ σ 1 2 ) k 2 2 S e 2 2 γ S e 2 ( 4 γ b β 2 σ 1 2 ) ( I k 3 I ) 2 + ( 4 γ b β 2 σ 1 2 ) k 3 2 I 2 + ( 4 γ + b β 2 ) I 2 + 1 2 σ 1 2 b I ( 2 γ + d γ σ 1 2 d σ 1 2 2 ) ( R k 4 R ) 2 + ( 2 γ + d γ σ 1 2 d σ 1 2 2 ) k 4 2 R 2 ( 2 γ + d γ ) R 2

这里取

k 1 = 8 γ S u + b β S u 4 γ S e 4 γ R ( 4 γ + b β 2 + σ 1 2 ) S u , k 2 = 2 γ 2 γ S u S e 2 γ + σ 1 2 , k 3 = 4 γ + b β 2 4 γ + b β 2 + σ 1 2 , k 4 = 2 γ d γ 2 γ S u R 2 γ d γ + σ 1 2 + d σ 1 2 2

由于

4 γ b β 2 4 γ b β 2 σ 1 2 > 0

所以

( 4 γ b β 2 ) < 0

那么就有

L V ( 4 γ b β 2 σ 1 2 ) ( S u k 1 2 S u ) 2 + k 1 2 4 ( 4 γ b β 2 σ 1 2 ) S u 2 ( 2 γ σ 1 2 ) ( S e k 2 S e ) 2 + ( 2 γ σ 1 2 ) k 2 2 S e 2 ( 4 γ b β 2 σ 1 2 ) ( I k 3 I ) 2 + ( 4 γ b β 2 σ 1 2 ) k 3 2 I 2 + 1 2 σ 1 2 b I ( 2 γ + d γ σ 1 2 d σ 1 2 2 ) ( R k 4 R ) 2 + ( 2 γ + d γ σ 1 2 d σ 1 2 2 ) k 4 2 R 2

N = k 1 2 4 ( 4 γ b β 2 σ 1 2 ) S u 2 + ( 2 γ σ 1 2 ) k 2 2 S e 2 + ( 4 γ b β 2 σ 1 2 ) k 3 2 I 2 + 1 2 σ 1 2 b I + ( 2 γ + d γ σ 1 2 d σ 1 2 2 ) k 4 2 R 2

k 5 = min { 4 γ b β 2 , 2 γ σ 1 2 , 2 γ + d γ σ 1 2 d σ 1 2 2 }

得到

L V k 5 [ ( S u k 1 2 S u ) 2 + ( S e k 2 S e ) 2 + ( I k 3 I ) 2 + ( R k 4 R ) 2 ] + N

进而可以有

0 E V ( S u , S e , I , R ) = V 0 + E 0 t L V ( S u ( s ) , S e ( s ) , I ( s ) , R ( s ) ) d s V 0 k 5 E 0 t [ ( S u k 1 2 S u ) 2 + ( S e k 2 S e ) 2 + ( I k 3 I ) 2 + ( R k 4 R ) 2 ] d s + N t

整理后得

E 0 t [ ( S u k 1 2 S u ) 2 + ( S e k 2 S e ) 2 + ( I k 3 I ) 2 + ( R k 4 R ) 2 ] d s V 0 + N t k 5

因此

lim t sup 1 t E 0 t [ ( S u k 1 2 S u ) 2 + ( S e k 2 S e ) 2 + ( I k 3 I ) 2 + ( R k 4 R ) 2 ] d s N k 5

定理得证。

5. 小结

本文介绍的是官方媒体影响下的随机谣言传播模型。在受官方媒体影响下的确定性谣言传播模型中,存在影响因子 γ ,受时间推移、社会认知等因素而使得模型中的四类人以 γ 比例的数量离开系统,其中文中仍然将影响因子 γ 作为固定的数值,但实际上影响因子很容易受到环境的改变而发生变化,因此用某个常数值来代替影响因子将不再适用于现实环境,于是本章建立新的随机微分方程,将影响因子考虑环境噪声的影响,即 γ γ + σ B ˙ ( t ) ,然后分析其动力学行为。由于确定性模型的平衡点不在满足随机微分方程,但这类平衡点在系统(2)中的性质仍值得研究。最后发现,当满足三点要求时,即1) 0 < 1 ,2) ( 1 + c 1 ) σ 1 2 < c 1 ( γ + μ ) + γ ;3) 2 σ 1 2 + 1 2 < 2 γ 其中 c 1 = 2 γ 2 β q Λ σ > 0 时,系统(2)的解在无谣言传播点附近是渐近稳定的。相比与确定性模型得到的条件时,即当 σ 1 0 ,随机模型还需要满足(2)、(3)两点,这两点是对白噪声的限制,因此控制白噪声 σ 1 的大小能够限制谣言传播的水平,甚至消失。当 σ 1 = 0 时,此时没有白噪声的影响,与确定性模型结论保持一致。在讨论正平衡点附近的渐近性态时,本章发现,当满足以下四点条件,即1) 0 > 1 ;2) 4 γ b β 2 > 0 ;3) 2 γ σ 1 2 > 0 ;4) 2 γ + d γ σ 1 2 d σ 1 2 2 > 0 时,系

统(2)的解在正平衡点是渐近稳定的,这说明相比于确定性模型的结论,随机模型增加了白噪声的影响,就需要满足另外的三个条件使得谣言传播在一个稳定的值附近波动。若没有白噪声的影响,就没有对 σ 1 的限制,则与确定性模型的结论保持一致。

参考文献

[1] Daley, D.J. and Kendall, D.G. (1965) Stochastic Rumours. IMA Journal of Applied Mathematics, 1, 42-55.
https://doi.org/10.1093/imamat/1.1.42
[2] Makidp, T. (1973) Mathematical Models and Applications. Prentice-Hall, New Jersey.
[3] Zanette, D.H. (2001) Critical Behavior of Propagation on Small-World Networks. Physical review E, 64, 050901.
https://doi.org/10.1103/PhysRevE.64.050901
[4] Zhao, L.J., Wang, J.J., Chen, Y.C., et al. (2012) SIHR Rumor Spresding Model in Social Networks. Physica A: Statistical Mechanics and Its Applications, 391, 2444-2453.
https://doi.org/10.1016/j.physa.2011.12.008
[5] Zhu, L. and Wang, Y.G. (2017) Rumor Spreading Model with Noise Interference in Complex Social Networks. Physica A: Statistical Mechanics and Its Applications, 469, 750-760.
https://doi.org/10.1016/j.physa.2016.11.119
[6] Liu, Q.M., Li, T. and Sun, M.C. (2017) The Analysis of an SEIR Rumor Propagation Model on Heterogeneous Network. Physica A: Statistical Mechanics and Its Applications, 469, 372-380.
https://doi.org/10.1016/j.physa.2016.11.067
[7] 陈安, 王子君, 陈樱花. 基于SEIR模型视角的重大公共卫生事件中伪科学网络谣言的传播治理: 以新冠肺炎疫情为例[J]. 科学导报, 2020, 38(4): 55-65.
[8] 宋雪, 武力兵. 一类带有媒体效应的谣言传播模型稳定性分析[J]. 辽宁科技大学学报, 2020, 43(4): 315-320.
[9] 陈华. 一类具有媒体播报效应的谣言传播模型的定性分析[J]. 先科技大学学报, 2016, 36(2): 255-264.
[10] 赵敏, 陈文霞. 基于微分方程的谣言传播模型建立与分析[J]. 四川师范大学学报, 2019, 42(5): 626-632.