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

循环噪声驱动下非对称双稳系统的驻留时间分布函数研究

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

摘要:提出了一种循环噪声驱动下非对称双稳系统驻留时间分布函数的理论计算方法. 利用具有分段逃逸速率的两态模型理论, 建立分段逃逸速率方程, 分段地推导出了驻留时间分布函数的解析表达式. 在此基础上, 从理论和数值模拟两方面阐明了在非对称性及循环噪声的影响下驻留时间分布函数呈现出反馈结构. 研究结果表明: 当非对称性、循环噪声的相关强度及循环滞后时间取适当值时均会造成驻留时间分布函数呈现出分段指数衰减现象且在循环滞后时间处出现骤然下降的趋势. 随着非对称性的减小, 驻留时间分布函数指数衰减的速率加快并出现单调指数衰减现象. 分别增大相关强度和循环滞后时间, 驻留时间分布函数在循环滞后时间处骤然下降的间距变小. 此外, 驻留时间分布函数在循环滞后时间处的值会随着噪声强度和相关强度的改变出现极大值, 说明系统发生了随机共振现象.
关键词: 循环噪声/
非对称双稳系统/
驻留时间分布函数/
随机共振

English Abstract


--> --> -->
在非线性系统中, 增加适量的噪声不但不会降低输入信号的周期性输出, 反而在一定程度上增加了系统的输出信噪比, 即输入信号与噪声之间产生了协同效应, 这就是经典的“随机共振”[1,2]现象. 随机共振的提出让人们不再认为噪声只有负面影响, 相反, 噪声在自然界中也是可以发挥积极作用的. 随后, 人们对随机共振及相关问题展开了研究并提出了量化随机共振现象的度量指标, 如信噪比[3-6]、线性响应[7-9]、驻留时间分布函数[10-16]等. 这些度量指标作为与噪声相关的函数, 随着噪声强度的改变, 若出现极大值则说明发生随机共振现象. 然而, 已有实验证明, 当把噪声注入到系统中时, 通过不同渠道的耦合噪声会发生分裂且分裂成多个含有时滞项的同源噪声, 此类噪声被称为循环噪声[17,18]. 用$\zeta \left( t \right)$表示该噪声, 考虑最简单的情形, 可表示为
$\zeta \left( t \right) = \xi \left( t \right) + \varepsilon \xi \left( {t - \tau } \right), $
其中$\xi \left( t \right)$是主噪声项, $\varepsilon \xi \left( {t - \tau } \right)$是噪声循环(第二噪声)项, $\varepsilon $表示相关强度, $\tau $表示循环滞后时间.
随着循环噪声的发现和兴起, 循环噪声驱动下系统的动力学行为研究引起了国内外****的广泛关注, 并取得了重大研究成果[19-29], 如诱导随机同步[19]、控制相干共振[20]、增强系统稳定性[22]、诱导随机分岔[28]等. 特别地, Sun等[23,24]建立了循环噪声驱动下双稳系统的两态模型理论并研究了该系统的共振动力学行为. 随后, 将该系统推广到了循环噪声驱动的时滞双稳系统上, 研究结果表明循环噪声能够调制时滞双稳系统的共振动力学行为[25]. 文献[2325]以对称双稳系统展开研究, 探讨了循环噪声对信噪比、线性响应及驻留时间分布函数结构的影响. 然而, 在许多实际的物理系统中对称性是难以保证的, 非对称系统更为普遍和具有实用价值, 因此研究非对称双稳系统的共振动力学现象应用前景更为广阔. 而如何计算循环噪声驱动下非对称性双稳系统的驻留时间分布函数以及非对称性对该系统驻留时间分布函数结构的影响目前还未被涉及.
鉴于此, 本文以循环噪声驱动下非对称双稳系统为研究对象, 运用两态模型理论推导出逃逸速率和驻留时间分布函数的解析表达式, 并基于此分析非对称性及循环噪声对驻留时间分布函数结构的影响. 本文的结构安排如下: 第二部分给出研究系统并推导出该系统逃逸速率的解析表达式; 第三部分基于逃逸速率方程推导出驻留时间分布函数的解析表达式; 第四部分讨论不同参数对驻留时间分布函数结构的影响; 第五部分总结全文.
考虑循环噪声驱动下的非对称双稳系统, 模型可用如下朗之万方程表示:
$\dot x = - V'\left( x \right) + \sqrt D \xi \left( t \right) + \varepsilon \sqrt D \xi \left( {t - \tau } \right),$

