view stat_pars_ddm.m @ 2:e2b116f3b69b

fixed reallocate and pinv_2_fast to work on Octave
added short example very quickly, might not be working properly

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