跳转至

第 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-KacCh 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 的解

\[dN_t = rN_t\,dt + \alpha N_t\,dB_t, \quad N_0 = N_0.\]
\[N_t = N_0 \exp\!\left(\left(r - \frac{1}{2}\alpha^2\right) t + \alpha B_t\right).\]
  • 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 的解

\[L Q'' + R Q' + \tfrac{1}{C} Q = G_t + \alpha W_t, \quad Q(0) = Q_0, \; Q'(0) = I_0.\]

定义 \(X_t = (Q_t, Q'_t)^T \in \mathbb R^2\)

\[dX = AX\,dt + H(t)\,dt + K\,dB_t, \quad A = \begin{pmatrix} 0 & 1 \\ -1/C & -R/L \end{pmatrix}, \; H(t) = \begin{pmatrix} 0 \\ G_t/L \end{pmatrix}, \; K = \begin{pmatrix} 0 \\ \alpha/L \end{pmatrix}.\]

求解:乘以 \(e^{-At}\)(矩阵指数),用多维 Itô 公式:

\[d(e^{-At} X(t)) = e^{-At} H(t)\,dt + e^{-At} K\,dB_t,\]
\[X(t) = e^{At} X(0) + \int_0^t e^{A(t-s)} H(s)\,ds + \int_0^t e^{A(t-s)} K\,dB_s.\]
  • 关键认识:线性 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ô 公式给出

\[dY_1 = -\sin(B_t)\,dB_t - \tfrac{1}{2}\cos(B_t)\,dt, \quad dY_2 = \cos(B_t)\,dB_t - \tfrac{1}{2}\sin(B_t)\,dt.\]

\(Y = (Y_1, Y_2) = e^{iB_t}\) 满足

\[dY = -\tfrac{1}{2} Y\,dt + K Y\,dB, \quad K = \begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}.\]
  • 几何解释\(Y_t\) 在单位圆 \(S^1\) 上走,\(K\) 是切向方向(\(K^2 = -I\),故 \(K\)\(90°\) 旋转)。
  • Ch 11 应用:单位圆上的几何 BM 是 Heston 随机波动率模型(Ch 12 扩展)的简化。

5.4 重对数律(Thm 5.1.2)

\[\limsup_{t \to \infty} \frac{B_t}{\sqrt{2t \log\log t}} = 1 \quad \text{a.s.}\]
  • 意义: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)—— 弱解但无强解的经典反例

\[dX_t = \text{sign}(X_t)\,dB_t, \quad X_0 = 0, \quad \text{sign}(x) = \begin{cases} +1 & x \ge 0 \\ -1 & x < 0 \end{cases}.\]
  • 无强解的论证:假设 \(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)—— 平稳分布的范例

\[dX_t = \mu X_t\,dt + \sigma\,dB_t, \quad \mu, \sigma \in \mathbb R, \; B_t \in \mathbb R.\]

显式解(用 integrating factor \(e^{-\mu t}\)):

\[X_t = e^{\mu t} X_0 + \sigma \int_0^t e^{\mu(t-s)}\,dB_s.\]
  • 均值\(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 类结构。

关键结论

  1. 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 显式求解
  2. 存在唯一性只需 Lipschitz + 线性增长(Thm 5.2.1)——与 ODE 情形完全平行(Picard-Lindelöf 定理)。但 SDE 情形下证明更复杂:要用 Itô 等距(Ch 3 核心定理)估计 Picard 迭代的 \(L^2\) 收敛。
  3. Gronwall 不等式是 SDE 估计的"瑞士军刀":唯一性证明、存在性证明、解的稳定性估计、爆炸性反例(Remark a/b)都依赖 Gronwall。这是 SDE 与 ODE 共用工具的关键例子
  4. 强解 vs 弱解的区分在金融工程中至关重要
  5. 强解假设:定价/对冲时"完美知道" BM 路径 → 唯一对冲组合(Ch 12 完备市场)。
  6. 弱解假设:市场无套利 → 存在等价鞅测度 → 唯一风险中性定价。
  7. 两者关系:强解存在 ⇒ 弱解存在 + 弱唯一;反之不真(Tanaka)。
  8. 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\) 差异是建模选择的物理体现,不是数学任意
  9. 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 多次用此。
  10. 线性 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 的自然推广。

挑战和开放性问题

  1. 存在唯一性的"必要条件"问题:Thm 5.2.1 的条件(Lipschitz + 线性增长)充分但不必要。例如:
  2. \(\sigma\) 局部 Lipschitz 但全局增长(\(dX = X^2 dB\))→ 存在局部解但可能爆炸到 \(\infty\) in 有限时间。
  3. \(\sigma\) 单边 Lipschitz(\(dX = \text{sign}(X) dB\))→ 弱解存在(Tanaka),强解不存在。
  4. 退化 SDE(\(\sigma = 0\) in some region)→ 多解或不可微。
  5. 局部时间依赖系数 \(b(t, x), \sigma(t, x)\) 的非一致 Lipschitz 情形:路径依赖 SDE(McKean-Vlasov 方程)。
  6. 高维 SDE 的"维数耦合":Thm 5.2.1 在 \(\mathbb R^n\) 中给出存在唯一性,但 \(\sigma \in \mathbb R^{n \times m}\) 的耦合效应使多维 SDE 的几何性质(如解的支撑集、长时间行为、平稳分布)远比线性情形复杂。
  7. Stratonovich vs Itô 的"建模"选择——非数学问题:
  8. 工程逼近(Wong-Zakai 极限)→ Stratonovich
  9. 金融无套利 → Itô
  10. 物理热噪声 / 电磁噪声 → 实际不严格归类,文献中"两个都有人用"
  11. 这意味着"用 SDE 解工程问题"时必须明确写出"用哪个积分",否则结果不同。
  12. Tanaka 方程的"病理"边界\(\sigma = \text{sign}\) 仅在 0 一点不连续。如果 \(\sigma\) 在更复杂的边界(\(|\sigma|\) 平方可积但不连续),弱解是否存在?Stroock-Varadhan 1979 给出"广义 SDE" 的存在性定理,要求 \(\sigma\) 一致椭圆 + 局部 Hölder;Zvonkin-Krylov 1981 用 Itô 公式 + 强解存在性的精巧论证处理 \(\sigma\) 退化情形。
  13. Picard 迭代的"收敛速度":Thm 5.2.1 证明用 \((k+1)!\) 因子的 \(L^2\) 估计,但实际数值算法(Euler-Maruyama、Milstein)收敛速度取决于 \(b, \sigma\) 的光滑性。理论存在性 \(\ne\) 数值高效——Ch 12 定价算法的精度问题(隐式 Euler, Runge-Kutta SDE 等)都是开放方向。
  14. SDE 解的"长时间行为":本章只给出 \([0, T]\) 上的局部存在唯一性。平稳分布、遍历性、混合时间等问题需要 ergodic theory 工具(Birkhoff, Oseledets 乘法遍历定理),Ch 7 / Ch 8 简述。
  15. SDE 解的"支撑集"问题:在 \(dX = b\,dt + \sigma\,dB\) 中,\(\sigma\) 的支集决定 \(X\) 能否到达某些区域。Stroock-Varadhan 支持定理(1969)给出"\(\sigma\) 的秩"控制解的可达性——是 Ch 9 边界可达性、Ch 12 障碍期权可达性的关键工具。
  16. 延迟 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)。