[理工] [線代] SVD

看板Grad-ProbAsk作者 (好崩潰)時間14年前 (2011/12/02 00:05), 編輯推噓7(7014)
留言21則, 7人參與, 最新討論串1/3 (看更多)
T 奇異值分解A=UΣV T 從AA 可得特徵根與特徵向量,其特徵向量做單範正交化 可得[u1,u2,...,uN]=U的行向量 T 從A A 可得特徵根與特徵向量,其特徵向量做單範正交化 可得[v1,v2,...,vN]=V的行向量 請問U跟V的行向量順序是唯一嗎?? 若不唯一,又要如從特徵向量來選取呢?? 也就是說萬一算出來的特徵根為重根的話 勢必span{}會有多個基底 那我該如何選取?? 有請高手們大心解答 鋼溫!! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 223.139.139.218

12/02 00:33, , 1F
奇異值擺的位置決定行向量位置
12/02 00:33, 1F

12/02 00:50, , 2F
其實多算幾題你就會發現,跟對角化很像
12/02 00:50, 2F

12/02 01:07, , 3F
通常都是先算出V 也就是(A^T)A 的eigenvector之後,再
12/02 01:07, 3F

12/02 01:08, , 4F
把V的行向量一一取出來, 假設 v0是V的其中一個行向量
12/02 01:08, 4F

12/02 01:10, , 5F
則對應到的 u0 (U中對應位置的行向量) 為
12/02 01:10, 5F

12/02 01:12, , 6F
u0 = 1/(奇異值)* A * v0 ,U中的行向量,跟你算A(A^T)
12/02 01:12, 6F

12/02 01:13, , 7F
的eigenvector一樣 ,小黃那本書裡面範例好像就這樣做的
12/02 01:13, 7F

12/02 01:14, , 8F
歹勢說的有點2266哈,不過大概做法就這樣
12/02 01:14, 8F

12/02 01:20, , 9F
v1 u1 對應λ1 多的特徵向量對到0
12/02 01:20, 9F

12/02 01:21, , 10F
一定可分解 因為他是對稱矩陣
12/02 01:21, 10F

12/02 01:21, , 11F
所以不用擔心重根找不到向量
12/02 01:21, 11F
我假設一個問題好了 ┌ ┐ A=│1 1│ │1 1│ │0 0│ └ ┘ T ┌ ┐ ┌ ┐ A A=│2 2│ 解得特徵根為4,0 ; 奇異值為2,0 Σ=│2 0│ │2 2│ │0 0│ └ ┘ │0 0│ └ ┘ 我知道奇異值需由大至小排,所以對應的V行向量一定會有順序 ┌ ┐ ┌ ┐ V(4)=span{│1│} V(0)=span{│ 1│} │1│ │-1│ └ ┘ └ ┘ ┌ ┐ V=│ 1 1 │ │--- ---│ │√2 √2│ │ │ │ 1 -1 │ │--- ---│ │√2 √2│ └ ┘ T ┌ ┐ AA =│2 2 0│ 解得特徵根為4,0,0 <---有重根 │2 2 0│ │0 0 0│ └ ┘ ┌ ┐ ┌ ┐ ┌ ┐ V(4)=span{│1│} V(0)=span{│ 1│,│0│} │1│ │-1│ │0│ │0│ │ 0│ │1│ └ ┘ └ ┘ └ ┘ ^^^^^^^^^^^^^^^^ 我的問題就是這個 我該如何決定我要先放哪個行向量在U中? ※ 編輯: jim055006 來自: 223.138.181.92 (12/02 13:27)

12/02 13:36, , 12F
隨便放嗎???!!!
12/02 13:36, 12F

12/02 14:18, , 13F
只要有對應到該對應的奇異值就可以,先後沒差
12/02 14:18, 13F

12/02 22:29, , 14F
我記得順序有差
12/02 22:29, 14F

12/02 23:06, , 15F
那到底要如何比較哩??
12/02 23:06, 15F

12/02 23:17, , 16F
AV=[AV1,AV2]=[2E1,0E2]=[E1,E2,E3]Σ=UΣ
12/02 23:17, 16F

12/02 23:20, , 17F
E2跟E3可換... 這是我的想法 有錯請指正了
12/02 23:20, 17F

12/02 23:37, , 18F
H大的意思是那兩個向量不管順序如何都會對應到兩個重根?
12/02 23:37, 18F

12/03 00:58, , 19F
恩 因為AVx=奇異值*Ex
12/03 00:58, 19F

12/03 12:51, , 20F
感謝!!!!!所有回答的人....鋼溫!!
12/03 12:51, 20F

09/11 14:38, , 21F
隨便放嗎???!!! https://daxiv.com
09/11 14:38, 21F
文章代碼(AID): #1ErwMtA5 (Grad-ProbAsk)
討論串 (同標題文章)
文章代碼(AID): #1ErwMtA5 (Grad-ProbAsk)