现代通信原理4.2:随机过程

  前面我们简单回顾了随机变量。这一部分我们来讨论随机过程,大家注意理解二者之间的区别与关系。

1、什么是随机过程?

  所谓随机过程,可以看成是随某参量(一般为时间t)变化的随机变量。举例来说,假定下图中所示随机噪声源,其输出电压为随机变量,且该变量岁时间变化。我们对其输出电压随时间变化的函数进行测量并记录,其中第 j j j次测试得到的结果用 X ( A j , t ) X(A_j,t) X(Aj,t)来表示,简写为 X j ( t ) X_j(t) Xj(t),我们称之为第 j j j个样本函数。显然,每个样本函数都是确定函数。如果一共进行了 N N N次测试,就可以得到 N N N个样本函数,如下图所示。所有可能的 N N N个样本函数,就构成样本函数空间 { X 1 ( t ) , X 2 ( t ) , … , X N ( t ) } \{X_1(t),X_2(t),\ldots,X_N(t)\} { X1(t),X2(t),,XN(t)}
在这里插入图片描述
   我们可以从两个不同的角度来理解随机过程。
  第一个角度,我们可以把随机过程看成样本函数的集合。从上面的例子来看,所有的 N N N个样本函数,构成样本函数空间。而每次测试得到的噪声源输出电压,一定是空间中的某个样本函数,但具体是哪个函数,则是随机的。

  注意这里的 N N N可以是有限值,也可以是无限大的。换句话说,样本空间中样本函数个数,可以是有限多个,也可以是无线多个。
  举例来说,如果我们要发送二进制数据的"0"或者“1”,我们可以用两种波形来表示它们,例如,用持续时间为1ms的+5V电平表示"1",持续时间为1ms的-5V电平表示"0"。因此我们的样本函数集中包含下图所示两个样本函数,即 N = 2 N=2 N=2,从接收机来看,会随机地收到两个函数中的其中一个,并对此进行检测。因此这是一个二进制的数字通信系统。
在这里插入图片描述
  考虑另外一种情况。如果我们用麦克风采集语音信号并发送给接收端。从接收端来看,可能收到的信号波形有无穷多种可能。因此样本函数集中的样本函数有无穷多种。这也就是我们前面介绍过的模拟通信系统。

一般来说,如果 N N N为无限值,即状态为连续的,我们称 X ( t ) X(t) X(t)为连续随机过程;如果 N N N为有限值,即状态为离散的,我们称 X ( t ) X(t) X(t)为离散随机过程。
  第二个角度,可以把随机过程看成一系列随机变量的集合。同样是上面的随机噪声源的例子,如下图所示,如果我们观察每次实验 t 1 t_1 t1时刻对应的电压值,会发现在第 j j j个样本函数上的值为 X ( A j , t 1 ) X(A_j,t_1) X(Aj,t1),而 N N N次不同的实验其取值都不相同。因此,噪声源在 t 1 t_1 t1时刻的输出 X ( t 1 ) X(t_1) X(t1)为一个随机变量。显然 t 2 t_2 t2时刻的输出 X ( t 2 ) X(t_2) X(t2)也为随机变量。就这样,在时间轴上的一系列随机变量 { X ( t 1 ) , X ( t 2 ) , … } \{X(t_1),X(t_2),\ldots\} { X(t1),X(t2),}构成了随机过程 X ( t ) X(t) X(t)
在这里插入图片描述
同样,随机过程所包含的随机变量的个数可以是有限个,也可以是无限多。如果是无限多,则在时间上是连续的;如果是有限多,则在时间上是离散的。时间上离散的随机过程,我们也常称之为随机序列。

2、平稳性与遍历性

  随机过程有两个重要的性质,平稳性和遍历性(也称各态历经性)。如果随机过程具有这两个性质,其分析过程就可以简化。

2.1 平稳性

  我们先来看 N N N阶平稳和严格平稳的定义。

如果随机过程 X ( t ) X(t) X(t)对于任意的时刻 t 1 , t 2 , … , t N t_1,t_2,\ldots,t_N t1,t2,,tN,都有其 N N N维联合概率密度函数
p X ( x 1 , x 2 , … , x N ; t 1 , t 2 , … , t N ) = p X ( x 1 , x 2 , … , x N ; t 1 + t 0 , t 2 + t 0 , … , t N + t 0 ) p_X(x_1,x_2,\ldots,x_N;t_1,t_2,\ldots,t_N)=p_X(x_1,x_2,\ldots,x_N;t_1+t_0,t_2+t_0,\ldots,t_N+t_0) pX(x1,x2,,xN;t1,t2,,tN)=pX(x1,x2,,xN;t1+t0,t2+t0,,tN+t0)其中 t 0 t_0 t0为任意时刻,则称 X ( t ) X(t) X(t) N N N阶平稳。此外,如果 X ( t ) X(t) X(t)对于 N → ∞ N\rightarrow \infty N阶平稳,则称 X ( t ) X(t) X(t)为严格平稳。

