基于Givens旋转完成QR分解进而求解实矩阵的逆矩阵-MATLAB代码

function [Q,R] = fx_GivensRotation_QR_decom(A) N=size(A,1); Q=eye(N); for j=1:N %列? for i=j+1:N %行? a=A(j,j); b=A(i,j); if(b~=0) %想消为0的位置不为0才进行消除? c=a/sqrt(a^2+b^2); s=b/sqrt(a^2+b^2); G=eye(N,N); %Givens矩阵,用第j列第j个元素去消第j列第i个元素 G(j,j)=c; G(j,i)=s; G(i,j)=-s; G(i,i)=c; A1=G*A; R=A1; Q=Q*G.'; A=A1; end end end end