필터 지우기
필터 지우기

Problem with frequency axis during Fourier Transform of Gaussian Pulse

조회 수: 3 (최근 30일)
Vivian
Vivian 2014년 2월 28일
댓글: Rick Rosson 2014년 3월 4일
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
Rick Rosson
Rick Rosson 2014년 3월 3일
편집: Rick Rosson 2014년 3월 3일
Please post screen shots of the spectrum for each of two different values of dt.

댓글을 달려면 로그인하십시오.

답변 (1개)

Rick Rosson
Rick Rosson 2014년 3월 3일
편집: Rick Rosson 2014년 3월 3일
A few minor tweaks:
t = 0:dt:2e-9-dt;
df = fs/nfft;
f = -fs/2:df:fs/2-df;
  댓글 수: 2
Vivian
Vivian 2014년 3월 3일
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 2014년 3월 4일
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.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Discrete Fourier and Cosine Transforms에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by