随机共振的概念是Benzi等^[1]在研究第四纪全球冰川期问题时提出的,揭示了噪声对非线性系统的动力学行为起到的积极有序的建设性作用的一面.近40年来, 随机共振及其相关问题的研究在理论、实验和应用方面都取得了丰硕成果^[2-6].然而,这些研究成果大多集中在经典的双稳系统和阈值系统,对于复杂非线性多稳态系统中噪声诱导共振的研究较少.许多实际的工程系统、电子电路、超导器件、控制器等,其基本的模型都是基于周期势系统建立的,比如:生物马达中的棘齿模型^[7]、物理领域的约瑟夫森结^[8-9]、工程力学中的单摆模型^[10]等.因此, 周期势系统中噪声诱导的共振研究成为非线性动力学关注的热点^[11-23].Fronzon和Mannella^[11]研究了倾斜的周期势系统中由噪声激励的布朗粒子的传输问题,并搭建了锁相环电路模型进行实验验证.Dan等^[13]从理论上研究了非均匀介质中质点在过阻尼周期势系统中的运动,取适当相位差时, 质点的运输会随着噪声强度的变化达到随机共振的效果.Zhang^[14]在恒力加白噪声驱动的欠阻尼单摆中研究了相干共振现象.Saikia等^[15-16]发现在驱动频率接近最小势阱的固有频率时,欠阻尼的周期势系统中存在随机共振现象,并把平均输入能量作为新的指标量引入到随机共振的研究中,发现两个稳定的动力学状态和随机共振的产生主要取决于系统的阻尼比和外激励的振幅.Jin等^[17-19]研究了不同噪声激励下周期势系统中的相干共振和随机共振现象,并探讨了有色噪声对系统响应和随机共振的影响.谢勇和陈若男^[21]在弱周期信号极限下,利用Floquet理论和非扰动展开法给出了一种非线性响应意义上的矩方法,对周期信号和高斯白噪声共同作用下过阻尼搓板势系统的随机共振进行了研究.Liu和Kang^[22]主要研究了非高斯列维噪声对欠阻尼周期势系统随机共振的影响,发现噪声参数对系统共振效应具有重要影响.Lou等^[23]研究了关联噪声激励下单自由度周期势系统的首次离出时间行为并讨论了系统中发生的共振激活现象.上述研究结果主要针对白噪声激励下的过阻尼系统展开,对欠阻尼周期势系统的研究较少.

噪声广泛存在于工程实际中^[24-26], 为了简化计算,常常将宽带或记忆时间很短的激励用高斯白噪声表示. 但事实上,真正的白噪声是不存在的,故实际问题中的随机激励需要用具有非零相关时间的噪声来描述.三值噪声是一种随机电报噪声, 是真实噪声的典型模型,能更好地描述自然界中的环境波动^[27-29]. 三值噪声包含了二值噪声的情形,并且在一定极限条件下能退化为高斯白噪声或者散粒白噪声. 特别地,生态系统中总存在随机环境扰动, 如, 环境(气温、湿度、光照等)的变化、气候变迁,这种扰动由于其多样性通常以三值或多值噪声形式表示. 例如:在描述N类种群的Lotka-Volterra随机系统中,需要考虑随机环境对种群数量的容纳能力的影响,该随机环境激励以三值噪声来刻画^[29]. 近年来,针对二值噪声作用下非线性系统的随机共振研究取得了一些成果^[30-35]. 如:Jin等^[30-31]研究了二值和三值噪声激励下捕食与被捕食模型的解矩稳定性,发现相比于高斯白噪声, 三值噪声更利于提高系统解的稳定性.Xu等^[32-33]利用数值方法研究了非对称二值噪声作用下双稳系统的随机共振和随机分岔,发现系统响应会呈现非对称性,且加性二值噪声的状态和强度能够诱导系统产生随机分岔现象.Fulinski^[34]研究了非马尔可夫二值噪声作用下系统中噪声诱导的迁移和随机共振现象.

本文研究了加性和乘性三值噪声及外周期信号驱动下的欠阻尼周期势系统的概率密度演化和随机共振.通过数值方法得到了系统的瞬态概率密度和联合稳态概率密度函数,分析了系统参数和噪声对概率密度演化的影响. 此外,利用系统平均输入能量来刻画随机共振,讨论了三值噪声强度、噪声转迁率、外周期信号振幅等对随机共振的影响.