$V\left( x \right) = \frac{1}{4}x{}^4 - \frac{1}{2}x{}^2 - \gamma x,$
其中$\gamma $表示势阱的非对称性, $\xi \left( t \right)$是高斯白噪声, 其统计性质满足$\left\langle {\xi \left( t \right)} \right\rangle \!=\! 0$, $\left\langle {\xi \left( {t + t'} \right)\xi \left( t \right)} \right\rangle \!= \! 2\delta \left( {t'} \right)$. 令$\zeta \left( t \right) = \xi \left( t \right) + \varepsilon \xi \left( {t - \tau } \right)$, 则$\zeta \left( t \right)$表示循环噪声, 此时$\zeta \left( t \right)$满足
$\begin{split}&\left\langle {\zeta \left( t \right)} \right\rangle = 0,\\ & \left\langle {\zeta \left( {t + t'} \right)\zeta \left( t \right)} \right\rangle \\={}& 2\left[\left( {1 + {\varepsilon ^2}} \right)\delta \left( {t'} \right) + \varepsilon \delta \left( {t' + \tau } \right) + \varepsilon \delta \left( {t' - \tau } \right) \right].\end{split}$
上式表明$\zeta \left( t \right)$在“0”和“$ \pm \tau $”处均是自相关的, 由此可知循环噪声与高斯白噪声的统计特性不同.
当非对称双稳系统只含主噪声项(即无噪声循环项)时, 可简化为
$\dot x = x - {x^3} + \gamma + \sqrt D \xi \left( t \right).$
$ - \sqrt {\dfrac{4}{{27}}} < \gamma < \sqrt {\dfrac{4}{{27}}} $时, 令$V'\left( x \right) = 0$可得该系统的两个稳定点和一个不稳定点. 记${x_ \mp }$表示处在左右势阱的两个稳定点, ${x_{\rm{u}}}$为系统的不稳定点, 则有
${x_ - } = - 2\sqrt {\frac{1}{3}} \cos \left[ {\frac{1}{3}{\rm{arc}}\cos \left( {\gamma \sqrt {\frac{{27}}{4}} } \right) - \frac{{\text{π}}}{{\rm{3}}}} \right],$
${x_ + } = 2\sqrt {\frac{1}{3}} \cos \left[ {\frac{1}{3}{\rm{arc}}\cos \left( {\gamma \sqrt {\frac{{27}}{4}} } \right)} \right], \qquad~~$
${x_{\rm{u}}} = - 2\sqrt {\frac{1}{3}} \cos \left[ {\frac{1}{3}{\rm{arc}}\cos \left( {\gamma \sqrt {\frac{{27}}{4}} } \right) + \frac{{\text{π}}}{3}} \right].$
根据克拉默法则可得由主噪声引起的粒子逃出左右势阱内的逃逸速率为
$\begin{split}W_ \pm ^{\rm{m}}\left( t \right) ={}& \frac{{\sqrt {\left| {V''\left( {{x_{\rm{u}}}} \right)V''\left( {{x_ \mp }} \right)} \right|} }}{{2\pi }}\\ & \times \exp \left[ { - \frac{{V\left( {{x_{\rm{u}}}} \right) - V\left( {{x_ \mp }} \right)}}{D}} \right]. \end{split}$
将(2)式、(4)—(6)式代入(7)式得$W_ \mp ^{\rm{m}}\left( t \right)$满足
$\begin{split} W_ - ^{\rm{m}}\left( t \right) =\; & \frac{1}{{2{\text{π}}}}\sqrt {\left| {\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - 1} \right\}\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right] - 1} \right\}} \right|} \\ &\times \exp \left[ { - \frac{1}{D}\left( {\frac{4}{9}\left\{ {{{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right]} \right\}} \right.} \right. \\ & - \frac{2}{3}\left\{ {{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{{\rm{3}}}} \right] - {{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right]} \right\} \\ & \left. {\left. { + \frac{2}{{\sqrt 3 }}\gamma \left\{ {\cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] + \cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right]} \right\}} \right)} \right],\end{split}$
$\begin{split} W_{+}^{\rm{m}}\left( t \right) =\; & \frac{1}{{2{\text{π}}}}\sqrt {\left| {\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - 1} \right\}\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right] - 1} \right\}} \right|} \\ & \times \exp \left[ { - \frac{1}{D}\left( {\frac{4}{9}\left\{ {{{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\}} \right.} \right. \\ & - \frac{2}{3}\left\{ {{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\} \\ & \left. {\left. { + \frac{2}{{\sqrt 3 }}\gamma \left\{ {\cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - \cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{{\rm{3}}}} \right]} \right\}} \right)} \right].\end{split}$
同理, 由噪声循环项引起的粒子逃出左右势阱的逃逸速率$W_{\rm{ \pm }}^{\rm{s}}\left( t \right)$
$\begin{split} W_{+}^{\rm{s}}\left( t \right) =\; & \frac{1}{{2{\text{π}}}}\sqrt {\left| {\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - 1} \right\}\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right] - 1} \right\}} \right|} \\ & \times \exp \left[ { - \frac{1}{{D{\varepsilon ^2}}}\left( {\frac{4}{9}\left\{ {{{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\}} \right.} \right. \\ & - \frac{2}{3}\left\{ {{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\} \\ & \left. {\left. { + \frac{2}{{\sqrt 3 }}\gamma \left\{ {\cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - \cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\}} \right)} \right],\end{split} $
$\begin{split} W_ - ^{\rm{s}}\left( t \right) =\; & \frac{1}{{2{\text{π}}}}\sqrt {\left| {\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - 1} \right\}\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right] - 1} \right\}} \right|} \\ & \times \exp \left[ { - \frac{1}{{D{\varepsilon ^2}}}\left( {\frac{4}{9}\left\{ {{{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right]} \right\}} \right.} \right. \\ & - \frac{2}{3}\left\{ {{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right]} \right\} \\ & \left. {\left. { + \frac{2}{{\sqrt 3 }}\gamma \left\{ {\cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] + \cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right)} \right]} \right\}} \right)} \right].\end{split} $
假设粒子在初始时刻${t_0}$时处在右势阱中, 令${W_ \pm }\left( t \right)$表示在系统(1)中粒子逃出左右势阱的逃逸速率. 当${t_0} \leqslant t < {t_0} + \tau $时, 在$\varepsilon \ne 0$及适当的时滞$\tau $的条件下, 可假设$\xi \left( t \right)$$\varepsilon \xi \left( {t - \tau } \right)$是完全不相关的, 则粒子此时逃出左右势阱的速率分别为
$\begin{split} {W_ + }\left( t \right) =\; & \frac{1}{{2{\text{π}}}}\sqrt {\left| {\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - 1} \right\}\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right] - 1} \right\}} \right|} \\ & \times \exp \left[ { - \frac{1}{{D\left( {1 + {\varepsilon ^2}} \right)}}\left( {\frac{4}{9}\left\{ {{{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\}} \right.} \right. \\ & - \frac{2}{3}\left\{ {{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\} \\ & \left. {\left. { + \frac{2}{{\sqrt 3 }}\gamma \left\{ {\cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - \cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) - \frac{{\text{π}}}{3}} \right]} \right\}} \right)} \right],\end{split} $
$\begin{split} {W_-}(t) =\; & \frac{1}{{2{\text{π}}}}\sqrt {\left| {\left\{ {4{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - 1} \right\}\left\{ {4{{\cos }^2}\left[ \frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) \right] - 1} \right\}} \right|} \\ & \times \exp \left[ { - \frac{1}{{D\left( {1 + {\varepsilon ^2}} \right)}}\left( {\frac{4}{9}\left\{ {{{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^4}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) }\right]} \right\}} \right.} \right. \\ & - \frac{2}{3}\left\{ {{{\cos }^2}\left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] - {{\cos }^2}\left[\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) \right]} \right\} \\ & \left. {\left. { + \frac{2}{{\sqrt 3 }}\gamma \left\{ {\cos \left[ {\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) + \frac{{\text{π}}}{3}} \right] + \cos \left[\frac{1}{3}\arccos \left( {\sqrt {\frac{{27}}{4}} \gamma } \right) \right]} \right\}} \right)} \right],\end{split} $
${r_ \pm } = {W_ \pm }\left( t \right)$, 则由两态模型理论[23,24]可得粒子逃逸速率为
${r_ \pm } = {W_ \pm }\left( t \right),$
${W_ - }\left( {t + n\tau } \right) = p_ - ^{\rm{m}} - p_ + ^{\rm{s}} + \left( {p_ + ^{\rm{s}} + p_ - ^{\rm{s}}} \right){n_ + }\left[ {t + \left( {n - 1} \right)\tau } \right],$
${W_ + }\left( {t + n\tau } \right) = p_ + ^{\rm{m}} + p_ + ^{\rm{s}} - \left( {p_ + ^{\rm{s}} + p_ - ^{\rm{s}}} \right){n_ + }\left[ {t + \left( {n - 1} \right)\tau } \right],$
其中$p_ \pm ^{\rm{m}} = W_ \pm ^{\rm{m}}\left( t \right),\; p_ \pm ^{\rm{s}} = W_ \pm ^{\rm{s}}\left( t \right)$, 以及
$\begin{split} {n_ + }\left( t \right) ={}& \frac{{{r_ - }}}{{{r_ + } + {r_ - }}}\exp \left[ { - \left( {{r_ + } + {r_ - }} \right)\left( {t - {t_0}} \right)} \right] \\& + \frac{{{r_ + }}}{{{r_ + } + {r_ - }}}, \end{split}$
$\begin{split}& {n_ + }\left( {t + n\tau } \right) \\=\; & \left[ {{n_ + }\left( {{t_0} + n\tau } \right) - \frac{{p_ + ^{\rm{m}} + p_ + ^{\rm{s}}}}{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}}} \right]\\ & \times\exp \left[ { - \left( {p_ + ^{\rm{m}} + p_ - ^{\rm{m}}} \right)\left( {t - {t_0}} \right)} \right] + \frac{{p_ + ^{\rm{m}} + p_ + ^{\rm{s}}}}{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}} \\ &- \left( {p_ + ^{\rm{s}} + p_ - ^{\rm{s}}} \right)\int_{{t_0}}^t {{n_ + }\left[ {\theta + \left( {n - 1} \right)\tau } \right]}\\ & \times \exp \left[ { - \left( {p_ + ^{\rm{m}} + p_ - ^{\rm{m}}} \right)\left( {t - \theta } \right)} \right]{\rm{d}}\theta,\end{split} $
其中$0 \leqslant t < \tau $以及$n = 1, 2, \cdots $.
${p_ + }\left( t \right)$ (${p_ - }\left( t \right)$)表示粒子处在右(左)势阱的概率, 在粒子逃逸后被吸收(也就是说忽略掉双势阱间的重复跳跃事件)的条件下, 粒子的瞬时逃逸速率方程满足:
${\dot p_ \pm }\left( t \right) = - {W_ \mp }\left( t \right){p_ \pm }\left( t \right).$
首先考虑粒子逃出右(左)势阱的逃逸时间分布函数${\rho _ \pm }\left( {t\left| {{t_0}} \right.} \right)$, 其中${\rho _ \pm }\left( {t\left| {{t_0}} \right.} \right)$表示粒子在${t_0}$时刻进入了右势阱而在t时刻逃离右势阱的概率密度函数. 由(19)式可得
${\rho _ \pm }\left( {t\left| {{t_0}} \right.} \right) = - \frac{1}{Z}\frac{{{\rm{d}}{p_ \pm }\left( t \right)}}{{{\rm{d}}t}},$
其中Z为归一化常数, 满足$\displaystyle\int_0^\infty {{\rho _ \pm }\left( {t\left| {{t_0}} \right.} \right){\rm{d}}} {t_0} = 1$.
${\kern 1 pt} T = t - {t_0}$, 则${\rho _ \pm }\left( T \right)$表示粒子在右(左)势阱的驻留时间分布函数. 对逃逸时间分布函数${\rho _ \pm }\left( {t\left| {{t_0}} \right.} \right)$取平均可得到驻留时间分布函数${\rho _ \pm }\left( T \right)$, 即
${\rho _ \pm }\left( T \right) = \frac{1}{Z}\int_0^\infty {{\rho _ \pm }\left( {t\left| {{t_0}} \right.} \right){\rho _ \pm }\left( {{t_0}} \right){\rm{d}}} {t_0},$
其中$\displaystyle\int_0^\infty \!{{\rho _ \pm }\!\left( {{t_0}} \right){\rm{d}}} {t_0} = 1$, ${\rho _ \pm }\!\left( {{t_0}} \right)$为粒子在$t_0\;({t_0} \geqslant 0)$时刻进入了右(左)势阱的概率密度函数. 此外, Z表示满足$\displaystyle\int_0^\infty {{\rho _ \pm }\left( T \right){\rm{d}}} T = 1$的归一化常数.
为便于讨论, 仅考虑计算粒子在右势阱的驻留时间分布函数. 由于粒子的逃逸速率具有分段的解析式, 因此可以推导出驻留时间分布函数的递归表达式. 对于${t_0} + n\tau \leqslant t < {t_0} + \left( {n + 1} \right)\tau, n = 0, 1, \cdots $, (19)式的解为
$\begin{split}& {p_ + }\left( {t + n\tau } \right) \\={}& {p_ + }\left( {{t_0} + n\tau } \right)\exp \left[ { - \int_{{t_0}}^t {{W_ - }\left( {\theta + n\tau } \right){\rm{d}}\theta } } \right], \end{split}$
其中
$\begin{split} &{p_ + }\left( {{t_0} + n\tau } \right) \\=\; & {p_ + }\left[ {{t_0} + \left( {n - 1} \right)\tau } \right]\\ & \times\exp \left\{ { - \int_{{t_0}}^{{t_0} + \tau } {{W_ - }\left[ {\theta + \left( {n - 1} \right)\tau } \right]{\rm{d}}\theta } } \right\}.\end{split}$
将(22)式代入(20)式可得粒子逃出右势阱的逃逸时间分布函数的递归表达式为
$\begin{split} & {\rho _ + }\left( {t + n\tau \left| {{t_0}} \right.} \right) \\=\;& \frac{1}{Z}{p_ + }\left( {{t_0} + n\tau } \right){W_ - }\left( {t + n\tau } \right)\\ & \times\exp \left[ { - \int_{{t_0}}^t {{W_ - }\left( {\theta + n\tau } \right){\rm{d}}\theta } } \right]. \end{split}$
由于驻留时间分布函数与粒子的初始时刻${t_0}$无关, 不妨取${t_0} = 0$. 令${T_{{n}}} = t + \left( {n - 1} \right)\tau - {t_0}$, 则$\left( {n - 1} \right)\tau \leqslant {T_{{n}}} < n\tau, n = 1, 2, \cdots $, 通过联合(21)式、(23)式和(24)式可推导出粒子在右势阱的驻留时间分布函数的递归表达式如下:
$\begin{split} & {\rho _ + }\left( {{T_{{n}}}} \right) = \frac{1}{Z}{p_ + }\left[ {\left( {n - 1} \right)\tau } \right]{W_ - }\left( {{T_{{n}}}} \right)\\ & ~~\times\exp \left\{ { - \int_0^{{T_{{n}}} - \left( {n - 1} \right)\tau } {{W_ - }\left[ {\theta + \left( {n - 1} \right)\tau } \right]{\rm{d}}\theta } } \right\},\end{split}$
其中Z是归一化常数且满足
$\displaystyle\sum\limits_{n = 1}^\infty \int_{\left( {n - 1} \right)\tau }^{n\tau } {\rho _ + }\left( {{T_{{n}}}} \right) {\rm{d}}{T_{{n}}} = 1.$
利用数值积分, 近似地计算出$Z = 1$.
$n = 1$时, 有$0 \leqslant {T_1} < \tau $, 再将初始条件${p_ + }\left( {{t_0} = 0} \right) = 1$和(14)式代入(25)式求得第一段驻留时间分布函数的解析式为
${\rho _ + }\left( {{T_1}} \right) = {r_ - }\exp \left( { - {r_ - }{T_1}} \right).$
$n = 2$时, 有$\tau \leqslant {T_2} < 2\tau $, 然后将(14)式、(15)式、(17)式和(23)式代入(25)式得第二段驻留时间分布函数的解析式为
$\begin{split} {\rho _ + }\left( {{T_2}} \right) =\; & \exp \left( { - \frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{{4{r_ - }}} - {r_ - }\tau } \right)\left\{\frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{2}\right.\\ & \left.\times\exp \left[ { - 2{r_ - }\left( {{T_2} - \tau } \right)} \right] + \frac{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}}{2} \right\} \\ & \times \exp \left\{\frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{{4{r_ - }}}\exp \left[ { - 2{r_ - }\left( {{T_2} - \tau } \right)} \right] \right.\\ &-\left. \frac{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}}{2}\left( {{T_2} - \tau } \right) \right\}.\\[-16pt]\end{split} $
$n = 3$时, 有$2\tau \leqslant {T_3} < 3\tau $, 根据(25)式, 通过一系列计算可得第三段驻留时间分布函数的解析式为
$\begin{split}{\rho _ + }\left( {{T_3}} \right) =\; & \exp \left[ - {r_ - }\tau + \frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{{4{r_ - }}}\exp \left( { - 2{r_ - }\tau } \right) \right. \\ & \left.- \frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{{4{r_ - }}}- \frac{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}}{2}\tau\right] \times E \times F,\end{split}$
其中
$\begin{split} E =\;& \frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{2}\left[ {\exp \left( { - 2{r_ - }\tau } \right) + \frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}} - 2{r_ - }}}} \right] \\&\times\exp \left[ { - \left( {p_ + ^{\rm{m}} + p_ - ^{\rm{m}}} \right)\left( {{T_3} - 2\tau } \right)} \right] \\&- \frac{{{{\left( {p_ + ^{\rm{s}} + p_ - ^{\rm{s}}} \right)}^2}}}{{2\left( {p_ + ^{\rm{m}} + p_ - ^{\rm{m}} - 2{r_ - }} \right)}}\exp \left[ { - 2{r_ - }\left( {{T_3} - 2\tau } \right)} \right] \\&+ \frac{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}}{2},\\[-16pt]\end{split} $
$\begin{split} F \!=\; & \exp \left\{\!{\frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{{2( {p_ + ^{\rm{m}} + p_ - ^{\rm{m}}} )}}\!\left[ {\exp ({ - 2{r_ - }\tau } ) \!+\! \frac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}} \!-\! 2{r_ - }}}} \right]} \right.\\& \times\left\{\exp \left[ { - \left( {p_ + ^{\rm{m}} + p_ - ^{\rm{m}}} \right)\left( {{T_3} - 2\tau } \right)} \right]- 1\right\} \\& - \frac{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}}{2}\left( {{T_3} - 2\tau } \right) - \frac{{{{\left( {p_ + ^{\rm{s}} + p_ - ^{\rm{s}}} \right)}^2}}}{{4{r_ - }\left( {p_ + ^{\rm{m}} + p_ - ^{\rm{m}} - 2{r_ - }} \right)}}\\& \bigg.\times\left\{\exp \left[ { - 2{r_ - }\left( {{T_3} - 2\tau } \right)} \right] - 1\right\}\bigg\}.\\[-18pt]\end{split} $
按照上述的计算方法, 容易推导出$n \geqslant 4$时驻留时间分布函数的解析表达式. 理论上讲, 根据(25)式, 当n取遍所有$n \geqslant 1$的正整数时, 便可求出粒子在右势阱的驻留时间分布函数的解析表达式. 同理, 也可求得粒子在左势阱的驻留时间分布函数的解析表达式.
由于非对称性的引入, 使得系统左右势阱差发生了改变, 且循环噪声又具有特殊的统计性质, 因此对粒子的跃迁行为会产生一定的影响. 为进一步探讨在非对称性及循环噪声影响下驻留时间分布函数的结构, 针对不同参数分别从理论和数值两方面进行分析.
2
4.1.非对称性$\gamma $的影响
-->当参数$D = 0.1,\; \varepsilon = 0.5,\; \tau = 100$保持不变时, 仅考虑非对称性$\gamma $对驻留时间分布函数的影响, 如图1所示: 当$\gamma = 0.07$时, 驻留时间分布函数呈现出分段的指数衰减现象, 且在$t = \tau $处呈现骤然下降趋势, 具有明显的间断现象. 当$\gamma = 0$时, 此时系统为对称双稳系统, 驻留时间分布函数仍呈现出分段的指数衰减现象, 且在$t = \tau $处间断的间距变小. 当$\gamma = - 0.07$时, 驻留时间分布函数指数衰减速率加快且在$t = \tau $处的间断性几乎消失, 呈现出单调衰减的趋势. 此外, 随着$\gamma $的减小, 驻留时间分布函数的衰减速率随之加快, 使得粒子在右势阱的驻留时间分布函数快速地趋于零. 也就是说, 随着$\gamma $的减小, 加剧了粒子的跃迁运动. 其原因在于当$\gamma = 0.07 > 0$时, 右势阱的势阱差较大, 而当$\gamma = - 0.07 < 0$时, 右势阱的势阱差较小, 如图2所示, $\gamma = - 0.07$时更有利于粒子跃迁, 使得粒子在右势阱的驻留时间变短乃至为零.
图 1 驻留时间分布函数随非对称性$\gamma $的变化(实线是从(25)式中得到的理论结果, 圈线表示数值模拟结果, 且$D = 0.1,\; \varepsilon = 0.5,\; \tau = 100$) (a) $\gamma = 0.07$; (b) $\gamma = 0$; (c) γ = –0.07
Figure1. Variation of residence-times distribution function (RTDF) with the asymmetry (the solid line denotes the theoretical results obtained from Eq. (25), and the circle line represents the numerical simulation results), where $D = 0.1, \varepsilon = 0.5, \tau = 100$: (a) $\gamma = 0.07$; (b) $\gamma = 0$; (c) γ = –0.07

