4.3.1 连续型随机变量

正态(高斯)分布

正态/高斯分布式最常用的分布之一。如果一个随机变量的概率密度函数为:
f x ( x ) = 1 2 π σ 2 exp ⁡ { − ( x − μ ) 2 σ 2 } , f_x(x) = \frac{1}{\sqrt{2 \pi \sigma^2}} \exp \{ \frac{-(x-\mu)}{2\sigma^2} \}, fx(x)=2πσ2 1exp{2σ2(xμ)} 这是关于参数 μ \mu μ(均值)为钟型的曲线。

图形特征

  • x = μ x=\mu x=μ处,取得最大值
  • 中间呈现凸形,两旁凹陷

它的分布函数由下面的变上限积分给出:
F x ( x ) = P ( X ≤ x ) = ∫ − ∞ x 1 2 π σ 2 exp ⁡ { − ( t − μ ) 2 2 σ 2 } d t ≡ G ( x − μ σ ) F_x(x) = P(X\leq x)=\int_{-\infty}^{x} \frac{1}{\sqrt{2 \pi \sigma^2}} \exp \{ - \frac{(t-\mu)^2}{2\sigma^2} \} \text{d} t \equiv G(\frac{x-\mu}{\sigma}) Fx(x)=P(Xx)=x2πσ2 1exp{2σ2(tμ)2}dtG(σxμ) 其中, G ( x ) = ∫ − ∞ x 1 2 π exp ⁡ { − y 2 2 } d y G(x) = \int_{-\infty}^{x} \frac{1}{2 \pi} \exp \{ -\frac{y^2}{2} \} \text{d} y G(x)=x2π1exp{2y2}dy 是最常使用的表达式。

因为 f x ( x ) f_x(x) fx(x)仅依赖两个参数 μ \mu μ σ 2 \sigma^2 σ2,常使用 x ∼ N ( μ , σ 2 ) \bm{x} \sim N(\mu,\sigma^2) xN(μ,σ2)来表达高斯概率密度函数,其中 ∼ \sim stands for “distributed as”。

标准正态分布的分布函数 Φ \bf{\Phi} Φ为:
Φ = P ( X ≤ x ) = 1 2 π ∫ − ∞ x exp ⁡ { − t 2 2 } d t \mathbf{\Phi} = P(X \leq x) = \frac{1}{\sqrt{2 \pi}} \int_{-\infty}^x \exp \{ -\frac{t^2}{2} \} \text{d} t Φ=P(Xx)=2π 1xexp{2t2}dt X ∼ N ( μ , σ 2 ) X \sim \mathcal{N}(\mu,\sigma^2) XN(μ,σ2),则 P ( x 1 < X < x 2 ) = Φ ( x 2 − μ σ ) − Φ ( x 1 − μ σ ) P(x_1<X<x_2) = \bf{\Phi}(\frac{x_2 - \mu}{\sigma}) - \bf{\Phi}(\frac{x_1 - \mu}{\sigma}) P(x1<X<x2)=Φ(σx2μ)Φ(σx1μ)

注意:常数 1 2 π σ 2 \frac{1}{\sqrt{2 \pi \sigma^2}} 2πσ2 1是一个归一化常数,它保证 f x ( x ) f_x(x) fx(x) ( − ∞ , ∞ ) (-\infty,\infty) (,)的积分等于1,这是因为 f ( x ) = ∫ − ∞ ∞ exp ⁡ { − x 2 2 σ 2 } = 2 π σ 2 f(x)=\int_{-\infty}^{\infty} \exp \{ -\frac{x^2}{2 \sigma^2} \} = \sqrt{2 \pi \sigma^2} f(x)=exp{2σ2x2}=2πσ2

特别地, x ∼ C N ( 0 , 1 ) \bm{x} \sim \mathcal{CN}(0,1) xCN(0,1)常被称为标准复高斯随机变量

性质

