参考连接:https://ww2.mathworks.cn/help/matlab/math/basic-spectral-analysis.html
注:生成的波形周期个数越多,得到的频谱越尖。增大tt可以使频谱更纯粹。
wx = 2pi f * x
%% matlab 采样,FFT,得到频率
FFSS = 4000;%采样频率
FF = 200; %200Hz
TT = 1/200;%0.005 s
QQ = 0; %相位
tt = 0: 1/FFSS: TT*20-1/FFSS;
Ssin = sin(2*pi*FF*tt+QQ);
figure
subplot(2,1,1);
plot(tt,Ssin);
fft_Ssin=fft(Ssin);
n =length(fft_Ssin);
f = (0:n-1)*(FFSS/n); % frequency range
power = abs(fft_Ssin).^2/n;
subplot(2,1,2);
plot(f,power);