作者alexndy (保少爺)
看板SMJH90308
標題Re: 呼叫保言~呼叫保言
時間Fri Oct 17 00:08:19 2008
※ 引述《sl613 (悠哉的午後)》之銘言:
: 唔...又有問題了>"<
: 如果要把A矩陣轉成L*(L轉) 有沒有標準步驟?
: 到底線代哪裡比較重要阿(泣)
: ※ 引述《alexndy (保少爺)》之銘言:
: : 對角線如果不是都1的話:
: : [2 3 4] [2 0 0][1 3/2 2]
: : [0 3 4] = [0 3 0][0 1 4/3]
: : [0 0 4] [0 0 4][0 0 1]
: : 如果是下三角的情況也是類似
: : 可以把diagonal的項再分出來
: : 然後使得LDU裡的L和U的對角線項都是1
ex:
[6 2 2] [1 0 0][6 2 2] [1 0 0][6 0 0][1 1/3 1/3]
[2 2 1] = [1/3 1 0][0 4/3 1/3] = [1/3 1 0][0 4/3 0][0 1 1/4]
[2 1 1] [1/3 1/4 1][0 0 1/4] [1/3 1/4 1][0 0 1/4][0 0 1]
A L U L D
U=L'
把D分成兩個根號D
然後再分別乘到L跟U就可以
A=(L*D^0.5)*(D^0.5*L')=(L*D^0.5)*(L*D^0.5)' 對角矩陣轉置還是對角
L L'
注意:討論是不是正定矩陣前提那個矩陣就是對稱矩陣了
--
◢█◤▏
█ ▎
◤˙˙◥
◣ ▼ ◢
迷 你 邀 武 八 在 你 後 面 , 他 非 常 火
Mini158 is behind you, she is furious.
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.211.228.36
→ sl613:跨嚨某..... 10/17 19:50
→ K0F:這矩陣似乎也太矩了吧!! 10/18 11:04
→ alexndy:就…把D分成兩個根號D 就是重點 10/19 23:49