具有Markov切换Poisson跳的随机微分方程的均方指数稳定性

Mean Square Exponential Stability of Stochastic Differential Equations with Markovian Switching and Poisson Jumps

摘要: 研究一类带Markov切换Poisson跳的随机微分方程的均方指数稳定性。运用Lyapunov稳定性理论、随机分析以及不等式技巧获得了该方程的平凡解是均方指数稳定性的充分条件。最后，给出一个例子说明所得的结果。

Abstract: This paper investigates the mean square exponential stability of stochastic differential equations with Markovian switching and Poisson jumps. By using Lyapunov stability theory, stochastic anal-ysis and inequality techniques, some sufficient conditions are derived to obtain the mean square exponential stability of the trivial solution. At last, an example is presented to illustrate the ob-tained results.

1. 引言

随机微分方程常被用来刻画随机动力系统，且已广泛地应用于物理、生物、医药、社会科学、经济、金融等诸多领域。在随机微分方程的研究中，稳定性是其研究的一个重要课题，关于随机微分方程稳定性的研究已取得一定的成果，见文献 [1] [2] [3] [4] 及其中的参考文献。实际中，许多物理系统在遭遇突变现象(如分支和内部联系紊乱、参数转移等)时，结构会发生随机改变。对于这样的物理系统，人们常用包含连续系统状态和离散系统状态的Markov切换系统来描述，关于带Markov切换的随机微分方程的稳定性研究也引起了学者们的极大兴趣，参阅文献 [5] [6] [7] [8]。

Brown运动是连续的随机过程，然而许多实际系统会遭受跳跃形式的随机突发扰动，如随机故障、地震、海啸等。在这些情形下，不能用Brown运动来刻画这些系统，因此将带跳过程引入到系统中来处理这些实际情形是合理的，见文献 [9] [10] [11]。近年来，关于带跳的随机微分方程的稳定性研究也获得了相应的成果，参阅文献 [12] [13] 及其中的参考文献。

本文研究一类不仅带有Markov切换还带有Poisson跳的随机微分方程的均方指数稳定性。该文结构如下：第2节给出所需的预备知识；第3节利用Lyapunov函数方法证明该类方程的平凡解是均方指数稳定的。第4节通过一个例子说明所得的结果。

2. 预备知识

记一个完备的概率空间为 $\left(\Omega ,\mathcal{F},P\right)$，其滤子 ${\mathcal{F}}_t$ 满足通常的条件，即 ${\left\{{\mathcal{F}}_t\right\}}_{t\ge 0}$ 是右连续的且 ${\mathcal{F}}_0$ 包含所

有的零测集。 ${\left\{r\left(t\right)\right\}}_{t\ge 0}$ 代表的是定义在此完备概率空间上的一个取值为 $S=\left\{1,2,\cdots ,N\right\}$ 的右连续时齐Markov链，其生成元(密度矩阵) $\Gamma ={\left({\gamma }_{ij}\right)}_{N\times N}$ 由转移概率矩阵确定，即

