poweror.m

来自「自己编写的程序」· M 代码 · 共 24 行

function [A1]=powerOR(A,N)
A1=A;
[m,n]=size(A1);
R=zeros(n);
for t=1:N
bk(:,1)=A1(:,1);
Q(:,1)=bk(:,1)/sqrt(sum(bk(:,1).^2));
for k=2:n
    c=zeros(n,1);
    for j=1:k-1
        c=c+(A1(:,k)'*Q(:,j))*Q(:,j);
    end
    bk(:,k)=A1(:,k)-c;
    Q(:,k)=bk(:,k)/sqrt(sum(bk(:,k).^2));
end
    for k=1:n
       R(k,k)=sqrt(sum(bk(:,k).^2));
       for i=k+1:n
           R(k,i)=A1(:,i)'*Q(:,k);
       end
    end
    A2=R*Q;
    A1=A2;
end