上传者: m0_52957036
|
上传时间: 2021-12-07 18:51:06
|
文件大小: 13KB
|
文件类型: -
function [Q,R,x]=QR(A,b) n=length(A;Q=eye(n;R=A; for i=1:n-1 x=R(i:n,i; if all(x(2:n-i+1)~=0) x=x/norm(x,inf;delta=norm(x,2*sign(x(1; v=[x(1)+delta;x(2:n-i+1)];beta=2/norm(v,2^2; else beta=0;v=0; end