隨機微分方程

隨機微分方程(英語:SDE, stochastic differential equation),是常微分方程的擴展,其項是隨機過程,解也是隨機過程。[1]其形容一個隨機變量的變動過程,也就是常微分方程加上一個白噪音[2]。一般微分方程的對象為可導函數,並以其建立等式。然而,隨機過程函數本身的導數不可定義,所以一般解微分方程的概念不適用於隨機微分方程。SDE在純數學中有許多應用,可用於模擬隨機模型的各種行為,如股價[3]、隨機增長模型[4]或受熱漲落影響的物理系統。

SDE具有隨機的微分,在最基本的情形下是以布朗運動導數計算的白噪聲,更一般地說是半鞅。不過,也可能存在其他隨機行為,如萊維過程跳躍過程[5]或有跳躍的半鞅。隨機微分方程還可擴展到微分流形[6][7][8][9]

隨機微分方程的概念最早以布朗運動的形式,由阿爾伯特·愛因斯坦在《熱的分子運動論所要求的靜液體中懸浮粒子的運動》論文中提出。這項研究隨後由保羅·朗之萬繼續。此後伊藤清魯斯蘭斯特拉托諾維奇英語Ruslan Stratonovich完善了隨機微分方程的數學基礎,使得這門領域更加的科學嚴謹。 一般而言,隨機微分方程的解是一隨機過程函數,但解方程需要先定義隨機過程函數的微分。最常見的定義為根據伊藤清所創,假設B布朗運動,則對於某函數H,作以下定積分之定義:

此稱為伊藤積分。伊藤式的隨機微分方程常用於在金融數學中。

背景

編輯

隨機微分方程源於愛因斯坦和Marian Smoluchowski提出的布朗運動理論(1905),不過Louis Bachelier是第一個建立布朗運動模型的人(1900),給出了一個非常早期的SDE實例,即現在所謂Bachelier模型。一些早期例子是線性SDE,也稱為郎之萬方程,得名於法國物理學家保羅·郎之萬,描述了諧振子在隨機力作用下的運動。1940年代,日本數學家伊藤清發展了SDE的數學理論,提出了隨機分析的概念,並開啟了非線性隨機微分方程的研究。後來,蘇聯物理學家魯斯蘭·斯特拉托諾維奇英語Ruslan Stratonovich提出了另一種方法,產生了類似於普通微積分的隨機積分

術語

編輯

文獻中,SDE最常見的形式是常微分方程,右式由一個取決於白噪音變量的擾動項。大多數時候SDE被理解為相應隨機差分方程的連續時間極限,這種SDE理解是模糊的,必須輔以相應積分的適當數學定義。[1][4]這種數學定義由伊藤清提出於1940年代,產生了伊藤積分。後來,蘇聯物理學家魯斯蘭·斯特拉托諾維奇英語Ruslan Stratonovich提出了另一種構造,即所謂隨機積分,與伊藤積分是相關但不同的物件,選擇取決於具體應用。伊藤積分以非預期性或因果性概念為基礎,這在以時間為變量的應用中很自然。而隨機積分的規則則與普通微積分相似,且具有內在的幾何特性,使它在處理流形上的隨機運動等問題時更自然。儘管通過伊藤SDE來模擬流形上的隨機運動也是可能的,且有時更可取[7],例如在試圖最佳化逼近子流形上的SDE時。[10]

SDE的另一種觀點是微分同胚的隨機流,這種理解十分確定,相當於隨機差分方程連續時間極限的斯特拉托諾維奇版本。與SDE相關的是Smoluchowski equation福克-普朗克方程,是描述概率分佈函數隨時間演化的方程。隨機演化算符的概念將福克-普朗克演化推廣為微分形式的時間演化。

在物理科學中,「郎之萬SDE」存在歧義:可以是更一般的形式,但通常指一類具有梯度流向量場的狹義SDE。這類SDE很受歡迎,是Parisi–Sourlas隨機量子化過程的起點,[11]產生了與超對稱量子力學密切相關的N=2超對稱模型。但從物理角度來看,這類SDE不怎麼有趣,因為從未表現出拓撲超對稱性的自發破缺:(過阻尼)郎之萬SDE永不混沌。

隨機分析

編輯

人們發現布朗運動維納過程在數學上異常複雜:維納過程幾乎肯定不可微,[1][4]因此要有自己的分析規則。隨機分析有伊藤積分隨機積分兩個版本,各有利弊,初學者往往搞不清楚特定問題用哪個更合適。有些指南(e.g. Øksendal, 2003)[4],人們可以很方便地將伊藤SDE轉換為等價的隨機SDE,反之亦然。[1][4]不過,最初寫下SDE時還是要決定使用哪種積分。

