矩阵的QR分解基于施密特正交化

function [qr]=qr（a）
% QR decomposition: A=QR.
% Decompose rectangqlar matrix A to a  orthogonal matrix Q a uppertriangle
% matrix R by means of Gram-Schmidt orthogonalization.
  [mn]=size（a）;  %计算矩阵A的大小
%------------给各矩阵赋初值------------------------------
   q=zeros（mn）; % q--正交基构成的矩阵
   r=zeros（mn）; % r--记录正交化过程中各元素的上三角矩阵
   v=zeros（n1）; % v--存储正交化过程中的列向量的范数
   vv=zeros（m1）;% vv 和 vs 均为过度列向量
   vs=zeros（m1）;
%----------Gram-Schmidt orthogonalization---------------
% 计算第一个正交基及其范数：q1和v1
   for j=1:m      
        v（1）=v（1）+a（j1）.^2;
   end
   r（11）=sqrt（v（1））;
   for j=1:m