随机时滞Duffing系统的13亚谐共振响应.docVIP

随机时滞Duffing系统的1/3亚谐共振响应 孙中奎* 徐 伟 (西北工业大学应用数学系，西安 710072) 摘要 本文应用随机多尺度方法研究了窄带激励下带有时滞反馈的Duffing系统的1/3亚谐共振响应。在分离系统的快慢变量的基础上讨论了系统的稳定性和分叉现象，系统分析了时滞和窄带噪声对Duffing系统稳态响应、稳定性及分叉的影响。并运用数值方法验证了文中所得结果的有效性。 关键词：时滞Duffing系统；窄带噪声；1/3亚谐共振响应；分叉 中图分类号：TB535；O322 引言 时滞是普遍存在的。事实上，由于反应的滞后性、传输速度的有限性等原因，时滞现象在现实世界中几乎是不可避免的，因此时滞广泛地存在于各种系统之中，例如，人口系统、神经网络系统、生物系统、经济系统、复杂网络系统等等。所以，时滞系统是对自然规律更本质、更真实地刻画，一般动力系统只是对现实世界规律的一个近似。 时滞系统的数学模型不再是一般的微分方程(组)，而是时滞微分方程(组)[1-4]，时滞系统的解空间是无穷维[1-4]的，这给时滞系统的解析求解带来了巨大的困难。目前国内外的广大研究人员在时滞系统解的存在性、唯一性以及稳定性等方面作了大量的工作：在文献[5]中，Massera证明了常微分方程中周期解的存在性；文献[6]推广了文献[5]的结果，使得Massera定理能够适用于时滞系统；Yoshizawa在文献[7]中证明了在时滞系统中当时滞小于或等于方程的周期时，如果系统的解一致有界、一致最终有界，则该时滞系统一定存在-周期解；Kaplan和Yorke在文献[8]中对一些具体类型的时滞系统，通过确立与时滞系统相应的常微分方程来判断此时滞系统周期解的存在性，该方法目前被称为Kaplan-Yorke法；文献[9,10]运用不动点定理证明了一类特殊形式的时滞系统周期解的存在性；Lyapunov函数法也被用来研究时滞系统周期解的存在性[11-13]。此外，常数变异法[1-4]、分步法[1, 3]、Laplace变换法[1, 3]等也被广泛地用在对时滞系统的求解过程中。特别地，李继彬等[14]将时滞系统和Hamilton理论结合起来以期得到时滞系统的周期解，给研究者提供了一个很好的思路和方法。然而，由于时滞系统的复杂性，目前只有极少数特殊的时滞系统可以求出精确解，因此，各种摄动方法得到了广泛地运用：Casal等将Poincare-Lindstedt方法用于对时滞系统周期解和分岔的研究[15]；McDonald等将谐波平衡法运用到对时滞系统的研究中[16]；多尺度方法也被广泛用于对时滞系统的分析，DAS等运用多尺度方法研究了一类自治和非自治的时滞系统在分岔点附近的周期解[17]，胡海岩在[18]中运用多尺度方法研究了含有时滞反馈的Duffing系统的主共振与亚谐共振并对研究时滞系统的其他方法进行了总结；此外，中心流型简化理论也被应用到对时滞系统地分析中，Xu等人在[19,20]中分别研究了带有时滞速度反馈的Duffing系统和带有时滞位移反馈的van der Pol-Duffing系统的动力学行为，借助中心流型和平均理论，得到了约化系统，讨论了时滞导致的分岔、拟周期解等现象，得到很好的结论。 本文在文献[21]已有工作的基础上，运用随机多尺度方法深入讨论带有时滞反馈的Duffing系统在谐和激励与窄带激励下的1/3亚谐共振响应，系统地分析了噪声、时滞对Duffing系统稳态响应、稳定性和分叉的影响。 1 基本公式推导 考虑如下窄带激励下的时滞Duffing系统 (1) 圆点表示对时间的导数， 是小参数， 和分别表示刚度项系数和自然频率，是非线性项的强度，是时滞项的系数，是随机过程，采用Wedig[22]形式 (2) 是随机激励强度，是中心频率，是标准的Wiener过程，是带宽。 在系统(1)中令，则上述系统转化为 . (3) 为了研究系统的1/3亚谐共振响应，我们令 则 . 将系统(1)的解展开成如下形式 (4) 其中分别是快慢时间尺度。记，则 (5) 将方程(4)和(5)代入方程(3)并且令的相同阶的系数相等，得到 , (6) , (7) 求解方程(6)有 ， (8) 这其中。 将(8)式代入方程(7)并消除久期项 (9) 将表示成极坐标形式 ， (10) 将式 (10)代入方程(9)中，并分离实部和虚部，得到 (11) 其中，。 因此，系统(1)的一阶近似解为