formq.m

来自「麻省理工ComputationalMethod」· M 代码 · 共 8 行

function Q=formQ(W)

[m,n]=size(W);
Q=eye(m);
for k=n:-1:1
  Q(k:m,:)=Q(k:m,:)-2*W(k:m,k)*(W(k:m,k)'*Q(k:m,:));
end