# HG changeset patch
# User smusevic
# Date 1375174587 -3600
# Node ID 72c011ed1977b5894b793a0ef64149035544e299
# Parent  dde5034365b48e8a0d30763e1deb0fafde5a5c06
more elaborate example with non-stat. estimate explanation

diff -r dde5034365b4 -r 72c011ed1977 ddm_example_1.m
--- a/ddm_example_1.m	Thu Jul 25 13:07:34 2013 +0100
+++ b/ddm_example_1.m	Tue Jul 30 09:56:27 2013 +0100
@@ -34,5 +34,21 @@
 [A, b, sig_fft] = ddm_lin_sys_fft(nr_bins, ddm_dgr, w', wd', poly_model', bffr, N, N_fft, fs);
 % fast pseudo-inverse
 ddm_estimates = flipud(pinv_2_fast(A,b,nr_bins,N-nr_bins+1).');
+
+% phase function estimates
+ddm_frequency_estimates_for_each_fft_bin = imag(ddm_estimates(2,:));
+ddm_linear_fm_estimates_for_each_fft_bin = imag(ddm_estimates(1,:));
+% amplitude funciton estimates
+ddm_linear_am_estimates_for_each_fft_bin = real(ddm_estimates(2,:));
+ddm_2nd_dgr_am_estimates_for_each_fft_bin = real(ddm_estimates(1,:));
+
 % re-assign the spectral magnitude 
-reallocate(abs(sig_fft), ddm_estimates, N_q, fs);
+ddm_spectrogram = reallocate(abs(sig_fft), ddm_estimates, N_q, fs);
+
+% find the peak
+[peak_val peak_bin] = max(ddm_spectrogram);
+
+% frequency trajectory of the peak bin
+frequency_trajectory = ddm_linear_fm_estimates_for_each_fft_bin(peak_bin) * t + ddm_frequency_estimates_for_each_fft_bin(peak_bin);
+
+