卡方分布

本文深入探讨了卡方分布,包括中心卡方分布的概率密度函数及其与标准正态分布的关系,伽玛函数的概念,以及当k=2时转化为指数分布的情况。此外,还介绍了非中心卡方分布,解释了非中心参数λ的意义,并提供了随机变量的概率密度函数表达式。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

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)z2k1e2z0,,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α1exdx,如果我们令 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α1et2dt进一步,我们有
Γ ( 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为卡方分布的概率密度函数曲线。
在这里插入图片描述

图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)=21e2z0,,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=1kXi2满足非中心卡方分布,它有两个参数,一个是自由度 k k k,另外一个参数 λ \lambda λ与随机变量 X i X_i Xi的均值有关,即
λ = ∑ i = 1 k μ i 2 . \lambda=\sum_{i=1}^{k}\mu_i^2. λ=i=1kμ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=022k+2iΓ(2k+2i)i!e2λz(2λ)iz2k+2i1,其均值和方差分别为
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λ).

评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值