Re: [理工] [工數] 矩陣代數運算

看板Grad-ProbAsk作者 (Bow)時間14年前 (2011/08/03 00:06), 編輯推噓1(102)
留言3則, 2人參與, 最新討論串2/2 (看更多)
綜合推文中大大們的想法,以及自己的算法 希望沒有漏氣XD 4 -1 1.AB=2I,由A 仍為方陣可看出A為方陣,得B=2A 4 [a b] 2.令C = A = [ ],(由於C為實對稱矩陣,因此可將C對角化) [b a] let (C-λI)u = 0 ,可得eigenvalues:(a+b)、(a-b) [1] [ 1] eigenvectors:[1]、[-1] 4 [1 1][a+b 0 ][1 1]-1 -1 C = A = [ ][ ][ ] = UDU [1 -1][ 0 a-b][1 -1] 1/4 1/4 1/4 1/4 1/4 1/4 -1 [(a+b) + (a-b) (a+b) - (a-b) ] then A=C =UD U =0.5*[ ] [ 1/4 1/4 1/4 1/4] [(a+b) - (a-b) (a+b) + (a-b) ] -1 1/4 -1 -1 4 -1 B=2A =2(UD U ) =2UD U [ 4 4 4 4] =[(a+b) + (a-b) (a+b) - (a-b) ] [ 4 4 4 4] [(a+b) - (a-b) (a+b) + (a-b) ] 呼,終於打完了,有錯麻煩指正:)) -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.114.207.10

08/03 00:12, , 1F
應該沒錯 但我之前是把U是以歸一正交向量組成
08/03 00:12, 1F

08/03 00:13, , 2F
U^(-1) = U^(T) 這樣不知道有無較快
08/03 00:13, 2F

08/03 00:40, , 3F
這個方法OK的,不過考慮到最後還要展開.所以就沒這樣做
08/03 00:40, 3F
文章代碼(AID): #1EE21zJd (Grad-ProbAsk)
文章代碼(AID): #1EE21zJd (Grad-ProbAsk)