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

news2024/11/23 21:11:05

通信原理(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.coloradmin.cn/o/1627160.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

python-opencv实现最近邻插值和双线性插值对图片上采样

使用背景 当我们需要把图像进行放大或者缩小的时候&#xff0c;第一反应是使用resize()实现。很多情况下&#xff0c;我们会调用最近邻插值和双线性插值去放大图片&#xff0c;当然要说没有分辨率的损失那是不可能的&#xff0c;只能说在放大图片的过程中尽可能增加了图片的分…

Hybrid Homomorphic Encryption:SE + HE

参考文献&#xff1a; [NLV11] Naehrig M, Lauter K, Vaikuntanathan V. Can homomorphic encryption be practical?[C]//Proceedings of the 3rd ACM workshop on Cloud computing security workshop. 2011: 113-124.[MJS16] Maux P, Journault A, Standaert F X, et al. To…

Matlab|交直流系统潮流计算(含5种控制模式)

目录 1 主要内容 程序参考流程图 2 部分代码 3 程序结果 4 下载链接 1 主要内容 该程序参考文献《交直流系统潮流计算及相互关联特性分析》&#xff0c;采用5种交直流潮流控制方式&#xff1a;1.定电流定电压 2.定电流定熄弧角 3.定功率定电压 4.定功率定熄弧角 5.定触发角…

Stable Diffusion 常用放大算法详解

常用放大算法 图像放大算法大致有两种&#xff1a; 传统图像放大算法&#xff08;Lantent、Lanczos、Nearest&#xff09;AI图像放大算法&#xff08;4x-UltraSharp、BSRGAN、ESRGAN等&#xff09; 传统图像放大算法是基于插值算法&#xff0c;计算出图像放大后新位置的像素…

用友政务财务系统 FileDownload 任意文件读取漏洞复现

0x01 产品简介 用友政务财务系统具有多项核心功能,旨在满足各类组织的财务管理需求。首先,它提供了财务核算功能,能够全面管理企业的总账、固定资产、现金、应付应收等模块,实时掌握企业的财务状况,并通过科目管理、凭证处理、报表分析等功能为决策提供有力支持。 0x02 …

Elcomsoft iOS Forensics Toolkit: iPhone/iPad/iPod 设备取证工具包

天津鸿萌科贸发展有限公司是 ElcomSoft 系列取证软件的授权代理商。 Elcomsoft iOS Forensics Toolkit 软件工具包适用于取证工作&#xff0c;对 iPhone、iPad 和 iPod Touch 设备执行完整文件系统和逻辑数据采集。对设备文件系统制作镜像&#xff0c;提取设备机密&#xff08…

13.电子产品拆解分析-插排带3USB

13.电子产品拆解分析-插排带3USB 一、功能介绍二、电路分析以及器件作用1、三个插座之间通过电线连接,总开关控制火线2、通过FSD3773低待机功耗原边反馈AC/DC驱动芯片控制5V的输出一、功能介绍 ①一键控制总电源开关;②带三路USB输出;③最大支持2500W输出,10A输出电流;④8…

STM32自己从零开始实操01:原理图

在听完老师关于 STM32 物联网项目的所有硬件课程之后&#xff0c;就是感觉自己云里雾里&#xff0c;明明课程都认真听完了&#xff0c;笔记也认真记录&#xff0c;但是就是感觉学到的知识还不是自己。 遂决定站在老师的肩膀上自己开始设计项目&#xff0c;将知识变成自己的&am…

【机器学习】集成学习:强化机器学习模型与创新能的利器

集成学习&#xff1a;强化机器学习模型预测性能的利器 一、集成学习的核心思想二、常用集成学习方法Bagging方法Boosting方法Stacking方法 三、集成学习代表模型与实现四、总结与展望 在大数据时代的浪潮下&#xff0c;机器学习模型的应用越来越广泛&#xff0c;而集成学习作为…

微信小程序简单实现购物车功能

微信小程序简单实现购物车结算和购物车列表展示功能 实现在微信小程序中对每一个购物车界面的商品订单&#xff0c;进行勾选结算和取消结算的功能&#xff0c;相关界面截图如下&#xff1a; 具体实现示例代码为&#xff1a; 1、js代码&#xff1a; Page({/*** 页面的初始数…

Windows下搭建Flutter开发环境

IDE:VS code Flutter官网:Flutter: 为所有屏幕创造精彩 - Flutter 中文开发者网站 - Flutter 下载&安装 下载Flutter SDK,如图,建议自行下载安装: SDK还是挺大的,近1G,使用迅雷下载会快不少。 下载完成,解压缩到指定目录即可! 设置Local SDK,按下面步骤操作即…

【Docker】Docker 实践(一):在 Docker 中部署第一个应用

Docker 实践&#xff08;一&#xff09;&#xff1a;在 Docker 中部署第一个应用 1.使用 YUM 方式安装 Docker2.验证 Docker 环境3.在 Docker 中部署第一个应用3.1 小插曲&#xff1a;docker pull 报 missing signature key 错误3.2 重新安装 Nginx 1.使用 YUM 方式安装 Docker…

git 基础知识(全能版)

文章目录 一 、git 有三个分区二、git 基本操作1、克隆—git clone2、拉取—git fetch / git pull3、查看—git status / git diff3.1 多人开发代码暂存技巧 本地代码4、提交—git add / git commit / git push5、日志—git log / git reflog6、删除—git rm ‘name’7、撤销恢…

云计算和边缘计算究竟有什么不同

在数据时代&#xff0c;无论是人的活动还是机器的运作都会产生各种各样海量的数据。在对数据梳理和筛选过程中&#xff0c;计算机的运算处理必不可少。为了减少本地计算机算力成本等限制&#xff0c;越来越多的企业选择了云计算和边缘计算。今天&#xff0c;德迅云安全就带您来…

中电金信:深度解析|数字化营销运营体系搭建

如何更好更快地梳理好体系搭建思路&#xff0c;稳步实现落地&#xff1f;下文将为大家明确搭建的推进步骤、执行要点&#xff0c;帮助商业银行理顺数字化营销运营体系的“点”“线”“面”~ 与所有转型的曲折、阵痛等特征一样&#xff0c;商业银行构建数字化营销运营体系过程中…

线性代数 --- 计算斐波那契数列第n项的快速算法(矩阵的n次幂)

计算斐波那契数列第n项的快速算法(矩阵的n次幂) The n-th term of Fibonacci Numbers&#xff1a; 斐波那契数列的是一个古老而又经典的数学数列&#xff0c;距今已经有800多年了。关于斐波那契数列的计算方法不难&#xff0c;只是当我们希望快速求出其数列中的第100&#xff0…

普通屏幕已过时?裸眼3D屏幕显示效果更胜一筹!

随着多媒体技术的迅猛进步&#xff0c;我们日常生活中的内容展现方式&#xff0c;已经经历了前所未有的变革。在这其中&#xff0c;裸眼3D屏幕的应用&#xff0c;无疑是最为引人注目的亮点&#xff0c;它相较于传统屏幕&#xff0c;在显示效果上展现出了鲜明的优势&#xff0c;…

走向大规模应用之前,DePIN 如何突破技术、数据与市场之网

近期&#xff0c;随着分布式物理基础设施网络&#xff08;DePIN&#xff09;的快速演变&#xff0c;一个旨在利用区块链技术彻底改造传统基础设施模型的新兴生态系统正在逐渐浮现。2024 年 4 月&#xff0c;以 peaq 为代表的 DePIN 项目成功筹集了 1500 万美元用于生态系统的扩…

Unreal Engine动态添加Button实例

在控件蓝图中添加容器&#xff0c;注意命名不要有中文 C代码中找到容器实例 1 2 3 4 5 6 7 8 UVerticalBox* verticalBox Cast<UVerticalBox>(CurrentWidget->GetWidgetFromName(TEXT("VerticalBox_0"))); if (verticalBox ! nullptr) { UScrollBox* …

【声网】实现web端与uniapp微信小程序端音视频互动

实现web端与uniapp微信小程序端音视频互动 利用声网实现音视频互动 开通声网服务 注册声网账号 进入Console 成功登录控制台后&#xff0c;按照以下步骤创建一个声网项目&#xff1a; 展开控制台左上角下拉框&#xff0c;点击创建项目按钮。 在弹出的对话框内&#xff0c;依…