通信原理——随机过程

通信原理——随机过程通信原理学习 通信原理随机过程

大家好,欢迎来到IT知识分享网。

1、随机过程的基本概念

  • 随机事件(简称事件、实现样本、记录):在随机试验中,可能出现也可能不出现,而在大量重复试验中具有某种规律性的事件
  • 随机变量:定义在样本空间上的实值函数,是不确定的
  • 随机过程:无数个随机变量构成的总体
  • 离散随机变量:随机变量x的取值个数是有限的或可数无穷个
  • 连续随机变量:随机变量x可能的取值充满某一有限或无限区间

1.1 随机过程的分布函数

ξ ( t ) \xi(t) ξ(t)表示一个随机过程,则在任意一个时刻 t 1 t_1 t1 ξ ( t 1 ) \xi(t_1) ξ(t1)是一个随机变量,定义随机过程 ξ ( t ) \xi(t) ξ(t)

一维概率分布函数: F 1 ( x 1 ; t 1 ) = P { ξ ( t 1 ) ⩽ x 1 } F_1(x_1;t_1) = P\left\{ \begin{aligned} \xi(t_1) \leqslant x_1 \end{aligned} \right\} F1(x1;t1)=P{
ξ(t1)x1}

一维概率密度函数: f 1 ( x 1 ; t 1 ) = ∂ F 1 ( x 1 ; t 1 ) ∂ x 1 f_1(x_1;t_1) = \frac {\partial F_1(x_1;t_1)}{\partial x_1} f1(x1;t1)=x1F1(x1;t1) F 1 ( x 1 ; t 1 ) = ∫ − ∞ x f 1 ( x 1 ; t 1 ) d x F_1(x_1;t_1) = \int_{-\infty}^{x}f_1(x_1;t_1)dx F1(x1;t1)=xf1(x1;t1)dx

n维概率分布函数: F n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) = P { ξ ( t 1 ) ≤ x 1 , ξ ( t 2 ) ≤ x 2 , ⋯   , ξ ( t n ) ≤ x n } Fn(x_1,x_2, \cdots ,x_n;t_1,t_2, \cdots ,t_n) = P\left\{ {\xi \left( {t_1} \right) \le x_1,\xi \left( {t_2} \right) \le x_2, \cdots ,\xi \left( {t_n} \right) \le x_n} \right\} Fn(x1,x2,,xn;t1,t2,,tn)=P{
ξ(t1)x1,ξ(t2)x2,,ξ(tn)xn}

n维概率密度函数: f n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) = ∂ F n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) ∂ x 1 ∂ x 2 ⋯ ∂ x n f_n(x_1,x_2, \cdots ,x_n;t_1,t_2, \cdots ,t_n) = \frac{
{\partial F_n(x_1,x_2, \cdots ,x_n;t_1,t_2, \cdots ,t_n)}}{
{\partial x_1\partial x_2 \cdots \partial x_n}}
fn(x1,x2,,xn;t1,t2,,tn)=x1x2xnFn(x1,x2,,xn;t1,t2,,tn)
F n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) = ∫ − ∞ x 1 ∫ − ∞ x 2 ⋯ ∫ − ∞ x n f n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) d x 1 d x 2 ⋯ d x n F_n(x_1,x_2, \cdots ,x_n;t_1,t_2, \cdots ,t_n) = \int\limits_{ – \infty }^{x_1} {\int\limits_{ – \infty }^{x_2} { \cdots \int\limits_{ – \infty }^{x_n} {f_n(x_1,x_2, \cdots ,x_n;t_1,t_2, \cdots ,t_n)} } } dx_1dx_2 \cdots dx_n Fn(x1,x2,,xn;t1,t2,,tn)=x1x2xnfn(x1,x2,,xn;t1,t2,,tn)dx1dx2dxn

1.2 随机过程的数字特征

  1. 均值(数学期望)

随机过程 ξ ( t ) \xi(t) ξ(t)的数学期望:

  • ξ ( t ) \xi(t) ξ(t)连续,定义 E [ ξ ( t ) ] = ∫ − ∞ ∞ x f 1 ( x , t ) d x = a ( t ) E\left[ {\xi \left( t \right)} \right] = \int_{ – \infty }^\infty x {f_1}\left( {x,t} \right)dx = a\left( t \right) E[ξ(t)]=xf1(x,t)dx=a(t)
  • ξ ( t ) \xi(t) ξ(t)离散,定义 E [ ξ i ( t ) ] = ∑ i = 1 K ξ i ( t ) P ( ξ i ( t ) ) = a ( t ) E\left[ {
    {\xi _{\rm{i}}}\left( t \right)} \right] = \sum\nolimits_{
    {\rm{i = }}1}^K {
    {\xi _{\rm{i}}}\left( t \right)P\left( {
    {\xi _{\rm{i}}}\left( t \right)} \right)} = a(t)
    E[ξi(t)]=i=1Kξi(t)P(ξi(t))=a(t)

ξ ( t ) \xi(t) ξ(t)的均值 E [ ξ ( t ) ] E[\xi(t)] E[ξ(t)]是时间的确定函数,常记为a(t),他表示随机过程的n个样本函数曲线的摆动中心。
2. 方差

随机过程 ξ ( t ) \xi(t) ξ(t)的方差:
D [ ξ ( t ) ] = E { [ ξ ( t ) − a ( t ) ] 2 } = E [ ξ 2 ( t ) ] − a 2 ( t ) = σ 2 ( t ) D\left[ {\xi \left( t \right)} \right] = E\left\{ {
{
{\left[ {\xi \left( t \right) – a\left( t \right)} \right]}^2}} \right\} = E[\xi^2(t)] – a^2(t) = {\sigma ^2}\left( t \right)
D[ξ(t)]=E{
[ξ(t)a(t)]2}
=
E[ξ2(t)]a2(t)=σ2(t)

方差等于均方值与均方平方之差,表示随机过程在时刻t相对于均值a(t)的偏离程度。

  1. 相关函数

