每週問題 November 14, 2016

證明一個常見於多變量統計學的矩陣 I_m-A(A^{\ast}A)^{-1}A^\ast 是半正定。

Let m\ge n and let A be a complex m\times n matrix of rank n. Show that the Hermitian matrix B=I_m-A(A^{\ast}A)^{-1}A^\ast is positive semidefinite.

 
參考解答:

因為 \hbox{rank}(A^\ast A)=\hbox{rank}A=nn\times n 階 Gramian 矩陣 A^\ast A 是可逆的。直接計算可得 B=B^\ast=B^2B 是一個 Hermitian 冪等 (idempotent) 矩陣,也就是正交投影矩陣。因此,B 的特徵值為 01,即得證。事實上,冪等矩陣 B 滿足 \hbox{rank}B=\hbox{trace}B,而

\begin{aligned}  \hbox{trace}B&=\hbox{trace}(I_m-A(A^{\ast}A)^{-1}A^\ast)\\  &=\hbox{trace}{I_m}-\hbox{trace}(A(A^{\ast}A)^{-1}A^\ast)\\  &=\hbox{trace}I_m-\hbox{trace}((A^{\ast}A)^{-1}A^\ast A)\\  &=\hbox{trace}I_m-\hbox{trace}I_n\\  &=m-n,  \end{aligned}

故知 Bm-n 個特徵值 1n 個特徵值 0

廣告
本篇發表於 pow 內積空間, 每週問題 並標籤為 , , , 。將永久鏈結加入書籤。

發表迴響

在下方填入你的資料或按右方圖示以社群網站登入:

WordPress.com Logo

您的留言將使用 WordPress.com 帳號。 登出 / 變更 )

Twitter picture

您的留言將使用 Twitter 帳號。 登出 / 變更 )

Facebook照片

您的留言將使用 Facebook 帳號。 登出 / 變更 )

Google+ photo

您的留言將使用 Google+ 帳號。 登出 / 變更 )

連結到 %s