Independent Gaussian ± \pm ± Gaussian
  • 如果 X ∼ ( μ , σ 2 ) X \sim (\mu, \sigma^2) X(μ,σ2) a a a b b b是实数,则 a X ± b ∼ N ( a μ ± b , ( a σ ) 2 ) aX \pm b \sim N(a\mu \pm b,(a\sigma)^2) aX±bN(aμ±b,()2)
  • 如果 X ∼ ( μ X , σ X 2 ) X \sim (\mu_X, \sigma_X^2) X(μX,σX2) Y ∼ ( μ Y , σ Y 2 ) Y \sim (\mu_Y, \sigma_Y^2) Y(μY,σY2)是统计独立(independent statistical)的正态随机变量 ,则它们的和/差也满足高斯分布 U = X + Y ∼ ( μ X + μ Y , σ X 2 + σ Y 2 ) U=X+Y \sim (\mu_X+\mu_Y, \sigma_X^2+\sigma_Y^2) U=X+Y(μX+μY,σX2+σY2) V = X − Y ∼ ( μ X − μ Y , σ X 2 + σ Y 2 ) V=X-Y \sim (\mu_X-\mu_Y, \sigma_X^2+\sigma_Y^2) V=XY(μXμY,σX2+σY2) 其中U和V相互独立。
Gaussian × / ÷ \times/ \div ×/÷ Gaussian
  • 如果 X ∼ ( 0 , σ X 2 ) X \sim (0, \sigma_X^2) X(0,σX2) Y ∼ ( 0 , σ Y 2 ) Y \sim (0, \sigma_Y^2) Y(0,σY2)是独立(Independent)的正态随机变量, 则积 X Y XY XY服从概率密度函数为 p p p的分布
    p ( z ) = 1 π σ X σ Y I 0 ( ∣ z ∣ σ X σ Y ) . p(z)=\frac{1}{\pi \sigma_{X} \sigma_{Y}} I_{0}\left(\frac{|z|}{\sigma_{X} \sigma_{Y}}\right). p(z)=πσXσY1I0(σXσYz). 其中, I 0 ( x ) I_0(x) I0(x)是零阶修正贝塞尔函数(Zeroth-order modified Bessel function): I n ( α ) = 1 π ∫ 0 π e α cos ⁡ x   d x I_{n}(\alpha)=\frac{1}{\pi} \int_{0}^{\pi}e^{\alpha \cos x} \mathrm{~d} x In(α)=π10πeαcosx dx
    注意: n n n阶修正贝塞尔函数为: I n ( α ) = 1 π ∫ 0 π cos ⁡ ( n x ) e α cos ⁡ x   d x I_{n}(\alpha)=\frac{1}{\pi} \int_{0}^{\pi} \cos (n x) e^{\alpha \cos x} \mathrm{~d} x In(α)=π10πcos(nx)eαcosx dx
Z = X 1 2 + X 2 2 + . . . + X n 2 Z = X_1^2 + X_2^2 +...+ X_n^2 Z=X12+X22+...+Xn2
  • 如果 X 1 , X 2 , . . . , X N X_1, X_2,...,X_N X1,X2,...,XN满足独立标准正态随机变量,即 X 1 , X 2 , . . . , X N ∼ C N ( 0 , 1 ) X_1, X_2,...,X_N \sim \mathcal{CN}(0,1) X1,X2,...,XNCN(0,1),则 X 1 2 + X 2 2 + . . . + X N 2 X_1^2 + X_2^2 + ... +X_N^2 X12+X22+...+XN2满足自由度为 n n n卡方分布

例如:

Complex Gaussian distriubution Z = X + i Y ∼ C N ( 0 , 1 ) Z = X + iY \sim \mathcal{CN}(0,1) Z=X+iYCN(0,1), X ∼ N ( 0 , 1 / 2 ) X \sim N(0,1/2) XN(0,1/2) and Y ∼ N ( 0 , 1 / 2 ) Y \sim N(0,1/2) YN(0,1/2).

∣ Z ∣ 2 = X 2 + Y 2 |Z|^2 = X^2 + Y^2 Z2=X2+Y2 where ∣ Z ∣ |Z| Z follows the Reyleigh distribution.

2 ∣ Z ∣ 2 = 2 X 2 + 2 Y 2 ∼ χ 2 2 2|Z|^2 = 2X^2 + 2Y^2 \sim \chi _{2}^{2} 2∣Z2=2X2+2Y2χ22,where 2 ∣ Z ∣ 2 2|Z|^2 2∣Z2 follows Chi-squared distribition, 2 X 2 ∼ N ( 0 , 1 / 2 ) 2X^2 \sim N(0,1/2) 2X2N(0,1/2), and 2 Y 2 ∼ N ( 0 , 1 / 2 ) 2Y^2\sim N(0,1/2) 2Y2N(0,1/2).

