QR分解解方程的MATLAB程序

function[Q,R,x]=QR(A,b)n=length(A);Q=eye(n);R=A;fori=1:n-1x=R(i:n,i);ifall(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;elsebeta=0;v=0;endh=eye(n-i+1)-beta*v*v';H=blkdiag(eye(i-1),h);Q=Q*H;R=H*R;endR=triu(R);y=Q'*b;fori=n:-1:1x(i)=(y(i)-R(i,i+1:n)*x(i+1:n))/R(i,i);end%注解：Ax=b，A=QR。