图 2 $\gamma = \pm 0.07$势阱图 (a) $\gamma = 0.07$; (b) $\gamma = - 0.07$
Figure2. Potential well of $\gamma = \pm 0.07$: (a) $\gamma = 0.07$; (b) γ = –0.07

特别地, 当$t = \tau $时, 根据(27)式可以得到驻留时间分布函数${\rho _ + }\left( \tau \right) = \left( {\dfrac{{p_ + ^{\rm{m}} + p_ - ^{\rm{m}}}}{2} + \dfrac{{p_ + ^{\rm{s}} + p_ - ^{\rm{s}}}}{2}} \right)\times \exp \left( { - {r_ - }\tau } \right)$. 为进一步说明非对称性$\gamma $对驻留时间分布函数结构的影响, 讨论${\rho _ + }\left( \tau \right)$分别作为噪声强度D、循环滞后时间$\tau $和相关强度$\varepsilon $的函数随非对称性$\gamma $的变化情况, 如图3所示. 由图3(a)可见, ${\rho _ + }\left( \tau \right)$随着D的增大先增大后减小, 故有峰值出现, 说明有共振现象发生. 随着$\gamma $的减小, ${\rho _ + }\left( \tau \right)$的峰值降低且向左偏移. 图3(b)显示, ${\rho _ + }\left( \tau \right)$随着$\tau $的增加单调衰减, 且随着$\gamma $的减小, 衰减速度加快. 在图3(c)中, ${\rho _ + }\left( \tau \right)$随着$\varepsilon $的增大先缓慢衰减后衰减加速, 且随着$\gamma $的减小, 衰减速度加快.
图 3 ${\rho _ + }\left( \tau \right)$分别作为噪声强度D、循环滞后时间$\tau $和相关强度$\varepsilon $的函数随非对称性$\gamma $的变化($\gamma = 0.1,\; \gamma = 0.04,\; \gamma = - 0.02$) (a) ε = 0.2, τ = 50; (b) D = 0.1, ε = 0.2; (c) D = 0.1, τ = 50
Figure3. Variation of ${\rho _ + }\left( \tau \right)$ which is a function of noise intensity D, relative strength $\varepsilon $ and recycling lag $\tau $ respectively with the asymmetry $\gamma $ ($\gamma \!= \! 0.1,\; \gamma = 0.04,\; \gamma \! = \! - 0.02$): (a) $\varepsilon = 0.2, \; \tau = 50$; (b) $D = 0.1, \; \varepsilon = 0.2$; (c) D = 0.1, τ = 50

