Computing and Ploting Fourier transform

3 次查看（过去 30 天）
omarahmed.m
omarahmed.m 2016-2-21
评论： omarahmed.m 2016-2-21
How can i use the fft command to transform a rectangular pulse to sinc function and plot the sinc function, i'm using a very traditional way to compute the fourier transform and plot it, but this way is too slow, when i use the fft command and try to plot, the ploting of the magnitude is only the pins or a spike, i need the plot to be a sinc function as the picture that i attached here, also i'm using the same method to compute the inverse fourier transform so i also need to use the ifft command, can anyone help?
if true
p = heaviside(t+1)-heaviside(t-1);
k = 0;
F=zeros(1,8001);
for f= -4:0.001:4
k=k+1;
F(k) = trapz(t,p.*exp(-1i*2*pi*f*t));
end
f = -4:0.001:4;
figure(1);
plot(f,F);
grid
end

请先登录，再回答此问题。

采纳的回答

Rick Rosson
Rick Rosson 2016-2-21
编辑：Rick Rosson 2016-2-21
Fs = 8;
dt = 1/Fs;
N = 8000;
t = dt*(-N/2:N/2-1)';
dF = Fs/N;
f = -Fs/2:dF:Fs/2-dF;
p = heaviside(t+1)-heaviside(t-1);
P = dt*fftshift(fft(ifftshift(p)));
figure;
plot(f,real(P));
grid on;
  1 个评论
omarahmed.m
omarahmed.m 2016-2-21
Thank you for your response, Could you just specify this command "P = dt*fftshift(fft(ifftshift(p)));", why it consist of fftshift of fft of ifftshift, and why multiplying by "dt". thank you very much.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by