从上面定义中我们可以看出,平稳性这个词很准确表述了概率密度函数不随时间变化的意思,即表现出来时间上的不变性。因为我们在求统计平均的时候都会用到概率密度函数,因此概率密度函数的平稳也就决定了均值、方差以及高阶矩等统计特性的平稳性。

2.2 遍历性

  如果一个随机过程的集合平均,等于任意样本函数的时间平均,则我们称该随机过程为遍历的(或各态历经的),并且

  • 如果一个随机过程是遍历的,则所有的时间平均和集合平均可以互换。
  • 遍历的随机过程一定是平稳的。
      如果说平稳性把对随机过程的分析简化成对于随机变量的分析,遍历性进一步把对随机变量的分析简化为对确定函数的时间分析。遍历性这个词也很准确描述了这样一个性质,即在时间维度上,任何一个样本函数都能够把该随机过程所有的随机状态历经到。特别注意的是,遍历的随机过程一定是平稳的。

  作为例子,我们先来看随机随机过程 X ( t ) X(t) X(t)的均值,即一阶原点矩:
m X ( t ) = E [ X ( t ) ] = ∫ − ∞ ∞ x p X ( x ; t ) d x m_X(t)={\rm E}[X(t)]=\int_{-\infty}^{\infty}xp_X(x;t)dx mX(t)=E[X(t)]=xpX(x;t)dx为关于时间的确定函数。如果 X ( t ) X(t) X(t)是一阶平稳的,则一阶概率密度函数不随时间 t t t变化,则其均值 E [ X ( t ) ] = m X {\rm E}[X(t)]=m_X E[X(t)]=mX也不随时间 t t t变化,我们称 X ( t ) X(t) X(t)为均值平稳。再进一步,若 X ( t ) X(t) X(t)为遍历的, X j ( t ) = X ( A j , t ) X_j(t)=X(A_j,t) Xj(t)=X(Aj,t)为其任一样本函数,则有
m X ( t ) = E [ X ( t ) ] = ∫ − ∞ ∞ x p X ( x ; t ) d x = X j ( t ) ‾ = lim ⁡ T → ∞ 1 T ∫ − T 2 T 2 X j ( t ) d t , \begin{aligned} m_X(t)&={\rm E}[X(t)]=\int_{-\infty}^{\infty}xp_X(x;t)dx\\ &=\overline{X_j(t)}=\lim_{T\rightarrow \infty}\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}X_j(t)dt, \end{aligned} mX(t)=E[X(t)]=xpX(x;t)dx=Xj(t)=TlimT12T2TXj(t)dt,即随机过程的统计平均等于样本函数的时间平均,注意这里的样本函数是确定函数。

3、自相关函数与广义平稳

3.1 自相关函数与自协方差函数

  对任意的两个时刻 t 1 , t 2 t_1,t_2 t1,t2,实随机过程 X ( t ) X(t) X(t)的自相关函数(auto-correlaton function, ACF)定义为
