qrfa.asv

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

function s=QRfa(a,N)
Q=eye(4);
for k=1:N
    for i=1:3
        for j=i+1:4
            d=sqrt(a(i,i).*a(i,i)+a(j,i).*a(j,i));
            disp(d);
            c=a(i,i)./d;
            s=a(j,i)./d;
            Q(i,i)=c;
            Q(j,j)=c;
            Q(i,j)=s;
            Q(j,i)=-s;
            %disp(Q);
           % R=inv(Q);
            s=Q*a;
            disp(s);
        end
    end
end
return