数值分析高斯塞德尔迭代MATLAB实现

%A为系数矩阵，b为等式右侧元素，x0为迭代初值，N代表最大迭代次数，p为误差允许范围
%其中需要输入x0、N、p
%需要输入nA，改变nA的值来改变方程组阶数
A=zeros(nA);
A(1,:)=-1;
A(:,1)=-1;
A(nA,:)=-1;
A(:,nA)=-1;
	for i=1:nA
	A(i,i)=30;
	end
A(1,nA)=0; 
A(nA,1)=0;
%以上构建方程组A的矩阵
nA=size(A,1);
b=zeros(nA,1);
for i=1:nA
    b(i,1)=2*i;
end
%以上构建方程组等式右侧矩阵b
n=length(A);
D=zeros(n);
x0=transpose(x0);  %转置
for i=1:n
    D(i,i)=A(i,i);
end
%构建D（A矩阵的对角线元素）矩阵
U=triu(-A)+D;
L=tril(-A)+D;
%构建对角线上部分元素矩阵U以及下部分元素矩阵L
k=0;
while k<=N
    x=inv(D-L)*U*x0+inv(D-L)*b;%迭代公式
    k=k+1;
    h=norm((x-x0),inf);
    if h<=p
        break
    end    
    x0=x;
end
if k>N
    disp(['迭代次数= ',num2str(k),',算法超过最大迭代次数。']);    
else
    disp(['迭代次数= ',num2str(k)]);    
end
x=x0;