FFT time delay, incorrect alighnment of the output wave

1 visualización (últimos 30 días)
Sophie
Sophie el 15 de Sept. de 2013
Hi,
I shift the FFT of a sampled signal by a certain number of samples and take the inverse FFT to get a delayed signal. I use the formula:
x[n+k] -> FFT(x[n]) * exp(-j*2*pi*k*m/N), where N = the length of the signal, k is the amount of time delay in time domain and m is the running index.
When I take the IFFT, I see that the signal starts correctly but then it gets misaligned (For example, in Figure 1 you can see there are 2 peaks in place of one peak at t~70 s). Why am I getting this result?
My code is below:
%for fs=2*10^8 Hz and f_signal=300 MHz
Fs=300000000; %fs=300MHz
T=1/Fs; %the sampling period
t=0:T:((0.7*10^-6)-T);
L=length(t);
f_signal=3*10^6;
x=cos(2*pi*f_signal*t); % signal_T=1/(3*10^6); %= 3.3333e-07
X=fft(x);
delayed_samples= 126; %delay the signal by 126 samples:
Y=X.*exp(1j*2*pi*delayed_samples*[0:L-1]/L );
y_1=real(ifft(Y));
figure, plot(Fs.*t, x, 'g'), hold on, plot(Fs.*t, y_1, 'r'), legend ('x','delayed x') , title 'signals vs. samples'
%%%end of the script%%%%%%%
Figure 1 is

Respuestas (1)

Matt J
Matt J el 16 de Sept. de 2013
Editada: Matt J el 16 de Sept. de 2013
Because it is a circulant shift. The discontinuity between x(1) and x(end) has been shifted circulantly to around frequency=83.
  2 comentarios
Sophie
Sophie el 18 de Sept. de 2013
If I need them to be aligned what do I do?
Matt J
Matt J el 19 de Sept. de 2013
Editada: Matt J el 19 de Sept. de 2013
Sample full periods of the sinusoid, so that there is no discontinuity between x(1) and x(end).

Iniciar sesión para comentar.

Categorías

Más información sobre Fourier Analysis and Filtering en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by