2.3 高斯變量

高斯分布
N(x|\mu,\sigma^2) = \frac{1}{(2\pi\sigma^2)^\frac{1}{2}} \exp\{-\frac{1}{2\sigma^2}(x-\mu)^2\}
多元高斯分布(D維)
N(x|\mu,\Sigma) = \frac{1}{(2\pi)^\frac{D}{2}}\frac{1}{|\Sigma|^\frac{1}{2}}\exp\{-\frac{1}{2}(x-\mu)^T\Sigma^{-1}(x-\mu)\}\\\Sigma是D*D的協(xié)方差矩陣,\mu是D維均值向量
本文旨在證明:\mu\Sigma為多元高斯分布的均值和方差


二次型\Delta^2
\Delta^2=(x-\mu)^T\Sigma^{-1}(x-\mu)\\\Delta叫做\mu和x之間的馬氏距離,\Sigma為單位矩陣時(shí)就變成歐式距離
矩陣\Sigma可以取對(duì)稱(chēng)矩陣,因?yàn)?strong>任何非對(duì)稱(chēng)項(xiàng)都會(huì)在指數(shù)中消失
考慮協(xié)方差特征向量方程\\ \Sigma\mu_i=\lambda\mu_i\\ 由于\Sigma為實(shí)對(duì)稱(chēng)矩陣,特征值也為實(shí)數(shù)猴誊,因此特征向量可以被選成單位正交\\ \mu_i^T\mu_j = I_{ij}\\協(xié)方差矩陣\Sigma可以展開(kāi)為\\\Sigma=\sum_{i=1}^D\lambda_i\mu_i\mu_i^T\\類(lèi)似的艰赞,\Sigma的逆矩陣\Sigma^{-1}可以寫(xiě)為\\\Sigma^{-1}=\sum_{i=1}^D\frac{1}{\lambda_i}\mu_i\mu_i^T
因此二次型\Delta^2可以寫(xiě)為
\Delta^2 = \sum^D_{i=1}\frac{y_i^2}{\lambda_i}\\其中y_i = \mu_i^T(x-\mu)
關(guān)于y
定義y=(y_1,...y_D)^T\\有\(zhòng)\y=U(x-\mu)\\可知道U時(shí)正交矩陣莱褒,滿足\\UU^T=U^TU=I
二次型\Delta^2y和原坐標(biāo)x的對(duì)應(yīng)關(guān)系


橢圓曲線表示二維空間
x=(x1,x2)
的高斯分布的常數(shù)概率密度的橢圓面浴栽,表示的概率密度為
e^{-\frac{1}{2}}
荒叼,值在
x=\mu
處計(jì)算。橢圓的軸由協(xié)方差矩陣的特征向量
\mu_i
定義典鸡,特征值(縮放因子)為
\lambda_i^\frac{1}{2}


x坐標(biāo)系到y坐標(biāo)系被廓,有jacobian矩陣J
J_{ij} = \frac{\delta x_i}{\delta y_i} = U_{ij}

U正交,因此
|J^2| = |U^T|^2 = |U^T||U| = |U||U^T| = |I| = 1\\ |J|=1

又行列式|\Sigma|可以寫(xiě)成特征值乘積|\Sigma|^{\frac{1}{2}} = \prod^D_{j=1}\lambda_j^{\frac{1}{2}}萝玷,故y坐標(biāo)系下嫁乘,高斯分布形式為
p(y) = p(x)|J|=\prod^D_{j=1}\frac{1}{(2\pi\lambda_j)^\frac{1}{2}}exp\{-\frac{y_j^2}{2\lambda_j}\}
以上公式是D個(gè)獨(dú)立一元高斯分布的乘積,特征向量定義了一個(gè)新的旋轉(zhuǎn)间护、平移的坐標(biāo)系亦渗,這個(gè)坐標(biāo)系下聯(lián)合概率分布可以分解成獨(dú)立分布的乘積

y坐標(biāo)系下的概率分布的積分為
\int p(y)d_y=\prod_{j=1}^D {\int}^{\infty}_{-\infty}\frac{1}{(2\pi\lambda_j)^\frac{1}{2}}\exp\{-\frac{y_j^2}{2\lambda_j}\}d{y_j}=1


