删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于真实信息传播者的谣言传播模型的动力学分析

本站小编 Free考研考试/2021-12-29

摘要:在谣言传播过程中加入真实信息的传播者, 考虑了人们对谣言的遗忘因素, 建立了SITR (susceptible-infective-true-removed)谣言传播模型. 利用下一代矩阵得到了谣言传播的阈值K0, 证明了K0 < 1时无谣言传播者无真实信息传播者平衡点的稳定性, 给出了边界平衡点(即有谣言传播者但无真实信息传播者, 及无谣言传播者但有真实信息传播者平衡点)存在的条件, 以及它们的稳定性, 发现了两个边界平衡点出现双稳的区域, 获得了不同条件下正平衡点的存在性, 及其局部稳定性. 最后, 通过数值模拟验证了理论结果, 模拟分析了真实信息传播者的初始值对谣言传播者的峰值及谣言的持续时间等的影响.
关键词: 谣言传播/
真实信息/
阈值/
平衡点

English Abstract


--> --> -->
谣言, 指的是没有相应事实基础, 却被捏造出来并通过一定手段推动传播的言论. 近年来, 随着移动互联网 + 时代的到来, 使信息化传播门槛进一步降低, 加之移动通信工具的普遍, 谣言传播变得速度更快, 作用力更强. 网络谣言滋生蔓延, 不仅败坏个人名誉, 给受害人造成极大的精神困扰, 更损害国家形象, 影响社会稳定. 因此, 对于谣言传播规律的研究十分必要. 谣言传播模型的研究始于20世纪60年代, Daley和Kendall[1]提出的D-K模型, 按照谣言传播把受众分为三类, 未曾听过谣言的人、传播谣言的人和听过谣言但不传播谣言的人, 利用随机过程的方法进行分析. Sudbury[2]认为可以借鉴传染病SIR (susceptible-infective-removed)模型研究谣言的传播. Zanette[3,4]和Nekovee等[5]分别在小世界网络和无标度网络上研究谣言的传播. Moreno等[6]基于SIR模型在非均匀网络上建立网络谣言传播模型. Zhao等[7,8]研究提出了遗忘机制、记忆机制对谣言传播的影响. Deng和Li[9]在研究谣言传播模型时考虑了遗忘率和记忆率. Zan等[10]基于SIR模型提出了具有反击机制的谣言传播模型. Wang等[11]基于SIR模型提出了带有信任机制的谣言传播模型. Yang等[12]和He等[13]研究了真理传播和辟谣策略对于谣言传播的影响. Huo等[14]考虑了传播节点的活跃度, 将谣言传播者分为高活跃状态下的传播者和低活跃状态下的传播者, 最后得出谣言传播者的活跃程度对于谣言传播有着非常大的影响. Zhang和Zhu[15]建立了I2S2R (ignorants-spreaders-stifler)的谣言传播动态模型. Xiao等[16]提出了一种基于进化博弈和反谣言信息的谣言传播动力学模型, 证明谣言不仅受到反谣言信息的影响, 而且还受到用户行为和心理因素的影响. 顾亦然和夏玲玲[17]基于在线社交网络建立SEIR (susceptible-exposed-infective-removed)谣言传播模型, 并提出重要熟人免疫策略. 王辉等[18]研究了社会加强效应对谣言传播的影响. 张亚明等[19]研究了双重社会强化对谣言传播的影响. 万贻平等[20]在谣言传播模型的基础上增加谣言清除者, 研究了谣言清除过程的网络谣言传播与抑制. 冉茂洁等[21]同时考虑个体兴趣度差异和辟谣机制问题, 提出建立了IWSR (ignorant-weak spreader-strong spreader-removal)谣言传播模型. 赵敏等[22]考虑媒体正面宣传报道和负面报道对谣言传播的影响, 基于社交网络建立了新的谣言传播模型.
在谣言传播过程中, 一些人已知事实真相或经过判断知道真相, 并传播有关谣言的真实信息, 但是已有的研究未考虑真实信息传播对于谣言传播的动力学影响. 因此加入真实信息的传播者会使得谣言传播更加符合实际, 而且现有的谣言模型对其进行数学分析的较少. 本文考虑真实信息传播的影响, 建立谣言传播模型并对其进行动力学分析, 第2节根据谣言传播的机理建立了含有真实信息传播者的谣言传播SITR模型; 第3节给出了谣言传播的阈值条件; 第4节讨论了谣言传播者存在真实信息传播者不存在平衡点的条件, 及谣言传播者不存在真实信息传播者存在平衡点的条件, 并分析了它们的动力学性态; 第5节分析了谣言传播者与真实信息传播者共同存在的平衡点存在性及其动力学性态; 第6节数值模拟验证理论结果; 第7节为总结.
在谣言传播过程中加入真实信息的传播者, 将人群分为四类: 未知者S (susceptible), 谣言的传播者I (infective), 真实信息的传播者T (true), 免疫者R (removed). 未知者与谣言传播者接触率为$\beta$, 接触后未知者以概率$\alpha_1$成为谣言传播者, 以概率$1-\alpha_1$成为免疫者; 未知者和谣言传播者与真实信息传播者接触, 未知者和谣言传播者以概率$\delta$成为真实信息传播者; 谣言传播者与谣言传播者接触后, 谣言传播者对谣言失去兴趣变为免疫者的概率为$\epsilon$; 谣言传播者的遗忘率为$\theta$.
$S(t),\; I(t),\; T(t)$$R(t)$分别表示在t时刻的未知者、谣言的传播者、真实信息的传播者和免疫者的数量, $N(t)$表示在t时刻四类人群的总数, 则有$N(t) = I(t)+S(t)+T(t)+R(t).$ 假设个体以恒定速率A进入人群且全部为未知者, 四类人群以速率$\mu$移出群体, 他们之间的关系如图1所示.
图 1 真实信息影响下谣言传播流程图
Figure1. Flow chart of rumor propagation under the influence of true information.

