隨機微分方程 (英語:SDE, stochastic differential equation ),是常微分方程 的擴展,其项是随机过程 ,解也是随机过程。[ 1] 其形容一個隨機變數的變動過程,也就是常微分方程加上一個白噪音 項[ 2] 。一般微分方程的對象為可導函數,並以其建立等式。然而,隨機過程函數本身的導數 不可定義,所以一般解微分方程的概念不適用於隨機微分方程。SDE在纯数学中有许多应用,可用于模拟随机模型 的各种行为,如股价 [ 3] 、随机增长模型[ 4] 或受热涨落 影响的物理系统。
SDE具有随机的微分,在最基本的情形下是以布朗运动 导数计算的白噪声 ,更一般地说是半鞅 。不过,也可能存在其他随机行为,如莱维过程 等跳跃过程 [ 5] 或有跳跃的半鞅。随机微分方程还可扩展到微分流形 。[ 6] [ 7] [ 8] [ 9]
隨機微分方程的概念最早以布朗運動 的形式,由阿爾伯特·愛因斯坦 在《热的分子运动论所要求的静液体中悬浮粒子的运动》論文中提出。这项研究隨後由保羅·朗之萬 继续。此後伊藤清 和鲁斯兰斯特拉托诺维奇 完善了隨機微分方程的數學基礎,使得這門領域更加的科學嚴謹。
一般而言,隨機微分方程的解是一隨機過程 函數,但解方程需要先定義隨機過程函數的微分 。最常見的定義為根據伊藤清 所創,假設B 為布朗運動 ,則對於某函數H ,作以下定積分之定義:
∫
0
t
H
d
B
=
lim
n
→
∞
∑
t
i
−
1
,
t
i
∈
π
n
H
t
i
−
1
(
B
t
i
−
B
t
i
−
1
)
.
{\displaystyle \int _{0}^{t}H\,dB=\lim _{n\rightarrow \infty }\sum _{t_{i-1},t_{i}\in \pi _{n}}H_{t_{i-1}}(B_{t_{i}}-B_{t_{i-1}}).}
此稱為伊藤積分 。伊藤式的隨機微分方程常用於在金融數學 中。
背景
随机微分方程源于爱因斯坦 和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具有广泛的适用性,从分子动力学到神经动力学,再到天体动力学,不一而足。更具体地说,SDE描述了所有动力系统,其中量子效应要么不重要,要么可以作为扰动加以考虑。SDE可被视为动力系统理论 对有噪模型的一种概括,这很重要,因为实际系统不可能与其环境完全隔离,总会受外部随机影响。
有一些标准化手段可通过引入新的未知数,将高阶方程转化为耦合的一阶方程组。以下是最常见的一类SDE:
d
x
(
t
)
d
t
=
F
(
x
(
t
)
)
+
∑
α
=
1
n
g
α
(
x
(
t
)
)
ξ
α
(
t
)
,
{\displaystyle {\frac {\mathrm {d} x(t)}{\mathrm {d} t}}=F(x(t))+\sum _{\alpha =1}^{n}g_{\alpha }(x(t))\xi ^{\alpha }(t),\,}
其中
x
∈
X
{\displaystyle x\in X}
是系统在相(或状态)空间中的位置,设
X
{\displaystyle X}
是可微流形,
F
∈
T
X
{\displaystyle F\in TX}
是表示确定演化规律的流向量场,
g
α
∈
T
X
{\displaystyle g_{\alpha }\in TX}
是一组向量场,定义了系统与高斯白噪声
ξ
α
{\displaystyle \xi ^{\alpha }}
的耦合。若
X
{\displaystyle X}
是线性空间,
g
{\displaystyle g}
是常数,则称系统受加性噪声影响,否则称系统受乘性噪声影响。这个术语有点误导性,它看上去像是有
g
(
x
)
∝
x
{\displaystyle g(x)\propto x}
的约束,实际上意味着一般情形。
对于固定的噪声,SDE存在唯一解,对初始条件可微。[ 13] 当将各种研究对象平均到噪声配置上时,随机情形的非平凡性便显露出来。乘性噪声SDE被理解为随机差分方程 的连续时间极限时,并不是唯一确定的实体,而必须辅以所谓“SDE诠释”,如伊藤或斯特拉托诺维奇解释。然而,当把SDE看做微分同胚的连续时间随机流时,则成了唯一确定的数学对象 ,相当于随机差分方程连续时间极限的斯特拉托诺维奇法。
物理学中,主要的求解方法是利用等效的福克-普朗克方程 (FPE)求出时间函数的概率分布 函数。福克-普朗克方程是确定的偏微分方程 ,给出了概率分布函数随时间的演变,类似于薛定谔方程 给出了量子波函数随时间的演变,或扩散方程 给出了化学浓度随时间的演变。此外还可用蒙特卡洛方法 获得数值解。其他手段还有利用统计物理学和量子力学 之间的类比关系进行路径积分 (例如,福克-普朗克方程可通过缩放几个变量,变换为薛定谔方程 ),或写下概率分布函数 矩 的常微分方程 。[來源請求]
在概率论和金融数学中的应用
概率论 (及概率论中的许多应用,如信号处理中的[滤波问题]]和金融数学 中的应用)使用的符号略有不同。它也是解决SDE数值方法的文献所用的符号。这种记法使物理公式中时间随机函数
η
m
{\displaystyle \eta _{m}}
的奇异特性更加明确。从严格的数学角度来说,
η
m
{\displaystyle \eta _{m}}
不能作为普通函数来选择,而只能作为广义函数 。数学公式在处理这问题时,比物理公式明确一些。
典型方程的形式是
d
X
t
=
μ
(
X
t
,
t
)
d
t
+
σ
(
X
t
,
t
)
d
B
t
,
{\displaystyle \mathrm {d} X_{t}=\mu (X_{t},t)\,\mathrm {d} t+\sigma (X_{t},t)\,\mathrm {d} B_{t},}
其中
B
{\displaystyle B}
表示维纳过程 (标准布朗运动)。方程被解释为表达相应积分方程 (下式)的一种非正式方式。
X
t
+
s
−
X
t
=
∫
t
t
+
s
μ
(
X
u
,
u
)
d
u
+
∫
t
t
+
s
σ
(
X
u
,
u
)
d
B
u
.
{\displaystyle X_{t+s}-X_{t}=\int _{t}^{t+s}\mu (X_{u},u)\mathrm {d} u+\int _{t}^{t+s}\sigma (X_{u},u)\,\mathrm {d} B_{u}.}
上式将连续时间随机过程 X t 的行为表为普通勒贝格积分 与伊藤积分 的和。对SDE的启发法 (但非常有用)的解释是,在长为δ 的微小时间区间内,随机过程X t 值的变化遵循期望 为μ (X t , t ) δ 、方差 为σ (X t , t )2 δ 的正态分布 ,且与过程过去的行为无关。这是因为,维纳过程的增量是独立的正态分布。函数μ 称为漂移系数,σ 称为扩散系数。随机过程X t 称为扩散过程 ,满足马尔可夫性质 。[ 1]
SDE的形式解释可据SDE解的构成给出。SDE的解主要有两种定义,有强解和弱解[ 1] ,都要求存在一个能解SDE积分方程形式的过程X t 。两者的差别在于依赖的概率空间 (
Ω
,
F
,
P
{\displaystyle \Omega ,\,{\mathcal {F}},\,P}
),弱解包括一个概率空间和满足积分方程的过程,强解则包含满足方程并定义在给定概率空间上的过程。
一个重要例子是几何布朗运动 方程
d
X
t
=
μ
X
t
d
t
+
σ
X
t
d
B
t
.
{\displaystyle \mathrm {d} X_{t}=\mu X_{t}\,\mathrm {d} t+\sigma X_{t}\,\mathrm {d} B_{t}.}
即金融数学中布莱克-舒尔斯模型 中的股价动态方程。[ 3]
推广几何布朗运动还可定义允许强解的SDE,分布是来自不同几何布朗运动或布莱克-舒尔斯模型的密度的凸组合,从而得到一个单一的SDE,其解的分布是不同布莱克-舒尔斯模型的对数正态分布的混合动力。[ 3] [ 14] [ 15] [ 16] 这就产生了可以处理金融数学中所谓波动性微笑 的模型。
更简单的SDE被称为算术布朗运动[ 4]
d
X
t
=
μ
d
t
+
σ
d
B
t
{\displaystyle \mathrm {d} X_{t}=\mu \,\mathrm {d} t+\sigma \,\mathrm {d} B_{t}}
Louis Bachelier在1900年将其作为第一个股价模型,即今天所谓Bachelier模型。
还有一些更一般的随机微分方程,其中的系数μ 、σ 不仅取决于过程X t 的现值,还取决于过程的前值,还可能取决于其他过程的现值或前值。这样,解过程X 便不是马尔可夫过程或扩散过程,而称为伊藤过程。当系数只取决于X 的现值和前值时,定义方程称为随机延迟微分方程。
将带斯特拉托诺维奇积分 的SDE推广到带跳跃的半鞅的是马库斯型SDE。马库斯积分是McShane随机积分的推广。[ 17]
奥恩斯坦-乌伦贝克过程 过程方程在随机金融学中有创新应用:
d
R
t
=
μ
R
t
d
t
+
σ
t
d
B
t
.
{\displaystyle \mathrm {d} R_{t}=\mu R_{t}\,\mathrm {d} t+\sigma _{t}\,\mathrm {d} B_{t}.}
是在收益率呈对数正态分布 的条件下,股价收益率的动态方程。在此假设下,Marcello Minenna开发的方法确定了预测区间,能识别可能隐藏市场滥用 现象的异常收益。[ 18] [ 19]
流形上的SDE
更一般地,可以将SDE理论扩展到可微流形 上,并使用斯特拉托诺维奇积分 。考虑流形
M
{\displaystyle M}
、某个有限维向量空间
E
{\displaystyle E}
、过滤概率空间
(
Ω
,
F
,
(
F
t
)
t
∈
R
+
,
P
)
{\displaystyle (\Omega ,{\mathcal {F}},({\mathcal {F}}_{t})_{t\in \mathbb {R} _{+}},P)}
,其中
(
F
t
)
t
∈
R
+
{\displaystyle ({\mathcal {F}}_{t})_{t\in \mathbb {R} _{+}}}
满足通常条件,并令
M
^
=
M
∪
{
∞
}
{\displaystyle {\widehat {M}}=M\cup \{\infty \}}
为单点紧化,且
x
0
{\displaystyle x_{0}}
为
F
0
{\displaystyle {\mathcal {F}}_{0}}
可测。则
M
{\displaystyle M}
上的随机微分方程可写为
d
X
=
A
(
X
)
∘
d
Z
{\displaystyle \mathrm {d} X=A(X)\circ dZ}
是一对
(
A
,
Z
)
{\displaystyle (A,Z)}
使得
Z
{\displaystyle Z}
是连续
E
{\displaystyle E}
值半鞅;
A
:
M
×
E
→
T
M
,
(
x
,
e
)
↦
A
(
x
)
e
{\displaystyle A:M\times E\to TM,(x,e)\mapsto A(x)e}
是
M
{\displaystyle M}
上向量丛 的同态。
∀
x
∈
M
{\displaystyle \forall x\in M}
,映射
A
(
x
)
:
E
→
T
x
M
{\displaystyle A(x):E\to T_{x}M}
都是线性的,
∀
e
∈
E
,
A
(
⋅
)
e
∈
Γ
(
T
M
)
{\displaystyle \forall e\in E,\ A(\cdot )e\in \Gamma (TM)}
for each
{\displaystyle }
。
初始条件为
X
0
=
x
0
{\displaystyle X_{0}=x_{0}}
的
M
{\displaystyle M}
上的SDE的解是连续的
{
F
t
}
{\displaystyle \{{\mathcal {F}}_{t}\}}
适应的
M
{\displaystyle M}
值过程
(
X
t
)
t
<
ζ
{\displaystyle (X_{t})_{t<\zeta }}
(寿命
ζ
{\displaystyle \zeta }
),且满足:对每个检验函数
f
∈
C
c
∞
(
M
)
{\displaystyle f\in C_{c}^{\infty }(M)}
,过程
f
(
X
)
{\displaystyle f(X)}
都是实值半鞅;对每个停止时间
τ
(
0
≤
τ
<
ζ
)
{\displaystyle \tau \ (0\leq \tau <\zeta )}
,方程
f
(
X
τ
)
=
f
(
x
0
)
+
∫
0
τ
(
d
f
)
X
A
(
X
)
∘
d
Z
{\displaystyle f(X_{\tau })=f(x_{0})+\int _{0}^{\tau }(\mathrm {d} f)_{X}A(X)\circ \mathrm {d} Z}
有
P
{\displaystyle P}
的把握成立,其中
(
d
f
)
X
:
T
x
M
→
T
f
(
x
)
M
{\displaystyle (df)_{X}:T_{x}M\to T_{f(x)}M}
是
X
{\displaystyle X}
处的微分形式 。若寿命最大,它就是最大解,即
{
ζ
<
∞
}
⊂
{
lim
t
↗
ζ
X
t
=
∞
in
M
^
}
{\displaystyle \{\zeta <\infty \}\subset \left\{\lim \limits _{t\nearrow \zeta }X_{t}=\infty {\text{ in }}{\widehat {M}}\right\}}
有
P
{\displaystyle P}
把握成立,由
f
(
X
)
{\displaystyle f(X)}
对每个检验函数
f
∈
C
c
∞
(
M
)
{\displaystyle f\in C_{c}^{\infty }(M)}
都是半鞅可知,
X
{\displaystyle X}
是
M
{\displaystyle M}
上的半鞅。给定最大解,可以将
X
{\displaystyle X}
的时间扩展到全部
R
+
{\displaystyle \mathbb {R} _{+}}
,然后在
M
^
{\displaystyle {\widehat {M}}}
上延拓
f
{\displaystyle f}
,可得
f
(
X
t
)
=
f
(
X
0
)
+
∫
0
t
(
d
f
)
X
A
(
X
)
∘
d
Z
,
t
≥
0
{\displaystyle f(X_{t})=f(X_{0})+\int _{0}^{t}(\mathrm {d} f)_{X}A(X)\circ \mathrm {d} Z,\quad t\geq 0}
(不可分过程)。[ 20]
虽然斯特拉托诺维奇 SDE满足链式法则,其漂移、扩散系数在坐标变化时表现为向量场,因此是流形上SDE的自然选择,但有时伊藤积分更可取。流形伊藤积分理论是首先由洛朗·施瓦兹 通过施瓦兹同态的概念提出,[ 7] 另见基于射流丛的流形上伊藤SDE的2-射流解释。[ 9] 当试图用给定空间上的SDE解与给定空间子曲面上的SDE解进行最佳近似时,这种解释很有帮助,[ 10] 因为基于斯特拉托诺维奇的投影达不到最佳效果。这已被应用于滤波问题 ,从而产生了最优投影滤波器。[ 10]
作为粗糙路径的SDE
SDE的求解需要概率设置,因为求解中隐含的积分是随机积分。若能逐路径处理微分方程,就不需要定义了,也就可以发展出独立于概率论的理论。
这就需要考虑SDE
d
X
t
(
ω
)
=
μ
(
X
t
(
ω
)
,
t
)
d
t
+
σ
(
X
t
(
ω
)
,
t
)
d
B
t
(
ω
)
{\displaystyle \mathrm {d} X_{t}(\omega )=\mu (X_{t}(\omega ),t)\,\mathrm {d} t+\sigma (X_{t}(\omega ),t)\,\mathrm {d} B_{t}(\omega )}
∀
ω
∈
Ω
{\displaystyle \forall \omega \in \Omega }
,都是唯一确定的微分方程,其中
Ω
{\displaystyle \Omega }
是给定概率空间(
Ω
,
F
,
P
{\displaystyle \Omega ,\,{\mathcal {F}},\,P}
)的样本空间。然而,从路径上直接解释SDE是不可能的,因为布朗运动路径的变化无界且无处可微的概率为1,因此没有简单方法赋予
d
B
t
(
ω
)
{\displaystyle \mathrm {d} B_{t}(\omega )}
之类的项以直观意义。这也排除了将随机积分定义为对每个
d
B
t
(
ω
)
{\displaystyle \mathrm {d} B_{t}(\omega )}
的积分的简单路径定义。不过,受Wong-Zakai结果[ 21] 对规则噪声的SDE解的极限的启发,利用粗糙路径理论,同时添加布朗运动迭代积分的选定定义,有可能为每个
ω
∈
Ω
{\displaystyle \omega \in \Omega }
定义确定的粗糙积分,例如,若特定选择迭代积分可实现与伊藤积分重合的概率为1。[ 21] 迭代积分的其他定义产生不同随机积分的确定性路径等价,如斯特拉托诺维奇积分。在金融数学中,这被用于无概率期权的定价。[ 22]
解的存在性与唯一性
与定微分方程一样,重要的是知道给定的SDE解的存在性和唯一性。下面是在n维欧氏空间R n 中取值,并由m维布朗运动B 驱动的伊藤SDE的典型存在性与唯一性定理;证明可见Øksendal (2003, §5.2)。[ 4]
令T > 0,并使
μ
:
R
n
×
[
0
,
T
]
→
R
n
;
{\displaystyle \mu :\mathbb {R} ^{n}\times [0,T]\to \mathbb {R} ^{n};}
σ
:
R
n
×
[
0
,
T
]
→
R
n
×
m
;
{\displaystyle \sigma :\mathbb {R} ^{n}\times [0,T]\to \mathbb {R} ^{n\times m};}
为可测函数 ,存在常数C 、D 使
|
μ
(
x
,
t
)
|
+
|
σ
(
x
,
t
)
|
≤
C
(
1
+
|
x
|
)
;
{\displaystyle {\big |}\mu (x,t){\big |}+{\big |}\sigma (x,t){\big |}\leq C{\big (}1+|x|{\big )};}
|
μ
(
x
,
t
)
−
μ
(
y
,
t
)
|
+
|
σ
(
x
,
t
)
−
σ
(
y
,
t
)
|
≤
D
|
x
−
y
|
;
{\displaystyle {\big |}\mu (x,t)-\mu (y,t){\big |}+{\big |}\sigma (x,t)-\sigma (y,t){\big |}\leq D|x-y|;}
对所有t ∈ [0, T ]、所有x 与y ∈ R n ,其中
|
σ
|
2
=
∑
i
,
j
=
1
n
|
σ
i
j
|
2
.
{\displaystyle |\sigma |^{2}=\sum _{i,j=1}^{n}|\sigma _{ij}|^{2}.}
令Z 是独立于由B s (s ≥ 0)生成的σ 代数的随机变量,且有有限二阶矩:
E
[
|
Z
|
2
]
<
+
∞
.
{\displaystyle \mathbb {E} {\big [}|Z|^{2}{\big ]}<+\infty .}
则随机微分方程/初值问题
d
X
t
=
μ
(
X
t
,
t
)
d
t
+
σ
(
X
t
,
t
)
d
B
t
for
t
∈
[
0
,
T
]
;
{\displaystyle \mathrm {d} X_{t}=\mu (X_{t},t)\,\mathrm {d} t+\sigma (X_{t},t)\,\mathrm {d} B_{t}{\mbox{ for }}t\in [0,T];}
X
0
=
Z
;
{\displaystyle X_{0}=Z;}
具有P-几乎必然 独特t 连续解(t , ω ) ↦ X t (ω ),使X 适应 由Z 、B s (s ≤ t )生成的滤子F t Z ;另外
E
[
∫
0
T
|
X
t
|
2
d
t
]
<
+
∞
.
{\displaystyle \mathbb {E} \left[\int _{0}^{T}|X_{t}|^{2}\,\mathrm {d} t\right]<+\infty .}
一般情形:局部利普希茨条件与最大解
上述SDE只是更一般形式的特例:
d
Y
t
=
α
(
t
,
Y
t
)
d
X
t
{\displaystyle \mathrm {d} Y_{t}=\alpha (t,Y_{t})\mathrm {d} X_{t}}
其中
X
{\displaystyle X}
是
R
n
{\displaystyle \mathbb {R} ^{n}}
上的连续半鞅;
Y
{\displaystyle Y}
是
R
d
{\displaystyle \mathbb {R} ^{d}}
上的连续半鞅
α
:
R
+
×
U
→
Lin
(
R
n
;
R
d
)
{\displaystyle \alpha :\mathbb {R} _{+}\times U\to \operatorname {Lin} (\mathbb {R} ^{n};\mathbb {R} ^{d})}
是从某非空开集
U
⊂
R
d
{\displaystyle U\subset \mathbb {R} ^{d}}
发出的映射,其中
Lin
(
R
n
;
R
d
)
{\displaystyle \operatorname {Lin} (\mathbb {R} ^{n};\mathbb {R} ^{d})}
R
n
{\displaystyle \mathbb {R} ^{n}}
到
R
d
{\displaystyle \mathbb {R} ^{d}}
的所有线性映射的空间
更一般地说,还可以研究流形 上的SDE。
方程的解收不收敛取决于
α
{\displaystyle \alpha }
的选择。假设
α
{\displaystyle \alpha }
满足某局部利普希茨条件,即对
t
≥
0
{\displaystyle t\geq 0}
和紧集
K
⊂
U
{\displaystyle K\subset U}
、常数
L
(
t
,
K
)
{\displaystyle L(t,K)}
,满足条件
|
α
(
s
,
y
)
−
α
(
s
,
x
)
|
≤
L
(
t
,
K
)
|
y
−
x
|
,
x
,
y
∈
K
,
0
≤
s
≤
t
,
{\displaystyle |\alpha (s,y)-\alpha (s,x)|\leq L(t,K)|y-x|,\quad x,y\in K,\;0\leq s\leq t,}
其中
|
⋅
|
{\displaystyle |\cdot |}
是欧氏范数。这一条件保证了所谓最大解的存在性和唯一性。
设
α
{\displaystyle \alpha }
连续、满足上述局部利普希茨条件,并设
F
:
Ω
→
U
{\displaystyle F:\Omega \to U}
为初始条件,即是关于初始σ-代数的可测函数。令
ζ
:
Ω
→
R
¯
+
{\displaystyle \zeta :\Omega \to {\overline {\mathbb {R} }}_{+}}
为可预测停时 ,
ζ
>
0
{\displaystyle \zeta >0}
几乎确定。
U
{\displaystyle U}
值半鞅
(
Y
t
)
t
<
ζ
{\displaystyle (Y_{t})_{t<\zeta }}
即下式方程的最大解
d
Y
t
=
α
(
t
,
Y
t
)
d
X
t
,
Y
0
=
F
{\displaystyle dY_{t}=\alpha (t,Y_{t})dX_{t},\quad Y_{0}=F}
若:
对
ζ
n
↗
ζ
{\displaystyle \zeta _{n}\nearrow \zeta }
,停止过程
Y
ζ
n
{\displaystyle Y^{\zeta _{n}}}
是下式停SDE的解:
d
Y
=
α
(
t
,
Y
)
d
X
ζ
n
{\displaystyle \mathrm {d} Y=\alpha (t,Y)\mathrm {d} X^{\zeta _{n}}}
在集合
{
ζ
<
∞
}
{\displaystyle \{\zeta <\infty \}}
上,几乎可以确定
Y
t
→
∂
U
(
t
→
ζ
)
{\displaystyle Y_{t}\to \partial U\ (t\to \zeta )}
[ 23]
则
ζ
{\displaystyle \zeta }
称为寿命或所谓“爆炸时间”。
線性SDE的一般解
d
X
t
=
(
a
(
t
)
X
t
+
c
(
t
)
)
d
t
+
(
b
(
t
)
X
t
+
d
(
t
)
)
d
W
t
{\displaystyle \mathrm {d} X_{t}=(a(t)X_{t}+c(t))\mathrm {d} t+(b(t)X_{t}+d(t))\mathrm {d} W_{t}}
X
t
=
Φ
t
,
t
0
(
X
t
0
+
∫
t
0
t
Φ
s
,
t
0
−
1
(
c
(
s
)
−
b
(
s
)
d
(
s
)
)
d
s
+
∫
t
0
t
Φ
s
,
t
0
−
1
d
(
s
)
d
W
s
)
{\displaystyle X_{t}=\Phi _{t,t_{0}}\left(X_{t_{0}}+\int _{t_{0}}^{t}\Phi _{s,t_{0}}^{-1}(c(s)-b(s)\mathrm {d} (s))\mathrm {d} s+\int _{t_{0}}^{t}\Phi _{s,t_{0}}^{-1}\mathrm {d} (s)\mathrm {d} W_{s}\right)}
其中
Φ
t
,
t
0
=
exp
(
∫
t
0
t
(
a
(
s
)
−
b
2
(
s
)
2
)
d
s
+
∫
t
0
t
b
(
s
)
d
W
s
)
{\displaystyle \Phi _{t,t_{0}}=\exp \left(\int _{t_{0}}^{t}\left(a(s)-{\frac {b^{2}(s)}{2}}\right)\mathrm {d} s+\int _{t_{0}}^{t}b(s)\mathrm {d} W_{s}\right)}
可約SDE: Case 1
d
X
t
=
1
2
f
(
X
t
)
f
′
(
X
t
)
d
t
+
f
(
X
t
)
d
W
t
{\displaystyle \mathrm {d} X_{t}={\frac {1}{2}}f(X_{t})f'(X_{t})\mathrm {d} t+f(X_{t})\mathrm {d} W_{t}}
對於一個可微函數
f
{\displaystyle f}
,
d
X
t
{\displaystyle \mathrm {d} X_{t}}
等價於 斯特拉托諾維奇 SDE
d
X
t
=
f
(
X
t
)
∘
W
t
{\displaystyle \mathrm {d} X_{t}=f(X_{t})\circ W_{t}}
其有一般解
X
t
=
h
−
1
(
W
t
+
h
(
X
0
)
)
{\displaystyle X_{t}=h^{-1}(W_{t}+h(X_{0}))}
其中的
h
(
x
)
{\displaystyle h(x)}
如下所示:
h
(
x
)
=
∫
x
d
s
f
(
s
)
{\displaystyle h(x)=\int ^{x}{\frac {\mathrm {d} s}{f(s)}}}
可約SDE: Case 2
d
X
t
=
(
α
f
(
X
t
)
+
1
2
f
(
X
t
)
f
′
(
X
t
)
)
d
t
+
f
(
X
t
)
d
W
t
{\displaystyle \mathrm {d} X_{t}=\left(\alpha f(X_{t})+{\frac {1}{2}}f(X_{t})f'(X_{t})\right)\mathrm {d} t+f(X_{t})\mathrm {d} W_{t}}
對於一個可微函數
f
{\displaystyle f}
,
d
X
t
{\displaystyle \mathrm {d} X_{t}}
等價於 斯特拉托諾維奇 SDE
d
X
t
=
α
f
(
X
t
)
d
t
+
f
(
X
t
)
∘
W
t
{\displaystyle \mathrm {d} X_{t}=\alpha f(X_{t})\mathrm {d} t+f(X_{t})\circ W_{t}}
可以簡化成下列形式:
d
Y
t
=
α
d
t
+
d
W
t
{\displaystyle \mathrm {d} Y_{t}=\alpha \mathrm {d} t+\mathrm {d} W_{t}}
其中
Y
t
=
h
(
X
t
)
{\displaystyle Y_{t}=h(X_{t})}
,這裡的
h
{\displaystyle h}
如上面所定義。他的一般解可以寫成
X
t
=
h
−
1
(
α
t
+
W
t
+
h
(
X
0
)
)
{\displaystyle X_{t}=h^{-1}(\alpha t+W_{t}+h(X_{0}))}
SDE与超对称
在SDE的超对称理论中,随机动力是通过作用于模型相空间微分形式 的随机演化算子定义的。在这一精确表述中,所有SDE都具有拓扑超对称性 ,即通过连续的时间流保持相空间连续性。这种超对称的自发破缺是混沌 、湍流 、自组织临界性 等诸多动力现象的数学本质,而南部定理 则解释了相关的长距动力行为,如蝴蝶效应 、粉红噪声 、爆裂声 ,以及地震、神经震荡、太阳耀斑等现象的无标统计等等。
相关条目
参考文献
^ 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 .
^ 陳, 松男. 金融數學與隨機微積分. ISBN 9789574145089 .
^ 3.0 3.1 3.2 Musiela, M., and Rutkowski, M. (2004), Martingale Methods in Financial Modelling, 2nd Edition, Springer Verlag, Berlin.
^ 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 .
^ 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
^ 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.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
^ 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.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.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 .
^ 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.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
^ 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 (英语) .
^ Fengler, M. R. (2005), Semiparametric modeling of implied volatility, Springer Verlag, Berlin. DOI https://doi.org/10.1007/3-540-30591-2
^ 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
^ Detecting Market Abuse . Risk Magazine. [2023-10-15 ] . (原始内容存档 于2023-05-18).
^ 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).
^ 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.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
^ 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).