[论文解读] Long time behavior of diffusions with Markov switching

Jean-Baptiste Bardet, Hélène Guérin+1 more|ArXiv.org|Dec 16, 2009

Stochastic processes and statistical mechanics参考文献 7被引用 37

一句话总结

本文建立了带马氏切换的 Ornstein-Uhlenbeck 扩散过程不变分布尾部行为的三支分类，根据模型参数识别出三种不同的尾部行为——重尾、类似指数尾或类似高斯尾。通过关联切换过程与漂移和扩散系数的矩阵的谱分析，提供了收敛到平衡态的定量速率。

ABSTRACT

Let $Y$ be an Ornstein-Uhlenbeck diffusion governed by an ergodic finite state Markov process $X$: $dY_t=-λ(X_t)Y_tdt+σ(X_t)dB_t$, $Y_0$ given. Under ergodicity condition, we get quantitative estimates for the long time behavior of $Y$. We also establish a trichotomy for the tail of the stationary distribution of $Y$: it can be heavy (only some moments are finite), exponential-like (only some exponential moments are finite) or Gaussian-like (its Laplace transform is bounded below and above by Gaussian ones). The critical moments are characterized by the parameters of the model.

研究动机与目标

• 表征由有限状态、遍历的马氏切换过程驱动的 Ornstein-Uhlenbeck 扩散过程的长期行为。
• 解决关于不变分布尾部行为是否严格为二分法（如先前所建议）还是允许第三种行为的开放问题。
• 通过耦合技术与由生成元导出的矩阵的谱分析，建立收敛到不变测度的定量速率。
• 根据模型参数，明确识别不变分布呈现重尾、类似指数尾或类似高斯尾的精确条件。
• 通过引入新参数 $\kappa$（由涉及切换生成元与漂移系数的矩阵的谱半径定义），拓展先前关于矩存在性与遍历性的结果。

提出的方法

• 将过程建模为跳跃扩散过程，其中漂移 $\lambda(X_t)$ 与扩散 $\sigma(X_t)$ 根据有限状态、不可约的连续时间马氏过程 $X_t$ 进行切换。
• 利用 SDE 的显式解，将 $Y_t$ 表示为具有随机时变系数的随机积分，从而实现对不变测度的分析。
• 基于切换过程的生成元 $A$ 与漂移 $\lambda$，定义矩阵 $A_p = A - p\Lambda$，并引入 $\eta_p = -\max\{\mathrm{Re}(\gamma) \mid \gamma \in \mathrm{Spec}(A_p)\}$ 以表征尾部行为。
• 引入 $\kappa = \sup\{p \geq 0 \mid \eta_p > 0\}$ 作为决定尾部行为区别的关键参数，满足 $\kappa \leq \infty$，且在 $(0,\kappa)$ 上 $\eta_p > 0$，在 $(\kappa,\infty)$ 上 $\eta_p < 0$，且当 $p \uparrow \kappa$ 时 $\eta_p \to 0$。
• 对具有不同初始分布的两个过程应用耦合技术，对切换过程使用粘性耦合，并通过共同的布朗运动耦合扩散路径。
• 利用谱参数 $\eta_p$ 与耦合时间，推导出 Wasserstein 距离 $W_p(\mathcal{L}(Y_t), \mathcal{L}(\tilde{Y}_t))$ 的指数上界，从而得到阶为 $\exp(-c t)$ 的收敛速率，其中 $c = \frac{\gamma \eta_p}{\gamma + s \eta_p}$。

实验结果

研究问题

• RQ1 扩散过程在马氏切换下，其不变分布呈现重尾、类似指数尾或类似高斯尾的精确条件是什么？
• RQ2 关键参数 $\kappa = \sup\{p \geq 0 \mid \eta_p > 0\}$ 如何与矩的存在性及不变测度的尾部行为相关联？
• RQ3 此类扩散过程收敛到平衡态的速率是否可以显式量化为指数速率，且这些速率如何依赖于切换动力学？
• RQ4 先前已知的重尾与轻尾二分法是否不足，若有，第三种行为是什么？
• RQ5 扩散系数 $\sigma(x)$ 与跳跃速率 $a(x)$ 如何影响尾部行为与收敛速度？

主要发现

• 扩散过程的不变分布呈现出尾部行为的三支分类：当 $\underline{\lambda} < 0$ 且 $\kappa < \infty$ 时为重尾；当 $\underline{\lambda} \geq 0$ 且 $\kappa < \infty$ 时为类似指数尾；当 $\underline{\lambda} \geq 0$ 且 $\kappa = \infty$ 时为类似高斯尾。
• $\kappa = \sup\{p \geq 0 \mid \eta_p > 0\}$ 是决定尾部行为区别的关键参数，其中 $\eta_p$ 在 $(0,\kappa)$ 上连续且为正，在 $(\kappa, \infty)$ 上为负，且当 $p \uparrow \kappa$ 时 $\eta_p \to 0$。
• 当 $\underline{\lambda} < 0$ 时，尾部呈多项式衰减：$t^\kappa \nu((t,\infty)) \to C > 0$ 当 $t \to \infty$，其中 $\kappa$ 是方程 $\rho(M_p) = 1$ 的唯一解，$M_p$ 由切换矩阵与漂移定义。
• 当 $\underline{\lambda} \geq 0$ 时，若 $\kappa = \infty$，则不变测度具有所有阶矩，且尾部为类似高斯尾：$\nu$ 的 Laplace 变换被高斯函数上下界控制。
• Wasserstein 距离 $W_p(\mathcal{L}(Y_t), \mathcal{L}(\tilde{Y}_t))$ 以速率 $\frac{\gamma \eta_p}{\gamma + s \eta_p}$ 指数衰减，其中 $\gamma$ 与 $s$ 分别与耦合时间及矩阶相关。
• 收敛到平衡态被定量控制：$W_p(\mathcal{L}(Y_t), \mathcal{L}(\tilde{Y}_t))^p \leq C_1 e^{-c t} + C_2 e^{-\eta_p t}$，其中 $c = \frac{\gamma \eta_p}{\gamma + s \eta_p}$，表明速率依赖于 $\eta_p$ 与耦合参数。