R X ( t 1 , t 2 ) = E [ X ( t 1 ) X ( t 2 ) ] = ∫ − ∞ ∞ x 1 x 2 p X ( x 1 , x 2 ; t 1 , t 2 ) d x 1 d x 2 , R_X(t_1,t_2)={\rm E}\left[X(t_1)X(t_2) \right]=\int_{-\infty}^{\infty}x_1x_2p_X(x_1,x_2;t_1,t_2)dx_1dx_2, RX(t1,t2)=E[X(t1)X(t2)]=x1x2pX(x1,x2;t1,t2)dx1dx2,它定义了随机变量 X ( t 1 ) X(t_1) X(t1) X ( t 2 ) X(t_2) X(t2)之间的二阶混合原点矩或称相关矩,而实随机过程 X ( t ) X(t) X(t)的自协方差函数(二阶混合中心距)定义为
C X ( t 1 , t 2 ) = E { [ X ( t 1 ) − m X ( t 1 ) ] [ X ( t 2 ) − m X ( t 2 ) ] } = ∫ − ∞ ∞ [ x 1 − m X ( t 1 ) ] [ x 2 − m X ( t 2 ) ] p X ( x 1 , x 2 ; t 1 , t 2 ) d x 1 d x 2 . \begin{aligned} C_X(t_1,t_2)&={\rm E}\left\{[X(t_1)-m_X(t_1)][X(t_2)-m_X(t_2)] \right\}\\ &=\int_{-\infty}^{\infty}[x_1-m_X(t_1)][x_2-m_X(t_2)]p_X(x_1,x_2;t_1,t_2)dx_1dx_2. \end{aligned} CX(t1,t2)=E{ [X(t1)mX(t1)][X(t2)mX(t2)]}=[x1mX(t1)][x2mX(t2)]pX(x1,x2;t1,t2)dx1dx2.显然,有
C X ( t 1 , t 2 ) = R X ( t 1 , t 2 ) − m X ( t 1 ) m X ( t 2 ) C_X(t_1,t_2)=R_X(t_1,t_2)-m_X(t_1)m_X(t_2) CX(t1,t2)=RX(t1,t2)mX(t1)mX(t2)成立。如果对于任意的 t 1 t_1 t1 t 2 t_2 t2都有 C X ( t 1 , t 2 ) = 0 C_X(t_1,t_2)=0 CX(t1,t2)=0,则称随机过程任意两个时刻是不相关的。

3.2 广义平稳

  如果随机过程 X ( t ) X(t) X(t)满足

  1. E [ X ( t ) ] = m X ( t ) = {\rm E}[X(t)]=m_X(t)= E[X(t)]=mX(t)=常数;
  2. R X ( t 1 , t 2 ) = R X ( τ ) ,   τ = t 1 − t 2 R_X(t_1,t_2)=R_X(\tau),\ \tau=t_1-t_2 RX(t1,t2)=RX(τ), τ=t1t2;

则称X(t)为广义平稳(wide-sense stationary,W.S.S.)。上面第一个条件意味着 X ( t ) X(t) X(t)是均值平稳的,第二个条件意味着 X ( t ) X(t) X(t)的自相关函数与 t 1 t_1 t1 t 2 t_2 t2的绝对时刻无关,只与两个时刻之间的间隔大小 τ \tau τ有关系。
  广义平稳随机过程的自相关函数有如下几个性质

  1. R X ( τ ) = R X ( − τ ) R_X(\tau)=R_X(-\tau) RX(τ)=RX(τ);
  2. R X ( 0 ) ≥ R X ( τ ) R_X(0)\ge R_X(\tau) RX(0)RX(τ);
  3. R X ( 0 ) = E [ X 2 ( t ) ] R_X(0)={\rm E}[X^2(t)] RX(0)=E[X2(t)].

性质3表明, R X ( 0 ) R_X(0) RX(0)等于 X ( t ) X(t) X(t)的均方值。

4、功率谱密度与Wiener-Khintchine 定理

4.1 随机过程功率谱密度函数定义

  随机过程 X ( t ) X(t) X(t)的功率谱密度函数定义为
P X ( f ) = lim ⁡ T → ∞ { E [ ∣ X T ( f ) ∣ ] T } , P_X(f)=\lim_{T\rightarrow\infty}\left\{\frac{ {\rm E}{\Large [}|X_T(f)|{\Large ]}}{T}\right\}, PX(f)=Tlim{ TE[XT(f)]},其中
X T ( f ) = ∫ − ∞ ∞ x T ( t ) e − j 2 π f t d t = ∫ − T 2 T 2 x ( t ) e − j 2 π f t d t , X_T(f)=\int_{-\infty}^{\infty}x_T(t)e^{-j2\pi ft}dt=\int_{-\frac{T}{2}}^{\frac{T}{2}}x(t)e^{-j2\pi ft}dt, XT(f)=xT(t)ej2πftdt=2T2Tx(t)ej2πftdt, x T ( t ) x_T(t) xT(t)为随机过程 X ( t ) X(t) X(t) [ − T 2 , T 2 ] [-\frac{T}{2},\frac{T}{2}] [2T,2T]上的截断函数。

4.2 Wiener-Khintchine 定理

  如果 X ( t ) X(t) X(t)为广义平稳随机过程,则其功率谱密度函数为自相关函数的傅立叶变换,即
R X ( τ ) ↔   F . T .   P X ( f ) . R_X(\tau)\xleftrightarrow{\ \mathcal{F.T.}\ } P_X(f). RX(τ) F.T.  PX(f).

