dan@0: function mfdi = multi_freq_ddm_int(frqs, mf_ders, sig, N, tol)
dan@0: %multi-frequency distribution derivative based estimator for
dan@0: % non-stationary sinusoidal analysis. Frequencies may not be at the bins of
dan@0: % FFT...
dan@0: %
dan@0: % Uses the Fourier kernel at specified frequencies. No assumptions about 
dan@0: %
dan@0: % [1] Michael Betser: Sinusoidal Polynomial Estimation Using The Distribution
dan@0: % Derivative, in IEEE Transactions on Signal Processing, Vol.57, Nr. 12,
dan@0: % December 2009
dan@0: %
dan@0: % 
dan@0: %
dan@0: 
dan@0: 
dan@0: R = size(frqs,2);
dan@0: t = [-(N-1)/2:(N-1)/2]'/fs;
dan@0: 
dan@0: t_mat = repmat(t,[1,R]);
dan@0: frq_mat = repmat(frqs,[N,1]);
dan@0: krnls = exp(1j .* t_mat .* frq_mat);
dan@0: krlns_ders = 1j * krnls .* frq_mat;
dan@0: 
dan@0: mfdi = gen_ddm_int(krnls, krlns_ders, mf_ders, sig, N, tol);
dan@0: 
dan@0: end