Processing math: 100%

搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

三态噪声激励下分数阶耦合系统的随机共振现象

彭皓 任芮彬 钟扬帆 蔚涛

Wang Cheng-Long, Wang Qing-Yu, Zhang Yue, Li Zhong-Yu, Hong Bing, Su Zhe, Dong Liang. Molecular dynamics study of cascade damage at SiC/C interface. Acta Phys. Sin., 2014, 63(15): 153402. doi: 10.7498/aps.63.153402
Citation: Wang Cheng-Long, Wang Qing-Yu, Zhang Yue, Li Zhong-Yu, Hong Bing, Su Zhe, Dong Liang. Molecular dynamics study of cascade damage at SiC/C interface. Acta Phys. Sin., 2014, 63(15): 153402. doi: 10.7498/aps.63.153402

三态噪声激励下分数阶耦合系统的随机共振现象

彭皓, 任芮彬, 钟扬帆, 蔚涛

Stochastic resonance of fractional-order coupled system excited by trichotomous noise

Peng Hao, Ren Rui-Bin, Zhong Yang-Fan, Yu Tao
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 为了刻画在黏弹性介质中具有质量涨落的耦合粒子的运动行为, 本文提出了相应模型, 即三态噪声激励下的分数阶耦合系统. 利用Shapiro-Loginov公式和Laplace变换, 发现了粒子间的统计同步性, 并得到了系统输出幅值增益的解析表达. 在此基础上, 针对模型涉及的关键要素, 即耦合系统、分数阶系统和三态噪声, 着重分析了耦合系数、系统阶数和噪声稳态转移概率对系统输出幅值增益的广义随机共振现象的影响, 并给出了合理解释. 具体地说, 1)随着耦合系数的增大, 共振现象将先增强后减弱, 直至收敛. 该现象表明适当的耦合作用能够促进系统共振现象的产生, 体现了研究耦合系统的重要性. 2)随着系统阶数的增大, 共振现象将逐渐减弱. 当系统阶数取值为1, 即系统退化为整数阶系统时, 其输出幅值增益的峰值最小, 该现象说明分数阶系统能比传统整数阶系统得到更大的输出幅值增益. 3)噪声稳态转移概率对系统输出幅值增益的影响会随着与之相关的其他参数的变化而变化. 在一定参数条件下, 三态噪声不仅能够使系统输出幅值获得比双态噪声激励时更大的增益, 还能改变系统的共振类型. 最后, 通过数值仿真验证了上述结果的正确性.
    In order to describe the motion behavior of coupled particles with mass fluctuations in a viscous medium, we propose a corresponding model, namely a fractional-order coupled system excited by trichotomous noise. By using the Shapiro-Loginov formula and the Laplace transform, we find the statistical synchronization of the system, then obtain analytical expression of the system output amplitude gain. On this basis, this paper focuses on the key points, which are the coupled system, the fractional order system and the trichotomous noise, analyzes the influences of coupling coefficient, system order and noise steady-state probability on the generalized stochastic resonance phenomenon of system’s output amplitude gain, and gives some reasonable explanations. Specifically, first, as the coupling coefficient increases, the generalized stochastic resonance phenomenon of the output amplitude gain of the system first increases and then weakens until it converges. This phenomenon shows that the appropriate coupling strength can promote the generation of system resonance, thereby reflecting the importance of studying coupled systems. Second, with the order of the system increases, the generalized stochastic resonance phenomenon of the system’s output amplitude gain weakens gradually. When the system order value is 1, that is, when the system degenerates into an integer order system, the peak value of its output amplitude gain is smallest. This phenomenon shows that the fractional order system can obtain a larger output amplitude gain than the traditional integer order system. Third, the effect of the steady-state probability of noise on the output amplitude gain of the system changes with other related parameters. Under certain parameter conditions, trichotomous noise can not only make the output amplitude of the system larger than that of the system excited by dichotomous noise, but also change the resonance type of the system. Finally, the correctness of the above results is verified by numerical simulation.
      PACS:
      71.20.Nr(Semiconductor compounds)
      72.20.Pa(Thermoelectric and thermomagnetic effects)
      通信作者: 蔚涛, scuyutao@163.com
    • 基金项目: 国家自然科学基金青年基金(批准号: 11501385)资助的课题.
      Corresponding author: Yu Tao, scuyutao@163.com
    • Funds: Project supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant No.11501385).

    随机共振(SR)[1]是一种重要的非线性现象, 过去的几十年中得到了学者们的广泛关注. 该现象表明, 适当的噪声可以增强系统对微弱信号的响应, 颠覆了人们以往对噪声只具有破坏性的认识, 引发了相关研究热潮[2-7]. SR的发生需要3个基本条件: 1)系统的非线性性; 2)微弱相干信号; 3)噪声. 因此, 早期研究[8-13]主要集中于受加性白噪声驱动的非线性系统. 然而, 大量近期研究[14-22]表明, 受乘性色噪声驱动的线性系统也可以产生SR, 在这类系统中, 系统的非线性是由乘性噪声提供的.

    考虑黏弹性介质中的Brown运动时, 环境分子会随机吸附或脱离于Brown粒子, 使Brown粒子质量发生随机涨落. 为此, 可通过在Langevin方程中引入乘性噪声ξ(t)—“质量涨落噪声”来描述该质量涨落[23]:

    (m+ξ(t))d2x(t)dt2+γdx(t)dt+ωx(t)=A0cos(Ωt)+η(t), (1)

    其中, m表示Brown粒子的质量, ξ(t)表示粒子质量的随机涨落, x(t)表示Brown粒子的位置, γω分别代表系统的阻尼系数和固有频率, A0cos(Ωt)η(t)分别代表系统受到的外部驱动力与内部噪声. 在早期研究[24-27]中, Brown粒子的质量涨落往往被建模成双态噪声[28]. 这种做法固然能简化对模型的求解, 但却将Brown粒子的质量限制在了两种状态上. 事实上, 真实粒子的质量往往存在3种状态: 增加、减小或不变. 如红细胞可将氧气输送至人体各部的原理是: 在氧含量高的区域, 红细胞易于加载氧分子; 而在氧含量低的区域, 红细胞更容易卸载氧分子. 如果将一个红细胞及其携带的氧分子视作一个整体(Brown粒子), 则在大多数时候, 该粒子的质量变化微弱, 不足以引起其运动行为的显著改变, 此时可认为粒子的质量保持不变; 但当人体行为模式发生较大变化时, 如开始或结束剧烈运动, 该粒子的质量也将随着大量氧分子的加载或卸载而显著增加或减少. 为此, 本文考虑将粒子质量涨落建模为三态噪声[29-32]. 一方面, 如上所述, 三态噪声能更好地描述客观实际, 使得相关研究工作具有更好的实际意义; 另一方面, 在极限条件下, 三态噪声可演变为双态噪声或高斯白噪声, 从而使得该研究具有较大理论价值[33-35].

    当Brown粒子处于黏弹性介质中时, 其运动行为不仅与系统当前状态有关, 还与系统的历史状态有关. 因此, 学者们提出将Langevin方程中的阻尼项由整数阶扩展为分数阶, 利用分数阶微积分的记忆性来体现系统的历史状态对Brown粒子运动行为的影响[20,36-39]. 分数阶阻尼项的表达式如下:

    γdαx(t)dtα=γt0(tu)αΓ(1α)˙x(u)du,  0<α1, (2)

    式中, α代表分数阶的阶数, 积分核函数hD(t)=tα/[Γ(1α)]表明该分数阶阻尼项服从幂律记忆性, 这与自然界中大量具有记忆性的材料和过程的性质相符. 因此, 具有质量涨落的分数阶系统可表示为

    (m+ξ(t))d2x(t)dt2+γdαx(t)dtα+mω2x(t)=A0cos(Ωt)+ηH(t), (3)

    式中, ηH(t)是依据涨落耗散定理确定的与分数阶阻尼项相关的内噪声.

    事实上, 在黏弹性介质中往往存在着大量相互关联的Brown粒子, 其耦合作用能为系统带来丰富的动力学行为. 已有工作将对SR的研究与耦合系统相结合, 得到了大量有意义的结果[20,40-44]. 如蔚涛等[20]针对双态噪声激励下的耦合系统的SR现象进行了研究, 发现粒子间的耦合作用对系统的共振行为有着显著影响. 具体表现为, 耦合作用不仅能影响系统共振强度的大小, 还能影响系统共振的类型(单峰共振和双峰共振).

    鉴于此, 为了推广前述工作, 本文对三态噪声激励下的分数阶耦合系统的SR现象进行了分析. 需要说明的是, 本文提及的SR现象是指Berdichevsky和Gitterman[14]提出的广义随机共振现象(GSR), 即系统输出的某些函数(如矩、自相关函数和功率谱等)随着系统参数发生非单调变化的现象. 本文的结构安排如下. 第2节给出系统模型, 即三态噪声激励下的分数阶耦合系统, 并解释了相关参数的含义. 第3节利用Shapiro-Loginov公式和Laplace变换, 得到了统计同步性及系统输出幅值增益的解析表达. 第4节基于所得的解析表达, 分析了耦合系数、系统阶数和噪声稳态转移概率对系统GSR现象的影响, 并给出了合理解释. 第5节通过数值仿真验证了上述结论的正确性. 第6节是对全文得到的一系列结果的讨论与总结.

    为考察黏弹性介质中的耦合Brown运动, 引入如下带质量涨落的分数阶耦合Langevin方程, 用以刻画三态噪声激励下的分数阶耦合系统:

    {(m+ξ1)d2x1dt2+γdαx1dtα+mω2x1=ε(x2x1)+A0cos(Ωt)+ηH1(t),(4a)(m+ξ2)d2x2dt2+γdαx2dtα+mω2x2=ε(x1x2)+A0cos(Ωt)+ηH2(t),(4b)

    其中, 下标1,2用于区分不同粒子.

    将粒子的质量涨落建模为对称三态噪声ξi(i=1,2), 其幅值在a,0,a中取值, 其稳态概率为

    Ps(a)=Ps(a)=q,  Ps(0)=12q, (5)

    其中, 0<q1/2. ξi (i=1,2)满足如下统计特性:

    ξi(t)=0,ξi(t)ξi(s)=2qa2eλ|ts|, (6)

    其中, 表示取统计平均. 噪声幅值a代表Brown粒子质量增加或减少的量, 由于减少的量不能超过粒子本身的质量, 且环境分子在不同Brown粒子上的吸附和脱离行为之间没有关系, 故对于对称三态噪声, 假设噪声幅值a<m, 并且ξ1ξ2相互独立.

    对于任意固定的时刻t, ηHi(t)(i=1,2)为零均值的随机变量, 即ηHi(t)=0, 且与ξi(t)(i=1,2)产生的机制完全无关, 故其与ξi(t)(i=1,2)相互独立. 因此, 在本文后续推导中, 所有包含ηHi(t)的项均为0, 表明该加性噪声对线性系统中的共振行为不产生任何影响, 故可将其略去不表, 以达到简化演算过程的目的.

    在该模型中, xi(i=1,2)ξj(j=1,2)的泛函, 从而由Shapiro-Loginov定理[45]可知:

    {ξjdnxidtn=(ddt+λ)nξjxi,i=1,2,j=1,2,ξ2jdxidt=dξ2jxidt+λξ2jxi2qa2λxi,i=1,2,j=1,2,ξ2jd2xidt2=(ddt+λ)2ξ2jxi4qa2λdxidt2qa2λ2xi,i=1,2,j=1,2,ξ1ξ2dxidt=ξ21ξ2dxidt=ξ1ξ22dxidt=(ddt+2λ)ξ1ξ2xi,i=1,2,ξ1ξ2d2xidt2=ξ21ξ2d2xidt2=ξ1ξ22d2xidt2=(ddt+2λ)2ξ1ξ2xi,i=1,2,ξ21ξ22dxidt=(ddt+2λ)ξ2iξ22xi8q2a4λxi,i=1,2,ξ21ξ22d2xidt2=(ddt+2λ)2ξ2iξ22xi16q2a4λdxidt16q2a4λ2xi,i=1,2,ξjdαxidtα=eλtdαdtα(ξjxieλt),i=1,2,j=1,2,ξ2jdαxidtα=eλtdαdtα[(ξ2jxi2qa2xi)eλt]+2qa2dαdtαxi,i=1,2,j=1,2,ξ1ξ2dαxidtα=e2λtdαdtα(ξ1ξ2e2λt)=ξ21ξ2dαxidtα=ξ1ξ22dαxidtα,i=1,2,ξ21ξ22dαxidtα=e2λtdαdtα[(ξ21ξ22xi4q2a4xi)e2λt]+4q2a4dαdtαxi,i=1,2. (7)

    (7)式中的各公式在后续的理论推导中将发挥重要作用.

    接下来, 通过计算x1x2来分析两个粒子运动轨迹间的统计同步性.

    首先, 将方程(1a)和方程(1b)作差可得

    (md2dt2+γdαdtα+mω2+2ε)(x1x2)+(ξ1d2x1dt2ξ2d2x2dt2)=0. (8)

    对方程(8)两端同时取平均并利用方程(7)可得

    (md2dt2+γdαdtα+mω2+2ε)x1x2+(ddt+λ)2ξ1x1ξ2x2=0. (9)

    将方程(8)乘以ξ1+ξ2并取平均可得

    [m(ddt+λ)2+mω2+2ε](ξ1x1ξ2x2+ξ2x1ξ1x2)+γeλtdαdtα[(ξ1x1ξ2x2+ξ2x1ξ1x2)eλt]+(ddt+2λ)2ξ1ξ2x1ξ1ξ2x2+(ddt+λ)2ξ21x1ξ22x24qa2λ(ddt+λ2)x1x2=0. (10)

    将方程(8)乘以ξ1ξ2并取平均可得

    [m(ddt+2λ)2+mω2+2ε]ξ1ξ2x1ξ1ξ2x2+γe2λtdαdtα[(ξ1ξ2x1ξ1ξ2x2)e2λt]+(ddt+2λ)2ξ21ξ2x1ξ1ξ22x2=0. (11)

    将方程(4a)和方程(4b)分别乘以ξ1ξ2后作差并取均值可得

    [m(ddt+λ)2+mω2+ε]ξ1x1ξ2x2+γeλtdαdtα(ξ1x1ξ2x2eλt)+(ddt+λ)2ξ21x1ξ22x24qa2λ(ddt+λ2)x1x2+εξ2x1ξ1x2=0. (12)

    将方程(4a)和方程(4b)分别乘以ξ21ξ22后作差并取均值可得

    [m(ddt+λ)2+mω2+ε]ξ21x1ξ22x2+γ[eλtdαdtα(ξ21x1ξ22x22qa2x1x2)eλt+2qa2dαdtαx1x2]+a2(ddt+λ)2ξ1x1ξ2x24mqa2λ(ddt+λ2)x1x2εξ21x2ξ22x1=0. (13)

    将方程(8)乘以ξ21+ξ22后取均值可得

    [m(ddt+λ)2+mω2+2ε](ξ21x1ξ22x2+ξ22x1ξ21x2)+γ{eλtdαdtα[(ξ21x1ξ22x2+ξ22x1ξ21x24qa2x1x2)eλt]+4qa2dαdtαx1x2}8mqa2λ(ddt+λ2)x1x2+a2(ddt+λ)2ξ1x1ξ2x2+(ddt+2λ)2ξ1ξ22x1ξ21ξ2x2=0. (14)

    将方程(4a)和方程(4b)分别乘以ξ21ξ2ξ1ξ22后作差并取均值可得

    [m(ddt+2λ)2+mω2+ε]ξ21ξ2x1ξ1ξ22x2+γe2λtdαdtα(ξ21ξ2x1ξ1ξ22x2e2λt)+a2(ddt+2λ)2ξ1ξ2x1ξ1ξ2x2+εξ1ξ22x1ξ21ξ2x2=0. (15)

    将方程(4a)和方程(4b)分别乘以ξ1ξ22ξ21ξ2后作差并取均值可得

    [m(ddt+2λ)2+mω2+ε]ξ1ξ22x1ξ21ξ2x2+γe2λtdαdtα(ξ1ξ22x1ξ21ξ2x2e2λt)+(ddt+2λ)2ξ21ξ22x1ξ21ξ22x216q2a4λ(ddt+λ)x1x2+εξ21ξ2x1ξ1ξ22x2=0. (16)

    将方程(8)乘以ξ21ξ22后取均值可得

    [m(ddt+2λ)2+mω2+2ε]ξ21ξ22x1ξ21ξ22x216mq2a4λ(ddt+λ)x1x2+γ[e2λtdαdtα(ξ21ξ22x1ξ21ξ22x24q2a4x1x2)e2λt+4q2a4dαdtαx1x2]+a2(ddt+2λ)2ξ1ξ22x1ξ21ξ2x2=0. (17)

    综上, 得到了一个包含以下9个变量的封闭线性方程组(方程(9)—方程(17)):

    {z1=x1x2,z2=ξ1x1ξ2x2,z3=ξ2x1ξ1x2,z4=ξ1ξ2x1ξ1ξ2x2,z5=ξ21x1ξ22x2,z6=ξ21x2ξ22x1,z7=ξ21ξ2x1ξ1ξ22x2,z8=ξ1ξ22x1ξ21ξ2x2,z9=ξ21ξ22x1ξ21ξ22x2. (18)

    根据常微分方程中的Picard存在唯一性定理, 该封闭线性方程组有唯一解. 经过长时间演化之后, 初值的影响将逐渐消失. 因而, 不失一般性, 本文假设系统具有零初始条件. 进而, 作Laplace变换后可得

    {b11Z1+b12Z2=0,  b21Z1+b22Z2+b23Z3+b24Z4+b25Z5=0,  b34Z4+b37Z7=0,b41Z1+b42Z2+b43Z3+b45Z5=0,  b51Z1+b52Z2+b55Z5+b56Z6=0,b61Z1+b62Z2+b65Z5+b66Z6+b68Z8=0,  b74Z4+b77Z7+b78Z8=0,b81Z1+b87Z7+b88Z8+b89Z9=0,  b91Z1+b98Z8+b99Z9=0, (19)

    其中, Zi(s)=L(zi)=0zi(t)estdt,i=1,2,3,,9.

    {b11=ms2+γsα+mω2+2ε,b12=(s+λ)2,b21=4qa2λ(s+λ2),b22=m(s+λ)2+γ(s+λ)α+mω2+2ε,b23=b22,b24=(s+2λ)2,b25=b12,b34=m(s+2λ)2+γ(s+2λ)α+mω2+2ε,b37=b24,b41=b21,b42=m(s+λ)2+γ(s+λ)α+mω2+ε,b43=ε,b45=b12,b51=4mqa2λ(s+λ2)2qa2γ(s+λ)α+2qa2γsα,b52=a2(s+λ)2,b55=b42,b56=ε,b61=2b51,b62=b52,b65=b66=b22,b68=b24,b74=a2(s+2λ)2,b77=m(s+2λ)2+γ(s+2λ)α+mω2+ε,b78=b43,b81=16q2a4λ(s+λ),b87=b43,b88=b77,b89=b24,b91=16mq2a4λ(s+λ)4q2a4γ(s+2λ)α+4q2a4γsα,b98=b74,b99=b34. (20)

    求解方程组(19), 可得Zi(s)(i=1,2,3,,9)的值为

    Zi(s)=0,i=1,2,3,,9. (21)

    对(21)式作Laplace逆变换可得相应的时域解如下:

    {x1x2=0,ξ1x1ξ2x2=0,ξ2x1ξ1x2=0,ξ1ξ2x1ξ1ξ2x2=0,ξ21x1ξ22x2=0,ξ21x2ξ22x1=0,ξ21ξ2x1ξ1ξ22x2=0,ξ1ξ22x1ξ21ξ2x2=0,ξ21ξ22x1ξ21ξ22x2=0. (22)

    解(22)可进一步等价表示为

    {x1=x2,ξ1x1=ξ2x2,ξ2x1=ξ1x2,ξ1ξ2x1=ξ1ξ2x2,ξ21x1=ξ22x2,ξ21x2=ξ22x1,ξ21ξ2x1=ξ1ξ22x2,ξ1ξ22x1=ξ21ξ2x2,ξ21ξ22x1=ξ21ξ22x2. (23)

    从(23)式可以看出: 两个粒子的行为在统计意义上是完全同步的. 这表明本文的系统也同样满足蔚涛等[26]提出的“统计同步性”, 也即, 系统的平均场行为与单粒子行为具有统计一致性. 在本文的后续工作中, 即通过研究单粒子行为来研究系统的平均场行为.

    本节利用3.1节所得结果推导一阶矩xi及系统输出幅值增益G的表达式.

    对方程(4a)取平均并利用(7)式可得

    (md2dt2+γdαdtα+mω2)x1+(ddt+λ)2ξ1x1=A0cos(Ωt). (24)

    将方程(4a)乘以ξ1后取平均并利用(7)式可得

    [m(ddt+λ)2+mω2+ε]ξ1x1+γeλtdαdtα(ξ1x1eλt)+(ddt+λ)2ξ21x14qa2λ(ddt+λ2)x1+(ddt+λ)2ξ21x1εξ2x1=0. (25)

    将方程(4a)乘以ξ2后取平均并利用(7)式可得

     [m(ddt+λ)2+mω2+ε]ξ2x1+γeλtdαdtα(ξ2x1eλt)εξ1x1+(ddt+2λ)2ξ1ξ2x1=0. (26)

    将方程(4a)乘以ξ1ξ2后取平均并利用(7)式可得

    [m(ddt+2λ)2+mω2]ξ1ξ2x1+(ddt+2λ)2ξ21ξ2x1+γe2λtdαdtα(ξ1ξ2x1e2λt)=0. (27)

    将方程(4a)乘以ξ21后取平均并利用(7)式可得

    [m(ddt+λ)2+mω2+ε]ξ21x14mqa2λ(ddt+λ2)x1+γ[eλtdαdtα(ξ21x12qa2x1)eλt+2qa2dαdtαx1]+a2(ddt+λ)2ξ1x1εξ22x1=2qa2A0cos(Ωt). (28)

    将方程(4a)乘以ξ22后取平均并利用(7)式可得

    [m(ddt+λ)2+mω2+ε]ξ22x14mqa2λ(ddt+λ2)x1+γ[eλtdαdtα(ξ22x12qa2x1)eλt+2qa2dαdtαx1]+(ddt+2λ)2ξ1ξ22x1εξ21x1=2qa2A0cos(Ωt). (29)

    将方程(4a)乘以ξ21ξ2后取平均并利用(7)式可得

    [m(ddt+2λ)2+mω2+ε]ξ21ξ2x1+γe2λtdαdtα(ξ21ξ2x1e2λt)+a2(ddt+2λ)2ξ1ξ2x1εξ1ξ22x1=0. (30)

    将方程(4a)乘以ξ1ξ22后取平均并利用(7)式可得

    [m(ddt+2λ)2+mω2+ε]ξ1ξ22x1+γe2λtdαdtα(ξ1ξ22x1e2λt)+(ddt+2λ)2ξ21ξ22x1εξ21ξ2x116q2a4λ(ddt+λ)x1=0. (31)

    将方程(4a)乘以ξ21ξ22后取平均并利用(7)式可得

    [m(ddt+2λ)2+mω2]ξ21ξ22x1+a2(ddt+2λ)2ξ1ξ22x116q2a4mλ(ddt+λ)x1+γ[e2λtdαdtα(ξ21ξ22x14q2a4x1)e2λt+4q2a4dαdtαx1]=4q2a4A0cos(Ωt). (32)

    综上, 得到了一个包含以下9个变量的封闭线性方程组(方程(24)—方程(32)):

    {y1=x1,y2=ξ1x1,y3=ξ2x1,y4=ξ1ξ2x1,y5=ξ21x1,y6=ξ22x1,y7=ξ21ξ2x1,y8=ξ1ξ22x1,y9=ξ21ξ22x1. (33)

    对封闭方程组(33)作Laplace变换可得

    {c11Y1+c12Y2=A0ss2+Ω2,  c21Y1+c22Y2+c23Y3+c25Y5=0,c32Y2+c33Y3+c34Y4=0,  c44Y4+c47Y7=0,c51Y1+c52Y2+c55Y5+c56Y6=2qa2A0ss2+Ω2,c61Y1+c65Y5+c66Y6+c68Y8=2qa2A0ss2+Ω2,c74Y4+c77Y7+c78Y8=0,c81Y1+c87Y7+c88Y8+c89Y9=0,c91Y1+c98Y8+c99Y9=4q2a4A0ss2+Ω2, (34)

    其中, Yi(s)=L(yi)=0yi(t)estdt,i=1,2,3,,9.

    c11=ms2+γsα+mω2,c12=(s+λ)2,
    c21=4qa2λ(s+λ2),c22=m(s+λ)2+γ(s+λ)α+mω2+ε,c23=ε,c25=c12,
    c32=c23,c33=c22,c34=(s+2λ)2,
    c44=m(s+2λ)2+γ(s+2λ)α+mω2,c47=c34,
    c51=4mqa2λ(s+λ2)2qa2γ(s+λ)α+2qa2γsα,c52=a2(s+λ)2,c55=c22,c56=c23,
    c61=c51,c65=c23,c66=c22,c68=c34,
    c74=a2(s+2λ)2,c77=m(s+2λ)2+γ(s+2λ)α+mω2+ε,c78=c23,
    c81=16q2a4λ(s+λ),c87=c23,c88=c77,c89=c34,
    c91=16mq2a4λ(s+λ)4q2a4γ(s+2λ)α+4q2a4γsα,c98=c74,c99=c44.

    求解方程组(34)可得Yi(s)(i=1,2,3,,9)的解析表达式. 特别地,

    Y1(s)=H1(s)A0ss2+Ω2, (35)

    其中,

    H1(s)=g(s)f(s),g(s)=(2qa2c212c222+2qa2c212c22c23+c12c222c52c422+2c222c223c423)(c223c244c234c274+2c34c44c74c77c244c277)+4q2a4c12c23c234(c12c22c34c74c12c22c44c77+c222c23c34c323c34),f(s)=(c11c12c222c52c11c422+2c11c222c223c11c423c212c222c51+c212c22c23c51+c12c21c322c12c21c22c223)×(c223c244c234c274+2c34c44c74c77c244c277)+c212c22c23c34(c234c74c91c34c44c74c81c34c44c77c91+c244c77c81)+c12c223c234(c222c34c91c222c44c81c223c34c91+c223c44c81). (36)

    对(35)式作Laplace逆变换可得x1的解析表达式如下:

    x1=A0t0h1(tτ)cos(Ωτ)dτ, (37)

    其中, h1(t)H1(s)作Laplace逆变换的结果.

    另一方面, 也可通过线性时不变系统的响应理论对x1进行求解. 方程组(34)可看作一个线性时不变系统, 因而, 其输出应该是与输入信号同频的余弦信号, 仅在幅值和相位上有所差别. 据此, x1可表示为

    x1=Acos(Ωt+ϕ), (38)

    其中Aϕ分别代表x1的幅值和相位, 和(35)式对比可得

    A=A0|H1(jΩ)|,ϕ=arg[H1(jΩ)], (39)

    其中, j代表虚数单位, 满足j2=1. 最终, 根据输出幅值增益的定义可得

    G=AA0=|H1(jΩ)|, (40)

    其中, H1(s)的具体表达式见(36)式.

    特别地, 当q=0.5时, 三态噪声将退化为双态噪声. 在这种情况下, (33)式和(34)式中的变量将满足以下关系:

    y5=y6=a2y1,y7=a2y3,y8=a2y2,y9=a4y1,c21=a2λ(2s+λ),c51=a2λ(2ms+mλγ),c81=4a4λ(s+λ),c91=a4λ(4ms+4mλ2γ). (41)

    将(41)式中的关系式代入方程组(34)的前4个方程中容易发现, 方程组(34)将退化为蔚涛等[27]已发表论文中的方程组. 因而, 当q=0.5时, 本文的研究结果可涵盖蔚涛等[27]之前的研究结果.

    本节将利用(40)式对系统输出幅值增益的GSR现象进行讨论, 以分析参数对该现象的影响.

    图1为q变化的GSR在α-ε平面的相图, 反映了GSR现象出现与否的参数范围. 在本文的相图中, 不同的灰度代表系统处于不同的共振状态, 具体地说:

    图 1 随q变化的GSR在$\alpha\text{-}\varepsilon$平面的相图, 其中$ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.4 $, $ \lambda=0.1 $, $ a=0.7 $\r\nFig. 1. Phase diagram for GSR versus q in the $\alpha\text{-}\varepsilon$ plane, with $ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.4 $, $ \lambda=0.1 $, $ a=0.7 $
    图 1  q变化的GSR在α-ε平面的相图, 其中m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7
    Fig. 1.  Phase diagram for GSR versus q in the α-ε plane, with m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7

    1)黑色区域(i)代表系统输出幅值增益未出现GSR现象;

    2)灰色区域(ii)代表系统输出幅值增益出现单峰共振现象;

    3)白色区域(iii)代表系统输出幅值增益出现双峰共振现象.

    图1所示, 在这一组参数设定下, 耦合系数和系统阶数均对系统输出幅值增益的共振现象有显著影响. 当系统阶数较小时, 系统输出幅值增益出现了单峰共振现象. 而当系统阶数较大时, 系统输出幅值增益没能出现共振现象. 与之相反, 当耦合系数较小时, 系统输出幅值增益没能出现共振现象. 当耦合系数足够大时, 系统输出幅值增益出现了单峰共振现象. 为了解释上述现象, 不失一般性, 在分别选择ε=4α=0.2的条件下, 进一步分析了耦合系数和系统阶数对系统GSR现象的影响, 结果如图2所示.

    图 2 系统输出幅值增益随噪声稳态转移概率q变化的GSR现象, 其中$ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.4 $, $ \lambda=0.1 $, $ a=0.7 $ (a) $ \varepsilon=4 $; (b) $ \alpha=0.2 $\r\nFig. 2. The OAG versus q with $ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.4 $, $ \lambda=0.1 $, $ a=0.7 $: (a) $ \varepsilon=4 $; (b) $ \alpha=0.2 $.
    图 2  系统输出幅值增益随噪声稳态转移概率q变化的GSR现象, 其中m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7 (a) ε=4; (b) α=0.2
    Fig. 2.  The OAG versus q with m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7: (a) ε=4; (b) α=0.2.

    图2(a)给出了系统输出幅值增益在不同系统阶数的条件下, 随着三态噪声稳态转移概率变化而发生的变化情况. 可以看出, 当系统阶数由0.2增大到0.4时, 系统GSR现象随之减弱. 当系统阶数由0.4进一步增大到1时, 系统GSR现象逐渐消失. 该现象与图1反映的结果一致. 为了解释这个现象, 本文分析了分数阶阻尼项的频率特征, 结果如图3[27]所示. 可以看出, 随着系统阶数的增加, 输出信号中的高频成分能量增强, 低频成分能量相应减弱. 根据绝热近似原理, 低频信号较高频信号更容易达到SR现象产生的条件.

    图 3 分数阶阻尼项的频率特征\r\nFig. 3. Frequency characteristics of the fractional damping term
    图 3  分数阶阻尼项的频率特征
    Fig. 3.  Frequency characteristics of the fractional damping term

    图2(b)给出了系统输出幅值增益在不同耦合系数的条件下, 随着三态噪声稳态转移概率变化而发生的变化情况. 能够发现, 当ε=0时, 系统退化为非耦合系统, 系统输出幅值增益没能产生GSR现象. 随着耦合系数的增大, 系统输出幅值增益开始出现广义共振现象. 具体地说, 当耦合系数由0.5增大到2时, 系统的共振现象随之增强. 当耦合系数进一步由2逐渐增大到100时, 系统的共振现象逐渐减弱, 直至收敛. 这是因为粒子间适当的耦合作用能够为系统提供产生共振现象所需的非线性条件. 当耦合作用过大时, 两个粒子将被强大的耦合力连接成一个刚性的整体, 从而使系统的非线性减弱, 进而减弱系统的共振现象.

    图2还可以看出, 系统输出幅值增益的共振峰均出现在q小于0.5时, 并在q等于0.5时取得最小值. 由于三态噪声将在q等于0.5时退化为双态噪声, 上述结果说明了在这一组参数条件下, 三态噪声与双态噪声相比, 能为系统输出创造更高的增益.

    图4可以看出, 噪声稳态转移概率和系统阶数均对系统输出幅值增益随耦合系数变化的GSR行为有着显著的影响, 即随着q的增大, 系统将由双峰共振演变为单峰共振. 并且, 当噪声稳态转移概率较大时, 随着系统阶数的增大, 系统也将由双峰共振演变为单峰共振. 不失一般性, 在分别选择q=0.3α=0.1的条件下, 进一步分析了系统阶数和噪声稳态转移概率对系统GSR现象的影响, 结果如图5 所示.

    图 4 随 ε变化的GSR现象在$q\text{-}\alpha$平面的相图, 其中$ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.6 $, $ \lambda=0.1 $, $ a=0.4 $\r\nFig. 4. Phase diagram for GSR versus ε in the $q\text{-}\alpha$ plane, with $ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.6 $, $ \lambda=0.1 $, $ a=0.4 $
    图 4  ε变化的GSR现象在q-α平面的相图, 其中m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4
    Fig. 4.  Phase diagram for GSR versus ε in the q-α plane, with m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4

    图5(a)给出了系统输出幅值增益在不同系统阶数的条件下, 随耦合系数的变化情况. 不难发现, 当系统阶数取值为0.2时, 系统输出幅值增益随耦合系数的变化产生双峰共振. 随着系统阶数的增大, 系统输出幅值增益逐渐由双峰共振演变为单峰共振, 且共振峰逐渐降低. 特别地, 当系统阶数取值为1, 即系统退化为整数阶系统时, 其输出幅值增益的峰值最小, 该现象说明分数阶系统能比传统整数阶系统得到更大的输出幅值增益.

    图5(b)给出了系统输出幅值增益在不同噪声稳态转移概率的条件下, 随耦合系数变化而变化的情况. 可以看出, 当噪声稳态转移概率取值为0时, 系统输出幅值增益没能产生广义共振现象. 其原因是此时三态噪声退化为常数, 随机系统也相应地退化为确定性系统. 粒子x1x2的运动行为完全一致, 致使系统中的耦合项ε(x1x2)ε(x2x1)消失, 进而消除了系统产生共振现象所必须的非线性条件. 当噪声稳态转移概率取值为0.25时, 系统输出幅值增益出现双峰共振. 当噪声稳态转移概率进一步增大到0.5, 即三态噪声退化为双态噪声时, 系统输出幅值增益出现单峰共振, 且共振峰的高度明显提高. 该结果表明, 在这一组参数条件下, 双态噪声能使系统输出幅值获得较三态噪声激励时更大的增益.

    图 5 系统输出幅值增益随耦合系数 ε变化的GSR现象, 其中$ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.6 $, $ \lambda=0.1 $, $ a=0.4 $ (a) $q= $$ 0.3$; (b) $ \alpha=0.1 $\r\nFig. 5. The OAG versus ε with $ m=1 $, $ \gamma=0.1 $, $ \omega=1 $, $ \varOmega=1.6 $, $ \lambda=0.1 $, $ a=0.4 $: (a) $ q=0.3 $; (b) $ \alpha=0.1 $.
    图 5  系统输出幅值增益随耦合系数 ε变化的GSR现象, 其中m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4 (a) q=0.3; (b) α=0.1
    Fig. 5.  The OAG versus ε with m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4: (a) q=0.3; (b) α=0.1.

    为验证理论结果的正确性, 下面通过数值仿真模拟模型(2)式所刻画的振子运动. 对于充分小的仿真时间间隔Δt, 该模型在离散时间下的一阶近似表达式为

    {xi(tn)=yi(tn1)Δtαmemo(xi,α),yi(tn)=[γyi(tn1)mω2xi(tn1)+ε(2k=1xk(tn1)2xi(tn1))+A0cos(Ωtn1)]Δt2αm+ξi(tn)memo(yi,2α), (42)

    其中,

    {memo(xi,α)=nk=1(1)k(αk)xi(tnk),(αk)=α(α1)(α2)(αk+1)k!. (43)

    在给定系统参数条件下, 取仿真时长t=300, 采样间隔Δt=102, 可得系统输出信号时域图和频域图如图6所示.

    图 6 系统输出信号 (a) 时域图; (b) 频域图.  $m = 1$, $\gamma = 0.1$, $\alpha = 0.5$, $\varepsilon  = 1$, $\omega  = 1$, $\varOmega  = 0.4\pi$, $a  = 0.2$, $\lambda  = 0.1$, $q  = 0.3$\r\nFig. 6. System output signal: (a) Time domain diagram; (b) frequency domain diagram. $ m=1 $, $ \gamma=0.1 $, $ \alpha=0.5 $, $ \varepsilon=1 $, $ \omega=1 $, $ \varOmega=0.4\pi $, $ a=0.2 $, $ \lambda=0.1 $, $ q=0.3 $.
    图 6  系统输出信号 (a) 时域图; (b) 频域图. m=1, γ=0.1, α=0.5, ε=1, ω=1, Ω=0.4π, a=0.2, λ=0.1, q=0.3
    Fig. 6.  System output signal: (a) Time domain diagram; (b) frequency domain diagram. m=1, γ=0.1, α=0.5, ε=1, ω=1, Ω=0.4π, a=0.2, λ=0.1, q=0.3.

    图6(a)为系统输出信号时域图, 可以看出, 受噪声的影响, 系统输出具有较大的随机性. 此时, 从时域图上已经无法判断出是否有正弦响应信号存在. 因此, 通过对时域输出信号进行傅里叶变换, 得到了与之对应的频域输出信号, 如图6(b)所示. 从图6(b)可清楚地看到, 系统输出在外部驱动频率(Ω=0.4π)处出现了明显的尖峰, 表明系统输出响应中含有与外部驱动信号(A0cos(Ωt))同频的正弦信号. 其峰值(2.149)代表系统对外部驱动信号的响应幅值, 噪声的存在使得该响应幅值与其真值之间存在随机偏差. 为此, 下面采用Monte-Carlo方法来消除该随机偏差.

    图6的参数条件下, 重复进行N次仿真实验, 取N次结果的平均值作为系统稳态响应振幅的仿真值, 以(40)式确定的理论结果为参考值, 可得仿真误差随仿真次数的变化, 如图7所示.

    图 7 仿真误差图 (a) 仿真与理论对比图; (b) 仿真误差与仿真次数关系图\r\nFig. 7. Simulation error graph: (a) Comparison between simulation and theory; (b) relationship between simulation error and simulation times
    图 7  仿真误差图 (a) 仿真与理论对比图; (b) 仿真误差与仿真次数关系图
    Fig. 7.  Simulation error graph: (a) Comparison between simulation and theory; (b) relationship between simulation error and simulation times

    图7(a)为仿真与理论结果对比图. 可以看出, 仿真结果始终在理论结果的附近随机波动, 表明仿真结果与理论结果之前确实存在随机偏差. 随着仿真次数N的增加, 随机波动的幅度逐渐减小. 为更加清楚地反映该随机偏差的收敛情况, 进一步给出了相应的绝对随机偏差随仿真次数N的变化图(图7(b)). 图7(b)表明, 随着仿真次数N的增加, 噪声带来的随机偏差将逐渐收敛. 特别地, 若以δ0=103为可接受的偏差门限, 则当N2881时, 即可满足该门限需求. 一般地, 噪声强度越大或者偏差门限越小, 则所需的仿真次数N越大.

    本文研究了三态噪声激励下的分数阶耦合系统的GSR现象. 通过对模型进行一系列理论推导, 发现了粒子间的统计同步性, 并进一步得到了系统输出幅值增益的解析表达. 为了体现本文提出的模型是对前人工作的延续和推广, 针对模型涉及的关键要素, 即耦合系统、分数阶系统和三态噪声, 重点分析了耦合系数、系统阶数和噪声稳态转移概率对系统输出幅值增益的GSR现象的影响. 具体结论如下: 1)耦合系数对GSR现象的影响. 随着耦合系数的增大, GSR现象先增强后减弱, 直至收敛. 该现象表明适当的耦合作用能够促进系统共振现象的产生. 其原因是适当的耦合作用增强了系统产生GSR所需的非线性条件. 当耦合作用过大时, 两个粒子将被强大的耦合力连接成一个刚性的整体, 从而使系统的非线性减弱, 进而减弱系统的共振现象. 2)系统阶数对GSR现象的影响. 随着系统阶数的增大, GSR现象将逐渐减弱. 从分数阶阻尼项的频率特征可以看出, 随着系统阶数的增加, 输出信号中的高频成分能量增强, 低频成分能量相应减弱. 根据绝热近似原理, 低频信号较高频信号更容易达到SR现象产生的条件. 特别地, 当系统阶数取值为1, 即系统退化为整数阶系统时, 其输出幅值增益的峰值最小, 该现象说明分数阶系统能比传统整数阶系统得到更大的输出幅值增益. 此外, 不同的系统阶数还能诱导系统出现不同类型的GSR, 即单峰共振和双峰共振. 3)噪声稳态转移概率对GSR现象的影响. 噪声稳态转移概率对系统输出幅值增益的影响会随着与之相关的其他参数的变化而变化. 在一定参数条件下, 三态噪声不仅能够使系统输出幅值获得比双态噪声激励时更大的增益, 还能改变系统的共振类型. 综上所述, 耦合系数、系统阶数和噪声稳态转移概率均对系统输出幅值增益的GSR现象有显著影响, 体现了对三态噪声激励下的分数阶耦合系统的动力学行为进行理论研究的重要意义.

    [1]

    Benzi R, Sutera A, Vulpiani A 1981 J. Phys. A 14 L453Google Scholar

    [2]

    Wiesenfeld K, Moss F 1995 Nature 373 33Google Scholar

    [3]

    Gammaitoni L, Hänggi P, Jung P, Marchesoni F 1998 Rev. Mod. Phys. 70 223Google Scholar

    [4]

    Hänggi P 2002 Chem. Phys. Chem. 3 285Google Scholar

    [5]

    Zhang W, Xiang B R 2006 Talanta 70 267Google Scholar

    [6]

    Chen H, Varshney P K, Kay S M, Michels J H 2007 IEEE Trans. Signal Process. 55 3172Google Scholar

    [7]

    Soika E, Mankin R, Ainsaar A 2010 Phys. Rev. E 81 011141Google Scholar

    [8]

    Gang H, Ditzinger T, Ning C Z, Haken H 1993 Phys. Rev. Lett. 71 807Google Scholar

    [9]

    Pikovsky A, Zaikin A, Ma D L C 2002 Phys. Rev. Lett. 88 050601Google Scholar

    [10]

    Tessone C J, Mirasso C R, Toral R, Gunton J D 2006 Phys. Rev. Lett. 97 194101Google Scholar

    [11]

    冷永刚, 王太勇, 郭焱, 吴振勇 2007 物理学报 56 30Google Scholar

    Leng Y G, Wang T Y, Guo Y, Wu Z Y 2007 Acta Phys. Sin. 56 30Google Scholar

    [12]

    林敏, 黄咏梅, 方利民 2008 物理学报 57 2041Google Scholar

    Lin M, Huang Y M, Fang L M 2008 Acta Phys. Sin. 57 2041Google Scholar

    [13]

    张莉, 刘立, 曹力 2010 物理学报 59 1494Google Scholar

    Zhang L, Liu L, Cao L 2010 Acta Phys. Sin. 59 1494Google Scholar

    [14]

    Berdichevsky V, Gitterman M 1999 Phys. Rev. E 60 1494Google Scholar

    [15]

    Yu T, Zhang L, Luo M K 2013 Phys. Scr. 88 045008Google Scholar

    [16]

    Zhong S C, Ma H, Peng H, Zhang L 2015 Nonlinear Dyn. 82 535Google Scholar

    [17]

    Yang B, Zhang X, Zhang L 2016 Phys. Rev. E 94 022119

    [18]

    Lin L F, Yu L, Wang H Q 2017 Commun. Nonlinear Sci. Numer. Simul. 43 171Google Scholar

    [19]

    Zhang L, Lai L, Peng H 2018 Phys. Rev. E 97 012147Google Scholar

    [20]

    Yu T, Zhang L, Ji Y D, Lai L 2019 Commun. Nonlinear Sci. Numer. Simul. 72 26Google Scholar

    [21]

    Lin L F, Wang H Q 2019 Nonlinear Dyn. 98 801Google Scholar

    [22]

    姜磊, 赖莉, 蔚涛, 罗懋康 2021 物理学报 70 130501Google Scholar

    Jiang L, Lai L, Yu T, Luo M K 2021 Acta Phys. Sin. 70 130501Google Scholar

    [23]

    Ausloos M, Lambiotte R 2006 Phys. Rev. E 73 011105Google Scholar

    [24]

    Gitterman M, Klyatskin V I 2010 Phys. Rev. E 81 051139Google Scholar

    [25]

    Gitterman M, Shapiro I 2011 J. Stat. Phys. 144 139Google Scholar

    [26]

    蔚涛, 张路, 罗懋康 2013 物理学报 62 120504Google Scholar

    Yu T, Zhang L, Luo M K 2013 Acta Phys. Sin. 62 120504Google Scholar

    [27]

    蔚涛, 罗懋康, 华云 2013 物理学报 62 210503Google Scholar

    Yu T, Luo M K, Hua Y 2013 Acta Phys. Sin. 62 210503Google Scholar

    [28]

    Broeck C V D 1983 J. Stat. Phys. 31 467Google Scholar

    [29]

    Zhong S C, Wei K, Gao S L, Ma H 2015 J. Stat. Phys. 159 195Google Scholar

    [30]

    Mankin R, Ainsaar A, Reiter E 1999 Phys. Rev. E 60 1374Google Scholar

    [31]

    Lang R L, Yang L, Qin H L, Di G H 2012 Nonlinear Dyn. 69 1423Google Scholar

    [32]

    Peng H, Ren R B, Li P F, Yu T 2020 Phys. Scr. 95 075214Google Scholar

    [33]

    Xu P F, Jin Y F 2020 Chaos, Solitons Fractals 138 109857Google Scholar

    [34]

    Jin Y F, Wang H Q 2020 Chaos, Solitons Fractals 133 109633Google Scholar

    [35]

    Jin Y F, Niu S Y 2016 Int. J. Mod. Phys. B 30 1650102Google Scholar

    [36]

    Bagley R L, Torvik P J 1984 Appl. Mech. 51 294Google Scholar

    [37]

    Lutz L 2001 Phys. Rev. E 64 051106Google Scholar

    [38]

    He G T, Tian Y, Wang Y 2013 J. Stat. Mech. 9 26Google Scholar

    [39]

    Jiang L, Lai L, Yu T, Luo M K 2021 J. Stat. Mech. 2021 063208Google Scholar

    [40]

    Cubero D 2008 Phys. Rev. E 77 021112Google Scholar

    [41]

    Tang Y, Zou W, Lu J, Kurths J 2012 Phys. Rev. E 85 046207Google Scholar

    [42]

    Atsumi Y, Hata S, Nakao H 2013 Phys. Rev. E 2013 88Google Scholar

    [43]

    Oyarzabal R S, Szezech Jr J D, Batista A M 2017 Commun. Nonlinear Sci. Numer. Simul. 54 62Google Scholar

    [44]

    Zhang L, Xu L, Yu T, Lai L, Zhong S C 2021 Commun. Nonlinear Sci. Numer. Simul. 93 105499Google Scholar

    [45]

    Shapiro V E, Loginov V M 1978 Phys. A 91 563Google Scholar

    期刊类型引用(3)

    1. Yu-Geng Li,Yong-Qiang Liu,Mo-Ran Wang,Wen-Qing Yao,Xiao-Huan Luo,Tu Lyu,Wei-Qin Ao,Chao-Hua Zhang,Fu-Sheng Liu,Li-Peng Hu. Leveraging crystal symmetry for thermoelectric performance optimization in cubic GeSe. Rare Metals. 2024(10): 5332-5345 . 必应学术
    2. 李强,陈硕,刘可可,鲁志强,胡芹,冯利萍,张清杰,吴劲松,苏贤礼,唐新峰. n型Bi_2Te_3基化合物的类施主效应和热电性能. 物理学报. 2023(09): 135-143 . 百度学术
    3. 陈上峰,孙乃坤,张宪民,王凯,李武,韩艳,吴丽君,岱钦. Mn_3As_2掺杂Cd_3As_2纳米结构的制备及热电性能. 物理学报. 2022(18): 237-245 . 百度学术

    其他类型引用(3)

  • 图 1  q变化的GSR在α-ε平面的相图, 其中m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7

    Fig. 1.  Phase diagram for GSR versus q in the α-ε plane, with m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7

    图 2  系统输出幅值增益随噪声稳态转移概率q变化的GSR现象, 其中m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7 (a) ε=4; (b) α=0.2

    Fig. 2.  The OAG versus q with m=1, γ=0.1, ω=1, Ω=1.4, λ=0.1, a=0.7: (a) ε=4; (b) α=0.2.

    图 3  分数阶阻尼项的频率特征

    Fig. 3.  Frequency characteristics of the fractional damping term

    图 4  ε变化的GSR现象在q-α平面的相图, 其中m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4

    Fig. 4.  Phase diagram for GSR versus ε in the q-α plane, with m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4

    图 5  系统输出幅值增益随耦合系数 ε变化的GSR现象, 其中m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4 (a) q=0.3; (b) α=0.1

    Fig. 5.  The OAG versus ε with m=1, γ=0.1, ω=1, Ω=1.6, λ=0.1, a=0.4: (a) q=0.3; (b) α=0.1.

    图 6  系统输出信号 (a) 时域图; (b) 频域图. m=1, γ=0.1, α=0.5, ε=1, ω=1, Ω=0.4π, a=0.2, λ=0.1, q=0.3

    Fig. 6.  System output signal: (a) Time domain diagram; (b) frequency domain diagram. m=1, γ=0.1, α=0.5, ε=1, ω=1, Ω=0.4π, a=0.2, λ=0.1, q=0.3.

    图 7  仿真误差图 (a) 仿真与理论对比图; (b) 仿真误差与仿真次数关系图

    Fig. 7.  Simulation error graph: (a) Comparison between simulation and theory; (b) relationship between simulation error and simulation times

  • [1]

    Benzi R, Sutera A, Vulpiani A 1981 J. Phys. A 14 L453Google Scholar

    [2]

    Wiesenfeld K, Moss F 1995 Nature 373 33Google Scholar

    [3]

    Gammaitoni L, Hänggi P, Jung P, Marchesoni F 1998 Rev. Mod. Phys. 70 223Google Scholar

    [4]

    Hänggi P 2002 Chem. Phys. Chem. 3 285Google Scholar

    [5]

    Zhang W, Xiang B R 2006 Talanta 70 267Google Scholar

    [6]

    Chen H, Varshney P K, Kay S M, Michels J H 2007 IEEE Trans. Signal Process. 55 3172Google Scholar

    [7]

    Soika E, Mankin R, Ainsaar A 2010 Phys. Rev. E 81 011141Google Scholar

    [8]

    Gang H, Ditzinger T, Ning C Z, Haken H 1993 Phys. Rev. Lett. 71 807Google Scholar

    [9]

    Pikovsky A, Zaikin A, Ma D L C 2002 Phys. Rev. Lett. 88 050601Google Scholar

    [10]

    Tessone C J, Mirasso C R, Toral R, Gunton J D 2006 Phys. Rev. Lett. 97 194101Google Scholar

    [11]

    冷永刚, 王太勇, 郭焱, 吴振勇 2007 物理学报 56 30Google Scholar

    Leng Y G, Wang T Y, Guo Y, Wu Z Y 2007 Acta Phys. Sin. 56 30Google Scholar

    [12]

    林敏, 黄咏梅, 方利民 2008 物理学报 57 2041Google Scholar

    Lin M, Huang Y M, Fang L M 2008 Acta Phys. Sin. 57 2041Google Scholar

    [13]

    张莉, 刘立, 曹力 2010 物理学报 59 1494Google Scholar

    Zhang L, Liu L, Cao L 2010 Acta Phys. Sin. 59 1494Google Scholar

    [14]

    Berdichevsky V, Gitterman M 1999 Phys. Rev. E 60 1494Google Scholar

    [15]

    Yu T, Zhang L, Luo M K 2013 Phys. Scr. 88 045008Google Scholar

    [16]

    Zhong S C, Ma H, Peng H, Zhang L 2015 Nonlinear Dyn. 82 535Google Scholar

    [17]

    Yang B, Zhang X, Zhang L 2016 Phys. Rev. E 94 022119

    [18]

    Lin L F, Yu L, Wang H Q 2017 Commun. Nonlinear Sci. Numer. Simul. 43 171Google Scholar

    [19]

    Zhang L, Lai L, Peng H 2018 Phys. Rev. E 97 012147Google Scholar

    [20]

    Yu T, Zhang L, Ji Y D, Lai L 2019 Commun. Nonlinear Sci. Numer. Simul. 72 26Google Scholar

    [21]

    Lin L F, Wang H Q 2019 Nonlinear Dyn. 98 801Google Scholar

    [22]

    姜磊, 赖莉, 蔚涛, 罗懋康 2021 物理学报 70 130501Google Scholar

    Jiang L, Lai L, Yu T, Luo M K 2021 Acta Phys. Sin. 70 130501Google Scholar

    [23]

    Ausloos M, Lambiotte R 2006 Phys. Rev. E 73 011105Google Scholar

    [24]

    Gitterman M, Klyatskin V I 2010 Phys. Rev. E 81 051139Google Scholar

    [25]

    Gitterman M, Shapiro I 2011 J. Stat. Phys. 144 139Google Scholar

    [26]

    蔚涛, 张路, 罗懋康 2013 物理学报 62 120504Google Scholar

    Yu T, Zhang L, Luo M K 2013 Acta Phys. Sin. 62 120504Google Scholar

    [27]

    蔚涛, 罗懋康, 华云 2013 物理学报 62 210503Google Scholar

    Yu T, Luo M K, Hua Y 2013 Acta Phys. Sin. 62 210503Google Scholar

    [28]

    Broeck C V D 1983 J. Stat. Phys. 31 467Google Scholar

    [29]

    Zhong S C, Wei K, Gao S L, Ma H 2015 J. Stat. Phys. 159 195Google Scholar

    [30]

    Mankin R, Ainsaar A, Reiter E 1999 Phys. Rev. E 60 1374Google Scholar

    [31]

    Lang R L, Yang L, Qin H L, Di G H 2012 Nonlinear Dyn. 69 1423Google Scholar

    [32]

    Peng H, Ren R B, Li P F, Yu T 2020 Phys. Scr. 95 075214Google Scholar

    [33]

    Xu P F, Jin Y F 2020 Chaos, Solitons Fractals 138 109857Google Scholar

    [34]

    Jin Y F, Wang H Q 2020 Chaos, Solitons Fractals 133 109633Google Scholar

    [35]

    Jin Y F, Niu S Y 2016 Int. J. Mod. Phys. B 30 1650102Google Scholar

    [36]

    Bagley R L, Torvik P J 1984 Appl. Mech. 51 294Google Scholar

    [37]

    Lutz L 2001 Phys. Rev. E 64 051106Google Scholar

    [38]

    He G T, Tian Y, Wang Y 2013 J. Stat. Mech. 9 26Google Scholar

    [39]

    Jiang L, Lai L, Yu T, Luo M K 2021 J. Stat. Mech. 2021 063208Google Scholar

    [40]

    Cubero D 2008 Phys. Rev. E 77 021112Google Scholar

    [41]

    Tang Y, Zou W, Lu J, Kurths J 2012 Phys. Rev. E 85 046207Google Scholar

    [42]

    Atsumi Y, Hata S, Nakao H 2013 Phys. Rev. E 2013 88Google Scholar

    [43]

    Oyarzabal R S, Szezech Jr J D, Batista A M 2017 Commun. Nonlinear Sci. Numer. Simul. 54 62Google Scholar

    [44]

    Zhang L, Xu L, Yu T, Lai L, Zhong S C 2021 Commun. Nonlinear Sci. Numer. Simul. 93 105499Google Scholar

    [45]

    Shapiro V E, Loginov V M 1978 Phys. A 91 563Google Scholar

  • [1] 王烨花, 何美娟. 高斯色噪声激励下非对称双稳耦合网络系统的随机共振. 物理学报, 2022, 71(19): 190501. doi: 10.7498/aps.71.20220909
    [2] 彭皓, 任芮彬, 蔚涛. 三态噪声激励下分数阶耦合系统的随机共振现象研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211272
    [3] 杨建华, 马强, 吴呈锦, 刘后广. 分数阶双稳系统中的非周期振动共振. 物理学报, 2018, 67(5): 054501. doi: 10.7498/aps.67.20172046
    [4] 马正木, 靳艳飞. 二值噪声激励下欠阻尼周期势系统的随机共振. 物理学报, 2015, 64(24): 240502. doi: 10.7498/aps.64.240502
    [5] 赖志慧, 冷永刚. 三稳系统的动态响应及随机共振. 物理学报, 2015, 64(20): 200503. doi: 10.7498/aps.64.200503
    [6] 屠浙, 赖莉, 罗懋康. 分数阶非对称耦合系统在对称周期势中的定向输运. 物理学报, 2014, 63(12): 120503. doi: 10.7498/aps.63.120503
    [7] 胡建兵, 赵灵冬. 分数阶系统稳定性理论与控制研究. 物理学报, 2013, 62(24): 240504. doi: 10.7498/aps.62.240504
    [8] 贾红艳, 陈增强, 薛薇. 分数阶Lorenz系统的分析及电路实现. 物理学报, 2013, 62(14): 140503. doi: 10.7498/aps.62.140503
    [9] 田艳, 黄丽, 罗懋康. 噪声交叉关联强度的时间周期调制对线性过阻尼系统的随机共振的影响. 物理学报, 2013, 62(5): 050502. doi: 10.7498/aps.62.050502
    [10] 张静静, 靳艳飞. 非高斯噪声激励下FitzHugh-Nagumo神经元系统的随机共振. 物理学报, 2012, 61(13): 130502. doi: 10.7498/aps.61.130502
    [11] 张广丽, 吕希路, 康艳梅. 稳定噪声环境下过阻尼系统中的参数诱导随机共振现象. 物理学报, 2012, 61(4): 040501. doi: 10.7498/aps.61.040501
    [12] 张静静, 靳艳飞. 非高斯噪声驱动下非对称双稳系统的平均首次穿越时间与随机共振研究. 物理学报, 2011, 60(12): 120501. doi: 10.7498/aps.60.120501
    [13] 胡建兵, 韩焱, 赵灵冬. 分数阶系统的一种稳定性判定定理及在分数阶统一混沌系统同步中的应用. 物理学报, 2009, 58(7): 4402-4407. doi: 10.7498/aps.58.4402
    [14] 王宝华, 陆启韶, 吕淑娟. 阈下激励与噪声联合作用下肝细胞系统的内钙时空随机共振问题. 物理学报, 2009, 58(11): 7458-7465. doi: 10.7498/aps.58.7458
    [15] 宁丽娟, 徐伟. 信号调制下分段噪声驱动的线性系统的随机共振. 物理学报, 2009, 58(5): 2889-2894. doi: 10.7498/aps.58.2889
    [16] 郭立敏, 徐 伟, 阮春蕾, 赵 燕. 二值噪声驱动下二阶线性系统的随机共振. 物理学报, 2008, 57(12): 7482-7486. doi: 10.7498/aps.57.7482
    [17] 周丙常, 徐 伟. 关联噪声驱动的非对称双稳系统的随机共振. 物理学报, 2008, 57(4): 2035-2040. doi: 10.7498/aps.57.2035
    [18] 林 敏, 方利民, 朱若谷. 双频信号作用下耦合双稳系统的双共振特性. 物理学报, 2008, 57(5): 2638-2642. doi: 10.7498/aps.57.2638
    [19] 周丙常, 徐 伟. 周期混合信号和噪声联合激励下的非对称双稳系统的随机共振. 物理学报, 2007, 56(10): 5623-5628. doi: 10.7498/aps.56.5623
    [20] 徐 伟, 靳艳飞, 徐 猛, 李 伟. 偏置信号调制下色关联噪声驱动的线性系统的随机共振. 物理学报, 2005, 54(11): 5027-5033. doi: 10.7498/aps.54.5027
  • 期刊类型引用(3)

    1. Yu-Geng Li,Yong-Qiang Liu,Mo-Ran Wang,Wen-Qing Yao,Xiao-Huan Luo,Tu Lyu,Wei-Qin Ao,Chao-Hua Zhang,Fu-Sheng Liu,Li-Peng Hu. Leveraging crystal symmetry for thermoelectric performance optimization in cubic GeSe. Rare Metals. 2024(10): 5332-5345 . 必应学术
    2. 李强,陈硕,刘可可,鲁志强,胡芹,冯利萍,张清杰,吴劲松,苏贤礼,唐新峰. n型Bi_2Te_3基化合物的类施主效应和热电性能. 物理学报. 2023(09): 135-143 . 百度学术
    3. 陈上峰,孙乃坤,张宪民,王凯,李武,韩艳,吴丽君,岱钦. Mn_3As_2掺杂Cd_3As_2纳米结构的制备及热电性能. 物理学报. 2022(18): 237-245 . 百度学术

    其他类型引用(3)

计量
  • 文章访问数:  4887
  • PDF下载量:  91
  • 被引次数: 6
出版历程
  • 收稿日期:  2021-07-08
  • 修回日期:  2021-09-10
  • 上网日期:  2022-01-19
  • 刊出日期:  2022-02-05

/

返回文章
返回