新建 文本文档 (2).txt

来自「高斯塞得迭代和雅可比迭代的程序,对计算矩阵有帮助」· 文本 代码 · 共 39 行

function [x,k,flag]=Jacobic(A,b,delta,max1)
% 求解线性方程组的迭代法
% A为方程组的系数矩阵
% b为方程的右端项
% delta为精度控制
% max1 为最大迭代次数

if nargin<4 
    max1=100; 
end
if nargin<3 
    delta<1e-5;
end
n=length(A) ;
k=0 ;
x=zeros(n,1); 
y=zeros(n,1) ;
flag='OK!' ;
while 1
    for i=1:n
        y(i)=b(i) ;
        for j=1:n
            if j~=i
                y(i)=y(i)-A(i,j)*x(j) ;
            end
        end
        if abs(A(i,i))<1e-10 | k==max1
            flag='Fail' ;
            return ;
        end
        y(i)=y(i)/A(i,i) ;
    end
    if norm(y-x,inf)<delta
        break ;
    end
    x=y ;
    k=k+1 ;
end