路徑積分法在隨機動力系統中的應用

莫協強

（梧州學院 數理系，廣西 梧州 543000）

路徑積分法在隨機動力系統中的應用

莫協強

（梧州學院 數理系，廣西 梧州 543000）

路徑積分法在隨機動力系統中有著重要作用，本文介紹了基于Guass-Lgenedre公式的路徑積分法的原理和算法，并研究了諧和激勵與隨機激勵作用下的Duffing一Ryaelihg振子，推導用于計算時間上的平均概率密度路徑積分數值表達式，探討隨機跳躍現象存在性和概率密度峰的個數之間的關系.

路徑積分法；隨機動力系統；Guass-Legenedre公式

它的基本思想是在空間和時間路徑積分的離散化,將取代路徑和積分、非線性隨機動力系統可能有復雜結構的確定性和隨機響應反應,Yim和Naess,等等…隨機響應分析混沌響應的概率密度,即通過不確定性系統中激勵引進隨機干擾,混沌吸引子的不確定性系統存在也能有效地使用相空間的隨機響應的概率密度演化來描述.本文主要利用路徑積分方法研究了非線性隨機動力系統混亂的反應.

1 路徑積分法的原理

在空間和時間是途徑和離散化,而不是整體,這是路徑積分方法的基本思想.即通過短期的轉移概率密度的連接形式的全球轉移概率密度,得到狀態向量的聯合概率密度函數.雖然每條路徑的貢獻的定性概念的總和是透明的,但它應該給予的總和更精確的數學定義,路徑是一個高階無限數目的音樂,是道路的間隔什么樣的措施應給與不明顯.目的是為了給這樣的數學定義,該定義是挺麻煩的.路徑積分法是最優越的特性可以得到負面,更準確的尾部概率密度.此外,路徑積分法也可以計算非平穩信號的系統暫態概率密度和問題,如第一次通過.最簡單的定義路徑積分方法是連續的過程離散化在空間和時間限制的小格子點通透.然而,一個連續的過程的離散并不是唯一的.因此,對于不同的離散規則,你有許多不同的路徑積分方法.為建立路徑積分協變的,需要選擇特定離散規則,導致許多研究人員并沒有被具體的離散規則,提出了各種各樣的推導的路徑積分法協變.Wehner和wolfer路徑積分方法研究首次提出以尋求“非線性隨機“FPK方程可得到精確穩態解形式解.給出不確定的道路路徑積分的近似解的初始值適當的點,路徑積分形式的解決方案可能在平衡熱力學領域熱力學平衡延續的概念.實際上大部分人致力于用路徑積分表示下列FPK方程：

其中P（q1,q2,…,qn,t）是滿足上述“非線性隨機”FPK方程的概率密度，是擴散張量，是漂移向量，重復下標表示求和.即將方程（1.1）的解表示為如下形式的泛函積分：

其中Dμ(q)是積分測度，通常看作是Onsgear-Machl up泛函.若使等式(1.2)是協變形式，那就是當來自行星光狀態矢量空間的非線性隨機點變換,Onsgaer-Machl up功能也根據通常的計算規則了相應的改造,那么問題將變得更為復雜.因此,困難的路徑積分是不可能Onsgaer-Machl up唯一肯定的功能.如果你想獲得基于路徑的非平衡熱力學理論一致的積分,那么這個屬性是必不可少的.最簡單的定義路徑積分方法是連續的過程離散在時間和空間充滿了小有限網格點.然而,一個連續的離散化過程不僅這樣的話,分別對應著不同的離散規則,你有許多不同的路徑積分方法.建立了路徑積分協變的,需要選擇特定離散規則.在所有可能的路徑,一個特定的路徑由適當的條件決定,最低的原則是表示,此次行動是最桶色條件的一種方式.換句話說.存在一定數量的年代,我們可以計算出,每條路徑,經典軌道;是最少的年代的軌道.事實上,真正的條件S只是一個極值.也就是說,如果這路徑;X稍作改變,S價值為一階來表達相同的,所以我們不給出具體如此在各種規則的建議是協變的沿程積分方法.在實際應用中,即使路徑積分法的概念具有一定的參考價值,而是解析解和數值,公式(1.2)的計算是不太可能推廣到一般的情況,wehner和wol f er首先討論基于方程(1.2)的數值計算方法,提出了等效離散形式,起點是路徑,路徑積分或離散網格.的目的是尋求FPK方程可得到精確穩態解(1.1)的解決方案,路徑,并能很好的為基礎的數值計算方法.路徑和通常寫成：