數值求解

編輯

解SDE的數值方法[12]歐拉-丸山法米爾斯坦法龍格-庫塔法

物理學中的應用

編輯

物理學中,SDE具有廣泛的適用性,從分子動力學到神經動力學,再到天體動力學,不一而足。更具體地說,SDE描述了所有動態系統,其中量子效應要麼不重要,要麼可以作為擾動加以考慮。SDE可被視為動態系統理論對有噪模型的一種概括,這很重要,因為實際系統不可能與其環境完全隔離,總會受外部隨機影響。

有一些標準化手段可通過引入新的未知數,將高階方程轉化為耦合的一階方程組。以下是最常見的一類SDE:

 

其中 是系統在相(或狀態)空間中的位置,設 是可微流形, 是表示確定演化規律的流向量場, 是一組向量場,定義了系統與高斯白噪聲 的耦合。若 是線性空間, 是常數,則稱系統受加性噪聲影響,否則稱系統受乘性噪聲影響。這個術語有點誤導性,它看上去像是有 的約束,實際上意味着一般情形。

對於固定的噪聲,SDE存在唯一解,對初始條件可微。[13]當將各種研究物件平均到噪聲配置上時,隨機情形的非平凡性便顯露出來。乘性噪聲SDE被理解為隨機差分方程的連續時間極限時,並不是唯一確定的實體,而必須輔以所謂「SDE詮釋」,如伊藤或斯特拉托諾維奇解釋。然而,當把SDE看做微分同胚的連續時間隨機流時,則成了唯一確定的數學物件,相當於隨機差分方程連續時間極限的斯特拉托諾維奇法。

物理學中,主要的求解方法是利用等效的福克-普朗克方程(FPE)求出時間函數的概率分佈函數。福克-普朗克方程是確定的偏微分方程,給出了概率分佈函數隨時間的演變,類似於薛定諤方程給出了量子波函數隨時間的演變,或擴散方程給出了化學濃度隨時間的演變。此外還可用蒙地卡羅方法獲得數值解。其他手段還有利用統計物理學和量子力學之間的類比關係進行路徑積分(例如,福克-普朗克方程可通過縮放幾個變量,轉換為薛定諤方程),或寫下概率分佈函數常微分方程[來源請求]

在概率論和金融數學中的應用

編輯

概率論(及概率論中的許多應用,如信號處理中的[濾波問題]]和金融數學中的應用)使用的符號略有不同。它也是解決SDE數值方法的文獻所用的符號。這種記法使物理公式中時間隨機函數 的奇異特性更加明確。從嚴格的數學角度來說, 不能作為普通函數來選擇,而只能作為廣義函數。數學公式在處理這問題時,比物理公式明確一些。

典型方程的形式是

 

其中 表示維納過程(標準布朗運動)。方程被解釋為表達相應積分方程(下式)的一種非正式方式。

 

