如何快速的将时域离散信号的频谱通过MATLAB的fft2函数表示出来:
clc
clear allclose all
fs = 500;
t = 0:1/fs:1.5;
f1 = 40;
f2 = 20;
x = sin(2 * pi * t * f1) + sin(2 * pi * t * f2); %input signalfigure(1)
plot(t, x);
title('input signal x[n]'); %input signal wave
figure(2)
f = (0:750) * fs / 751 - fs / 2; %750 derives from t, which (1.5 - 0) / (1 / fs) = 750, 0:1:750 there are 751 values
plot(f, abs(fftshift(fft(x)))); % plot input signal's frequency spectrum
%利用fftshift函数是必要的,fftshift函数的作用正是让正负半轴的图像分别关于各自的中心对称,