jacobi.m

来自「jacobi迭代法程序和矩阵的QR分解算法」· M 代码 · 共 20 行

function x=jacobi(A,b)
[m,m]=size(A);
x=zeros(m,1);
D=diag(diag(A));
N=100;
for i=1:m
    L=zeros(m,m);
    U=zeros(m,m);
    for j=1:i
        L(i,j)=D(i,j)-A(i,j);
    end
    for j=i:m
        U(i,j)=D(i,j)-A(i,j);
    end
end
while i<=N
    for k=2:m
        x(k,1)=inv(D)*(L+U)*x(k-1,1)+inv(D)*b;
    end
end