根据图1, 建立SITR谣言传播模型如下:
$\left\{\begin{aligned} & \dfrac{{\rm d} S(t)}{{\rm d}t} = A-\beta SI-\delta ST-\mu S,\\ & \dfrac{{\rm d} I(t)}{{\rm d}t} = \alpha_1\beta SI-\delta IT-\epsilon I^2-\theta I-\mu I,\\ &\dfrac{{\rm d} T(t)}{{\rm d}t} = \delta T(S+I)-\mu T,\\ &\dfrac{{\rm d} R(t)}{{\rm d}t} = (1-\alpha_1)\beta SI+\epsilon I^2+\theta I-\mu R.\end{aligned} \right. $
将系统(1)的四个方程相加, 得到$\dfrac{{\rm d} N(t)}{{\rm d}t} = A-$$\mu N(t),$ 解得$N(t) = N_0{\rm e}^{-t}+\dfrac{A}{\mu},$ 其中$N(0) = N_0$表示初始时刻四类人群的总数. 有
${\lim_{t \to +\infty}N(t) = \dfrac{A}{\mu}}.$
因为系统(1)的前三个方程不含R, 因此只考虑系统(1)的前三个方程, 得
$\left\{\begin{aligned} &\dfrac{{\rm d} S(t)}{{\rm d}t} = A-\beta SI-\delta ST-\mu S,\\ &\dfrac{{\rm d} I(t)}{{\rm d}t} = \alpha_1\beta SI-\delta IT-\epsilon I^2-\theta I-\mu I,\\ &\dfrac{{\rm d} T(t)}{{\rm d}t} = \delta T(S+I)-\mu T.\end{aligned}\right. $
容易验证系统(2)的正向不变集是
$\begin{split} D = & \left\{(S(t),\;I(t),\;T(t))|S(t),\;I(t),\;T(t)\geqslant 0, \right. \\ & \left. 0\leqslant S(t)+I(t)+T(t)\leqslant \dfrac{A}{\mu}\right\}. \end{split}$

系统(2)有一个无谣言传播者无真实信息传播者平衡点$E_0\left({A}/{\mu},\;0,\;0\right),$ 下面用下一代矩阵法[23]求系统(2)的阈值.
${ X} = (I,T)^\top$, $\top$表示向量$(I,T)$的转置, 则
${{\cal{ F}}(X)} \!=\! \left(\!\!\!\begin{array}{*{20}{c}} \alpha_1\beta SI\\ \delta TS \end{array}\!\!\!\right), \; {{\cal{ V}}(X)} \!=\! \left(\!\!\!\begin{array}{*{20}{c}} \delta IT+\epsilon I^2+\theta I+\mu I\\ \mu T-\delta IT \end{array}\!\!\!\right).$
在无谣言传播者无真实信息传播者平衡点$E_0$处有
${ F}\left( {{E_0}} \right) = \left(\!\!{\begin{array}{*{20}{c}}{{\alpha _1}\beta \dfrac{A}{\mu }}&0\\0&{\delta \dfrac{A}{\mu }}\end{array}}\!\!\right),\;\;\;{ V}\left( {{E_0}} \right) = \left(\!\!{\begin{array}{*{20}{c}}{\theta + \mu }&0\\0&\mu \end{array}}\!\!\right),$
${{ {FV}}^{-1}} = \left(\!\!\begin{array}{*{20}{c}} \dfrac{\alpha_1\beta A}{\mu(\theta+\mu)} & 0\\ 0 & {\delta A}/{\mu^2} \end{array}\!\!\right),$ 则系统(2)的阈值是${ {FV}}^{-1}$的最大谱半径, 即
$K_0 = \max\{\rho({ {FV}}^{-1})\} = \max\{K_1,K_2\},$
其中$K_1 = \dfrac{\delta A}{\mu^2},$ $K_2 = \dfrac{\alpha_1\beta A}{\mu(\theta+\mu)}.$
定理1  当$K_0 < 1$时, 系统(2)的无谣言传播者无真实信息传播者平衡点$E_0$是局部渐近稳定的.
2
4.1.无谣言传播者但真实信息传播者存在的平衡点
-->在谣言传播的过程中, 由于真实信息传播者的存在, 人群中最终没有了谣言传播者, 即无谣言传播者但有真实信息传播者平衡点存在性的问题. 通过对系统(2)求平衡点得, 如果$K_1 > 1$, 那么系统(2)存在无谣言传播者但有真实信息传播者平衡点$E_1\left(\dfrac{\mu}{\delta},0,\dfrac{\mu}{\delta}(K_1-1)\right)$, 简称为系统(2)的边界平衡点$E_1$.
定理2  当$K_1 > 1$, 且$K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 系统(2)的边界平衡点$E_1$是局部渐近稳定的.
证明  系统(2)在$E_1$处的雅可比矩阵为
${{ J}(E_1)} = \left( \begin{array}{*{20}{c}} -\mu K_1 & - \dfrac{\beta \mu}{\delta} & -\mu\\ 0 & \dfrac{\alpha_1\beta \mu}{\delta}-\mu K_1-\theta & 0 \\ \mu K_1-\mu & \mu K_1-\mu & 0 \end{array} \right),$
特征方程为: $(\lambda-\dfrac{\alpha_1\beta \mu}{\delta}+\mu K_1+\theta)[\lambda^2+\mu K_1\lambda+$$\mu^2(K_1-1)] = 0$. 因此, 特征方程有一个特征根为$\lambda_1 = \dfrac{\alpha_1\beta \mu}{\delta}-\mu K_1-\theta,$ 另外两个特征根满足方程:
$\begin{array}{*{20}{c}} \lambda^2+\mu K_1\lambda+\mu^2(K_1-1) = 0, \end{array} $
$K_1 > 1$时, 方程(3)的两个根都具有负实部. 又当$K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, $\lambda_1 < 0$. 因此, 当$K_1 > 1$, 且$K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, $E_1$是局部渐近稳定的.
2
4.2.谣言传播者存在但无真实信息传播者的平衡点
-->在谣言传播的过程中, 可能会出现谣言传播者存在但没有真实信息传播者, 即$I\neq0,$ $T = 0$. 令系统(2)右端为零, 当$I\neq0,$ $T = 0$时, 得到方程
$\begin{array}{*{20}{c}} aI^2+bI+c = 0, \end{array} $
其中$a = \epsilon\beta \!>\! 0,$ $b = \epsilon\mu+\!\beta(\theta\!+\!\mu) \!>\! 0,$ $c = \mu(\theta\!+\!\mu)\times$$(1-K_2).$
$K_2 < 1$时, 方程(4)无正实根. 当$K_2 > 1$时, 方程(4)存在唯一正根$I_2 = \dfrac{-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta}{2\epsilon\beta},$ 其中, $\varDelta = [\epsilon\mu+\beta(\theta+\mu)]^2-4\epsilon\beta\mu(\theta+\mu)(1-K_2).$ 故当$K_2 > 1$时, 系统(2)存在有谣言传播者但无真实信息传播者的平衡点$E_2\left(\dfrac{-\epsilon\mu+\beta(\theta+\mu)+\sqrt\varDelta}{2\alpha_1\beta^2}\right.,$$\left.\dfrac{-\epsilon\mu-\beta(\theta\!+\!\mu)\!+\!\sqrt\varDelta}{2\epsilon \beta},\;0\right)$, 我们称为边界平衡点$E_2$.
定理3  当$K_2 > 1$, 且
$K_1 < \dfrac{2\epsilon\beta(\theta+\mu)}{(\epsilon\beta-\alpha_1\beta^2)(\theta+\mu)+(\sqrt\varDelta-\epsilon\mu)(\epsilon+\alpha_1\beta)}K_2$时, 系统(2)的边界平衡点$E_2$是局部渐近稳定的.
证明  系统(2)在$E_2$处的雅可比矩阵为

$ {{ J}(E_2)}\!=\!\left(\!\!\!\!\begin{array}{*{20}{c}} \dfrac{-\epsilon\mu+\beta(\theta+\mu)-\sqrt\varDelta}{2\epsilon} & \dfrac{\epsilon\mu-\beta(\theta+\mu)-\sqrt\varDelta}{2\alpha_1\beta} & \dfrac{\delta[\epsilon\mu-\beta(\theta+\mu)-\sqrt\varDelta]}{2\alpha_1\beta^2}\\ \dfrac{\alpha_1[-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta]}{2\epsilon} & \dfrac{\epsilon\mu+\beta(\theta+\mu)-\sqrt\varDelta}{2\beta} & \dfrac{\delta[\epsilon\mu+\beta(\theta+\mu)-\sqrt\varDelta]}{2\epsilon\beta} \\ 0 & 0 & \delta\left[\dfrac{-\epsilon\mu+\beta(\theta+\mu)+\sqrt\varDelta}{2\alpha_1\beta^2}+\dfrac{-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta}{2\epsilon \beta}\right]-\mu \end{array}\!\!\!\!\right)\!,$
相应的特征方程为
$\begin{split} & \left\{\lambda+\mu-\delta\left[\dfrac{-\epsilon\mu+\beta(\theta+\mu)+\sqrt\varDelta}{2\alpha_1\beta^2}\right.\right.\\ + &\left.\left.\dfrac{-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta}{2\epsilon \beta}\right]\right\}(\lambda^2+p_0\lambda+q_0) = 0,\end{split}$
其中,
$p_0 = \dfrac{(\beta+\epsilon)[-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta]}{2\epsilon\beta}+\mu,$
$q_0 = \dfrac{\sqrt\varDelta[-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta]}{2\epsilon\beta}.$ 特征方程有一个特征根为
$\begin{split}\lambda_1 = & \;\delta\left[\dfrac{-\epsilon\mu+\beta(\theta+\mu)+\sqrt\varDelta}{2\alpha_1\beta^2}\right.\\ & \left.+\dfrac{-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta}{2\epsilon \beta}\right]-\mu,\end{split}$
另外两个特征根满足方程
$\begin{array}{*{20}{c}} \lambda^2+p_0\lambda+q_0 = 0, \end{array} $
因为$-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta > 0$, 从而有$p_0 > 0,$ $q_0 > 0$. 故方程(5)的两个根都具有负实部. 因此, 当
$\begin{aligned}\lambda_1 = & \; \delta\left[\dfrac{-\epsilon\mu+\beta(\theta+\mu)+\sqrt\varDelta}{2\alpha_1\beta^2}\right. \\ &+ \left.\dfrac{-\epsilon\mu-\beta(\theta+\mu)+\sqrt\varDelta}{2\epsilon \beta}\right]-\mu < 0\end{aligned}$
时, 即
$K_1 < \dfrac{2\epsilon\beta(\theta+\mu)}{(\epsilon\beta-\alpha_1\beta^2)(\theta+\mu)+(\sqrt\varDelta-\epsilon\mu)(\epsilon+\alpha_1\beta)}K_2,$
系统(2)的边界平衡点$E_2$是局部渐近稳定的. 令
$m = \dfrac{(\epsilon\beta-\alpha_1\beta^2)(\theta+\mu)+(\sqrt\varDelta-\epsilon\mu)(\epsilon+\alpha_1\beta)-2\epsilon\beta\theta}{2\epsilon\beta\mu},$ $n = \dfrac{2\epsilon\beta(\theta+\mu)}{(\epsilon\beta-\alpha_1\beta^2)(\theta+\mu)+(\sqrt\varDelta-\epsilon\mu)(\epsilon+\alpha_1\beta)}$. 直线$K_2 = \dfrac{1}{n}K_1$的斜率为
$\dfrac{1}{n} = \dfrac{(\epsilon\beta-\alpha_1\beta^2)(\theta+\mu)+(\sqrt\varDelta-\epsilon\mu)(\epsilon+\alpha_1\beta)}{2\epsilon\beta(\theta+\mu)}$, 因为$\sqrt\varDelta > \epsilon\mu+\beta(\theta+\mu)$, 所以${1}/{n} > 1$. 因此$K_2 = $$K_2 = \dfrac{1}{n}K_1$$ {K_1(\mu K_1+\theta)}/({\mu+\theta})$始终有一个交点, 且$K_2 = \dfrac{1}{n}K_1$$K_2 = {K_1(\mu K_1+\theta)}/({\mu+\theta})$的交点的横坐标为$m > 1$.
综上所述, 当$K_1 > 1,$$K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 系统(2)存在边界平衡点$E_1$, 且$E_1$是局部渐近稳定的(见图2). 当$K_2 > 1$, 且$K_1 < nK_2$时, 系统(2)存在边界平衡点$E_2$, 且$E_2$是局部渐近稳定的(见图2). 当$K_1 > m,$$\dfrac{1}{n}K_1 < K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 系统(2)存在两个边界平衡点$E_1$$E_2$, 且$E_1$$E_2$是局部渐近稳定的, 系统(2)出现双稳现象(见图2红色区域). 即在谣言传播的过程中, 初始值处于不同范围时, 系统(2)的解或者趋于有真实信息传播者但无谣言传播者的平衡点, 或者趋于有谣言传播者但无真实信息传播者的平衡点.
图 2 $K_1$$K_2$ 的关系图, 参数取值为$\alpha_1=0.6, $$\beta=0.6,\; A=0.2, \;\mu=0.2, \;\theta=0.01, \;\epsilon=0.01$
Figure2. Diagram of $K_1$ and $K_2$. The parameter values are $\alpha_1= 0.6,\; \beta=0.6,\; A=0.2,\; \mu =0.2, \;\theta =0.01,\; \epsilon=0.01 $

下面求在谣言传播中, 谣言传播者和真实信息传播者共存的平衡点, 为了简便起见称其为正平衡点.
令系统(2)右端为零, 当$I\neq0,$ $T\neq0$时, 解得
$S^* = \dfrac{\mu}{\delta}-I^*,\quad T^* = \dfrac{(\theta+\mu)\left(\dfrac{K_2}{K_1}-1\right)-(\epsilon+\alpha_1\beta)I^*}{\delta},$
其中, $I^*$满足方程
$\begin{array}{*{20}{c}} a_0I^*{^2}+b_0I^*+c_0 = 0, \end{array} $
$a_0 = \beta-\epsilon-\alpha_1\beta,$ $b_0 = (\theta+\mu)\dfrac{K_2}{K_1}-\theta-\dfrac{\mu}{\delta}a_0,$ $c_0 = $ $ \dfrac{\mu(\theta+\mu)}{\delta K_1}\left[\dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}-K_2\right].$ 方程(6)的判别式为
$\varDelta_1 = \left[(\theta+\mu)\dfrac{K_2}{K_1}-\theta+\dfrac{\mu}{\delta}a_0\right]^2-4a_0A.$
若正平衡点$E^*(S^*,\;I^*,\;T^*)$存在, 则要求$S^* > 0,$$ I^* > 0,\; T^* > 0$. 可得到$0 < I^* < {\mu}/{\delta},$ $0 < I^* <$ $ \dfrac{{(\theta\!+\!\mu)\left({K_2}/{K_1} \!-\! 1\right)}}{\epsilon+\alpha_1\beta}$. 因此, 要想使得正平衡点$E^*$存在, 只需$K_2 > K_1$. 又因为$\dfrac{{(\theta+\mu)\left(\dfrac{K_2}{K_1}-1\right)}}{\epsilon+\alpha_1\beta} < \dfrac{\mu}{\delta}$, 故$0 < I^* < \dfrac{{(\theta+\mu)\left(\dfrac{K_2}{K_1}-1\right)}}{\epsilon+\alpha_1\beta}$.
下面分情况讨论方程(6)正根的存在性.
情形1  当$K_2 > \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 有$c_0 < 0$, $(\theta+\mu)\dfrac{K_2}{K_1}-\theta > \dfrac{\delta}{\mu}A.$ 所以$\varDelta_1 \! > \!\left(\dfrac{\delta}{\mu}A-\dfrac{\mu}{\delta}a_0\right)^2 \!\geqslant \! 0$.
(i)当$a_0\geqslant 0,$$\epsilon\leqslant (1-\alpha_1)\beta$时, 方程(6)有唯一正实根$I_1{^*} = \dfrac{-b_0+\sqrt\varDelta_1}{2a_0}$. 因此, 系统(2)存在唯一正平衡点$E_1{^*}(S_1{^*},\; I_1{^*},\; T_1{^*})$.
(ii)当$a_0 < 0$, 即$\epsilon > (1-\alpha_1)\beta$时, 可由$(\theta+\mu)\left(\dfrac{K_2}{K_1}-1\right) > 0$, 得到$(\theta+\mu)\dfrac{K_2}{K_1}-\theta > 0$, 故有$b_0 > 0$, 此时方程(6)有两个正实根$I_2{^*} = \dfrac{-b_0+\sqrt\varDelta_1}{2a_0}$$I_3{^*} = \dfrac{-b_0-\sqrt\varDelta_1}{2a_0} > \dfrac{\mu}{\delta}$(舍去). 因此, 系统(2)存在唯一的正平衡点$E_2{^*}(S_2{^*}, \;I_2{^*}, \;T_2{^*})$.
情形2 当$K_2 = \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 有$c_0 = 0$, $\varDelta_1 = b_0{^2} > 0.$
(i)当$a_0 > 0,$ $b_0 < 0$, 即$\epsilon < (1-2\alpha_1)\beta+\dfrac{\delta\theta}{\mu}$时, 方程(6)有唯一正实根$I_4{^*} = \dfrac{-b_0}{a_0}$. 因此, 系统(2)存在唯一正平衡点$E_4{^*}(S_4{^*}, I_4{^*}, T_4{^*})$.
(ii)当$a_0 < 0$, 即$\epsilon > (1-\alpha_1)\beta$时, 有$b_0 > 0$, 此时方程(6)有唯一正实根$I_5{^*} = \dfrac{-b_0}{a_0} > \dfrac{\mu}{\delta}$(舍去).
情形3 当$K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 有$c_0 > 0$.
(i)当$a_0 > 0$, $b_0 < 0$, $\varDelta_1\geqslant 0$时, 即$\epsilon < $$ (1-2\alpha_1)\beta+\dfrac{\delta\theta}{\mu}$$\dfrac{\mu(\beta-\epsilon)}{\delta(\theta+2\sqrt {a_0A})}\;\geqslant\; 1$, 此时方程(6)存在两个正实根$I_6{^*} = \dfrac{-b_0+\sqrt\varDelta_1}{2a_0}$$I_7{^*} = $$ \dfrac{-b_0-\sqrt\varDelta_1}{2a_0}$. 因此, 系统(2)存在两个正平衡点$E_6{^*}(S_6{^*},\; I_6{^*},\; T_6{^*})$$E_7{^*}(S_7{^*}, \;I_7{^*},\; T_7{^*})$.
(ii)当$a_0 < 0$, 即$\epsilon > (1-\alpha_1)\beta$时, 有$b_0 > 0$, $\varDelta_1 > 0$. 此时方程(6)有唯一正实根$I_8{^*} = $$\dfrac{-b_0-\sqrt\varDelta_1}{2a_0} > \dfrac{\mu}{\delta}$(舍去). 因此, 系统(2)的正平衡点存在性仅在$K_2 > K_1$范围内考虑.
综上所述, 系统(2)平衡点的存在性分下面两种情形:
1) $\epsilon < (1-\alpha_1)\beta$: 当$K_2 > \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 若$K_2 < 1$, $E_1{^*}$存在, 此时系统(2)存在两个平衡点$E_0,$ $E_1{^*}$(见图3(a)黄色区域). 若$K_2 > 1,K_1 < 1$时, $E_1{^*}$存在, 此时系统(2)存在三个平衡点$E_0,$ $E_2,$ $E_1{^*}$(见图3(a)绿色区域). 若$K_2 > 1,$ $K_1 > 1$时, $E_1{^*}$存在, 此时系统(2)存在四个平衡点$E_0,$ $E_1,$ $E_2,$ $E_1{^*}$(见图3(a)紫色区域). 当$K_2 < $$ \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$, 且满足$\epsilon < (1-2\alpha_1)\beta+\dfrac{\delta\theta}{\mu}$, $\varDelta_1\geqslant 0$时, $E_6{^*}$$E_7{^*}$存在, 此时系统(2)存在五个平衡点$E_0,$ $E_1,$ $E_2,$ $E_6{^*},$ $E_7{^*}$(见图3(a)蓝色区域). 若$1 < K_2 < K_1$, 此时系统(2)存在三个平衡点$E_0,$ $E_1,$ $E_2$(见图3(a)白色区域). 若$K_1 > 1 > $$K_2$, 此时系统(2)存在两个平衡点$E_0,$ $E_1$(见图3(a)红色区域). 若$1 > K_1 > K_2$时, 此时系统(2)存在一个平衡点$E_0$存在(见图3(a)橙色区域).
图 3 $K_1$$K_2$关系图($\alpha_1=0.7, \beta=0.9, A=0.2,\ \mu=$$0.2,\; \theta=0.01)$ (a) $\epsilon=0.01$, $\epsilon<(1-\alpha_1)\beta$; (b) $\epsilon=0.3$, $\epsilon>(1-\alpha_1)\beta$
Figure3. Diagram of $K_1$ and $K_2$ ($\alpha_1\!= 0.7,\; \beta\!=0.9,\; A\!=0.2,$ $\mu=0.2,\; \theta=0.01$): (a) $\epsilon=0.01$, $\epsilon<(1-\alpha_1)\beta$; (b) $\epsilon=$ $0.3,\; \epsilon > (1-\alpha_1)\beta$

