pwitrt.m

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

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