qrfac.m

来自「Numerical analysis （Gauss、Bisection、jaco」· M 代码 · 共 25 行

function s=QRfac(a,N)
for k=1:N
    for i=1:5
        for j=i+1:6
            d=sqrt(a(i,i).*a(i,i)+a(j,i).*a(j,i));
            disp(d);
            if d==0
                continue;
            end
            c=a(i,i)./d;
            s=a(j,i)./d;
            Q=eye(6);
            Q(i,i)=c;
            Q(j,j)=c;
            Q(i,j)=s;
            Q(j,i)=-s;
            %disp(Q);
            R=inv(Q);
            a=Q*a*R;
            %disp(a);
        end
    end
end
s=a;
return