概率论与数理统计(6)——随机过程、Markov链与平稳随机过程

概率论与数理统计第十至十二章——随机过程、Markov链与平稳随机过程

第十章 随机过程

10.1 随机过程的概念

随机过程从动态的角度考察概率问题,研究对象随时间演变。

  • 随机过程

$\{X(e,t),e\in S,t\in T \}$是定义在$S,T$上的二元函数。若此函数对任意固定的$t\in T$,$X(e,t)$是个随机变量, 则称$\{X(e,t),e\in S,t\in T \}$是随机过程

  • 样本函数

对随机过程$\{X(e,t),e\in S,t\in T \}$进行一次试验,即给定e,那么它是t的函数,称为随机过程的样本函数

  • 状态

对任意固定的e和t,$X(e,t)$称为过程的状态

  • 状态空间

$X(e,t)$所有可能的取值称为状态空间。之后将$X(e,t)$简记为$X(t)$。

随机过程可以根据参数集T和任一时刻的状态状态S分为四类:

连续参数 连续型

的随机过程、连续参数离散型的随机过程、离散参数连续型的随机过程、离散参数离散型的随机过程。

10.2 随机过程的统计描述

10.2.1 分布函数

  • 分布函数族

在随机过程$\{X(t),t\in T\}$中,对每一固定的$t\in T$,称

\[F_x(x,t)=P\left\{ X(t)\leq x\right\},x\in \mathbb{R}\]

为该随机过程的一维分布函数,而

\[\left\{ F_x(x,t),t\in T\right\}\]

称为一维函数分布族

  • 推广到n维

n维分布函数

\[F_x(x_1,x_2,...,x_n;t_1,t_2,...,t_n)=P\left\{ X(t_1)\leq x_1,X(t_2)\leq x_2,...,X(t_n)\leq x_n\right\}, x_i\in \mathbb{R}\]

n维分布函数族

\[\left\{F_x(x_1,x_2,...,x_n;t_1,t_2,...,t_n),t_i\in T\right\}\]

10.2.2 数字特征

  • 某一时间点t的数字特征

均值函数

\[\mu_X(t)=E[X(t)]\]

均方值函数

\[\Psi_X^2(t)=E[X^2(t)]\]

方差函数

\[\sigma_X^2(t)=D_X(t)=E\left\{[X(t)-\mu_X(t)]^2\right\}\]

标准差函数

\[\sigma_X(t)=\sqrt{D_X(t)}\]
  • 不同时间点的数字特征

对任意$t_1,t_2\in T$,

(自)相关函数

\[R_{XX}(t_1,t_2)=E[X(t_1)X(t_2)]\]

简记为$R_X(t_1,t_2)$。

(自)协方差函数

\[\begin{aligned} C_{XX}(t_1,t_2)& =Cov[X(t_1)X(t_2)]\\ & =E\left\{[X(t_1)-\mu_X(t_1)][X(t_2)-\mu_X(t_2)]\right\} \end{aligned}\]

简记为$C_X(t_1,t_2)$。

  • 各数字特征之间的关系
\[\Psi_X^2(t)=R_X(t,t)\] \[C_X(t_1,t_2)=R_X(t_1,t_2)-\mu_X(t_1)\mu_X(t_2)\] \[\sigma_X^2(t)=C_X(t,t)=R_X(t,t)-\mu_X^2(t)\]
  • 二阶矩过程

在随机过程$\{X(t),t\in T\}$中,对每一固定的$t\in T$,$E[X^2(t)]$都存在,则称$X(t)$为二阶矩过程

  • 正态过程

若对随机过程$\{X(t),t\in T\}$,任意整数$n$和任意$t_1,t_2,…,t_n\in T$,$(X(t_1),X(t_2),…,X(t_n))$服从$n$维正态分布,则称$\{X(t),t\in T\}$是正态过程。

10.2.3 二维随机过程的分布函数和数字特征

  • 二维随机过程

设$X(t),Y(t)$是依赖于同一参数$t\in T$的随机过程,对于不同的$t\in T$,$(X(t),Y(t))$是不同的二维随机变量,称$\{X(t),Y(t)\}$为二维随机过程。

  • 二维随机过程的n+m维分布函数

n+m维随机变量$(X(t_1),X(t_2),…,X(t_n);Y(t_1’),Y(t_2’),…,Y(t_m’))$的分布函数

