推 sl613:簡單來說就是故意把上三角拆成對角矩陣嗎? 10/09 19:55
→ alexndy:對的 10/09 20:49
推 sl613:感謝你的指導XD 本版有史以來最認真的一篇文 10/09 21:09
→ alexndy:XD 10/09 21:28
※ 引述《alexndy (保少爺)》之銘言：
: ※ 引述《sl613 (悠哉的午後)》之銘言：
: : 阿保 可以幫我解釋一下LDU分解的步驟嗎= =
: : 卡住了
: : 偽第三行
: 就先把他上三角化 (用高斯消去)
: 高斯消去法就等於對原來矩陣乘一個elementary matrix
: 所以 假設原來的是A
: (…)(E2)(E1)A=U
: U是上三角後的矩陣 Ek是elementary matrix 這些都會是下三角
: 然後再乘反矩陣過去 就是LU
: 再把對角線的再分出一個矩陣就可以
對角線如果不是都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
--
◢█◤▏
█ ▎
◤˙˙◥
◣ ▼ ◢
迷 你 邀 武 八 在 你 後 面 ， 他 非 常 火
Mini158 is behind you, she is furious.
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.211.228.36