和大家分享一下另一種算法
(1).先將A列運算至對角方塊矩陣出現0矩陣(行運算亦可,看個人喜好)
R12(1) [ 1 2 0 0 ] [X O]
A---------->[-4 -4 0 0 ]= A_R =[ ]
R31(0.5) [ 0 0 2 0 ] [ ]
R42(-2) [ 2 2 0 1 ] [Y Z]
(2).
A_R=ER42(-2)*ER31(0.5)*ER12(1)*A
稍微整理之後 => A=ER12(-1)*ER31(-0.5)*ER42(2)*A_R
-1 -1
等號兩邊取inverse => A =A_R *ER42(-2)*ER31(0.5)*ER12(1)
(3).利用方塊矩陣求A_R^-1,詳細推導恕不詳述
-1 [ -1 ] [-1 -0.5 0 1 ]
A_R =[ X O ] =[ 1 0.25 0 -0.5 ]
[ -1 -1 -1 ] [ 0 0 0.5 0 ]
[-Z Y X Z ] [ 0 0.5 0 0 ]
-1
計算後可得A ,算了兩次終於跟原PO一樣XD
提出來參考參考,這方法也沒有比較快..囧
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.171.28.73