以上是證明\mu\Sigma為高斯分布的均值和方差的前備條件,接下來(lái)進(jìn)行證明

  1. E[x]=\mu
    對(duì)連續(xù)概率密度函數(shù)求期望(積分)
    E[x] = \frac{1}{(2\pi)^\frac{D}{2}}\frac{1}{|\Sigma|^\frac{1}{2}}\int\exp\{-\frac{1}{2}(x-\mu)^T\Sigma^{-1}(x-\mu)\}xdx
    z=x-\mu
    E[x] = \frac{1}{(2\pi)^\frac{D}{2}}\frac{1}{|\Sigma|^\frac{1}{2}}\int\exp\{-\frac{1}{2}z^T\Sigma^{-1}z\}(z+\mu)dz
    由于積分區(qū)域是(-\infty,\infty),根據(jù)對(duì)稱(chēng)性可得(z+\mu)中的z項(xiàng)為零汁尺,因此:
    E[x]=\mu
    也就是證明了文章一開(kāi)始的D維均值向量就是多元高斯分布的均值(應(yīng)該是這樣??)

  2. var[x] = \Sigma
    求高斯分布的二階矩(PS:二階(非中心)矩是對(duì)變量的平方求期望法精,一階矩就是對(duì)變量求期望)
    一元變量下,二階矩由E[x^2]給出痴突;對(duì)于多元高斯分布搂蜓,有D^2個(gè)由E[x_ix_j]給出的二階矩,也就是矩陣E[xx^T]
    E[xx^T]=\frac{1}{(2\pi)^{\frac{D}{2}}}\frac{1}{|\Sigma|^{\frac{1}{2}}}\int\exp\{-\frac{1}{2}(x-\mu)^T\Sigma^{-1}(x-\mu)\}xx^Tdx
    z=x-\mu
    E[xx^T] = \frac{1}{(2\pi)^\frac{D}{2}}\frac{1}{|\Sigma|^\frac{1}{2}}\int\exp\{-\frac{1}{2}z^T\Sigma^{-1}z\}(z+\mu)(z+\mu)^Tdz\\=\frac{1}{(2\pi)^\frac{D}{2}}\frac{1}{|\Sigma|^\frac{1}{2}}\int\exp\{-\frac{1}{2}z^T\Sigma^{-1}z\}(zz^T+\mu\mu^T+\mu z^T+z\mu^T)dz
    由于對(duì)稱(chēng)性\mu z^Tz\mu^T項(xiàng)互相抵消辽装,\mu\mu^T為常數(shù)帮碰,因此我們先計(jì)算zz^T項(xiàng)
    這里的證明沒(méi)有看懂,最后可以寫(xiě)成\\ \frac{1}{(2\pi)^\frac{D}{2}}\frac{1}{|\Sigma|^\frac{1}{2}}\int\exp\{-\frac{1}{2}z^T\Sigma^{-1}z\} zz^Tdz\\ =\frac{1}{(2\pi)^\frac{D}{2}}\frac{1}{|\Sigma|^\frac{1}{2}}\sum^D_{i=1}\sum^D_{j=1}u_iu_j^T\int \exp\{-\sum^D_{k=1}\frac{y_k^2}{2\lambda_k}\}y_iy_jdy\\=\sum^D_{i=1} u_i u_i^T \lambda_i=\Sigma\\ 因此 E[xx^T]=\mu\mu^T+\Sigma
    定義協(xié)方差var[x]
    var[x] = E[(x-E[x])(x-E[x])^T]
    由于高斯分布E[x]=\mu拾积,結(jié)合E[xx^T]殉挽,得到
    var[x]=\Sigma
    也就是文首的D*D的協(xié)方差矩陣

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市拓巧,隨后出現(xiàn)的幾起案子斯碌,更是在濱河造成了極大的恐慌,老刑警劉巖肛度,帶你破解...
    沈念sama閱讀 211,290評(píng)論 6 491
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件傻唾,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡承耿,警方通過(guò)查閱死者的電腦和手機(jī)冠骄,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 90,107評(píng)論 2 385
  • 文/潘曉璐 我一進(jìn)店門(mén),熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)加袋,“玉大人惊楼,你說(shuō)我怎么就攤上這事显设〈辔辏” “怎么了黎做?”我有些...
    開(kāi)封第一講書(shū)人閱讀 156,872評(píng)論 0 347
  • 文/不壞的土叔 我叫張陵蝙砌,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我跋理,道長(zhǎng)择克,這世上最難降的妖魔是什么? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 56,415評(píng)論 1 283
  • 正文 為了忘掉前任前普,我火速辦了婚禮肚邢,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘拭卿。我一直安慰自己骡湖,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 65,453評(píng)論 6 385
  • 文/花漫 我一把揭開(kāi)白布峻厚。 她就那樣靜靜地躺著响蕴,像睡著了一般。 火紅的嫁衣襯著肌膚如雪惠桃。 梳的紋絲不亂的頭發(fā)上浦夷,一...
    開(kāi)封第一講書(shū)人閱讀 49,784評(píng)論 1 290
  • 那天,我揣著相機(jī)與錄音辜王,去河邊找鬼劈狐。 笑死,一個(gè)胖子當(dāng)著我的面吹牛呐馆,可吹牛的內(nèi)容都是我干的肥缔。 我是一名探鬼主播,決...
    沈念sama閱讀 38,927評(píng)論 3 406
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼汹来,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼续膳!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起收班,我...
    開(kāi)封第一講書(shū)人閱讀 37,691評(píng)論 0 266
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤姑宽,失蹤者是張志新(化名)和其女友劉穎,沒(méi)想到半個(gè)月后闺阱,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 44,137評(píng)論 1 303
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡舵变,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 36,472評(píng)論 2 326
  • 正文 我和宋清朗相戀三年酣溃,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片纪隙。...
    茶點(diǎn)故事閱讀 38,622評(píng)論 1 340
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡赊豌,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出绵咱,到底是詐尸還是另有隱情碘饼,我是刑警寧澤,帶...
    沈念sama閱讀 34,289評(píng)論 4 329
  • 正文 年R本政府宣布,位于F島的核電站艾恼,受9級(jí)特大地震影響住涉,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜钠绍,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,887評(píng)論 3 312
  • 文/蒙蒙 一舆声、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧柳爽,春花似錦媳握、人聲如沸。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,741評(píng)論 0 21
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至赵誓,卻和暖如春打毛,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背架曹。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,977評(píng)論 1 265
  • 我被黑心中介騙來(lái)泰國(guó)打工隘冲, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人绑雄。 一個(gè)月前我還...
    沈念sama閱讀 46,316評(píng)論 2 360
  • 正文 我出身青樓展辞,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親万牺。 傳聞我的和親對(duì)象是個(gè)殘疾皇子罗珍,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 43,490評(píng)論 2 348

