作者shyfang (shyfang)
看板NCTU-STAT94G
标题Re: 不会是这样的吧??? orz
时间Mon Aug 1 13:02:52 2005
※ 引述《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