一類(lèi)具有媒體播報(bào)效應(yīng)的謠言傳播模型的定性分析*

陳　華(1．西安科技大學(xué)能源學(xué)院，陜西西安710054; 2．延安大學(xué)校辦，陜西延安716000)

?

一類(lèi)具有媒體播報(bào)效應(yīng)的謠言傳播模型的定性分析*

陳華1，2
(1．西安科技大學(xué)能源學(xué)院，陜西西安710054; 2．延安大學(xué)校辦，陜西延安716000)

摘要:為了準(zhǔn)確刻畫(huà)謠言傳播機(jī)理，抑制謠言傳播，文中研究了一類(lèi)具有媒體播報(bào)效應(yīng)的謠言傳播模型。利用常微分方程穩(wěn)定性理論，詳細(xì)分析了系統(tǒng)平衡點(diǎn)的存在性及局部漸近穩(wěn)定性;通過(guò)比較定理證明了邊界平衡點(diǎn)E0(rumor-free equilibrium)的全局漸近穩(wěn)定性，利用Lyapunov-La-Salle不變性原理給出了正平衡點(diǎn)E+(rumor-endemic equilibrium)全局漸近穩(wěn)定的充分條件。結(jié)論表明，隨著時(shí)間的推移，要么傳謠者的數(shù)量趨于0，要么謠言持續(xù)，媒體報(bào)道對(duì)謠言傳播的規(guī)模有一定影響。最后，通過(guò)數(shù)值模擬驗(yàn)證了媒體播報(bào)對(duì)謠言傳播和控制的影響。

關(guān)鍵詞:謠言傳播;穩(wěn)定性; Lyapunov函數(shù);媒體播報(bào)

0　引言

突發(fā)公共事件應(yīng)急管理是公共安全保障的核心問(wèn)題，包含著豐富而深刻的復(fù)雜性科學(xué)問(wèn)題。突發(fā)公共事件應(yīng)急管理的研究已經(jīng)成為不同領(lǐng)域?qū)＜?、學(xué)者共同關(guān)注的研究課題［1－5］。突發(fā)公共事件會(huì)引發(fā)高度的社會(huì)關(guān)注，在信息傳播過(guò)程中容易出現(xiàn)消息失真、輿論失實(shí)的情況，加之部分民眾的主觀猜測(cè)與惡意關(guān)聯(lián)，致使謠言比平時(shí)更容易滋生、感染和傳播。謠言是作為一種普遍的社會(huì)輿論現(xiàn)象，具有信息性、傳播性和未知性等特點(diǎn)。在突發(fā)事件乃至各種危機(jī)中，一些關(guān)于國(guó)家發(fā)展、社會(huì)生活及個(gè)人利益的敏感信息、丑聞信息或謠言信息等，在經(jīng)過(guò)社交網(wǎng)絡(luò)媒介的發(fā)酵之后，可迅速發(fā)展為網(wǎng)絡(luò)突發(fā)事件，引起巨大的負(fù)面效應(yīng)，造成社會(huì)心理恐慌和非理性反應(yīng)，其后果嚴(yán)重影響民眾生活與社會(huì)穩(wěn)定［6－8］。例如，2011年，日本核泄漏引發(fā)的搶購(gòu)食鹽風(fēng)波極大地影響了正常的社會(huì)秩序;“六翅肯德基怪雞”、“康師傅地溝油”、“娃哈哈肉毒桿菌”等謠言信息引起的網(wǎng)絡(luò)突發(fā)事件不僅給涉事品牌造成了巨大傷害，也引起了各方的極度恐慌。因此，如何找到行之有效的方法，分析謠言傳播的機(jī)理，抑制謠言的傳播擴(kuò)散，成為了社會(huì)科學(xué)和自然科學(xué)不同領(lǐng)域共同關(guān)注的課題。