考虑乘性和加性三值噪声及外周期力驱动下的欠阻尼周期势系统, 其动力学方程如下

其中, $\gamma $为阻尼系数, 势函数$U(x)=-\sin x-bx$, $b$是一个偏置常数,外周期力$F(t)=F_{0} \cos (\omega_{0} t)$, 这里$F_{0} $和$\omega_{0}$分别是其振幅和频率. $\xi (t)$和$\eta (t)$是两个独立的三值噪声,它们的统计特性为

其中, 噪声自相关时间$\tau_{i} $, 噪声强度$D_{i} (i=1,2)$.

对于三值噪声(2),假设$\xi (t)$和$\eta (t)$分别取值为$\left\{ {-A_{i} ,0,A_{i} } \right\}$ $(i=1,2)$, 从$\pm A_{i} $到0的转迁率为$\alpha_{i} $, 其他态之间的转迁率为$\beta_{i} $. 对于对称的三值噪声, $\alpha_{i} =\beta_{i} (i=1,2)$. 则有转迁概率满足^[29]

根据方程(2)和(3), 可以得到以下对应关系: $\tau_{i} =(\alpha_{i} +2\beta_{i})^{-1}$, $q_{i} ={{\beta_{i} }/{(\alpha_{i} +2\beta_{i} )}}$, $D_{i}=4q_{i} A_{i}^{2} \tau_{i} (i=1,2)$. 现以加性噪声为例, 给出三值噪声的时间历程图(见图1).

针对系统(1), 分析系统参数对随机响应的影响. 在图2和图3中,固定噪声强度$D_{1} =0.4$, $D_{2} =0.2$, $b=0.2$, $\omega_{0}={\pi /4}$, $\alpha_{1} =\alpha_{2} =10$.通过数值方法对方程式(1)进行数值积分, 所取的初值为$x_{0} \in[ -\pi :\pi /50:\pi ]$, $\dot{{x}}_{0} =0$,对不同初始值得到的响应轨线进行平均,给出系统响应的时间历程图、稳态概率密度函数. 从图2来看,当阻尼系数和外周期力振幅均较小时,系统响应在$x=1.5$附近振荡,即系统粒子仅在一个稳态势阱内运动, 其对应的稳态概率密度是单峰结构. 在图3中,当$F_{0} $增大到 0.4时, 可以看到系统粒子开始在多个稳态之间跃迁,其对应的稳态概率密度是多峰结构. 上述现象可以解释为, 对于固定的噪声强度,通过增加外周期力的振幅能够使系统在多个势阱间做跃迁运动.

下面为了研究系统的联合概率密度函数的演化,对方程(1)采用蒙特卡洛方法进行数值模拟, 令$D_{1} =0.2$, $D_{2} =0.1$,系统的其他参数选取如上.在计算区域$\{(x,\dot{{x}}):(-20,20)\times (-20,20)\}$和$\{(x,\dot{{x}}):(-2\pi ,2\pi )\times (-2\pi ,2\pi )\}$上分别取200$\times$200组初始值, 每个初始值计算组样例,计算系统平均稳态联合概率密度函数^[36-37]和瞬态概率密度函数.图4和图5显示了平均稳态联合概率密度函数在外周期力作用下的演变情况. 在图4中,对于$\gamma =0.1$, 一个周期内的平均稳态联合概率密度具有多个高度基本相等的峰. 随着$F_{0}$的增大, 峰之间的跃迁更加频繁, 多个峰的高度逐渐变得不相等,且随着$x$由$-$20增大到20, 峰值逐渐增大, 即粒子由于$F_{0} $的增加,有更大的几率运动到离稳定点更远的地方, 见图4(b).图5中固定$\gamma=0.9$, 可以看出, 此时的平均稳态联合概率密度函数呈现3个独立的尖峰结构,说明随着阻尼系数的增大, 势阱之间的跃迁运动变得困难, 符合物理直观.对比图4和图5, 可以发现阻尼系数和外周期力振幅的变化可以引起平均稳态联合概率密度函数的形态发生拓扑结构的变化,这种现象类似于随机P-分岔现象.图6显示了在一个周期内的不同时刻,系统瞬态联合概率密度函数在$(x,y)$平面内的投影, 反映了联合概率密度函数的演化过程. 可以清晰地看出,概率密度函数峰的位置之间相互连接, 其界限比较模糊, 随着时间的延长,概率密度函数峰的位置逐步演化成3个相互独立的“吸引子”,显示了瞬态概率密度函数在一个周期内的演化过程.