复正态(高斯)分布

若复高斯分布 Z = X + i Y Z=X+iY Z=X+iY,且满足 X ∼ ( μ x , σ x 2 ) X \sim (\mu_x,\sigma^2_x) X(μx,σx2), Y ∼ ( μ y , σ y 2 ) Y \sim (\mu_y,\sigma^2_y) Y(μy,σy2), μ = μ x = μ y \mu=\mu_x=\mu_y μ=μx=μy, σ = σ x 2 = σ y 2 \sigma = \sigma^2_x = \sigma^2_y σ=σx2=σy2,则 μ z = μ x + i μ y , σ z 2 = 2 σ x 2 = 2 σ y 2 \textcolor{red}{\mu_z=\mu_x + i \mu_y,\sigma^2_z =2 \sigma^2_x = 2\sigma^2_y} μz=μx+iμy,σz2=2σx2=2σy2.
The probability density function of the 2-D random variable ( x , y ) (x,y) (x,y) is
p x y = 1 2 π σ 2 exp ⁡ { − ( x − μ z ) 2 + ( y − μ y ) 2 2 σ 2 } . p_{x y}=\frac{1}{2 \pi \sigma^{2}} \exp \left\{ -\frac{\left(x-\mu_{z}\right)^{2}+\left(y-\mu_{y}\right)^{2}}{2 \sigma^{2}} \right\}. pxy=2πσ21exp{2σ2(xμz)2+(yμy)2}. Noting that Z = X + i Y Z=X+iY Z=X+iY, the probability density function can be represented in terms of z z z: p z = 1 2 π σ 2 exp ⁡ { − ( z − μ z ) 2 2 σ 2 } . p_{z}=\frac{1}{2 \pi \textcolor{red}{\sigma^{2}}} \exp \left\{-\frac{\left(z-\textcolor{red}{\mu_{z}}\right)^{2}}{\textcolor{red}{2 \sigma^{2}}} \right\}. pz=2πσ21exp{2σ2(zμz)2}. p z = 1 π σ z 2 exp ⁡ { − ( z − μ z ) 2 σ z 2 } . p_{z}=\frac{1}{ \pi \textcolor{red}{\sigma_{z}^{2}}} \exp \left\{-\frac{\left(z-\textcolor{red}{\mu_{z}}\right)^{2}}{\textcolor{red}{ \sigma_{z}^{2}}} \right\}. pz=πσz21exp{σz2(zμz)2}.
注意:复高斯随机变量的密度函数,分母无根号

零均值循环对称负高斯随机变量
μ = μ x = μ y = 0 \mu=\mu_x=\mu_y=0 μ=μx=μy=0 Z Z Z称为零均值循环对称复高斯随机变量(Zero Mean Circle Symmetric Complex Gaussian, ZMCSCG), σ 2 \sigma^2 σ2称为每个实数维度上的方差。
X ∼ C N ( 0 , 1 / 2 ) ,   Y ∼ C N ( 1 / 2 ) X \sim \mathcal{CN}(0,1/2), ~Y \sim \mathcal{CN}(1/2) XCN(0,1/2), YCN(1/2),则 Z ∼ C N ( 0 , 1 ) . Z\sim \mathcal{CN}(0,1). ZCN(0,1).

与正态分布相关的函数

1. Q函数

Q函数又称标准正态分布的右尾函数
Q ( x ) = P ( X > x ) = ∫ x ∞ 1 2 π exp ⁡ { − t 2 2 } d t = 1 − Φ ( x ) . Q(x) =P(X>x) =\int_{x}^{\infty} \frac{1}{\sqrt{2 \pi}} \exp \{ -\frac{t^2}{2} \} \text{d} t =1-\bf{\Phi}(x). Q(x)=P(X>x)=x2π 1exp{2t2}dt=1Φ(x).

  • Q ( x ) Q(x) Q(x)与标准正态分布函数: Q ( x ) = 1 − Φ ( x ) Q(x) = 1-\mathbf{\Phi}(x) Q(x)=1Φ(x)

2. 误差函数(Error Function)