上式將連續時間隨機過程Xt的行為表為普通勒貝格積分伊藤積分的和。對SDE的啟發法(但非常有用)的解釋是,在長為δ的微小時間區間內,隨機過程Xt值的變化遵循期望值μ(Xttδ方差σ(Xtt)2 δ正態分佈,且與過程過去的行為無關。這是因為,維納過程的增量是獨立的正態分佈。函數μ稱為漂移系數,σ稱為擴散系數。隨機過程Xt稱為擴散過程,滿足馬可夫性質[1]

SDE的形式解釋可據SDE解的構成給出。SDE的解主要有兩種定義,有強解和弱解[1],都要求存在一個能解SDE積分方程形式的過程Xt。兩者的差別在於依賴的概率空間( ),弱解包括一個概率空間和滿足積分方程的過程,強解則包含滿足方程並定義在給定概率空間上的過程。

一個重要例子是幾何布朗運動方程

 

即金融數學中布萊克-舒爾斯模型中的股價動態方程。[3] 推廣幾何布朗運動還可定義允許強解的SDE,分佈是來自不同幾何布朗運動或布萊克-舒爾斯模型的密度的凸組合,從而得到一個單一的SDE,其解的分佈是不同布萊克-舒爾斯模型的對數正態分佈的混合動力。[3][14][15][16]這就產生了可以處理金融數學中所謂波動性微笑的模型。

更簡單的SDE被稱為算術布朗運動[4]

 

Louis Bachelier在1900年將其作為第一個股價模型,即今天所謂Bachelier模型。

還有一些更一般的隨機微分方程,其中的系數μσ不僅取決於過程Xt的現值,還取決於過程的前值,還可能取決於其他過程的現值或前值。這樣,解過程X便不是馬可夫過程或擴散過程,而稱為伊藤過程。當系數只取決於X的現值和前值時,定義方程稱為隨機延遲微分方程。

將帶斯特拉托諾維奇積分的SDE推廣到帶跳躍的半鞅的是馬庫斯型SDE。馬庫斯積分是McShane隨機積分的推廣。[17]

奧恩斯坦-烏倫貝克過程過程方程在隨機金融學中有創新應用:

 

是在收益率呈對數正態分佈的條件下,股價收益率的動態方程。在此假設下,Marcello Minenna開發的方法確定了預測區間,能識別可能隱藏市場濫用現象的異常收益。[18] [19]

流形上的SDE

編輯

更一般地,可以將SDE理論擴展到可微流形上,並使用斯特拉托諾維奇積分。考慮流形 、某個有限維向量空間 、過濾概率空間 ,其中 滿足通常條件,並令 為單點緊化,且  可測。則 上的隨機微分方程可寫為

 

是一對 使得

  •  是連續 值半鞅;
  •   向量叢的同態。

 ,映射 都是線性的,  for each

初始條件為  上的SDE的解是連續的 適應的 值過程 (壽命 ),且滿足:對每個檢驗函數 ,過程 都是實值半鞅;對每個停止時間 ,方程

 

 的把握成立,其中  處的微分形式。若壽命最大,它就是最大解,即

 

 把握成立,由 對每個檢驗函數 都是半鞅可知,  上的半鞅。給定最大解,可以將 的時間擴展到全部 ,然後在 上延拓 ,可得

 

(不可分過程)。[20]

雖然斯特拉托諾維奇 SDE滿足連鎖律,其漂移、擴散系數在坐標變化時表現為向量場,因此是流形上SDE的自然選擇,但有時伊藤積分更可取。流形伊藤積分理論是首先由洛朗·施瓦次通過施瓦次同態的概念提出,[7]另見基於射流叢的流形上伊藤SDE的2-射流解釋。[9]當試圖用給定空間上的SDE解與給定空間子曲面上的SDE解進行最佳近似時,這種解釋很有幫助,[10]因為基於斯特拉托諾維奇的投影達不到最佳效果。這已被應用於濾波問題,從而產生了最優投影濾波器。[10]

作為粗糙路徑的SDE

編輯

SDE的求解需要概率設置,因為求解中隱含的積分是隨機積分。若能逐路徑處理微分方程,就不需要定義了,也就可以發展出獨立於概率論的理論。 這就需要考慮SDE

 

 ,都是唯一確定的微分方程,其中 是給定概率空間( )的樣本空間。然而,從路徑上直接解釋SDE是不可能的,因為布朗運動路徑的變化無界且無處可微的概率為1,因此沒有簡單方法賦予 之類的項以直觀意義。這也排除了將隨機積分定義為對每個 的積分的簡單路徑定義。不過,受Wong-Zakai結果[21]對規則噪聲的SDE解的極限的啟發,利用粗糙路徑理論,同時添加布朗運動迭代積分的選定定義,有可能為每個 定義確定的粗糙積分,例如,若特定選擇迭代積分可實現與伊藤積分重合的概率為1。[21]迭代積分的其他定義產生不同隨機積分的確定性路徑等價,如斯特拉托諾維奇積分。在金融數學中,這被用於無概率期權的定價。[22]

解的存在性與唯一性

編輯

與定微分方程一樣,重要的是知道給定的SDE解的存在性和唯一性。下面是在n維歐氏空間Rn中取值,並由m維布朗運動B驅動的伊藤SDE的典型存在性與唯一性定理;證明可見Øksendal (2003, §5.2)。[4]

T > 0,並使

 
 

可測函數,存在常數CD使

 
 

對所有t ∈ [0, T]、所有xy ∈ Rn,其中

 

Z是獨立於由Bss ≥ 0)生成的σ代數的隨機變量,且有有限二階矩:

 

則隨機微分方程/初值問題

 
 

具有P-幾乎必然獨特t連續解(tω) ↦ Xt(ω),使X適應ZBs(s ≤ t)生成的濾子FtZ;另外

 

一般情形:局部利普希茨條件與最大解

編輯

上述SDE只是更一般形式的特例:

 

其中

  •   上的連續半鞅;  上的連續半鞅
  •  是從某非空開集 發出的映射,其中   的所有線性映射的空間

更一般地說,還可以研究流形上的SDE。

方程的解收不收斂取決於 的選擇。假設 滿足某局部利普希茨條件,即對 和緊集 、常數 ,滿足條件

 

