pwitrt.m

数值线性代数的Matlab应用程序包 共13个程序函数

function lmd=PwItrt(A,N)
[m,n]=size(A);
v(:,1)=eye(m,1);
for k=2:N+1
    w=A*v(:,k-1);
    v(:,k)=w/sqrt(w'*w);
    lmd(k-1)=v(:,k)'*A*v(:,k);
end