2) $\epsilon > (1-\alpha_1)\beta$: 具体见图3(b).
因此, 我们进一步表示系统(2)平衡点之间的关系见表1.
系统(2)的平衡点
$ \epsilon< \left( {1 - {\alpha _1}} \right){\rm{\beta }}$${K_2} > {K_1}$${K_2} > \dfrac{{{K_1}\left( {\mu {R_1} + \theta } \right)}}{{\mu + \theta }}$${K_2} < 1$图3(a)中黄色区域${E_0},\;E_1^*$
${K_2} > 1 > {K_1}$图3(a)中绿色区域${E_0},\;{E_2},\;E_1^*$
${K_1} > 1$图3(a)中紫色区域${E_0},\;{E_1},\;{E_2},\;E_1^*$
${K_2} < \dfrac{{{K_1}\left( {\mu {K_1} + \theta } \right)}}{{\mu + \theta }}$$\epsilon < \left( {1 - 2{\alpha _1} } \right)\beta + \dfrac{ {\delta \theta } }{\mu },{\varDelta _1} \geqslant 0$图3(a)中蓝色区域${E_0},\;{E_1},\;{E_2},\;E_6^*,\;E_7^*$
${K_2} < {K_1}$${K_1} < 1$图3(a)中橙色区域$E_0$
${K_1} > 1 > {K_2}$图3(a)中红色区域$E_0,\;E_1$
${K_1} > {K_2} > 1$图3(a)中空白区域$E_0,\;E_1,\;E_2$
$\epsilon > \left( {1 - {\alpha _1}} \right){\rm{\beta }}$${K_2} > {K_{1,}}$${K_2} > \dfrac{{{K_1}\left( {\mu {K_1} + \theta } \right)}}{{\mu + \theta }}$${K_2} < 1$图3(b)中黄色区域$E_0,\;E_2^*$
${K_2} > 1 > {K_1}$图3(b)中绿色区域$E_0,\;E_2,\;E_2^*$
${K_1} > 1$图3(b)中紫色区域$E_0,\;E_1,\;E_2,\;E_2^*$
$K_2<K_1$$K_1<1$图3(b)中橙色区域$E_0$
$K_1>1>K_2$图3(b)中红色区域$E_0,\;E_1$
$1 < {K_2} < \dfrac{{{K_1}\left( {\mu {K_1} + \theta } \right)}}{{\mu + \theta }}$图3(b)中空白区域$E_0,\;E_1,\;E_2$


