anld.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 19 行

function [h,q]=Anld(A)
[m,n]=size(A);
q=zeros(m,n);
h=zeros(m,n);
q(1,1)=1;
for i=1:m
    v=A*q(:,i);
    for j=1:n
        h(j,i)=q(:,j)'*v;
        v=v-h(j,i)*q(:,j);
    end
    h(i+1,i)=sqrt(v'*v);
%     if h(i+1,i)<1E-10&&h(i+1,i)>-1E-10
%         break
%     end
    q(:,i+1)=v/h(i+1,i);
end
% h=h(1:i,1:i);