1、中心卡方分布
如果我们有
k
k
k个独立的标准正态分布随机变量
x
1
,
x
2
,
…
,
x
k
x_1,x_2,\ldots,x_k
x1,x2,…,xk,则其平方和
Z
=
∑
i
=
1
k
x
i
2
Z=\sum_{i=1}^kx_i^2
Z=∑i=1kxi2满足自由度为
k
k
k的卡方分布,概率密度函数为
(1)
p
Z
(
z
)
=
{
z
k
2
−
1
e
−
z
2
2
k
2
Γ
(
k
2
)
,
z
>
0
0
,
o
t
h
e
r
w
i
s
e
.
\tag{1} p_Z(z)=\left\{\begin{aligned} \frac{z^{\frac{k}{2}-1}e^{-\frac{z}{2}}}{2^{\frac{k}{2}}\Gamma(\frac{k}{2})}&,&z>0\\ 0&,&{\rm otherwise}. \end{aligned} \right.
pZ(z)=⎩⎪⎨⎪⎧22kΓ(2k)z2k−1e−2z0,,z>0otherwise.(1)
其均值为自由度
k
k
k,方差为
2
k
2k
2k。
下面我们首先来看看伽玛分布。
2、伽玛函数
下面我们来看看Gamma函数的定义。
Γ
(
α
)
=
∫
0
+
∞
x
α
−
1
e
−
x
d
x
,
\Gamma(\alpha)=\int_{0}^{+\infty}x^{\alpha-1}e^{-x}dx,
Γ(α)=∫0+∞xα−1e−xdx,如果我们令
x
=
t
2
x=t^2
x=t2,则有
Γ
(
α
)
=
∫
0
+
∞
2
t
2
α
−
1
e
−
t
2
d
t
。
\Gamma(\alpha)=\int_{0}^{+\infty}2t^{2\alpha-1}e^{-t^2}dt。
Γ(α)=∫0+∞2t2α−1e−t2dt。进一步,我们有
Γ
(
1
)
=
1
Γ
(
1
2
)
=
π
Γ
(
α
+
1
)
=
α
Γ
(
α
)
Γ
(
n
+
1
)
=
n
Γ
(
n
)
=
n
!
\begin{aligned} &\Gamma(1)=1\\ &\Gamma(\frac{1}{2})=\sqrt{\pi}\\ &\Gamma(\alpha+1)=\alpha\Gamma(\alpha)\\ &\Gamma(n+1)=n\Gamma(n)=n! \end{aligned}
Γ(1)=1Γ(21)=πΓ(α+1)=αΓ(α)Γ(n+1)=nΓ(n)=n!推荐一个关于Gamma函数的视频(https://www.youtube.com/watch?v=ixaz8-q90L8)。
3、中心卡方分布的概率密度函数曲线
//下面我们用MATLAB看下卡方分布的概率密度函数曲线。
clear;
for k=1:12
for z1=1:1400
z=z1/100;
pdf(k,z1)=power(z,k/2-1)*exp(-z/2)/power(2,k/2)/gamma(k/2);
end
end
x=0.01:0.01:14;
plot(x,pdf(1,:),'r')
hold on;
plot(x,pdf(2,:),'y')
plot(x,pdf(4,:),'g')
plot(x,pdf(6,:),'b')
plot(x,pdf(11,:),'k')
当然也可以调用MATLAB里面的函数chi2pdf(Z,k)。图1为卡方分布的概率密度函数曲线。
4、指数分布
若
k
=
2
k=2
k=2,则
Z
Z
Z满足指数分布,即
(2)
p
Z
(
z
)
=
{
1
2
e
−
z
2
,
z
>
0
0
,
o
t
h
e
r
w
i
s
e
.
\tag{2} p_Z(z)=\left\{\begin{aligned} \frac{1}{2}e^{-\frac{z}{2}}&,&z>0\\ 0&,&{\rm otherwise}. \end{aligned} \right.
pZ(z)=⎩⎨⎧21e−2z0,,z>0otherwise.(2)
5、非中心卡方分布
令
X
1
,
X
2
,
…
,
X
i
,
…
,
X
k
X_1,X_2,\ldots,X_i,\ldots,X_k
X1,X2,…,Xi,…,Xk为
k
k
k个独立且正态分布的随机变量,其中
X
i
X_i
Xi的均值为
μ
i
\mu_i
μi,方差为1,这里
i
=
1
,
2
,
…
,
k
i=1,2,\ldots,k
i=1,2,…,k,则随机变量
Z
=
∑
i
=
1
k
X
i
2
Z=\sum_{i=1}^{k}X_i^2
Z=i=1∑kXi2满足非中心卡方分布,它有两个参数,一个是自由度
k
k
k,另外一个参数
λ
\lambda
λ与随机变量
X
i
X_i
Xi的均值有关,即
λ
=
∑
i
=
1
k
μ
i
2
.
\lambda=\sum_{i=1}^{k}\mu_i^2.
λ=i=1∑kμi2.
λ
\lambda
λ有时也被称为非中心参数。非中心卡方分布随机变量
Z
Z
Z的概率密度函数为
p
Z
(
z
)
=
∑
i
=
0
∞
e
−
λ
−
z
2
(
λ
2
)
i
z
k
+
2
i
2
−
1
2
k
+
2
i
2
Γ
(
k
+
2
i
2
)
i
!
,
p_Z(z)=\sum_{i=0}^{\infty}\frac{e^{\frac{-\lambda-z}{2}}(\frac{\lambda}{2})^iz^{\frac{k+2i}{2}-1}}{2^{\frac{k+2i}{2}}\Gamma(\frac{k+2i}{2})i!},
pZ(z)=i=0∑∞22k+2iΓ(2k+2i)i!e2−λ−z(2λ)iz2k+2i−1,其均值和方差分别为
E
(
Z
)
=
λ
+
k
V
a
r
(
Z
)
=
2
(
k
+
2
λ
)
.
\begin{aligned} {\rm E}(Z)&=\lambda+k\\ {\rm Var}(Z)&=2(k+2\lambda). \end{aligned}
E(Z)Var(Z)=λ+k=2(k+2λ).