how to extend time after fft ?

15 vues (au cours des 30 derniers jours)
omar thamer
omar thamer le 8 Nov 2013
Commenté : omar thamer le 13 Nov 2013
I have a signal runs for 500 seconds ( 1:1:500 sec) after i did the fft, i want to go back to time domain and see how the signal looks like for 1000 Seconds (1:1:1000) ? Any advice how to extend the time ? Is it better to do it in frequency domain ? Is it simply to add the sines and cosines with loops ?
Thank you

Réponse acceptée

Vilnis Liepins
Vilnis Liepins le 11 Nov 2013
ifft(edft(YourData,1500),1500)
Signal for 1000 sec will be ifft output 1:1000.
  2 commentaires
omar thamer
omar thamer le 12 Nov 2013
Thank you Vilnis, really helped a lot.Though, results not as expected for real signal but absolutely true for periodic computer generated one. I am wondering if you can advice on my signal which i attached here. As you can see a lot o discontinuities and no clear center frequency.
Any advice how to get periodicity to my signal ? I tried with STFT it also yields no results.
omar thamer
omar thamer le 13 Nov 2013
Its not working properly .. Have a look to my code:
t=1:1:400;
x=sin(2*pi*t*1/20);
X=edft(x,500);
Y=ifft(X,500);
the outcome for extended signal has different magnitude and frequency.

Connectez-vous pour commenter.

Plus de réponses (1)

Walter Roberson
Walter Roberson le 8 Nov 2013
Try
ifft(fft(YourData), 1000)
  2 commentaires
omar thamer
omar thamer le 8 Nov 2013
it will be padded with zeros then and return a complex vector. Its not working
Walter Roberson
Walter Roberson le 11 Nov 2013
fft() assumes that its data is periodic, so to get 1000 seconds of output, ifft() to get the 500 second output and then replicate it.

Connectez-vous pour commenter.

Catégories

En savoir plus sur Fourier Analysis and Filtering dans Help Center et File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by