Matlab里实现施密特正交化

%斯密特正交化，其中a是一个191*191的矩阵
b = zeros(191,191);%创建一个容器，用于容纳a的每一列
b(:,1) = a(:,1);
for i = 2:191
    sum = zeros(191,1);
    for j = 1:i-1%求投影
        touyin = (dot(a(:,i),b(:,j))/dot(b(:,j),b(:,j)))*b(:,j);
        sum = sum + touyin;
    end
    b(:,i) = a(:,i) - sum;
end
%b就是a正交化后的矩阵