其中 是歐氏範數。這一條件保證了所謂最大解的存在性和唯一性。

 連續、滿足上述局部利普希茨條件,並設 為初始條件,即是關於初始σ-代數的可測函數。令 可預測停時 幾乎確定。 值半鞅 即下式方程的最大解

 

若:

  •  ,停止過程 是下式停SDE的解:
 
  • 在集合 上,幾乎可以確定 [23]

 稱為壽命或所謂「爆炸時間」。

能寫出解析解的SDE[12]

編輯

線性SDE的一般解

編輯
 
 

其中

 

可約SDE: Case 1

編輯
 

對於一個可微函數    等價於 斯特拉托諾維奇 SDE

 

其有一般解

 

其中的 如下所示:

 

可約SDE: Case 2

編輯
 

對於一個可微函數    等價於 斯特拉托諾維奇 SDE

 

可以簡化成下列形式:

 

其中   ,這裏的  如上面所定義。他的一般解可以寫成

 

SDE與超對稱

編輯

在SDE的超對稱理論中,隨機動力是通過作用於模型相空間微分形式的隨機演化算子定義的。在這一精確表述中,所有SDE都具有拓撲超對稱性,即通過連續的時間流保持相空間連續性。這種超對稱的自發破缺是混沌湍流自組織臨界性等諸多動力現象的數學本質,而南部定理則解釋了相關的長距動力行為,如蝴蝶效應粉紅噪聲爆裂聲,以及地震、神經震盪、太陽耀斑等現象的無標統計等等。

相關條目

編輯

參考文獻

編輯
  1. ^ 1.0 1.1 1.2 1.3 1.4 1.5 Rogers, L.C.G.; Williams, David. Diffusions, Markov Processes and Martingales, Vol 2: Ito Calculus 2nd ed., Cambridge Mathematical Library. Cambridge University Press. 2000. ISBN 0-521-77594-9. OCLC 42874839. doi:10.1017/CBO9780511805141. 
  2. ^ 陳, 松男. 金融數學與隨機微積分. ISBN 9789574145089. 
  3. ^ 3.0 3.1 3.2 Musiela, M., and Rutkowski, M. (2004), Martingale Methods in Financial Modelling, 2nd Edition, Springer Verlag, Berlin.
  4. ^ 4.0 4.1 4.2 4.3 4.4 4.5 4.6 Øksendal, Bernt K. Stochastic Differential Equations: An Introduction with Applications. Berlin: Springer. 2003. ISBN 3-540-04758-1. 
  5. ^ Kunita, H. (2004). Stochastic Differential Equations Based on Lévy Processes and Stochastic Flows of Diffeomorphisms. In: Rao, M.M. (eds) Real and Stochastic Analysis. Trends in Mathematics. Birkhäuser Boston. https://doi.org/10.1007/978-1-4612-2054-1_6
  6. ^ Imkeller, Peter; Schmalfuss, Björn. The Conjugacy of Stochastic and Random Differential Equations and the Existence of Global Attractors. Journal of Dynamics and Differential Equations. 2001, 13 (2): 215–249 [2023-10-15]. ISSN 1040-7294. S2CID 3120200. doi:10.1023/a:1016673307045. (原始內容存檔於2023-03-07). 
  7. ^ 7.0 7.1 7.2 Michel Emery (1989). Stochastic calculus in manifolds. Springer Berlin, Heidelberg. Doi https://doi.org/10.1007/978-3-642-75051-9
  8. ^ Zdzisław Brzeźniak and K. D. Elworthy, Stochastic differential equations on Banach manifolds, Methods Funct. Anal. Topology 6 (2000), no. 1, 43-84.
  9. ^ 9.0 9.1 Armstrong J. and Brigo D. (2018). Intrinsic stochastic differential equations as jets. Proc. R. Soc. A., 474: 20170559, http://doi.org/10.1098/rspa.2017.0559
  10. ^ 10.0 10.1 10.2 Armstrong, J., Brigo, D. and Rossi Ferrucci, E. (2019), Optimal approximation of SDEs on submanifolds: the Itô-vector and Itô-jet projections. Proc. London Math. Soc., 119: 176-213. https://doi.org/10.1112/plms.12226.
  11. ^ Parisi, G.; Sourlas, N. Random Magnetic Fields, Supersymmetry, and Negative Dimensions. Physical Review Letters. 1979, 43 (11): 744–745. Bibcode:1979PhRvL..43..744P. doi:10.1103/PhysRevLett.43.744. 
  12. ^ 12.0 12.1 Peter E. Kloeden, Eckhard Platen (1992). Numerical Solution of Stochastic Differential Equations. Springer Berlin, Heidelberg. DOI: https://doi.org/10.1007/978-3-662-12616-5
  13. ^ Slavík, A. Generalized differential equations: Differentiability of solutions with respect to initial conditions and parameters. Journal of Mathematical Analysis and Applications. 2013, 402 (1): 261–274. doi:10.1016/j.jmaa.2013.01.027  (英語). 
  14. ^ Fengler, M. R. (2005), Semiparametric modeling of implied volatility, Springer Verlag, Berlin. DOI https://doi.org/10.1007/3-540-30591-2
  15. ^ Brigo, Damiano; Mercurio, Fabio. Lognormal-mixture dynamics and calibration to market volatility smiles. International Journal of Theoretical and Applied Finance. 2002, 5 (4): 427–446. doi:10.1142/S0219024902001511. 
  16. ^ Brigo, D, Mercurio, F, Sartorelli, G. (2003). Alternative asset-price dynamics and volatility smile, QUANT FINANC, 2003, Vol: 3, Pages: 173 - 183, ISSN 1469-7688
  17. ^ Steven Marcus, Modeling and approximation of stochastic differential equation driven by semimartigales, Stochastics 4, 1981, 4: 223–245 
  18. ^ Detecting Market Abuse. Risk Magazine. [2023-10-15]. (原始內容存檔於2023-05-18). 
  19. ^ The detection of Market Abuse on financial markets: a quantitative approach (PDF). Consob – The Italian Securities and Exchange Commission. [2023-10-15]. (原始內容存檔 (PDF)於2023-03-21). 
  20. ^ Hackenbroch, Wolfgang; Thalmaier, Anton. Stochastische Analysis: Eine Einführung in die Theorie der stetigen Semimartingale. Vieweg+Teubner Verlag Wiesbaden. 1994: 364-365. ISBN 978-3-519-02229-9 (德語). 
  21. ^ 21.0 21.1 Friz, P. and Hairer, M. (2020). A Course on Rough Paths with an Introduction to Regularity Structures, 2nd ed., Springer-Verlag, Heidelberg, DOI https://doi.org/10.1007/978-3-030-41556-3
  22. ^ Armstrong, J., Bellani, C., Brigo, D. and Cass, T. (2021). Option pricing models without probability: a rough paths approach. Mathematical Finance, vol. 31, pages 1494–1521.
  23. ^ Hackenbroch, Wolfgang; Thalmaier, Anton. Stochastische Analysis: Eine Einführung in die Theorie der stetigen Semimartingale. Vieweg+Teubner Verlag Wiesbaden. 1994: 297–299. ISBN 978-3-519-02229-9 (德語). 