ξ ( t ) \xi(t) ξ(t)的协方差函数:
B [ t 1 , t 2 ] = E { [ ξ ( t 1 ) − a ( t 1 ) ] [ ξ ( t 2 ) − a ( t 2 ) ] } = ∫ − ∞ ∞ ∫ − ∞ ∞ [ x 1 − a ( t ) ] [ x 2 − a ( t 2 ) ] f 2 ( x 1 , x 2 ; t 1 , t 2 ) d x 1 d x 2 B\left[ {
{t_1},{t_2}} \right] = E\left\{ {\left[ {\xi \left( {
{t_1}} \right) – a\left( {
{t_1}} \right)} \right]\left[ {\xi \left( {
{t_2}} \right) – a\left( {
{t_2}} \right)} \right]} \right\} \\= \int_{-\infty}^{\infty}\int_{-\infty}^{\infty}[x_1 – a(t)][x_2 – a(t_2)]f_2(x_1,x_2;t_1,t_2)dx_1dx_2
B[t1,t2]=E{
[ξ(t1)a(t1)][ξ(t2)a(t2)]}
=[x1a(t)][x2a(t2)]f2(x1,x2;t1,t2)dx1dx2

式中: a ( t 1 ) a(t_1) a(t1) a ( t 2 ) a(t_2) a(t2)分别是在 t 1 t_1 t1 t 2 t_2 t2时刻得到的 ξ ( t ) \xi(t) ξ(t)的均值; f 2 ( x 1 , x 2 ; t 1 , t 2 ) f_2(x_1,x_2;t_1,t_2) f2(x1,x2;t1,t2) ξ ( t ) \xi(t) ξ(t)的二维概率密度函数。

ξ ( t ) \xi(t) ξ(t)的相关函数:
R [ t 1 , t 2 ] = E [ ξ ( t 1 ) ξ ( t 2 ) ] = ∫ − ∞ ∞ ∫ − ∞ ∞ x 1 x 2 f 2 ( x 1 , x 2 ; t 1 , t 2 ) d x 1 d x 2 R[t1,t2] = E[\xi \left( {t1} \right)\xi \left( {t2} \right)] = \int_{-\infty}^{\infty}\int_{-\infty}^{\infty}x_1x_2f_2(x_1,x_2;t_1,t_2)dx_1dx_2 R[t1,t2]=E[ξ(t1)ξ(t2)]=x1x2f2(x1,x2;t1,t2)dx1dx2
式中: ξ ( t 1 ) \xi(t_1) ξ(t1) ξ ( t 2 ) \xi(t_2) ξ(t2)分别是在 t 1 t_1 t1 t 2 t_2 t2时刻观测 ξ ( t ) \xi(t) ξ(t)得到的随机变量
协方差函数和相关函数之间有着如下确定的关系:
B [ t 1 , t 2 ] = R [ t 1 , t 2 ] − E [ ξ ( t 1 ) ] E [ ξ ( t 2 ) ] B[t1,t2] = R[t1,t2] – E[\xi \left( {t1} \right)]E[\xi \left( {t2} \right)] B[t1,t2]=R[t1,t2]E[ξ(t1)]E[ξ(t2)]

ξ ( t ) , η ( t ) \xi(t),\eta(t) ξ(t),η(t)
互协方差函数: B ξ η [ t 1 , t 2 ] = E [ ξ ( t 1 ) − a ξ ( t 1 ) ] [ η ( t 2 ) − a η ( t 2 ) ] B_{\xi \eta} [t1,t2] = E[\xi \left( {t1} \right) – a\xi (t1)][\eta \left( {t2} \right) – a\eta (t2)] Bξη[t1,t2]=E[ξ(t1)aξ(t1)][η(t2)aη(t2)]

互相关函数: R ξ η [ t 1 , t 2 ] = E [ ξ ( t 1 ) η ( t 2 ) ] R_{\xi \eta} [t1,t2] = E[\xi \left( {t1} \right)\eta \left( {t2} \right)] Rξη[t1,t2]=E[ξ(t1)η(t2)]

互相关系数: ρ = E [ ξ ( t 1 ) − a ξ ( t 1 ) ] [ η ( t 2 ) − a η ( t 2 ) ] E { [ ξ ( t 1 ) − a ξ ( t 1 ) ] 2 } E { [ η ( t 2 ) − a η ( t 2 ) ] 2 } = B ξ η [ t 1 , t 2 ] σ ξ ( t 1 ) σ η ( t 2 ) \rho = \frac{
{E[\xi \left( {t_1} \right) – {a_\xi }(t1)][\eta \left( {t_2} \right) – a_\eta (t_2)]}}{
{\sqrt {E\{ {
{[\xi \left( {t_1} \right) – {a_\xi }(t_1)]}^2}\} E\{ {
{[\eta \left( {t_2} \right) – a_\eta (t_2)]}^2}\} } }} = \frac{
{B_{\xi \eta} [t_1,t_2]}}{
{\sigma_\xi (t_1)\sigma_\eta (t_2)}}
ρ=E{
[ξ(t1)aξ(t1)]2}E{
[η(t2)aη(t2)]2}

E[ξ(t1)aξ(t1)][η(t2)aη(t2)]
=
σξ(t1)ση(t2)Bξη[t1,t2]

(1) ∣ ρ ∣ ⩽ 1 |\rho|\leqslant1 ρ1
(2)相关性:若相关系数 ρ = 0 \rho = 0 ρ=0,则 ξ ( t ) , η ( t ) \xi(t),\eta(t) ξ(t),η(t)是线性不相关的
(3)独立与相关性:若 ξ ( t ) , η ( t ) \xi(t),\eta(t) ξ(t),η(t)是独立的,则线性不相关;反之不然。
结论:一般随机过程的统计特性,原则上都与时刻 t 1 、 t 2 ( t 2 = t 1 + t a u ) t_1 、t_2 (t_2 =t_1+tau) t1t2(t2=t1+tau)有关或者说与时间起点 t 1 t_1 t1及时间间隔 τ \tau τ有关。

2、平稳随机过程

2.1 定义

