Mercurial > hg > ddm
annotate stat_pars_ddm.m @ 0:a4a7e3405062
Import DDM code by Sašo Muševič
| author | Dan Stowell <dan.stowell@elec.qmul.ac.uk> |
|---|---|
| date | Thu, 09 May 2013 20:04:15 +0100 |
| parents | |
| children |
| rev | line source |
|---|---|
| dan@0 | 1 function [ests r0] = stat_pars_ddm(sig,gdi,t,N,w) |
| dan@0 | 2 frq_est = imag(gdi(1)); |
| dan@0 | 3 ply_est = fliplr(gdi.'); |
| dan@0 | 4 ply_est(end) = real(ply_est(end)); % hm.... |
| dan@0 | 5 wmp = win_mod_poly(N, t, [ply_est 0], w); |
| dan@0 | 6 krnl_dmd = exp(1j * frq_est * t'); |
| dan@0 | 7 dtft_dmd = sum(w .* sig .* conj(krnl_dmd),1); |
| dan@0 | 8 r0 = dtft_dmd / wmp * 2; |
| dan@0 | 9 %ply_est(end) = ply_est(end) + 1j*frq_est; |
| dan@0 | 10 ests = [fliplr(gdi.') abs(r0) + 1j*angle(r0)]; |
| dan@0 | 11 end |