查看“︁斯特拉托诺维奇积分”︁的源代码

在[[随机过程]]中，'''斯特拉托诺维奇积分'''或'''菲斯克-斯特拉托诺维奇积分'''是一种[[随机分析|随机积分]]，冠名于同时发展出该技术的{{Le|鲁斯兰·斯特拉托诺维奇|Ruslan Stratonovich}}和[[唐纳德·菲斯克]]。它是[[伊藤积分]]最常见的替代品。应用数学中通常选择伊藤积分，而物理学中则频繁使用斯特拉托诺维奇积分。

在某些场景下，斯特拉托诺维奇积分更易于操作。斯特拉托诺维奇积分的定义使其满足[[链式法则]]，而[[伊藤积分]]则并非如此。

遇到斯特拉托诺维奇积分的一类主要场景是斯特拉托诺维奇[[隨機微分方程|随机微分方程]](SDE)的解。其与伊藤随机微分方程是等价的，并且可以根据何种更为方便而在两种定义之间进行切换。

== 定义 ==
斯特拉托诺维奇积分可以用类似于[[黎曼积分]]的方式（[[黎曼和]]的[[极限 (数学)|极限]]）来定义。设 <math>W : [0, T] \times \Omega \to \mathbb{R}</math> 是一个[[维纳过程]]， <math>X : [0, T] \times \Omega \to \mathbb{R}</math> 是[[适应过程|适应]]于维纳过程的自然{{Le|滤过 (概率论)|Filtration (probability_theory)}}的一个[[半鞅]]。

类似[[黎曼-斯蒂尔杰斯积分]]中的策略，我们考虑区间 <math>[0, T]</math> 的一个分割 <math>0 = t_{0} < t_{1} < \dots < t_{k} = T</math> ，对于每个这样的分割都可定义下列和式

<math display="block">\sum_{i = 0}^{k - 1} {X_{t_{i+1}} + X_{t_i}\over 2} \left( W_{t_{i+1}} - W_{t_i} \right),</math>

随着区间分割的精细化，该和式可[[随机变量的收敛|均方收敛]]于一个 <math>\Omega \to \mathbb{R}</math> 的随机变量，称为 <math>X</math> 关于 <math>W</math> 在 <math>[0,T]</math> 上的'''斯特拉托诺维奇积分'''<ref>参见 Gardiner (2004) 的第98页，以及第101页的注释</ref>，记作

<math display="block">\int_0^T X_{t} \circ \mathrm{d} W_t.</math>

== 计算 ==
普通微积分的许多积分技术都可用于斯特拉托诺维奇积分，例如：设 <math>f : \mathbb{R} \to \mathbb{R}</math> 是一个光滑函数，那么

: <math>\int_0^T f'(W_t) \circ \mathrm{d} W_t = f(W_T)-f(W_0)</math>

更一般地，对于光滑函数 <math>f : \mathbb{R} \times \mathbb{R} \to \mathbb{R}</math> 有

: <math>\int_0^T {\partial f\over\partial W}(W_t,t) \circ \mathrm{d} W_t + \int_0^T {\partial f\over\partial t}(W_t,t)\, \mathrm{d}t = f(W_T,T)-f(W_0,0).</math>

后一条规则类似于普通微积分的链式法则。

