-
为了刻画在黏弹性介质中具有质量涨落的耦合粒子的运动行为, 本文提出了相应模型, 即三态噪声激励下的分数阶耦合系统. 利用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.
-
Keywords:
- stochastic resonance /
- fractional system /
- coupled system /
- trichotomous noise
1. 引 言
随机共振(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(t−u)−αΓ(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节是对全文得到的一系列结果的讨论与总结.
2. 模 型
为考察黏弹性介质中的耦合Brown运动, 引入如下带质量涨落的分数阶耦合Langevin方程, 用以刻画三态噪声激励下的分数阶耦合系统:
{(m+ξ1)d2x1dt2+γdαx1dtα+mω2x1=ε(x2−x1)+A0cos(Ωt)+ηH1(t),(4a)(m+ξ2)d2x2dt2+γdαx2dtα+mω2x2=ε(x1−x2)+A0cos(Ωt)+ηH2(t),(4b) 其中, 下标
1,2 用于区分不同粒子.将粒子的质量涨落建模为对称三态噪声
ξi (i=1,2) , 其幅值在−a,0,a 中取值, 其稳态概率为Ps(−a)=Ps(a)=q, Ps(0)=1−2q, (5) 其中,
0<q⩽1/2 .ξi (i=1,2) 满足如下统计特性:⟨ξi(t)⟩=0,⟨ξi(t)ξi(s)⟩=2qa2e−λ|t−s|, (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⟨ξ2jxi⟩dt+λ⟨ξ2jxi⟩−2qa2λ⟨xi⟩,i=1,2,j=1,2,⟨ξ2jd2xidt2⟩=(ddt+λ)2⟨ξ2jxi⟩−4qa2λd⟨xi⟩dt−2qa2λ2⟨xi⟩,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ξ22xi⟩−8q2a4λ⟨xi⟩,i=1,2,⟨ξ21ξ22d2xidt2⟩=(ddt+2λ)2⟨ξ2iξ22xi⟩−16q2a4λd⟨xi⟩dt−16q2a4λ2⟨xi⟩,i=1,2,⟨ξjdαxidtα⟩=e−λtdαdtα(⟨ξjxi⟩eλt),i=1,2,j=1,2,⟨ξ2jdαxidtα⟩=e−λtdαdtα[(⟨ξ2jxi⟩−2qa2⟨xi⟩)eλt]+2qa2dαdtα⟨xi⟩,i=1,2,j=1,2,⟨ξ1ξ2dαxidtα⟩=e−2λtdαdtα(⟨ξ1ξ2⟩e2λt)=⟨ξ21ξ2dαxidtα⟩=⟨ξ1ξ22dαxidtα⟩,i=1,2,⟨ξ21ξ22dαxidtα⟩=e−2λtdαdtα[(⟨ξ21ξ22xi⟩−4q2a4⟨xi⟩)e2λt]+4q2a4dαdtα⟨xi⟩,i=1,2. (7) (7)式中的各公式在后续的理论推导中将发挥重要作用.
3. 理论结果
3.1 同步性
接下来, 通过计算
⟨x1−x2⟩ 来分析两个粒子运动轨迹间的统计同步性.(md2dt2+γdαdtα+mω2+2ε)(x1−x2)+(ξ1d2x1dt2−ξ2d2x2dt2)=0. (8) (md2dt2+γdαdtα+mω2+2ε)⟨x1−x2⟩+(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−ξ22x2⟩−4qa2λ(ddt+λ2)⟨x1−x2⟩=0. (10) 将方程(8)乘以
ξ1ξ2 并取平均可得[m(ddt+2λ)2+mω2+2ε]⟨ξ1ξ2x1−ξ1ξ2x2⟩+γe−2λ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−ξ2x2⟩eλt)+(ddt+λ)2⟨ξ21x1−ξ22x2⟩−4qa2λ(ddt+λ2)⟨x1−x2⟩+ε⟨ξ2x1−ξ1x2⟩=0. (12) 将方程(4a)和方程(4b)分别乘以
ξ21 与ξ22 后作差并取均值可得[m(ddt+λ)2+mω2+ε]⟨ξ21x1−ξ22x2⟩+γ[e−λtdαdtα(⟨ξ21x1−ξ22x2⟩−2qa2⟨x1−x2⟩)eλt+2qa2dαdtα⟨x1−x2⟩]+a2(ddt+λ)2⟨ξ1x1−ξ2x2⟩−4mqa2λ(ddt+λ2)⟨x1−x2⟩−ε⟨ξ21x2−ξ22x1⟩=0. (13) 将方程(8)乘以
ξ21+ξ22 后取均值可得[m(ddt+λ)2+mω2+2ε](⟨ξ21x1−ξ22x2⟩+⟨ξ22x1−ξ21x2⟩)+γ{e−λtdαdtα[(⟨ξ21x1−ξ22x2⟩+⟨ξ22x1−ξ21x2⟩−4qa2⟨x1−x2⟩)eλt]+4qa2dαdtα⟨x1−x2⟩}−8mqa2λ(ddt+λ2)⟨x1−x2⟩+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⟩+γe−2λtdαdtα(⟨ξ21ξ2x1−ξ1ξ22x2⟩e2λ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⟩+γe−2λtdαdtα(⟨ξ1ξ22x1−ξ21ξ2x2⟩e2λt)+(ddt+2λ)2⟨ξ21ξ22x1−ξ21ξ22x2⟩−16q2a4λ(ddt+λ)⟨x1−x2⟩+ε⟨ξ21ξ2x1−ξ1ξ22x2⟩=0. (16) 将方程(8)乘以
ξ21ξ22 后取均值可得[m(ddt+2λ)2+mω2+2ε]⟨ξ21ξ22x1−ξ21ξ22x2⟩−16mq2a4λ(ddt+λ)⟨x1−x2⟩+γ[e−2λtdαdtα(⟨ξ21ξ22x1−ξ21ξ22x2⟩−4q2a4⟨x1−x2⟩)e2λt+4q2a4dαdtα⟨x1−x2⟩]+a2(ddt+2λ)2⟨ξ1ξ22x1−ξ21ξ2x2⟩=0. (17) 综上, 得到了一个包含以下9个变量的封闭线性方程组(方程(9)—方程(17)):
{z1=⟨x1−x2⟩,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)e−stdt,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逆变换可得相应的时域解如下:
{⟨x1−x2⟩=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.2 输出幅值增益
本节利用3.1节所得结果推导一阶矩
⟨xi⟩ 及系统输出幅值增益G的表达式.(md2dt2+γdαdtα+mω2)⟨x1⟩+(ddt+λ)2⟨ξ1x1⟩=A0cos(Ωt). (24) [m(ddt+λ)2+mω2+ε]⟨ξ1x1⟩+γe−λtdαdtα(⟨ξ1x1⟩eλt)+(ddt+λ)2⟨ξ21x1⟩−4qa2λ(ddt+λ2)⟨x1⟩+(ddt+λ)2⟨ξ21x1⟩−ε⟨ξ2x1⟩=0. (25) [m(ddt+λ)2+mω2+ε]⟨ξ2x1⟩+γe−λtdαdtα(⟨ξ2x1⟩eλt)−ε⟨ξ1x1⟩+(ddt+2λ)2⟨ξ1ξ2x1⟩=0. (26) [m(ddt+2λ)2+mω2]⟨ξ1ξ2x1⟩+(ddt+2λ)2⟨ξ21ξ2x1⟩+γe−2λtdαdtα(⟨ξ1ξ2x1⟩e2λt)=0. (27) [m(ddt+λ)2+mω2+ε]⟨ξ21x1⟩−4mqa2λ(ddt+λ2)⟨x1⟩+γ[e−λtdαdtα(⟨ξ21x1⟩−2qa2⟨x1⟩)eλt+2qa2dαdtα⟨x1⟩]+a2(ddt+λ)2⟨ξ1x1⟩−ε⟨ξ22x1⟩=2qa2A0cos(Ωt). (28) [m(ddt+λ)2+mω2+ε]⟨ξ22x1⟩−4mqa2λ(ddt+λ2)⟨x1⟩+γ[e−λtdαdtα(⟨ξ22x1⟩−2qa2⟨x1⟩)eλt+2qa2dαdtα⟨x1⟩]+(ddt+2λ)2⟨ξ1ξ22x1⟩−ε⟨ξ21x1⟩=2qa2A0cos(Ωt). (29) [m(ddt+2λ)2+mω2+ε]⟨ξ21ξ2x1⟩+γe−2λtdαdtα(⟨ξ21ξ2x1⟩e2λt)+a2(ddt+2λ)2⟨ξ1ξ2x1⟩−ε⟨ξ1ξ22x1⟩=0. (30) [m(ddt+2λ)2+mω2+ε]⟨ξ1ξ22x1⟩+γe−2λtdαdtα(⟨ξ1ξ22x1⟩e2λt)+(ddt+2λ)2⟨ξ21ξ22x1⟩−ε⟨ξ21ξ2x1⟩−16q2a4λ(ddt+λ)⟨x1⟩=0. (31) [m(ddt+2λ)2+mω2]⟨ξ21ξ22x1⟩+a2(ddt+2λ)2⟨ξ1ξ22x1⟩−16q2a4mλ(ddt+λ)⟨x1⟩+γ[e−2λtdαdtα(⟨ξ21ξ22x1⟩−4q2a4⟨x1⟩)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)e−stdt,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+c12c222c52−c422+2c222c223−c423)(c223c244−c234c274+2c34c44c74c77−c244c277)+4q2a4c12c23c234(c12c22c34c74−c12c22c44c77+c222c23c34−c323c34),f(s)=(c11c12c222c52−c11c422+2c11c222c223−c11c423−c212c222c51+c212c22c23c51+c12c21c322−c12c21c22c223)×(c223c244−c234c274+2c34c44c74c77−c244c277)+c212c22c23c34(c234c74c91−c34c44c74c81−c34c44c77c91+c244c77c81)+c12c223c234(c222c34c91−c222c44c81−c223c34c91+c223c44c81). (36) 对(35)式作Laplace逆变换可得
⟨x1⟩ 的解析表达式如下:⟨x1⟩=A0∫t0h1(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]之前的研究结果.4. 系统输出幅值增益的GSR现象
本节将利用(40)式对系统输出幅值增益的GSR现象进行讨论, 以分析参数对该现象的影响.
4.1 系统输出幅值增益随噪声稳态转移概率q变化的GSR
图1为随q变化的GSR在
α-ε 平面的相图, 反映了GSR现象出现与否的参数范围. 在本文的相图中, 不同的灰度代表系统处于不同的共振状态, 具体地说:1)黑色区域(i)代表系统输出幅值增益未出现GSR现象;
2)灰色区域(ii)代表系统输出幅值增益出现单峰共振现象;
3)白色区域(iii)代表系统输出幅值增益出现双峰共振现象.
如图1所示, 在这一组参数设定下, 耦合系数和系统阶数均对系统输出幅值增益的共振现象有显著影响. 当系统阶数较小时, 系统输出幅值增益出现了单峰共振现象. 而当系统阶数较大时, 系统输出幅值增益没能出现共振现象. 与之相反, 当耦合系数较小时, 系统输出幅值增益没能出现共振现象. 当耦合系数足够大时, 系统输出幅值增益出现了单峰共振现象. 为了解释上述现象, 不失一般性, 在分别选择
ε=4 和α=0.2 的条件下, 进一步分析了耦合系数和系统阶数对系统GSR现象的影响, 结果如图2所示.图2(a)给出了系统输出幅值增益在不同系统阶数的条件下, 随着三态噪声稳态转移概率变化而发生的变化情况. 可以看出, 当系统阶数由0.2增大到0.4时, 系统GSR现象随之减弱. 当系统阶数由0.4进一步增大到1时, 系统GSR现象逐渐消失. 该现象与图1反映的结果一致. 为了解释这个现象, 本文分析了分数阶阻尼项的频率特征, 结果如图3[27]所示. 可以看出, 随着系统阶数的增加, 输出信号中的高频成分能量增强, 低频成分能量相应减弱. 根据绝热近似原理, 低频信号较高频信号更容易达到SR现象产生的条件.
图2(b)给出了系统输出幅值增益在不同耦合系数的条件下, 随着三态噪声稳态转移概率变化而发生的变化情况. 能够发现, 当
ε=0 时, 系统退化为非耦合系统, 系统输出幅值增益没能产生GSR现象. 随着耦合系数的增大, 系统输出幅值增益开始出现广义共振现象. 具体地说, 当耦合系数由0.5增大到2时, 系统的共振现象随之增强. 当耦合系数进一步由2逐渐增大到100时, 系统的共振现象逐渐减弱, 直至收敛. 这是因为粒子间适当的耦合作用能够为系统提供产生共振现象所需的非线性条件. 当耦合作用过大时, 两个粒子将被强大的耦合力连接成一个刚性的整体, 从而使系统的非线性减弱, 进而减弱系统的共振现象.从图2还可以看出, 系统输出幅值增益的共振峰均出现在q小于0.5时, 并在q等于0.5时取得最小值. 由于三态噪声将在q等于0.5时退化为双态噪声, 上述结果说明了在这一组参数条件下, 三态噪声与双态噪声相比, 能为系统输出创造更高的增益.
4.2 系统输出幅值增益随耦合系数
ε 变化的GSR从图4可以看出, 噪声稳态转移概率和系统阶数均对系统输出幅值增益随耦合系数变化的GSR行为有着显著的影响, 即随着q的增大, 系统将由双峰共振演变为单峰共振. 并且, 当噪声稳态转移概率较大时, 随着系统阶数的增大, 系统也将由双峰共振演变为单峰共振. 不失一般性, 在分别选择
q=0.3 和α=0.1 的条件下, 进一步分析了系统阶数和噪声稳态转移概率对系统GSR现象的影响, 结果如图5 所示.图5(a)给出了系统输出幅值增益在不同系统阶数的条件下, 随耦合系数的变化情况. 不难发现, 当系统阶数取值为0.2时, 系统输出幅值增益随耦合系数的变化产生双峰共振. 随着系统阶数的增大, 系统输出幅值增益逐渐由双峰共振演变为单峰共振, 且共振峰逐渐降低. 特别地, 当系统阶数取值为1, 即系统退化为整数阶系统时, 其输出幅值增益的峰值最小, 该现象说明分数阶系统能比传统整数阶系统得到更大的输出幅值增益.
图5(b)给出了系统输出幅值增益在不同噪声稳态转移概率的条件下, 随耦合系数变化而变化的情况. 可以看出, 当噪声稳态转移概率取值为0时, 系统输出幅值增益没能产生广义共振现象. 其原因是此时三态噪声退化为常数, 随机系统也相应地退化为确定性系统. 粒子
x1 和x2 的运动行为完全一致, 致使系统中的耦合项ε(x1−x2) 和ε(x2−x1) 消失, 进而消除了系统产生共振现象所必须的非线性条件. 当噪声稳态转移概率取值为0.25时, 系统输出幅值增益出现双峰共振. 当噪声稳态转移概率进一步增大到0.5, 即三态噪声退化为双态噪声时, 系统输出幅值增益出现单峰共振, 且共振峰的高度明显提高. 该结果表明, 在这一组参数条件下, 双态噪声能使系统输出幅值获得较三态噪声激励时更大的增益.5. 数值仿真验证
为验证理论结果的正确性, 下面通过数值仿真模拟模型(2)式所刻画的振子运动. 对于充分小的仿真时间间隔
Δt , 该模型在离散时间下的一阶近似表达式为{xi(tn)=yi(tn−1)Δtα−memo(xi,α),yi(tn)=[−γyi(tn−1)−mω2xi(tn−1)+ε(2∑k=1xk(tn−1)−2xi(tn−1))+A0cos(Ωtn−1)]Δt2−αm+ξi(tn)−memo(yi,2−α), (42) 其中,
{memo(xi,α)=n∑k=1(−1)k(αk)xi(tn−k),(αk)=α(α−1)(α−2)⋯(α−k+1)k!. (43) 在给定系统参数条件下, 取仿真时长
t=300 , 采样间隔Δt=10−2 , 可得系统输出信号时域图和频域图如图6所示.图6(a)为系统输出信号时域图, 可以看出, 受噪声的影响, 系统输出具有较大的随机性. 此时, 从时域图上已经无法判断出是否有正弦响应信号存在. 因此, 通过对时域输出信号进行傅里叶变换, 得到了与之对应的频域输出信号, 如图6(b)所示. 从图6(b)可清楚地看到, 系统输出在外部驱动频率(
Ω=0.4π )处出现了明显的尖峰, 表明系统输出响应中含有与外部驱动信号(A0cos(Ωt ))同频的正弦信号. 其峰值(2.149)代表系统对外部驱动信号的响应幅值, 噪声的存在使得该响应幅值与其真值之间存在随机偏差. 为此, 下面采用Monte-Carlo方法来消除该随机偏差.在图6的参数条件下, 重复进行N次仿真实验, 取N次结果的平均值作为系统稳态响应振幅的仿真值, 以(40)式确定的理论结果为参考值, 可得仿真误差随仿真次数的变化, 如图7所示.
图7(a)为仿真与理论结果对比图. 可以看出, 仿真结果始终在理论结果的附近随机波动, 表明仿真结果与理论结果之前确实存在随机偏差. 随着仿真次数N的增加, 随机波动的幅度逐渐减小. 为更加清楚地反映该随机偏差的收敛情况, 进一步给出了相应的绝对随机偏差随仿真次数N的变化图(图7(b)). 图7(b)表明, 随着仿真次数N的增加, 噪声带来的随机偏差将逐渐收敛. 特别地, 若以
δ0=10−3 为可接受的偏差门限, 则当N⩾2881 时, 即可满足该门限需求. 一般地, 噪声强度越大或者偏差门限越小, 则所需的仿真次数N越大.6. 结 论
本文研究了三态噪声激励下的分数阶耦合系统的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 L453
Google Scholar
[2] Wiesenfeld K, Moss F 1995 Nature 373 33
Google Scholar
[3] Gammaitoni L, Hänggi P, Jung P, Marchesoni F 1998 Rev. Mod. Phys. 70 223
Google Scholar
[4] Hänggi P 2002 Chem. Phys. Chem. 3 285
Google Scholar
[5] Zhang W, Xiang B R 2006 Talanta 70 267
Google Scholar
[6] Chen H, Varshney P K, Kay S M, Michels J H 2007 IEEE Trans. Signal Process. 55 3172
Google Scholar
[7] Soika E, Mankin R, Ainsaar A 2010 Phys. Rev. E 81 011141
Google Scholar
[8] Gang H, Ditzinger T, Ning C Z, Haken H 1993 Phys. Rev. Lett. 71 807
Google Scholar
[9] Pikovsky A, Zaikin A, Ma D L C 2002 Phys. Rev. Lett. 88 050601
Google Scholar
[10] Tessone C J, Mirasso C R, Toral R, Gunton J D 2006 Phys. Rev. Lett. 97 194101
Google Scholar
[11] 冷永刚, 王太勇, 郭焱, 吴振勇 2007 物理学报 56 30
Google Scholar
Leng Y G, Wang T Y, Guo Y, Wu Z Y 2007 Acta Phys. Sin. 56 30
Google Scholar
[12] 林敏, 黄咏梅, 方利民 2008 物理学报 57 2041
Google Scholar
Lin M, Huang Y M, Fang L M 2008 Acta Phys. Sin. 57 2041
Google Scholar
[13] 张莉, 刘立, 曹力 2010 物理学报 59 1494
Google Scholar
Zhang L, Liu L, Cao L 2010 Acta Phys. Sin. 59 1494
Google Scholar
[14] Berdichevsky V, Gitterman M 1999 Phys. Rev. E 60 1494
Google Scholar
[15] Yu T, Zhang L, Luo M K 2013 Phys. Scr. 88 045008
Google Scholar
[16] Zhong S C, Ma H, Peng H, Zhang L 2015 Nonlinear Dyn. 82 535
Google 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 171
Google Scholar
[19] Zhang L, Lai L, Peng H 2018 Phys. Rev. E 97 012147
Google Scholar
[20] Yu T, Zhang L, Ji Y D, Lai L 2019 Commun. Nonlinear Sci. Numer. Simul. 72 26
Google Scholar
[21] Lin L F, Wang H Q 2019 Nonlinear Dyn. 98 801
Google Scholar
[22] 姜磊, 赖莉, 蔚涛, 罗懋康 2021 物理学报 70 130501
Google Scholar
Jiang L, Lai L, Yu T, Luo M K 2021 Acta Phys. Sin. 70 130501
Google Scholar
[23] Ausloos M, Lambiotte R 2006 Phys. Rev. E 73 011105
Google Scholar
[24] Gitterman M, Klyatskin V I 2010 Phys. Rev. E 81 051139
Google Scholar
[25] Gitterman M, Shapiro I 2011 J. Stat. Phys. 144 139
Google Scholar
[26] 蔚涛, 张路, 罗懋康 2013 物理学报 62 120504
Google Scholar
Yu T, Zhang L, Luo M K 2013 Acta Phys. Sin. 62 120504
Google Scholar
[27] 蔚涛, 罗懋康, 华云 2013 物理学报 62 210503
Google Scholar
Yu T, Luo M K, Hua Y 2013 Acta Phys. Sin. 62 210503
Google Scholar
[28] Broeck C V D 1983 J. Stat. Phys. 31 467
Google Scholar
[29] Zhong S C, Wei K, Gao S L, Ma H 2015 J. Stat. Phys. 159 195
Google Scholar
[30] Mankin R, Ainsaar A, Reiter E 1999 Phys. Rev. E 60 1374
Google Scholar
[31] Lang R L, Yang L, Qin H L, Di G H 2012 Nonlinear Dyn. 69 1423
Google Scholar
[32] Peng H, Ren R B, Li P F, Yu T 2020 Phys. Scr. 95 075214
Google Scholar
[33] Xu P F, Jin Y F 2020 Chaos, Solitons Fractals 138 109857
Google Scholar
[34] Jin Y F, Wang H Q 2020 Chaos, Solitons Fractals 133 109633
Google Scholar
[35] Jin Y F, Niu S Y 2016 Int. J. Mod. Phys. B 30 1650102
Google Scholar
[36] Bagley R L, Torvik P J 1984 Appl. Mech. 51 294
Google Scholar
[37] Lutz L 2001 Phys. Rev. E 64 051106
Google Scholar
[38] He G T, Tian Y, Wang Y 2013 J. Stat. Mech. 9 26
Google Scholar
[39] Jiang L, Lai L, Yu T, Luo M K 2021 J. Stat. Mech. 2021 063208
Google Scholar
[40] Cubero D 2008 Phys. Rev. E 77 021112
Google Scholar
[41] Tang Y, Zou W, Lu J, Kurths J 2012 Phys. Rev. E 85 046207
Google Scholar
[42] Atsumi Y, Hata S, Nakao H 2013 Phys. Rev. E 2013 88
Google Scholar
[43] Oyarzabal R S, Szezech Jr J D, Batista A M 2017 Commun. Nonlinear Sci. Numer. Simul. 54 62
Google Scholar
[44] Zhang L, Xu L, Yu T, Lai L, Zhong S C 2021 Commun. Nonlinear Sci. Numer. Simul. 93 105499
Google Scholar
[45] Shapiro V E, Loginov V M 1978 Phys. A 91 563
Google 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] Benzi R, Sutera A, Vulpiani A 1981 J. Phys. A 14 L453
Google Scholar
[2] Wiesenfeld K, Moss F 1995 Nature 373 33
Google Scholar
[3] Gammaitoni L, Hänggi P, Jung P, Marchesoni F 1998 Rev. Mod. Phys. 70 223
Google Scholar
[4] Hänggi P 2002 Chem. Phys. Chem. 3 285
Google Scholar
[5] Zhang W, Xiang B R 2006 Talanta 70 267
Google Scholar
[6] Chen H, Varshney P K, Kay S M, Michels J H 2007 IEEE Trans. Signal Process. 55 3172
Google Scholar
[7] Soika E, Mankin R, Ainsaar A 2010 Phys. Rev. E 81 011141
Google Scholar
[8] Gang H, Ditzinger T, Ning C Z, Haken H 1993 Phys. Rev. Lett. 71 807
Google Scholar
[9] Pikovsky A, Zaikin A, Ma D L C 2002 Phys. Rev. Lett. 88 050601
Google Scholar
[10] Tessone C J, Mirasso C R, Toral R, Gunton J D 2006 Phys. Rev. Lett. 97 194101
Google Scholar
[11] 冷永刚, 王太勇, 郭焱, 吴振勇 2007 物理学报 56 30
Google Scholar
Leng Y G, Wang T Y, Guo Y, Wu Z Y 2007 Acta Phys. Sin. 56 30
Google Scholar
[12] 林敏, 黄咏梅, 方利民 2008 物理学报 57 2041
Google Scholar
Lin M, Huang Y M, Fang L M 2008 Acta Phys. Sin. 57 2041
Google Scholar
[13] 张莉, 刘立, 曹力 2010 物理学报 59 1494
Google Scholar
Zhang L, Liu L, Cao L 2010 Acta Phys. Sin. 59 1494
Google Scholar
[14] Berdichevsky V, Gitterman M 1999 Phys. Rev. E 60 1494
Google Scholar
[15] Yu T, Zhang L, Luo M K 2013 Phys. Scr. 88 045008
Google Scholar
[16] Zhong S C, Ma H, Peng H, Zhang L 2015 Nonlinear Dyn. 82 535
Google 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 171
Google Scholar
[19] Zhang L, Lai L, Peng H 2018 Phys. Rev. E 97 012147
Google Scholar
[20] Yu T, Zhang L, Ji Y D, Lai L 2019 Commun. Nonlinear Sci. Numer. Simul. 72 26
Google Scholar
[21] Lin L F, Wang H Q 2019 Nonlinear Dyn. 98 801
Google Scholar
[22] 姜磊, 赖莉, 蔚涛, 罗懋康 2021 物理学报 70 130501
Google Scholar
Jiang L, Lai L, Yu T, Luo M K 2021 Acta Phys. Sin. 70 130501
Google Scholar
[23] Ausloos M, Lambiotte R 2006 Phys. Rev. E 73 011105
Google Scholar
[24] Gitterman M, Klyatskin V I 2010 Phys. Rev. E 81 051139
Google Scholar
[25] Gitterman M, Shapiro I 2011 J. Stat. Phys. 144 139
Google Scholar
[26] 蔚涛, 张路, 罗懋康 2013 物理学报 62 120504
Google Scholar
Yu T, Zhang L, Luo M K 2013 Acta Phys. Sin. 62 120504
Google Scholar
[27] 蔚涛, 罗懋康, 华云 2013 物理学报 62 210503
Google Scholar
Yu T, Luo M K, Hua Y 2013 Acta Phys. Sin. 62 210503
Google Scholar
[28] Broeck C V D 1983 J. Stat. Phys. 31 467
Google Scholar
[29] Zhong S C, Wei K, Gao S L, Ma H 2015 J. Stat. Phys. 159 195
Google Scholar
[30] Mankin R, Ainsaar A, Reiter E 1999 Phys. Rev. E 60 1374
Google Scholar
[31] Lang R L, Yang L, Qin H L, Di G H 2012 Nonlinear Dyn. 69 1423
Google Scholar
[32] Peng H, Ren R B, Li P F, Yu T 2020 Phys. Scr. 95 075214
Google Scholar
[33] Xu P F, Jin Y F 2020 Chaos, Solitons Fractals 138 109857
Google Scholar
[34] Jin Y F, Wang H Q 2020 Chaos, Solitons Fractals 133 109633
Google Scholar
[35] Jin Y F, Niu S Y 2016 Int. J. Mod. Phys. B 30 1650102
Google Scholar
[36] Bagley R L, Torvik P J 1984 Appl. Mech. 51 294
Google Scholar
[37] Lutz L 2001 Phys. Rev. E 64 051106
Google Scholar
[38] He G T, Tian Y, Wang Y 2013 J. Stat. Mech. 9 26
Google Scholar
[39] Jiang L, Lai L, Yu T, Luo M K 2021 J. Stat. Mech. 2021 063208
Google Scholar
[40] Cubero D 2008 Phys. Rev. E 77 021112
Google Scholar
[41] Tang Y, Zou W, Lu J, Kurths J 2012 Phys. Rev. E 85 046207
Google Scholar
[42] Atsumi Y, Hata S, Nakao H 2013 Phys. Rev. E 2013 88
Google Scholar
[43] Oyarzabal R S, Szezech Jr J D, Batista A M 2017 Commun. Nonlinear Sci. Numer. Simul. 54 62
Google Scholar
[44] Zhang L, Xu L, Yu T, Lai L, Zhong S C 2021 Commun. Nonlinear Sci. Numer. Simul. 93 105499
Google Scholar
[45] Shapiro V E, Loginov V M 1978 Phys. A 91 563
Google 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)
计量
- 文章访问数: 4887
- PDF下载量: 91
- 被引次数: 6