※ 引述《lalalarc (初來乍到)》之銘言：
: ※ 引述《roussas (QQ)》之銘言：
: : http://www.cyut.edu.tw/~ckhung/b/la/gje.shtml
: : 這裡也有高斯jordan的做法...
: 如果是用 Gauss-Seidel 法 必須要考慮到數列收斂的問題
: 因為用這方法的前提是可以找到一初始值使得n至無窮時x會收斂到真正的值
: 所以就算方程組有解 用此法可能會成數列發散而找不到解
: 如果掉換一下各方程式的位置 重新寫一組矩陣 或許就可收斂了
: 或是改變初始值
: ex:
: x+2y+4z=1
: 3y+z=1
: 2x+y+3z=1 (x,y,z)=(1/3,1/3,0)
: 用Gauss-Seidel法用不出來
: 但將第一列第三列交換數列就收斂了........
好像有一個條件
fabs(ajj) > sum( fabs(aij) ,for any i != j )
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.7.249