通信原理(2)--随机过程

通信原理(2)–随机过程

3.1随机过程的基本概念

随机过程{x(t)}由一族时间函数 x i ( t ) x_i(t) xi(t),i=1,2.3…组成,每一个时间函数 x i ( t ) x_i(t) xi(t)称为随机过程{x(t)}的一个样本函数(一个实现)
每个样本函数在时间上,在幅度取值上都是连续变化的波形。
若固定时刻,随机过程在该时刻所有样本函数的取值则为一随机变量
为了描述随机过程在不同时刻的相互关系,用n维联合分布函数(概率密度函数族)来描述n个不同时刻对应的n个随机变量。

image-20240421161957462

image-20240421162306892

随机过程其实是随机变量在连续时间上的一次升维

3.1.1随机过程的分布函数

一维分布函数
F 1 ( x 1 , t 1 ) = P [ ξ ( t 1 ) ≤ x 1 ] F_1(x_1,t_1)=P\Big[\xi(t_1)\leq x_1\Big] F1(x1,t1)=P[ξ(t1)x1]

{ 0 ≤ F 1 ( x 1 , t 1 ) ≤ 1 F 1 ( − ∞ , t 1 ) = 0 , F 1 ( + ∞ , t 1 ) = 1 x 2 > x 1 , F 1 ( x 2 , t 1 ) > F 1 ( x 1 , t 1 ) \begin{cases}0\leq F_1(x_1,t_1)\leq1\\F_1(-\infty,t_1)=0,F_1(+\infty,t_1)=1\\x_2>x_1,F_1(x_2,t_1)>F_1(x_1,t_1)\end{cases} 0F1(x1,t1)1F1(,t1)=0,F1(+,t1)=1x2>x1,F1(x2,t1)>F1(x1,t1)

一维概率密度函数
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 ) ≥ 0 F 1 ( x 1 , t 1 ) = ∫ − ∞ x 1 f 1 ( x 1 ′ , t 1 ) d x 1 ′ P [ a 1 ≤ ξ ( t 1 ) ≤ a 1 ] = ∫ a 1 a 2 f 1 ( x 1 , t 1 ) d x 1 ∫ − ∞ + ∞ f 1 ( x 1 , t 1 ) d x 1 = 1 \begin{cases}f_1(x_1,t_1)\geq0\\[2ex]F_1(x_1,t_1)=\int_{-\infty}^{x_1}f_1(x_1^{\prime},t_1)\mathrm{d}x_1^{\prime}\\[2ex]P\Big[a_1\leq\xi(t_1)\leq a_1\Big]=\int_{a_1}^{a_2}f_1(x_1,t_1)\mathrm{d}x_1\\[2ex]\int_{-\infty}^{+\infty}f_1(x_1,t_1)\mathrm{d}x_1=1\end{cases} f1(x1,t1)0F1(x1,t1)=x1f1(x1,t1)dx1P[a1ξ(t1)a1]=a1a2f1(x1,t1)dx1+f1(x1,t1)dx1=1

n维分布函数与n维概率密度函数

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 ] F_n(x_1,x_2,\cdots,x_n;t_1,t_2,\cdots,t_n)=P\Big[\xi(t_1)\leq x_1,\xi(t_2)\leq x_2,\cdots,\xi(t_n)\leq x_n\Big] 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)

数学期望

a ( t ) = E [ ξ ( t ) ] = ∫ − ∞ + ∞ x f 1 ( x , t ) d x a(t)=E\Big[\xi(t)\Big]=\int_{-\infty}^{+\infty}xf_1(x,t)\mathrm{d}x a(t)=E[ξ(t)]=+xf1(x,t)dx

道机过程的数学期望与t有关,表示n个样本函数曲线的摆动中心。

E[C]=C常数的数学期望是常数本身:
E [ a X 1 + b X 2 ] = a E [ X 1 ] + b E [ X 2 ] E\left[aX_1+bX_2\right]=aE\left[X_1\right]+bE\left[X_2\right] E[aX1+bX2]=aE[X1]+bE[X2]
对于两个统计独立的随机变量 X 1 , X 2 ,有 E [ X 1 X 2 ] = E [ X 1 ] E [ X 2 ] X_1,X_2\text{,有}E[X_1X_2]=E[X_1]E[X_2] X1,X2,E[X1X2]=E[X1]E[X2]

方差

D [ ξ ( t ) ] = E { [ ξ ( t ) − a ( t ) ] 2 } = E [ ξ 2 ( t ) ] − a 2 ( t ) D\Big[\xi(t)\Big]=E\Big\{\Big[\xi(t)-a(t)\Big]^2\Big\}=E\Big[\xi^2(t)\Big]-a^2(t) D[ξ(t)]=E{[ξ(t)a(t)]2}=E[ξ2(t)]a2(t)

随机过程的方差表示了随机过程 ξ ( t ) \xi(t) ξ(t)在t时刻偏离数学期望的程度。

D [ C ] = 0 D[C]=0 D[C]=0 常数的方差是0;

D [ C X ] = C 2 D [ X ] ,   D [ X + C ] = D [ X ] D[CX]=C^{2}D[X],\:D[X+C]=D[X] D[CX]=C2D[X],D[X+C]=D[X]
对于两个随机变量 X 1 , X 2 X_1,X_2 X1,X2, 有
D [ a X 1 + b X 2 ] = a 2 D [ X 1 ] + b 2 D [ X 2 ] + 2 a b ⋅ C o v [ X 1 , X 2 ] D\left[aX_1+bX_2\right]=a^2D\left[X_1\right]+b^2D\left[X_2\right]+2ab\cdot\mathrm{Cov}\left[X_1,X_2\right] D[aX1+bX2]=a2D[X1]+b2D[X2]+2abCov[X1,X2]
对于两个不相关的随机变量 X 1 , X 2 X_1,X_2 X1,X2, 有

D [ a X 1 + b X 2 ] = a 2 D [ X 1 ] + b 2 D [ X 2 ] D\Big[aX_1+bX_2\Big]=a^2D\Big[X_1\Big]+b^2D\Big[X_2\Big] D[aX1+bX2]=a2D[X1]+b2D[X2]

3.1.2随机过程的统计特性

数学期望和方差只能反映随机变量 ξ ( 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(t_{1},t_{2})=E\Big[\xi(t_{1})\xi(t_{2})\Big]=\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}x_{1}x_{2}f_{2}(x_{1},x_{2};t_{1},t_{2})\mathrm{d}x_{1}\mathrm{d}x_{2} R(t1,t2)=E[ξ(t1)ξ(t2)]=++x1x2f2(x1,x2;t1,t2)dx1dx2

反映不同时刻随机过程取值的相关性。

  • 自协方差函数

B ( t 1 , t 2 ) = E { [ ξ ( t 1 ) − a ( t 1 ) ] [ ξ ( t 2 ) − a ( t 2 ) ] } B(t_1,t_2)=E\left\{\left[\xi(t_1)-a(t_1)\right]\left[\xi(t_2)-a(t_2)\right]\right\} B(t1,t2)=E{[ξ(t1)a(t1)][ξ(t2)a(t2)]}

= ∫ − ∞ + ∞ ∫ − ∞ + ∞ [ x 1 − a ( t 1 ) ] [ x 2 − a ( t 2 ) ] f 2 ( x 1 , x 2 ; t 1 , t 2 ) d x 1 d x 2 =\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}\Big[x_{1}-a(t_{1})\Big]\Big[x_{2}-a(t_{2})\Big]f_{2}(x_{1},x_{2};t_{1},t_{2})\mathrm{d}x_{1}\mathrm{d}x_{2} =++[x1a(t1)][x2a(t2)]f2(x1,x2;t1,t2)dx1dx2

自协方差函数反映不同时刻随机过程的统计相关性。

  • 自相关函数与自协方差函数的关系

B ( t 1 , t 2 ) = R ( t 1 , t 2 ) − a ( t 1 ) a ( t 2 ) B(t_1,t_2)=R(t_1,t_2)-a(t_1)a(t_2) B(t1,t2)=R(t1,t2)a(t1)a(t2)

对于零均值随机过程: B ( t 1 , t 2 ) = R ( t 1 , t 2 ) B(t_1,t_2)=R(t_1,t_2) B(t1,t2)=R(t1,t2)
假定 t 2 > t 1 t_2>t_1 t2>t1,令 t 2 = t 1 + τ t_2=t_1+\tau t2=t1+τ,则 R ( t 1 , t 2 ) = R ( t 1 , t 1 + τ )   ; R(t_1,t_2)=R(t_1,t_1+\tau)\:; R(t1,t2)=R(t1,t1+τ);

