作者PanJC ((#‵Д′)f〒﹌﹌﹌﹌﹌ꄩ
看板NCTU-STAT100
标题[闲聊] gaus
时间Thu Jul 28 15:24:26 2011
Step 1: for k= 0 : n-1
do Step 2 ~ Step 6
Step 2: (is, js[k]) = arg max{ |a(i,j)| | i = k:n-1, j = k:n-1 }
(i,j)
Step 3: if( |a(is,js[k])| < eps ) then stop
else do Step 4 ~ Step 6
Step 4: if( is != k or js[k] != k ) then
a(k,:) <=> a(is,:)
a(:,k) <=> a(:,js[k])
b(k) <=> b(is)
Step 5: a(k,k) = 1.0 / a(k,k)
for( j= k+1 : n-1)
do a(k,j) = a(k,j) * a(k,k)
b(k) = b(k) * a(k,k)
Step 6: for i = k+1 : n-1
for j = k+1 : n-1
a(i,j) = a(i,j) - a(i,k) * a(k,j)
for i = k+1 : n-1
b(i) = b(i) - a(i,k) * b(k)
Step 7: x(n-1) = b(n-1)
for i= n-2 : 0
x(i) = b(i) - sum_{j=i+1}^{n-1} a(i,j) * x(j)
Step 8: for( k = n-1 : 0 )
x(k) <=> x(js[k])
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.7.248
※ 编辑: PanJC 来自: 140.113.7.248 (07/28 15:34)
※ 编辑: PanJC 来自: 118.233.7.208 (07/28 20:23)
※ 编辑: PanJC 来自: 118.233.7.208 (07/28 20:24)
1F：→ c76068:有点不懂...step4 为什麽会有is[k]? 07/29 01:37
打错, 已改过
※ 编辑: PanJC 来自: 118.233.7.208 (07/29 06:14)
2F：→ c76068:最後只要印出解就好了吗? 07/30 00:18
3F：→ PanJC:yes. 07/30 10:32
4F：推 gentlecat:step1 是不是应该改成做 step2~5阿?? 07/30 17:40
5F：推 gentlecat:喔喔没事我蠢了= = " 07/30 17:47
6F：推 gentlecat:请问一下有正确答案吗@@? 07/31 20:21
7F：推 WhywhyTseng:我直接手动带回去算耶XDDD 07/31 20:39
8F：→ PanJC:手算啊?! 写个程式代回去验算可能比较快.... 07/31 20:49
9F：推 c76068:为什麽不直接用matlab算? 或者跑一个没有截距项的回归..也 07/31 21:05
10F：→ c76068:可以得到答案....你有统计套装软体吧.... 07/31 21:05
11F：→ PanJC:当然可以, 但前提假设是同学们已经至少会一样统计软体. 07/31 21:12
12F：→ PanJC:如果都不会, 至少目前还有 C 可以写个程式验算. 07/31 21:12
13F：推 smallstrong2:推荐那个用手动的XDD 08/01 10:59
14F：推 gentlecat:手动+1XDD 08/01 11:17