由于謠言在社交網(wǎng)絡(luò)中的散布和病毒傳播、擴(kuò)散很相似，Daley和Kendall［4］借鑒傳染病“倉(cāng)室”模型的概念提出了謠言傳播的數(shù)學(xué)模型(簡(jiǎn)稱DK模型［9］)，指出雖然表面上看傳染病和謠言的傳播是相似的，但是二者在傳播機(jī)理上的細(xì)微差異卻導(dǎo)致傳播結(jié)果的明顯不同。隨后，學(xué)者們對(duì)D －K模型進(jìn)行了各種不同的推廣，現(xiàn)有的謠言傳播模型大都借鑒了傳染病模型。Dietz［10］在總結(jié)傳染病傳播和謠言擴(kuò)散機(jī)理之間異同的基礎(chǔ)上，重點(diǎn)闡述了利用傳染病模型解釋謠言傳播的合理性。Maki和Thompson［11］對(duì)DK模型的傳播機(jī)制進(jìn)行了修改，認(rèn)為謠言是通過(guò)傳播者與人群中其他人的雙向接觸進(jìn)行傳播的，并由此得到了新的著名的Maki-Thompson謠言傳播模型。Kawachi［12］等考慮各種不同的接觸方式，建立了謠言傳播的SIR模型，但是由于參數(shù)過(guò)多導(dǎo)致難以對(duì)數(shù)學(xué)模型進(jìn)行完整的定性分析，通過(guò)數(shù)學(xué)分析與數(shù)值模擬相結(jié)合確定了謠言是否傳播的閾值?；袅及驳龋?3］基于D－K模型的建?？蚣埽紤]個(gè)體心理因素對(duì)于不實(shí)信息的傳播影響，采用非線性函數(shù)刻畫(huà)不實(shí)信息傳播力度的變化趨勢(shì)，通過(guò)比較與數(shù)據(jù)的模擬說(shuō)明了非線性傳播率的合理性，探討了不實(shí)信息控制的閾值。趙來(lái)軍等［14］在考慮遺忘率的基礎(chǔ)上，在D－K模型基礎(chǔ)上增加一類(lèi)Hibernators人群(表示由于遺忘機(jī)制不能立即傳播謠言，而又由于記憶機(jī)制后來(lái)又傳播謠言的人群)，根據(jù)對(duì)遺忘機(jī)制和記憶機(jī)制的研究，提出SIHR謠言傳播模型，研究表明由于遺忘機(jī)制和記憶機(jī)制的存在推遲了謠言傳播終止的時(shí)間，降低了謠言的最大影響力。上述研究表明基于傳染病模型和D－K模型的謠言傳播模型是描述、研究謠言傳播規(guī)律的有效手段。

媒體報(bào)道在傳播理論中最早被稱為:“皮下注射”、“魔幻子彈”理論，媒體報(bào)道直接影響個(gè)體對(duì)于事物的認(rèn)識(shí)，對(duì)主體價(jià)值觀的形成具有重要影響［15］。研究表明Ma在研究中表明媒體報(bào)道對(duì)于民眾理解能力、應(yīng)急管理都存在巨大的影響力［16］。在大型傳染病爆發(fā)時(shí)，媒體報(bào)道對(duì)疾病傳播的影響是復(fù)雜的，會(huì)影響疾病的最終傳播狀態(tài)，適度報(bào)道可以限制傳染病的傳播，減少傳染病傳播帶來(lái)的損失，但是過(guò)度報(bào)道會(huì)引起恐慌，也會(huì)帶來(lái)額外的災(zāi)難后果［17］。劉玉英和肖燕妮［18］指出媒體報(bào)道的作用會(huì)極大的減小感染率，通常采用媒體影響因子函數(shù)來(lái)刻畫(huà)媒體報(bào)道對(duì)感染率的影響。文獻(xiàn)［16］～［24］及相關(guān)參考文獻(xiàn)討論了不同媒體影響因子函數(shù)對(duì)傳染病模型的影響，結(jié)論表明媒體影響并不是決定疾病爆發(fā)與否的最重要因素，但媒體對(duì)疾病傳播的規(guī)模有很大的影響。上述文獻(xiàn)的研究焦點(diǎn)主要集中于媒體播報(bào)對(duì)傳染率的影響，并沒(méi)有對(duì)媒體播報(bào)給予量化處理。考慮到媒體播報(bào)數(shù)量，Misra等［25］建立了一類(lèi)媒體播報(bào)被量化的傳染病模型，潘艷雪等［26］建立了一類(lèi)具有時(shí)滯和媒體播報(bào)雙重效應(yīng)的流感模型，研究表明媒體播報(bào)數(shù)量嚴(yán)重影響傳染者的數(shù)量，適當(dāng)?shù)拿襟w播報(bào)對(duì)防控具有重要意義。事實(shí)上，謠言傳播跟傳染病傳播有類(lèi)似之處，也會(huì)受到媒體報(bào)道的影響，媒體的報(bào)道使民眾認(rèn)識(shí)了解不實(shí)信息，辨識(shí)虛假信息，不信謠、不傳謠。霍良安和黃培清［15］基于系統(tǒng)動(dòng)力學(xué)的思想，提出了不實(shí)信息的動(dòng)態(tài)傳播模型，刻畫(huà)了科普教育以及媒體報(bào)道對(duì)于不實(shí)信息傳播的影響。趙洪涌等［7］結(jié)合時(shí)空滯后、空間擴(kuò)散、媒體報(bào)道等因素對(duì)社交網(wǎng)絡(luò)中謠言傳播的影響，建立偏微分方程描述的社交網(wǎng)絡(luò)謠言時(shí)空傳播模型，通過(guò)引入概率函數(shù)來(lái)描述隨著媒體報(bào)道的力度增大，真正進(jìn)入網(wǎng)絡(luò)中的謠言無(wú)知者密度會(huì)逐漸減少這一特性，研究表明媒體對(duì)社交網(wǎng)絡(luò)謠言的深度報(bào)道可以極大地減小網(wǎng)絡(luò)中謠言傳播者密度，并擴(kuò)大系統(tǒng)的穩(wěn)定區(qū)域。由于信息傳播與疾病傳播的本質(zhì)差異，在謠言傳播過(guò)程中，媒體播報(bào)除了會(huì)降低未知者向傳謠者轉(zhuǎn)移的幾率，未知者還會(huì)有幾率直接成為事實(shí)真相的知曉者從而成為免疫者。然而，現(xiàn)有文獻(xiàn)并未考慮這種差異，文中的目的是在理論上對(duì)具有媒體播報(bào)效應(yīng)的謠言傳播模型進(jìn)行動(dòng)力學(xué)分析，并由此來(lái)研究媒體報(bào)道對(duì)謠言傳播的影響。

