Normal-Distribution

x’Ax 和 x’AAx 在單位球面上的預期比率?

  • August 10, 2021

認為 A 是對稱正定矩陣。以下量的第一時刻是否有一個很好的表達方式?

xTAxxTA2x

在哪裡 x 分佈為 Normal(0,In) . 這是在球體表面上計算的兩個二次型的比率。

什麼時候 A 有特徵值 1,12 ,這個期望等於 43 ,如下圖(筆記本

在此處輸入圖像描述

編輯 8 月 18 日 為了擴展超平面的答案,我們可以採取 A 在不失一般性的情況下成為對角線並將解寫為

E[xTAxxTA2x]=a,z

在哪裡 ai=1Aii zi=EyN(0,A)y2i|y|2

對於對稱矩陣 A,B , 數量 RA(x)=xTAxxTxRA,B(x)=xTAxxTBx 被稱為 (廣義)瑞利商。此處已提出有關此問題的問題:瑞利商的分佈或此處瑞利商的預期值

公認的答案是指Mathai 和 Provost於 1992 年出版的《隨機變量中的二次形式》一書。

在第 144 頁,我們參考了 Gurland 於 1956 年發表的論文正態分佈隨機變量中的二次形式,其中討論了廣義瑞利商的分佈和期望。除其他外,作者表明:

E[xTAxxTBx]=n1j=0k=0(1)j+12j+2vj+k+1cjgkB(j+k+1,3n2j1)

這裡, B(x,y)Beta 函數v , cjgk 是與特徵值/特徵多項式相關的係數 AB .

還有許多其他參考文獻給出了不同的系列/積分展開/表示的時刻 RA,B(x)

這些都沒有表明存在一個一般簡單的“封閉形式” E[RA,B(x)] .


給定特徵值分解,我們在任何情況下都可以做一些簡化步驟 B=UTΛU

ExN(0,𝕀)[xTAxxTBx]=EyN(0,𝕀)[yTUTAUyyTΛy]=EzN(0,Λ)[zTΛ1/2UTAUΛ1/2zzTz]

C=Λ1/2UTAUΛ1/2 並使用線性我們有:

EzN(0,Λ)[zTCzzTz]=EzN(0,Λ)[C,;zzTzTz]=C,;EzN(0,Λ)[zzTzTz]

在這裡,兩者 EzN(0,Λ)[zzT]=ΛEzN(0,Λ)[zTz]=tr(Λ) 是微不足道的,但是數值模擬表明 EzN(0,Λ)[zzTzTz] 是一個對角矩陣,其對角線具有一些非平凡的非線性關係 wrt Λ .

引用自:https://stats.stackexchange.com/questions/539720