[線代] 高維向量夾角

看板Math作者 (真。盤古大洋)時間6年前 (2019/11/11 12:35), 6年前編輯推噓0(0022)
留言22則, 3人參與, 6年前最新討論串1/1
數學版的大家好,有一題線性代數想請教 考慮k 個N維隨機的向量 u1, u2, ... uk 滿足兩兩正交,且長度為 1 今有一個 N*N 矩陣 A (可能給定也可能隨機) 則問: (1) 若Au1 和 Au2 的夾角為 θ 2,則E[cos^2(θ 2)]=? 註:也就是 (u1’A’Au2)^2/(u1’A’Au1)(u2’A’Au2) 的期望值 (2) 推廣:若 Auk 與 span{Au1, ..., Au(k-1)} 之間的夾角為 θ k,則E[cos^2(θ k)]= 註:設 U=[Au1 ... Au(k-1)] 定義 P=U(U’U)^(-1)U’ 則所求就是 (uk’A’PAuk)^2/(uk’A’Auk)^2 的期望值 我想知道的是,(1) 和 (2) 的結果和 A 的什麼性質有關,是 singular values 嗎,或是 determinant 呢,還是其他的? 什麼關鍵字比較好找到類似的題目呢? 目前知道的兩個較 trivial 的情況: (i) A 是正交矩陣 那兩個結果都會是 0 (ii) A 的各 entries是 zero-mean i.i.d. Gaussian distributed r.v. 那(2)的答案似乎是 (k-1)/N 特別讓我在意的是這兩種 A: (iii) A 是由兩個矩陣相加組成,一個是 identity matrix I, 另一個則是各 entries是 zero-mean i.i.d. Gaussian distributed r.v. 的矩陣(給定 va riance 為 σ^2) (iv) A 是任意一個 Toeplitz 矩陣或circulant matrix 謝謝大家 ----- Sent from JPTT on my iPhone -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 111.71.127.185 (臺灣) ※ 文章網址: https://www.ptt.cc/bbs/Math/M.1573446929.A.600.html ※ 編輯: Panthalassa (111.71.127.185 臺灣), 11/11/2019 12:43:37 ※ 編輯: Panthalassa (111.71.127.185 臺灣), 11/11/2019 12:45:39

11/11 17:08, 6年前 , 1F
有沒有Ker(A) = {0}?
11/11 17:08, 1F

11/11 17:09, 6年前 , 2F
沒有的話期望值定義可能有零在分母
11/11 17:09, 2F

11/11 17:13, 6年前 , 3F
其實你的問題裡面只有 (A^T A)
11/11 17:13, 3F

11/11 17:14, 6年前 , 4F
diagonalizable by orthonormal transformation
11/11 17:14, 4F

11/11 17:15, 6年前 , 5F
eigenvalue 是 A 的 singular value 平方
11/11 17:15, 5F

11/11 17:16, 6年前 , 6F
反正先diagonalize,一切都應該好做很多
11/11 17:16, 6F

11/11 23:06, 6年前 , 7F
謝謝提醒 :) 不過有一部分是 A^T PA 比較難做對角化
11/11 23:06, 7F

11/11 23:06, 6年前 , 8F
P 我只知道 eigenvalues 為 0 or 1
11/11 23:06, 8F

11/11 23:06, 6年前 , 9F
這方面有沒有什麼領域在探討呢..隨機矩陣?
11/11 23:06, 9F

11/12 02:11, 6年前 , 10F
哪來的P?你下一步要推廣到中間夾projection嗎?
11/12 02:11, 10F

11/12 02:16, 6年前 , 11F
喔你的(2)
11/12 02:16, 11F

11/12 02:22, 6年前 , 12F
我現在頭腦不是很清楚,但是我覺得你做完SVD以後,
11/12 02:22, 12F

11/12 02:22, 6年前 , 13F
那個P在中間應該也是對角的
11/12 02:22, 13F

11/12 02:23, 6年前 , 14F
所以整個 A^t P A 也可以弄成對角
11/12 02:23, 14F

11/12 02:25, 6年前 , 15F
啊一想就知道果然錯了 XD 睡覺去zzz...
11/12 02:25, 15F

11/12 10:06, 6年前 , 16F
沒關係 謝謝XD
11/12 10:06, 16F

11/12 13:19, 6年前 , 17F
P 可以寫成 \sum (A uk)(uk' A')
11/12 13:19, 17F

11/12 13:20, 6年前 , 18F
把這個形式插進你的(2),你還是可以先對角化A'A再說
11/12 13:20, 18F

11/12 13:21, 6年前 , 19F
喔我忘了normalize,你看著辦吧
11/12 13:21, 19F

11/12 19:22, 6年前 , 20F
標題看成高雄...
11/12 19:22, 20F

11/13 11:39, 6年前 , 21F
謝謝 w大~
11/13 11:39, 21F

11/13 11:40, 6年前 , 22F
a大 真的有點像哈哈
11/13 11:40, 22F
文章代碼(AID): #1ToEKHO0 (Math)