%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)*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;
数值分析雅可比迭代MATLAB代码
于 2023-10-11 10:10:26 首次发布