高斯迭代法例題 怎樣用高斯-賽德?tīng)柕ㄇ蠼饩仃嚪匠蹋?/h1>
怎樣用高斯-賽德?tīng)柕ㄇ蠼饩仃嚪匠蹋?Gauss-Seidel函數(shù)%a是未知數(shù)的系數(shù)矩陣;B是方程的右常數(shù)列向量;x0取1(m,1)m作為未知數(shù);EPS是精度;如果不輸入EPS,則默認(rèn)為1.0e-6
怎樣用高斯-賽德?tīng)柕ㄇ蠼饩仃嚪匠蹋?/h2>
%Gauss-Seidel函數(shù)%a是未知數(shù)的系數(shù)矩陣;B是方程的右常數(shù)列向量;x0取1(m,1)m作為未知數(shù);EPS是精度;如果不輸入EPS,則默認(rèn)為1.0e-6函數(shù)[x,n]=Gauss-Seidel(a,B,x0,EPS)如果nargin==3 EPS=1.0e-6elseif nargin<3 error Return d=diag(diag(a))U
在Gauss-Seidel迭代中,a=d-l-U,但是AX=B(d-l)x=UX,bx=(d-l)^{-1}UX(d-l)^{-1}B
so B=(d-l)^{-1}U