1　模型的建立

文中基于D－K模型的框架，建立具有媒體播報(bào)效應(yīng)的謠言傳播動(dòng)力學(xué)模型。首先，根據(jù)社交網(wǎng)絡(luò)中謠言傳播的現(xiàn)實(shí)背景和傳播機(jī)制，作如下幾點(diǎn)說(shuō)明

1)根據(jù)倉(cāng)室模型方法，主要將社交網(wǎng)絡(luò)中的人群分成3大類(lèi)，S(t) :表示易染類(lèi)個(gè)體，個(gè)體從未聽(tīng)過(guò)謠言，也稱為易染者(或未知者，Ignorants，類(lèi)似于Susceptile individuals) ; I(t) :表示感染類(lèi)個(gè)體，個(gè)體知道謠言并傳播謠言，也稱為傳播者(Spreaders，類(lèi)似于Infective individuals) ; R(t) :表示移出類(lèi)個(gè)體，個(gè)體知道謠言但不傳播謠言也稱為免疫者(或緘默者，Stiflers，Removed individuals)。

2)利用傳染病模型中研究媒體報(bào)道因素影響的方法，考慮媒體報(bào)道影響下的謠言傳播模型，M (t)表示媒體播報(bào)的數(shù)量。媒體播報(bào)數(shù)量(M)與傳播者(I)數(shù)量相關(guān)，記γ為播報(bào)率，即由于傳播者(I)數(shù)量而引起的媒體播報(bào)數(shù)量的變化率。隨著事件的發(fā)展媒體欄目對(duì)某一新聞事件報(bào)道的次數(shù)會(huì)逐漸減少，假設(shè)減少率為μ0．

3)假設(shè)個(gè)體以恒定的比率(B)進(jìn)入特定的社交網(wǎng)絡(luò);每一類(lèi)個(gè)體由于某種原因以一定的概率(μ)移出該類(lèi)人群。

4)基于D－K模型對(duì)于傳播規(guī)律的描述，進(jìn)一步假設(shè):傳播者(I)與未知者(S)相遇時(shí)以一定的概率β傳播謠言;傳播者(I)與傳播者(I)或傳播者(I)與免疫者(R)相遇時(shí)，傳播者(I)因失去繼續(xù)傳播的興趣以概率α變?yōu)槊庖哒?R) ;受媒體報(bào)道的影響，未知者(S)和傳播者(I)會(huì)以一定的概率變?yōu)槊庖哒?R)，變化值正比于未知者(S) (或傳播者(I) )與播報(bào)量(M)的乘積，即λSM，ηIM．其中λ和η分別表示受媒體播報(bào)影響，未知者(S)與傳播者(I)向免疫者(R)轉(zhuǎn)變的概率。(λ，η分別表示什么)如圖1所示。

圖1　媒體報(bào)道影響下謠言傳播示意圖Fig．1 Transfer diagram of the rumor transmission model with media coverage

基于以上分析，建立如下謠言傳播模型

假設(shè)各參數(shù)都是正數(shù)，初始條件為

2　模型分析

本節(jié)將對(duì)系統(tǒng)(1)進(jìn)行定性分析，關(guān)于解的存在性、非負(fù)性、有界性有如下定理。

定理1系統(tǒng)(1)滿足初始條件(2)的解(S (t)，I(t)，R(t)，M(t) )在［0，+∞)上存在、非負(fù)、有界并且關(guān)于Ω是正向不變的，其中

證明容易證明系統(tǒng)(1)滿足初始條件(2)的解(S(t)，I(t)，R(t)，M(t) )在［0，+∞)上存在、非負(fù)。將系統(tǒng)(1)的前3個(gè)方程相加，可以得到如下的方程

(S(t) + I(t) + R(t) ) ' = B－μ(S(t) + I(t) + R(t) )，

且有S(0) +I(0) +R(0)≥0．顯然，對(duì)所有的t≥0有

由此，可以得到系統(tǒng)(1)的耗散性和

又由系統(tǒng)(1)的第4個(gè)方程，

從而

從而，系統(tǒng)(1)滿足初始條件(2)的解(S(t)，I (t)，R(t)，M(t) )在［0，+∞)上最終有界。

下面證明系統(tǒng)(1)關(guān)于Ω是正向不變的。

