斯特拉托諾維奇積分

隨機過程中,斯特拉托諾維奇積分菲斯克-斯特拉托諾維奇積分是一種隨機積分,冠名於同時發展出該技術的魯斯蘭·斯特拉托諾維奇英語Ruslan Stratonovich唐納德·菲斯克。它是伊藤積分最常見的替代品。應用數學中通常選擇伊藤積分,而物理學中則頻繁使用斯特拉托諾維奇積分。

在某些場景下,斯特拉托諾維奇積分更易於操作。斯特拉托諾維奇積分的定義使其滿足鏈式法則,而伊藤積分則並非如此。

遇到斯特拉托諾維奇積分的一類主要場景是斯特拉托諾維奇隨機微分方程(SDE)的解。其與伊藤隨機微分方程是等價的,並且可以根據何種更為方便而在兩種定義之間進行切換。

定義

斯特拉托諾維奇積分可以用類似於黎曼積分的方式(黎曼和極限)來定義。設   是一個維納過程 適應於維納過程的自然濾過 (概率論)英語Filtration (probability_theory)的一個半鞅

類似黎曼-斯蒂爾傑斯積分中的策略,我們考慮區間   的一個分割   ,對於每個這樣的分割都可定義下列和式

 

隨著區間分割的精細化,該和式可均方收斂於一個   的隨機變量,稱為   關於    上的斯特拉托諾維奇積分[1],記作

 

計算

普通微積分的許多積分技術都可用於斯特拉托諾維奇積分,例如:設   是一個光滑函數,那麼

 

更一般地,對於光滑函數  

 

後一條規則類似於普通微積分的鏈式法則。

數值方法

隨機積分很少能以解析形式求解,這使得隨機數值積分成為隨機積分的各種運用中的一個重要話題。有各種各樣的數值近似都收斂到斯特拉托諾維奇積分;它們的變體用於求解斯特拉托諾維奇隨機微分方程(Kloeden & Platen 1992) 。但注意,求解朗之萬方程的數值解時最為廣泛使用的歐拉-丸山方法要求方程為伊藤形式。[2]

微分記號

若有隨機過程   滿足

 

可以寫出如下的式子

 

這種記號通常用於表述隨機微分方程,它們實際上是關於隨機積分的方程。它兼容於普通微積分的符號,例如

 

與伊藤積分的比較

隨機過程   關於維納過程  伊藤積分記作  (不帶圓圈)。在定義它時,用的是和上文中定義斯特拉托諾維奇積分時相同的步驟,只不過是在每個子區間的左端點取   的值,也就是說

  取代  

與斯特拉托諾維奇積分不同,該積分不遵循普通的鏈式法則,而是須改用稍複雜些的伊藤引理

用以下公式便可在伊藤積分和斯特拉托諾維奇積分之間的進行轉換

 

其中   可以是任意的二元連續可微函數,後一個積分是伊藤積分(Kloeden & Platen 1992,第101頁)。

指明問題是斯特拉托諾維奇還是伊藤表述的重要性可從朗之萬方程的例子看出。假定   是一個時間均勻伊藤擴散,其擴散係數   連續可微,也就是說   滿足隨機微分方程   。為了得到該方程的斯特拉托諾維奇版本,項   應當用   表達出來:

 

顯然,若   是一個常函數(   獨立於   ),朗之萬方程的兩種表述將給出相同的形式。在這種情況下,噪聲項被稱為是「加性」的(因為乘在噪聲項   上的僅是一個固定係數)。否則,伊藤形式的朗之萬方程通常可能與斯特拉托諾維奇形式的朗之萬方程不同,後者的噪聲項稱為乘性的(即,噪聲項中的   乘上了一個  的函數   )。

更一般地,對於任何兩個半鞅   

 

其中  協變差的連續部分。

斯特拉托諾維奇積分的應用

伊藤積分有一項重要性質:不涉及未來的信息。而斯特拉托諾維奇積分則不是這樣。在現實世界的許多應用中(如股票價格建模),人們僅擁有有關過去事件的信息,因此伊藤表述更為自然。在金融數學中通常使用伊藤表述。

然而,在物理學中,隨機積分是作為朗之萬方程的解出現的。朗之萬方程是更微觀模型(Risken 1996)粗粒度版本;關於斯特拉托諾維奇或伊藤表述乃至其他更奇特的解釋(例如等溫表述)何者更為合適,要看所考慮的問題是什麼。物理科學中最常用的是斯特拉托諾維奇表述。


由於斯特拉托諾維奇微積分滿足常規的鏈式法則,斯特拉托諾維奇表述的隨機微分方程可以更直接地定義在可微流形上(而不僅是在   上)。伊藤微積分的複雜鏈式法則則對於微分流形而言更為尷尬。

文內引注

  1. ^ 參見 Gardiner (2004) 的第98頁,以及第101頁的注釋
  2. ^ Perez-Carrasco R.; Sancho J.M. Stochastic algorithms for discontinuous multiplicative white noise (PDF). Phys. Rev. E. 2010, 81 (3): 032104 [2024-08-06]. Bibcode:2010PhRvE..81c2104P. PMID 20365796. doi:10.1103/PhysRevE.81.032104. (原始內容存檔 (PDF)於2018-07-22). 

參考資料