精选优质文档-倾情为你奉上四.实验代码: 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
Copyright © 2018-2021 Wenke99.com All rights reserved
工信部备案号:浙ICP备20026746号-2
公安局备案号:浙公网安备33038302330469号
本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。