jacobi_g.asv

来自「数值分析常用几种方法」· ASV 代码 · 共 34 行

function X=jacobi_G(A,n)
%A为非奇异矩阵

D=zeros(n);
L=zeros(n);
U=zeros(n);
for i=1:n
    for j=1:n
        if i==j
            D(i,j)=A(i,j);          
        end
        if j<i
        L(i,j)=A(i,j);
    end
    end
end
for i=2:n
   for j=1:n
    if j<i
        L(i,j)=A(i,j);
    end
   end
end
for i=1:n-1
   for j=1:n
    if j>i
        U(i,j)=A(i,j);
    end
   end
end
G=-inv(D)*(L+U);
G