作业帮Matlab实现雅各比矩阵
来源：学生作业帮 编辑：拍题作业网作业帮 分类：综合作业 时间：2024/04/28 07:10:19
Matlab实现雅各比矩阵
function X=jacobi(A,b,P,delta,max1)
% A是n维非奇异阵
% B是n维向量
% P是初值
% delta是误差界
% X为所求的方程组AX=B的近似解
N=length(b);
for k=1:max1
for j=1:N
X(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);
end
err=abs(norm(X'-P));
P=X';
if (err<delta)
break
end
end
X=X';k,err
>> A=[4,1,-1;1,-5,-1;2,-1,-6]
>> b=[13;-8;-2]
>> P=[0;0;0]
>> X=jacobi(A,b,P,10^(-4),20)
k =
9
err =
2.5713e-005
X =
3.0000
2.0000
1.0000
% A是n维非奇异阵
% B是n维向量
% P是初值
% delta是误差界
% X为所求的方程组AX=B的近似解
N=length(b);
for k=1:max1
for j=1:N
X(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j);
end
err=abs(norm(X'-P));
P=X';
if (err<delta)
break
end
end
X=X';k,err
>> A=[4,1,-1;1,-5,-1;2,-1,-6]
>> b=[13;-8;-2]
>> P=[0;0;0]
>> X=jacobi(A,b,P,10^(-4),20)
k =
9
err =
2.5713e-005
X =
3.0000
2.0000
1.0000