jacdd.m

数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序

function X=jacdd(A,b,x0,P,wucha,max1)
[n,m]=size(A);
for j=1:m
    a(j)=sum(abs(A(:,j)))-2*(abs(A(j,j)));
end
for i=1:n
    if a(i)>=0
        disp('请注意：系数矩阵A不是严格对角占优的，此雅可比迭代不一定收敛')
        return
    end
end
if a(i)<0
    disp('请注意：系数矩阵A是严格对角占优的，此方程组有唯一解，且雅可比迭代收敛')
end
for k=1:max1
    k
    for j=1:m
        X(j)=(b(j)-A(j,[1:j-1,j+1:m])*x0([1:j-1,j+1:m]))/A(j,j);
    end
    X,djwcx=norm(X'-x0,P);xdwcx=djwcx/(norm(X',P)+eps);x0=X';X1=A\b;
    if (djwcx<wucha)&(xdwcx<wucha)
        disp('请注意：雅可比迭代收敛,此方程组的精确解jX和近似解X如下：')
        return
    end
end
if(djwcx>wucha)&(xdwcx>wucha)
    disp('请注意：雅可比迭代次数已经超过最大迭代次数max1')
end
a,jX=X1',X=X,