表1系统(2)平衡点之间的关系
Table1.Relationship table between the equilibria of system (2).

定理4  当$K_2 > \max\{K_1, \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}\}$, 且$\epsilon < (1-\alpha_1)\beta$时, 如果
$\dfrac{[(a_0\!+\!\epsilon)I_1{^*}\!+\!g]\left[a_0\epsilon I_1{^*}{^2}\!+\!\left(g\epsilon \!+\!\alpha_1\beta^2\dfrac{\mu}{\delta}\right)I_1{^*}\!-\!\mu h\right]}{\alpha_1\beta^2(a_0+\epsilon)I_1{^*}{^3}+\alpha_1\beta^2gI_1{^*}{^2}-\delta h\sqrt\varDelta_1} \! >\!1,$ 那么系统(2)的正平衡点$E_1{^*}$是局部渐近稳定的.
证明  系统(2)在正平衡点$E^*(S^*,I^*,T^*)$处的雅可比矩阵为

${{ J}(E^*)} = \left( \begin{array}{*{20}{c}} -\beta I^*-\delta T^*-\mu & -\beta S^* & -\delta S^*\\ \alpha_1\beta I^* & \alpha_1\beta S^*-\delta T^*-2\epsilon I^*-\theta-\mu & -\delta I^* \\ \delta T^* & \delta T^* & \delta (S^*+I^*)-\mu \end{array} \right),$
其中$S^* = \dfrac{\mu}{\delta}-I^*,\; T^* = \dfrac{(\theta+\mu)\left({K_2}/{K_1}-1\right)-(\epsilon+\alpha_1\beta)I^*}{\delta},$ 代入化简得
${{ J}(E^*)} = \left( \begin{array}{*{20}{c}} -a_0I^*-(\theta+\mu)\dfrac{K_2}{K_1}+\theta & -\beta \dfrac{\mu}{\delta}+\beta I^* & -\mu+\delta I^*\\ \alpha_1\beta I^* & -\epsilon I^* & -\delta I^* \\ (\theta+\mu)(\dfrac{K_2}{K_1}-1)-(\epsilon+\alpha_1\beta) I^* & (\theta+\mu)\left(\dfrac{K_2}{K_1}-1\right)-(\epsilon+\alpha_1\beta) I^* & 0 \end{array} \right),$
相应的特征方程为
$\begin{array}{*{20}{c}} \lambda^3+b_1\lambda^2+b_2\lambda+b_3 = 0, \end{array}$
其中, $b_1 = (a_0+\epsilon)I^*+g,$ $b_2 = (a_0\epsilon-\alpha_1\beta^2)I^*{^2}+\left(g\epsilon+\alpha_1\beta^2\dfrac{\mu}{\delta}\right)I^*-\mu h,$ $b_3 = -\delta h(2 a_0I^*{^2}+b_{0}I^*),$ $g = $$ (\theta+\mu)\dfrac{K_2}{K_1}-\theta > 0,$ $h = (\epsilon+\alpha_1\beta) I^*-(\theta+\mu)\left(\dfrac{K_2}{K_1}-1\right) < 0$. 计算得
$W_1 = b_1 = (a_0+\epsilon)I^*+g > 0,$
$\begin{split}{W_2} = & \left| \begin{array}{*{20}{c}} b_1 & 1 \\ b_3 & b_2 \end{array} \right| = [(a_0+\epsilon)I^*+g]\left[a_0\epsilon I^*{^2}+\left(g\epsilon+\alpha_1\beta^2\dfrac{\mu}{\delta}\right)I^*-\mu h\right]\\ & +h[2\delta a_0I^*{^2}-(\mu a_0-g\delta)I^*]-\alpha_1\beta^2(a_0+\epsilon)I^*{^3}-\alpha_1\beta^2gI^*{^2}\end{split}$