對(duì)于任意解(S(t)，I(t)，R(t)，M(t) )∈Ω， I(t) = 0是系統(tǒng)(1)的軌線且I(0) = 0，因此不可能有如任何一個(gè)解在有限時(shí)間內(nèi)與I(t) =0相交，所以在第一卦限內(nèi)的解必保持在第一卦限內(nèi)。因此，系統(tǒng)(2)關(guān)于Ω是正向不變的。

下面研究系統(tǒng)(1)的平衡點(diǎn)的存在性和穩(wěn)定性。關(guān)于平衡點(diǎn)的存在性，有如下定理。

定理2 1)系統(tǒng)(1)總存在邊界平衡點(diǎn)(Ru-

mor-free equilibrium，RFE) E0=;

2)如果基本再生數(shù)R0=＞1，系統(tǒng)(1)存在唯一正平衡點(diǎn)(Rumor-endemic equilibrium，REE) E+=，I*，－I*，)，其中I*是方程βB =0的唯一正根。

證明令系統(tǒng)(1)方程的右邊等于零，即

．

顯然，邊界平衡點(diǎn)E0=總是存在的。

當(dāng)S≠0，I≠0，R≠0，M≠0時(shí)，得

其中

由文獻(xiàn)［22］，定義基本再生數(shù)為

由方程組(3)的第4個(gè)方程得M*=．將I*， M*代入方程組(3)的第一個(gè)方程得

關(guān)于平衡點(diǎn)的穩(wěn)定性，有如下定理。

定理3當(dāng)R0＜1時(shí)，E0是局部漸近穩(wěn)定的; 當(dāng)R0＞1時(shí)，E0是不穩(wěn)定的。

證明系統(tǒng)(1)在邊界平衡點(diǎn)E0的雅可比矩陣為

矩陣JE0的特征值為

由于系統(tǒng)(1)中所有系數(shù)均非負(fù)，因此，ξ1，2＜0，ξ3＜0，當(dāng)μ2－βB＞0，即R0＜1時(shí)，ξ4＜0．由Routh-Hurwitz穩(wěn)定判據(jù)［27］可知，當(dāng)R0＜1時(shí)，E0是局部漸近穩(wěn)定的;當(dāng)R0＞1時(shí)，E0是不穩(wěn)定的。

定理4當(dāng)R0＞1且珚ω＞0時(shí)，E+是局部漸近穩(wěn)定的;當(dāng)R0＞1且珚ω＜0時(shí)，E+是不穩(wěn)定的，其中

證明系統(tǒng)(1)在正平衡點(diǎn)E+的雅可比矩陣為

相應(yīng)的特征方程為

其中

顯然，特征方程(5)的一個(gè)特征根為x1=－μ＜0，從而E+的穩(wěn)定性取決于方程

的根的符號(hào)。

由于

而a3=γλ(α+β) S*I*+μ0(α+β)βS*I*+γηI*(μ+λM*+βI*)＞0，從而H3的符號(hào)與H2相同。

由Routh-Hurwitz穩(wěn)定判據(jù)［27］可知，當(dāng)＞0時(shí)，E+是局部漸近穩(wěn)定的;當(dāng)＜0時(shí)，E+是不穩(wěn)定的。

定理5當(dāng)R0＜1時(shí)，系統(tǒng)(1)的邊界平衡點(diǎn)E0關(guān)于Ω是全局漸近穩(wěn)定的。

證明下面利用比較定理［28－30］證明邊界平衡點(diǎn)E0關(guān)于Ω是全局漸近穩(wěn)定的。由于參數(shù)α，β，η，μ均為正數(shù)且在Ω內(nèi)系統(tǒng)(1)滿足初始條件的任一解(S(t)，I(t)，R(t)，M(t) )都是非負(fù)的。因此，由系統(tǒng)(1)的第2個(gè)方程可以得到

當(dāng)R0＜1時(shí)，≤0．

因此，由比較定理［28－30］，當(dāng)t→∞，I(t)→0．再由系統(tǒng)(1)可得當(dāng)I(t) = 0時(shí)，S(t)→，R(t)→0，M (t)→0．

因此，當(dāng)R0＜1時(shí)，系統(tǒng)(1)的邊界平衡點(diǎn)E0關(guān)于Ω是全局漸近穩(wěn)定的。

定理6當(dāng)R0＞1，珚ω＞0且Υ≥0時(shí)，系統(tǒng)(1)的正平衡點(diǎn)E+是全局漸近穩(wěn)定的。

由定理1，系統(tǒng)(7)所有滿足初始條件的解的ω極限集均在Z = 0上，即t．在解平面S + I + R =上，系統(tǒng)(7)的解滿足系統(tǒng)(1)的正平衡點(diǎn)E+(S*，I*，R*，M*)對(duì)應(yīng)于系統(tǒng)(8)的正平衡點(diǎn)E+(I*，R*，M*)其中，S*=Bμ －I*－R*．

由定理4，當(dāng)R0＞1＞0時(shí)，E+是局部漸近穩(wěn)定的。下面只需證明系統(tǒng)(8)的正平衡點(diǎn)是全局漸近穩(wěn)定的。下面利用Lyapunov-LaSalle不變性原理［24］證明E+關(guān)于Ω是全局吸引的。

