
概率论-分布函数(高斯分布、复高斯分布、瑞丽分布、Nakagami-m分布、均匀分布、卡方分布)
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πσ21exp{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(X≤x)=∫−∞x2πσ21exp{−2σ2(t−μ)2}dt≡G(σ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) x∼N(μ,σ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(X≤x)=2π1∫−∞xexp{−2t2}dt 令
X
∼
N
(
μ
,
σ
2
)
X \sim \mathcal{N}(\mu,\sigma^2)
X∼N(μ,σ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πσ21是一个归一化常数,它保证 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) x∼CN(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±b∼N(aμ±b,(aσ)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=X−Y∼(μ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σY∣z∣). 其中, 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(α)=π1∫0π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(α)=π1∫0π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,...,XN∼CN(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+iY∼CN(0,1), X ∼ N ( 0 , 1 / 2 ) X \sim N(0,1/2) X∼N(0,1/2) and Y ∼ N ( 0 , 1 / 2 ) Y \sim N(0,1/2) Y∼N(0,1/2).
∣ Z ∣ 2 = X 2 + Y 2 |Z|^2 = X^2 + Y^2 ∣Z∣2=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∣Z∣2=2X2+2Y2∼χ22,where 2 ∣ Z ∣ 2 2|Z|^2 2∣Z∣2 follows Chi-squared distribition, 2 X 2 ∼ N ( 0 , 1 / 2 ) 2X^2 \sim N(0,1/2) 2X2∼N(0,1/2), and 2 Y 2 ∼ N ( 0 , 1 / 2 ) 2Y^2\sim N(0,1/2) 2Y2∼N(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)
X∼CN(0,1/2), Y∼CN(1/2),则
Z
∼
C
N
(
0
,
1
)
.
Z\sim \mathcal{CN}(0,1).
Z∼CN(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)=∫x∞2π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)=π2∫0xexp{−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)=1−2Q(2x)
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)=π2∫x∞exp{−t2}dt=1−erf(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(2x)=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,x≥0x<0
注意:
- 在通信系统中,随机接收信号的幅度通常用瑞利分布描述。
- 服从复高斯分布的变量的模用瑞丽分布描述,例如, X ∼ C N ( 0 , σ 2 ) X\sim \mathcal{CN}(0,\sigma^2) X∼CN(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},x≥0, 累积分布函数为 F x ( x ) = 1 − exp ( − x 2 2 σ 2 ) . F_x(x)=1-\exp\left( -\frac{x^2}{2\sigma^2}\right). Fx(x)=1−exp(−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+Y∼CN(0,σ2)服从复高斯分布,其中 X , Y ∼ N ( 0 , σ 2 2 ) X, Y \sim \mathcal{N}(0,\frac{\sigma^2}{2}) X,Y∼N(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)mx2m−1exp{−Ωmx2}, 0,x>0x≤0
与瑞丽分布相比,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)={b−a1, 0,a≤x≤bothers 则称它服从均匀分布,纪作
x
∼
U
(
a
,
b
)
\bm{x}\sim U(a,b)
x∼U(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,b−ax−a, 0,x≥ba≤x<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,...,XN∼CN(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)={2Γ(n/2)1(2x)2n−1e−2x,0,x>0x≤0 其中 Γ ( α ) = ∫ 0 + ∞ x α − 1 e − x d x \Gamma(\alpha) = \int_{0}^{+\infty} x^{\alpha-1} e^{-x} \text{d} x Γ(α)=∫0+∞xα−1e−xdx. -
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

更多推荐



所有评论(0)