※ 引述《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