yesf.m

来自「Matlab语言常用算法程序集.rar 里面有很多有用的Matlab数值运」· M 代码 · 共 27 行

function [x,IA]= Yesf(A,b)
N = size(A);
n = N(1);
A1 = A -eye(n,n);
A2 = zeros(n,n);
A3 = zeros(n,n);

for m=1:n
    for i=1:n
        for j=1:n
            if(i ~= m)
                A2(i,j) = A1(i,j);
            else
                if(j == m)
                    A2(i,j) = 1;
                else
                    A2(i,j) = 0;
                end
            end
        end
    end                                         %第一步递推
    A3 = A2 - A2(1:n,m)*A1(m,1:n)/(1+A1(m,m));  %第二步递推
    A1 = A3;
end
IA=A3;
x = A3*b;         %求解方程