Mercurial > hg > ddm
view stat_pars_ddm.m @ 1:462de60d7e94
missing files
author | Saso Musevic musevic@gmail.com |
---|---|
date | Thu, 25 Jul 2013 10:40:55 +0100 |
parents | a4a7e3405062 |
children |
line wrap: on
line source
function [ests r0] = stat_pars_ddm(sig,gdi,t,N,w) frq_est = imag(gdi(1)); ply_est = fliplr(gdi.'); ply_est(end) = real(ply_est(end)); % hm.... wmp = win_mod_poly(N, t, [ply_est 0], w); krnl_dmd = exp(1j * frq_est * t'); dtft_dmd = sum(w .* sig .* conj(krnl_dmd),1); r0 = dtft_dmd / wmp * 2; %ply_est(end) = ply_est(end) + 1j*frq_est; ests = [fliplr(gdi.') abs(r0) + 1j*angle(r0)]; end