2
4.2.相关强度$\varepsilon $的影响
-->固定参数$D = 0.1,\; \gamma = 0.1,\; \tau = 150$, 考虑相关强度$\varepsilon $对驻留时间分布函数结构的影响, 如图4所示: 当$\varepsilon $适当改变时, 驻留时间分布函数呈现出指数衰减现象, 且在$t = \tau $处出现间断; 随着$\varepsilon $的增大, 驻留时间分布函数指数衰减的速率加快且在$t = \tau $处间断的间距变小, 这说明了噪声循环项对粒子在双势阱的跃迁起到促进作用.
图 4 驻留时间分布函数随相关强度$\varepsilon $的变化(实线是从(25)式中得到的理论结果, 圈线表示数值模拟结果, 且$D = 0.1,\; \gamma = 0.1,\; \tau = 150$) (a) $\varepsilon = 0.2$; (b) $\varepsilon = 0.5$; (c) $\varepsilon = 0.7$
Figure4. Variation of the RTDF with the relative strength ε (the solid line denotes the theoretical results obtained from Eq. (25), and the circle represents the numerical simulation results), where $D = 0.1,\; \gamma = 0.1,\; \tau = 150$: (a) $\varepsilon = 0.2$; (b) $\varepsilon = 0.5$; (c) $\varepsilon = 0.7$.