4.3 随机过程功率谱密度函数性质

  随机过程 X ( t ) X(t) X(t)的功率谱密度函数 P X ( t ) P_X(t) PX(t)具有如下性质

  • P X ( f ) P_X(f) PX(f)总为实数;
  • P X ( f ) ≥ 0 P_X(f)\ge 0 PX(f)0
  • 如果 X ( t ) X(t) X(t)为实随机过程,则 P X ( f ) = P X ( − f ) P_X(f)=P_X(-f) PX(f)=PX(f)
  • 随机信号 X ( t ) X(t) X(t)的平均功率 P = ∫ − ∞ ∞ P X ( f ) d f P=\int_{-\infty}^{\infty}P_X(f)df P=PX(f)df,如果 X ( t ) X(t) X(t)为W.S.S.,则 P = ∫ − ∞ ∞ P X ( f ) d f = R X ( 0 ) = E [ X 2 ( t ) ] P=\int_{-\infty}^{\infty}P_X(f)df=R_X(0)={\rm E}[X^2(t)] P=PX(f)df=RX(0)=E[X2(t)]

【证明】根据Wiener-Khintchine 定理,有
R X ( τ ) = ∫ − ∞ ∞ P X ( f ) e j 2 π f τ d f , R_X(\tau)=\int_{-\infty}^{\infty}P_X(f)e^{j2\pi f\tau}df, RX(τ)=PX(f)ej2πfτdf, τ = 0 \tau=0 τ=0,有
R X ( 0 ) = ∫ − ∞ ∞ P X ( f ) d f = P . R_X(0)=\int_{-\infty}^{\infty}P_X(f)df=P. RX(0)=PX(f)df=P.

  • P X ( 0 ) = ∫ − ∞ ∞ R X ( τ ) d τ P_X(0)=\int_{-\infty}^{\infty}R_X(\tau)d\tau PX(0)=RX(τ)dτ.

5、遍历随机信号的物理参数

  若随机信号 X ( t ) X(t) X(t)为遍历的,显然样本函数的时间平均与其统计平均是等效的。因此,《现代通信原理2.2:信号时间平均算子与信号物理参数》中对确定信号物理参数的定义,就可以推广到随机信号。用 ξ ( t ) \xi(t) ξ(t)表示 X ( t ) X(t) X(t)的样本函数,则有

  • 直流
    X D C = ξ ( t ) ‾ = E [ X ( t ) ] = m X X_{\rm DC}=\overline{\xi(t)}={\rm E}[X(t)]=m_X XDC=ξ(t)=E[X(t)]=mX
  • 归一化平均功率
    P = △ ξ 2 ( t ) ‾ = E [ X 2 ( t ) ] = R X ( 0 ) = R ξ ( 0 ) P\overset{\triangle}{=}\overline{\xi^2 (t)}={\rm E}[X^2(t)]=R_X(0)=R_{\xi}(0) P=ξ2(t)=E[X2(t)]=RX(0)=Rξ(0)这里的 R X ( τ ) R_X(\tau) RX(τ)为随机过程 X ( t ) X(t) X(t)的自相关函数, R ξ ( τ ) R_{\xi}(\tau) Rξ(τ)为样本函数 ξ ( t ) \xi (t) ξ(t)的自相关函数。
  • 归一化直流功率
    P D C = △ ξ 2 ( t ) ‾ = { E [ X ( t ) ] } 2 = m X 2 P_{\rm DC}\overset{\triangle}{=}\overline{\xi^2 (t)}=\{ {\rm E}[X(t)]\}^2=m_X^2 PDC=ξ2(t)={ E[X(t)]}2=mX2
  • 归一化交流功率
    P A C = △ [ ξ ( t ) − X D C ] 2 ‾ = E [ X ( t ) − m X ] 2 = P − P D C = σ X 2 P_{\rm AC}\overset{\triangle}{=}\overline{[\xi (t)-X_{\rm DC}]^2}={\rm E}[X(t)-m_X]^2=P-P_{\rm DC}=\sigma_X^2 PAC=[ξ(t)XDC]2=E[X(t)mX]2=PPDC=σX2
  • 均方根值(RMS)
    X r m s = △ ξ 2 ( t ) ‾ = E [ X 2 ( t ) ] = R X ( 0 ) = ∫ − ∞ ∞ P X ( f ) d f X_{\rm rms}\overset{\triangle}{=}\sqrt{\overline{\xi^2 (t)}}=\sqrt{ {\rm E}[X^2(t)]}=\sqrt{R_X(0)}=\sqrt{\int_{-\infty}^{\infty}P_X(f)df} Xrms=ξ2(t) =E[X2(t)] =RX(0) =PX(f)df

猜你喜欢

转载自blog.csdn.net/tanghonghanhaoli/article/details/100185439