${W_3} = \left| \begin{array}{*{20}{c}} b_1 & 1 &0\\ b_3 & b_2 &b_1\\ 0 & 0 &b_3\\ \end{array} \right| = b_3W_2.$
$I_1{^*}$代入得$b_3 = -\delta h\sqrt\varDelta_1 > 0,$ 从而知$W_3$$W_2$同号. 进一步把$I_1{^*}$代入$W_2$, 由Hurwitz 判据[24]可知, 当$W_2 > 0$, 即
$\dfrac{[(a_0+\epsilon)I_1{^*}+g]\left[a_0\epsilon I_1{^*}{^2}+\left(g\epsilon+\alpha_1\beta^2\dfrac{\mu}{\delta}\right)I_1{^*}-\mu h\right]}{\alpha_1\beta^2(a_0+\epsilon)I_1{^*}{^3}+\alpha_1\beta^2gI_1{^*}{^2}-\delta h\sqrt\varDelta_1} > 1$时, $E{_1}^{*}$是局部渐近稳定的.
注: 1)可类似证明其他正平衡点的稳定性.
2)当$\epsilon = 0$时, 阈值$K_0$, 无谣言传播者无真实信息传播者平衡点$E_0$, 及无谣言传播者但真实信息传播者存在的平衡点$E_1$没有变化; 谣言传播者存在但无真实信息传播者的平衡点$E_2$的大小发生变化, 影响着双稳区域; 谣言传播者真实信息传播者共存的平衡点(即正平衡点)存在的条件及个数发生变化, 影响着其动力学性态.
为了验证理论结果, 下面对系统(2)进行数值模拟.
参数取值为$\alpha_1 = 0.3$, $\beta = 0.6$, $A = 0.2$, $\mu = $$ 0.2$, $\theta = 0.01$, $\epsilon = 0.01$, $\delta = 0.15$时, $K_1 = 0.75 < 1,$ $K_2 = 0.857 < 1$, 则$E_0$渐近稳定(见图4).
图 4 $\alpha_1 = 0.3,\;\beta = 0.6,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.15$时, (a) S, I, T的时间序列图和(b) I-T相平面图
Figure4. When $\alpha_1 = 0.3,\;\beta = 0.6,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.15$, (a) time series graph of S, I, and T; (b) phase plan of I-T.