=== 数值方法 ===
随机积分很少能以解析形式求解，这使得[[随机分析|随机]][[數值積分|数值积分]]成为随机积分的各种运用中的一个重要话题。有各种各样的数值近似都收敛到斯特拉托诺维奇积分；它们的变体用于求解斯特拉托诺维奇随机微分方程{{Harvard citation|Kloeden|Platen|1992}} 。但注意，[[朗之万方程|求解朗之万方程]]的数值解时最为广泛使用的[[欧拉-丸山法|欧拉-丸山方法]]要求方程为伊藤形式。<ref>{{Cite journal |last=Perez-Carrasco R. |last2=Sancho J.M. |title=Stochastic algorithms for discontinuous multiplicative white noise |url=http://discovery.ucl.ac.uk/10025194/9/PhysRevE.81.032104%20%281%29.pdf |journal=Phys. Rev. E |year=2010 |volume=81 |issue=3 |page=032104 |bibcode=2010PhRvE..81c2104P |doi=10.1103/PhysRevE.81.032104 |pmid=20365796 |access-date=2024-08-06 |archive-date=2018-07-22 |archive-url=https://web.archive.org/web/20180722175203/http://discovery.ucl.ac.uk/10025194/9/PhysRevE.81.032104%20%281%29.pdf |dead-url=no }}</ref>

== 微分记号 ==
若有随机过程 <math>X_t, Y_t, Z_t</math> 满足

<math display="block">\forall T>0,\quad X_T-X_0=\int_0^T Y_{t} \circ \mathrm{d} W_t + \int_0^T Z_{t} \,\mathrm{d}t,</math>

可以写出如下的式子

<math display="block">\mathrm{d}X=Y\circ\mathrm{d}W + Z\,\mathrm{d}t.</math>

这种记号通常用于表述[[隨機微分方程|随机微分方程]]，它们实际上是关于随机积分的方程。它兼容于普通微积分的符号，例如

<math display="block">\mathrm{d}(t^2\,W^3)=3 t^2 W^2\circ\mathrm{d}W + 2t W^3\,\mathrm{d}t.</math>

== 与伊藤积分的比较 ==
随机过程 <math>X</math> 关于维纳过程 <math>W</math> 的[[伊藤积分]]记作<math display="block">\int_0^T X_{t} \,\mathrm{d} W_t</math> （不带圆圈）。在定义它时，用的是和上文中定义斯特拉托诺维奇积分时相同的步骤，只不过是在每个子区间的左端点取 <math>X</math> 的值，也就是说

: 以 <math>X_{t_{i}}</math> 取代 <math>(X_{t_{i+1}}+ X_{t_{i}})/ 2</math>

与斯特拉托诺维奇积分不同，该积分不遵循普通的链式法则，而是须改用稍复杂些的[[伊藤引理]]。

用以下公式便可在伊藤积分和斯特拉托诺维奇积分之间的进行转换

: <math>\int_{0}^{T} f(W_{t},t) \circ \mathrm{d} W_{t} = \frac{1}{2} \int_{0}^{T} {\partial f\over\partial W}(W_{t},t) \, \mathrm{d} t + \int_{0}^{T} f(W_{t},t) \, \mathrm{d} W_{t},</math>

其中 <math>f</math> 可以是任意的二元[[连续可微]]函数，后一个积分是伊藤积分{{Harvard citation|Kloeden|Platen|1992|p=101}}。

指明问题是斯特拉托诺维奇还是伊藤表述的重要性可从朗之万方程的例子看出。假定 <math>X_t</math> 是一个时间均匀伊藤扩散，其扩散系数 <math>\sigma</math> 连续可微，也就是说 <math>X_t</math> 满足随机微分方程 <math>\mathrm{d} X_t = \mu(X_t)\,\mathrm{d} t + \sigma(X_t)\,\mathrm{d} W_t</math> 。为了得到该方程的斯特拉托诺维奇版本，项 <math>\sigma(X_t)\,\mathrm{d} W_t</math> 应当用 <math> \sigma (X_{t}) \circ \mathrm{d} W_{t}</math> 表达出来：

: <math>\int_{0}^{T} \sigma (X_{t}) \circ \mathrm{d} W_{t} = \frac{1}{2} \int_{0}^{T} \sigma'(X_{t}) \sigma(X_{t}) \, \mathrm{d} t + \int_{0}^{T} \sigma (X_{t}) \, \mathrm{d} W_{t}.</math>

显然，若 <math> \sigma </math> 是一个常函数（ <math> \sigma(X_t) </math> 独立于 <math>X_t </math> ），朗之万方程的两种表述将给出相同的形式。在这种情况下，噪声项被称为是“加性”的（因为乘在噪声项 <math> dW_t </math> 上的仅是一个固定系数）。否则，伊藤形式的朗之万方程通常可能与斯特拉托诺维奇形式的朗之万方程不同，后者的噪声项称为乘性的（即，噪声项中的 <math> dW_t </math> 乘上了一个<math> X_t </math> 的函数 <math> \sigma(X_t) </math> ）。

更一般地，对于任何两个[[半鞅]] <math>X</math> 和 <math>Y</math>

: <math>\int_{0}^{T} X_{s-} \circ \mathrm{d} Y_s = \int_0^T X_{s-}\,\mathrm{d}Y_s+ \frac{1}{2} [X,Y]_T^c,</math>

其中 <math> [X,Y]_T^c</math> 是[[二次变差|协变差]]的连续部分。

== 斯特拉托诺维奇积分的应用 ==
伊藤积分有一项重要性质：不涉及未来的信息。而斯特拉托诺维奇积分则不是这样。在现实世界的许多应用中（如股票价格建模），人们仅拥有有关过去事件的信息，因此伊藤表述更为自然。在金融数学中通常使用伊藤表述。

然而，在物理学中，随机积分是作为[[朗之万方程]]的解出现的。朗之万方程是更微观模型{{Harvard citation|Risken|1996}}粗粒度版本；关于斯特拉托诺维奇或伊藤表述乃至其他更奇特的解释（例如等温表述）何者更为合适，要看所考虑的问题是什么。物理科学中最常用的是斯特拉托诺维奇表述。


由于斯特拉托诺维奇微积分满足常规的链式法则，斯特拉托诺维奇表述的随机微分方程可以更直接地定义在[[可微流形]]上（而不仅是在 <math>\mathbb{R}^n</math> 上）。伊藤微积分的复杂链式法则则对于微分流形而言更为尴尬。

== 文内引注 ==
{{Reflist}}{{Reflist}}

== 参考资料 ==

* {{Cite book|last=Øksendal, Bernt K.|authorlink=Bernt Øksendal|title=Stochastic Differential Equations: An Introduction with Applications|publisher=Springer, Berlin|year=2003|isbn=3-540-04758-1}}
* {{Cite book|last=Gardiner, Crispin W.|title=Handbook of Stochastic Methods|publisher=Springer, Berlin Heidelberg|year=2004|edition=3|isbn=3-540-20882-8}}
* {{Cite journal |last=Jarrow |first=Robert |last2=Protter |first2=Philip |title=A short history of stochastic integration and mathematical finance: The early years, 1880–1970 |journal=IMS Lecture Notes Monograph |year=2004 |volume=45 |page=1–17 |citeseerx=10.1.1.114.632}}
* {{Cite book|last=Kloeden|first=Peter E.|last2=Platen|first2=Eckhard|title=Numerical solution of stochastic differential equations|publisher=[[Springer-Verlag]]|location=Berlin, New York|series=Applications of Mathematics|isbn=978-3-540-54062-5|year=1992}}
* {{Cite book|last1=Risken|first1=Hannes|title=The Fokker-Planck Equation|url=https://archive.org/details/fokkerplanckequa0000risk|publisher=[[Springer-Verlag]]|location=Berlin, Heidelberg|series=Springer Series in Synergetics|isbn=978-3-540-61530-9|year=1996}}.
[[Category:积分的定义]]