mgsqr.m

来自「一个二阶中心插值滤波器（DDF2）的目标跟踪仿真实例。」· M 代码 · 共 14 行

function [R]=MGSQR(A)
[m,n]=size(A);
Q=zeros(m,n);
R=zeros(n,n);
for i=1:n
    R(i,i)=norm(A(:,i));
    Q(:,i)=A(:,i)/R(i,i);
    if i<n
    for j=i+1:n
        R(i,j)=Q(:,i)'*A(:,j);
        A(:,j)=A(:,j)-Q(:,i)*R(i,j);
    end
    end
end