参数取值为$\alpha_1 = 0.7$, $\beta = 0.9$, $A = 0.2$, $\mu = $$ 0.2$, $\theta = 0.01$, $\epsilon = 0.01$. 当取$\delta = 0.36$时, $ K_1 =$$ 1.8 > 1$, $K_2 = 3 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta} = 3.17$, 则$E_1$渐近稳定(见图5). 当取$\delta = 0.22$时, $K_2 = 3 > 1$, $K_1 = $$ 1.1 < nK_2 = 1.296$, 则$E_2$渐近稳定(见图6).
图 5 $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.36$时, (a) S, I, T的时间序列图和(b) I-T相平面图
Figure5. When $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.36$, (a) time series graph of S, I, and T; (b) phase plan of I-T.

图 6 $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.22$时, (a) S, I, T的时间序列图和(b) I-T相平面图
Figure6. When $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.22$, (a) time series graph of S, I, and T; (b) phase plan of I-T.

图5可以看出, 随着谣言的传播, 谣言传播者的数量会达到一个峰值, 之后不断下降, 真实信息传播者数量随着时间逐渐增大, 最后达到稳定状态. 即系统(2)达到没有谣言传播者, 只有真实信息传播者和易感者的稳定状态. 从图6可以看出, 随着谣言的传播, 谣言传播者的数量逐渐增加, 真实信息传播者数量会逐渐减少变为零, 最后达到稳定状态. 即系统(2)达到没有真实信息传播者, 只有谣言传播者和易感者的稳定状态.
参数取值为$\alpha = 0.4,\; \beta = 0.7,\; A = 0.2,\; \mu = 0.2,$$ \epsilon = 0.01,\; \theta = 0.01,\; \delta = 0.234$, 且满足$\dfrac{1}{n}K_{1} < K_{2} < $$\dfrac{K_{1}(K_{1}\mu+\theta)}{\mu+\theta}$, 当$T(0)$值不相同时, 系统(2)的解会趋于不同的边界平衡点(见图7).
图 7$T(0)$不相同时, 系统(2)的解趋于不同的边界平衡点 (a) I的时间序列图; (b) 时间增加后I 的时间序列图
Figure7. When $T(0)$ is different, the solution of system (2) tends to the different boundary equilibrium: (a) Time series graph of I; (b) time series graph of I after time increment.

