Subsample Delay Estimation: y=fourier_series(Y,t); - File Exchange

Code covered by the BSD License

Highlights from
Subsample Delay Estimation

[e f_e G P]=zzb(s,n1,n2, f_s,... calculates the modified Ziv-Zakai lower bound for time delay estimation.
[x_max y_max A]=crit_interp_p...
delay=delayest_3point(u2,u1,m... delay=delayest_3point(u2,u1,method,estimator,parameter);
delay=delayest_fft(y,u); delay=delayest_fft(y,u);
delayest_iterative(u2,u1,N_i_...
x=rand_white(N);
y=fft_circshift(u,d) y=fft_circshift(u,d)
y=fourier_series(Y,t);
compare_delay_methods_6_03.m
simple_example.m
View all files

from Subsample Delay Estimation by Travis Wiens
Demonstrates a number of methods of estimating the delay between two signals to subsample accuracy.

y=fourier_series(Y,t);

function y=fourier_series(Y,t);
%y=fourier_series(Y,t);
%Calculates the fourier series (i.e. the inverse fft for a single point,
%not necessarily on an integer sample).
%Assumes the ifft(Y) is real and has an even number of points.

%    Copyright Travis Wiens 2009 travis.mlfx@nutaksas.com

N_Y=numel(Y);
C_n=1/N_Y;

y=C_n*real(Y(1))*ones(1,numel(t));%mean
y=y+2*C_n*real(Y(2:(N_Y/2))*exp(2j*pi*C_n*(1:(N_Y/2-1))'*t));
%for real valued output, Y above Nyquist is complex conjugate
y=y+C_n*real(Y(N_Y/2+1)*exp(2j*pi*C_n*N_Y/2'*t));%Nyquist frequency

Highlights from Subsample Delay Estimation

Highlights from
Subsample Delay Estimation