施密特正交化QR分解matlab仿真

算法原理参考：https://blog.youkuaiyun.com/m0_66360845/article/details/137006605

function [Q, R] = QR_MGS(A)
%% A = (A1,A2,A3)，A1是列向量
[m,n] = size(A);
for j=1:n
   v=A(:,j);
   for i=1:j-1
       R(i,j)=Q(:,i)'*v;
       v=v-R(i,j)*Q(:,i);
   end
   R(j,j)=norm(v);
   Q(:,j)=v/R(j,j);
end

end