※ 引述《aliking (aliking)》之銘言:
: 我想請問一個vb的問題,
: 也就是求解一個Ax=b的問題,這個解矩陣的問題,
: 應該會有高手已經寫好模組吧?
: 請問有人知道要去那找這個模組呢?
: 一般好像是用高斯消去…
: 不過我的問題中的A應該是幾百乘幾百的矩陣!
: 用高斯好像…不太適合!
: 請各位幫忙一下! ALIKING
Problem Set:(我只會當A為方陣的case)
如果A是n乘n的矩陣, x與b分別是n乘1的向量, A與b為已知, 試求向量x。
------------------------------------------------------------------
Step1. 將A化成上三角矩陣(利用row operator)
Step2. 之後由下而上, 求得x向量的每一個element的值
--------------------------------------------------------------------
For example:
1 2 3 x1 1
Ax=b為 ﹝4 5 6﹞﹝x2﹞=﹝2﹞
7 8 9 x3 3
1 2 3 1 2 3
Step1: 把A化為上三角矩陣 => A = ﹝4 5 6﹞= ﹝0 -3 -6﹞
7 8 9 0 0 -3
1 2 3 x1 1
Step2: ﹝0 -3 -6﹞﹝x2﹞=﹝2﹞
0 0 -3 x3 3
=> x3=-1,
x2=(-3)*(x2)+(-6)*(x3)=2 => x2=4/3
x1=4/3 (同x2之算法)
如果你要Step1.的src, 下一篇會post
--
﹏﹏﹏﹏﹏ ▁▆▅▄▆█▂ ︴︴︴︴︴ ◎▁▁█▁▁◎
﹌ ﹌ \ / ╭ ╮ ︵ ︵
∩ ∩ ●─● ︶ ︶ ◎─◎
﹤ ︷ 卅 ︰︿: ╭ˍˍ╮
╰╯ ∩ ︹ ╲╱
喜 怒 哀 樂 collegeman.
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.59.31.47
※ 編輯: collegeman 來自: 61.59.31.47 (05/20 01:51)
※ 編輯: collegeman 來自: 61.59.31.47 (05/20 01:52)
※ 編輯: collegeman 來自: 61.59.31.47 (05/20 01:54)
※ 編輯: collegeman 來自: 61.59.31.47 (05/20 01:54)