随机过程

随机过程

1 概率论基础

(1)随机变量函数的概率密度

对于任意的单调函数 g ( x ) g(x) g(x),都有
f Y ( y ) = f X ( x ) ∣ J ∣ x = g − 1 ( y ) ( J = d x d y ) f_Y(y)=f_X(x)|J|_{x=g^{-1}(y)} (J=\frac{dx}{dy}) fY(y)=fX(x)Jx=g1(y)(J=dydx)
对于非单调函数,可以根据单调性分段。


例: 考虑一个平方律检波的例子,假定输入输出的关系为
Y = b X 2 Y=bX^2 Y=bX2
求Y的概率密度
解: 由于 Y Y Y的值不可能为负,故 y &lt; 0 y&lt;0 y<0时, f Y ( y ) = 0 f_Y(y)=0 fY(y)=0。若 y &gt; 0 y&gt;0 y>0,这是对于任意的 y y y,有两个 x x x值与之对应,即
x 1 = y / b , x 2 = − y / b x_{1}=\sqrt{y / b}, \quad x_{2}=-\sqrt{y / b} x1=y/b ,x2=y/b \ 由于 J 1 = d x 1 d y = 1 2 b y , J 2 = d x 2 d y = − 1 2 b y J_{1}=\frac{\mathrm{d} x_{1}}{\mathrm{d} y}=\frac{1}{2 \sqrt{b y}}, J_{2}=\frac{\mathrm{d} x_{2}}{\mathrm{d} y}=-\frac{1}{2 \sqrt{b y}} J1=dydx1=2by 1,J2=dydx2=2by 1,因此
f Y ( y ) = 1 2 b y [ f x ( y / b ) + f X ( − y / b ) ] ( y &gt; 0 ) f_{Y}(y)=\frac{1}{2 \sqrt{b y}}\left[f_{x}(\sqrt{y / b})+f_{X}(-\sqrt{y / b})\right] \quad(y&gt;0) fY(y)=2by 1[fx(y/b )+fX(y/b )](y>0)


(2)基本公式

E X = ∫ − ∞ + ∞ x f ( x ) d x EX=\int_{-\infty}^{+\infty} x f(x) \mathrm{d} x EX=+xf(x)dx
D X = E ( X − E X ) = E X 2 − E 2 X DX=E(X-EX)=EX^2-E^2X DX=E(XEX)=EX2E2X
Cov ⁡ ( X , Y ) = E [ ( X − E X ) ( Y − E Y ) ] = E [ X Y ] − E X E Y \operatorname{Cov}(X, Y)=E[(X-EX)(Y-EY)]=E[X Y]-EXEY Cov(X,Y)=E[(XEX)(YEY)]=E[XY]EXEY
ρ X Y = Cov ⁡ ( X , Y ) D ( X ) D ( Y ) \rho_{X Y}=\frac{\operatorname{Cov}(X, Y)}{\sqrt{D(X)} \sqrt{D(Y)}} ρXY=D(X) D(Y) Cov(X,Y)

2 随机过程

(1)一维概率分布

F x ( x , t ) = P ( X ( t ) ⩽ x } F_{x}(x, t)=P(X(t) \leqslant x\} Fx(x,t)=P(X(t)x}
如果 F X ( x , t ) F_X(x,t) FX(x,t)的一阶导数存在,则定义
f x ( x , t ) = ∂ F X ( x , t ) ∂ x f_{x}(x, t)=\frac{\partial F_{X}(x, t)}{\partial x} fx(x,t)=xFX(x,t)
为随机过程 X ( t ) X(t) X(t)的一维概率密度

(2)二维概率分布

