Matlab实现雅各比矩阵

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/22 18:27:03
Matlab实现雅各比矩阵Matlab实现雅各比矩阵Matlab实现雅各比矩阵functionX=jacobi(A,b,P,delta,max1)%A是n维非奇异阵%B是n维向量%P是初值%delta

Matlab实现雅各比矩阵
Matlab实现雅各比矩阵

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