通信原理(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.rhkb.cn/news/316658.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

http1.1和http2.0的同源请求数限制

判断协议版本 :scheme: 在请求头中表示使用的是HTTP/2协议。即 出现 :开头的请求头Chrome 只支持查看 HTTP/1.x 的 Raw Headers&#xff0c;对这种请求&#xff0c;会给出 view source 选项。HTTP2.0不给出。可继续学习 https://www.cnblogs.com/kirito-c/p/10360868.html抓包…

笔记:编写程序,分别采用面向对象和 pyplot 快捷函数的方式绘制正弦曲线 和余弦曲线。 提示:使用 sin()或 cos()函数生成正弦值或余弦值。

文章目录 前言一、面向对象和 pyplot 快捷函数的方式是什么&#xff1f;二、编写代码面向对象的方法&#xff1a;使用 pyplot 快捷函数的方法&#xff1a; 总结 前言 本文将探讨如何使用编程语言编写程序&#xff0c;通过两种不同的方法绘制正弦曲线和余弦曲线。我们将分别采用…

DSP开发实战教程--EPWM模块的影子寄存器详细讲解原理和代码实例

EPWM模块影子寄存器的原理 在TI&#xff08;Texas Instruments&#xff09;的DSP28335中&#xff0c;EPWM&#xff08;Enhanced Pulse Width Modulator&#xff09;模块提供了高精度、高灵活性的PWM信号生成功能。为了能在不影响当前PWM波形输出的情况下预装载新的PWM参数&…

cocos-lua资源管理

本文介绍cocos-lua项目的资源管理和工作流&#xff0c;适用人群包括初学者和有经验开发者&#xff0c;故读者可根据自己的需要有选择性的查阅自己需要的内容&#xff0c;下文以ccs代指Cocos Studio 一.简单案例解析 下文通过介绍一个简单demo&#xff0c;介绍合图和资源目录结…

[C++][算法基础]最大不相交区间数量(贪心 + 区间问题2)

给定 &#x1d441; 个闭区间 [&#x1d44e;&#x1d456;,&#x1d44f;&#x1d456;]&#xff0c;请你在数轴上选择若干区间&#xff0c;使得选中的区间之间互不相交&#xff08;包括端点&#xff09;。 输出可选取区间的最大数量。 输入格式 第一行包含整数 &#x1d4…

【算法基础实验】图论-UnionFind连通性检测之quick-union

Union-Find连通性检测之quick-union 理论基础 在图论和计算机科学中&#xff0c;Union-Find 或并查集是一种用于处理一组元素分成的多个不相交集合&#xff08;即连通分量&#xff09;的情况&#xff0c;并能快速回答这组元素中任意两个元素是否在同一集合中的问题。Union-Fi…

PHP源码_最新在线工具箱网站系统源码

项目运行截图 源码贡献 https://githubs.xyz/boot?app41 部分数据库表 SET NAMES utf8mb4; SET FOREIGN_KEY_CHECKS 0;-- ---------------------------- -- Table structure for toolbox_category -- ---------------------------- DROP TABLE IF EXISTS toolbox_category…

StarRocks x Paimon 构建极速实时湖仓分析架构实践

Paimon 介绍 Apache Paimon 是新一代的湖格式&#xff0c;可以使用 Flink 和 Spark 构建实时 Lakehouse 架构&#xff0c;以进行流式处理和批处理操作。Paimon 创新性地使用 LSM&#xff08;日志结构合并树&#xff09;结构&#xff0c;将实时流式更新引入 Lakehouse 架构中。 …

医学vr虚拟仿真综合实验教学平台为科研教学提供了坚实的基础

在兽医专业的广袤领域中&#xff0c;动物解剖学作为基石学科&#xff0c;为组织胚胎学、生理学、病理解剖学、外科手术学、临床诊断学等科研教学提供了坚实的基础。而如今&#xff0c;随着科技的飞速发展&#xff0c;我们迎来了一个全新的学习时代——3D数字动物解刨虚拟仿真实…

[iOS]使用CocoaPods发布公开库

1.检查库名是否已被占用 选择库名时&#xff0c;尽量选择具有描述性并且独特的名字&#xff0c;这不仅可以避免命名冲突&#xff0c;还可以帮助用户更好地理解库的用途和功能。 在实际创建和发布 CocoaPods 库之前&#xff0c;确实应该检查库名是否已经被占用&#xff0c;以避…

AutoCAD 2025 for mac/win:设计未来,触手可及

在数字化时代&#xff0c;设计不再局限于纸笔之间&#xff0c;而是跃然于屏幕之上&#xff0c;AutoCAD 2025正是这一变革的杰出代表。无论是Mac用户还是Windows用户&#xff0c;AutoCAD 2025都以其卓越的性能和出色的用户体验&#xff0c;成为了CAD设计绘图领域的佼佼者。 Aut…

Vuforia AR篇(三)— AR模型出场效果

目录 前言一、AR模型出场二、AR出场特效三、添加过渡效果四、效果 前言 在这个数字化日益增长的时代&#xff0c;增强现实&#xff08;AR&#xff09;技术正以前所未有的速度发展。AR模型&#xff0c;作为这一技术的核心组成部分&#xff0c;不仅改变了我们与数字世界的互动方…

【MATLAB源码-第201期】基于matlab的黏菌群优化算法(SMA)无人机三维路径规划,输出做短路径图和适应度曲线

操作环境&#xff1a; MATLAB 2022a 1、算法描述 黏菌优化算法&#xff08;Slime Mould Algorithm, SMA&#xff09;是一种新颖的启发式优化方法&#xff0c;其灵感来源于自然界中的真菌——黏菌。这种算法模拟了黏菌在寻找食物时的行为和网络形成策略。在本文中&#xff0c…

Python | Leetcode Python题解之第58题最后一个单词的长度

题目&#xff1a; 题解&#xff1a; class Solution:def lengthOfLastWord(self, s: str) -> int:ls[]for i in s.split():ls.append(i)return len(ls[-1])

【华为】NAT的分类和实验配置

【华为】NAT的分类和实验配置 NAT产生的技术背景IP地址分类NAT技术原理NAT分类静态NAT动态NATNAPTEasy IP&#xff08;PAT&#xff09;NAT Server 配置拓扑静态NAT测试抓包 动态NAT测试抓包 NAPT测试抓包 PAT测试抓包 NAT Server检测抓包 PC1PC2服务器 NAT产生的技术背景 随着…

排序算法:插入、希尔、选择、推排、冒泡、快速、归并排序

排序算法 目录 前言 一、排序的概念 1.1排序的概念 1.2 常见的排序算法 二、常见排序算法的实现 2.1 插入排序 2.2 希尔排序 2.3 选择排序 2.4 堆排序 2.5 冒泡排序 2.6 快速排序 2.6.1 hoare版本 2.6.2 前后指针版本 2.6.3 非递归版本 2.7 归并排序 归并排序 2.8 计数排序 三、…

Android视角看鸿蒙第十二课-鸿蒙的布局之相对布局RelativeContainer

Android视角看鸿蒙第十二课-鸿蒙的布局之相对布局RelativeContainer 导读 相对布局和线性、层叠布局一样都是类似于Android布局的&#xff0c;之前两篇文章已经了解线性、层叠布局的使用方法&#xff0c;这篇文章一起来学习下鸿蒙中的相对布局。 之前的文章中&#xff0c;我偶…

10.MMD 室内场景导入背景视频和灯光

导入背景视频 1. 导入人物和场景 场景是Akali’s room&#xff0c;可以在墙壁上添加视频 先添加主场景 2. 修改视频文件格式 在背景里选择导入背景视频文件 需要将mp4视频格式转化为AVI格式 方法一 先将视频导入格式工厂 点击配置 将视频编码改成DivX 再开始处理 …

数据结构八:线性表之循环队列的设计

上篇博客&#xff0c;学习了栈&#xff0c;我们可以知道他也是一种线性表&#xff0c;遵从先进后出的原则&#xff0c;在本节&#xff0c;我们进一步学习另一种线性表—队列。就像饭堂里排队打饭的的队伍&#xff0c;作为一种先进先出的线性表&#xff0c;他又有哪些特别之处呢…

敏捷之Scrum开发

目录 一、什么是 Scrum 1.1 Scrum 的定义 二、Scrum 迭代开发过程 2.1 迭代开发过程说明 2.1.1 开发方法 2.1.1.1 增量模型 2.1.1.1.1 定义 2.1.1.1.2 模型方法说明 2.1.1.2 迭代模型 2.1.1.2.1 定义 2.1.1.2.2 模型方法说明 2.1.2 迭代过程 2.1.2.1 产品需求Produ…