在外周期力和三值噪声共同作用下的周期势系统(1)中,布朗粒子存在有序的运动和无序的热运动,系统的随机共振体现在周期力做功随着噪声强度变化所表现出的非单调上^[38],该方法也称为随机能量法.

定义外周期力在一个周期$T_{0} ={{2\pi }/{\omega_{0} }}$内对系统所做的功为

将不同初始位置下的所有外加周期力输入能量式(4)进行平均, 得到系统的平均输入能量

其中$N$是计算的周期数.

根据方程式(1)、式(4)和式(5),图7$\sim \!$图11分别讨论了加性噪声强度、加性噪声转迁率和乘性噪声转迁率对系统的平均输入能量的影响.在数值计算中, 取$\gamma =0.3$, $\omega_{0} ={\pi /4}$, $N=2000$.图7和图8描述了系统分别在加性噪声激励、加性和乘性噪声共同激励下,$\bar{{W}}$随加性三值噪声强度$D_{1} $的变化规律. 在图7中, 当$F_{0}$分别取0.1, 0.2, 0.3时, $\bar{{W}}$随$D_{1} $的增大呈现非单调变化,有共振峰出现, 即出现随机共振. 随着$F_{0} $的进一步增大,$\bar{{W}}$随$D_{1} $的增大出现单调递减变化, 此时随机共振消失.类似地, 在图8中固定$F_{0} =0.3$, 系统(1)在乘性和加性噪声共同作用下,$\bar{{W}}$随$D_{1} $的增加出现共振峰,但是$\bar{{W}}$和最优的$D_{1} $的值较图7变小. 总之,当噪声强度和外周期力振幅取适当值时, 系统(1)会出现随机共振现象.

图9和图10分别给出了不同周期力振幅下,$\bar{{W}}$随加性和乘性噪声转迁率变化的曲线. 在图9中, 当$F_{0} =0.3$时,周期外力对粒子的做功较小, 粒子不能克服势垒只能做阱内运动, 当$F_{0}$分别取0.4和0.5时曲线出现了极小值, $\bar{{W}}$随着加性噪声转迁率$\alpha_{1}$的增大先减小再增大, 表明$\alpha_{1} $对系统有抑制共振的作用. 当$F_{0}=0.7$ 时曲线出现了极大值, $\bar{{W}}$随$\alpha_{1} $先增大再减小,此时系统出现了随机共振. 在图10中, 当$F_{0} =0.1$时,$\bar{{W}}$随乘性噪声转迁率$\alpha_{2} $的增大而单调递减, 但幅度非常小,未发现随机共振. 当 $F_{0} $分别取0.6和0.7时, $\bar{{W}}$作为$\alpha_{2}$的函数出现了极大值, 表明系统出现了随机共振. 与图9相比发现,在同样的周期外力作用下($F_{0} =0.7)$, $\bar{{W}}$随$\alpha_{1}$和$\alpha_{2} $的变化均出现随机共振,即加性噪声和乘性噪声转迁率对系统的影响是依赖于外周期力振幅的大小. 在图11中,同时考虑了加性和乘性噪声的影响, 发现$\bar{{W}}$作为$\alpha_{1}$的函数呈现单调变化趋势, 在此情形下, 系统(1)无随机共振现象出现.

本文主要研究了由加性和乘性三值噪声激励下,二阶欠阻尼周期势系统中的非线性随机动力学.通过计算系统的联合稳态概率密度函数和瞬态概率密度函数,讨论了系统参数和外周期力振幅对联合概率密度函数的影响,显示了瞬态概率密度函数在一个周期内的演化过程.利用随机能量法研究了周期力和噪声对系统的平均输入能量的影响,揭示了系统的随机共振现象,发现存在适当的噪声强度和外周期力振幅使得平均输入能量曲线存在极大值,导致系统出现随机共振现象.在平均输入能量随加性噪声转迁率的变化规律中发现当周期力幅值较小时,出现了抑制共振的现象;平均输入能量随乘性噪声转迁率的变化中也出现了随机共振现象,但相比于加性噪声转迁率, 未见抑制共振现象,说明乘性三值噪声和加性三值噪声对系统随机共振的影响是不同的.