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