dan@0: function [fk fk_ders] = fourier_kernel_win(frqs, N,fs, win, win_der) dan@0: R = size(frqs,2); dan@0: t = [- (N - 1) / 2:(N - 1) / 2]' / fs; dan@0: win_mat = repmat(win, [1,R]); dan@0: win_der_mat = repmat(win_der, [1,R]); 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: fk = win_mat .* krnls; dan@0: fk_ders = (1j * win_mat .* frq_mat + win_der_mat) .* krnls; dan@0: end