※ 引述《roussas (QQ)》之銘言：
: : http://www.na.ntu.edu.tw/chinese/course/online/num_anal/chap2.htm
: : 這裡有演算法 不知道有沒有幫助
: : 不然~~去圖書館借
: : 數值方法入門-使用C語言- 全華科技圖書
: : 一共有兩本
: : 裡面有C寫的程式跟講解 Jacobi Method的作法
: : 星期三我借的時候還有一本 一本我借走了
: : 希望對你有幫助
: 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法用不出來
但將第一列第三列交換數列就收斂了........
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.185.114