Schmidt

A为矩阵
[m,n]=size(A);

Q=zeros(m,n);

tQ=zeros(m,n);

R=zeros(n,m);

for j=1:n
    v=A(:,j);
    for i=1:j-1
        R(i,j)=(tQ(:,i)'*A(:,j))/(norm(tQ(:,i))^2)
        v=v-R(i,j)*tQ(:,i);
    end
    R(j,j)=norm(v);
    tQ(:,j)=v;
    Q(:,j)=v/R(j,j);
end

SMatrix = Q;

 

 大家试试 怎么样的  嘻嘻