function [correlator] = gen_mtchd_fltr(Ns, m, h, QT)
phase = zeros(1,Ns);
correlator = zeros(m,m,m,Ns);
v = zeros(1,3);
%   phase(ak-2,ak-1,ak,t)
for ak2 = 1:m
    for ak1 = 1:m
        for ak = 1:m
            v = [(2*(ak2-1)-(m-1)),(2*(ak1-1)-(m-1)),(2*(ak-1)-(m-1))];
            phase = mod(  2*pi*h*( v(1)*QT(2*Ns+1:3*Ns) + v(2)*QT(Ns+1:2*Ns) + v(3)*QT(1:Ns) ), 2*pi );
            correlator(ak2,ak1,ak,:) = exp(j*phase);
        end;
    end;
end;