第二章 随机过程
1. 随机过程基本概念及其统计特性
具有随机相位 ϕ \phi ϕ的交流电压波形: V ( t ) = A c o s ( ω t + ϕ ) V(t)=Acos(\omega t+\phi) V(t)=Acos(ωt+ϕ).
定义:对于每一个特定时间 t i t_i ti, X ( t i , s ) X(t_i,s) X(ti,s)都是随机变量,则称 X ( t i , s ) X(t_i,s) X(ti,s)为随机函数。
1.1. 随机事件和概率
- 我们把任意可以重复进行的过程称为一个试验(experiment)。如果试验的结果不止一个,成为随机试验(random experiment)。
- 随机试验的所有可能的结果(outcom)构成的集合,称为样本空间(sample space)。样本空间的子集称为随机事件(random event)。只有一个元素的子集称为基本事件(elementary event)。
- 刻画时间发生可能性大小的数量指标,称为事件的概率。
1.2. 随机过程的分类
- 连续型随机:时间连续,状态连续
- 离散型随机:时间连续,状态离散
- 连续随机:时间离散,状态连续
- 离散随机(随机数字序列,数字信号):状态量化后的连续随机序列,即时间离散,状态离散
1.3. 随机过程的概率分布
1.一维概率分布
F
X
(
x
;
t
1
)
=
P
{
X
(
t
1
)
≤
x
}
F_X(x;t_1)=P\{X(t_1)\leq x\}
FX(x;t1)=P{X(t1)≤x}
f X ( x ; t 1 ) = ∂ F X ( x ; t 1 ) ∂ x f_{X}\left(x ; t_{1}\right)=\frac{\partial F_{X}\left(x ; t_{1}\right)}{\partial x} fX(x;t1)=∂x∂FX(x;t1)
连续性随机变量的分布函数为连续函数,因此 P ( x = a ) = 0 P(x=a)=0 P(x=a)=0,一件事情的概率为0,并不代表这件事情不会发生。
2.二维概率分布
F
X
(
x
1
,
x
2
;
t
1
,
t
2
)
=
P
{
X
(
t
1
)
≤
x
1
,
X
(
t
2
)
≤
x
2
}
F_X(x_1,x_2;t_1,t_2)=P\{X(t_1)\leq x_1,X(t_2)\leq x_2\}
FX(x1,x2;t1,t2)=P{X(t1)≤x1,X(t2)≤x2}
f X ( x 1 , x 2 ; t 1 , t 2 ) = ∂ 2 F X ( x 1 , x 2 ; t 1 , t 2 ) ∂ x 1 ∂ x 2 f_X(x_1,x_2;t_1,t_2)=\frac{\partial^2F_X(x_1,x_2;t_1,t_2)}{\partial x_1\partial x_2} fX(x1,x2;t1,t2)=∂x1∂x2∂2FX(x1,x2;t1,t2)
1.4. 随机过程的数字特征
对于一个随机过程,要研究它的变化规律,常常需要建立起它的“函数关系”,也就是建立随机过程的多维分布。因为堆积过程的多维分布可以比较全面地描述随机过程的整个变化规律的统计特性,弹药建立过程的分布函数一般比较复杂,使用也不便,甚至不可能。
因此,在许多实际应用中,当随机过程的“函数关系”不好确定时,我们往往可以退而求其次,像引入随机变量的数则特征一样,引入随机过程的数字特征
1.4.1. 数学期望(均值)
m X ( T ) = E [ X ( t ) ] = ∫ − ∞ ∞ x f X ( x ; t ) d x m_X(T)=E[X(t)]=\int^{\infty}_{-\infty}xf_X(x;t)dx mX(T)=E[X(t)]=∫−∞∞xfX(x;t)dx
1.4.2. 均方值和方差
二阶原点矩:
Ψ
X
2
(
t
)
=
E
[
X
2
(
t
)
]
=
∫
−
∞
∞
x
2
f
X
(
x
;
t
)
d
x
\varPsi_X^2(t)=E[X^2(t)]= \int^{\infty}_{-\infty}x^2f_X(x;t)dx
ΨX2(t)=E[X2(t)]=∫−∞∞x2fX(x;t)dx
二阶中心矩:
σ
X
2
(
t
)
=
D
[
X
(
t
)
]
=
E
{
[
X
(
t
)
]
2
}
=
E
{
[
X
(
t
)
−
m
X
(
t
)
]
2
}
=
E
[
X
2
(
t
)
]
−
[
m
X
(
t
)
]
2
\sigma^2_X(t)=D[X(t)]=E\{[X(t)]^2\}\\=E\{[X(t)-m_X(t)]^2\}=E[X^2(t)]-[m_X(t)]^2
σX2(t)=D[X(t)]=E{[X(t)]2}=E{[X(t)−mX(t)]2}=E[X2(t)]−[mX(t)]2
1.4.2.1. 数字特征的性质
- η = g ( ξ ) : E ( η ) = E ( g ( ξ ) ) = ∫ − ∞ ∞ g ( x ) f ( x ) d x \eta=g(\xi):E(\eta)=E(g(\xi))=\int^{\infty}_{-\infty}g(x)f(x)dx η=g(ξ):E(η)=E(g(ξ))=∫−∞∞g(x)f(x)dx特别的,若 η = a ξ + b , 则 E ( η ) = a E ( ξ ) + b , 且 D ( η ) = a 2 D ( ξ ) . \eta=a\xi+b,则E(\eta)=aE(\xi)+b,且D(\eta)=a^2D(\xi). η=aξ+b,则E(η)=aE(ξ)+b,且D(η)=a2D(ξ).
- μ n = E [ ξ − E ( ξ ) ] 2 = ∑ k = 1 n ( − 1 ) n − k C n k ν 1 n − k ν k \mu_n=E[\xi-E(\xi)]^2= \sum_{k = 1}^{n} (-1)^{n-k}C^k_n\nu_1^{n-k}\nu_k μn=E[ξ−E(ξ)]2=k=1∑n(−1)n−kCnkν1n−kνk特别的, D ( ξ ) = E ( ξ 2 ) − [ E ( ξ ) ] 2 . 又 D ( ξ ) ≥ 0 D(\xi)=E(\xi^2)-[E(\xi)]^2.又D(\xi)\geq 0 D(ξ)=E(ξ2)−[E(ξ)]2.又D(ξ)≥0,由此得 ( ∫ − ∞ ∞ x f ( x ) d x ) 2 ≤ ∫ − ∞ ∞ x 2 f ( x ) d x (\int^{\infty}_{-\infty}xf(x)dx)^2\leq \int^{\infty}_{-\infty}x^2f(x)dx (∫−∞∞xf(x)dx)2≤∫−∞∞x2f(x)dx
- E [ ξ − E ( ξ ) ] 2 ≤ E [ ( ξ − x ) 2 ] E[\xi-E(\xi)]^2\leq E[(\xi-x)^2] E[ξ−E(ξ)]2≤E[(ξ−x)2]
1.4.3. 自相关、协方差函数
随机过程的均值、方差描述了随机过程在各个鼓励时刻的重要数字特征值,但他们不能反映随机过程的内在联系。因此便引入了自相关和协方差函数。
1.4.3.1. 自相关函数(二阶联合原点矩)
R X ( t 1 , t 2 ) = E [ X ( t 1 ) X ( t 2 ) ] = ∫ − ∞ ∞ ∫ − ∞ ∞ x 1 x 2 f ( x 1 , x 2 ; t 1 , t 2 ) d x 1 d x 2 R_X(t_1,t_2)=E[X(t_1)X(t_2)]\\=\int^{\infty}_{-\infty}\int^{\infty}_{-\infty}x_1x_2f(x_1,x_2;t_1,t_2)dx_1dx_2 RX(t1,t2)=E[X(t1)X(t2)]=∫−∞∞∫−∞∞x1x2f(x1,x2;t1,t2)dx1dx2
- 自相关函数反映随机过程在两个不同时刻的相关程度。
- 自相关函数可正可负,其绝对值越大,表示相关性越强。一般来说,时间相隔越远,相关性越弱,自相关函数的绝对值越小。
- 当两个时刻重合时,其相关性是最强的。
1.4.3.2. 协方差(二阶联合中心距)
K X ( t 1 , t 2 ) = E [ ( X ( t 1 ) − m X ( t 1 ) ) ( X ( t 2 − m X ( t 2 ) ) ] = R X ( t 1 , t 2 ) − m X ( t 1 ) m X ( t 2 ) K_X(t_1,t_2)=E[(X(t_1)-m_X(t_1))(X(t_2-m_X(t_2))]\\=R_X(t_1,t_2)-m_X(t_1)m_X(t_2) KX(t1,t2)=E[(X(t1)−mX(t1))(X(t2−mX(t2))]=RX(t1,t2)−mX(t1)mX(t2)
- 协方差函数反应随机过程两个时刻的起伏值之间的相关程度
当 t 1 = t 2 t_1=t_2 t1=t2时,
R X ( t 1 , t 2 ) = R X ( t , t ) = E [ X 2 ( t ) ] R_X(t_1,t_2)=R_X(t,t)=E[X^2(t)] RX(t1,t2)=RX(t,t)=E[X2(t)]
K X ( t 1 , t 2 ) = D [ X ( t ) ] = E ( X 2 ( t ) ) − m X 2 ( t ) K_X(t_1,t_2)=D[X(t)]=E(X^2(t))-m_X^2(t) KX(t1,t2)=D[X(t)]=E(X2(t))−mX2(t)
若 X ( t 1 ) , X ( t 2 ) X(t_1),X(t_2) X(t1),X(t2)独立
- R X ( t 1 , t 2 ) = E [ X ( t 1 ) X ( t 2 ) ] = E [ X ( t 1 ) ] E [ X ( t 2 ) ] R_X(t_1,t_2)=E[X(t_1)X(t_2)]=E[X(t_1)]E[X(t_2)] RX(t1,t2)=E[X(t1)X(t2)]=E[X(t1)]E[X(t2)]
- K X ( t 1 , t 2 ) = 0 , 即 X ( t 1 ) , X ( t 2 ) 不 相 关 K_X(t_1,t_2)=0,即X(t_1),X(t_2)不相关 KX(t1,t2)=0,即X(t1),X(t2)不相关
1.4.3.3. 相关系数(normalized covariance)
为去除起伏值强度的影响,对协方差函数做归一化处理
r X ( t 1 , t 2 ) = K X ( t 1 , t 2 ) σ X ( t 1 ) σ X ( t 2 ) r_X(t_1,t_2)=\frac{K_X(t_1,t_2)}{\sigma_X(t_1)\sigma_X(t_2)} rX(t1,t2)=σX(t1)σX(t2)KX(t1,t2)
用来描述 X ( t 1 ) , X ( t 2 ) X(t_1),X(t_2) X(t1),X(t2)的线性相关程度。measure the degree of correlation
- ∣ r X ( t 1 , t 2 ) ∣ ≤ 1 |r_X(t_1,t_2)|\leq 1 ∣rX(t1,t2)∣≤1
- r X ( t 1 , t 2 ) = ± 1 r_X(t_1,t_2)=\pm 1 rX(t1,t2)=±1,表示二者线性相关
- r X ( t 1 , t 2 ) = 0 , r_X(t_1,t_2)=0, rX(t1,t2)=0,表明二者不相关
r X Y = 0 ⟺ K X Y = 0 ⟺ E ( X Y ) = E ( X ) E ( Y ) ⟺ D ( X ± Y ) = D ( X ) + D ( Y ) r_{XY}=0\\\Longleftrightarrow K_{XY}=0\\\Longleftrightarrow E(XY)=E(X)E(Y)\\\Longleftrightarrow D(X\pm Y)=D(X)+D(Y) rXY=0⟺KXY=0⟺E(XY)=E(X)E(Y)⟺D(X±Y)=D(X)+D(Y)
不相关与相互独立:
不相关
⇏
\nRightarrow
⇏ 相互独立
相互独立
⇒
\Rightarrow
⇒ 不相关
正交随机变量:
定义: 当 E [ X Y ] = 0 E[XY]=0 E[XY]=0,则称随机变量 X , Y X,Y X,Y是正交的。
命题: 若
E
[
X
]
=
0
E[X]=0
E[X]=0或
E
[
Y
]
=
0
E[Y]=0
E[Y]=0至少有一成立,则
X
,
Y
X,Y
X,Y不相关与
X
,
Y
X,Y
X,Y正交等价。
K
X
Y
=
E
[
X
Y
]
−
E
[
X
]
E
[
Y
]
=
E
[
X
Y
]
K_{XY}=E[XY]-E[X]E[Y]=E[XY]
KXY=E[XY]−E[X]E[Y]=E[XY]
1.5. 随机信号的特征函数
1.5.1. 一维特征函数
C X ( u 1 ; t 1 ) = E [ e j u 1 X ( t 1 ) ] = ∫ − ∞ ∞ f X ( x 1 ; t 1 ) e j u 1 x 1 d x 1 C_X(u_1; t_1)=E[e^{ju_1X(t_1)}]=\int^{\infin}_{-\infty}f_X(x_1;t_1)e^{ju_1x_1}dx_1 CX(u1;t1)=E[eju1X(t1)]=∫−∞∞fX(x1;t1)eju1x1dx1
f X ( x 1 ; t 1 ) = 1 2 π ∫ − ∞ ∞ C x ( u 1 ; t 1 ) e − j u 1 x d u 1 f_X(x_1;t_1)=\frac{1}{2\pi}\int^{\infin}_{-\infin}C_x(u_1;t_1)e^{-ju_1x}du_1 fX(x1;t1)=2π1∫−∞∞Cx(u1;t1)e−ju1xdu1
1.5.1.1. 特征函数的性质
- ∣ C ( u ) ∣ ≤ C ( 0 ) = 1 ; C X ( − u ) = C X ( u ) ‾ |C(u)|\leq C(0)=1; C_X(-u)=\overline{C_X(u)} ∣C(u)∣≤C(0)=1;CX(−u)=CX(u)
- 若 η = a ξ + b \eta =a\xi +b η=aξ+b,则 C η = e j u b C ξ ( a u ) C_{\eta}=e^{jub}C_{\xi}(au) Cη=ejubCξ(au)
- E ( ξ k ) = j − k C ( k ) ( 0 ) E(\xi^{k})=j^{-k}C^{(k)}(0) E(ξk)=j−kC(k)(0)
- 密度函数与特征函数具有一一对应的关系
C ( u ) = ∫ − ∞ ∞ f ( x ) e j u x d x 1 ⟺ f ( x ) = 1 2 π ∫ − ∞ ∞ C ( u ) e − j u x d u C(u)=\int^{\infin}_{-\infty}f(x)e^{jux}dx_1\Longleftrightarrow f(x)=\frac{1}{2\pi}\int^{\infin}_{-\infin}C(u)e^{-jux}du C(u)=∫−∞∞f(x)ejuxdx1⟺f(x)=2π1∫−∞∞C(u)e−juxdu
1.5.1.2. 二维特征函数
C X ( u 1 , u 2 ; t 1 , t 2 ) = E [ e j [ u 1 X ( t 1 ) + u 2 X ( t 2 ) ] ] = ∫ − ∞ ∞ f X ( x 1 , x 2 ; t 1 , t 2 ) d x 1 e j [ u 1 x 1 + u 2 x 2 ] d x 1 d x 2 C_X(u_1,u_2; t_1,t_2)=E[e^{j[u_1X(t_1)+u_2X(t_2)]}]=\int^{\infin}_{-\infty}f_X(x_1,x_2;t_1,t_2)dx_1e^{j[u_1x_1+u_2x_2]}dx_1dx_2 CX(u1,u2;t1,t2)=E[ej[u1X(t1)+u2X(t2)]]=∫−∞∞fX(x1,x2;t1,t2)dx1ej[u1x1+u2x2]dx1dx2
X 1 , X 2 独立 ⟺ C X 1 , X 2 ( u 1 , u 2 ) = C X 1 ( u 1 ) C X 2 ( u 2 ) X_1, X_2\text{独立}\Longleftrightarrow C_{X_1,X_2}(u_1,u_2)=C_{X_1}(u_1)C_{X_2}(u_2) X1,X2独立⟺CX1,X2(u1,u2)=CX1(u1)CX2(u2)
2. 随机过程的微分和积分
- 具有二阶矩的随机过程可以建立随机微积分的概念,对于研究常用的平稳过程很有用。
- 为定义二阶矩过程的连续性、导数和积分概念,需要利用随机变量的极限,下面主要研究二阶矩过程在均方极限意义下的随机微积分。
- 在普通的微积分中,连续、导数和积分等概念都是建立在极限概念的基础上。以随机序列极限为基础,研究分析随机过程的连续、导数和积分等概念和性质。
2.1. 随机变量序列的均方极限
2.1.1. 定义1(二阶矩变量空间)
称定义在概率空间
(
Ω
,
F
,
P
)
(\Omega,F,P)
(Ω,F,P)上的具有二阶矩的随机变量的全体所组成的集合
H
=
{
X
∣
E
[
∣
X
∣
2
]
<
+
∞
}
H=\{X|E[|X|^2]<+\infty\}
H={X∣E[∣X∣2]<+∞}
为二阶矩空间变量。
引理: 设 X , Y ∈ H , X,Y\in H, X,Y∈H,则对任意的复数 a , b a,b a,b 有 a X + b Y ∈ H aX+bY\in H aX+bY∈H
2.1.2. 定义2(均方极限)
设
X
,
X
n
∈
H
,
n
=
1
,
2
,
⋯
。
X,X_n\in H, n=1,2,\cdots。
X,Xn∈H,n=1,2,⋯。如果
lim
n
→
∞
E
∣
X
n
−
X
∣
2
=
0
\lim_{n\rightarrow \infty}E|X_n-X|^2=0
n→∞limE∣Xn−X∣2=0
则称
{
X
n
,
n
=
1
,
2
,
⋯
}
\{X_n,n=1,2,\cdots\}
{Xn,n=1,2,⋯}均收敛于
X
X
X。或称
X
X
X为
{
X
n
,
n
=
1
,
2
,
⋯
}
\{X_n,n=1,2,\cdots\}
{Xn,n=1,2,⋯}的均方极限,记为:
l
.
i
.
m
n
→
∞
X
n
=
X
l.i.m_{n\rightarrow\infty}X_n=X
l.i.mn→∞Xn=X
l . i . m l.i.m l.i.m是英文 limit in mean square 的缩写,表示均方意义下的极限
2.1.2.1. 性质
( 1 ) l . i . m n → ∞ c n = lim n → ∞ c n = c (1)\quad l.i.m_{n\rightarrow\infty}c_n=\lim_{n \rightarrow \infty }c_n=c (1)l.i.mn→∞cn=limn→∞cn=c
( 2 ) l . i . m n → ∞ Z = Z (2)\quad l.i.m_{n\rightarrow\infty}Z=Z (2)l.i.mn→∞Z=Z
( 3 ) l . i . m n → ∞ ( c n Z ) = c Z (3)\quad l.i.m_{n\rightarrow\infty}(c_nZ)=cZ (3)l.i.mn→∞(cnZ)=cZ
2.1.2.2. 定理
设${X_n,n=1,2\cdots},{Y_n,n=1,2,\cdots} \subset H. 其 中 其中 其中X,Y\in H$.
且
l
.
i
.
m
n
→
∞
(
X
n
)
=
X
,
l
.
i
.
m
→
∞
(
Y
n
)
=
Y
,
a
,
b
l.i.m_{n\rightarrow\infty}(X_n)=X,l.i.m_{\rightarrow\infty}(Y_n)=Y,a,b
l.i.mn→∞(Xn)=X,l.i.m→∞(Yn)=Y,a,b为常数,则
l
.
i
.
m
n
→
∞
(
a
X
n
+
b
Y
n
)
=
a
X
+
b
Y
;
★
l.i.m_{n\rightarrow\infty}(aX_n+bY_n)=aX+bY;\bigstar
l.i.mn→∞(aXn+bYn)=aX+bY;★
lim m , n → ∞ E ( X m ‾ Y n ) = E [ X ‾ Y ] ; ★ \lim_{m,n\rightarrow\infty}E(\overline{X_m}Y_n )=E[\overline XY];\bigstar m,n→∞limE(XmYn)=E[XY];★
2.1.2.3. 推论
设 { X n , n = 1 , 2 ⋯ } ⊂ H , X ∈ H , 且 \{X_n,n=1,2\cdots\}\subset H,X\in H,且 {Xn,n=1,2⋯}⊂H,X∈H,且 l . i . m n → ∞ X n = X l.i.m_{n\rightarrow \infty}X_n=X l.i.mn→∞Xn=X,则
- l i m n → ∞ E X n = E X = E [ l . i . m n → ∞ X n ] lim_{n\rightarrow\infty}EX_n=EX=E[l.i.m_{n\rightarrow\infty}X_n] limn→∞EXn=EX=E[l.i.mn→∞Xn]
- l i m n → ∞ E ∣ X n ∣ 2 = E ∣ X ∣ 2 = E [ ∣ l . i . m n → ∞ X n ∣ 2 ] lim_{n\rightarrow\infty}E|X_n|^2=E|X|^2=E[|l.i.m_{n\rightarrow\infty}X_n|^2] limn→∞E∣Xn∣2=E∣X∣2=E[∣l.i.mn→∞Xn∣2]
- l i m n → ∞ D X n = D X = D [ l . i . m n → ∞ X n ] lim_{n\rightarrow\infty}DX_n=DX=D[l.i.m_{n\rightarrow\infty}X_n] limn→∞DXn=DX=D[l.i.mn→∞Xn]
2.1.3. 均方收敛判定准则
定理(Cauchy准则)设
{
X
n
,
n
=
1
,
2
⋯
}
⊂
H
,
\{X_n,n=1,2\cdots\}\subset H,
{Xn,n=1,2⋯}⊂H,则
{
X
n
}
\{X_n\}
{Xn}均方收敛的充要条件是:
lim
n
,
m
→
∞
E
∣
X
m
−
X
n
∣
2
=
0
\lim_{n,m\rightarrow\infty}E|X_m-X_n|^2=0
n,m→∞limE∣Xm−Xn∣2=0
定理(Loeve准则)设
{
X
n
,
n
=
1
,
2
⋯
}
⊂
H
,
\{X_n,n=1,2\cdots\}\subset H,
{Xn,n=1,2⋯}⊂H,则
{
X
n
}
\{X_n\}
{Xn}均方收敛的充要条件是:
lim
n
,
m
→
∞
E
[
X
‾
m
X
n
]
=
c
c
<
+
∞
,
是
常
数
\lim_{n,m\rightarrow\infty}E[\overline X_{m}X_n]=c\quad c<+\infty,是常数
n,m→∞limE[XmXn]=cc<+∞,是常数
2.2. 随机过程的微分及积分
2.2.1. 随机过程的微分
数字特征:
随机过程导数的均值等于其均值的导数
E
[
X
˙
(
t
)
]
=
d
d
t
E
[
X
(
t
)
]
E[\dot{X}(t)]=\frac{d}{dt}E[X(t)]
E[X˙(t)]=dtdE[X(t)]
随机过程导数的相关函数等于其相关函数的混合偏导数
R
X
˙
(
t
1
,
t
2
)
=
∂
2
R
X
(
t
1
,
t
2
)
∂
t
1
∂
t
2
R_{\dot{X}}\left(t_{1}, t_{2}\right)=\frac{\partial^{2} R_{X}\left(t_{1}, t_{2}\right)}{\partial t_{1} \partial t_{2}}
RX˙(t1,t2)=∂t1∂t2∂2RX(t1,t2)
2.2.2. 随机过程的积分
Y = ∫ a b X ( t ) d t Y=\int^b_aX(t)dt Y=∫abX(t)dt
数字特征:
随机过程积分的均值等于其均值的积分
E
[
Y
]
=
∫
a
b
m
X
(
t
)
d
t
E[Y]=\int^b_am_X(t)dt
E[Y]=∫abmX(t)dt
随机过程积分的均方值等于其相关函数的二重积分
E
[
Y
2
]
=
∫
a
b
∫
a
b
R
X
(
t
1
,
t
2
)
d
t
1
d
t
2
E[Y^2]=\int^b_a\int^b_aR_X(t_1,t_2)dt_1dt_2
E[Y2]=∫ab∫abRX(t1,t2)dt1dt2
随机过程积分的方差等于其协方差函数的二重积分
σ
Y
2
=
∫
a
b
∫
a
b
K
X
(
t
1
,
t
2
)
d
t
1
d
t
2
\sigma^2_Y=\int^b_a\int^b_aK_X(t_1,t_2)dt_1dt_2
σY2=∫ab∫abKX(t1,t2)dt1dt2
变上限积分:
随机过程: Y ( t ) = ∫ 0 t X ( λ ) d λ Y(t)=\int^t_0X(\lambda)d\lambda Y(t)=∫0tX(λ)dλ
均值:
E [ Y ( t ) ] = ∫ 0 t m X ( λ ) d λ E[Y(t)]=\int^t_0m_X(\lambda)d\lambda E[Y(t)]=∫0tmX(λ)dλ
相关系数:
R
Y
(
t
1
,
t
2
)
=
∫
0
t
1
∫
0
t
2
R
X
(
λ
,
ν
)
d
λ
d
ν
R_Y(t_1,t_2)=\int^{t_1}_0\int^{t_2}_0R_X(\lambda,\nu)d\lambda d\nu
RY(t1,t2)=∫0t1∫0t2RX(λ,ν)dλdν
3. 平稳随机过程和遍历性过程
stationarity:平稳性
ergodicity:遍历性
3.1. 平稳随机过程
平稳随机过程是指统计特性不随时间的推移而变化的随机过程,分为严平稳和宽平稳。
3.1.1. 严平稳随机过程定义
严平稳随机过程: Strict-Sense Stationary:SSS
- SSS一维分布、密度函数和均值都与时间无关
- SSS二维分布、密度函数和自相关函数都与两个时刻的绝对值无关,而仅与两者之差有关
3.1.2. 严平稳随机过程的低阶矩(一、二维概率密度及数字特征)
这只是严平稳随机过程的性质,并不能作为判定其是否为严平稳随机过程的条件。
- E [ X ( t ) ] = m X E[X(t)]=m_X E[X(t)]=mX
- E [ X 2 ( t ) ] = Ψ X 2 E[X^2(t)]=\varPsi_X^2 E[X2(t)]=ΨX2
- D [ X ( t ) ] = σ X 2 D[X(t)]=\sigma^2_X D[X(t)]=σX2
- R X ( t , t + τ ) = R X ( τ ) R_X(t,t+\tau)=R_X(\tau) RX(t,t+τ)=RX(τ)
- K X ( t 1 , t 2 ) = K X ( τ ) = R X ( τ ) − m X 2 K_X(t_1,t_2)=K_X(\tau)=R_X(\tau)-m_X^2 KX(t1,t2)=KX(τ)=RX(τ)−mX2
- K X ( 0 ) = σ X 2 = R X ( 0 ) − m X 2 K_X(0)=\sigma_X^2=R_X(0)-m_X^2 KX(0)=σX2=RX(0)−mX2
3.1.3. 宽平稳随机过程
Wide-Sense Stationary:WSS
3.1.3.1. 定义
若随机过程 X ( t ) X(t) X(t)的数学期望为一常数,其相关函数只与时间间隔 τ = t 2 − t 1 \tau=t_2-t_1 τ=t2−t1有关,且它的均方值有限。即满足:
- E [ X ( t ) ] = m X ( t ) = m X : c o n s t a n t E[X(t)]=m_X(t)=m_X:constant E[X(t)]=mX(t)=mX:constant
- R X ( t , t + τ ) = R X ( τ ) : ∀ t , τ R_X(t,t+\tau )=R_X(\tau):\forall t,\tau RX(t,t+τ)=RX(τ):∀t,τ
- E [ X 2 ( t ) ] < ∞ E[X^2(t)]<\infty E[X2(t)]<∞
则称 X ( t ) X(t) X(t)为宽平稳随机过程(或称广义平稳随机过程)。
Remark:
- 严平稳过程只要二阶矩存在,则它必定也是宽平稳的,反之不成立。
- 宽平稳的正态过程必定也是严平稳的。
- 由于宽平稳随机过程的定义只涉及与一、二维概率密度有关的数字特征,所以,一个严平稳过程只要均方有界,则其必然是广义平稳的。但反之不一定成立。
上述三个公式可作为判定宽平稳随机过程的依据。
3.1.4. 循环平稳随机过程
CycloStationary:CS
- E [ X ( t ) ] = m X ( t ) = m X ( t + k T ) E[X(t)]=m_X(t)=m_X(t+kT) E[X(t)]=mX(t)=mX(t+kT)
- R X ( t , t + τ ) = R X ( t , τ ) = R X ( t + k T , τ ) R_X(t,t+\tau)=R_X(t,\tau)=R_X(t+kT,\tau) RX(t,t+τ)=RX(t,τ)=RX(t+kT,τ)
3.2. 平稳随机过程相关函数的性质
相关函数反映随机过程在两个不同时刻的相关度
R X ( t , t + τ ) = E [ X ( t ) X ( t + τ ) ] = R X ( τ ) R_X(t,t+\tau)=E[X(t)X(t+\tau)]=R_X(\tau) RX(t,t+τ)=E[X(t)X(t+τ)]=RX(τ)
3.2.1. 基本性质
-
R X ( t , t + τ ) = R X ( τ ) = R X ( − τ ) R_X(t,t+\tau)=R_X(\tau)=R_X(-\tau) RX(t,t+τ)=RX(τ)=RX(−τ)表示是一个偶函数。
-
R X ( 0 ) ≥ 0 R_X(0)\geq 0 RX(0)≥0相关函数在 t = 0 t=0 t=0处的值是大于等于0的。
-
K X ( τ ) = R X ( τ ) − m X 2 , σ 2 = R X ( 0 ) − m X 2 K_X(\tau)=R_X(\tau)-m_X^2,\\\sigma^2=R_X(0)-m_X^2 KX(τ)=RX(τ)−mX2,σ2=RX(0)−mX2方差和协方差的计算方法,注意!!!$\bigstar $
-
∣ R X ( τ ) ∣ ≤ R X ( 0 ) |R_X(\tau)|\leq R_X(0) ∣RX(τ)∣≤RX(0)在时间为0处 R X ( t ) R_X(t) RX(t)最大
-
∃ T ≠ 0 → R X ( T ) = R X ( 0 ) ⇒ R X ( τ ) = R X ( τ + k T ) \exists T\neq 0\rightarrow R_X(T)=R_X(0)\\\Rightarrow R_X(\tau)=R_X(\tau+kT) ∃T=0→RX(T)=RX(0)⇒RX(τ)=RX(τ+kT)周期性的问题。
-
若严平稳过程中不含有任何周期分量,则有: lim ∣ τ ∣ → ∞ R X ( τ ) = R X ( ∞ ) = m X 2 \lim_{|\tau|\rightarrow\infty}R_X(\tau)=R_X(\infty)=m_X^2 ∣τ∣→∞limRX(τ)=RX(∞)=mX2
3.2.2. 相关系数 & 相关时间
相关系数:
r
X
(
τ
)
=
K
X
(
τ
)
K
X
(
0
)
=
R
X
(
τ
)
−
m
X
2
σ
X
2
r_X(\tau)=\frac{K_X(\tau)}{K_X(0)}=\frac{R_X(\tau)-m_X^2}{\sigma_X^2}
rX(τ)=KX(0)KX(τ)=σX2RX(τ)−mX2
相关时间:
对于一般的随机国车过而言,随着时间间隔
τ
\tau
τ正大相关程度减弱,因此相关系数也随着减弱,当间隔大到一定程度(假定为
τ
0
\tau_0
τ0),相关系数很小可以认为起伏值不相关了,这个时间就称为相关时间。
3.3. 遍历性过程
遍历性(ergodicity),也称为各态历经性过程。
任何一条样本函数所包含的取值状态与随机过程(任意时刻)所有的状态相同,而且出现的频率与随机过程各状态的概率相同。
辛钦证明:
在具备一定的补充条件下,对平稳随机过程的一个样本函数取时间均值(观察时间足够长),就从概率意义上趋近于此过程的统计(集合)均值。对于这样的随机过程,我们说它具备遍历性。
随机过程的遍历性,可以理解为随机过程的各个样本函数都同样地经历了随机过程的各种可能状态,因此,从随机过程的任何一个样本函数就能得到随机过程的全部统计信息,任何一个样本函数的特性都唔那个充分地代表整个随机过程的特性。
3.3.1. 遍历性过程的定义
- 随机过程的时间平均
- 严遍历性过程的定义: 随机过程的各种时间平均值,在平均时间足够长的条件下都依概率1等于相应的集合平均
3.3.1.1. 严遍历性过程的定义
如果一个随机过程 X ( t ) X(t) X(t),它的各种时间平均(时间足够长)依概率1收敛于相应的集合平均,则称过程 X ( t ) X(t) X(t)具有严格(或狭义)遍历性,并称此过程为严格(或狭义)遍历性过程,简称严遍历过程。
3.3.1.2. 随机过程的时间平均
A ⟨ X ( t ) ⟩ = X ( t ) ‾ = lim T → ∞ 1 2 T ∫ − T T X ( t ) d t A\langle X(t) \rangle =\overline{X(t)}=\lim_{T\rightarrow\infty}\frac{1}{2T}\int^T_{-T}X(t){\rm d}t A⟨X(t)⟩=X(t)=T→∞lim2T1∫−TTX(t)dt
R X ( t , t + τ ) = X ( t ) X ( t + τ ) ‾ = lim T → ∞ 1 2 T ∫ − T T X ( t ) X ( t + τ ) d t \mathfrak{R}_X(t,t+\tau)=\overline{X(t)X(t+\tau)}=\lim_{T\rightarrow\infty}\frac{1}{2T}\int^T_{-T}X(t)X(t+\tau){\rm d}t RX(t,t+τ)=X(t)X(t+τ)=T→∞lim2T1∫−TTX(t)X(t+τ)dt
分别称之为过程 X ( t ) X(t) X(t)的 时间均值 和 时间相关函数 。
3.3.1.3. 宽遍历性过程的定义
设 X ( t ) X(t) X(t)是一个平稳随机过程
(1) 若
A
⟨
X
(
t
)
⟩
=
X
(
t
)
‾
=
E
[
X
(
t
)
]
=
m
X
A\langle X(t) \rangle =\overline{X(t)}=E[X(t)]=m_X
A⟨X(t)⟩=X(t)=E[X(t)]=mX
则称过程 X ( t ) X(t) X(t)的均值具有遍历性。
(2) 若
R X ( t , t + τ ) = X ( t ) X ( t + τ ) ‾ = E [ X ( t ) X ( t + τ ) ] = R X ( τ ) \mathfrak{R}_X(t,t+\tau)=\overline{X(t)X(t+\tau)}=E[X(t)X(t+\tau)]=R_X(\tau) RX(t,t+τ)=X(t)X(t+τ)=E[X(t)X(t+τ)]=RX(τ)
则称过程 X ( t ) X(t) X(t)的自相关函数具有遍历性。若在 τ = 0 \tau=0 τ=0时上式成立,则称过程 X ( t ) X(t) X(t)的均方值具有遍历性。
若(1)和(2)同时满足,则 X ( t ) X(t) X(t)称为宽遍历随机过程。(判定方法)
3.3.2. 随机过程具备遍历性的条件
随机过程均值遍历的充要条件: E [ X ( t ) ] = A [ X ( t ) ] E[X(t)]=A[X(t)] E[X(t)]=A[X(t)]
E [ X ( t ) ] E[X(t)] E[X(t)]与样本无关,是时间的函数,是任意时刻上的概率平均。 E [ X ( t ) ] = m ( c o n s t ) E[X(t)]=m(const) E[X(t)]=m(const)
A [ X ( t ) ] A[X(t)] A[X(t)]与实践无关,与取哪条样本函数有关,是随机变量。 D { A [ X ( t ) ] } = 0 D\{A[X(t)]\}=0 D{A[X(t)]}=0
自相关函数遍历的充分条件: E [ X ( t ) X ( t + τ ) ] = A [ X ( t ) X ( t + τ ) ] E[X(t)X(t+\tau)]=A[X(t)X(t+\tau)] E[X(t)X(t+τ)]=A[X(t)X(t+τ)]