\[F(x_1,x_2,...,x_n;t_1,t_2,...,t_n;y_1,y_2,...,y_m;t_1',t_2',...,t_m')\]

称为二维随机过程的n+m维分布函数。

  • 二维随机过程的数字特征

互相关函数

\[R_{XY}(t_1,t_2)=E[X(t_1)Y(t_2)],\quad t_1,t_2\in T\] \[R_{YX}(t_1,t_2)=E[Y(t_1)X(t_2)],\quad t_1,t_2\in T\]

互协方差函数

\[C_{XY}(t_1,t_2)=R_{XY}(t_1,t_2)-\mu_{X}(t_1)\mu_Y(t_2)\] \[C_{YX}(t_1,t_2)=R_{YX}(t_1,t_2)-\mu_Y(t_1)\mu_X(t_2)\]

若对任意$t_1,t_2\in T$,恒有$C_{XY}(t_1,t_2)=0$,则$X(t)$和$Y(t)$是不相关的。

10.3 泊松过程及维纳过程

10.3.1 独立增量过程

  • 增量

给定二阶矩过程$\{X(t),t\geq 0\}$,对$s,t(0\leq s<t)$,称随机变量

\[X(t)-X(s)\]

为随机过程在区间$(s,t]$上的增量

  • 独立增量过程

如果对任意选定的正整数n和$0\leq t_0<t_1<t_2<\cdots<t_n$,n个增量

\[\left\{X(t_i)-X(t_{i-1})\right\},i=1,2,...,n\]

则称$\{X(t),t\geq 0\}$为独立增量过程

  • 平稳性

如果任给实数h和$0\leq s+h<t+h$,$X(t+h)-X(s+h)$和$X(t)-X(s)$具有相同的分布,则称增量具有平稳性

  • 独立增量过程的性质

1、 $X(t)$的有限维分布函数族可以由增量$X(t)-X(s)(0\leq s<t)$的分布所确定。

2、 设方差函数$D_X(t)$已知,则协方差函数

\[C_X(s,t)=D_X(\min (s,t))\]

10.3.2 泊松过程

以$N(t_0),t\geq 0$表示在时间间隔$(0,t]$内出现的质点数。记$N(t_0,t)=N(t)-N(t_0),0\leq t_0<t$,它表示时间间隔$(t_0,t]$内出现的质点数,其概率记为

\[P_k=(t_0,t)=P\left\{N(t_0,t)=k\right\},\quad k=0,1,2,...\]
  • 泊松过程

计数过程$N(t)$满足如下条件,称作强度为$\lambda$的泊松过程

1、 在不相重叠的区间上的增量具有独立性

2、 对于充分小的$\Delta t$,

\[P_1(t,t+\Delta t)=P\left\{(N(t,t+\Delta t))=1\right\}=\lambda \Delta t+o(\Delta t)\]

其中常数$\lambda$称为$N(t)$的强度。

3、 对于充分小的$\Delta t$,

\[\sum_{j=2}^{\infty}P_j(t,t+\Delta t)=\sum_{j=2}^{\infty}P\left\{N(t,t+\Delta t)=j\right\}=o(\Delta t)\]

4、 $N(0)=0$

  • 分布律

若$N(t)$是强度为$\lambda$的泊松过程,则

\[P_k(t_0,t)=P\left\{N(t_0,t)=k\right\}=\frac{[\lambda (t-t_0)]^k}{k!}e^{-\lambda (t-t_0)},t>t_0,k=0,1,2,...\]

即$N(t_0,t)\sim \pi[\lambda (t-t_0)]$。

因此有泊松过程的另一等价定义

1、 它是独立增量过程

2、 对任意的$t>t_0\geq 0$,增量

\[N(t)-N(t_0)\sim\pi[\lambda (t-t_0)]\]

3、 $N(0)=0$

  • 数字特征

均值函数

\[E[N(t_0,t)]=E\left\{N(t)-N(t_0)\right\}=\lambda (t-t_0)\]

方差函数

\[D[N(t_0,t)]=D\left\{N(t)-N(t_0)\right\}=\lambda (t-t_0)\]

特别地,当$t_0=0$时,$\mu_N(t)=D_N(t)=\lambda t$

协方差函数

\[C_N(s,t)=D_N(\min(s,t))=\lambda \min(s,t),s,t\geq 0\]

自相关函数

\[R_N(s,t)=C_N(s,t)+\mu_N(s)\mu_N(t)=\lambda \min(s,t)+\lambda^2st\]
  • 应用:粒子统计

1、 记$W_n$是第$n$个质点出现的时间,则其概率密度为

\[f_{W_n}(t)= \left\{ \begin{aligned} & \frac{\lambda (\lambda t)^{n-1}e^{-\lambda t}}{(n-1)!},& t>0\\ & 0,& t\leq 0 \end{aligned} \right.\]

服从$\Gamma(n,\lambda)$分布。

特别地,质点首次出现的等待时间服从指数分布

\[f_{W_1}(t)= \left\{ \begin{aligned} & \lambda e^{-\lambda t},& t>0\\ & 0,& t\leq 0 \end{aligned} \right.\]

2、 记$T_i=W_i-W_{i-1},W_0=0$,称之为相继出现的第$i-1$个质点和第$i$个质点的点间间距,则其概率密度为

\[f_{T_i}(t)= \left\{ \begin{aligned} & \lambda e^{-\lambda t},& t>0\\ & 0,& t\leq 0 \end{aligned} \right.\]

两个定理:

定理1 强度为$\lambda$的泊松过程的点间间距是相互独立的随机变量,且服从同一指数分布。

定理2 如果任意相继出现的两个质点的点间间距相互独立,且服从同一个指数分布:

\[f(t)= \left\{ \begin{aligned} & \lambda e^{-\lambda t},& t>0\\ & 0,& t\leq 0 \end{aligned} \right.\]

则质点流构成强度为$\lambda$的泊松过程。

10.3.3 维纳过程

  • 维纳过程

对于二阶矩过程$\{W(t),t\geq 0\}$,如果它满足:

1、 具有独立增量

2、 对任意$t>s\geq 0$,增量$W(t)-W(s)\sim N(0,\sigma^2(t-s))$且$\sigma>0$

3、 $W(0)=0$

则称之为维纳过程

  • 性质

1、 维纳过程是齐次的独立增量过程(当增量具有平稳性时,称为齐次的独立增量过程)

2、 维纳过程是正态过程,故其分布完全由它的均值函数和自协方差函数(即自相关函数)所确定

3、 维纳过程的数字特征

\[\mu_W(t)=E(W(t))=0\] \[D_W(t)=D(W(t))=\sigma^2t\] \[C_W(s,t)=R_W(s,t)=D_W[\min (s,t)]=\sigma^2\min (s,t),s,t>0\]

第十一章 Markov链

11.1 Markov过程及其概率分布

11.1.1 Markov性与Markov过程

  • Markov性

过程在时刻$t_0$所处的状态已知,若它在时刻$t>t_0$所处的状态的条件分布与过程在时刻$t_0$之前所处的状态无关。

  • Markov过程

具有Markov性的过程。即对随机过程$\{X(t),t\in T\}$,对于参数集T中的任意n个数值$t_1<t_2<\cdots<t_n,n\geq 3,t\in T$,

\[P\left\{X(t_n)\leq x_n\mid X(t_1)=x_1,...,X(t_{n-1})=x_{n-1}\right\}=P\left\{X(t_n)\leq x_n\mid X(t_{n-1})=x_{n-1}\right\}\]
  • Markov链

时间和状态都离散的Markov过程称为Markov链,简称马氏链,记为

\[\left\{X_n=X(n),n=0,1,2,...\right\}\]

用条件分布律表示为

\[P\left\{X_{m+n}=a_j\mid X_m=a_i\right\} =: P_{ij}(m,m+n)\]
  • 转移概率

上述$P_{ij}=(m,m+n)$称为马氏链在时间m处于状态$a_i$的条件下,在时间$m+n$转移到状态$a_j$的转移概率。转移概率有一个显然的性质:

\[\sum_{j=1}^{\infty}P_{ij}(m,m+n)=1,j=1,2,...\]
  • 转移概率矩阵
\[P(m,m+n)= \left( \begin{matrix} & P_{11}(m,m+n) & P_{12}(m,m+n) & \cdots\\ & P_{21}(m,m+n) & P_{22}(m,m+n) & \cdots\\ & \vdots & \vdots & \ddots \end{matrix} \right)\]

此矩阵的每一行元素之和为1。

当$P_{ij}(m,m+n)$只与$i,j$及$n$有关时,把它记为$P_{ij}(n)$,称此概率为马氏链的n步转移概率

当转移概率具有平稳性时,称此链是齐次马氏链

  • 初始分布
\[P_j(0)=P\left\{X_0=a_j\right\},a_j\in I,j=1,2,...\]
  • 一维分布
\[P_j(n)=P\left\{X_n=a_j\right\},a_j\in I,n\in T\]

性质

\[\sum_{j=1}^{\infty}P_j(n)=1\]

一个重要的公式

\[P_j(n)=\sum_{i=1}^{\infty}P_i(0)P_{ij}(n),j=1,2,...\]

11.2 多步转移概率的确定

11.2.1 C-K方程

设$\{X_n,n=1,2,…\}$是一齐次马氏链,则对任意的$u,v\in T_1=\{0,1,2,…\}$,有

\[P_{ij}(u+v)=P_{ik}(u)P_{kj}(v),i,j=1,2,...\]

根据矩阵乘法的性质,C-K方程可以写成矩阵乘法形式

\[P(u+v)=P(u)P(v)\]

设$P(n)$是n步转移概率矩阵,则

\[P(n)=P^n\]

齐次马氏链的有限维分布可由初始分布与一步转移概率完全确定。

11.3 遍历性

  • 遍历性

对固定的状态j,不管链在某一时刻从什么状态i出发,经过长时间的转移,到达状态j的概率都趋近于$\pi_j$,这就是遍历性

即满足下式

\[\lim_{n\to \infty}P(n)= \left( \begin{matrix} & \pi_1 & \pi_2 & \cdots \\ & \pi_1 & \pi_2 & \cdots \\ & \vdots & \vdots & \ddots \end{matrix} \right)\]
  • 链的极限分布

由于$\pi_0+\pi_1=1$,所以$(\pi_0,\pi_1)$记为$\pi$,构成一分布律,称为链的极限分布

即满足下式

\[\sum_j\pi_j=1\]
  • 有限链的遍历性的充分条件

设齐次马氏链$\{X_n\geq 1\}$的状态空间为$I=\{a_1,a_2,…,a_N\}$,P是它的一步转移概率矩阵,如果存在正整数m,使对任意的$a_i,a_j\in I$,都有$P_{ij}(m)>0,i,j=1,2,…,N$,则此链具有遍历性,且有极限分布$\pi=(\pi_1,\pi_2,…,\pi_N)$,它是方程组

\[\pi=\pi P\]

的满足条件

\[\pi_j>0,\sum_{j=1}^N\pi_j=1\]

的唯一解。

第十二章 平稳随机过程

12.1 平稳随机过程的概念

  • 平稳性

对任意的$n(n=1,2,…),t_1,t_2,…,t_n\in T$和任意实数$h$,当$t_1+h,t_2+h,…,t_n+h\in T$时,

\[(X(t_1),X(t_2),...,X(t_n))\]

\[(X(t_1+h),X(t_2+h),...,X(t_n+h))\]

具有相同的分布函数,则称随机过程具有平稳性,称此过程为严平稳随机过程,简称严平稳过程

  • 数字特征

$X(t)$与$X(t+h)$同分布,取$h=-t$,则$X(t)$与$X(0)$同分布,从而有相同的数学期望(均值函数)

\[\mu_X(t)=E[X(t)]=E[X(0)]=:\mu_X为常数\]

取$h=-t_1$,则$(X(t_1),X(t_2))$与$(X(0),X(t_2-t_1))$同分布,从而自相关函数仅是时间差$t_2-t_1=\tau$的函数。

\[R_X(t_1,t_2)=E[X(t_1)X(t_2)]=R_X(0,t_2-t_1)=:R_X(t_2-t_1)=R_X(\tau)\]

协方差函数

\[C_X(\tau)=R_X(\tau)-\mu_X^2\]

方差函数:这里$t_2=t_1=t$,因此$\tau=0$

\[D_X(t)=C_X(0)=R_X(0)-\mu_X^2\]
  • 宽平稳过程

给定二阶矩过程$\{X(t),t\in T\}$,如果对任意的$t,t+\tau\in T$,都有

\[E[X(t)]=\mu_X为常数\] \[E[X(t)X(t+\tau)]=R_X(\tau)\]

则称$\{X(t),t\in T\}$为宽平稳过程。

  • 平稳相关(联合平稳)

若两个平稳过程的互相关函数也只是时间差的函数,记为$R_{XY}(\tau)$,即

\[R_{XY}(t,t+\tau)=E[X(t)Y(t+\tau)]=R_{XY}(\tau)\]

则称$X(t)$和$Y(t)$是平稳相关的,或称这两个过程是联合(宽)平稳的。

12.2 各态历经性

考虑到平稳过程的统计特性不随时间的推移而变化,可以定义随机积分。

老师突然告诉我这一章不考了,摆!