Question

怎样实现 sinc函数的频谱图

Answer 1

t=-15:0.001:15;
fa=sinc(t/pi); %Sa函数的原函数
subplot(211);
plot(t,fa); %在（-15,15）区间内构造Sa函数图像
xlabel('t');
ylabel('Sa（t）');
title('在区间（-15，15）内构建Sa(t)函数');
grid on;
%下面是想要实现对sinc函数的傅里叶变换并且构图。。但是不知道fft函数要怎么用。
subplot(212);
Fa=log(1+abs(fftshift(fft(Sa)))); %对Sa函数进行傅里叶变换
plot(Fa); %在（-3,3）内构造频谱图
xlabel('w');
ylabel('Fa（t）');
title('Sa(t)的频谱函数');
grid on;