matlab做施密特正交

for i=1:tiqvduan
if i==1
shimi=[shimi,duan(:,i)]
else
tou=duan(:,i)'*shimi
mo=shimi.*shimi
mo=sum(mo,1)
%mo=sqrt(mo)
k=tou./mo
k=repmat(k,[h 1])
alltou=k.*shimi
alltou=sum(alltou,2)
shimi=[shimi,duan(:,i)-alltou]
end
end
shimo=repmat(sqrt(sum(shimi.*shimi,1)),[h 1])
shimi=shimi./shimo