Effacer les filtres
Effacer les filtres

Problem with frequency axis during Fourier Transform of Gaussian Pulse

3 vues (au cours des 30 derniers jours)
Vivian
Vivian le 28 Fév 2014
Commenté : Rick Rosson le 4 Mar 2014
Hello, I noticed that whenever I change dt(sampling time) the bandwidth of FFT Gaussian changes. Smaller dt will have wider bandwidth.
  • Why is this happening? Doesn't smaller dt just allows more points for better signal sampling in time domain?
clc
clear all
close all
dt=15e-12;
t=0:dt:2e-9;
tao = .5e-9;
nfft = length(t);
%frequency spectrum
fs = 1/dt;
df = fs/(nfft-1);
f=(-fs/2):df:(fs/2);
sigma1=.1E-9;
x1=((t-tao)/sigma1).*((t-tao)/sigma1);
G1=(1/(sqrt(6.28)*sigma1)).*exp(.5.*(-x1));
%Fourier Transform
G7 = fftshift(abs(fft(G1)));
plot(f/1E9,G7)
xlim([-50 50]);
Thanks,
  1 commentaire
Rick Rosson
Rick Rosson le 3 Mar 2014
Modifié(e) : Rick Rosson le 3 Mar 2014
Please post screen shots of the spectrum for each of two different values of dt.

Connectez-vous pour commenter.

Réponses (1)

Rick Rosson
Rick Rosson le 3 Mar 2014
Modifié(e) : Rick Rosson le 3 Mar 2014
A few minor tweaks:
t = 0:dt:2e-9-dt;
df = fs/nfft;
f = -fs/2:df:fs/2-df;
  2 commentaires
Vivian
Vivian le 3 Mar 2014
Thank you for looking into this. Actually, I wasn't quite clear on what I was asking. I was wondering about why if changed:
t=0:dt:2e-9;
to
t=0:dt*4:2e-9;
The bandwidth of G7 widened. I think it is incorrect to change dt to dt*4, as this changes the time spacing between the number of samples in time domain without reflecting the changes of it in the frequency domain. The proper way is to:
dt=4*15e-12;
To show the correct fourier transform of the sampling frequency. Please correct me if I'm wrong.
Thank you!
Rick Rosson
Rick Rosson le 4 Mar 2014
Yes, that is correct. If you change dt, then that also changes fs, and that will change the range of frequencies over which the Fourier spectrum is defined. So it is crucial that when you change dt, you recompute fs and everything else that depends on dt.

Connectez-vous pour commenter.

Community Treasure Hunt

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

Start Hunting!

Translated by