我们一般研究时间间隔 τ \tau τ 对相关性的影响;

一般情况下, τ \tau τ越大,相关性越小

对于两个随机过程 ξ ( t ) \xi(t) ξ(t) η ( t ) \eta(t) η(t), 定义:

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

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

3.2平稳随机过程

3.2.1定义

平稳随机过程的统计特性不随时间的推移而变化,其一维分布与时间 t t t 无关,二维分布只与时间间隔 τ \tau τ 有关。即:

一维分布与时间t无关
f 1 ( x , t ) = f 1 ( x , t + h ) = f 1 ( x ) ∀ h f_1(x,t)=f_1(x,t+h)=f_1(x)\quad\forall h f1(x,t)=f1(x,t+h)=f1(x)h

二维分布只与时间间隔 τ \tau τ有关
f 2 ( x 1 , x 2 ; t 1 , t 2 ) = f 2 ( x 1 , x 2 ; t 1 + h , t 2 + h ) f_2(x_1,x_2;t_1,t_2)=f_2(x_1,x_2;t_1+h,t_2+h) f2(x1,x2;t1,t2)=f2(x1,x2;t1+h,t2+h)

= f 2 ( x 1 , x 2 ; t 2 − t 1 ) = f 2 ( x 1 , x 2 ; τ ) ∀ h =f_2(x_1,x_2;t_2-t_1)=f_2(x_1,x_2;\tau)\quad\forall h =f2(x1,x2;t2t1)=f2(x1,x2;τ)h

3.2.2 平稳随机过程的数字特征

均值: E [ ξ ( t ) ] = ∫ − ∞ + ∞ x f 1 ( x ) E\left [ \xi( t) \right ] = \int _{- \infty}^{+ \infty}xf_{1}( x) E[ξ(t)]=+xf1(x)d x = α x= \alpha x=α t t t无关,样本函数围绕一水平线起伏

方差: D [ ξ ( t ) ] = ∫ − α + ∞ ( x − a ) 2 f 1 ( x ) D\left [ \xi( t) \right ] = \int _{- \alpha}^{+ \infty}( x- a) ^{2}f_{1}( x) D[ξ(t)]=α+(xa)2f1(x)d x = E [ ξ 2 ( t ) ] − a 2 = σ 2 x= E\left [ \xi^2( t) \right ] - a^{2}= \sigma^{2} x=E[ξ2(t)]a2=σ2 与== t t t无关==,为常数

自相关函数:

R ( t , t + τ ) = E [ ξ ( t ) ξ ( t + τ ) ] R(t,t+\tau)=E\bigl[\xi(t)\xi(t+\tau)\bigr] R(t,t+τ)=E[ξ(t)ξ(t+τ)]

= ∫ − ∞ + ∞ ∫ − ∞ + ∞ x 1 x 2 f 2 ( x 1 , x 2 ; τ ) = \int _{- \infty}^{+ \infty}\int _{- \infty}^{+ \infty}x_{1}x_{2}f_{2}( x_{1}, x_{2}; \tau) =++x1x2f2(x1,x2;τ)d x 1 x_1 x1d x 2 = R ( τ ) x_2= R( \tau) x2=R(τ)
自相关只与时间间隔 τ \tau τ有关,而与时间起点无关。

3.2.3 狭义(严)平稳和广义(宽)平稳

  • 狭义平稳随机过程(严平稳)

随机过程的各维概率密度函数都不随时间的推移而变化。

f n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) = f n ( x 1 , x 2 , ⋯   , x n ; t 1 + h , t 2 + h , ⋯   , t n + h ) ,   ∀ n , h 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+h,t_2+h,\cdots,t_n+h),\:\forall n,h fn(x1,x2,,xn;t1,t2,,tn)=fn(x1,x2,,xn;t1+h,t2+h,,tn+h),n,h
狭义平稳过程又称严平稳过程,条件很强,对其作放宽,可得广义平稳、即宽平稳过程的概念。

  • 广义平稳随机过程(宽平稳)

一个二阶矩随机过程,均值为常数自相关函数仅仅是时间间隔的函数,称这样的过程为广义平稳随机过程或宽平稳随机过程。

image-20240421172255107

后文的平稳随机过程均为宽平稳(除特指)

3.2.3各态历经性

  • 平稳随机过程的各态经历性
    假定 x ( t ) x(t) 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 \overline{a}=\overline{x(t)}=\lim_{T\to\infty}\frac{1}{T}\int_{-T/2}^{T/2}x(t)\mathrm{d}t\quad\overline{R(\tau)}=\overline{x(t)x(t+\tau)}=\lim_{T\to\infty}\frac{1}{T}\int_{-T/2}^{T/2}x(t)x(t+\tau)\mathrm{d}t a=x(t)=TlimT1T/2T/2x(t)dtR(τ)=x(t)x(t+τ)=TlimT1T/2T/2x(t)x(t+τ)dt
若平稳随机过程依概率1满足 a ˉ = a \bar{a}=a aˉ=a, R ( τ ) ‾ = R ( τ ) \overline{R(\tau)}=R(\tau) R(τ)=R(τ)
则称该平稳过程具有各态历经性。

  • 含义

随机过程中任一实现都经历了随机过程的所有可能状态。

化“统计平均”为“时间平均”, 用任意一个样本函数刻画整个随机过程的所有特征,简化实际的测量和计算。

3.2.4平稳过程的自相关函数

3.2.4.1平稳随机过程自相关函数的性质

R ( 0 ) = E [ ξ 2 ( t ) ] = S R(0)=E\left[\xi^2(t)\right]=S R(0)=E[ξ2(t)]=S 物理意义:随机过程的平均功率;

R ( ∞ ) = E 2 [ ξ ( t ) ] = a 2 R(\infty)=E^2\left[\xi(t)\right]=a^2 R()=E2[ξ(t)]=a2 物理意义:随机过程的直流功率;

时间间隔无限大时, ξ ( t ) \xi(t) ξ(t) ξ ( t + τ ) \xi(t+\tau) ξ(t+τ)趋于独立

由于平稳, E [ ξ ( t + τ ) ] = E [ ξ ( t ) ] E\Big[\xi(t+\tau)\Big]=E\Big[\xi(t)\Big] E[ξ(t+τ)]=E[ξ(t)]

R ( ∞ ) = lim ⁡ τ → α E [ ξ ( t ) ξ ( t + τ ) ] = lim ⁡ τ → ∞ E [ ξ ( t ) ] E [ ξ ( t + τ ) ] = E [ ξ ( t ) ] E [ ξ ( t ) ] = E 2 [ ξ ( t ) ] R(\infty)=\lim_{\tau\to\alpha}E\left[\xi(t)\xi(t+\tau)\right]=\lim_{\tau\to\infty}E\left[\xi(t)\right]E\left[\xi(t+\tau)\right]=E\left[\xi(t)\right]E\left[\xi(t)\right]=E^2\left[\xi(t)\right] R()=limταE[ξ(t)ξ(t+τ)]=limτE[ξ(t)]E[ξ(t+τ)]=E[ξ(t)]E[ξ(t)]=E2[ξ(t)]

R ( 0 ) − R ( ∞ ) = E [ ξ 2 ( t ) ] − α 2 = σ 2 R(0)-R(\infty)=E\Bigl[\xi^2(t)\Bigr]-\alpha^2=\sigma^2 R(0)R()=E[ξ2(t)]α2=σ2 物理意义:随机过程的交流功率;对应方差公式 D [ ξ ( t ) ] = E [ ξ 2 ( t ) ] − E 2 [ ξ ( t ) ] D[\xi(t)]=E\Big[\xi^2(t)\Big]-E^2\Big[\xi(t)\Big] D[ξ(t)]=E[ξ2(t)]E2[ξ(t)]

R ( τ ) = R ( − τ ) R(\tau)=R(-\tau) R(τ)=R(τ) 自相关函数为偶函数;

∣ R ( τ ) ∣ ≤ R ( 0 ) \begin{vmatrix}R(\tau)\end{vmatrix}\leq R(0)\quad R(τ) R(0)给出了自相关函数的上界,与自身时刻相关性最大。

3.2.5平稳过程的功率谱密度

  • 定义

假定 f ( t ) f(t) f(t)为随机过程 ξ ( t ) \xi(t) ξ(t)的任一实现,对其进行 T T T 长度的截断,记

f r ( t ) f_r(t) fr(t),其傅里叶变换为 F T ( ω ) F_T(\omega) FT(ω),则任一实现的功率谱为:

P f ( ω ) = lim ⁡ T → ∞ ∣ F T ( ω ) ∣ 2 T P_f(\omega)=\lim_{T\to\infty}\frac{\left|F_T(\omega)\right|^2}T Pf(ω)=TlimTFT(ω)2

ξ ( t ) \xi(t) ξ(t) 的功率谱密度为:

  • 维纳-辛钦定理

平稳随机过程 ξ ( t ) \xi(t) ξ(t)功率谱密度函数 P ξ ( ω ) P_\xi(\omega) Pξ(ω)自相关函数 R ( τ ) R(\tau) R(τ)为一对傅里叶变换

{ P ξ ( ω ) = ∫ − ∞ + ∞ R ( τ ) e − j ω τ d τ R ( τ ) = 1 2 π ∫ − ∞ + ∞ P ξ ( ω ) e j ω τ d ω 或 { P ξ ( f ) = ∫ − ∞ + ∞ R ( τ ) e − j 2 π f τ d τ R ( τ ) = ∫ − ∞ + ∞ P ξ ( f ) e j 2 π f τ d f R ( τ ) ⇔ P ξ ( f ) \begin{cases}P_{\xi}(\omega)=\int_{-\infty}^{+\infty}R(\tau)e^{-j\omega\tau}\mathrm{d}\tau\\R(\tau)=\dfrac{1}{2\pi}\int_{-\infty}^{+\infty}P_{\xi}(\omega)e^{j\omega\tau}\mathrm{d}\omega\end{cases}\text{或}\begin{cases}P_{\xi}(f)=\int_{-\infty}^{+\infty}R(\tau)e^{-j2\pi f\tau}\mathrm{d}\tau\\R(\tau)=\int_{-\infty}^{+\infty}P_{\xi}(f)e^{j2\pi f\tau}\mathrm{d}f\end{cases}R(\tau)\Leftrightarrow P_{\xi}(f) Pξ(ω)=+R(τ)eτdτR(τ)=2π1+Pξ(ω)eτdω{Pξ(f)=+R(τ)ej2πfτdτR(τ)=+Pξ(f)ej2πfτdfR(τ)Pξ(f)

  • 功率谱密度的性质

功率谱密度具有非负性: P ξ ( f ) ≥ 0 P_\xi(f)\geq0 Pξ(f)0

功率谱密度是偶函数 : P ξ ( − f ) = P ξ ( f ) :P_\xi(-f)=P_\xi(f) :Pξ(f)=Pξ(f) ·

单边、双边功率谱密度互换: P ξ 单边 ( f ) = { 2 P ξ 双边 ( f ) f ≥ 0 0 f < 0 P_{\xi\text{单边}}(f)=\begin{cases}2P_{\xi\text{双边}}(f)&f\geq0\\0&f<0&\end{cases} Pξ单边(f)={2Pξ双边(f)0f0f<0

平均功率的计算方法
S = R ( 0 ) = E [ ξ 2 ( t ) ] S = ∫ − ∞ + ∞ P ξ ( f ) d f = 1 2 π ∫ − ∞ + ∞ P ξ ( ω ) d ω S=R(0)=E\bigg[\xi^2(t)\bigg]\quad S=\int_{-\infty}^{+\infty}P_\xi(f)\mathrm{d}f=\frac{1}{2\pi}\int_{-\infty}^{+\infty}P_\xi(\omega)\mathrm{d}\omega S=R(0)=E[ξ2(t)]S=+Pξ(f)df=2π1+Pξ(ω)dω

3.3高斯随机过程

3.3.1定义

若随机过程的任意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 exp ⁡ [ − 1 2 ∣ B ∣ ∑ j = 1 n ∑ k = 1 n ∣ B j k ∣ ( x j − a j σ j ) ( x k − a k σ k ) ] \begin{aligned}&f_{n}(x_{1},x_{2},\cdots,x_{n};t_{1},t_{2},\cdots,t_{n})\\&=\frac{1}{\left(2\pi\right)^{n/2}\sigma_{1}\sigma_{2}\cdots\sigma_{n}\left|B\right|^{1/2}}\exp\left[-\frac{1}{2\left|B\right|}\sum_{j=1}^{n}\sum_{k=1}^{n}\left|B_{jk}\right|\left(\frac{x_{j}-a_{j}}{\sigma_{j}}\right)\left(\frac{x_{k}-a_{k}}{\sigma_{k}}\right)\right]\end{aligned} fn(x1,x2,,xn;t1,t2,,tn)=(2π)n/2σ1σ2σnB1/21exp[2B1j=1nk=1nBjk(σjxjaj)(σkxkak)]

3.3.2重要性质

  • 高斯过程的n维分布完全由n个随机变量 x 1 , x 2 , ⋯   , x n x_1,x_2,\cdots,x_n x1,x2,,xn的数学期望,方差以及两两变量之间的归一化协方差函数决定;
  • 广义平稳的高斯过程也狭义平稳;
  • 高斯过程在不同时刻取值不相关,则它们也统计独立

f n ( x 1 , x 2 , ⋯   , x n ; t 1 , t 2 , ⋯   , t n ) = f ( x 1 , t 1 ) f ( x 2 , t 2 ) ⋯ f ( x n , t n ) f_n(x_1,x_2,\cdots,x_n;t_1,t_2,\cdots,t_n)=f(x_1,t_1)f(x_2,t_2)\cdots f(x_n,t_n) fn(x1,x2,,xn;t1,t2,,tn)=f(x1,t1)f(x2,t2)f(xn,tn)

对于高斯变量来说,不相关和独立是等价的

  • 高斯过程经过线性变换(或线性系统)后仍为高斯过程,但数字特征发生改变:
  • 若干个高斯过程的代数和仍为高斯过程,但数字特征发生改变。

X 1 ∼ N ( μ 1 , σ 1 2 ) , X 2 ∼ N ( μ 2 , σ 2 2 ) ⟶ a X 1 + b X 2 ∼ N ( a μ 1 + b μ 2 , a 2 σ 1 2 + b 2 σ 2 2 ) X_1\sim N(\mu_1,\sigma_1^2),X_2\sim N(\mu_2,\sigma_2^2)\boxed{\longrightarrow}aX_1+bX_2\sim N(a\mu_1+b\mu_2,a^2\sigma_1^2+b^2\sigma_2^2) X1N(μ1,σ12),X2N(μ2,σ22)aX1+bX2N(aμ1+bμ2,a2σ12+b2σ22)

3.3.3一维高斯过程

3.3.3.1概率密度函数与数字特征

f ( x ) = 1 2 π σ exp ⁡ [ − ( x − a ) 2 2 σ 2 ] , − ∞ < x < + ∞ ⟹ X ∼ N ( a , σ 2 ) f(x)=\frac1{\sqrt{2\pi}\sigma}\exp\biggl[-\frac{(x-a)^2}{2\sigma^2}\biggr],-\infty<x<+\infty\quad\boxed{\Longrightarrow}X\sim N({a},{\sigma^2}) f(x)=2π σ1exp[2σ2(xa)2],<x<+XN(a,σ2)

3.3.3.2概率密度曲线及其性质
  • 对称轴:均值 x = a x=a x=a

  • ∫ − ∞ + ∞ f ( x ) d x = 1 ∫ − ∞ a f ( x ) d x = ∫ a + ∞ f ( x ) d x = 1 2 \int_{-\infty}^{+\infty}f(x)\mathrm{d}x=1\int_{-\infty}^{a}f(x)\mathrm{d}x=\int_{a}^{+\infty}f(x)\mathrm{d}x=\frac12 +f(x)dx=1af(x)dx=a+f(x)dx=21

  • 均值 α \alpha α反映变量的分布中心, 方差 σ 2 \sigma^{2} σ2反映变量的集中程度;

  • 当均值 α = 0 \alpha=0 α=0,方差 σ 2 = 1 \sigma^2=1 σ2=1时,高斯分布为标准正态分布N(0,1)

image-20240425184319277

3.3.3.3 高斯变量的标准化/归一化

X ∼ N ( a , σ 2 ) ⟶ Y = X − a σ Y ∼ N ( 0 , 1 ) P ( X ≤ x ) = P ( Y ≤ x − a σ ) X\sim N(a,\sigma^2)\quad \overset{Y=\frac{X-a}\sigma}{\longrightarrow } \quad Y\sim N(0,1)\\P(X\leq x)=P(Y\leq\frac{x-a}\sigma) XN(a,σ2)Y=σXaYN(0,1)P(Xx)=P(Yσxa)

给定任意高斯变量,都可以经过归一化变为标准正态分布Y。

3.3.4 正态分布函数、Q函数、误差函数、互补误差函数

3.3.4.1 正态分布
  • 正态分布函数

X ∼ N ( a , σ 2 ) X\sim N(a,\sigma^2) XN(a,σ2) F ( x ) = P ( X ≤ x ) = ∫ − ∞ x 1 2 π σ exp ⁡ ∣ − ( t − a ) 2 2 σ 2 ∣ F( x) = P( X\leq x) = \int _{- \infty}^x\frac 1{\sqrt {2\pi}\sigma}\exp \left | - \frac {\left ( t- a\right ) ^2}{2\sigma^2}\right | F(x)=P(Xx)=x2π σ1exp 2σ2(ta)2 d t t t

  • 标准正态分布函数

F ( x ) = ∫ − ∞ x 1 2 π e − t 2 2 d t F(x)=\int_{-\infty}^x\frac{1}{\sqrt{2\pi}}e^{-\frac{t^2}{2}}\mathrm{d}t F(x)=x2π 1e2t2dt

积分式无闭合形式,只有数值解,需要经过处理以后通过查表求解概率。

3.3.4.2 Q函数

X ∼ N ( 0 , 1 ) X\sim N(0,1) XN(0,1)

Q ( x ) = ∫ x + ∞ 1 2 π e − t 2 2 d t Q(x)=\int_x^{+\infty}\frac{1}{\sqrt{2\pi}}e^{-\frac{t^2}{2}}\mathrm{d}t Q(x)=x+2π 1e2t2dt
X ∼ N ( a , σ 2 ) → Y = X − a σ Y ∼ N ( 0 , 1 ) X\sim N(a,\sigma^2)\quad\xrightarrow{\begin{array}{c}Y=\frac{X-a}\sigma\\\end{array}}\quad Y\sim N(0,1) XN(a,σ2)Y=σXa YN(0,1)

P ( X ≤ x ) = P ( Y ≤ x − a σ ) Q ( x ) = 1 − F ( x − a σ ) P(X\leq x)=P(Y\leq\frac{x-a}{\sigma})\quad\boxed{Q(x)=1-F\left(\frac{x-a}{\sigma}\right)} P(Xx)=P(Yσxa)Q(x)=1F(σxa)

给定任意高斯变量 X X X,经过归一化变为标准正态分布 Y Y Y后,可以利用Q函数查表求有关的概率。
Q ( 0 ) = 0.5 , Q ( x ) = 1 − Q ( − x ) Q(0)=0.5,\quad Q(x)=1-Q(-x) Q(0)=0.5,Q(x)=1Q(x)
image-20240426152530771

3.4平稳随机过程通过线性系统

随机机过程通过线性系统

只考虑平稳随机过程通过线性时不变系统的情形。

在这里插入图片描述

对确知信号(视为随机过程的样本函数):

v o ( t ) = v i ( t ) ∗ h ( t ) = ∫ 0 + ∞ h ( τ ) v i ( t − τ ) d τ V o ( ω ) V i ( ω ) H ( ω ) V o ( ω ) = V i ( ω ) ⋅ H ( ω ) \begin{aligned}&v_{o}(t)=v_{i}(t)*h(t)=\int_{0}^{+\infty}h(\tau)v_{i}(t-\tau)\mathrm{d}\tau\\&V_{o}(\omega)\quad V_{i}(\omega)\quad H(\omega)\quad V_{o}(\omega)=V_{i}(\omega)\cdot H(\omega)\end{aligned} vo(t)=vi(t)h(t)=0+h(τ)vi(tτ)dτVo(ω)Vi(ω)H(ω)Vo(ω)=Vi(ω)H(ω)
对整个随机过程: ξ o ( t ) = ξ i ( t ) ∗ h ( t ) = ∫ 0 + ∞ h ( τ ) ξ i ( t − τ ) d τ \xi_o(t)=\xi_i(t)*h(t)=\int_0^{+\infty}h(\tau)\xi_i(t-\tau)\mathrm{d}\tau ξo(t)=ξi(t)h(t)=0+h(τ)ξi(tτ)dτ

3.4.1 随机过程通过线性系统

  • 数学期望

E [ ξ i ( t ) ] = a ⟶ E [ ξ o ( t ) ] = a ⋅ H ( 0 ) E\left[\xi_i(t)\right]=a\quad\longrightarrow\quad E\left[\xi_o(t)\right]=a\cdot H(0) E[ξi(t)]=aE[ξo(t)]=aH(0)

输出过程的数学期望为常数,即输入直流分量a与系统直流增益H(0)的乘积。
R i ( t , t + τ ) = R i ( τ ) ⟶ R o ( t , t + τ ) = R o ( τ ) R_i(t,t+\tau)=R_i(\tau)\quad\longrightarrow\quad R_o(t,t+\tau)=R_o(\tau) Ri(t,t+τ)=Ri(τ)Ro(t,t+τ)=Ro(τ)
输入平稳随机过程,则输出过程自相关函数只与时间间隔有关。

平稳随机过程通过线性系统,输出也是平稳随机过程。

  • 证明

E [ ξ o ( t ) ] = a ⋅ H ( 0 ) E[\xi_o(t)]=a\cdot H(0) E[ξo(t)]=aH(0)

E [ ξ o ( t ) ] = E [ ∫ 0 + ∞ h ( τ ) ξ i ( t − τ ) d τ ] = ∫ 0 + ∞ h ( τ ) E [ ξ i ( t − τ ) ] d τ = E [ ξ i ( t ) ] ⋅ ∫ 0 + ∞ h ( τ ) d τ = a ⋅ ∫ 0 + ∞ h ( τ ) d τ ⟶ E [ ξ o ( t ) ] = a ⋅ H ( 0 ) \begin{aligned} E[\xi_{o}(t)]& =E\bigg[\int_{0}^{+\infty}h(\tau)\xi_{i}(t-\tau)\mathrm{d}\tau\bigg]= \int_{0}^{+\infty}h(\tau)E\bigg[\xi_{i}(t-\tau)\bigg]\mathrm{d}\tau \\ &=E\left[\xi_i(t)\right]\cdot\int_0^{+\infty}h(\tau)\mathrm{d}\tau \\ &=a\cdot\int_0^{+\infty}h(\tau)\mathrm{d}\tau\quad\longrightarrow\quad E\left[\xi_o(t)\right]=a\cdot H(0) \end{aligned} E[ξo(t)]=E[0+h(τ)ξi(tτ)dτ]=0+h(τ)E[ξi(tτ)]dτ=E[ξi(t)]0+h(τ)dτ=a0+h(τ)dτE[ξo(t)]=aH(0)

H ( f ) = ∫ 0 + ∞ h ( τ ) e − j 2 π f τ d τ H ( 0 ) = ∫ 0 + ∞ h ( τ ) e − j 2 π f τ ⏟ 1 d τ ∣ f = 0 = ∫ 0 + ∞ h ( τ ) d τ \begin{aligned}&H(f)=\int_{0}^{+\infty}h(\tau)e^{-j2\pi f\tau}\mathrm{d}\tau\\&H(0)=\int_{0}^{+\infty}h(\tau)\underbrace{e^{-j2\pi f\tau}}_{1}\mathrm{d}\tau\Bigg|_{f=0}=\int_{0}^{+\infty}h(\tau)\mathrm{d}\tau\end{aligned} H(f)=0+h(τ)ej2πfτdτH(0)=0+h(τ)1 ej2πfτdτ f=0=0+h(τ)dτ

  • 证明

R o ( t , t + τ ) = R o ( τ ) R_o(t,t+\tau)=R_o(\tau) Ro(t,t+τ)=Ro(τ)

R o ( t , t + τ ) = E [ ξ o ( t ) ξ o ( t + τ ) ] = E [ ∫ 0 + ∞ h ( α ) ξ i ( t − α ) d α ∫ 0 + ∞ h ( β ) ξ i ( t + τ − β ) d β ] = ∫ 0 + ∞ ∫ 0 + ∞ h ( α ) h ( β ) E [ ξ i ( t − α ) ξ i ( t + τ − β ) ] d α d β \begin{aligned} R_{o}(t,t+\tau)& =E\Big[\xi_{o}(t)\xi_{o}(t+\tau)\Big] \\ &=E\left[\int_0^{+\infty}h(\alpha)\xi_i(t-\alpha)\mathrm{d}\alpha\int_0^{+\infty}h(\beta)\xi_i(t+\tau-\beta)\mathrm{d}\beta\right] \\ &=\int_{0}^{+\infty}\int_{0}^{+\infty}h(\alpha)h(\beta)E\Big[\xi_{i}(t-\alpha)\xi_{i}(t+\tau-\beta)\Big]\mathrm{d}\alpha\mathrm{d}\beta \end{aligned} Ro(t,t+τ)=E[ξo(t)ξo(t+τ)]=E[0+h(α)ξi(tα)dα0+h(β)ξi(t+τβ)dβ]=0+0+h(α)h(β)E[ξi(tα)ξi(t+τβ)]dαdβ

E [ ξ i ( t − α ) ξ i ( t + τ − β ) ] = R i ( τ + α − β ) R o ( t , t + τ ) = ∫ 0 + ∞ ∫ 0 + ∞ h ( α ) h ( β ) R i ( τ + α − β ) d α d β = R o ( τ ) \begin{aligned} &E\Big[\xi_{i}(t-\alpha)\xi_{i}(t+\tau-\beta)\Big]=R_{i}(\tau+\alpha-\beta) \\ &R_{o}(t,t+\tau)=\int_{0}^{+\infty}\int_{0}^{+\infty}h(\alpha)h(\beta)R_{i}(\tau+\alpha-\beta)\mathrm{d}\alpha\mathrm{d}\beta=R_{o}(\tau) \end{aligned} E[ξi(tα)ξi(t+τβ)]=Ri(τ+αβ)Ro(t,t+τ)=0+0+h(α)h(β)Ri(τ+αβ)dαdβ=Ro(τ)

因此输出过程自相关函数只依赖时间间隔,而与时间起点无关。

  • 功率谱密度

P i ( ω ) ⟹ P 0 ( ω ) = P i ( ω ) ∣ H ( ω ) ∣ 2 P_i(\omega)\quad\Longrightarrow\quad P_0(\omega)=P_i(\omega)\left|H(\omega)\right|^2 Pi(ω)P0(ω)=Pi(ω)H(ω)2

输出过程的功率谱为输入功率谱 P i ( ω ) P_i(\omega) Pi(ω)与系统功率增益 ∣ H ( ω ) ∣ 2 \left|H(\omega)\right|^2 H(ω)2的乘积。

研究一般随机过程经过线性系统的输出较复杂,但特殊地,线性系统输入高斯过程,则输出也为高斯过程,但数字特征会发生改变

  • 随机过程通过线性系统

P 0 ( ω ) = P i ( ω ) ∣ H ( ω ) ∣ 2 P_0(\omega)=P_i(\omega){\left|H(\omega)\right|}^2 P0(ω)=Pi(ω)H(ω)2

R o ( t , t + τ ) = ∫ 0 + ∞ ∫ 0 + ∞ h ( α ) h ( β ) R i ( τ + α − β ) d α d β R_o(t,t+\tau)=\int_0^{+\infty}\int_0^{+\infty}h(\alpha)h(\beta)R_i(\tau+\alpha-\beta)\mathrm{d}\alpha\mathrm{d}\beta Ro(t,t+τ)=0+0+h(α)h(β)Ri(τ+αβ)dαdβ

P o ( ω ) = ∫ − ∞ + ∞ R o ( τ ) e − j ω τ d τ = ∫ − ∞ + ∞ ∫ 0 + ∞ ∫ 0 + ∞ [ h ( α ) h ( β ) R i ( τ + α − β ‾ ) d α d β ] e − j ω τ d τ τ ′ = τ + α − β τ = τ ′ − α + β , d τ = d τ ′ = ∫ 0 + ∞ h ( α ) e j ω α ‾ d α ⋅ ∫ 0 + ∞ h ( β ) e − j ω β d β ⋅ ∫ − ∞ + ∞ R i ( τ ′ ) e − j ω τ ′ d τ ′ = H ∗ ( ω ) ⋅ H ( ω ) ⋅ P i ( ω ) = P i ( ω ) ∣ H ( ω ) ∣ 2 \begin{aligned} P_{o}(\omega)& =\int_{-\infty}^{+\infty}R_{o}(\tau)e^{-j\omega\tau}\mathrm{d}\tau \\ &=\int_{-\infty}^{+\infty}\int_{0}^{+\infty}\int_{0}^{+\infty}\left[h(\alpha)h(\beta)R_{i}(\underline{\tau+\alpha-\beta})\mathrm{d}\alpha\mathrm{d}\beta\right]e^{-j\omega\tau}\mathrm{d}\tau \\ &\tau^{\prime}=\tau+\alpha-\beta\quad\tau=\tau^{\prime}-\alpha+\beta,\mathrm{d}\tau=\mathrm{d}\tau^{\prime} \\ &=\int_{0}^{+\infty}h(\alpha)e^{\underline{j\omega\alpha}}\mathrm{d}\alpha\cdot\int_{0}^{+\infty}h(\beta)e^{-j\omega\beta}\mathrm{d}\beta\cdot\int_{-\infty}^{+\infty}R_{i}(\tau^{\prime})e^{-j\omega\tau^{\prime}}\mathrm{d}\tau^{\prime} \\ &=H^{*}(\omega)\cdot H(\omega)\cdot P_{i}(\omega)=P_{i}(\omega)\big|H(\omega)\big|^{2} \end{aligned} Po(ω)=+Ro(τ)eτdτ=+0+0+[h(α)h(β)Ri(τ+αβ)dαdβ]eτdττ=τ+αβτ=τα+β,dτ=dτ=0+h(α)eαdα0+h(β)eβdβ+Ri(τ)eτdτ=H(ω)H(ω)Pi(ω)=Pi(ω) H(ω) 2

image-20240426162500843

3.5窄带随机过程

image-20240426172832422
P ξ ( f ) = { P ξ ( f ) ∣ f ∣ ≤ f 0 0 ∣ f ∣ > f 0 P ξ ( f ) = { P ξ ( f ) f c − f 0 ≤ ∣ f ∣ ≤ f c + f 0 0 其他 P_\xi(f)=\begin{cases}P_\xi(f)&\left|f\right|\leq f_0\\\\0&\left|f\right|>f_0\end{cases}\quad \quad\quad P_\xi(f)=\begin{cases}P_\xi(f)&f_c-f_0\leq\left|f\right|\leq f_c+f_0\\\\0&\text{其他}\end{cases} Pξ(f)= Pξ(f)0ff0f>f0Pξ(f)= Pξ(f)0fcf0ffc+f0其他

3.5.1 窄带随机过程的定义

通带宽度: Δ f ≪ f c \Delta f\ll f_c Δffc, 且中心频率 f c f_c fc 远离零频率的系统。

窄带过程:随机过程通过以 f c f_c fc为中心频率的窄带系统的输出过程。

  • 窄带过程时、频示意图

image-20240426180359350

窄带随机过程的数学表示
ξ ( t ) = a ξ ( t ) ‾ cos ⁡ [ ω c t + φ ξ ( t ) ‾ ] \xi(t)=\underline{a_\xi(t)}\cos\left[\omega_ct+\underline{\varphi_\xi(t)}\right] ξ(t)=aξ(t)cos[ωct+φξ(t)]
a ξ ( t ) {a_\xi(t)} aξ(t):随机包络

φ ξ ( t ) {\varphi_\xi(t)} φξ(t):随机相位

image-20240426180726416

  • 同相分量-正交分量表达式

ξ ( t ) = a ξ ( t ) cos ⁡ [ ω c t + φ ξ ( t ) ] = a ξ ( t ) cos ⁡ φ ξ ( t ) ‾ cos ⁡ ω c t − a ξ ( t ) sin ⁡ φ ξ ( t ) ‾ sin ⁡ ω c t = ξ c ( t ) ‾ cos ⁡ ω c t − ξ s ( t ) ‾ sin ⁡ ω c t \begin{aligned} \xi(t)& =a_{\xi}(t)\cos\Bigl[\omega_{c}t+\varphi_{\xi}(t)\Bigr] \\ &=\underline{a_\xi(t)\cos\varphi_\xi(t)}\cos\omega_ct-\underline{a_\xi(t)\sin\varphi_\xi(t)}\sin\omega_ct \\ &=\underline{\xi_c(t)}\cos\omega_ct-\underline{\xi_s(t)}\sin\omega_ct \end{aligned} ξ(t)=aξ(t)cos[ωct+φξ(t)]=aξ(t)cosφξ(t)cosωctaξ(t)sinφξ(t)sinωct=ξc(t)cosωctξs(t)sinωct

同相分量: ξ c ( t ) = a ξ ( t ) cos ⁡ φ ξ ( t ) \xi_c(t)=a_\xi(t)\cos\varphi_\xi(t) ξc(t)=aξ(t)cosφξ(t)

正交分量: ξ s ( t ) = a ξ ( t ) sin ⁡ φ ξ ( t ) \xi_s(t)=a_\xi(t)\sin\varphi_\xi(t) ξs(t)=aξ(t)sinφξ(t)
{ a ξ ( t ) = ξ c 2 ( t ) + ξ s 2 ( t ) φ ξ ( t ) = arctan ⁡ ξ s ( t ) ξ c ( t ) \begin{cases}a_\xi(t)=\sqrt{\xi_c^2(t)+\xi_s^2(t)}\\\varphi_\xi(t)=\arctan\frac{\xi_s(t)}{\xi_c(t)}\end{cases} {aξ(t)=ξc2(t)+ξs2(t) φξ(t)=arctanξc(t)ξs(t)

3.5.2 ξ c ( t ) 和 ξ s ( t ) \xi_{c}(t)\text{和}\xi_{s}(t) ξc(t)ξs(t)的统计特性

ξ ( t ) \xi(t) ξ(t) 为平稳高斯窄带过程,均值为0,方差为 σ ξ 2 \sigma_\mathrm{\xi}^2 σξ2

ξ ( 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 ) \xi_c(t) ξc(t) ξ s ( t ) \xi_s(t) ξs(t)数学期望均为0;

E [ ξ c ( t ) ] = E [ ξ s ( t ) ] = 0 E[\xi_c(t)]=E[\xi_s(t)]=0 E[ξc(t)]=E[ξs(t)]=0 ξ c ( t ) \xi_c(t) ξc(t) ξ s ( t ) \xi_s(t) ξs(t)相关函数均只与时间间隔有关,故必然平稳

ξ c ( t ) \xi_c(t) ξc(t) ξ s ( t ) \xi_s(t) ξs(t)自相关函数相同;

ξ c ( t ) \xi_c(t) ξc(t) ξ s ( t ) \xi_s(t) ξs(t)方差相同; σ ξ 2 = σ c 2 = σ s 2 \sigma_\xi^2=\sigma_c^2=\sigma_s^2 σξ2=σc2=σs2

ξ c ( t ) \xi_c(t) ξc(t) ξ s ( t ) \xi_s(t) ξs(t)互不相关统计独立高斯变量

ξ c ( t ) 和  ξ s ( t ) 为零均值,等方差,不相关,独立的平稳高斯过程。 \xi_c(t)\text{和 }\xi_s(t)\text{为零均值,等方差,不相关,独立的平稳高斯过程。} ξc(t) ξs(t)为零均值,等方差,不相关,独立的平稳高斯过程。

3.5.3 a ξ ( t )  和  φ ξ ( t ) 的统计特性 a_\xi(t)\text{ 和 }\varphi_\xi(t)\text{的统计特性} aξ(t)  φξ(t)的统计特性的统计特性

  • 瑞利分布

f ( a ξ ) = a ξ σ ξ 2 e x p ( − a ξ 2 2 σ ξ 2 ) , a ξ ≥ 0 f(a_\xi)=\frac{a_\xi}{\sigma_\xi^2}\mathrm{exp}\biggl(-\frac{a_\xi^2}{2\sigma_\xi^2}\biggr),a_\xi\geq0 f(aξ)=σξ2aξexp(2σξ2aξ2),aξ0

  • 均匀分布

f ( φ ξ ) = 1 2 π , 0 ≤ φ ξ ≤ 2 π f(\varphi_\xi)=\frac1{2\pi},0\leq\varphi_\xi\leq2\pi f(φξ)=2π1,0φξ2π

  • 随机相位和随机包络统计独立

f ( a ξ , φ ξ ) = f ( a ξ ) ⋅ f ( φ ξ ) f(a_\xi,\varphi_\xi)=f(a_\xi)\cdot f(\varphi_\xi) f(aξ,φξ)=f(aξ)f(φξ)

3.6正弦波加窄带高斯噪声

3.6.1 数学模型

  • 信号模型

r ( t ) = A cos ⁡ ( ω c t + θ ‾ ) + n ( t ) ‾ r(t)=A\cos(\omega_ct+\underline{\theta})+\underline{n(t)} r(t)=Acos(ωct+θ)+n(t)

θ \theta θ:随机相位

n ( t ) {n(t)} n(t):均值为0,方差为c的窄带高斯噪声

  • 同相分量-正交分量表达式

n ( t ) = n c ( t ) cos ⁡ ω c t − n s ( t ) sin ⁡ ω c t r ( t ) = A cos ⁡ ( ω c t + θ ) + n c ( t ) cos ⁡ ω c t − n s ( t ) sin ⁡ ω c 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 \begin{aligned} &n(t) =n_c(t)\cos\omega_ct-n_s(t)\sin\omega_ct \\ &r(t) =A\cos(\omega_ct+\theta)+n_c(t)\cos\omega_ct-n_s(t)\sin\omega_ct \\ &=\underline{\begin{bmatrix}A\cos\theta+n_c(t)\end{bmatrix}}\cos\omega_ct-\underline{\begin{bmatrix}A\sin\theta+n_s(t)\end{bmatrix}}\sin\omega_ct \\ &=\underline{z_c(t)}\cos\omega_ct-\underline{z_s(t)}\sin\omega_ct \end{aligned} n(t)=nc(t)cosωctns(t)sinωctr(t)=Acos(ωct+θ)+nc(t)cosωctns(t)sinωct=[Acosθ+nc(t)]cosωct[Asinθ+ns(t)]sinωct=zc(t)cosωctzs(t)sinωct

同相分量: z c ( t ) = A cos ⁡ θ + n c ( t ) z_c(t)=A\cos\theta+n_c(t) zc(t)=Acosθ+nc(t)

正交分量: z s ( t ) = A sin ⁡ θ + n s ( t ) z_s(t)=A\sin\theta+n_s(t) zs(t)=Asinθ+ns(t)

  • 幅度-相位表达式

r ( t ) = z c ( t ) cos ⁡ ω c t − z s ( t ) sin ⁡ ω c t = z ( t ) ‾ cos ⁡ [ ω c t + φ ( t ) ‾ ] r(t)=z_c(t)\cos\omega_ct-z_s(t)\sin\omega_ct=\underline{z(t)}\cos{}[ \omega_ct+\underline{\varphi(t)} ] r(t)=zc(t)cosωctzs(t)sinωct=z(t)cos[ωct+φ(t)]

随机包络 z ( t ) = z c 2 ( t ) + z s 2 ( t ) \text{随机包络}\quad z(t)=\sqrt{z_c^2(t)+z_s^2(t)}\quad 随机包络z(t)=zc2(t)+zs2(t)

随机相位 φ ( t ) = arctan ⁡ z s ( t ) z c ( t ) \text{随机相位}\quad\varphi(t)=\arctan\frac{z_s(t)}{z_c(t)} 随机相位φ(t)=arctanzc(t)zs(t)

  • 同相分量-正交分量与幅度-相位的关系

image-20240426183523682

3.6.2 同相分量和正交分量和统计特性

z c ( t ) = A cos ⁡ θ + n c ( t ) z_c(t)=A\cos\theta+\boxed{n_c(t)} zc(t)=Acosθ+nc(t)

z s ( t ) = A sin ⁡ θ + n s ( t ) z_s(t)=A\sin\theta+\boxed{n_s(t)} zs(t)=Asinθ+ns(t)

高斯窄带过程的同相分量与正交分量零均值,方差相同的高斯变量

 

随机相位 θ \theta θ给定时,$z_c( t) $ 和 $z_s( t) $为取值不相关,相互独立,且 同分布的高斯变量,其数字特征为:
{ E [ z c ( t ) ] = A cos ⁡ θ E [ z s ( t ) ] = A sin ⁡ θ σ n 2 = σ c 2 = σ s 2 \begin{cases}E[z_c(t)]=A\cos\theta\\E[z_s(t)]=A\sin\theta\\\sigma_n^2=\sigma_c^2=\sigma_s^2\end{cases} E[zc(t)]=AcosθE[zs(t)]=Asinθσn2=σc2=σs2

3.6.3 随机包络和随机相位的统计特性

随机包络 z ( t ) 的统计特性 \text{随机包络}z(t)\text{的统计特性} 随机包络z(t)的统计特性

莱斯分布/广义瑞利分布
f ( z ) = z σ n 2 exp ⁡ ( − z 2 + A 2 2 σ n 2 ) I 0 ( A z σ n 2 ) , z ≥ 0 f(z)=\frac{z}{\sigma_n^2}\exp\biggl(-\frac{z^2+A^2}{2\sigma_n^2}\biggr)I_0\biggl(\frac{Az}{\sigma_n^2}\biggr),z\geq0 f(z)=σn2zexp(2σn2z2+A2)I0(σn2Az),z0

I 0 ( x ) = 1 2 π ∫ 0 2 π exp ⁡ [ x cos ⁡ φ ] d φ 第一类零阶修正贝塞尔函数 I_0(x)=\frac1{2\pi}\int_0^{2\pi}\exp[x\cos\varphi]\mathrm{d}\varphi\quad\text{第一类零阶修正贝塞尔函数} I0(x)=2π102πexp[xcosφ]dφ第一类零阶修正贝塞尔函数
大信噪比环境: r = A 2 2 σ 2 ≫ 1 r=\frac{A^2}{2\sigma^2}\gg1 r=2σ2A21 近似为高斯分布

小信噪比环境: A → 0 A\to0 A0 退化为瑞利分布

随机相位 φ ( t ) \varphi(t) φ(t)的概率密度 f ( φ ) f(\varphi) f(φ)与信噪比环境有关。

  • 小信噪比环境: f ( φ ) f(\varphi) f(φ) 服从均匀分布
  • 大信噪比环境:集中在有用信号相位附近,在0相位附近取值集中。

image-20240426184421510
f ( a ξ ) = a ξ σ ξ 2 exp ⁡ ( − a ξ 2 2 σ ξ 2 ) , a ξ ≥ 0 f(a_\xi)=\frac{a_\xi}{\sigma_\xi^2}\exp\biggl(-\frac{a_\xi^2}{2\sigma_\xi^2}\biggr),a_\xi\geq0 f(aξ)=σξ2aξexp(2σξ2aξ2),aξ0
退化后
f ( z ) = z σ n 2 exp ⁡ ( − z 2 + A 2 2 σ n 2 ) I 0 ( A z σ n 2 ) , z ≥ 0 f ( φ ξ ) = 1 2 π , 0 ≤ φ ξ ≤ 2 π f(z)=\frac{z}{\sigma_{n}^{2}}\exp\biggl(-\frac{z^{2}+A^{2}}{2\sigma_{n}^{2}}\biggr)I_{0}\biggl(\frac{Az}{\sigma_{n}^{2}}\biggr),z\geq0\\f(\varphi_{\xi})=\frac{1}{2\pi},0\leq\varphi_{\xi}\leq2\pi f(z)=σn2zexp(2σn2z2+A2)I0(σn2Az),z0f(φξ)=2π1,0φξ2π
窄带过程”和“正弦波+高斯窄带过程”两类模型的有关结论常用于分析系统抗噪声性能
结合误差函数、互补误差函数的定义可以求解系统的总误码率。

3.7高斯白噪声和带限白噪声

3.7.1 高斯白噪声

  • 白噪声

功率谱密度均匀分布在整个频率范围内的随机过程称为白噪声。

  • 功率谱密度

P n ( ω ) = n 0 2 , − ∞ < ω < + ∞ P n ( ω ) = n 0 , ω ≥ 0 P_n(\omega)=\frac{n_0}2,-\infty<\omega<+\infty\qquad P_n(\omega)=n_0,\omega\geq0 Pn(ω)=2n0,<ω<+Pn(ω)=n0,ω0

n 0  ⁣ : n_0\colon n0: 单边噪声功率谱密度

n 0 2 : \frac{n_0}2: 2n0: 双边噪声功率谱密度 2

  • 自相关函数

P n ( ω ) = n 0 2 , − ∞ < ω < + ∞ ⟺ R ( τ ) = n 0 2 δ ( τ ) P_{n}(\omega)=\frac{n_{0}}{2},-\infty<\omega<+\infty\quad\Longleftrightarrow\quad R(\tau)=\frac{n_{0}}{2}\delta(\tau) Pn(ω)=2n0,<ω<+R(τ)=2n0δ(τ)

在这里插入图片描述

白噪声在任意两个不同时刻上取值互不相关,只有在间隔 τ = 0 \tau=0 τ=0 时才相关。

  • 高斯白噪声

服从高斯分布白噪声称为高斯白噪声。

  • 性质

高斯白噪声的特点是:在任意两个不同时刻上随机变量互不相关且统计独立

  • 实际应用

实际中不存在理想白噪声,因为其功率为无穷大,若噪声功率谱均匀分布的频率范围远大于通信系统的工作频带,则视为白噪声。

3.7.2 低通白噪声

  • 功率谱密度和自相关函数

P n ( ω ) = { n 0 2 ∣ ω ∣ ≤ ω H 0 ∣ ω ∣ > ω H ⟶ R ( τ ) = n 0 f H sin ⁡ ω H τ ω H τ P_n(\omega)=\begin{cases}\frac{n_0}{2}&\left|\omega\right|\leq\omega_H\\\\0&\left|\omega\right|>\omega_H\end{cases}\quad{\longrightarrow}\quad R(\tau)=n_0f_H\frac{\sin\omega_H\tau}{\omega_H\tau} Pn(ω)= 2n00ωωHω>ωHR(τ)=n0fHωHτsinωHτ
image-20240426185419742

  • 平均功率

S = R ( 0 ) = n 0 f H S = ∫ − ∞ + ∞ P n ( f ) d f = 1 2 π ∫ − ∞ + ∞ P n ( ω ) d ω = 1 2 π ∫ − ω H ω H n 0 2 d ω = n 0 f H \begin{aligned}&S=R(0)=n_0f_H\\&S=\int_{-\infty}^{+\infty}P_n(f)\mathrm{d}f=\frac1{2\pi}\int_{-\infty}^{+\infty}P_n(\omega)\mathrm{d}\omega\\&=\frac1{2\pi}\int_{-\omega_H}^{\omega_H}\frac{n_0}2\mathrm{d}\omega=n_0f_H\end{aligned} S=R(0)=n0fHS=+Pn(f)df=2π1+Pn(ω)dω=2π1ωHωH2n0dω=n0fH

对低通白噪声按照抽样定理抽样时,各抽样值为互不相关的随机变量。

3.7.3 带通白噪声

  • 功率谱密度和自相关函数

P n ( f ) = { n 0 2 f c − B 2 ≤ ∣ f ∣ ≤ f c + B 2 0 其他 ⟺ R ( τ ) = n 0 B S a ( π B τ ) cos ⁡ 2 π f c τ P_n(f)=\begin{cases}\frac{n_0}{2}&f_c-\frac{B}{2}\leq\left|f\right|\leq f_c+\frac{B}{2}\\0&\text{其他}\end{cases}\Longleftrightarrow\quad R(\tau)=n_0BSa(\pi B\tau)\cos2\pi f_c\tau Pn(f)={2n00fc2Bffc+2B其他R(τ)=n0BSa(πBτ)cos2πfcτ

image-20240426185842841

  • 平均功率

S = R ( 0 ) = n 0 B S = ∫ − ∞ + ∞ P n ( f ) d f = 2 B ⋅ n 0 2 = n 0 B \begin{aligned}&S=R(0)=n_0B\\&S=\int_{-\infty}^{+\infty}P_{n}(f)\mathrm{d}f=2B\cdot\frac{n_{0}}{2}=n_{0}B\end{aligned} S=R(0)=n0BS=+Pn(f)df=2B2n0=n0B

  • 低通白噪声

P n ( ω ) = { n 0 2 ∣ ω ∣ ≤ ω H 0 ∣ ω ∣ > ω H ⟺ R ( τ ) = n 0 f H sin ⁡ ω H τ ω H τ P_n(\omega)=\begin{cases}\frac{n_0}{2}&\left|\omega\right|\leq\omega_H\\0&\left|\omega\right|>\omega_H\end{cases}\quad\Longleftrightarrow\quad R(\tau)=n_0f_H\frac{\sin\omega_H\tau}{\omega_H\tau} Pn(ω)={2n00ωωHω>ωHR(τ)=n0fHωHτsinωHτ

  • 带通白噪声

P n ( f ) = { n 0 2 f c − B 2 ≤ ∣ f ∣ ≤ f c + B 2 0 其他 ⟺ R ( τ ) = n 0 B S a ( π B τ ) cos ⁡ 2 π f c τ P_n(f)=\begin{cases}\frac{n_0}{2}&f_c-\frac{B}{2}\leq\left|f\right|\leq f_c+\frac{B}{2}\\0&\text{其他}\end{cases}\quad\Longleftrightarrow\quad R(\tau)=n_0BSa(\pi B\tau)\cos2\pi f_c\tau Pn(f)={2n00fc2Bffc+2B其他R(τ)=n0BSa(πBτ)cos2πfcτ

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mfbz.cn/a/576783.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

使用composer开发自己的扩展包

前言 日常的开发中我们经常用到composer去安装其他人封装好的扩展包&#xff0c;如果你有好的功能代码想分享给其他人使用&#xff0c;就可以使用composer打包成扩展包。其他人用composer安装后就可以使用你的扩展包了。这篇文章教你如何打包自己的composer扩展包。 1.新建仓…

d10(104-112)-勇敢开始Java,咖啡拯救人生

目录 多态[面向对象的三大特征&#xff1a;封装 继承 多态]polymorphism 多态的好处 多态下的类型转换问题 注意事项 final关键字 注意事项 常量 抽象类abstract 好处 运用场景 接口interface 好处 接口中新增的三种方法&#xff1a; 接口的多继承 内部类[类中五…

江苏开放大学2024年春《机电设备安装与调试 050095》第三次形成性考核作业参考答案

电大搜题 多的用不完的题库&#xff0c;支持文字、图片搜题&#xff0c;包含国家开放大学、广东开放大学、超星等等多个平台题库&#xff0c;考试作业必备神器。 公众号 答案&#xff1a;更多答案&#xff0c;请关注【电大搜题】微信公众号 答案&#xff1a;更多答案&#…

boot2docker实践@FreeBSD

boot2docker是专为Docker而设计的轻量极Linux发型包&#xff0c;解决Windows、OS X、FreeBSD不能安装Docker的问题。Boot2Docker完全运行于内存中&#xff0c;24M大小&#xff0c;启动仅5-6秒。Boot2Docker需要运行在VirtualBox中&#xff0c;同时意味着有VirtualBox就能运行Bo…

万业企业发布23年年报、24年一季报,集成电路业务同比大增近七成 转型成效显著

4月26日晚间&#xff0c;万业企业&#xff08;600641&#xff09;发布2023年年度报告及2024年一季度报告。2023年&#xff0c;公司实现营业收入9.65亿元&#xff0c;归母净利润1.51亿元。其中&#xff0c;公司集成电路设备制造业务收入较上年同期大幅增长67.53%&#xff0c;公司…

CYCLE:学习自我完善代码生成

目录 IntriductionOverview of the Approach 预训练的代码语言模型在代码生成方面取得了可喜的性能&#xff0c;并提高了人类开发人员的编程效率。然而&#xff0c;现有的代码 LM 评估通常忽略了它们的 自我求精能力&#xff0c;这些评估仅关注一次性预测的准确性。对于代码 L…

Leetcode - 周赛394

目录 一&#xff0c;3120. 统计特殊字母的数量 I 二&#xff0c;3121. 统计特殊字母的数量 II 三&#xff0c;3122. 使矩阵满足条件的最少操作次数 四&#xff0c;3123. 最短路径中的边 一&#xff0c;3120. 统计特殊字母的数量 I 本题就是统计有多少个字母的大小写同时出现…

AIGC学习步骤

目录 AIGC学习步骤 步骤一&#xff1a;理解基本概念 步骤二&#xff1a;学习资源 步骤三&#xff1a;深入研究 步骤四&#xff1a;联系专家 步骤五&#xff1a;实践应用 步骤六&#xff1a;持续学习 AIGC学习步骤 我们先来说说什么是AIGC&#xff1f; 生成式人工智能—…

第二篇:Python环境搭建:从初学者到专家

Python环境搭建&#xff1a;从初学者到专家 在编程的世界里&#xff0c;准备好一个高效而舒适的开发环境是走向成功的第一步。在这篇博客文章中&#xff0c;我们将一起探索如何为Python编程搭建一个理想的环境。无论你是完全的新手还是希望提升现有的技能&#xff0c;本文都会…

IDEA快速入门

目录 1. 概述 2. 安装 3. 激活 4. 关闭自动更新 5. 创建Java项目 5.1 配置JRE 5.2 创建项目 6. 配置设置 6.1 主题 6.2 设置字体默认大小 6.3 鼠标滚轮改变字体大小 6.4 设置自动导入 6.5 项目选择 7. lombok插件 7.1 安装插件 7.2 启用注解 8. 安装包及插件…

SecuPress Pro 专业级WordPress网站安全防护插件优化版

下载地址&#xff1a;SecuPress Pro 专业版.zip SecuPress Pro&#xff1a;专业的WordPress安全解决方案 如果您没有时间进行每周扫描&#xff0c;SecuPress Pro将是您的理想选择。SecuPress Pro提供了所有SecuPress Free的功能&#xff0c;同时还增加了一些高级选项&#xff…

优维全新力作:统一采控平台

在本月&#xff0c;优维新一代核心系统「EasyOps」7.0大版本重磅上线&#xff0c;为广大用户带来了“更核心、更智能、更开放、更客制”的产品能力。&#xff08;点击回看&#xff1a;重磅&#xff01;优维科技发布EasyOps7.0大版本&#xff09;在本次版本能力分享上&#xff0…

VBA技术资料MF145:清空回收站

我给VBA的定义&#xff1a;VBA是个人小型自动化处理的有效工具。利用好了&#xff0c;可以大大提高自己的工作效率&#xff0c;而且可以提高数据的准确度。“VBA语言専攻”提供的教程一共九套&#xff0c;分为初级、中级、高级三大部分&#xff0c;教程是对VBA的系统讲解&#…

各种CSS导航代码集合!CV即可!水平导航,垂直导航,jd粘性导航…

你好&#xff0c;我是云桃桃。 一个希望帮助更多朋友快速入门 WEB 前端的程序媛。 云桃桃-大专生&#xff0c;一枚程序媛&#xff0c;感谢关注。回复 “前端基础题”&#xff0c;可免费获得前端基础 100 题汇总&#xff0c;回复 “前端工具”&#xff0c;可获取 Web 开发工具合…

mysql的多表查询和子查询

多表查询&#xff1a;查询数据时&#xff0c;需要使用多张表来查询 多表查询分类&#xff1a; 1.内连接查询 2.外连接查询 3.子查询 笛卡尔积&#xff1a; create table class (id int primary key auto_increment,name varchar(10) ); create table student (id int primar…

【TCP:可靠数据传输,快速重传,流量控制,TCP流量控制】

文章目录 可靠数据传输TCP&#xff1a;可靠数据传输TCP发送方事件快速重传流量控制TCP流量控制 可靠数据传输 TCP&#xff1a;可靠数据传输 TCP在IP不可靠服务的基础上建立了rdt 管道化的报文段 GBN or SR 累计确认&#xff08;像GBN&#xff09;单个重传定时器&#xff08;像…

搞懂特殊的引用类型接口、枚举、注解、记录、包装类!

面向对象系列——特殊的引用类型&#xff08;与class并列地位&#xff09; 文章目录 面向对象系列——特殊的引用类型&#xff08;与class并列地位&#xff09;特殊的&#xff1a;接口、枚举、注解、记录&#xff0c;包装类一、interface接口(不是类)—提供一种规范(具有多态性…

ionic 中对Input输入框、searchbar进行solr检索

一、概述 Ionic 是一个用于开发跨平台应用程序的开源工具&#xff0c;可以使用 Angular、React 或 Vue 等前端框架。要在 Ionic 应用程序中实现实时与 Solr 通信&#xff0c;可以使用 HTTP 客户端&#xff08;如 Angular 的 HttpClient 或 Ionic 的 Native HTTP&#xff09;…

MySQL—MySQL的存储引擎之InnoDB

MySQL—MySQL的存储引擎之InnoDB 存储引擎及种类 存储引擎说明MyISAM高速引擎&#xff0c;拥有较高的插入&#xff0c;查询速度&#xff0c;但不支持事务InnoDB5.5版本后MySQL的默认数据库存储引擎&#xff0c;支持事务和行级锁&#xff0c;比MyISAM处理速度稍慢ISAMMyISAM的…

一分钟教你下载小程序视频到手机#下载高手

本文就教你们如何下载小程序的视频到手机 首先将小程序视频下载到电脑上&#xff0c;然后再传送到手机上,这样就实现了下载小程序到手机上 可是怎么将小程序下载到电脑上呢&#xff1f;这里要借用一个工具&#xff1a;下载高手 下载高手工具我已经打包好了&#xff0c;有需要…
最新文章