我是靠谱客的博主 笑点低招牌,最近开发中收集的这篇文章主要介绍随机过程,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

随机过程

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)=Pleft{Xleft(t_{1}right) leqslant x_{1}, Xleft(t_{2}right) leqslant x_{2}rightrangle 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)=Eleft{left[X(t)-m_{X}(t)right]^{2}right}=Eleft[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)=Eleft[Xleft(t_{1}right) Xleft(t_{2}right)right]=int_{-infty}^{+infty} int_{-infty}^{+infty} x_{1} x_{2} fleft(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)=Eleft{left[Xleft(t_{1}right)-m_{X}left(t_{1}right)right]left[Xleft(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;=Eleft[Xleft(t_{1}right) Xleft(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)=Eleft[Xleft(t_{1}right) Yleft(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)=Eleft{left[Xleft(t_{1}right)-m_{x}left(t_{1}right)right]left[Yleft(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 随机过程功率谱

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

最后

以上就是笑点低招牌为你收集整理的随机过程的全部内容,希望文章能够帮你解决随机过程所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(35)

评论列表共有 0 条评论

立即
投稿
返回
顶部