参数取值为$\alpha_1 = 0.7$, $\beta = 0.9$, $\delta = 0.3$, $A = $0.2, $\mu = 0.2$, $\theta = 0.01$, 则$K_2 = 3 > \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta} = $ 2.21429. 当取$\epsilon = 0.01$时, $\epsilon = 0.01 < (1-\alpha_1)\beta =$0.27, $E_1{^*}$是渐近稳定的(见图8). 当取$\epsilon = 0.3$时, $\epsilon = 0.3 > (1-\alpha_1)\beta = 0.27$, $E_2{^*}$是渐近稳定的(见图9).
图 8 $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.3 $时, (a) S, I, T的时间序列图和(b) I-T相平面图
Figure8. When $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.3 $, (a) time series graph of S, I, and T; (b) phase plan of I-T.

图 9 $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.3,\;\delta = 0.3 $时, (a) S, I, T的时间序列图和(b) I-T相平面图
Figure9. When $\alpha_1 = 0.7,\;\beta = 0.9,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.3,\;\delta = 0.3 $, (a) time series graph of S, I, and T; (b) phase plan of I-T.

$\alpha_1 \!= 0.5$, $\beta = 0.6$, $\delta \!= 0.24$, $A \!= 0.2$, $\mu \!= 0.2,$ $\theta = 0.01$, 则$K_2 = \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta} = 1.42857$, 当取$\epsilon = 0.01 < (1-2\alpha_1)\beta+ {\delta\theta}/{\mu} = 0.012$, $E_4^*$是渐近稳定的(见图10).
图 10 $ \alpha_1 = 0.5,\;\beta = 0.6,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.24 $时, (a) S, I, T的时间序列图和(b) I-T相平面图
Figure10. When $\alpha_1 = 0.5,\;\beta = 0.6,\; A = 0.2,\;\mu =0.2,\;\theta = 0.01,\;\epsilon = 0.01,\;\delta = 0.24 $, (a) time series graph of S, I, and T; (b) phase plan of I-T.

图8图10可看出, 随着谣言的传播, 易感者数量先减小后增大, 谣言传播者的数量会先增大后减小, 真实信息传播者数量逐渐增大, 最后达到稳定状态. 即系统(2)达到同时存在易感者、言传播者和真实信息传播者的稳定状态.
图11(a)可以看出, 真实信息传播者的初始值影响着谣言传播者的峰值以及谣言传播的时间, 真实信息传播者的初始值越大谣言传播者的峰值越小, 峰值到达的时间越早, 谣言灭绝得越快. 从图11(b)可以看出, 谣言传播者的初始值影响着谣言传播者的峰值以及谣言传播者峰值到达的时间, 谣言传播者的初始值越大谣言传播者的峰值越大, 峰值到达的时间越早, 但不影响谣言持续的时间.
图 11 (a)不同真实信息传播者的初始值对谣言传播的影响; (b)不同谣言传播者的初始值对谣言传播的影响; 参数取值 $\alpha_1=0.7$, $\beta=0.9$, $\delta=0.36$, $A=0.2$, $\mu=0.2$, $\theta=0.01$, $\epsilon=0.01$
Figure11. (a) Influence of different initial values of true information spreaders; (b) influence of different initial values of rumor spreaders on the rumor propagation. The parameter values are $\alpha_1=0.7$, $\beta=0.9$, $\delta=0.36$, $A=0.2$, $\mu=0.2$, $\theta=0.01$, $\epsilon=0.01$.

图11(a)得到谣言传播者峰值$I_{\rm {max}}$的散点图(见图12(a)), 通过MATLAB拟合得到函数表达式: $I_{\rm {max}} \!=\! x_1{\rm e}^{x_2T(0)}\!+\!y_1{\rm e}^{y_2T(0)}$, 其中, $x_1 \!=\! 0.2881,$ $x_2 = -20.11,$ $y_1 = 0.1296,$ $y_2 = -6.828$ (结果见图12(b)), 故T(0)与$I_{\rm {max}}$之间的函数关系为$ I_{\rm {max}} =$$ 0.2881{\rm e}^{-20.11T(0)}+0.1296{\rm e}^{-6.828T(0)}$, 可以看出谣言传播者峰值$I_{\rm {max}}$随着T(0)的增加按上述函数关系衰减.
图 12 (a) T(0)与谣言传播者峰值$I_{\rm {max}}$的散点图; (b) 根据数据点拟合得到的曲线图
Figure12. (a) Scatter plot of T(0) and the peak value $I_{\rm {max}}$ of rumor spreaders; (b) the curve graph of fitting data points.

