怎樣用高斯-賽德爾迭代法求解矩陣方程？%Gauss-Seidel函數(shù)%a是未知數(shù)的系數(shù)矩陣；B是方程的右常數(shù)列向量；x0取1（m，1）m作為未知數(shù)；EPS是精度；如果不輸入EPS，則默認(rèn)為1.0e-6

怎樣用高斯-賽德爾迭代法求解矩陣方程？

%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

如果使用Jacobi或Gauss-Seidel型迭代方法，則需要重新排列系數(shù)矩陣將非零元素排列在對(duì)角線上