其中 G(qi+1,qi,τ)=μiexp[-τv(qi+1,qi,τ)]稱作短時傳播子，不同離散規則有不同的G，傳播子的唯一條件是它滿足方程（1.1）達到O(τ2).確定數值方法的簡單的、最準確的路徑和公式是：

2 基于Gauss-Legendre公式的路徑積分法

在許多喬條件出現的特色功能是一種高斯分布,泊松分布假設隨機脈沖,從而介紹其二階近似的指數關系,并得到高斯功能,事實上,一些本質的物理過程看起來像個他們分布.在傳統的風險理論,高斯分布符合物理現象.這是因為很多隨機出現的綜合結果的獨立的事情.這就是那個機遇去中心極限定理理論的結論,同樣的結論適用于分布的功能.因為許多路徑積分可以劃分為若干路徑積分因素、制度是在一個特定國家和功能,可以寫每種模式下的是產品的一個因素.最基本的理念是路徑積分在空間和時間離散化、道路、而不是整體,即通過連接短期轉移概率密度形成全球轉移概率密度,得到狀態向量的聯合概率密度函數.設(X)t是n維狀態向量，其演化概率密度為

其中q(x,t|x(0),t0)是轉移概率密度，由非線性隨機動力學系統相應的FPK方程確定，p(x(0),t0)是X（t）在t=t 0處初始概率分布，R 是X（t）的 n維狀態空間.將[t 0，t]分成 N 個子區間可得：

實際上，路徑積分是在縮減的狀態空間Rx內積分，Rx之外的區域轉移概率密度因為充分小而忽略不計.即：

利用不同的插值公式就形成不同的路徑積分數值方法.目前主要有兩類廣泛應用的路徑積分數值方法.一種是，Naess等利用三次B樣條插值來離散積分式(2.3)，可以計算尾部低于10-6數量級的平穩概率密度.進一步,Naess提出了概率密度的B樣條插值的對數函數代替直接在概率密度自身的插值,能明顯地提高計算精度和效率.另一個是,玉將(2.3)根據Guass-Lgenedre公式離散化,得到Guass-Legnedre基于路徑的積分公式數值計算方法.高斯積分點其他的點概率密度能通過插值.所以,只有我i時刻高斯積分點計算相應的轉移概率密度,大大減少了計算量.最終可得到適合編程實現的表達式.一維情形即：

其中K是子區間數,Lk是第k子區間的高斯積分點數，δk是第k子區間的長度，xkl是高斯積分點，ckl是相應的權值.這種方法的優點是沒有重新構造插值函數數值積分團聚的,可以直接利用Guass-Lgenedre離散積分公式的類型(2.4).對時變系統,特別是膨脹系數漂移或包含一個傳播周期函數的系統,它的時刻就在時間迫使一個騎表達期刊,用簡單的方法是通過求解一個周期的時間步ChangShang力矩方程出發.每一刻的解方程周期保持不變,為下星期期間短期轉移概率密度的均值和方差的估計.路徑積分的數值方法,主要計算是一種轉移概率密度的近似,因為相同的電力系統的時候,正好方程表達不依賴時間的,只要計算一下方程;不為高斯白噪聲的動態激勵機制,一般沒有辦法直接轉矩方程d只有通過開拓系統d系統推力矩方程出發,再用路徑積分擴張維數值對應的FPK方程可得到精確穩態解高維系統的分析.噪音強度在一定的非線性隨機系統狀態估計的概率密度演化可以用來描繪了混沌算子的吸引結構特點.

3 路徑積分法基于Gauss-Legenedre公式求解時間上平均概率密度

對時變系統,特別是漂移或包含在擴散系數的周期函數的系統,介紹了概率密度定義的平均時間.