本文在谣言传播模型中加入真实信息的传播者T, 构建了SITR谣言传播模型, 利用下一代矩阵得到了谣言传播的阈值$K_0$, 当$K_0 < 1$时, 无谣言传播者无真实信息传播者平衡点存在, 且是局部渐近稳定的. 当$K_1 > 1$时, 边界平衡点$E_1$存在, 若$K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$成立, 则边界平衡点$E_1$是局部渐近稳定的. 当$K_2 > 1$时, 边界平衡点$E_2$存在, 若$K_1 < nK_2$成立, 则边界平衡点$E_2$是局部渐近稳定的. 当$K_1 > m,$$\dfrac{1}{n}K_1 < K_2 < \dfrac{K_1(\mu K_1+\theta)}{\mu+\theta}$时, 两个边界平衡点同时稳定, 即系统出现双稳. 进一步得到了不同情形下正平衡点的存在性及其条件(见表1), 分析了正平衡点的稳定性. 最后通过数值模拟验证了理论结果, 且发现当真实信息传播者的初始值不同时, 谣言传播者的峰值是不相同的, 真实信息传播者的初始值越大谣言传播者的峰值越小, 谣言灭绝时间越短等. 因此, 在谣言传播过程中, 通过对谣言传播模型的数学分析, 可以为有效降低谣言负面效应提供理论依据. 在谣言出现时, 当传播真实信息的人越多, 谣言的危害性会越低. 我们国家和政府要通过各种途径传播谣言的危害性, 加大宣传力度, 加强对民众的教育. 政府和相关部门要随时做好应急预案, 确保当谣言出现时, 以“最真信息、最快速度、最佳效果”为目标做好工作, 减少谣言传播给社会带来的危害.
系统(2)中考虑了人口的移入和移出, 所以人群构成的网络是一个动态网络, 即网络中的节点和连边都是实时变化的, 在动态网络中, 新进入网络的个体的连边数目的分布(新个体的度分布)以及新节点与网络中已有节点的连边机制(随机连接、度优先连接、自适应连接等)等都影响着网络结构的变化, 而网络结构又影响着网络上的传播过程, 目前, 网络动态的演化对网络上传染病的传播影响已有一些研究[25-27]. 如果考虑网络结构, 需要考虑不同的新节点度分布, 和新节点与旧节点的连边机制对网络结构的影响, 然后再考虑不同情况下谣言传播的动力学行为, 这是值得深入研究的内容, 下一步我们将做这些方面的研究.
相关话题/传播 信息 系统 网络 传播者

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于时间透镜系统的冲击脉冲产生与特性研究
    摘要:冲击点火方案具备低点火能量阈值、高增益以及更好的流体力学稳定性等优势,已成为实现惯性约束聚变点火的核心方案之一.在冲击点火方案中,高质量的冲击脉冲是实现成功点火的必要条件.本文基于光纤环相位调制时间透镜系统,提出一种利用时域非对称相位调制结合频域线性色散补偿的方案产生对脉宽和峰值功率对比度高精 ...
    本站小编 Free考研考试 2021-12-29
  • 关于非均匀系统局部平均压力张量的推导及对均匀流体的分析
    摘要:由维里定理导出的适用于均匀系统的平衡态压力张量表达式可以分成两部分:动压力张量和位形压力张量.人们进而对平衡的非均匀系统进行物理分析得到了局部平均压力张量表达式.本文用更为简洁的方法推导出这一表达式.给出以原子直径为长度单位的局部平均尺寸L*>8条件下均匀流体系统平均位形压力中的三部分贡献项( ...
    本站小编 Free考研考试 2021-12-29
  • 表面等离激元传播的调制
    摘要:光学衍射极限极大地制约了纳米光子器件的发展和应用.基于表面等离激元的纳米器件由于具有突破衍射极限传播和处理光信号的特性而受到广泛关注.通过调控表面等离激元的传播,可以实现纳米尺度上光信号的调制,对片上集成纳米光子回路和光信息处理技术的发展具有重要意义.本文主要介绍了表面等离激元传播调制的基本原 ...
    本站小编 Free考研考试 2021-12-29
  • 双气泡振子系统的非线性声响应特性分析
    摘要:对初始半径不同的双气泡振子系统在声波作用下的共振行为和声响应特征进行了分析.利用微扰法分析了双泡系统的非线性共振频率,由于气泡间耦合振动的非线性影响,双泡系统存在双非线性共振频率.倍频共振和分频共振现象的存在使得双泡系统振幅-频率响应曲线有多共振峰,且随着非线性增强,共振区向低频区移动.通过对 ...
    本站小编 Free考研考试 2021-12-29
  • 矢量光共焦扫描显微系统纳米标准样品的制备与物理测量精度
    摘要:针对超分辨领域分辨率测量标准的缺失情况,本文介绍了一种用于纳米尺度分辨率测试的标准样品的设计方案和制备方法,该样品适用于矢量光共聚焦激光扫描显微系统.该设计方案包含一系列测量图案和明确的指示标记,具有测量范围宽、线宽梯级序列分布合理、制备精度高等特点.首先在非晶硅片上实现硅纳米标准样品的制备, ...
    本站小编 Free考研考试 2021-12-29
  • 光纤偏振编码量子密钥分发系统荧光边信道攻击与防御
    摘要:实际安全性是目前量子密钥分发系统中最大的挑战.在实际实现中,接收单元的单光子探测器在雪崩过程的二次光子发射(反向荧光)会导致信息泄露.目前,已有研究表明该反向荧光会泄露时间和偏振信息并且窃听行为不会在通信过程中产生额外误码率,在自由空间量子密钥分发系统中提出了利用反向荧光获取偏振信息的攻击方案 ...
    本站小编 Free考研考试 2021-12-29
  • 基于人工神经网络在线学习方法优化磁屏蔽特性参数
    摘要:磁屏蔽在磁场敏感的装置如原子钟、原子干涉仪等精密设备中发挥重要的作用,在变化外磁场下的某个磁屏蔽内部剩余磁场,可以通过Jiles-Atherton磁滞模型和磁屏蔽系数计算得出,根据计算结果可以进行主动补偿来抵消内部磁场的改变.然而实际应用中磁滞模型中五个与磁屏蔽相关的参数以及磁场衰减的两个参数 ...
    本站小编 Free考研考试 2021-12-29
  • 基于矢量像差理论的离轴反射光学系统初始结构设计
    摘要:传统的离轴反射光学系统初始结构设计方法是先求取轴对称反射光学系统结构,然后通过光瞳离轴、视场离轴或二者结合的方法实现无遮拦设计.由于同轴光学系统像差分布规律不适用于离轴光学系统,因此离轴后的反射光学系统结构像差较大,而且系统无遮拦设计过程复杂.本文提出了一种基于矢量像差理论的离轴反射光学系统初 ...
    本站小编 Free考研考试 2021-12-29
  • Fabry-Perot腔与光学参量放大复合系统中实现可调谐的非常规光子阻塞
    摘要:本文提出在Fabry-Perot腔和光学参量放大复合系统中实现非常规光子阻塞效应.此系统包含可调谐的复合型驱动强度相位,用二阶关联函数描述光子统计性质,数值模拟不同参数下的光子阻塞效应,研究发现通过调节复合型驱动强度相位可以控制非常规光子阻塞.在弱驱动条件下,计算得到了强光子反聚束的最优化条件 ...
    本站小编 Free考研考试 2021-12-29
  • 二维含多孔介质周期复合结构声传播分析
    摘要:为取得理想的隔声性能,本文结合多孔介质和周期结构两类声振调控方案,讨论了一种新型含多孔介质周期复合结构;采用等效模型描述振子系统,利用薄板理论和Biot理论建立了相应的声振耦合理论模型.利用此模型计算得出的结果与文献中数据吻合良好.研究结果表明,采用简单振子系统或组合振子系统都可以在其特征频率 ...
    本站小编 Free考研考试 2021-12-29