若一个随机过程 ξ ( t ) \xi(t) ξ(t)的统计特性与时间起点无关,即时间平移不影响其任何统计特性,则称该随机过程是在严格意义下的平稳随机过程,简称严平稳随机过程(狭义平稳随机过程)。
因此,平稳随机过程 ξ ( t ) \xi(t) ξ(t)的任意有限概率密度 函数与时间起点有关,也就是说,对于任意的正整数n和所有实数,有
f n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) = f n ( x 1 , x 2 , ⋯   , x n ; t 1 + τ , t 2 + τ , ⋯   , t n + τ ) f_n(x_1,x_2, \cdots ,x_n;t_1,t_2, \cdots ,t_n)\\ = f_n(x_1,x_2, \cdots ,x_n;t_1 + \tau ,t_2 + \tau , \cdots ,t_n + \tau ) fn(x1,x2,,xn;t1,t2,,tn)=fn(x1,x2,,xn;t1+τ,t2+τ,,tn+τ)
他的一维概率密度函数与时间t无关,即
f 1 ( x 1 , t 1 ) = f 1 ( x 1 ) f_1(x_1,t_1) = f_1(x_1) f1(x1,t1)=f1(x1)
而二维分布函数只与时间间隔有关,即
f 2 ( x 1 , x 2 ; t 1 , t 2 ) = f 2 ( x 1 , x 2 ; τ ) f_2(x_1,x_2;t_1,t_2) = f_2(x_1,x_2;\tau) f2(x1,x2;t1,t2)=f2(x1,x2;τ)
数学特征:
a ( t ) = a σ 2 ( t ) = σ 2 R ( t 1 , t 1 + τ ) = R ( τ ) a(t) = a\\{\sigma ^2}(t) = {\sigma ^2}\\R(t_1,t_1 + \tau ) = R(\tau ) a(t)=aσ2(t)=σ2R(t1,t1+τ)=R(τ)
可见,平稳随机过程 ξ ( t ) \xi(t) ξ(t)具有简明的数学特征:①均值与t无关,为常数a②自相关函数只与时间间隔 τ = t 2 − t 1 \tau = t_2 – t_1 τ=t2t1有关,即 R ( t 1 , t 1 + τ ) = R ( τ ) R(t_1,t_1+\tau) = R(\tau) R(t1,t1+τ)=R(τ),同时满足①和②的过程定义为广义平稳随机过程。

2.2 各态历经性