其中T=2π/ω,ω一般是諧和激勵或周期系數的角頻率.當（3.1）式中t'充分大，在整個諧和激勵周期T上平均的概率密度可逼近系統響應的近似平穩概率密度.將Yu的路徑積分數值方法推廣到時間上平均的概率密度的計算，不難導出利用Gauss-Legnedre公式計算平均概率密度的表達式.即積分表達式(3.1)可離散成如下的Guass-Legendre積分表達式

其中K是[t'，t'+T]內的子區間數，Lk是第k子區間的高斯積分點，δK是第k子區間的長度，每一tkl是一個周期內的高斯積分時刻，而ckl是相應的權重，而p(x1,x2,tkl)的計算可直接參考上面的計算方法.Guass-Legnedre于使用路徑積分公式法研究了高斯白噪聲的沖擊只有幾個隨機非線性動態系統,結果表明該方法在計算盡頭的一個更小的概率密度是歷史上其他數值方法精度高,本研究是首次通過更重要.合理假設短期轉移概率密度為高斯分布,利用高斯三種方法已關閉的力矩方程出發,推導出能算出這道短的轉移概率密度的一階矩和二階矩.在此基礎上,利用路徑積分法求解系統的平滑概率密度.不同于簡單形狀的激情,基于路徑積分方法Guass-Lgenedre公式的隨機性和沖突,促進共同作用的非線性隨機動力學系統,利用非線性動態系統參與的隨機研究更為復雜.這些研究結果將和其他數值方法的結果進行了比較,得出了計算尾概率密度具有較高的計算精度.

4 結語

在路徑積分對幾串理論的這些應用中，如果被積函數是高斯型的.那么我們就可以大量地應用這個技術.這些問題正好都是一些用其他方法也可以解決的問題，不一定要用路徑積分.人們可以有理由懷疑路徑積分有無實際用途.如果問題不是高所型的，則至少可以用路徑積分研究它并將其公式化.然而，路徑積分方法使一個問題的各種表述方式之間可以迅速地相互轉換，并且常常清楚而迅速地得到關于某個關系式的提示，而用更通常的方式推導它會要慢一些.

〔1〕李杰.隨機結構分析的擴階系統方法(Ⅱ)——結構動力分析[J].地震工程與工程振動,1995(04).

〔2〕李杰.隨機結構動力分析的擴階系統方法[J].工程力學,1996(01).

〔3〕廖松濤,李杰.隨機結構正交展開分析的Ritz動力聚縮法[J].計算力學學報,2002(01).

〔4〕李杰.復合隨機振動分析的擴階系統方法[J].力學學報,1996(01).

〔5〕馬少娟,徐偉,李偉,靳艷飛.基于 Chebyshev多項式逼近的隨機vanderPol系統的倍周期分岔分析[J].物理學報,2005(08).

〔6〕孫曉娟,徐偉,馬少娟.含有界隨機參數的雙勢阱DuffingvanderPol系統的倍周期分岔[J].物理學報,2006(02).

〔7〕馬少娟,徐偉,李偉.基于Laguerre多項式逼近法的隨機雙勢阱Duffing系統的分岔和混沌研究 [J].物理學報,2006(08).

〔8〕徐偉,方同,戎海武.有界窄帶激勵下具有黏彈項的 Duffing 振子[J].力學學報,2002（05）.

〔9〕戎海武,徐偉,孟光,方同.諧和與隨機噪聲聯合作用下非線性非線性隨機系統的響應[J].應用力學學報,2001（04）.

〔10〕戎海武,徐偉,王向東,孟光,方同.窄帶噪聲作用下二自由度非線性非線性隨機系統的響應[J].力學學報,2001（06）.

〔11〕朱位秋,黃志龍,應祖光.非線性隨機動力學與控制的哈密頓理論框架[J].力學與實踐,2002（03）.

〔12〕沈仙華,李璋靜,趙玉林.非線性隨機的單自由度系統相加型聯合共振[J].科學技術與工程,2009（18）.

O172.2

A

1673-260X（2012）05-0009-03