Mercurial > hg > ddm

diff stat_pars_ddm.m @ 0:a4a7e3405062

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/stat_pars_ddm.m	Thu May 09 20:04:15 2013 +0100
@@ -0,0 +1,11 @@
+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
\ No newline at end of file