考慮緊集Ω上的非負(fù)函數(shù)

顯然V1在Ω上連續(xù)，并且沿著系統(tǒng)(8)的解的導(dǎo)數(shù)滿足

考慮緊集Ω上的非負(fù)函數(shù)

顯然V2在Ω上連續(xù)，并且沿著系統(tǒng)(8)的解的導(dǎo)數(shù)滿足

考慮緊集Ω上的非負(fù)函數(shù)

其中Λ為待定常數(shù)。顯然V3在Ω上連續(xù)，并且沿著系統(tǒng)(8)的解的導(dǎo)數(shù)滿足

定義W = V1+ V2+ V3．顯然W在Ω上連續(xù)，并且沿著系統(tǒng)(8)的解的導(dǎo)數(shù)滿足

由定理1，I(t)，M(t)，R(t)均是有界的，從而可記

則

因此，W是系統(tǒng)(8)在Ω上的一個(gè)Lyapunov函數(shù)。定義Ω的子集E'為

E' = { (I，R，M) |(I，R，M)∈Ω，W' =0}．同時(shí)，令M'為系統(tǒng)(8)在E'上的最大不變子集。

因此，

E' = { (I，R，M) | (I，R，M)∈Ω，I = I*，R = R*，M = M*}．

由M'的不變性和系統(tǒng)(1)與(8)，容易證明E' = E+．

所以，由Lyapunov-LaSalle不變性原理得到E+是全局吸引的。又由定理4，當(dāng)R0＞1，珚ω＞0時(shí)，E+是局部漸近穩(wěn)定的。因此，當(dāng)R0＞1，珚ω＞0且Υ ≥0時(shí)，E+是全局漸近穩(wěn)定的。

3　數(shù)值仿真分析

本節(jié)根據(jù)系統(tǒng)(1)的實(shí)際意義，取不同參數(shù)，利用Matlab軟件進(jìn)行數(shù)值模擬，討論媒體報(bào)道對(duì)謠言傳播的影響。

3．1邊界平衡點(diǎn)的穩(wěn)定性

為保證R0＜1，在系統(tǒng)(1)中取參數(shù)(說(shuō)明選取的依據(jù)是什么) B = 0．9，μ= 0．2，β= 0．02，α= 0．1，η=0．05，λ= 0．04，γ= 0．3，μ0= 0. 1.經(jīng)計(jì)算E0= (4．5，0，0，0)，R0=0．45＜1，滿足定理5的條件。圖2和圖3顯示了在上述參數(shù)下，邊界平衡點(diǎn)E0的全局漸近穩(wěn)定性。

為探討媒體播報(bào)對(duì)謠言傳播的影響，在其他參數(shù)不變的條件下，改變播報(bào)率γ的取值，觀察γ分別取0. 01，0. 4，0. 7時(shí)S(t)，I(t)，R(t)和M(t)的變化曲線，如圖4～圖7．從圖4可以得出，γ越小S(t)越小，且播報(bào)率越大未知者的人數(shù)越小，說(shuō)明媒體播報(bào)使更多的人了解謠言。從圖5可以得出，γ越大I(t)趨于0的速度越快，說(shuō)明媒體播報(bào)對(duì)謠言傳播的具有限制作用;從圖6可以得出，γ越大R(t)的峰值越大，說(shuō)明媒體的播報(bào)使得聽(tīng)說(shuō)過(guò)謠言的人數(shù)增多;從圖7可以得出，γ越大M(t)的峰值越大，說(shuō)明媒體播報(bào)數(shù)量的增大。

圖2　S－I相平面，顯示E0全局漸近穩(wěn)定Fig．2　S－I plane，E0is globally asymptotically stable

圖3　I－M相平面，顯示E0全局漸近穩(wěn)定Fig．3　I－M plane，E0is globally asymptotically stable

圖4　S(t)與播報(bào)率γ的關(guān)系圖Fig．4　Variation of Ignorant S(t) with time for different values of γ

圖5　I(t)與播報(bào)率γ的關(guān)系圖Fig．5　Variation of Spreader I(t) with time for different values of γ

3. 2正平衡點(diǎn)的穩(wěn)定性

為保證R0＞1，在系統(tǒng)(1)中取參數(shù)B =0. 9，μ =0. 1，β= 0. 02，α= 0. 01，η= 0. 05，λ= 0. 04，γ= 0. 05，μ0= 0. 1.經(jīng)計(jì)算E+= (6．948 5，0．738 3，1．313 2，0．369 2)，R0= 1．8＞1．圖8和圖9顯示了在上述參數(shù)下，正平衡點(diǎn)E+的全局漸近穩(wěn)定性。

圖6　R(t)與播報(bào)率γ的關(guān)系圖Fig．6　Variation of Stifler R(t) with time for different values of γ

圖7　M(t)與播報(bào)率γ的關(guān)系圖Fig．7　Variation of Media M(t) with time for different values of γ

