王國(guó)威，程慶華，徐大海 (長(zhǎng)江大學(xué)物理科學(xué)與技術(shù)學(xué)院，湖北 荊州 434023

王淑平 (鎮(zhèn)賚第四中學(xué)，吉林 鎮(zhèn)賚 137300)

周期信號(hào)和關(guān)聯(lián)噪聲作用下非對(duì)稱雙穩(wěn)系統(tǒng)的統(tǒng)計(jì)特性研究

王國(guó)威，程慶華，徐大海 (長(zhǎng)江大學(xué)物理科學(xué)與技術(shù)學(xué)院，湖北 荊州 434023

王淑平 (鎮(zhèn)賚第四中學(xué)，吉林 鎮(zhèn)賚 137300)

研究了關(guān)聯(lián)噪聲和周期信號(hào)共同驅(qū)動(dòng)下非對(duì)稱雙穩(wěn)系統(tǒng)的穩(wěn)態(tài)性質(zhì)和各階矩。根據(jù)劉維方程、諾維科夫原理等方法，得到了系統(tǒng)相應(yīng)的?？?普朗克方程，并計(jì)算了相應(yīng)的穩(wěn)態(tài)概率分布函數(shù)Pst(x)的表達(dá)式。定義R=D/Q(其中,D是乘性噪聲強(qiáng)度，Q是加性噪聲強(qiáng)度)，分3種情況(R=1,R>1,R<1)討論了周期信號(hào)振幅A和頻率Ω對(duì)系統(tǒng)定態(tài)性質(zhì)和各階矩的影響。結(jié)果表明：①當(dāng)R=1時(shí)，Pst(x)隨態(tài)變量x變化的三維圖呈現(xiàn)1條等高的雙峰曲線，且態(tài)變量的各階矩均處于3種情況所對(duì)應(yīng)的3條曲線的最上方；②當(dāng)R>1時(shí)，Pst(x)隨態(tài)變量x變化的三維圖呈現(xiàn)1條左低右高的不對(duì)稱雙峰曲線，且態(tài)變量的各階矩均處于3種情況所對(duì)應(yīng)的3條曲線的最下方；③當(dāng)R<1時(shí)，Pst(x)隨態(tài)變量x變化的三維圖呈現(xiàn)1條左低右高的不對(duì)稱雙峰曲線，且態(tài)變量的各階矩均處于3種情況所對(duì)應(yīng)的3條曲線的中間位置上。

非對(duì)稱雙穩(wěn)系統(tǒng)；穩(wěn)態(tài)；矩；關(guān)聯(lián)噪聲；郎之萬(wàn)方程

近年來(lái)噪聲對(duì)雙穩(wěn)系統(tǒng)的穩(wěn)態(tài)和瞬態(tài)性質(zhì)的影響一直是人們關(guān)注的一個(gè)熱點(diǎn)。研究發(fā)現(xiàn)[1]，在某種特定情況下，外部噪聲(乘性噪聲)和內(nèi)部噪聲(加性噪聲)可能來(lái)自同一個(gè)噪聲源，并且他們之間彼此關(guān)聯(lián)。關(guān)聯(lián)噪聲對(duì)非線性系統(tǒng)的動(dòng)力學(xué)研究帶來(lái)了很大的影響，并且被應(yīng)用于很多領(lǐng)域,如生物系統(tǒng)、激光系統(tǒng)、雙穩(wěn)系統(tǒng)等。在關(guān)于關(guān)聯(lián)噪聲和周期信號(hào)共同驅(qū)動(dòng)下的雙穩(wěn)系統(tǒng)的穩(wěn)態(tài)性質(zhì)方面的文獻(xiàn)中，研究對(duì)象一般是對(duì)稱的雙穩(wěn)系統(tǒng)[2-6]，而對(duì)非對(duì)稱的雙穩(wěn)系統(tǒng)研究較少，且沒(méi)有考慮周期信號(hào)的作用。為此，筆者對(duì)周期信號(hào)對(duì)關(guān)聯(lián)噪聲作用下非對(duì)稱雙穩(wěn)系統(tǒng)的定態(tài)性質(zhì)和矩的影響進(jìn)行了研究。

1 穩(wěn)態(tài)概率密度函數(shù)

關(guān)聯(lián)噪聲和周期信號(hào)共同驅(qū)動(dòng)的非對(duì)稱雙穩(wěn)系統(tǒng)[7]，可由以下一維無(wú)量綱郎之萬(wàn)方程給出[8]:

(1)

式中，x為態(tài)變量；t為時(shí)間；r為系統(tǒng)非對(duì)稱參數(shù)；A為周期信號(hào)振幅；Ω為信號(hào)頻率；ξ(t)和η(t)是具有零均值和非零自相關(guān)時(shí)間的白噪聲，其統(tǒng)計(jì)性質(zhì)如下[9]:

〈ξ(t)〉=〈η(t)〉=0 〈ξ(t)ξ(t′)〉=2Dδ(t-t′)

(2)

式中，D和Q分別為乘性噪聲強(qiáng)度和加性噪聲強(qiáng)度；λ為ξ(t)和η(t)之間的關(guān)聯(lián)系數(shù)。

根據(jù)文獻(xiàn)[10]，作如下代換：

f(x)=-r+x-x3+AcosΩtg1(x)=xg2(x)=1

(3)

則式(1)可變形為:

(4)

由劉維方程[11]得:

(5)

對(duì)式(5)作系綜平均，且由P(x,t)=〈δ(x(t)-x)〉，可得:

(6)

根據(jù)諾維科夫原理[12]:

(7)

對(duì)式(6)進(jìn)行整理，得:

(8)

對(duì)式(8)進(jìn)行整理，得:

(9)

(10)

(11)

在定態(tài)條件下，式(9)的穩(wěn)態(tài)解為:

(12)

對(duì)式(12)積分得穩(wěn)態(tài)概率分布函數(shù):

(13)

態(tài)變量x的各階矩和標(biāo)準(zhǔn)化方差分別為：

(14)

其中，態(tài)變量x的均值表示式(14)中n=1時(shí)的矩。

2 乘性噪聲強(qiáng)度D和加性噪聲強(qiáng)度Q的比率對(duì)定態(tài)性質(zhì)和各階矩的影響

定義R=D/Q，即乘性噪聲強(qiáng)度D和加性噪聲強(qiáng)度Q的比率。根據(jù)周期信號(hào)和關(guān)聯(lián)噪聲共同驅(qū)動(dòng)下的非對(duì)稱雙穩(wěn)系統(tǒng)的穩(wěn)態(tài)概率分布函數(shù)、各階矩和標(biāo)準(zhǔn)化方差的表達(dá)式，討論3種情況下(R=1,R<1,R>1)乘性噪聲強(qiáng)度和加性噪聲強(qiáng)度的比率R對(duì)系統(tǒng)的定態(tài)性質(zhì)、各階矩和標(biāo)準(zhǔn)化方差的影響。

2.1R對(duì)穩(wěn)態(tài)概率密度分布的影響

r=7.4;λ=0.6;A=0.5;Ω=1,t=0.01

圖1所示為穩(wěn)態(tài)概率分布函數(shù)Pst(x)隨態(tài)變量x和周期信號(hào)振幅A在R=1、R>1和R<1時(shí)的變化圖。從圖1(a)可以看出，當(dāng)R=1時(shí)，Pst(x)在沿著態(tài)變量x變大的方向上，出現(xiàn)了2個(gè)峰值，且左右2個(gè)峰值高度基本一致，但左邊的峰比右邊的峰要寬。隨著周期信號(hào)振幅A的增加，2峰出現(xiàn)不同的變化，即左邊的峰值高度基本無(wú)變化，而右邊的峰值高度出現(xiàn)明顯的波動(dòng)。從圖1(b)和圖1(c)可以看出，在R<1和R>1的情況下，左邊的峰值高度均變小，2峰均不等高，這和R=1時(shí)明顯不同。但是隨著周期信號(hào)振幅A的增加，2峰峰值高度均會(huì)出現(xiàn)波動(dòng)。

r=-10;λ=0.6;Ω=10,t=0.01 r=-10,λ=0.6,A=100,t=0.01

r=10;λ=0.6;Ω=100,t=0.01 r=10,λ=0.6,A=100,t=0.01

2.2R對(duì)一階矩〈x〉的影響

圖2所示為一階矩〈x〉(態(tài)變量均值)作為周期信號(hào)振幅A的函數(shù)在R=1、R>1和R<1時(shí)的變化曲線圖。從圖2可以看出，一階矩〈x〉是周期信號(hào)振幅A的單調(diào)增函數(shù)，且在R=1時(shí)，對(duì)于輸入系統(tǒng)的某一確定的振幅，其均值大于R<1和R>1時(shí)系統(tǒng)的態(tài)變量均值。圖3給出了一階矩〈x〉作為周期信號(hào)頻率Ω的函數(shù)的變化曲線，從圖3可以看出，一階矩〈x〉隨著周期信號(hào)頻率Ω增加而單調(diào)減小，而且當(dāng)周期信號(hào)頻率Ω取值較小時(shí)，〈x〉隨Ω的增加而減小的趨勢(shì)較緩。當(dāng)Ω的取值大于20時(shí)，一階矩〈x〉隨著Ω增加而減小的較快。而且，當(dāng)R=1時(shí)，一階矩〈x〉曲線整體位于R<1和R>1時(shí)〈x〉曲線的上方。

2.3R對(duì)二階矩〈x2〉的影響

圖4所示為態(tài)變量二階矩〈x2〉隨周期信號(hào)振幅A變化曲線圖。從圖4可以看出，當(dāng)輸入系統(tǒng)信號(hào)振幅較小時(shí)，增加A的取值對(duì)二階矩〈x2〉的影響不明顯，但是當(dāng)A的取值大于25時(shí)，隨著A的增加，二階矩〈x2〉明顯變大，而且在R=1時(shí)〈x2〉隨著A的增加而變大的趨勢(shì)要比R>1和R<1時(shí)要快。圖5所示為二階矩作為周期信號(hào)頻率Ω的函數(shù)隨R變化曲線圖。從圖5中可以看出，〈x2〉是Ω的單調(diào)減函數(shù)，且R=1時(shí)，二階矩〈x2〉曲線整體位于R>1和R<1時(shí)〈x2〉曲線的上方。

r=10;λ=0.6;Ω=100,t=0.01 r=10,λ=0.6,A=100,t=0.01

r=10;λ=0.6;Ω=100,t=0.01 r=10,λ=0.6,A=100,t=0.01

2.4R對(duì)三階矩〈x3〉的影響

圖6所示為三階矩〈x3〉作為周期信號(hào)振幅A的函數(shù)隨R變化曲線圖。從圖6可以看出，〈x3〉是A的單調(diào)增函數(shù)。當(dāng)A小于20時(shí)，隨著A的增加，三階矩〈x3〉增加趨勢(shì)較緩；但是當(dāng)A大于20時(shí)，三階矩〈x3〉則隨著A的增加而迅速增加，而且當(dāng)R=1時(shí)〈x3〉隨著A的增加而增加的趨勢(shì)要比R>1和R<1時(shí)要快。圖7所示為三階矩〈x3〉作為輸入周期信號(hào)頻率Ω的函數(shù)隨R變化曲線圖。從圖7可以看出，三階矩〈x3〉是Ω的單調(diào)減函數(shù)，且當(dāng)Ω的取值小于20時(shí)，〈x3〉隨Ω的增加而緩慢減小；當(dāng)Ω的取值大于20時(shí)， 〈x3〉隨Ω的增加而迅速減?。欢耶?dāng)R=1時(shí)，其〈x3〉曲線整體位于R>1和R<1時(shí)〈x3〉曲線的上方。

2.5R對(duì)態(tài)變量的標(biāo)準(zhǔn)化方差κ2的影響

以輸入周期信號(hào)振幅A和頻率Ω為變量，態(tài)變量x的標(biāo)準(zhǔn)化方差κ2在R=1、R>1和R<1時(shí)的變化曲線圖如圖8和圖9所示。從圖8可以看出，在R>1和R<1時(shí)，隨著輸入周期信號(hào)振幅A增加，κ2迅速減小至趨于0附近；但是當(dāng)R=1時(shí)，隨著A增加，κ2先緩慢減小，然后迅速減小，最后緩慢減小至趨于0。從圖9可以看出，κ2隨著周期信號(hào)頻率Ω增加在R=1、R>1和R<1時(shí)的變化趨勢(shì)基本一致：當(dāng)Ω較小時(shí)，隨著Ω增加，κ2基本無(wú)變化；當(dāng)Ω的取值大于60時(shí)，κ2則隨著Ω的增加而變大。另外，當(dāng)R>1時(shí)，κ2曲線整體高于R<1和R=1時(shí)的κ2曲線。

3 結(jié) 語(yǔ)

通過(guò)理論計(jì)算，得到了關(guān)聯(lián)噪聲和周期信號(hào)驅(qū)動(dòng)下非對(duì)稱雙穩(wěn)系統(tǒng)的近似?？?普朗克方程，并在定態(tài)條件下推導(dǎo)了相應(yīng)的穩(wěn)態(tài)概率分布函數(shù)和各階矩的表達(dá)式。在定義乘性噪聲強(qiáng)度和加性噪聲強(qiáng)度比率 的情況下，對(duì)計(jì)算結(jié)果進(jìn)行了數(shù)值分析，討論了 對(duì)系統(tǒng)的統(tǒng)計(jì)性質(zhì)的影響,對(duì)繼續(xù)研究系統(tǒng)的平均首次通過(guò)時(shí)間、信噪比等具有一定理論指導(dǎo)意義。

[1]Jia Y, Li J R. Stochastic system with colored correlation between white noise and colored noise[J].Physica A,1998,252:417-427.

[2] 寧麗娟，徐偉，楊曉麗.色關(guān)聯(lián)噪聲驅(qū)動(dòng)的非對(duì)稱雙穩(wěn)系統(tǒng)中平均首次穿越時(shí)間的研究[J].物理學(xué)報(bào), 2007, 56 (1)：25-29.

[3] 周丙常，徐偉.關(guān)聯(lián)噪聲驅(qū)動(dòng)的非對(duì)稱雙穩(wěn)系統(tǒng)的隨機(jī)共振[J].物理學(xué)報(bào), 2008, 57 (4)：2035-2040.

[4] Wang C J, Dai Z C, Mei D C. Asymmetric effects on escape rates of bistable system [J].Commun Theor Phys,2011,56(6):1041-1044.

[5] 靳艷飛，徐偉，李偉，等. 色關(guān)聯(lián)的色噪聲驅(qū)動(dòng)的雙穩(wěn)杜芬模型的穩(wěn)態(tài)分析[J].動(dòng)力學(xué)與控制學(xué)報(bào)，2005,3(2):60-65.

[6] 王兵，吳秀清，邵繼紅.雙色噪聲驅(qū)動(dòng)非對(duì)稱雙穩(wěn)系統(tǒng)平均第一穿越時(shí)間研究[J].物理學(xué)報(bào)，2009, 58(3):1391-1395.

[7] 張靜靜，靳艷飛. 非高斯噪聲驅(qū)動(dòng)下非對(duì)稱雙穩(wěn)系統(tǒng)的平均首次穿越時(shí)間與隨機(jī)共振研究[J].物理學(xué)報(bào), 2011, 60 (12): 120501-5.

[8] 周丙常，徐偉，張曉燕. 2種不同類型噪聲驅(qū)動(dòng)下的非對(duì)稱雙穩(wěn)系統(tǒng)的隨機(jī)共振[J].應(yīng)用力學(xué)學(xué)報(bào), 2009, 26(1):27-30.

[9] 胡崗. 隨機(jī)力與非線性系統(tǒng)[M].上海：上海科技教育出版社，1994.

[10] Liang G Y, Cao L, Wu D J. Approximate Fokker-Planck equation of system driven by multiplicative colored noises with colored cross-correlation[J].Physica A,2004,335:371-384.

[11] Jia Y, Zheng X P, Hu X M, et al. Effects of colored noise on stochastic resonance in a bistable system subject to multiplicative and additive noise[J]. Physical Review E, 2001, 63: 031107.

[12] Luo X Q, Wu D, Zhu S Q. Stochastic resonance in a time-delayed bistable system with colored coupling between noise terms[J].International Journal of Modern Physics B, 2012，26(30):1250149.

O431.2

A

1673-1409(2013)28-0018-04

[編輯] 李啟棟