这是代码:
关闭所有;
%定义样本的数量
fs = 100;
f = 400;%赫兹
% N =长度(信号);
%定义信号
t = 0:1 / fs: 1 - 1 / f;
% =罪信号(2 *π* f * t);
信号= xlsread (“testdata.xlsx”);
%的情节来说明它是一个正弦波
情节(t,信号)
标题(时域信号的);
%进行傅里叶变换
fftSignal = fft(信号);
% fftshift适用于把它放在我们习惯于形式(见文档)
fftSignal2 = fftshift (fftSignal);
% xdft = xdft(1:长度(s) / 2 + 1);
%,计算频率轴,定义的采样率
f = f / 2 * linspace (1, 1, fs);
%因为信号复杂,我们需要情节让它的大小
%看起来正确,所以我们用abs(绝对值)
图;
情节(f, abs (fftSignal2));
标题(“正弦的FFT级”);
包含(的频率(赫兹));
ylabel (“级”);