正态函数的积分是一个非基本函数(即不是初等函数), 称为误差函数(高斯误差函数,Gauss error function)。
e r f ( x ) = 2 π ∫ 0 x exp ⁡ { − t 2 } d t erf(x) = \frac{2}{\sqrt{\pi}} \int_{0}^{x} \exp \{-t^2 \} \text{d} t erf(x)=π 20xexp{t2}dt

Q ( x ) Q(x) Q(x) e r f erf erf函数:

  • e r f ( x ) = 1 − 2 Q ( 2 x ) {erf}(x)=1-2 Q(\sqrt{2} x) erf(x)=12Q(2 x)
3. 互补误差函数(Complementary Error Function)

e r f c ( x ) = 2 π ∫ x ∞ exp ⁡ { − t 2 } d t = 1 − e r f ( x ) erfc(x) = \frac{2}{\sqrt{\pi}} \int_{x}^{\infty} \exp \{-t^2 \} \text{d} t = 1 - erf(x) erfc(x)=π 2xexp{t2}dt=1erf(x) e r f ( x ) + e r f c ( x ) = 1 erf(x) + erfc(x) = 1 erf(x)+erfc(x)=1

Q ( x ) Q(x) Q(x) e r f c erfc erfc函数:

  • Q ( x ) = 1 2 e r f c ( x / 2 ) Q(x) = \frac{1}{2} erfc(x/\sqrt{2}) Q(x)=21erfc(x/2 );
  • 2 Q ( 2 x ) = e r f c ( x ) 2Q(\sqrt{2}x) = erfc(x) 2Q(2 x)=erfc(x);

瑞丽(Rayleigh)分布

