[線代] 如何將一個nonsymmetric矩陣拆解成兩個symmetric?

[letibe (remember the fate)]

因為程式需求要找出一個non symmetric matrix 的 Jordan form

雖然matlab可以很快的找出來

不過老師要求要用SAS/IML跑

可是SAS並沒有相關的指令

他唯一的功能是

一個nonsymmetric matrix A 可以拆成 symmetric 的B，C

且A= inverse(B) * C

這時候藉由B和C可以找出generalized vector

我想用這樣來組成Jordan form

可是問題來了....要怎樣才能拆解B和C呢?

希望哪位前輩給點意見

只有數學式也好~ 我想辦法弄成code看看

感謝

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.119.144.53
※ 編輯: letibe          來自: 140.119.144.53       (04/07 18:36)
推 powerkshs   :只會用手XD