3.0000 2.0000 1.0000 3.0000 2.0000 1.0000 3.0000 2.0000 1.0000 3.0000 2.0000 1.0000 ans = 3.0000 2.0000
1.0000
解的迭代情况
1.5
0.5
图一。解的迭代过程
程序说明:xhis为迭代过程,ans为最终迭代结果.本次迭代次数为10.图一中红点标记每次迭代的解的轨迹。从图中可以看出,初始值选择在(0,0,0),经过几次就应经到达很高的精度,表明Gauss—Seidel的迭代法收敛很快。
三.MATLAB程序
function y=GS(A,b) if(nargin==0)
disp('你没有输入任何参数');
A=ceil(10*rand(3,3));b=ceil(10*rand(3,1)); A,b GS(A,b)