如果随机变量 x \bm{x} x的概率密度函数是: f x ( x ) = { x σ 2 exp ⁡ { − x 2 2 σ 2 } , x ≥ 0 0 , x < 0 f_x(x) = \begin{cases} \frac{x}{\sigma^2} \exp\{ - \frac{x^2}{2\sigma^2} \},& x \geq 0 \\ \qquad 0, & x<0 \end{cases} fx(x)={σ2xexp{2σ2x2},0,x0x<0
注意:

  • 在通信系统中,随机接收信号的幅度通常用瑞利分布描述。
  • 服从复高斯分布的变量的模用瑞丽分布描述,例如, X ∼ C N ( 0 , σ 2 ) X\sim \mathcal{CN}(0,\sigma^2) XCN(0,σ2), ∣ X ∣ |X| X服从瑞利分布,该分布的概率密度函数 f x ( x ) = x σ 2 exp ⁡ { − x 2 2 σ 2 } , x ≥ 0 f_x(x) = \frac{x}{\sigma^2} \exp\{ - \frac{x^2}{2\sigma^2} \}, x \geq 0 fx(x)=σ2xexp{2σ2x2},x0, 累积分布函数 F x ( x ) = 1 − exp ⁡ ( − x 2 2 σ 2 ) . F_x(x)=1-\exp\left( -\frac{x^2}{2\sigma^2}\right). Fx(x)=1exp(2σ2x2).
  • 如果 A A A服从瑞丽分布,则有 E ( A ) = π 2 σ \mathbb{E} (A) = \sqrt{\frac{\pi}{2}} \sigma E(A)=2π σ D ( A ) = 4 − π 2 σ \mathbb{D} (A) = \frac{4-\pi}{2}\sigma D(A)=24πσ
  • 信道 h = X + ȷ Y ∼ C N ( 0 , σ 2 ) h=X+ \jmath Y \sim \mathcal{CN}(0,\sigma^2) h=X+YCN(0,σ2)服从复高斯分布,其中 X , Y ∼ N ( 0 , σ 2 2 ) X, Y \sim \mathcal{N}(0,\frac{\sigma^2}{2}) X,YN(0,2σ2). 则 ∣ h ∣ = X 2 + Y 2 |h|= \sqrt{X^2 + Y^2} h=X2+Y2 服从瑞丽分布,其均值为 E ( ∣ h ∣ ) = π 2 ∗ σ 2 2 = π 2 σ \mathbb{E}(|h|) = \sqrt{\frac{\pi}{2}} *\sqrt{ \frac{\sigma^2}{2} } = \frac{\pi}{2}\sigma E(h)=2π 2σ2 =2πσ.

Nakagami-m分布

通过引入参数 m m m,从瑞利分布可以导出一类更一般的分布——Nakagami- m m m分布,其概率密度函数为:
f x ( x ) = { 2 Γ ( m ) ( m Ω ) m x 2 m − 1 exp ⁡ { − m x 2 Ω } , x > 0   0 , x ≤ 0 f_x(x) = \begin{cases} \frac{2}{\Gamma(m)} \left( \frac{m}{\Omega} \right)^m x^{2m-1} \exp\{ - \frac{mx^2}{\Omega} \},& x > 0 \\ \qquad \qquad ~ 0, & x \leq 0 \end{cases} fx(x)={Γ(m)2(Ωm)mx2m1exp{Ωmx2}, 0,x>0x0
与瑞丽分布相比,Nakagami- m m m分布在模拟通信理论中衰落信道时具有更大的灵活性。在上式中, m = 1 m=1 m=1对应瑞丽分布,并且通过调整 m m m可控制密度函数的拖尾。如下图所示,当 m < 1 m<1 m<1时,密度函数的拖尾衰减比瑞利分布慢, m > 1 m>1 m>1,拖尾衰减比瑞利分布快。

均匀分布

若随机变量 x \bm{x} x的密度函数在区间 ( a , b ) (a,b) (a,b)是常数, − ∞ < a < b < + ∞ -\infty<a<b<+\infty <a<b<+,而其他地方为零,即 f x ( x ) = { 1 b − a , a ≤ x ≤ b    0 , o t h e r s f_x(x) = \begin{cases} \frac{1}{b-a} ,& a \leq x \leq b \\ ~~0, & others \end{cases} fx(x)={ba1,  0,axbothers 则称它服从均匀分布,纪作 x ∼ U ( a , b ) \bm{x}\sim U(a,b) xU(a,b)。密度函数如下图所所示,其分布函数为:
F x ( x ) = { 1 , x ≥ b x − a b − a , a ≤ x < b    0 , x < a . F_x(x) = \begin{cases} 1, &x \geq b \\ \frac{x-a}{b-a} ,& a \leq x < b \\ ~~0, & x<a\end{cases}. Fx(x)= 1,baxa,  0,xbax<bx<a.


卡方分布(Chi-Square Distribution)

  • 如果 X 1 , X 2 , . . . , X n X_1, X_2,...,X_n X1,X2,...,Xn满足独立标准正态随机变量, X 1 , X 2 , . . . , X N ∼ C N ( 0 , 1 ) X_1, X_2,...,X_N \sim \mathcal{CN}(0,1) X1,X2,...,XNCN(0,1),则 X 2 = X 1 2 + X 2 2 + . . . + X n 2 X^2=X_1^2 + X_2^2 + ... +X_n^2 X2=X12+X22+...+Xn2满足自由度为 n n n卡方分布, 即 X 2 ∼ χ 2 ( 2 ) X^2 \sim \chi^2(2) X2χ2(2), 其中自由度为独立变量的个数

  • 卡方分布的概率密度为
    f n ( x ) = { 1 2 Γ ( n / 2 ) ( x 2 ) n 2 − 1 e − x 2 , x > 0 0 , x ≤ 0 f_{n}(x)=\left\{\begin{array}{cc} \frac{1}{2 \Gamma(n / 2)}\left(\frac{x}{2}\right)^{\frac{n}{2}-1} e^{-\frac{x}{2}}, & x>0 \\ 0, & x \leq 0 \end{array}\right. fn(x)={(n/2)1(2x)2n1e2x,0,x>0x0 其中 Γ ( α ) = ∫ 0 + ∞ x α − 1 e − x d x \Gamma(\alpha) = \int_{0}^{+\infty} x^{\alpha-1} e^{-x} \text{d} x Γ(α)=0+xα1exdx.

  • E ( X 2 ) = n , D ( X 2 ) = 2 n E(X^2)=n, D(X^2)=2n E(X2)=n,D(X2)=2n

  • 卡方分布是由正态分布构造而成的一个新的分布,当自由度 n n n很大时,卡方分布近似为正态分布。

  • n n n为2时,分布的CDF为: F ( x ) = ∫ − ∞ x exp ⁡ ( − 1 2 t ) d t F(x) = \int_{-\infty}^{x} \exp (- \frac{1}{2}t) \text{d} t F(x)=xexp(21t)dt


Logo

旨在为数千万中国开发者提供一个无缝且高效的云端环境,以支持学习、使用和贡献开源项目。

更多推荐