发明名称 一种复杂网络中的谣言传播控制方法
摘要 本发明公开了一种复杂网络中的谣言传播控制方法,在SIR模型为基础上建立一个新的谣言传播模型,考虑到感染率随着传播谣言节点数量增加反而下降,引入一个分段函数来描述感染率,这样就能更准确得描述复杂网络中的谣言传播行为。其次,在该新模型的基础上,引入一个最优控制变量,通过数学方法可以求出最优控制变量,将最优控制变量加入到该模型中,可以使得尽量多的健康节点转变为免疫节点,从而使网络中传播谣言的节点最少,达到控制谣言传播的目的。
申请公布号 CN104361231A 申请公布日期 2015.02.18
申请号 CN201410633493.2 申请日期 2014.11.11
申请人 电子科技大学 发明人 徐杰;余雅红;高成毅;孙健
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 一种复杂网络中的谣言传播控制方法,其特征在于,包括以下步骤:(1)、建立谣言传播模型将logistic增长规律引入到SIR谣言传播模型上,建立谣言传播模型:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mi>dS</mi><mi>dt</mi></mfrac><mo>=</mo><mi>&gamma;S</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>S</mi><mi>k</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mi>&alpha;SI</mi><mo>-</mo><mi>&mu;S</mi></mtd></mtr><mtr><mtd><mfrac><mi>dI</mi><mi>dt</mi></mfrac><mo>=</mo><mi>&alpha;SI</mi><mo>-</mo><mi>&beta;I</mi><mo>-</mo><mi>&mu;I</mi></mtd></mtr><mtr><mtd><mfrac><mi>dR</mi><mi>dt</mi></mfrac><mo>=</mo><mi>&beta;I</mi><mo>-</mo><mi>&mu;R</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000605396950000011.GIF" wi="729" he="421" /></maths>其中,S表示时间t内健康节点数量,I表示时间t内传播节点数量,R表示时间t内免疫节点数量,α表示当谣言传播节点与健康节点接触时健康节点以概率α变为传播节点I,β表示当谣言传播节点与免疫节点接触时传播节点以概率β变为免疫节点,μ表示单位时间移出现有网络的用户数量,将网络中所有的节点看成一个群体,则γ为这个群体的内在增长率且γ>μ,k为环境容量;设α是固定的,则健康节点S转变为传播节点I用分段函数T(I)表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mi>&alpha;I</mi></mtd><mtd><mi>I</mi><mo>&lt;</mo><mo>=</mo><msub><mi>I</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mi>&alpha;</mi><msub><mi>I</mi><mn>0</mn></msub></mtd><mtd><mi>I</mi><mo>></mo><mo>=</mo><msub><mi>I</mi><mn>0</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000605396950000012.GIF" wi="461" he="156" /></maths>其中,I<sub>0</sub>表示当健康节点S转变为传播节点I的数量达到饱和时的传播节点I的个数,用T(I)来替代谣言传播模型中的αI,则谣言传播模型可以表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mi>dS</mi><mi>dt</mi></mfrac><mo>=</mo><mi>&gamma;S</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mi>S</mi><mi>k</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mi>ST</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mo>-</mo><mi>&mu;S</mi></mtd></mtr><mtr><mtd><mfrac><mi>dI</mi><mi>dt</mi></mfrac><mo>=</mo><mi>ST</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mo>-</mo><mi>&beta;I</mi><mo>-</mo><mi>&mu;I</mi></mtd></mtr><mtr><mtd><mfrac><mi>dR</mi><mi>dt</mi></mfrac><mo>=</mo><mi>&beta;I</mi><mo>-</mo><mi>&mu;R</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi>b</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000605396950000013.GIF" wi="770" he="420" /></maths>(2)、求解最优控制变量引入控制变量u(t),u(t)∈U<sub>ad</sub>,U<sub>ad</sub>表示控制空间且U<sub>ad</sub>={u(t):0≤u(t)≤1,t∈[0,t<sub>end</sub>]};建立目标泛函数J(u),<img file="FDA0000605396950000021.GIF" wi="548" he="156" />τ表示控制预算参数,且<img file="FDA0000605396950000022.GIF" wi="390" he="205" />选取合适的控制变量u(t),使得目标泛函数达到最小,即求解出最优控制变量,记为:u<sup>*</sup>(t);(3)、控制谣言传播节点I的数量将最优控制变量u<sup>*</sup>(t)加入到谣言传播模型(b)中,得到受控下的微分方程组,如下:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msup><mi>dS</mi><mo>*</mo></msup><mi>dt</mi></mfrac><mo>=</mo><msup><mi>&gamma;S</mi><mo>*</mo></msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><msup><mi>S</mi><mo>*</mo></msup><mi>k</mi></mfrac><mo>)</mo></mrow><mo>-</mo><msup><mi>S</mi><mo>*</mo></msup><mi>T</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>&mu;S</mi><mo>*</mo></msup><mo>-</mo><msup><mi>u</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>S</mi><mo>*</mo></msup></mtd></mtr><mtr><mtd><mfrac><msup><mi>dI</mi><mo>*</mo></msup><mi>dt</mi></mfrac><mo>=</mo><msup><mi>S</mi><mo>*</mo></msup><mi>T</mi><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>&beta;I</mi><mo>*</mo></msup><mo>-</mo><msup><mi>&mu;I</mi><mo>*</mo></msup></mtd></mtr><mtr><mtd><mfrac><msup><mi>dR</mi><mo>*</mo></msup><mi>dt</mi></mfrac><msup><mrow><mo>=</mo><mi>&beta;I</mi></mrow><mo>*</mo></msup><mo>-</mo><msup><mi>&mu;R</mi><mo>*</mo></msup><mo>+</mo><msup><mi>u</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msup><mi>S</mi><mo>*</mo></msup></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000605396950000023.GIF" wi="892" he="439" /></maths>其中,S<sup>*</sup>、I<sup>*</sup>、R<sup>*</sup>分别表示在时间t内受控下的健康节点、传播节点、免疫节点数量,再通过求解微分方程组可以看到加入最优控制变量u<sup>*</sup>(t)之后,谣言传播节点I的数量得到控制。
地址 611731 四川省成都市高新区(西区)西源大道2006号