作者PanJC ((#‵Д′)f〒﹌﹌﹌﹌﹌ꄩ
看板NCTU-STAT100
标题[闲聊] HINT: inverse matrix
时间Thu Jul 28 21:18:52 2011
Step 1: for k= 0 : n-1
do Step 2 ~ Step 6
Step 2: (is[k], js[k]) = arg max{ |a(i,j)| | i = k:n-1, j = k:n-1 }
(i,j)
Step 3: if( | a(is[k], js[k]) | < eps ) then stop
else do Step 4 ~ Step 6
Step 4: if( is[k] != k ) then
a(k,:) <=> a(is[i],:)
if( js[k] != k ) then
a(:,k) <=> a(:,js[k])
Step 5: a(k,k) = 1.0 / a(k,k)
for j = 0 : n-1 ; j != k // 注意这里的 j != k 意指 j 要跳过 k
do a(k,j) = a(k,j) * a(k,k)
Step 6: for i = 0 : n-1 ; i != k
for j = 0 : n-1 ; j != k
a(i,j) = a(i,j) - a(i,k) * a(k,j)
for i = 0 : n-1 ; i != k
a(i,k) = - a(i,k) * a(k,k)
Step 7: for k = n-1 : 0
for l = 0 : n-1
a(k,l) <=> a(js[k],l)
a(l,k) <=> a(l,is[k])
请同学们多多了解为什麽这个程式能 work, 在使用它时, 信心就会多一点.
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 118.233.7.208
1F：推 gentlecat:推荐这个认真的学长XD 07/28 22:58
2F：→ PanJC:大家都很认真, 加油 07/29 06:28
3F：推 sukotoko:谢谢学长:) 07/30 20:25
4F：推 gentlecat:可以问一下STEP5~7的意思吗? 有点看不太懂@@ 08/01 00:56