施密特正交化matlab代码

以下是施密特正交化的 MATLAB 代码：

function [Q,R] = schmidt(A)
% 施密特正交化
% A: m*n 矩阵，m>n
% Q: m*n 正交矩阵
% R: n*n 上三角矩阵

[m,n] = size(A);
Q = zeros(m,n);
R = zeros(n,n);

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

end

希望能对你有所帮助！