如果您无法下载资料，请参考说明：
1、部分资料下载需要金币，请确保您的账户上有足够的金币
2、已购买过的文档，再次下载不重复扣费
3、资料包下载后请先用软件解压，在使用对应软件打开
function[Q,R]=CQR(A)%对矩阵A进行QR分解[m,n]=size(A);Q=eye(n);fori=1:na=A(i:n,i);[~,s]=max(abs(a));e=zeros(length(a),1);e(1)=1;I=eye(length(a));b=sign(a(s))*(norm(a));u=a+b.*e;w=u/(norm(u));h=I-2*w*(w');H=zeros(n,n);fork=1:iH(k,k)=1;endforj=i:nfort=i:nH(j,t)=h(j-k+1,t-k+1);endendA=H*A;Q=H*Q;endQ=Q';R=A;
立即下载
