annotate stat_pars_ddm.m @ 4:72c011ed1977
tip
more elaborate example with non-stat. estimate explanation
author |
smusevic |
date |
Tue, 30 Jul 2013 09:56:27 +0100 |
parents |
a4a7e3405062 |
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 |