Symmetries and exact solutions of stochastic differential equations: the CIR model

技术摘要：随机微分方程的对称性与精确解：CIR 模型

问题陈述
Cox-Ingell-Ross (CIR) 模型是用于描述利率期限结构的经典随机微分方程 (SDE)。在模型的参数（$\kappa$, $\theta$, 和 $s$）之间存在一个著名的临界条件，即 $4\kappa\theta = s^2$。在此条件下，该模型在强解的存在性以及似然估计量的行为方面表现出特定的性质。虽然近期文献通过对应前向 Kolmogorov 方程的对称性视角分析了这一约束，但参数约束与 SDE 本身内部对称性的存在之间的直接关系，以及此类对称性如何促进精确解的求解，仍需进行显式论证。

方法论
本文采用了 Itô 型随机微分方程的对称性分析方法。其研究步骤如下：

1. 对称性定义： 作者定义了一个无穷小对称生成元 $V = \phi(x, t; w)\partial_x$，用于标量 Itô 方程 $dx = f(x,t)dt + \sigma(x,t)dw$。
2. 确定方程： $V$ 成为对称性的条件导致了一个关于漂移项 $f$、噪声项 $\sigma$ 和生成元函数 $\phi$ 的超定偏微分方程组（即确定方程）。这些方程利用了 Itô 拉普拉斯算子，并可以通过 Stratonovich 漂移项 $b(x,t)$ 进行重构。
3. 分类： 作者将此框架应用于具有泊松型噪声（$\sigma = s\sqrt{x}$）的一般自回归 Itô 方程，这涵盖了 CIR 模型。分析根据确定方程的相容性将其分为三种情况：
   • 情况 1：对漂移项无约束（导致平方贝塞尔过程）。
   • 情况 2：存在产生无对称性的约束。
   • 情况 3：对漂移项和噪声参数存在非平凡约束。
4. 解的构造： 当识别出对称性时，应用 Kozlov 变量代换法。该变换将原始 SDE 映射为一个新的方程，其中漂移项和噪声项与新的状态变量无关，从而实现显式积分。

主要贡献与结果

• 推导临界关系： 通过求解 CIR 模型的确定方程，本文证明了当参数满足关系式 $4\kappa\theta = s^2$ 时，存在非平凡对称性。这证实了此前在解的存在性和统计估计背景下观察到的临界条件，与该特定模型具有此类对称性的情形相一致。作者指出，该结果适用于 CIR 模型的一个特定情形，而非一般形式。
• 显式对称生成元： 对于 $4\kappa\theta = s^2$ 的特定情况，本文推导出了显式对称生成元的形式：
  $$\phi(x, t; w) = \sqrt{x} e^{-\kappa t/2}$$
• 通过对称性获得精确解： 利用 Kozlov 变换 $y(t) = \int \frac{dx}{\phi}$，作者将 CIR SDE 转化为一个更简单的方程 $dy = s e^{\kappa t/2} dw$。对该方程进行积分，可得到利率 $r(t)$（在变换变量中记为 $x(t)$）的精确强解：
  $$x(t) = \frac{1}{4} e^{-\kappa t} \left[ 2\sqrt{x_0} + s \int_0^t e^{\kappa t'/2} dw(t') \right]^2$$
  该结果重现了文献（如 [19]）中已有的已知解，但直接从方程的对称性质出发进行了推导。
• 推广： 该分析还识别了其他具有泊松噪声的漂移函数下的对称性，包括一种特定的非线性漂移情形（命题 2），尽管其主要焦点仍在于 CIR 模型的线性漂移。

意义与主张
本文声称其主要贡献在于建立了参数约束 $4\kappa\theta = s^2$ 与 CIR SDE 中对称性存在之间的直接理论联系。作者认为，这种对称性允许在该特定参数范围内构造精确解。

此外，这项工作强调了对称性方法在随机微积分中的实用价值。虽然 SDE 的精确解非常罕见，但本文指出，识别对称性可以通过变量变换（Kozlov 方法）来构建精确解。作者认为，这种方法对于改进数值方案具有重要意义；在数值近似中保持方程的对称性可以提高精度和效率，这一概念在随机文献中虽有提及，但在确定性微分方程领域已取得成功，但在随机领域仍未得到充分利用。

文章最后指出，它重现并确认了近期研究 [7] 中关于 CIR 模型的结果，但是在 SDE 对称性的特定框架内进行的，而非通过 Kolmogorov 方程，并且利用对称性显式地构造了解。作者还提到，在推导显式解的过程中，与局部时间相关的项（Skorokhod 引理）已被舍去。