圖8　S－I相平面，顯示E+全局漸近穩(wěn)定Fig．8　S－Iplane，E+is globally asymptotically stable

圖9　I－M相平面，顯示E+全局漸近穩(wěn)定Fig．9 I－M plane，E+is globally asymptotically stable

圖10　S(t)與播報(bào)率γ的關(guān)系圖Fig．10　Variation of Ignorant S(t) with time for different values of γ

圖11　I(t)與播報(bào)率γ的關(guān)系圖Fig．11　Variation of Spreader I(t) with time for different values of γ

圖12　R(t)與播報(bào)率γ的關(guān)系圖Fig．12　Variation of Stifler R(t) with time for different values of γ

圖13　M(t)與播報(bào)率γ的關(guān)系圖Fig．13 Variation of Media M(t) with time for different values of γ

為探討媒體播報(bào)對(duì)謠言傳播的影響，在其他參數(shù)不變的條件下，改變播報(bào)率γ的取值，觀察γ分別取0. 05，0. 4，0. 7時(shí)S(t)，I(t)，R(t)和M(t)的變化曲線，如圖10～圖13所示。從圖10可以得出，γ越小S(t)的峰值越大，趨于穩(wěn)定的值越小，說(shuō)明媒體播報(bào)使得未聽(tīng)過(guò)謠言的人減少;從圖11可以得出，γ越大I(t)的峰值越小且穩(wěn)定值也越小，說(shuō)明媒體播報(bào)對(duì)謠言傳播具有限制作用;從圖12可以得出，γ越大R(t)的穩(wěn)定值越大，說(shuō)明媒體的播報(bào)使得聽(tīng)過(guò)謠言但不傳謠的人增多，然而從圖中可以看出當(dāng)γ= 0．4時(shí)R(t)的峰值大于γ=0．7時(shí)，說(shuō)明過(guò)度播報(bào)可能使得聽(tīng)過(guò)謠言但不傳謠的人數(shù)減少，恰是適得其反;從圖13可以得出，γ越大M(t)的峰值越大，說(shuō)明媒體播報(bào)數(shù)量越大。

4　結(jié)論

1)在疾病傳播過(guò)程中，受媒體報(bào)道影響，易感人群只會(huì)降低被感染的概率，不會(huì)直接變成免疫人群。但是，在謠言傳播過(guò)程中，由于信息傳播特點(diǎn)與疾病傳播特點(diǎn)的本質(zhì)差異，媒體播報(bào)除了會(huì)降低未知者向傳謠者轉(zhuǎn)移的幾率，未知者還會(huì)有幾率直接成為事實(shí)真相的知曉者從而成為免疫者。文中基于上述考慮，建立了一類(lèi)具有媒體播報(bào)效應(yīng)影響的謠言傳播模型。利用常微分方程穩(wěn)定性分析方法，得到了平衡點(diǎn)全局漸近穩(wěn)定的充分條件;

2)數(shù)值實(shí)驗(yàn)結(jié)果表明，媒體報(bào)道可使未聽(tīng)過(guò)謠言的人數(shù)減少;當(dāng)播報(bào)率較小時(shí)，媒體報(bào)道可使聽(tīng)過(guò)謠言但不傳謠的人增多;但是，當(dāng)播報(bào)率較大時(shí)，過(guò)度播報(bào)也可能使得聽(tīng)過(guò)謠言但不傳謠的人數(shù)減少適度的媒體報(bào)道可以有效控制謠言的擴(kuò)散，但是過(guò)度的媒體報(bào)道反而會(huì)適得其反。因此，媒體報(bào)道是把雙刃劍，在謠言控制過(guò)程中應(yīng)合理規(guī)劃媒體播報(bào)次數(shù);

3)由于個(gè)體的傳播行為會(huì)受到個(gè)體獲得的信息量以及個(gè)體的偏好決定的影響，而謠言傳播過(guò)程中，媒體、政府的應(yīng)對(duì)措施等因素都會(huì)對(duì)信息的傳播與擴(kuò)散產(chǎn)生影響，文中的結(jié)論表明，政府在應(yīng)急管理中積極采取應(yīng)對(duì)措施，媒體及時(shí)公布權(quán)威信息、加大媒體覆蓋面，引導(dǎo)公眾輿論，是應(yīng)急管理中信息管理的關(guān)鍵所在。

參考文獻(xiàn)References

［1］袁曉芳，李紅霞，田水承，等．基于BioLAB的煤礦工人突發(fā)事件場(chǎng)景下應(yīng)急反應(yīng)的評(píng)估研究［J］．西安科技大學(xué)學(xué)報(bào)，2014，34(5) : 546－552. YUAN Xiao-fang，LI Hong-xia，TIAN Shui-cheng，et al． Evaluation of coal miners’emergency response capacity in emergency scenarios based on BioLAB［J］．Journal of Xi’an University of Science and Technology，2014，34 (5) : 546－552.