閱讀更多

編輯
  • Adomian, George. Stochastic systems. Mathematics in Science and Engineering (169). Orlando, FL: Academic Press Inc. 1983. 
  • Adomian, George. Nonlinear stochastic operator equations . Orlando, FL: Academic Press Inc. 1986. 
  • Adomian, George. Nonlinear stochastic systems theory and applications to physics. Mathematics and its Applications (46). Dordrecht: Kluwer Academic Publishers Group. 1989. 
  • Calin, Ovidiu. An Informal Introduction to Stochastic Calculus with Applications. Singapore: World Scientific Publishing. 2015: 315. ISBN 978-981-4678-93-3. 
  • Teugels, J. and Sund B. (eds.). Encyclopedia of Actuarial Science. Chichester: Wiley. 2004: 523–527. 
  • C. W. Gardiner. Handbook of Stochastic Methods: for Physics, Chemistry and the Natural Sciences. Springer. 2004: 415. 
  • Thomas Mikosch. Elementary Stochastic Calculus: with Finance in View. Singapore: World Scientific Publishing. 1998: 212. ISBN 981-02-3543-7. 
  • Seifedine Kadry. A Solution of Linear Stochastic Differential Equation. Wseas Transactions on Mathematics (USA: WSEAS TRANSACTIONS on MATHEMATICS, April 2007.). 2007: 618. ISSN 1109-2769. 
  • Higham., Desmond J. An Algorithmic Introduction to Numerical Simulation of Stochastic Differential Equations. SIAM Review. January 2001, 43 (3): 525–546. Bibcode:2001SIAMR..43..525H. CiteSeerX 10.1.1.137.6375 . doi:10.1137/S0036144500378302. 
  • Desmond Higham and Peter Kloeden: "An Introduction to the Numerical Simulation of Stochastic Differential Equations", SIAM, ISBN 978-1-611976-42-7 (2021).