推薦閱讀更多精彩內(nèi)容

  • 二元變量 伯努利分布 似然函數(shù)為 對(duì)數(shù)似然函數(shù)為 如果我們令關(guān)于u的導(dǎo)數(shù)等于零,我們就得到了最大似然的估計(jì)值 現(xiàn)在...
    初七123閱讀 1,548評(píng)論 0 4
  • 考試科目:高等數(shù)學(xué)脚粟、線性代數(shù)覆旱、概率論與數(shù)理統(tǒng)計(jì) 考試形式和試卷結(jié)構(gòu) 一、試卷滿分及考試時(shí)間 試卷滿分為150分核无,考...
    Saudade_lh閱讀 1,074評(píng)論 0 0
  • 2017年考研數(shù)學(xué)一大綱原文 考試科目:高等數(shù)學(xué)扣唱、線性代數(shù)、概率論與數(shù)理統(tǒng)計(jì) 考試形式和試卷結(jié)構(gòu) 一团南、試卷滿分及考...
    SheBang_閱讀 620評(píng)論 0 7
  • 考試形式和試卷結(jié)構(gòu)一噪沙、試卷滿分及考試時(shí)間 試卷滿分為150分,考試時(shí)間為180分鐘 二吐根、答題方式 答題方式為閉卷正歼、...
    幻無(wú)名閱讀 749評(píng)論 0 3
  • (1)孟凡剛 90天踐行總結(jié) 孟凡剛_65bc 作者 2017.07.17 11:46 打開(kāi)App 沒(méi)有記錄就沒(méi)有...
    珠寶女神肖肖閱讀 123評(píng)論 0 0