第 5 章 随机微分方程(Stochastic Differential Equations)
作者
Bernt Øksendal(奥斯陆大学)。本章是 SDE 理论的核心:把 Ch 1 提出的 7 个 Problem 转化为严格的 SDE 问题,给出存在唯一性定理,并把"如何显式求解线性 SDE"作为例证。
Itô Kiyosi 1951 论文已经给出 Lipschitz 条件下的存在唯一性(Picard 迭代 + Gronwall),本章沿用这一框架,但同时强调工程应用:用 4 个例子(种群增长 / LRC 电路 / 圆上 BM / Ornstein-Uhlenbeck)展示 SDE 的物理 / 工程含义。"强解 vs 弱解" 的区分来自 Stroock-Varadhan 1979,是 SDE 理论的精细化——本章用 Tanaka 方程(Ex 5.3.2)作为"无强解但有弱解"的经典反例。
内容概述
本章正式"回答" Ch 1 提出的 7 个 Problem 1–2:含噪种群增长 + 含噪 LRC 电路。Ch 6+ 将继续回答 Problem 3–7。
本章中心任务:给定 SDE \(dX_t = b(t, X_t)\,dt + \sigma(t, X_t)\,dB_t\),(A) 何时存在唯一解?(B) 如何显式求解?
主要内容: 1. 4 个例子(Ch 5.1): - Example 5.1.1:几何布朗运动 \(dN_t = rN_t dt + \alpha N_t dB_t\) → \(N_t = N_0 \exp((r - \tfrac12 \alpha^2)t + \alpha B_t)\)。这是 Ch 1 Problem 1 的解。 - Example 5.1.3:含噪 LRC 电路 → 2 维线性 SDE → 矩阵指数 \(e^{At}\)。Ch 1 Problem 2 的解。 - Example 5.1.4:\(B\) 在单位圆上 → \(dY = -\tfrac12 Y dt + KY dB\)。 - 重述 Itô vs Stratonovich:同一 \(\sigma\) 系数 \(\alpha\) 下 Itô 解和 Stratonovich 解差 \(\tfrac12 \alpha^2 t\) 漂移修正——两种解的渐近行为不同(LIL 5.1.2:Itô 解 \(N_t \to 0\) a.s. 当 \(r < \tfrac12 \alpha^2\);Stratonovich 解 \(N_t \to 0\) a.s. 当 \(r < 0\))。 2. 存在唯一性定理(Thm 5.2.1):Lipschitz + 线性增长 → SDE 有唯一 \(t\)-连续解。证明用 Picard 迭代 + Gronwall 不等式。 3. 强解 vs 弱解(Ch 5.3): - 强解:\(B_t\) 固定,\(X_t\) 是 \(\mathcal F^Z_t\)-适应的。 - 弱解:\(B_t\) 和 \(X_t\) 同时构造。 - 反例(Ex 5.3.2 Tanaka 方程):\(dX_t = \text{sign}(X_t)\,dB_t\), \(X_0 = 0\),\(\sigma = \text{sign}\) 不满足 Lipschitz,故无强解;但有弱解(\(X_t\) 取任何 BM,再定义 \(B_t\))。 4. Tanaka 方程的意义:\(\sigma\) 仅在 \(0\) 一点不连续 / 不可微,强解可不存在;这意味着"\(\sigma\) 整体 Lipschitz"是"强解存在"的关键。 5. 弱唯一性(Lem 5.3.1):两个弱解的"分布"相同。 6. Ornstein-Uhlenbeck 过程(Ex 5.5):\(dX_t = \mu X_t dt + \sigma dB_t\) → 显式解 \(X_t = e^{\mu t} X_0 + \sigma \int_0^t e^{\mu(t-s)} dB_s\)。均值回归、平稳分布等性质是 Merton 投资组合(Ch 11)和利率模型(Ch 12 扩展)的标准工具。
在全书中位置:核心方法论章。Ch 6 滤波(线性 SDE)、Ch 7 扩散的 Markov 性质、Ch 8 Feynman-Kac、Ch 10 停时最优停、Ch 11 HJB 方程、Ch 12 Black-Scholes 都依赖本章的存在唯一性 + Itô 公式。
前置知识:Ch 3(Itô 积分)+ Ch 4(Itô 公式、鞅表示、二次变差)。Picard 迭代的证明仅需 Gronwall 不等式(标准 ODE 工具),但 Itô 等距是 SDE 情形比 ODE 情形多的关键。
核心方程与概念
本章是 SDE 显式求解 + 存在唯一性理论的"工作章节"。下面列出最重要的 9 个对象。
5.1 几何布朗运动(Geometric BM, Example 5.1.1)— Ch 1 Problem 1 的解
- Itô 公式推导:取 \(g(t, x) = \log x\),\(d(\log N_t) = (1/N_t)\,dN_t - \tfrac12 (1/N_t^2)(dN_t)^2 = dN_t/N_t - \tfrac12 \alpha^2\,dt\)。两边积分得 \(\log(N_t/N_0) = rt + \alpha B_t - \tfrac12 \alpha^2 t\)。
- Itô 修正项 \(-\tfrac12 \alpha^2 t\) 的来源:\(\tfrac12\) 来自 Itô 公式的二阶项,\(\alpha^2\) 来自 \(\sigma = \alpha N_t\)。
- 渐近行为(LIL Thm 5.1.2):\(\limsup_{t\to\infty} B_t / \sqrt{2t \log\log t} = 1\) a.s.,故
- \(r > \tfrac12 \alpha^2\):\(N_t \to \infty\) a.s.(指数增长胜出对数扰动)
- \(r < \tfrac12 \alpha^2\):\(N_t \to 0\) a.s.
- \(r = \tfrac12 \alpha^2\):\(N_t\) 在 0 和 \(\infty\) 之间反复波动
- 与 Stratonovich 解的对比:Stratonovich \(N_t^{(S)} = N_0 e^{rt + \alpha B_t}\),故 \(E[N_t^{(S)}] = N_0 e^{(r + \tfrac12 \alpha^2)t}\)(多了 \(\tfrac12 \alpha^2 t\) 漂移)。Itô 解和 Stratonovich 解对相同的 SDE 形式给出不同的过程。
- Ch 12 应用:几何 BM 是 Black-Scholes 股票价格模型的标准选择。
5.2 LRC 电路(Example 5.1.3)— Ch 1 Problem 2 的解
定义 \(X_t = (Q_t, Q'_t)^T \in \mathbb R^2\):
求解:乘以 \(e^{-At}\)(矩阵指数),用多维 Itô 公式:
- 关键认识:线性 SDE 可用"矩阵指数 \(e^{At}\)"作 integrating factor,把 SDE 化为可直接积分的形式。这与 ODE 情形 \(dX/dt = AX\) 的解 \(X = e^{At}X_0\) 完全平行。
- Ch 6 应用:Kalman-Bucy 滤波的"信号过程" \(dX = AX\,dt + K\,dB\) 即线性 SDE。
5.3 圆上 Brown Motion(Example 5.1.4)
取 \(g(t, x) = e^{ix} = (\cos x, \sin x)\),\(X_t = B_t\)。Itô 公式给出
故 \(Y = (Y_1, Y_2) = e^{iB_t}\) 满足
- 几何解释:\(Y_t\) 在单位圆 \(S^1\) 上走,\(K\) 是切向方向(\(K^2 = -I\),故 \(K\) 是 \(90°\) 旋转)。
- Ch 11 应用:单位圆上的几何 BM 是 Heston 随机波动率模型(Ch 12 扩展)的简化。
5.4 重对数律(Thm 5.1.2)
- 意义:BM 路径的精确增长率——比 \(\sqrt{t}\) 慢(系数 \(1/\sqrt{\log\log t}\) 衰减),但比任何 \(t^\alpha\)(\(\alpha < 1/2\))快。
- Ch 7 应用:用 LIL 证明扩散过程的"常返性 / 逃逸性"——对 \(n = 2\) 维 BM,\(|B_t| \to \infty\) 实际上比 \(\sqrt{t}\) 慢(\(B_t/\sqrt{t \log\log t} \to 0\) a.s.),故 \(B_t\) 在平面上是常返的(Ch 9 Dirichlet 边值的可解性条件)。
- 历史:Khintchine 1924 首次证明 LIL;Lévy 1937、Strassen 1964 给出强化版本。
5.5 Itô vs Stratonovich 的渐近差(Ch 5.1 末)
对种群增长模型 \(dN_t = rN_t dt + \alpha N_t dW_t\):
| Itô 解 | Stratonovich 解 | |
|---|---|---|
| 表达式 | \(N_t = N_0 e^{(r - \alpha^2/2)t + \alpha B_t}\) | \(N_t^{(S)} = N_0 e^{rt + \alpha B_t}\) |
| \(E[N_t]\) | \(N_0 e^{rt}\) | \(N_0 e^{(r + \alpha^2/2)t}\) |
| 渐近 | \(r > \alpha^2/2 \Rightarrow N_t \to \infty\) a.s. | \(r > 0 \Rightarrow N_t^{(S)} \to \infty\) a.s. |
- 解读:Stratonovich 解在均值上多 \(\tfrac12 \alpha^2 t\) 的"漂移"——这是 Itô 公式 \(\tfrac12 \sigma^2\) 项的 Stratonovich 版本(无修正)。
- 建模选择:若物理系统的扰动是"真实随机微分"(如电磁噪声),Stratonovich 极限更自然(Wong-Zakai 1969);若扰动是"对数价格波动"(金融 tick data),Itô 假设下无套利(Ch 12)。
5.6 存在唯一性定理(Thm 5.2.1)
条件:\(b, \sigma\) 满足 - (5.2.1) 线性增长:\(|b(t, x)| + |\sigma(t, x)| \le C(1 + |x|)\) - (5.2.2) Lipschitz:\(|b(t, x) - b(t, y)| + |\sigma(t, x) - \sigma(t, y)| \le D|x - y|\)
结论:SDE (5.2.3) 有唯一 \(t\)-连续解 \(X_t\),适应于 \(\mathcal F^Z_t\)(\(Z\) = 初始条件),且 \(E[\int_0^T |X_t|^2 dt] < \infty\)。
唯一性证明(Picard 距离估计):设 \(X^1, X^2\) 是两个解, $\(E[|X^1_t - X^2_t|^2] \le 3E[|Z - \tilde Z|^2] + 3(1 + T) D^2 \int_0^t E[|X^1_s - X^2_s|^2] ds.\)$ 用 Gronwall 不等式 \(v(t) \le F + A\int_0^t v(s)\,ds \Rightarrow v(t) \le F e^{At}\)。\(Z = \tilde Z\) 时 \(F = 0\),\(v(t) = 0\)。
存在性证明(Picard 迭代): $\(Y^{(k+1)}_t = Z + \int_0^t b(s, Y^{(k)}_s)\,ds + \int_0^t \sigma(s, Y^{(k)}_s)\,dB_s, \quad Y^{(0)}_t = Z.\)$ 归纳证 \(E[|Y^{(k+1)}_t - Y^{(k)}_t|^2] \le A_2 t^{k+1} / (k+1)!\)(关键因子 \((k+1)!\))。\(L^2\) Cauchy 序列 → 极限 \(X_t\) 满足 SDE。
- Gronwall 不等式:ODE 中 \(v(t) \le \int_0^t \beta(s) v(s)\,ds \Rightarrow v(t) = 0\)。本章用"差 $v(t) - $ 常数估计"形式。
5.7 强解 vs 弱解(Ch 5.3)
| 强解 | 弱解 | |
|---|---|---|
| \(B_t\) 何时给 | 先给(背景 \(\sigma\)-代数固定) | 同时构造(与 \(X_t\) 一起) |
| \(X_t\) 适应于 | \(\mathcal F^Z_t\)(含 \(B\) 历史) | 某个 \(\mathcal H_t\)(由 \((X, B)\) 共同生成) |
| 存在性 | 需要 \(\sigma\) Lipschitz | 仅需 \(\sigma\) 连续性 + 弱唯一性 |
| Tanaka 方程 \(dX = \text{sign}(X)\,dB\) | 无 | 有 |
- 强解的"路径唯一性":\(X_t^1(\omega) = X_t^2(\omega)\) a.s. for all \(t\)。
- 弱解的"分布唯一性":\(X_t^1\) 和 \(X_t^2\) 有相同有限维分布。
- 强解存在 \(\Rightarrow\) 弱解存在 + 弱唯一。反之不真(Tanaka 反例)。
- Ch 12 应用:金融工程用"弱解"——因为市场无套利 \(\Rightarrow\) 折现价格鞅存在(等价鞅测度),但不指定 BM 路径。
5.8 Tanaka 方程(Ex 5.3.2)—— 弱解但无强解的经典反例
- 无强解的论证:假设 \(X_t\) 是强解。由 Tanaka 公式 (4.3.12),\(B_t = \int_0^t \text{sign}(X_s)\,dX_s = |X_t| - L_t\),故 \(B_t\) 由 \((|X_s|: s \le t)\) 决定——\(\mathcal F_t \subset \mathcal M_t^{|X|}\)。但 \(B_t\) 也是 Itô 积分 \(\int_0^t \text{sign}(X_s)\,dX_s\)——反之 \(\mathcal M_t^{|X|} \subset \mathcal F_t\)。矛盾。
- 有弱解:取 \(X_t = \tilde B_t\)(任何 BM),定义 \(B_t = \int_0^t \text{sign}(\tilde B_s)\,d\tilde B_s\)——直接满足 \(dX_t = \text{sign}(X_t)\,dB_t\)。
- 意义:在 \(\sigma\) 仅在 0 不连续(\(\text{sign}\) 在 0 跳变)的情形,强解不存在;弱解存在且唯一。这警告"\(\sigma\) 整体 Lipschitz"不是 SDE 可解性的必要条件——只是强解存在的充分条件。
5.9 Ornstein-Uhlenbeck 过程(Ex 5.5)—— 平稳分布的范例
显式解(用 integrating factor \(e^{-\mu t}\)):
- 均值:\(E[X_t] = e^{\mu t} X_0\)(\(\mu > 0\) 指数发散,\(\mu < 0\) 指数回归 0)。
- 方差:\(\text{Var}(X_t) = \sigma^2 \int_0^t e^{2\mu(t-s)}\,ds = \sigma^2 (e^{2\mu t} - 1)/(2\mu)\)。\(\mu < 0\) 时收敛到 \(-\sigma^2/(2\mu)\),平稳分布 \(N(0, -\sigma^2/(2\mu))\)。
- OU 过程 vs BM:BM 无平稳分布,OU 有(高斯)。这是"均值回归"动力学的标准模型——Ch 11 利率模型、Ch 12 Heston 波动率模型都用 OU 类结构。
关键结论
- Itô 公式是显式求解 SDE 的核心工具:4 个 Example(5.1.1, 5.1.3, 5.1.4, Ex 5.5)都用 Itô 公式(Ch 4)作 integrating factor 或对数变换。Itô 公式 + 适当 \(g\) 变换 = 80% 的 SDE 显式求解。
- 存在唯一性只需 Lipschitz + 线性增长(Thm 5.2.1)——与 ODE 情形完全平行(Picard-Lindelöf 定理)。但 SDE 情形下证明更复杂:要用 Itô 等距(Ch 3 核心定理)估计 Picard 迭代的 \(L^2\) 收敛。
- Gronwall 不等式是 SDE 估计的"瑞士军刀":唯一性证明、存在性证明、解的稳定性估计、爆炸性反例(Remark a/b)都依赖 Gronwall。这是 SDE 与 ODE 共用工具的关键例子。
- 强解 vs 弱解的区分在金融工程中至关重要:
- 强解假设:定价/对冲时"完美知道" BM 路径 → 唯一对冲组合(Ch 12 完备市场)。
- 弱解假设:市场无套利 → 存在等价鞅测度 → 唯一风险中性定价。
- 两者关系:强解存在 ⇒ 弱解存在 + 弱唯一;反之不真(Tanaka)。
- Itô vs Stratonovich 的渐近差(Ch 5.1 末表)——同一 \(\alpha\) 系数下 Itô 解"均值"= \(N_0 e^{rt}\)(无 \(\alpha^2\) 项),Stratonovich 解"均值"= \(N_0 e^{(r + \alpha^2/2) t}\)。这种 \(\tfrac12 \alpha^2 t\) 差异是建模选择的物理体现,不是数学任意。
- LIL (Thm 5.1.2) 是 BM 渐近分析的"显微镜":\(\limsup B_t / \sqrt{2t \log\log t} = 1\) a.s. 给出 BM 路径的精确增长率,比"\(\sqrt{t}\)"慢(系数 \(1/\sqrt{\log\log t}\) 衰减),但比 \(t^\alpha\)(\(\alpha < 1/2\))快。Ch 7 / Ch 9 多次用此。
- 线性 SDE 都有"矩阵指数 \(e^{At}\)"解(Ex 5.1.3, 5.5):\(X(t) = e^{At} X_0 + \int_0^t e^{A(t-s)} (\sigma\,dB_s + b\,ds)\)。这与 ODE 情形完全平行——线性 SDE = ODE 的自然推广。
挑战和开放性问题
- 存在唯一性的"必要条件"问题:Thm 5.2.1 的条件(Lipschitz + 线性增长)充分但不必要。例如:
- \(\sigma\) 局部 Lipschitz 但全局增长(\(dX = X^2 dB\))→ 存在局部解但可能爆炸到 \(\infty\) in 有限时间。
- \(\sigma\) 单边 Lipschitz(\(dX = \text{sign}(X) dB\))→ 弱解存在(Tanaka),强解不存在。
- 退化 SDE(\(\sigma = 0\) in some region)→ 多解或不可微。
- 局部时间依赖系数 \(b(t, x), \sigma(t, x)\) 的非一致 Lipschitz 情形:路径依赖 SDE(McKean-Vlasov 方程)。
- 高维 SDE 的"维数耦合":Thm 5.2.1 在 \(\mathbb R^n\) 中给出存在唯一性,但 \(\sigma \in \mathbb R^{n \times m}\) 的耦合效应使多维 SDE 的几何性质(如解的支撑集、长时间行为、平稳分布)远比线性情形复杂。
- Stratonovich vs Itô 的"建模"选择——非数学问题:
- 工程逼近(Wong-Zakai 极限)→ Stratonovich
- 金融无套利 → Itô
- 物理热噪声 / 电磁噪声 → 实际不严格归类,文献中"两个都有人用"
- 这意味着"用 SDE 解工程问题"时必须明确写出"用哪个积分",否则结果不同。
- Tanaka 方程的"病理"边界:\(\sigma = \text{sign}\) 仅在 0 一点不连续。如果 \(\sigma\) 在更复杂的边界(\(|\sigma|\) 平方可积但不连续),弱解是否存在?Stroock-Varadhan 1979 给出"广义 SDE" 的存在性定理,要求 \(\sigma\) 一致椭圆 + 局部 Hölder;Zvonkin-Krylov 1981 用 Itô 公式 + 强解存在性的精巧论证处理 \(\sigma\) 退化情形。
- Picard 迭代的"收敛速度":Thm 5.2.1 证明用 \((k+1)!\) 因子的 \(L^2\) 估计,但实际数值算法(Euler-Maruyama、Milstein)收敛速度取决于 \(b, \sigma\) 的光滑性。理论存在性 \(\ne\) 数值高效——Ch 12 定价算法的精度问题(隐式 Euler, Runge-Kutta SDE 等)都是开放方向。
- SDE 解的"长时间行为":本章只给出 \([0, T]\) 上的局部存在唯一性。平稳分布、遍历性、混合时间等问题需要 ergodic theory 工具(Birkhoff, Oseledets 乘法遍历定理),Ch 7 / Ch 8 简述。
- SDE 解的"支撑集"问题:在 \(dX = b\,dt + \sigma\,dB\) 中,\(\sigma\) 的支集决定 \(X\) 能否到达某些区域。Stroock-Varadhan 支持定理(1969)给出"\(\sigma\) 的秩"控制解的可达性——是 Ch 9 边界可达性、Ch 12 障碍期权可达性的关键工具。
- 延迟 SDE / 路径依赖 SDE:\(dX_t = b(t, X_t, X_{t-\tau})\,dt + \sigma(t, X_t, X_{t-\tau})\,dB_t\)(\(\tau > 0\) 延迟)——在金融(利率期限结构、volatility surface 拟合)、生物(记忆效应)有重要应用,但 Thm 5.2.1 不直接适用。需要新理论(Mohammed 1984, Reiß 2002)。
个人反思与批判性分析
本章是 Itô 理论"对 ODE"的胜利——Picard 迭代 + Gronwall + Lipschitz 的标准 ODE 工具 + Itô 等距的 SDE 修正 = SDE 的存在唯一性。这是 SDE 在 1960-70 年代成为成熟数学分支的标志——之前 SDE 常被视为"带噪声的 ODE",无独立理论;之后 SDE 有自己的"基本定理三元组"(Itô 公式 + 鞅表示 + SDE 存在唯一性)。
强解 vs 弱解的"哲学"对立: - 强解视角:"物理噪声"是固定的(Wong-Zakai 工程逼近),SDE 解是唯一确定的。 - 弱解视角:"物理噪声"是"等价类"(无套利的金融建模),SDE 解只需分布正确。 - 两种视角在不同应用中各有优势。强解更"经典"——给出"路径";弱解更"灵活"——不预设路径唯一性。Tanaka 方程表明弱解可以存在而强解不存在,提示"强解假设"在某些 SDE 中过强。
Itô vs Stratonovich 的工程含义——同一 \(\alpha\) 系数下 Itô 解和 Stratonovich 解的渐近行为不同(LIL 5.1.2 应用给出临界 \(r\) 不同:\(r < \alpha^2/2\) vs \(r < 0\))。这意味着"用哪个积分"不是数学细节,是建模决策——直接影响预测。物理学家一般用 Stratonovich(Wong-Zakai 极限),金融工程师一般用 Itô(无套利)。两种选择背后的"哲学"分歧: - Itô 强调"信息流":过去 → 现在的因果链。 - Stratonovich 强调"测度论":Wiener 过程的连续逼近(用光滑过程逼近 BM)。
对 ODE 经验的迁移价值——存在唯一性定理的形式与 ODE 完全平行(Lipschitz + 线性增长 → Picard 迭代 + Gronwall)。这给 ODE 背景的研究者一个"快速入口":不需要重新学测度论/鞅论,只需"加上 Itô 等距估计"即可。这是 Itô 理论的"工程师友好"特性。
与 Karatzas-Shreve (1998) 的比较——KS Ch 5 给出 SDE 存在唯一性,但条件较弱(局部 Lipschitz + 局部增长)。Øksendal 用全局 Lipschitz + 线性增长简化条件——对初学者友好,但丢失了"解的爆炸性"分析(局部 Lipschitz 时解可能在有限时间爆炸到 \(\infty\))。
对初学者的建议——本章的 Exercise 5.1(5 个具体 SDE 的解)、Exercise 5.3(多维线性 SDE)、Exercise 5.5(OU 过程)都是 Ch 6 滤波 / Ch 12 金融工程的"必备工具"。亲手做完后,Ch 6 几乎不需要重读 Picard 迭代证明。
重要参考文献
[X1] K. Itô. On stochastic differential equations. Memoirs of the American Mathematical Society 4: 1–51, 1951. — Itô SDE 存在唯一性定理的首次严格证明(Picard + Gronwall)。本章 Thm 5.2.1 的来源。
[X2] A. Friedman. Stochastic Differential Equations and Applications, Vol. 1. Academic Press, 1975. ISBN 978-0122684012. — SDE 存在唯一性的不同条件(局部 Lipschitz、爆破时间)研究(Ch 4)。
[X3] D. W. Stroock, S. R. S. Varadhan. Multidimensional Diffusion Processes. Springer, 1979. ISBN 978-3540902983. — 弱解存在性 + 弱唯一性(Stroock-Varadhan 支持定理)。
[X4] N. V. Krylov, A. K. Zvonkin. On a certain class of stochastic equations. Theory of Probability and Its Applications 25: 790–793, 1981. — 强解存在性对 \(\sigma\) 退化的处理(Krylov-Zvonkin 论证)。
[X5] A. Khintchine. Über dyadische Reihen. Mathematische Zeitschrift 24: 309–325, 1926. — 重对数律 LIL 的首次证明。
[X6] E. Wong, M. Zakai. On the relation between ordinary and stochastic differential equations. International Journal of Engineering Science 3: 213–229, 1965. — Wong-Zakai 逼近:光滑过程逼近 BM 时 ODE 极限 → Stratonovich SDE("建模选 Itô 还是 Stratonovich"的物理基础)。
[X7] J. Lamperti. Stochastic Processes: A Survey of the Mathematical Theory. Springer, 1977. ISBN 978-0387902753. — LIL 证明的标准参考(§22, Thm 5.1.2 引用)。
[X8] T. C. Gard. Introduction to Stochastic Differential Equations. Marcel Dekker, 1988. ISBN 978-0824776974. — 1 维 SDE 化简方法(Ch 4, 引用在 Ch 5.1 末)。
[X9] G. E. Uhlenbeck, L. S. Ornstein. On the theory of the Brownian motion. Physical Review 36: 823–841, 1930. — OU 过程首次引入(Ex 5.5 引用)。
[X10] P. Langevin. Sur la théorie du mouvement brownien. Comptes Rendus de l'Académie des Sciences (Paris) 146: 530–533, 1908. — Langevin 方程首次引入(OU 过程的物理原型)。
[X11] D. Revuz, M. Yor. Continuous Martingales and Brownian Motion, 3rd ed. Springer, 1999. ISBN 978-3540643253. — OU 过程平稳分布、遍历性的标准参考。
[X12] I. Karatzas, S. E. Shreve. Brownian Motion and Stochastic Calculus, 2nd ed. Springer, 1991. ISBN 978-0387976556. — 强解 vs 弱解的完整理论(KS Ch 5)、Picard 迭代的更精细估计(KS Ch 5.2D)。