媒体干预下三分意见群体网络舆情传播模型的制作方法

本发明涉及网络安全领域，具体涉及到媒体作用下的网络舆情传播安全。

背景技术：

媒体作用下的网络舆情传播引起了学者们的广泛关注，但大部分研究均局限于通过加入媒体节点，对经典复杂网络传播模型进行改进，关注于媒体节点本身在模型上的体现。

近年来，对媒体作用下互联网舆情传播定量研究工作现已成为热点问题，musselle考察了社会网络中两个持对立观点媒体对个体意见传播演化的影响，发现当对立媒体连接个体数量越多，随着时间变化个体意见演化受媒体意见影响越大。crokidakis把大众媒体作为外部影响因素引入二维sznajd模型中，发现在媒体因素影响下，即使初始持赞同观点的个体密度很小，但受媒体传播影响，最终个体意见也会达成统一。

互联网舆情传播研究工作通常借鉴传染病模型，如利用sir(susceptible,infected,recovered)模型和sis(susceptible-infected-susceptible)模型来模拟互联网舆情传播。自sudbury首次将传染病sir模型用于解释谣言传播以来，不断有学者提出改进模型来描述网络舆情传播。如：顾亦然提出一个新的基于在线社交网络的谣言传播seir(susceptible,exposedinfected,recovered)模型，并给出了抑制谣言传播的免疫策略。钱颖等基于sir模型对微博上的舆情传播进行了研究，建立了微博舆情传播模型。xu等通过将信息价值与用户行为引入社交网络信息传播模型中，改进和修改了传统的seir模型，构建了基于社交网络的s-seir(socialsusceptible,exposedinfected,recovered)模型。李林等利用网络传播的sir模型，考虑两种冲突信息在社会网络上同时传播，传播过程中引入个体处理冲突信息的三种行为方式，建立了信息的传播演化模型。张立凡等通过构建带有讨论机制的舆情传播模型,研究媒体作用下网络舆情演变的内在规律,把媒体对舆情传播的干预作用抽象为强化度及分歧度,构建新的siaibr模型。本发明在siaibr模型的基础上，引入中立群体构建新的si³r模型。

技术实现要素：

在考虑网络媒体传播实际情况后，提出如下带有媒体作用讨论机制的网络舆情话题传播模型。该模型中网络群体处于5种状态：未知态(s)、支持传播态(ia)、中立传播态(ib)、反对传播态(ic)、移除态(r)。用s(t)、ia(t)、ib(t)、ic(t)、r(t)分别表示t时刻处于5种状态的网民占网民总数的比率。本发明中还用到了强化度、分歧度和渗透率三个概念，下面给出了相关的三个定义。

定义1强化度ε。在话题出现的时候，媒体的导向会使得未感染人群加入话题行列，表现为未知状态往移除状态转变的概率减小，并完整叠加到未知状态转变为转播状态的概率。

定义2分歧度η。当媒体观点与原有观点发生碰撞时，会导致传播过程中激烈的争论，从而产生意见分歧，改变传播态势。分歧度定义为一个概率，表示一部分传播者因受到媒体的影响而持另一种观点的概率。

定义3渗透率γ。表示持不同观点的个体在交互过程中，因受到他方的影响而改变之前对时间的看法，进而转变到他方的阵营中。这里假设有三方：支持方ia,中立方ib,反对方ic。这种因接触而发生状态转移的概率定义为渗透度，这里假设支持方与反对方的转换必须经过中立方，即不能直接转换，相应的渗透率分别为：γab，γba，γbc，γcb。

本发明解决方案主要是通过在支持群体和反对群体之间引入了中立群体，其主要特征为：在媒体传播过程中，主要影响对象是中立传播个体，因为这些个体容易随意改变自己的观点，将这种中立状态的个体考虑进传播模型中，能够更好反应网络舆论网民的分布规律。本发明通过引入中立群体提出的si³r模型，对应的动力学方程为：

