Mercurial > hg > ddm
changeset 4:72c011ed1977 tip
more elaborate example with non-stat. estimate explanation
| author | smusevic |
|---|---|
| date | Tue, 30 Jul 2013 09:56:27 +0100 |
| parents | dde5034365b4 |
| children | |
| files | ddm_example_1.m |
| diffstat | 1 files changed, 17 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- 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); + +