2
4.3.循环滞后时间$\tau $的影响
-->$D = 0.09,\; \gamma = 0.08,\; \varepsilon = 0.5$保持不变时, 仅改变循环滞后时间$\tau $的值, 驻留时间分布函数呈现出分段指数衰减的现象, 且$t = \tau $时出现间断, 如图5所示. 此外, 随着$\tau $的增加, 驻留时间分布函数在$t = \tau $处的间断变得不太明显, 甚至消失, 这说明了对于较大的$\tau $($\tau \gg {1 / {{r_ - }}}$), 在循环噪声的驱动下非对称双稳系统将会进入似稳状态, 此时, 驻留时间分布函数满足${\rho _ + }\left( t \right) \propto \exp \left( { - {r_ - }t} \right)$.
图 5 驻留时间分布函数随循环滞后时间$\tau $的变化(实线是从(25)式中得到的理论结果, 圈线表示数值模拟结果, 且 $D \!= \!0.09,\; \gamma \!= \! 0.08,\; \varepsilon = 0.5$) (a) $\tau \!=\! 150$; (b) $\tau \!=\! 200$; (c) $\tau = 300$
Figure5. Variation of the RTDF with the recycling lag $\tau $ (the solid line denotes the theoretical results obtained from Eq. (25), and the circle represents the numerical simulation results), where $D = 0.09, \;\gamma = 0.08,\; \varepsilon = 0.5$: (a) $\tau = 150$; (b) $\tau = 200$; (c) $\tau = 300$.