其中，αi(i＝0,1,2,3)、βj(j＝1,2,3)分别为状态转换概率；ia、ib、ic和r分别为支持层、中间层、反对层和移除状态；α0、α1、α2、α3表示未知个体向移除个体、支持个体、中立个体以及反对个体转变的概率；β1、β2、β3表示支持个体、中立个体和反对个体变为移除状态的概率；γab、γba、γbc、γcb分别表示支持个体转变为中立状态、中立个体转变为支持状态、中立个体转变为反对状态、反对状态转变为中立状态的概率；a为新网民输入率。

另外，本发明中所用到的传播参数为α、β、γ，它们之间的转换关系如下(其中ε∈(0,1)，η∈(0,0.5))：

下面对模型进行具体分析，求解平衡点。

上述模型(1)中前4个变量不含变量r，因此仅需考虑前4个方程构成的系统。得到如下模型：

传播动力学理论表面，网络传播动力系统中都会存在一个传播阈值r0，又叫基本再生数，它决定了网络舆情传播与否。当r0≤1时，传播将逐渐消失，当r0>1时，舆情信息将会在一定区域内传播。

令x＝(ia,ib,ic,s)^t，则模型(3)可以表示为dx＝f(x)-v(x),计算基本再生数，如下：

则

由a＝fv^-1,可得谱半径(基本再生数r0)如下：

经使用matlab的solve符号演算可知，上述齐次方程不存在精确解，即无法得到精确平衡点。但可以通过简化方程得到它的4个平衡点：

可以看出p1、p2、p3属于对称的3个点,参数p1中的α1,β1在p2、p3中分别以α2,β2或α3,β3替换。p1、p2、p3中ia、ib、ic分别有2个为0,在实际传播情况中代表只有一种讨论声音。

上述分析给出了动力学方程组出现的4组平衡点精确解。然而若结合媒体作用下的舆论传播模型分析可发现,由于η＜1–η,α1＞α2、α3,β1＜β2、β3,均衡状态下,ia不会比ib、ic小,因此p2、p3解不会出现。因此,当r0＞1时,基于媒体干预的带有讨论机制的舆情传播模型只有一个平衡点p1。

附图说明

结合附图和下面的具体实施方式，将更加清楚地了解本发明所提出的媒体干预下三分意见群体网络舆情传播模型。图1是本发明的流程框图，图2是本发明的模型示意图。

具体实施方式

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图1和具体实施方式对本发明作进一步详细的说明，这些附图均为简化的示意图，仅以示意方式说明本方面的基本结构，因此其仅显示与本发明有关的构成。

本发明的实验平台为matlab2012a，计算机主机配置：cpu为i3_2310m,内存6g；操作系统为microsoftwindows10professional。

请参考图1，展示出了本发明的媒体干预下三分意见群体网络舆情传播模型算法100的具体实施例子的流程图。所述方法100包括：

步骤102，初始化变量：s(未知态)、ia(支持传播态)、ib(中立传播态)、ic(反对传播态)、r(移除态)，给定这五个变量的初始值；

步骤104，利用常微分方程数值解的欧拉法，设置迭代次数step＝150，针对下面的常微分方程组中的s(未知态)、ia(支持传播态)、ib(中立传播态)、ic(反对传播态)、r(移除态)这五个变量，进行循环求解；

步骤106，对求解出的五个变量进行输出。

本发明能够更好的反映传播群体内部的层次结构，符合实际谣传传播系统的特点。整个系统是在三个层次的传播(支持传播、中立传播和反对传播)影响下运行的，更能够反映网络舆论网民的分布规律。

上述说明已经充分表述了本发明的具体实施方式。需要指出的是，熟悉该领域的技术人员对本发明的具体实施方式所做的任何改动均不脱离本发明的权利要求书的范围。相应地，本发明的权利要求的范围也并不仅仅局限于所述具体实施方式。