• 正态分布密度函数的系数


    正态分布的密度函数,可以一般化地写为

    [f(x) = k expleft[-dfrac{1}{2}(x-b)' A (x-b) ight] ]

    事实上,如果某个多维随机变量的密度函数可以写成该形式,那么它就服从正态分布。其中(b)是均值,正定矩阵(A)是协方差矩阵的逆,它们共同决定的正态分布的形式。而另外一个字母(k),仅仅是归一化系数,它是满足整个密度函数的积分等于(1)的那个值。

    如果有人背过公式,会发现这个系数的形式比较复杂。本文具体来看看,它是怎么计算出来的。

    由于(A)是正定的,必有分解(A=CC')。先做个变换,令(x-b=(C')^{-1}y),那么

    [(x-b)' A (x-b) = y' C^{-1}A(C')^{-1}y = y'y ]

    同时,该变换的Jacobian matrix为(J = det [(C')^{-1}]=1/det(C))

    假设(x)(d)维,则(y)也是(d)维,将其各维写出,有(y = (y_1,cdots,y_d)')。接下来,对密度函数进行积分:

    [egin{aligned} &int_{-infty}^{infty} cdots int_{-infty}^{infty} f(x) dx_1 cdots dx_d\ =& int_{-infty}^{infty}cdotsint_{-infty}^{infty} kexp(-dfrac{1}{2}y'y) |J| dy_1 cdots dy_d\ =& dfrac{k}{|det(C)|} int_{-infty}^{infty}cdotsint_{-infty}^{infty} prod_{i=1}^{d}exp(-dfrac{1}{2} y_i^2) dy_1 cdots dy_d\ =& dfrac{k}{|det(C)|} prod_{i=1}^{d} int_{-infty}^{infty} exp(-dfrac{1}{2} y_i^2) dy_i \ =& dfrac{k}{|det(C)|} prod_{i=1}^{d} sqrt{2pi} \ =& dfrac{k}{|det(C)|} (2pi)^{d/2} \ =& k [det(A)]^{-1/2} (2pi)^{d/2} end{aligned} ]

    上述积分必定等于(1),因此,

    [k=(2pi)^{-d/2} [det(A)]^{1/2} ]

    同名公众号:分析101
  • 相关阅读:
    美国大学排名之本科中最用功的学校top15
    PhpStorm (强大的PHP开发环境)2017.3.2 附注册方法
    获取地址栏的URL: PHP JS
    怎么给php下拉框默认选中
    在JS中使用全局变量
    原生和jQuery的ajax用法
    XAMPP重要文件目录及配置
    select获取下拉框的值 下拉框默认选中
    h5 时间控件问题,怎么设置type =datetime-local 的值
    JS截取字符串常用方法详细整理
  • 原文地址:https://www.cnblogs.com/analysis101/p/14873278.html
Copyright © 2020-2023  润新知