［2］董丁穩(wěn)，陳贊．基于貝葉斯分類(lèi)的煤礦系統(tǒng)安全動(dòng)態(tài)定量評(píng)價(jià)［J］．西安科技大學(xué)學(xué)報(bào)，2015，35(3) : 320－324，330. DONG Ding-wen，CHEN Zan．Dynamic and quantitative safety evaluation of mine system based on Bayesian classification［J］．Journal of Xi’an University of Science and Technology，2015，35(3) : 320－324，330.

［3］張儉讓，王江艷，董丁穩(wěn)，等．基于優(yōu)化定權(quán)與灰色聚類(lèi)的人流密集場(chǎng)所應(yīng)急疏散能力評(píng)價(jià)［J］．西安科技大學(xué)學(xué)報(bào)，2014，34(4) : 451－455. ZHANG Jian-rang，WANG Jiang-yan，DONG Ding-wen，et al．Evaluation of emergency evacuation capacity based on weight optimization and grey clustering in assembly occupancies［J］．Journal of Xi’an University of Science and Technology，2014，34(4) : 451－455.

［4］周星．大型教學(xué)樓火災(zāi)性能指標(biāo)的數(shù)值分析［J］．西安科技大學(xué)學(xué)報(bào)，2013，33(4) : 404－409. ZHOU Xing．Numerical analysis of fire performance index of teaching building［J］．Journal of Xi’an University of Science and Technology，2013，33(4) : 404－409.

［5］王秋紅，張娟，鄧軍，等．基于火災(zāi)風(fēng)險(xiǎn)評(píng)價(jià)理論的城鎮(zhèn)消防站選址優(yōu)化［J］．西安科技大學(xué)學(xué)報(bào)，2014，34(6) : 681－685. WANG Qiu-hong，ZHANG Juan，DENG Jun，et al．Optimization location of fire fighting stations based on fire risk assessment theory［J］．Journal of Xi’an University of Science and Technology，2014，34(6) : 681－685.

［6］張芳，司光亞，羅批．謠言傳播模型綜述［J］．復(fù)雜系統(tǒng)與復(fù)雜性科學(xué)，2009，6(4) : 1－11. ZHANG Fang，SI Guang-ya，LUO Pi．A survey for rumor propagation models［J］．Complex Systems and Complexity Science，2009，6(4) : 1－11.

［7］趙洪勇，朱霖河．社交網(wǎng)絡(luò)中謠言傳播動(dòng)力學(xué)研究［J］．南京航空航天大學(xué)學(xué)報(bào)，2015，47(3) : 332－342. ZHAO Hong-yong，ZHU Lin-he．Research on rumor spreading dynamics in social networks［J］．Journal of Nanjing University of Aeronautics＆Astronautic，2015，47(3) : 332－342.

［8］劉怡君，李倩倩，牛文元．輿論動(dòng)力學(xué)模型綜述［J］．管理評(píng)論，2013，25(1) : 167－176. LIU Yi-jun，LI Qian-qian，NIU Wen-yuan．Review on the models of opinion dynamics［J］．Management Review，2013，25(1) : 167－176.

［9］Daley D J，Kendall D G．Stochastic rumours［J］．Ima Journal of Applied Mathematics，1965，1(1) : 42－55.

［10］Dietz Klaus．Epidemics and rumours: a survey［J］．Journal of the Royal Statistical Society，1967，130(4) : 505 －528.

［11］Maki Daniel P，Thompson Maynard．Mathematical models and applications，with emphasis on social，life，and management sciences［M］．Englewood Cliffs，New Jersey: Prentice Hall，1973.

［12］Kazuki Kawachi，Motohide Seki，Hiraku Yoshida，et al．A rumor transmission model with various contact interactions［J］．Journal of Theoretical Biology，2008，253(1) : 55－60.

［13］HUO Liang-an，HUANG Pei-qing，GUO Chun-xiang．Analyzing the dynamics of a rumor transmission model with incubation［J］．Discrete Dynamics in Nature＆Society，2012，65(2) : 267－278.

［14］ZHAO Lai-jun，WANG Jia-jia，CHEN Yu-cheng，et al．SIHR rumor spreading model in social networks［J］．Physica A，2012，391(7) : 2 444－2 453.

［15］霍良安，黃培清．科普教育及媒體報(bào)道對(duì)于不實(shí)信息傳播的影響［J］．系統(tǒng)工程理論與實(shí)踐，2014，34(2) : 365－375. HUO Liang-an，HUANG Pei-qing．Study of the impact of science popularization and media coverage on the transmission of unconfirmed information［J］．Systems Engineering-Theory＆Practice，2014，34(2) : 365－375．

［16］MA Rin-go．Media，crisis，and SARS: an introduction ［J］．Asian Journal of Communication，2005，15(3) : 241 －246.

［17］郭淑利，張萍．媒體宣傳影響下的脈沖控制雙菌株SIR傳染病模型［J］．信陽(yáng)師范學(xué)院學(xué)報(bào)，2015，28 (1) : 21－24．GUO Shu-li，ZHANG Ping．A two-strain SIR epidemic model with impulsive control and media coverage［J］．Journal of Xinyang Normal University，2015，28(1) : 21 －24．

