隨機微分方程式(英語:SDE, stochastic differential equation),是常微分方程式的擴展,其項是隨機過程,解也是隨機過程。[1]其形容一個隨機變數的變動過程,也就是常微分方程式加上一個白噪音項[2]。一般微分方程式的對象為可導函數,並以其建立等式。然而,隨機過程函數本身的導數不可定義,所以一般解微分方程式的概念不適用於隨機微分方程式。SDE在純數學中有許多應用,可用於模擬隨機模型的各種行為,如股價[3]、隨機增長模型[4]或受熱漲落影響的物理系統。
SDE具有隨機的微分,在最基本的情形下是以布朗運動導數計算的白噪聲,更一般地說是半鞅。不過,也可能存在其他隨機行為,如萊維過程等跳躍過程[5]或有跳躍的半鞅。隨機微分方程式還可擴展到微分流形。[6][7][8][9]
隨機微分方程式的概念最早以布朗運動的形式,由阿爾伯特·愛因斯坦在《熱的分子運動論所要求的靜液體中懸浮粒子的運動》論文中提出。這項研究隨後由保羅·朗之萬繼續。此後伊藤清和魯斯蘭斯特拉托諾維奇完善了隨機微分方程式的數學基礎,使得這門領域更加的科學嚴謹。
一般而言,隨機微分方程式的解是一隨機過程函數,但解方程式需要先定義隨機過程函數的微分。最常見的定義為根據伊藤清所創,假設B為布朗運動,則對於某函數H,作以下定積分之定義:
此稱為伊藤積分。伊藤式的隨機微分方程式常用於在金融數學中。
隨機微分方程式源於愛因斯坦和Marian Smoluchowski提出的布朗運動理論(1905),不過Louis Bachelier是第一個建立布朗運動模型的人(1900),給出了一個非常早期的SDE實例,即現在所謂Bachelier模型。一些早期例子是線性SDE,也稱為郎之萬方程式,得名於法國物理學家保羅·郎之萬,描述了諧振子在隨機力作用下的運動。1940年代,日本數學家伊藤清發展了SDE的數學理論,提出了隨機分析的概念,並開啟了非線性隨機微分方程式的研究。後來,蘇聯物理學家魯斯蘭·斯特拉托諾維奇提出了另一種方法,產生了類似於普通微積分的隨機積分。
文獻中,SDE最常見的形式是常微分方程式,右式由一個取決於白噪音變量的擾動項。大多數時候SDE被理解為相應隨機差分方程式的連續時間極限,這種SDE理解是模糊的,必須輔以相應積分的適當數學定義。[1][4]這種數學定義由伊藤清提出於1940年代,產生了伊藤積分。後來,蘇聯物理學家魯斯蘭·斯特拉托諾維奇提出了另一種構造,即所謂隨機積分,與伊藤積分是相關但不同的物件,選擇取決於具體應用。伊藤積分以非預期性或因果性概念為基礎,這在以時間為變量的應用中很自然。而隨機積分的規則則與普通微積分相似,且具有內在的幾何特性,使它在處理流形上的隨機運動等問題時更自然。儘管通過伊藤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數值方法的文獻所用的符號。這種記法使物理公式中時間隨機函數的奇異特性更加明確。從嚴格的數學角度來說,不能作為普通函數來選擇,而只能作為廣義函數。數學公式在處理這問題時,比物理公式明確一些。
典型方程式的形式是
其中表示維納過程(標準布朗運動)。方程式被解釋為表達相應積分方程式(下式)的一種非正式方式。
上式將連續時間隨機過程Xt的行為表為普通勒貝格積分與伊藤積分的和。對SDE的啟發法(但非常有用)的解釋是,在長為δ的微小時間區間內,隨機過程Xt值的變化遵循期望值為μ(Xt, t) δ 、變異數為σ(Xt, t)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理論擴展到可微流形上,並使用斯特拉托諾維奇積分。考慮流形、某個有限維向量空間、過濾機率空間,其中滿足通常條件,並令為單點緊化,且為可測。則上的隨機微分方程式可寫為
是一對使得
- 是連續值半鞅;
- 是上向量叢的同態。
,映射都是線性的, for each 。
初始條件為的上的SDE的解是連續的適應的值過程(壽命),且滿足:對每個檢定函數,過程都是實值半鞅;對每個停止時間,方程式
有的把握成立,其中是處的微分形式。若壽命最大,它就是最大解,即
有把握成立,由對每個檢定函數都是半鞅可知,是上的半鞅。給定最大解,可以將的時間擴展到全部,然後在上延拓,可得
(不可分過程)。[20]
雖然斯特拉托諾維奇 SDE滿足連鎖律,其漂移、擴散係數在坐標變化時表現為向量場,因此是流形上SDE的自然選擇,但有時伊藤積分更可取。流形伊藤積分理論是首先由洛朗·施瓦次通過施瓦次同態的概念提出,[7]另見基於射流叢的流形上伊藤SDE的2-射流解釋。[9]當試圖用給定空間上的SDE解與給定空間子曲面上的SDE解進行最佳近似時,這種解釋很有幫助,[10]因為基於斯特拉托諾維奇的投影達不到最佳效果。這已被應用於濾波問題,從而產生了最優投影濾波器。[10]
SDE的求解需要機率設置,因為求解中隱含的積分是隨機積分。若能逐路徑處理微分方程式,就不需要定義了,也就可以發展出獨立於機率論的理論。
這就需要考慮SDE
,都是唯一確定的微分方程式,其中是給定機率空間()的樣本空間。然而,從路徑上直接解釋SDE是不可能的,因為布朗運動路徑的變化無界且無處可微的機率為1,因此沒有簡單方法賦予之類的項以直觀意義。這也排除了將隨機積分定義為對每個的積分的簡單路徑定義。不過,受Wong-Zakai結果[21]對規則噪聲的SDE解的極限的啟發,利用粗糙路徑理論,同時添加布朗運動迭代積分的選定定義,有可能為每個定義確定的粗糙積分,例如,若特定選擇迭代積分可實現與伊藤積分重合的機率為1。[21]迭代積分的其他定義產生不同隨機積分的確定性路徑等價,如斯特拉托諾維奇積分。在金融數學中,這被用於無機率期權的定價。[22]
與定微分方程式一樣,重要的是知道給定的SDE解的存在性和唯一性。下面是在n維歐氏空間Rn中取值,並由m維布朗運動B驅動的伊藤SDE的典型存在性與唯一性定理;證明可見Øksendal (2003, §5.2)。[4]
令T > 0,並使
為可測函數,存在常數C、D使
對所有t ∈ [0, T]、所有x與y ∈ Rn,其中
令Z是獨立於由Bs(s ≥ 0)生成的σ代數的隨機變數,且有有限二階矩:
則隨機微分方程式/初值問題
具有P-幾乎必然獨特t連續解(t, ω) ↦ Xt(ω),使X適應由Z、Bs(s ≤ t)生成的濾子FtZ;另外
上述SDE只是更一般形式的特例:
其中
- 是上的連續半鞅;是上的連續半鞅
- 是從某非空開集發出的映射,其中到的所有線性映射的空間
更一般地說,還可以研究流形上的SDE。
方程式的解收不收斂取決於的選擇。假設滿足某局部利普希茨條件,即對和緊集、常數,滿足條件
其中是歐氏範數。這一條件保證了所謂最大解的存在性和唯一性。
設連續、滿足上述局部利普希茨條件,並設為初始條件,即是關於初始σ-代數的可測函數。令為可預測停時,幾乎確定。值半鞅即下式方程式的最大解
若:
- 對,停止過程是下式停SDE的解:
- 在集合上,幾乎可以確定[23]
則稱為壽命或所謂「爆炸時間」。
主條目:隨機動力的超對稱理論
在SDE的超對稱理論中,隨機動力是通過作用於模型相空間微分形式的隨機演化算子定義的。在這一精確表述中,所有SDE都具有拓撲超對稱性,即通過連續的時間流保持相空間連續性。這種超對稱的自發破缺是混沌、湍流、自組織臨界性等諸多動力現象的數學本質,而南部定理則解釋了相關的長距動力行為,如蝴蝶效應、粉紅噪聲、爆裂聲,以及地震、神經震盪、太陽耀斑等現象的無標統計等等。
Musiela, M., and Rutkowski, M. (2004), Martingale Methods in Financial Modelling, 2nd Edition, Springer Verlag, Berlin.
Øksendal, Bernt K. Stochastic Differential Equations: An Introduction with Applications. Berlin: Springer. 2003. ISBN 3-540-04758-1.
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
Zdzisław Brzeźniak and K. D. Elworthy, Stochastic differential equations on Banach manifolds, Methods Funct. Anal. Topology 6 (2000), no. 1, 43-84.
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.
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 (英語).
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.
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
Steven Marcus, Modeling and approximation of stochastic differential equation driven by semimartigales, Stochastics 4, 1981, 4: 223–245
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 (德語).
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.
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).