假设x(t)是平稳过程 ξ ( t ) \xi(t) ξ(t)的任意一次实现(样本),由于它是时间的确定函数,可以求得它的时间平均值。其时间均值和时间相关函数分别定义为:
{ a ‾ = x ( t ) ‾ = lim ⁡ T → ∞ 1 T ∫ − T / 2 T / 2 x ( t ) d t R ( τ ) ‾ = x ( t ) x ( t + τ ) ‾ = lim ⁡ T → ∞ 1 T ∫ − T / 2 T / 2 x ( t ) x ( t + τ ) d t \left\{\begin{matrix} \overline{a} = \overline{x(t)} = \lim_{T\rightarrow\infty} \frac 1 T \int_{-T/2}^{T/2}x(t)dt\\ \overline{R(\tau)} = \overline{x(t)x(t+\tau)} = \lim_{T\rightarrow\infty}\frac 1T\int_{-T/2}^{T/2}x(t)x(t+\tau)dt \end{matrix} \right. {
a=x(t)=limTT1T/2T/2x(t)dtR(τ)=x(t)x(t+τ)=limTT1T/2T/2x(t)x(t+τ)dt

如果平稳过程使
{ a = a ‾ R ( τ ) = R ( τ ) ‾ \left\{\begin{matrix} a = \overline{a}\\ R(\tau) = \overline{R(\tau)} \end{matrix} \right. {
a=aR(τ)=R(τ)

成立。也就是说,平稳过程的统计平均值等于它的任一次实现的时间平均值,则称该平稳过程具有各态历经性。
具有各态历经性的随机过程一定使平稳过程,反之不成立。

2.3 平稳过程的自相关函数

ξ ( t ) \xi(t) ξ(t)为实平稳过程,则它的自相关函数
R ( τ ) = E [ ξ ( t ) ξ ( t + τ ) ] R(\tau) = E[\xi(t)\xi(t+\tau)] R(τ)=E[ξ(t)ξ(t+τ)]
具有如下主要性质:
(1) R ( 0 ) = R [ ξ 2 ( t ) ] R(0) = R[\xi^2(t)] R(0)=R[ξ2(t)],表示 ξ ( t ) \xi(t) ξ(t)的平均功率。
(2) R ( τ ) = R ( − τ ) R(\tau) = R(-\tau) R(τ)=R(τ),表示 τ \tau τ的偶函数。
(3) ∣ R ( τ ) ∣ ⩽ R ( 0 ) |R(\tau)|\leqslant R(0) R(τ)R(0),表示 R ( τ ) R(\tau) R(τ)的上界。
(4) R ( ∞ ) = lim ⁡ τ → ∞ R ( τ ) = E [ ξ ( t ) ] E [ ξ ( t + τ ) ] = E 2 [ ξ ( t ) ] = a 2 R\left( \infty \right) = \mathop {\lim }\limits_{\tau \to \infty } R\left( \tau \right) = E\left[ {\xi \left( t \right)} \right]E\left[ {\xi \left( {t + \tau } \right)} \right] = {E^2}\left[ {\xi \left( t \right)} \right] = {a^2} R()=τlimR(τ)=E[ξ(t)]E[ξ(t+τ)]=E2[ξ(t)]=a2,表示 ξ ( t ) \xi(t) ξ(t)的直流功率。
(5) R ( 0 ) − R ( ∞ ) = σ 2 R\left( 0 \right) – R\left( \infty \right) = {\sigma ^2} R(0)R()=σ2,表示平稳过程的交流功率。

2.4 平稳过程的功率谱密度

对于任意的确定功率信号 f ( t ) f(t) f(t),它的功率谱密度定义为:
P f ( f ) = lim ⁡ T → ∞ ∣ F T ( f ) ∣ 2 T P_f(f) = \lim_{T\rightarrow\infty} \frac {|F_T(f)|^2}T Pf(f)=TlimTFT(f)2
不妨把 f ( t ) f(t) f(t)看成是平稳过程 ξ ( t ) \xi(t) ξ(t)的任一样本,一般而言,不同样本的谱密度不同,因此,过程的功率谱密度应看成是对所有样本的功率谱的统计平均,即
P ξ ( f ) = E [ P f ( f ) ] = lim ⁡ T → ∞ E [ ∣ F T ( f ) ∣ 2 ] T P_{\xi}(f) = E[P_f(f)] = \lim_{T\rightarrow\infty} \frac {E[|F_T(f)|^2]}T Pξ(f)=E[Pf(f)]=TlimTE[FT(f)2]
维纳-辛钦定理:
平稳过程的功率谱密度 P ξ ( f ) P_{\xi}(f) Pξ(f)与其自相关函数 R ( τ ) R(\tau) R(τ)也是一对傅里叶变换关系,即
P s ( ω ) = lim ⁡ T → ∞ E [ ∣ F T ( ω ) ∣ 2 ] T = ∫ − ∞ ∞ R ( τ ) e − j ω τ d τ R ( τ ) = 1 2 π ∫ − ∞ ∞ P s ( ω ) e − j ω τ d ω S = R ( 0 ) = 1 2 π ∫ − ∞ ∞ P s ( ω ) d ω P s ( ω ) ⇔ R ( τ ) {P_s}\left( \omega \right) = \mathop {\lim }\limits_{T \to \infty } \frac{
{E\left[ {
{
{\left| {
{F_T}\left( \omega \right)} \right|}^2}} \right]}}{
{\rm{T}}}{\rm{ = }}\int_{ – \infty }^\infty {R\left( \tau \right){e^{ – j\omega \tau }}d\tau }\\R\left( \tau \right) = \frac{1}{
{2\pi }}\int_{ – \infty }^\infty {
{P_s}\left( \omega \right){e^{ – j\omega \tau }}d\omega } \\S = R\left( 0 \right) = \frac{1}{
{2\pi }}\int_{ – \infty }^\infty {
{P_s}\left( \omega \right)d\omega } \\ {P_s}\left( \omega \right) \Leftrightarrow R\left( \tau \right)
Ps(ω)=TlimTE[FT(ω)2]=R(τ)eτdτR(τ)=2π1Ps(ω)eτdωS=R(0)=2π1Ps(ω)dωPs(ω)R(τ)

在维纳-辛钦定理的基础上,可以得出以下结论
(1)对功率谱密度进行积分,可以得到平稳过程的平均功率:
R ( 0 ) = ∫ − ∞ ∞ P ξ ( f ) d f R(0) = \int_{-\infty}^{\infty}P_{\xi}(f)df R(0)=Pξ(f)df
(2)各态历经过程的任一样本函数的功率谱密度等于过程的功率谱密度。
(3)功率谱密度 P ξ ( f ) P_{\xi}(f) Pξ(f)具有非负性和实偶性。

3、高斯随机过程

3.1 定义

如果随机过程 ξ ( t ) \xi(t) ξ(t)的任意n维分布均服从正态分布,则称它为正态过程或高斯过程。其n维正态概率密度函数表示如下:
f n ( x 1 , x 2 , . . . x n ; t 1 , t 2 , . . . t n ) = 1 ( 2 π ) n / 2 σ 1 σ 2 ⋯ σ n ∣ B ∣ 1 / 2 e [ − 1 2 ∣ B ∣ ∑ j = 1 n ∑ k = 1 n ∣ B ∣ j k ( x j − a j σ j ) ( x k − a k σ k ) ] \begin{array}{l} {f_n}\left( {
{x_1},{x_2},…{x_n};{t_1},{t_2},…{t_n}} \right) = \frac{1}{
{
{
{\left( {2\pi } \right)}^{n/2}}{\sigma _1}{\sigma _2} \cdots {\sigma _n}{
{\left| B \right|}^{1/2}}}} e^{\left[ {\frac{
{ – 1}}{
{2\left| B \right|}}\sum\nolimits_{j = 1}^n {
{
{\sum\nolimits_{k = 1}^n {\left| B \right|} }_{jk}}\left( {\frac{
{
{x_j} – {a_j}}}{
{
{\sigma _j}}}} \right)\left( {\frac{
{
{x_k} – {a_k}}}{
{
{\sigma _k}}}} \right)} } \right]} \end{array}
fn(x1,x2,xn;t1,t2,tn)=(2π)n/2σ1σ2σnB1/21e[2B1j=1nk=1nBjk(σjxjaj)(σkxkak)]

式中: a k = E [ ξ ( t k ) ] , σ k 2 = E [ ξ ( t k ) − a k ] ; ∣ B ∣ a_k = E[\xi(t_k)],\sigma^2_k = E[\xi(t_k) – a_k];|B| ak=E[ξ(tk)],σk2=E[ξ(tk)ak];B为归一化协方差矩阵的行列式,即:
∣ B ∣ = ∣ 1 b 12 ⋯ b 1 n b 21 1 ⋯ b 2 n ⋮ ⋮ ⋯ ⋮ b n 1 b n 2 ⋯ 1 ∣ |B| = \left|\begin{matrix} 1&b_{12}&\cdots&b_{1n}\\b_{21}&1&\cdots&b_{2n}\\ \vdots&\vdots&\cdots&\vdots\\b_{n1}&b_{n2}&\cdots&1 \end{matrix}\right| B=
1b21bn1b121bn2b1nb2n1

∣ B ∣ j k |B|_{jk} Bjk是行列式 ∣ B ∣ |B| B中元素 b j k b_{jk} bjk的代数余子式; b j k b_{jk} bjk为归一化协方差函数,即:
b j k = E ∣ [ ξ ( t j ) − a k ] [ ξ ( t K ) − a k ] ∣ σ j σ k b_{jk} = \frac {E|[\xi(t_j) – a_k][\xi(t_K) – a_k]|}{\sigma_j\sigma_k} bjk=σjσkE[ξ(tj)ak][ξ(tK)ak]

3.2 重要性质

(1)高斯分布的n维分布只依赖各个随机变量的均值、方差和归一化协方差。
(2)广义平稳的高斯过程也是严平稳的。
(3)如果高斯过程在不同时刻的取值是不相关的,则:
f n ( x 1 , x 2 , . . . x n ; t 1 , t 2 , . . . t n ) = f n ( x 1 ; t 1 ) f n ( x 2 ; t 2 ) ⋯ f n ( x n ; t n ) {f_n}\left( {
{x_1},{x_2},…{x_n};{t_1},{t_2},…{t_n}} \right) = {f_n}\left( {
{x_1};{t_1}} \right){f_n}\left( {
{x_2};{t_2}} \right) \cdots {f_n}\left( {
{x_n};{t_n}} \right)
fn(x1,x2,xn;t1,t2,tn)=fn(x1;t1)fn(x2;t2)fn(xn;tn)

这表明,如果高斯过程在不同时刻的取值是不相关的,那么他们也是统计独立的。
(4)高斯过程经过线性变换后生成的过程仍是高斯过程。

3.3 高斯随机变量

高斯过程在任一时刻的取值是一个正态分布的随机变量,也称高斯随机变量,其以为概率密度函数为:
f ( x ) = 1 2 π σ exp ⁡ [ − ( x − a ) 2 2 σ 2 ] f\left( x \right) = \frac{1}{
{\sqrt {2\pi } \sigma }}\exp \left[ { – \frac{
{
{
{\left( {x – a} \right)}^2}}}{
{2{\sigma ^2}}}} \right]
f(x)=2π
σ
1
exp[2σ2(xa)2]

特性:

  • 对称于 x = a x=a xa这条直线,即有 f ( a + x ) = f ( a − x ) f(a+x) = f(a-x) f(a+x)=f(ax)
  • ∫ − ∞ ∞ f ( x ) d x = 1 \int_{ – \infty }^\infty {f\left( x \right)dx = 1} f(x)dx=1 ∫ − ∞ a f ( x ) d x = ∫ a ∞ f ( x ) d x = 1 / 2 \int_{ – \infty }^a {f\left( x \right)dx} = \int_a^\infty {f\left( x \right)dx = 1/2} af(x)dx=af(x)dx=1/2
  • 在点a处达到极大值
  • σ \sigma σ不变,增加a,f(x)图形左移;减小a,f(x)图形右移;若a不变,增加 σ \sigma σ,f(x)图形变低;减小 σ \sigma σ,f(x)图形变高;

把正态分布的概率密度的积分定义为正态分布函数:
F ( x ) = ∫ − ∞ x f ( x ) d x = 1 2 π σ ∫ − ∞ x exp ⁡ [ − ( z − a ) 2 2 σ 2 ] d z = 1 2 + e r f ( x − a 2 a ) F\left( x \right) = \int_{ – \infty }^x {f\left( x \right)dx} = \frac{1}{
{\sqrt {2\pi } \sigma }}\int_{ – \infty }^x {\exp \left[ { – \frac{
{
{
{\left( {z – a} \right)}^2}}}{
{2{\sigma ^2}}}} \right]dz} = \frac{1}{2} + erf\left( {\frac{
{x – a}}{
{\sqrt 2 a}}} \right)
F(x)=xf(x)dx=2π
σ
1
xexp[2σ2(za)2]dz=
21+erf(2
a
xa
)

e r f ( x ) erf(x) erf(x)表示误差函数,定义为:
e r f ( x ) = 2 π ∫ 0 x e − t 2 d t erf(x) = \frac 2{\sqrt{\pi}}\int_0^xe^{-t^2}dt erf(x)=π
2
0xet2dt

F(x)也可以用互补误差函数erfc(x)表示,即
F ( x ) = 1 − 0.5 ∗ e r f c ( x − a 2 σ ) F(x) = 1 – 0.5*erfc(\frac {x – a}{\sqrt2\sigma}) F(x)=10.5erfc(2
σ
xa
)

e r f c ( x ) erfc(x) erfc(x)定义为:
e r f c ( x ) = 1 − e r f ( x ) = 2 π ∫ x ∞ e t 2 d t erfc(x) = 1 – erf(x) = \frac 2 {\sqrt \pi}\int_x^\infty e^{t^2}dt erfc(x)=1erf(x)=π
2
xet2dt

当x很大时(实际应用中只要x>2),互补误差函数可以近似为
e r f c ( x ) = 1 x π e − x 2 erfc(x) = \frac 1 {x\sqrt \pi}e^{-x^2} erfc(x)=xπ
1
ex2

4、平稳随机过程通过线性系统

假设输入过程 ξ i ( t ) \xi_i(t) ξi(t)是平稳的,其均值为a,自相关函数为 R i ( τ ) R_i(\tau) Ri(τ),功率谱密度为 P i ( ω ) P_i(\omega) Pi(ω)

  • 输出过程 ξ o ( t ) \xi_o(t) ξo(t)的均值

E [ ξ o ( t ) ] = a ∫ − ∞ ∞ h ( τ ) d τ = a H ( 0 ) E[\xi_o(t)] = a\int_{-\infty}^\infty h(\tau)d\tau = aH(0) E[ξo(t)]=ah(τ)dτ=aH(0)

  • 输出过程 ξ o ( t ) \xi_o(t) ξo(t)的自相关函数

R o ( t 1 , t 1 + τ ) = ∫ − ∞ ∞ ∫ − ∞ ∞ h ( α ) h ( β ) R i ( τ + α − β ) d α d β R_o(t_1,t_1+\tau) = \int_{-\infty}^{\infty}\int_{-\infty}^{\infty}h(\alpha)h(\beta)R_i(\tau+\alpha-\beta)d\alpha d\beta Ro(t1,t1+τ)=h(α)h(β)Ri(τ+αβ)dαdβ
若线性系统的输入过程是平稳的,那么输出过程也是平稳的。

  • 输出过程 ξ o ( t ) \xi_o(t) ξo(t)的功率谱

P o ( f ) = H ∗ ( f ) H ( f ) P i ( f ) = ∣ H ( f ) ∣ 2 P i ( f ) P_o(f) = H^*(f)H(f)P_i(f) = |H(f)|^2P_i(f) Po(f)=H(f)H(f)Pi(f)=H(f)2Pi(f)
输出过程的功率谱密度是输入过程的功率谱密度乘以系统频率响应模值的平方。

  • 输出过程 ξ o ( t ) \xi_o(t) ξo(t)的概率分布

高斯过程经线性变换后的过程仍为高斯过程

5、窄带随机过程

若随机过程 ξ ( t ) \xi(t) ξ(t)的谱密度集中在中心频率 f c f_c fc附近相对窄的频带范围 Δ f \Delta f Δf内,且 f c f_c fc远离零频率,则称该 ξ ( t ) \xi(t) ξ(t)为窄带随机过程。
窄带随机过程的一个样本的波形如同包络和相位随机缓变的正弦波。可表示为:
ξ ( t ) = a ξ ( t ) cos ⁡ [ ω c t + ϕ ξ ] \xi(t) = a_{\xi}(t)\cos[\omega_ct + \phi_{\xi}] ξ(t)=aξ(t)cos[ωct+ϕξ]
上式可以进行三角函数展开:
ξ ( t ) = ξ c ( t ) cos ⁡ ( ω c t ) − ξ s ( t ) sin ⁡ ( ω c t ) \xi(t) = \xi_c(t)\cos(\omega_ct) – \xi_s(t)\sin(\omega_ct) ξ(t)=ξc(t)cos(ωct)ξs(t)sin(ωct)
其中:
ξ c ( t ) = a ξ ( t ) cos ⁡ ( ϕ ξ ( t ) ) ξ s ( t ) = a ξ ( t ) sin ⁡ ( ϕ ξ ( t ) ) \xi_c(t) = a_{\xi}(t)\cos(\phi_{\xi}(t))\\ \xi_s(t) = a_{\xi}(t)\sin(\phi_{\xi}(t)) ξc(t)=aξ(t)cos(ϕξ(t))ξs(t)=aξ(t)sin(ϕξ(t))
这里的 ξ c ( t ) \xi_c(t) ξc(t) ξ s ( t ) \xi_s(t) ξs(t)分别被称为同相分量和正交分量。

假设 ξ ( t ) \xi(t) ξ(t)是一个均值为0,方差为 σ ξ 2 \sigma_{\xi}^2 σξ2的平稳高斯窄带过程。

5.1 ξ c ( t ) \xi_c(t) ξc(t) ξ s ( t ) \xi_s(t) ξs(t)的统计特性

ξ ( t ) \xi(t) ξ(t)的数学期望:
E [ ξ ( t ) ] = E [ ξ c ( t ) ] cos ⁡ ( ω c t ) − E [ ξ s ( t ) ] sin ⁡ ( ω c t ) E[\xi(t)] = E[\xi_c(t)]\cos(\omega_ct) – E[\xi_s(t)]\sin(\omega_ct) E[ξ(t)]=E[ξc(t)]cos(ωct)E[ξs(t)]sin(ωct)
因为 ξ ( t ) \xi(t) ξ(t)平稳且均值为0,那么对于任意的时间t,都有 E [ ξ ( t ) ] = 0 E[\xi(t)] = 0 E[ξ(t)]=0,所以可得:
E [ ξ c ( t ) ] = 0 ,      E [ ξ s ( t ) ] = 0 E[\xi_c(t)] = 0,~~~~E[\xi_s(t)] = 0 E[ξc(t)]=0,    E[ξs(t)]=0
ξ ( t ) \xi(t) ξ(t)的自相关函数:
R ξ ( t , t + τ ) = E [ ξ ( t ) ξ ( t + τ ) ] = R c ( t , t + τ ) cos ⁡ ( ω c t ) cos ⁡ ( ω c ( t + τ ) ) − R c s ( t , t + τ ) cos ⁡ ( ω c t ) sin ⁡ ( ω c ( t + τ ) ) − R s c ( t , t + τ ) sin ⁡ ( ω c t ) cos ⁡ ( ω c ( t + τ ) ) + R s ( t + τ ) sin ⁡ ( ω c t ) sin ⁡ ( ω c ( t + τ ) ) R_{\xi}(t,t+\tau) = E[\xi(t)\xi(t+\tau)]\\ = R_c(t,t+\tau)\cos(\omega_ct)\cos(\omega_c(t+\tau)) – R_{cs}(t,t+\tau)\cos(\omega_ct)\sin(\omega_c(t+\tau)) – \\ R_{sc}(t,t+\tau)\sin(\omega_ct)\cos(\omega_c(t+\tau)) + R_s(t+\tau)\sin(\omega_ct)\sin(\omega_c(t+\tau)) Rξ(t,t+τ)=E[ξ(t)ξ(t+τ)]=Rc(t,t+τ)cos(ωct)cos(ωc(t+τ))Rcs(t,t+τ)cos(ωct)sin(ωc(t+τ))Rsc(t,t+τ)sin(ωct)cos(ωc(t+τ))+Rs(t+τ)sin(ωct)sin(ωc(t+τ))
其中:
R c ( t , t + τ ) = E [ ξ c ( t ) ξ c ( t + τ ) ] R c s ( t , t + τ ) = E [ ξ c ( t ) ξ s ( t + τ ) ] R s c ( t , t + τ ) = E [ ξ s ( t ) ξ c ( t + τ ) ] R s ( t , t + τ ) = E [ ξ s ( t ) ξ s ( t + τ ) ] R_c(t,t+\tau) = E[\xi_c(t)\xi_c(t+\tau)]\\ R_{cs}(t,t+\tau) = E[\xi_c(t)\xi_s(t+\tau)]\\ R_{sc}(t,t+\tau) = E[\xi_s(t)\xi_c(t+\tau)]\\ R_s(t,t+\tau) = E[\xi_s(t)\xi_s(t+\tau)] Rc(t,t+τ)=E[ξc(t)ξc(t+τ)]Rcs(t,t+τ)=E[ξc(t)ξs(t+τ)]Rsc(t,t+τ)=E[ξs(t)ξc(t+τ)]Rs(t,t+τ)=E[ξs(t)ξs(t+τ)]
因为 ξ ( t ) \xi(t) ξ(t)是平稳的,故有:
R ξ ( t + τ ) = R ( τ ) R_{\xi}(t+\tau) = R(\tau) Rξ(t+τ)=R(τ)
因此,若令 t = 0 t = 0 t=0,有:
R ξ ( τ ) = R c ( τ ) cos ⁡ ( ω c τ ) − R c s ( τ ) sin ⁡ ( ω c τ ) R_{\xi}(\tau) = R_c(\tau)\cos(\omega_c\tau) – R_{cs}(\tau)\sin(\omega_c\tau) Rξ(τ)=Rc(τ)cos(ωcτ)Rcs(τ)sin(ωcτ)
同理,若令 t = π 2 ω c t = \frac {\pi}{2\omega_c} t=2ωcπ,有:
R ξ ( τ ) = R s ( τ ) cos ⁡ ( ω c τ ) + R s c ( τ ) sin ⁡ ( ω c τ ) R_{\xi}(\tau) = R_s(\tau)\cos(\omega_c\tau) + R_{sc}(\tau)\sin(\omega_c\tau) Rξ(τ)=Rs(τ)cos(ωcτ)+Rsc(τ)sin(ωcτ)
由此可知,若窄带随机过程是平稳的,则其同相分量和正交分量也必然是平稳的。
一些性质:
R c ( τ ) = R s ( τ ) R c s ( τ ) = − R s c ( τ ) R c s ( t a u = R s c ( − τ ) R s c ( t a u = − R s c ( − τ ) R s c ( 0 ) = R c s ( 0 ) = 0 R ξ ( 0 ) = R c ( 0 ) = R s ( 0 ) σ ξ 2 = σ c 2 = σ s 2 R_c(\tau) = R_s(\tau)\\ R_{cs}(\tau) = – R_{sc}(\tau)\\ R_{cs}(tau = R_{sc}( – \tau)\\ R_{sc}(tau = – R_{sc}( – \tau)\\ R_{sc}(0) = R_{cs}(0) = 0\\ R_{\xi}(0) = R_c(0) = R_s(0)\\ \sigma_{\xi}^2 = \sigma_c^2 = \sigma_s^2 Rc(τ)=Rs(τ)Rcs(τ)=Rsc(τ)Rcs(tau=Rsc(τ)Rsc(tau=Rsc(τ)Rsc(0)=Rcs(0)=0Rξ(0)=Rc(0)=Rs(0)σξ2=σc2=σs2

5.2 a ξ ( t ) a_{\xi}(t) aξ(t) ϕ ξ ( t ) \phi_{\xi}(t) ϕξ(t)的统计特性

由上面的分析可知, ξ c , ξ s \xi_c,\xi_s ξc,ξs的联合概率密度函数为:
f ( ξ c , ξ s ) = f ( ξ c ) f ( ξ s ) = 1 2 π σ ξ 2 e x p [ − ξ c 2 + ξ s 2 2 σ ξ 2 ] f(\xi_c,\xi_s) = f(\xi_c)f(\xi_s) = \frac {1}{2\pi \sigma^2_{\xi}}exp{[-\frac {\xi_c^2 + \xi_s^2}{2\sigma^2_{\xi}}]} f(ξc,ξs)=f(ξc)f(ξs)=2πσξ21exp[2σξ2ξc2+ξs2]
a ξ , ϕ ξ a_{\xi},\phi_{\xi} aξ,ϕξ的联合概率密度函数为 f ( a ξ , ϕ ξ ) f(a_{\xi},\phi_{\xi}) f(aξ,ϕξ),则有:
f ( a ξ , ϕ ξ ) = f ( ξ c , ξ s ) ∣ ∂ ( ξ c , ξ s ) ∂ ( a ξ , ϕ ξ ) ∣ f(a_{\xi},\phi_{\xi}) = f(\xi_c,\xi_s)|\frac {\partial(\xi_c,\xi_s)}{\partial(a_{\xi},\phi_{\xi})}| f(aξ,ϕξ)=f(ξc,ξs)(aξ,ϕξ)(ξc,ξs)
可以求得:
f ( a ξ , ϕ ξ ) = a ξ f ( ξ c , ξ s ) = a ξ 2 π σ ξ 2 e x p [ − ( a ξ cos ⁡ ( ϕ ξ ) ) 2 + ( a ξ sin ⁡ ( ϕ ξ ) ) 2 2 σ ξ ] = a ξ 2 π σ ξ 2 e x p [ − a ξ 2 2 σ ξ 2 ] f(a_{\xi},\phi_{\xi}) = a_{\xi}f(\xi_c,\xi_s) = \frac{a_{\xi}}{2\pi\sigma_{\xi}^2}exp{[-\frac {(a_{\xi}\cos(\phi_{\xi}))^2 + (a_{\xi}\sin(\phi_{\xi}))^2}{2\sigma_{\xi}}]}\\ = \frac{a_{\xi}}{2\pi\sigma_{\xi}^2}exp[-\frac {a_{\xi}^2}{2\sigma_{\xi}^2}] f(aξ,ϕξ)=aξf(ξc,ξs)=2πσξ2aξexp[2σξ(aξcos(ϕξ))2+(aξsin(ϕξ))2]=2πσξ2aξexp[2σξ2aξ2]
分别对 a ξ , ϕ ξ a_{\xi},\phi_{\xi} aξϕξ求积分求得一维概率密度函数为:
f ( a ξ ) = a ξ σ ξ 2 e x p [ − a ξ 2 2 σ ξ 2 ] f ( ϕ ξ ) = 1 2 π f(a_{\xi}) = \frac{a_{\xi}}{\sigma_{\xi}^2}exp[-\frac {a_{\xi}^2}{2\sigma_{\xi}^2}]\\ f(\phi_{\xi}) = \frac 1{2\pi} f(aξ)=σξ2aξexp[2σξ2aξ2]f(ϕξ)=2π1
可见, a ξ ( t ) a_{\xi}(t) aξ(t)服从瑞利分布, ϕ ξ ( t ) \phi_{\xi}(t) ϕξ(t)服从均匀分布。且 a ξ ( t ) 和 ϕ ξ ( t ) a_{\xi}(t)和\phi_{\xi}(t) aξ(t)ϕξ(t)是统计独立的,即:
f ( a ξ , ϕ ξ ) = f ( a ξ ) f ( ϕ ξ ) f(a_{\xi},\phi_{\xi}) = f(a_{\xi})f(\phi_{\xi}) f(aξ,ϕξ)=f(aξ)f(ϕξ)

6、正弦波加窄带高斯噪声

在许多调制系统中,传输的信号是用正弦波作为载波的已调信号。并且经过信道传输时会有噪声影响,为了减少噪声的影响,通常在解调器前端设置一个带通滤波器,其输出是正弦波已调信号与窄带高斯噪声的混合波形。
设正弦波加窄带高斯噪声的混合信号为:
r ( t ) = A cos ⁡ ( ω c t + θ ) + n ( t ) r(t) = A\cos(\omega_ct + \theta) + n(t) r(t)=Acos(ωct+θ)+n(t)
式中: n ( t ) = n c ( t ) cos ⁡ ( ω c t ) − n s ( t ) sin ⁡ ( ω c t ) ) n(t) = n_c(t)\cos(\omega_ct) – n_s(t)\sin(\omega_ct)) n(t)=nc(t)cos(ωct)ns(t)sin(ωct)),为窄带高斯噪声,其均值为0,方差为 σ n 2 \sigma_n^2 σn2 θ \theta θ为正弦波的随机相位,在 ( 0 , 2 π ) (0,2\pi) (0,2π)上均与分布;整幅A和角频率 ω c \omega_c ωc均假定为确知量。
于是:
r ( t ) = [ A cos ⁡ ( θ ) + n c ( t ) ] cos ⁡ ( ω c t ) − [ A sin ⁡ ( θ ) + n s ( t ) ] sin ⁡ ( ω c t ) = z c ( t ) cos ⁡ ( ω c t ) − z s ( t ) sin ⁡ ( ω c t ) = z ( t ) cos ⁡ [ ω c t + ϕ ( t ) ] r(t) = [A\cos(\theta) + n_c(t)]\cos(\omega_ct) – [A\sin(\theta) + n_s(t)]\sin(\omega_ct)\\ = z_c(t)\cos(\omega_ct) – z_s(t)\sin(\omega_ct) = z(t)\cos[\omega_ct+\phi(t)] r(t)=[Acos(θ)+nc(t)]cos(ωct)[Asin(θ)+ns(t)]sin(ωct)=zc(t)cos(ωct)zs(t)sin(ωct)=z(t)cos[ωct+ϕ(t)]
其中:
z c ( t ) = A cos ⁡ ( θ ) + n c ( t ) z s ( t ) = A sin ⁡ ( θ ) + n s ( t ) z_c(t) = A\cos(\theta) + n_c(t)\\ z_s(t) = A\sin(\theta) + n_s(t) zc(t)=Acos(θ)+nc(t)zs(t)=Asin(θ)+ns(t)
可以求得包络z的概率密度函数为:
f ( z ) = z σ n 2 e x p [ − 1 2 σ n 2 ( z 2 + A 2 ) ] I O ( A z σ n 2 ) f(z) = \frac z{\sigma_n^2}exp[-\frac 1{2\sigma_n^2}(z^2 + A^2)]I_O(\frac {Az}{\sigma_n^2}) f(z)=σn2zexp[2σn21(z2+A2)]IO(σn2Az)
这个概率密度函数称为广义瑞利分布,又称莱斯分布
性质:
当信噪比 A 2 2 σ 2 很大时, f ( z ) 近似为高斯分布 当 A = 0 时,信噪比 A 2 2 σ 2 = 0 , f ( z ) 退化为瑞利分布 当信噪比\frac {A^2}{2\sigma^2}很大时,f(z)近似为高斯分布\\ 当A = 0时,信噪比\frac {A^2}{2\sigma^2} = 0,f(z)退化为瑞利分布 当信噪比2σ2A2很大时,f(z)近似为高斯分布A=0时,信噪比2σ2A2=0f(z)退化为瑞利分布

7、高斯白噪声和带限白噪声

白噪声定义:功率谱密度在整个频域内都是均匀分布的噪声。是理想的宽带过程。即双边功率谱密度为:
P n ( f ) = n 0 2 ( W / H z ) , ( − ∞ < f < ∞ ) {P_n}\left( f \right) = \frac{
{
{n_0}}}{2}\left( {W/Hz} \right),\left( { – \infty < f < \infty } \right)
Pn(f)=2n0(W/Hz),(<f<)

白噪声的自相关函数: R n ( τ ) ⇔ P n ( f ) ⇒ R n ( τ ) = n 0 2 δ ( t ) {R_n}\left( \tau \right) \Leftrightarrow {P_n}\left( f \right) \Rightarrow {R_n}\left( \tau \right) = \frac{
{
{n_0}}}{2}\delta (t)
Rn(τ)Pn(f)Rn(τ)=2n0δ(t)

带限白噪声:当白噪声通过信道时,频带将受到限制,称之为带限白噪声。分为:

  1. 理想低通白噪声:白噪声通过理想矩形的低通信道时就产生这类限带白噪声: R n ( τ ) = ∫ − f 0 f 0 n 0 2 e j 2 π f τ d f = f 0 n 0 sin ⁡ ( ω 0 τ ) ( ω 0 τ ) {R_n}\left(\tau\right) = \int_{ – {f_0}}^{
    {f_0}}{\frac{
    {
    {n_0}}}{2}{e^{j2\pi f\tau }}df = {f_0}} {n_0}\frac{
    {\sin\left( {
    {\omega _0}\tau } \right)}}{
    {\left({
    {\omega _0}\tau } \right)}}
    Rn(τ)=f0f02n0ej2πfτdf=f0n0(ω0τ)sin(ω0τ)
  2. 理想带通白噪声(窄带高斯白噪声): R n ( τ ) = ∫ − f 0 − B / 2 − f 0 + B / 2 n 0 2 e j 2 π f τ d f + ∫ f 0 − B / 2 f 0 + B / 2 n 0 2 e j 2 π f τ d f = n 0 B sin ⁡ ( π B τ ) ( π B τ ) cos ⁡ ( 2 π f 0 τ ) \begin{array}{l}{R_n}\left( \tau \right) = \int_{ – {f_0} – B/2}^{ – {f_0} + B/2} {\frac{
    {
    {n_0}}}{2}{e^{j2\pi f\tau }}df} + \int_{
    {f_0} – B/2}^{
    {f_0} + B/2} {\frac{
    {
    {n_0}}}{2}{e^{j2\pi f\tau }}df} \\= {n_0}B\frac{
    {\sin \left( {\pi B\tau } \right)}}{
    {\left( {\pi B\tau } \right)}}\cos (2\pi {f_0}\tau )\end{array}
    Rn(τ)=f0B/2f0+B/22n0ej2πfτdf+f0B/2f0+B/22n0ej2πfτdf=n0B(πBτ)sin(πBτ)cos(2πf0τ)

免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://haidsoft.com/144198.html

(0)
上一篇 2025-04-28 14:33
下一篇 2025-04-28 14:45

相关推荐

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注微信