为进一步研究循环噪声驱动下非对称双稳系统的共振行为, 采用驻留时间分布函数在$t = \tau $处的值${\rho _ + }\left( \tau \right)$进行具体分析. 如图6所示, 随着噪声强度D和相关强度$\varepsilon $的改变, ${\rho _ + }\left( \tau \right)$出现了极大值, 这说明了在循环噪声的作用下非对称双稳系统发生了随机共振现象.
图 6 $ {\rho _ + }\left( \tau \right)$D$ \varepsilon $的变化, 其中$ \gamma = 0.05, \;\tau = 150$, $ 0.1 \leqslant \varepsilon \leqslant 1$, $ 0.03 \leqslant D \leqslant 0.3$
Figure6. Variation of ${\rho _ + }\left( \tau \right)$ with D and $\varepsilon $, where $\gamma = 0.05, \tau = 150$, $0.1 \leqslant \varepsilon \leqslant 1$, $0.03 \leqslant D \leqslant 0.3$.

本文主要研究了循环噪声驱动下非对称双稳系统驻留时间分布函数的计算方法, 数值模拟结果验证了该方法的可行性. 利用具有分段逃逸速率的两态模型理论, 推导出了驻留时间分布函数的分段解析式. 基于此, 分别从理论和数值模拟的角度分析了不同参数对驻留时间分布函数结构产生的影响, 着重说明了在非对称性$\gamma $的影响下驻留时间分布函数产生的反馈结构: 当$\gamma > 0$且取适当值时, 驻留时间分布函数呈现出分段指数衰减现象, 且当$t = \tau $时, 表现出明显的间断现象; 当$\gamma < 0$时, 驻留时间分布函数在$t = \tau $处的间断性几乎消失且指数衰减速率加快, 呈现出单调衰减现象. 以驻留时间分布函数在$t = \tau $处的值进行具体分析, 进一步说明随着非对称性的减小, 驻留时间分布函数的衰减速度加快. 此外, 驻留时间分布函数在$t = \tau $处的值随着噪声强度和相关强度的改变出现极大值, 说明在非对称双稳系统中, 循环噪声能够诱导随机共振现象的发生. 本文研究方法考虑了系统的非对称性, 更具有实用价值, 该方法可用于计算循环噪声驱动下时滞非对称双稳系统的驻留时间分布函数.
相关话题/系统 结构 计算 函数 噪声

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 超强耦合电路量子电动力学系统中反旋波效应对量子比特频率移动的影响
    摘要:从实验上研究了四结磁通量子比特与多模共面波导谐振腔构成的超强耦合电路量子电动力学系统.通过传输谱测量和数值拟合,确定量子比特与腔第一模式的耦合强度已达到0.1倍谐振腔频率,进入了超强耦合区域;通过色散读出方法得到了系统的能谱,并通过增加探测场光子,从能谱上得到了量子比特频率随探测光子的位移.这 ...
    本站小编 Free考研考试 2021-12-29
  • 基于希尔伯特变换的结构光照明快速三维彩色显微成像方法
    摘要:结构光照明显微是一种宽场显微技术,可以实现超分辨成像和三维光切片成像.基于HSV(色相、饱和度、明度)彩色空间的结构光照明全彩色三维光切片成像技术可以复原样品表面的真彩色信息,但每一层光切片都需要采集3幅固定相移差的原始图像,这对于需要多视场拼接的大尺寸样品而言,图像采集数据量大、图像重构时间 ...
    本站小编 Free考研考试 2021-12-29
  • 马约拉纳零能模的非阿贝尔统计及其在拓扑量子计算的应用
    摘要:自1937年被预言以来,马约拉纳费米子在粒子物理领域和暗物质领域就广受关注.它们在凝聚态物理中的“副本”,马约拉纳零能模(Majoranazeromode,MZM),被指出可以通过拓扑超导实现,并由于满足非阿贝尔统计及可以用来实现容错的量子计算机而成为凝聚态领域最受关注的研究方向之一.尤其在近 ...
    本站小编 Free考研考试 2021-12-29
  • 三模腔-原子闭环系统中可控的量子干涉和光子传输
    摘要:通过构造一个由相互垂直的两腔和一个二能级原子组成的光学腔-原子系统,研究可控的量子干涉引起的非传统光子传输现象.该系统中,两个正交腔之间通过光纤直接耦合和通过放在两腔交叉处的二能级原子间接耦合.该三模系统支持两个相互垂直的传播方向,即两探测场相互垂直.在考虑原子弛豫速率的情况下,该闭环系统中的 ...
    本站小编 Free考研考试 2021-12-29
  • 两种Ge-Sb-Se薄膜的光学性质及微观结构
    摘要:提出一种综合利用区域逼近法和柯西拟合法精确获取Ge20Sb15Se65薄膜和Ge28Sb12Se60薄膜透射光谱范围内任意波长处折射率与色散的多点柯西法,并从理论上证明了该方法的准确性.实验上,采用磁控溅射法制备了这两种Ge—Sb—Se薄膜,利用傅里叶红外光谱仪测得了透射光谱曲线,运用分段滤波 ...
    本站小编 Free考研考试 2021-12-29
  • 深海海底反射区声场角谱域分布结构分析及在声纳波束俯仰上的应用
    摘要:深海海底反射区的声场干涉导致能量起伏,存在不连续的若干声纳可探测区.主动声纳探测海底反射区目标时,必须建立起声纳可探测区与波束俯仰角间的量化关系,通过合理选择最优发射波束俯仰角,才能使其对准声纳可探测区.本文通过理论分析和数值仿真,指出海底反射区离散的声纳可探测区的形成与不同掠射角声线能量周期 ...
    本站小编 Free考研考试 2021-12-29
  • 分子动力学模拟冷却速率对非晶合金结构与变形行为的影响
    摘要:非晶合金因具有独特的无序结构、优异或独特的各种性能以及良好的应用前景,而受到专家****的广泛关注.其中,制备过程中的冷却速率对非晶的结构与性能起着非常重要的调控作用.本文采用分子动力学的模拟方法,分别以4种冷却速率获得相同尺寸的Zr48Cu45Al7三元非晶合金的制备态原子结构模型,并模拟了 ...
    本站小编 Free考研考试 2021-12-29
  • 受激辐射损耗超分辨显微成像系统研究的新进展
    摘要:由于受到衍射极限的影响,传统光学显微镜的分辨率被限制在半个波长左右.近二十年来出现了许多通过不同方法绕过光学衍射极限的超分辨成像技术,其中,受激辐射损耗显微(stimulatedemissiondepletionmicroscopy,STED)通过引入一束环形损耗光来抑制荧光光斑外围荧光分子的 ...
    本站小编 Free考研考试 2021-12-29
  • Re对Ni<sub>3</sub>Al微观结构及力学性质影响的第一原理研究
    摘要:应用基于密度泛函理论和广义梯度近似的第一原理方法,探究了Re元素掺杂镍铝合金中${\gamma^{\prime}}$相Ni3Al之后微观结构和力学性质的变化.结果表明,在大部分化学计量比范围内,Re原子在Ni3Al中易于替代Al位.Re的掺杂引起Ni3Al晶格常数小幅度地增大,且不会引起严重的 ...
    本站小编 Free考研考试 2021-12-29
  • 磁电势垒结构中光场辅助电子自旋输运特性
    摘要:基于Floquet理论和传输矩阵方法,理论研究了光场对电子隧穿两类磁电垒结构的自旋极化输运特性的影响,计算结果表明光场对两类磁电垒结构中电子的输运有显著影响:首先,原来不存在自旋过滤特性的结构应用光场后会产生低能区域明显的自旋过滤效应;其次,原来存在自旋过滤特性的结构应用光场后自旋过滤明显增强 ...
    本站小编 Free考研考试 2021-12-29