数值分析QR方法求矩阵特征值和特征向量(共4页).doc

精选优质文档-倾情为你奉上四.实验代码： function H,B=Hessenberg(A)n=length(A);B=eye(n);for k=1:n-2 X=zeros(n-k,1); H=eye(n); for i=1:n-k X(i)=A(i+k,k); end a=max(abs(X); if a=0.0 break end X=X/a;c=X(1); b1=sqrt(sum(X.2); if X(1)=0 b1=-b1;end X(1)=X(1)-b1; b=b12-b1*c; H0=eye(n-k)-X*X/b; for i=1:n-k for j=1:n-k H(i+k,j+k)=H0(i,j); end end A=H*A*H; B=B*H;en