对于任意的时刻 t 1 , t 2 t_1,t_2 t1,t2,以及任意两个实数 x 1 , x 2 x_1,x_2 x1,x2,定义
F x ( x 1 , x 2 , t 1 , t 2 ) = P { X ( t 1 ) ⩽ x 1 , X ( t 2 ) ⩽ x 2 ⟩ F_{x}\left(x_{1}, x_{2}, t_{1}, t_{2}\right)=P\left\{X\left(t_{1}\right) \leqslant x_{1}, X\left(t_{2}\right) \leqslant x_{2}\right\rangle Fx(x1,x2,t1,t2)=P{X(t1)x1,X(t2)x2
为随机过程 X ( t ) X(t) X(t)的二维概率分布,如果 F X ( x 1 , x 2 , t 1 , t 2 ) F_{X}\left(x_{1}, x_{2}, t_{1}, t_{2}\right) FX(x1,x2,t1,t2) x 1 , x 2 x_1,x_2 x1,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 z f_{x}\left(x_{1}, x_{2}, t_{1}, t_{2}\right)=\frac{\partial^{2} F_{x}\left(x_{1}, x_{2}, t_{1}, t_{2}\right)}{\partial x_{1} \partial x_{z}} fx(x1,x2,t1,t2)=x1xz2Fx(x1,x2,t1,t2)
为随机过程 X ( t ) X(t) X(t)的二维概率密度。

(3)随机过程的数字特征
均值

m X ( t ) = E [ X ( t ) ] = ∫ − ∞ + ∞ x f x ( x , t ) d x m_{X}(t)=E[X(t)]=\int_{-\infty}^{+\infty} x f_{x}(x, t) \mathrm{d} x mX(t)=E[X(t)]=+xfx(x,t)dx

方差

σ X 2 ( t ) = E { [ X ( t ) − m X ( t ) ] 2 } = E [ X 2 ( t ) ] − m X 2 ( t ) \sigma_{X}^{2}(t)=E\left\{\left[X(t)-m_{X}(t)\right]^{2}\right\}=E\left[X^{2}(t)\right]-m_{X}^{2}(t) σX2(t)=E{[X(t)mX(t)]2}=E[X2(t)]mX2(t)

自相关函数与协方差函数

对于任意两个时刻 t 1 , t 2 t_1,t_2 t1,t2,定义

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}\left(t_{1}, t_{2}\right)=E\left[X\left(t_{1}\right) X\left(t_{2}\right)\right]=\int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} x_{1} x_{2} f\left(x_{1}, x_{2}, t_{1}, t_{2}\right) \mathrm{d} x_{1} \mathrm{d} x_{2} RX(t1,t2)=E[X(t1)X(t2)]=++x1x2f(x1,x2,t1,t2)dx1dx2
为随机过程 X ( t ) X(t) X(t)的自相关函数。自相关就是函数和函数本身的相关性,当函数中有周期性分量的时候,自相关函数的极大值能够很好的体现这种周期性
相关性的描述除了用相关函数外,有时也用协方差函数。定义
K x ( t 1 , t 2 ) = E { [ X ( t 1 ) − m X ( t 1 ) ] [ X ( t 2 ) − m X ( t 2 ) ] } K_{x}\left(t_{1}, t_{2}\right)=E\left\{\left[X\left(t_{1}\right)-m_{X}\left(t_{1}\right)\right]\left[X\left(t_{2}\right)-m_{X}\left(t_{2}\right)\right]\right\} Kx(t1,t2)=E{[X(t1)mX(t1)][X(t2)mX(t2)]}
为随机过程 X ( t ) X(t) X(t)的协方差函数,协方差函数也可以表示为
K X ( t 1 , t 2 ) = E [ X ( t 1 ) X ( t 2 ) ] − m X ( t 1 ) m X ( t 2 ) = R X ( t 1 , t 2 ) − m X ( t 1 ) m X ( t 2 ) \begin{aligned} K_{X}\left(t_{1}, t_{2}\right) &amp;=E\left[X\left(t_{1}\right) X\left(t_{2}\right)\right]-m_{X}\left(t_{1}\right) m_{X}\left(t_{2}\right) \\ &amp;=R_{X}\left(t_{1}, t_{2}\right)-m_{X}\left(t_{1}\right) m_{X}\left(t_{2}\right) \end{aligned} KX(t1,t2)=E[X(t1)X(t2)]mX(t1)mX(t2)=RX(t1,t2)mX(t1)mX(t2)
相关运算从线性空间的角度看其实是内积运算,而两个向量的内积在线性空间中表示一个向量向另一个向量的投影,表示两个向量的相似程度,所以相关运算就体现了这种相似程度。

(4)广义平稳随机过程

如果随机过程 X ( t ) X(t) X(t)的均值为常数,自相关函数只与$\tau=t_{1}-t_{2}$有关,即
m X ( t ) = m X R X ( t 1 , t 2 ) = R X ( τ ) , τ = t 1 − t 2 \begin{array}{l}{m_{X}(t)=m_{X}} \\ {R_{X}\left(t_{1}, t_{2}\right)=R_{X}(\tau), \quad \tau=t_{1}-t_{2}}\end{array} mX(t)=mXRX(t1,t2)=RX(τ),τ=t1t2
则称随机过程 X ( t ) X(t) X(t)是广义平稳的

广义平稳随机过程的自相关函数
在这里插入图片描述

(5)平稳随机过程的相关系数和相关时间

相关系数

r X ( τ ) = K x ( τ ) σ X 2 = R X ( τ ) − m X 2 σ X 2 r_{X}(\tau)=\frac{K_{x}(\tau)}{\sigma_{X}^{2}}=\frac{R_{X}(\tau)-m_{X}^{2}}{\sigma_{X}^{2}} rX(τ)=σX2Kx(τ)=σX2RX(τ)mX2

相关时间

∣ r X ( τ 0 ) ∣ ⩽ 0.05 \left|r_{X}\left(\tau_{0}\right)\right| \leqslant 0.05 rX(τ0)0.05

(6)各态历经过程

   各态历经随机过程一个样本函数经历了随机过程所有可能的状态,通过对一条样本函数的观测就可以估计出随机过程的均值、方差和相关函数
  对大多数的平稳随机过程而言,都具有各态历经性
  在这里插入图片描述

3 随机过程的联合分布和互相关函数

互相关函数

R X Y ( t 1 , t 2 ) = E [ X ( t 1 ) Y ( t 2 ) ] = ∫ − ∞ + ∞ ∫ − ∞ + ∞ x y f x y ( x , y , t 1 , t 2 ) d x d y R_{X Y}\left(t_{1}, t_{2}\right)=E\left[X\left(t_{1}\right) Y\left(t_{2}\right)\right]=\int_{-\infty}^{+\infty} \int_{-\infty}^{+\infty} x y f_{x y}\left(x, y, t_{1}, t_{2}\right) \mathrm{d} x \mathrm{d} y RXY(t1,t2)=E[X(t1)Y(t2)]=++xyfxy(x,y,t1,t2)dxdy

互协方差函数

K X Y ( t 1 , t 2 ) = E { [ X ( t 1 ) − m x ( t 1 ) ] [ Y ( t 2 ) − m X ( t 2 ) ] } K_{X Y}\left(t_{1}, t_{2}\right)=E\left\{\left[X\left(t_{1}\right)-m_{x}\left(t_{1}\right)\right]\left[Y\left(t_{2}\right)-m_{X}\left(t_{2}\right)\right]\right\} KXY(t1,t2)=E{[X(t1)mx(t1)][Y(t2)mX(t2)]}

互协方差函数与互相关函数的关系

K X Y ( t 1 , t 2 ) = R X Y ( t 1 , t z ) − m X ( t 1 ) m X ( t z ) K_{X Y}\left(t_{1}, t_{2}\right)=R_{X Y}\left(t_{1}, t_{z}\right)-m_{X}\left(t_{1}\right) m_{X}\left(t_{z}\right) KXY(t1,t2)=RXY(t1,tz)mX(t1)mX(tz)

广义联合平稳

如果
m X ( t ) = m x m Y ( t ) = m Y R X Y ( t 1 , t 2 ) = R X Y ( τ ) , τ = t 1 − t 2 \begin{array}{l}{m_{X}(t)=m_{x}} \\ {m_{Y}(t)=m_{Y}} \\ {R_{X Y}\left(t_{1}, t_{2}\right)=R_{X Y}(\tau), \quad \tau=t_{1}-t_{2}}\end{array} mX(t)=mxmY(t)=mYRXY(t1,t2)=RXY(τ),τ=t1t2
则称 X ( t ) X(t) X(t) Y ( t ) Y(t) Y(t)是广义联合平稳的

联合平稳随机过程互相关函数性质

(1)

R X Y ( − τ ) = R Y X ( τ ) K X Y ( − τ ) = K Y X ( τ ) \begin{array}{l}{R_{X Y}(-\tau)=R_{Y X}(\tau)} \\ {K_{X Y}(-\tau)=K_{Y X}(\tau)}\end{array} RXY(τ)=RYX(τ)KXY(τ)=KYX(τ)
因为
R X Y ( − τ ) = E [ X ( t − τ ) Y ( t ) ] = E [ Y ( t ) X ( t − τ ) ] = R Y X ( τ ) R_{X Y}(-\tau)=E[X(t-\tau) Y(t)]=E[Y(t) X(t-\tau)]=R_{Y X}(\tau) RXY(τ)=E[X(tτ)Y(t)]=E[Y(t)X(tτ)]=RYX(τ)

(2)

∣ R X Y ( τ ) ∣ 2 ⩽ R X ( 0 ) R Y ( 0 ) 2 R X Y ( τ ) ⩽ R X ( 0 ) + R Y ( 0 ) ∣ K X Y ( τ ) ∣ 2 ⩽ σ X 2 σ Y 2 \begin{array}{l}{\left|R_{X Y}(\tau)\right|^{2} \leqslant R_{X}(0) R_{Y}(0)} \\ {2 R_{X Y}(\tau) \leqslant R_{X}(0)+R_{Y}(0)} \\ {\left|K_{X Y}(\tau)\right|^{2} \leqslant \sigma_{X}^2 \sigma_{Y}^{2}}\end{array} RXY(τ)2RX(0)RY(0)2RXY(τ)RX(0)+RY(0)KXY(τ)2σX2σY2

(3)

X ( t ) X(t) X(t) Y ( t ) Y(t) Y(t)是联合平稳的,则 Z ( t ) = X ( t ) + Y ( t ) Z(t)=X(t)+Y(t) Z(t)=X(t)+Y(t)是平稳的,且
R Z ( τ ) = R X ( τ ) + R Y ( τ ) + R X Y ( τ ) + R Y X ( τ ) R_{Z}(\tau)=R_{X}(\tau)+R_{Y}(\tau)+R_{X Y}(\tau)+R_{Y X}(\tau) RZ(τ)=RX(τ)+RY(τ)+RXY(τ)+RYX(τ)

4 随机过程功率谱

实质是推广的频谱分析,使用了截尾函数,使得信号满足绝对可积的狄里克雷条件

  • 4
    点赞
  • 13
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值