$P\left\{r\left(t+\Delta \right)=j|r\left(t\right)=i\right\}=\{\begin{array}{ll}{\gamma }_{ij}\Delta +o\left(\Delta \right),\hfill & 当i\ne j,\hfill \\ 1+{\gamma }_{ii}\Delta +o\left(\Delta \right),\hfill & 当i=j,\hfill \end{array}$

这里， $\Delta >0$， $\underset{\Delta \to 0}{\lim }o\left(\Delta \right)/\Delta =0$，对于 $i\ne j$， ${\gamma }_{ij}\ge 0$ 表示从状态 $i$ 到状态 $j$ 的转移速率且 ${\gamma }_{ii}=-{\displaystyle \underset{j\ne i}{\sum }{\gamma }_{ij}}$。

作为一个事实， $r(\cdot )$ 的每个样本轨道是一个右连续的阶梯函数，且在 $R^{+}=\left[0,+\infty \right)$ 上的任何一个有限区间里至多存在有限个跳跃点(参阅 [14] )。

考虑如下形式的带Markov切换Poisson跳的随机微分方程：

$\begin{array}{l}\text{d}x\left(t\right)=f\left(x\left(t\right),t,r\left(t\right)\right)\text{d}t+g\left(x\left(t\right),t,r\left(t\right)\right)\text{d}B\left(t\right)\\ +h\left(x\left(t\right),t,r\left(t\right)\right)\text{d}N\left(t\right),t\ge 0,\end{array}$ (1)

初始值 $x\left(0\right)=x_0\in R^m$， $r\left(0\right)=i_0\in S$，其中 $f,g,h:R^m\times R^{+}\times S\to R^m$。 $B\left(t\right)$ 和 $N\left(t\right)$ 分别是定义在该概率空间上的1-维Brown运动和强度为 $\lambda$ 的Poisson过程，这里假设Markov链 $r\left(t\right)$ 与Brown运动 $B\left(t\right)$ 和Poisson过程 $N\left(t\right)$ 是相互独立的。 $\tilde{N}\left(t\right)=N\left(t\right)-\lambda t$ 表示 $N\left(t\right)$ 的补偿Poisson过程。本文假设函数 $f,g,h$ 满足方程(1)解的全局存在唯一性的必需条件。而且，对 $\forall t\ge 0$ 和 $i\in S$，假设 $f\left(0,t,i\right)=g\left(0,t,i\right)=h\left(0,t,\text{i}\right)=0$。因此，当方程(1)的初始值 $x_0=0$ 时，该方程存在平凡解 $x\left(t\right)\equiv 0$。

定义 若存在正常数 $\gamma$ 和 $K$ 使得对于任意给定的初始值 $x\left(0\right)=x_0$ 有

$E{\left|x\left(t\right)\right|}^2\le K{\text{e}}^{-\gamma t}E{\left|x_0\right|}^2,t\ge 0$

成立，则称方程(1)的平凡解是均方指数稳定的。

令 $C^{2,1}\left(R^m\times R^{+};R^{+}\right)$ 表示关于变量 $x$ 二阶连续可导且关于变量 $t$ 一阶连续可导的全体非负函数 $V\left(x,t\right)$ 的集合。给定任意的 $V\left(x,t\right)\in C^{2,1}\left(R^m\times R^{+};R^{+}\right)$，定义算子

$\begin{array}{c}LV\left(x,t,i\right)=V_t\left(x,t\right)+V_x\left(x,t\right)f\left(x,t,i\right)\\ +\frac{1}{2}\text{trace}\left[g^{\text{T}}\left(x,t,i\right)V_{xx}\left(x,t\right)g\left(x,t,i\right)\right]\\ +\lambda \left[V\left(x+h\left(x,t,i\right),t\right)-V\left(x,t\right)\right],\end{array}$

其中， $V_t\left(x,t\right)=\frac{\partial v\left(x,t\right)}{\partial t}$， $V_{xx}\left(x,t\right)={\left(\frac{{\partial }^{2}v\left(x,t\right)}{\partial x_i\partial x_j}\right)}_{m\times m}$， $V_x\left(x,t\right)=\left(\frac{\partial v\left(x,t\right)}{\partial x_1},\frac{\partial v\left(x,t\right)}{\partial x_2},\cdots ,\frac{\partial v\left(x,t\right)}{\partial x_m}\right)$。

3. 主要结果

本节利用Lyapunov函数证明方程(1)的平凡解是均方指数稳定的。为了讨论方程(1)的平凡解的稳定性，对 $\forall i\in S$，假设存在常数 ${\alpha }_{i1}$ 和 ${\alpha }_{i2}$，以及正常数 ${\alpha }_{i3}$ 满足如下条件：

$x^{\text{T}}\left(t\right)f\left(x\left(t\right),t,i\right)+\frac{1}{2}{\left|g\left(x\left(t\right),t,i\right)\right|}^2\le {\alpha }_{i1}{\left|x\left(t\right)\right|}^2,$ (2)

$x^{\text{T}}\left(t\right)h\left(x\left(t\right),t,i\right)\le {\alpha }_{i2}{\left|x\left(t\right)\right|}^2,$ (3)

${\left|h\left(x\left(t\right),t,i\right)\right|}^2\le {\alpha }_{i3}{\left|x\left(t\right)\right|}^2.$ (4)

定理 若条件(2)~(4)成立，且存在正常数 ${\alpha }_0$ 使得

${\alpha }_0+2{\alpha }_{i1}+\lambda \left(2a_{i2}+{\alpha }_{i3}\right)<0$ (5)

成立，则方程(1)的平凡解是均方指数稳定的。

证明 取Lyapunov函数 $V\left(x\left(t\right),t\right)={\text{e}}^{{\alpha }_{0}t}{\left|x\left(t\right)\right|}^2$。则对 $\forall i\in S$，

$\begin{array}{c}LV\left(x\left(t\right),t,i\right)={\alpha }_0{\text{e}}^{{\alpha }_{0}t}{\left|x\left(t\right)\right|}^2+2{\text{e}}^{{\alpha }_{0}t}{x}^{\text{T}}\left(t\right)f\left(x\left(t\right),t,i\right)+{\text{e}}^{{\alpha }_{0}t}{\left|g\left(x\left(t\right),t,i\right)\right|}^2\\ +{\text{e}}^{{\alpha }_{0}t}\lambda \left(2x^{\text{T}}\left(t\right)h\left(x\left(t\right),t,i\right)+{\left|h\left(x\left(t\right),t,i\right)\right|}^2\right)\\ ={\text{e}}^{{\alpha }_{0}t}[{\alpha }_0{\left|x\left(t\right)\right|}^2+2x^{\text{T}}\left(t\right)f\left(x\left(t\right),t,i\right)+{\left|g\left(x\left(t\right),t,i\right)\right|}^2\\ +\lambda \left(2x^{\text{T}}\left(t\right)h\left(x\left(t\right),t,i\right)+{\left|h\left(x\left(t\right),t,i\right)\right|}^2\right)].\end{array}$

由条件(2)~(4)可得

$LV\left(x\left(t\right),t,i\right)\le {\text{e}}^{{\alpha }_{0}t}\left[{\alpha }_0+2{\alpha }_{i1}+\lambda \left(2a_{i2}+{\alpha }_{i3}\right)\right]{\left|x\left(t\right)\right|}^2.$

利用条件(5)进一步可得

$ELV\left(x\left(t\right),t,i\right)\le 0.$

因

$EV\left(x\left(t\right),t\right)-EV\left(x\left(0\right),0\right)\le {\displaystyle {\int }_0^{t}ELV}\left(x\left(s\right),s,r\left(s\right)\right)\text{d}s\le 0,$

所以

$EV\left(x\left(t\right),t\right)\le EV\left(x\left(0\right),0\right).$

注意 $V\left(x\left(t\right),t\right)={\text{e}}^{{\alpha }_{0}t}{\left|x\left(t\right)\right|}^2$，从而

${\text{e}}^{{\alpha }_{0}t}E{\left|x\left(t\right)\right|}^2\le E{\left|x_0\right|}^2,$

即

$E{\left|x\left(t\right)\right|}^2\le {\text{e}}^{-{\alpha }_{0}t}E{\left|x_0\right|}^2.$

证毕。

4. 例子

考虑如下带Markov切换Poisson跳的随机微分方程：

$\text{d}x\left(t\right)=f\left(x\left(t\right),t,r\left(t\right)\right)\text{d}t+g\left(x\left(t\right),t,r\left(t\right)\right)\text{d}B\left(t\right)+h\left(x\left(t\right),t,r\left(t\right)\right)\text{d}N\left(t\right),$ (6)

其中， $B\left(t\right)$ 和 $N\left(t\right)$ 分别是1-维Brown运动和强度 $\lambda =0.5$ 的Poisson过程， $r\left(t\right)$ 是取值为 $S=\left\{1,2\right\}$ 的Markov链且其生成元如下：

$\text{Γ}=\left(\begin{array}{cc}-1& 1\\ 2& -2\end{array}\right).$

方程(6)中的系数参数如下：

当 $i=1$ 时， $f\left(x,t,1\right)=-2x,g\left(x,t,1\right)=x,h\left(x,t,1\right)=x$ ；

当 $i=2$ 时， $f\left(x,t,1\right)=-3x,g\left(x,t,1\right)=2x,h\left(x,t,1\right)=-\frac{1}{2}x$。

对 $\forall i\in S$，可知 $f\left(0,t,i\right)=g\left(0,t,i\right)=h\left(0,t,i\right)=0$。经计算得，当 $i=1$ 时，

$x^{\text{T}}\left(t\right)f\left(x\left(t\right),t,1\right)+\frac{1}{2}{\left|g\left(x\left(t\right),t,1\right)\right|}^2=-\frac{3}{2}{\left|x\left(t\right)\right|}^2,$

$x^{\text{T}}\left(t\right)h\left(x\left(t\right),t,1\right)={\left|x\left(t\right)\right|}^2,$

${\left|h\left(x\left(t\right),t,1\right)\right|}^2={\left|x\left(t\right)\right|}^2.$

当 $i=2$ 时，

$x^{\text{T}}\left(t\right)f\left(x\left(t\right),t,2\right)+\frac{1}{2}{\left|g\left(x\left(t\right),t,2\right)\right|}^2=-{\left|x\left(t\right)\right|}^2,$

$x^{\text{T}}\left(t\right)h\left(x\left(t\right),t,2\right)=-\frac{1}{2}{\left|x\left(t\right)\right|}^2,$

${\left|h\left(x\left(t\right),t,1\right)\right|}^2=\frac{1}{4}{\left|x\left(t\right)\right|}^2.$

取 ${\alpha }_0=1$，由此可知对 $\forall i\in S$，有 ${\alpha }_0+2{\alpha }_{i1}+\lambda \left(2a_{i2}+{\alpha }_{i3}\right)<0$ 成立。根据第三节中的定理可得方程(6)的平凡解是均方指数稳定的。

NOTES

^*通讯作者。