［18］劉玉英，肖燕妮．一類(lèi)受媒體影響的傳染病模型的研究［J］．應(yīng)用數(shù)學(xué)和力學(xué)，2013，34(4) : 399－407. LIU Yu-ying，XIAO Yan-ni．An epidemic model with saturated media/psychological impact［J］．Applied Mathematics and Mechanics，2013，34(4) : 399－407．

［19］JEAN M Tchuenche，NOTHABO Dube，CLAVER P Bhunu，et al．The impact of media coverage on the transmission dynamics of human influenza［J］．Bmc Public Health，2011，11 (suppl) : S5．

［20］LIU Rong-song，WU Jian-hong，ZHU Huai-ping．Media/ Psychological impact on multiple outbreaks of emerging infectious diseases［J］．Computational＆Mathematical Methods in Medicine，2007，8(3) : 153－164．

［21］XIAO Dong-mei，RUAN Shi-gui．Global analysis of an epidemic model with nonmonotone incidence rate［J］．Mathematical Biosciences，2007，208(2) : 419－429．

［22］LI Bo，YUAN San-ling，ZHANG Wei-guo．Analysis on an epidemic model with a ratio-dependent nonlinear incidence rate［J］．International Journal of Biomathematics，2011，4(2) : 227－239．

［23］SUN Cheng-jun，YANG Wei，ARINO Julien，et al．Effect of media-induced social distancing on disease transmission in a two patch setting［J］．Mathematical Biosciences，2011，230(2) : 87－95．

［24］CUI Jing-an，SUN Yong-hong，ZHU Huai-ping．The impact of media on the control of infectious diseases［J］．Journal of Dynamics and Differential Equations，2008，20 (1) : 31－53．

［25］Misra A K，Sharma Anupama，Shukla J B．Modeling and analysis of effects of awareness programs by media on the spread of infectious diseases［J］．Mathematical＆Computer Modelling，2011，53(5－6) : 1 221－1 228．

［26］潘艷雪，王曉靜，許傳青．具有時(shí)滯和媒體播報(bào)效應(yīng)的流感模型［J］．北京建筑大學(xué)學(xué)報(bào)，2014，30(3) : 64 －69. PAN Yan-xue，WANG Xiao-jing，XU Chuan-qing．Model of influenza disease with delay and media coverage［J］．Journal of Beijing University of Civil Engineering and Architecture，2014，30(3) : 64－69．

［27］馬知恩，周義倉(cāng)．常微分方程定性與穩(wěn)定性方法［M］．北京:科學(xué)出版社，2001. MA Zhi-en，ZHOU Yi-cang．The qualitative analysis and stability methods for ordinary differential equations ［M］．Beijing: Science Press，2001.

［28］Lakshmikantham V，Leela S，Martynyuk A A．Stability analysis of nonlinear systems［M］．New York: Marcel Dekker Inc．，1989.

［29］廖曉昕．穩(wěn)定性的理論、方法和應(yīng)用［M］．武漢:華中科技大學(xué)出版社，1990. LIAO Xiao-xin．Theory，methods and application of stability［M］．Wuhan: Huazhong University of Science and Technology Press，1990．

［30］陳蘭蓀，陳鍵．非線性生物動(dòng)力系統(tǒng)［M］．北京:科學(xué)出版社，1993. CHEN Lan-sun，CHEN Jian．Nonlinear biology dynamics ［M］．Beijing: Science Press，1993．

Qualitative analysis of a rumor transmission model with media coverage

CHEN Hua
(1．College of Energy Science and Engineering，Xi’an University of Science and Technology，Xi’an 710054，China; 2．President’s Office，Yan’an University，Yan’an 716000，China)

Abstract:To understand the mechanism of the rumor propagation and inhibit the spread of rumors．We consider a rumor transmission model with media coverage．Based on the stability of ordinary differential equation，detailed qualitative analysis about the existence and local asymptotic stability of its equilibria is carried out．Using comparison theorem，we prove the globally asymptotic stability of the equilibrium E0(rumor-free equilibrium)．By use of Lyapunov-LaSalle invariant principle，we give some sufficient conditions which lead to the globally asymptotic stability of the equilibrium E+(rumor-endemic equilibrium)．We show that the number of infective individuals tends to zero as time evolves or the rumor persists．Furthermore，the numerical simulations show the impact of media coverage on the spread and control of rumor spreading．

Key words:rumor propagation; stability; lyapunov function; media broadcast

作者簡(jiǎn)介:陳華(1967－)，男，山東曲阜人，副教授，E-mail: chenhua1130@163．com

*收稿日期:2015－10－10責(zé)任編輯:高佳

DOI:10．13800/j．cnki．xakjdxxb．2016．0217

文章編號(hào):1672－9315(2016) 02－0255－10

中圖分類(lèi)號(hào):O 175; N 945

文獻(xiàn)標(biāo)志碼:A