%Data 为一维采样数组
% Fs 为采样频率
Data_length=length(Data);
% DFT需要的采样点数为2的幂指数,但是输入的点数有可能是一个一般整数,于是找一个离指定点数最近的2的幂指数用来做DFT
NFFT = 2^nextpow2(Data_length);
% 生成结果的x轴,即频域。正常出来的频谱是左右对称的,为便于处理只要前半个域。
F_domain = (Fs/2)*linspace(0,1,NFFT/2);
% FFT变换
Y = fft(Data,NFFT)/Data_length;
% 幅值
Amp =2*abs(Y(1:NFFT/2));
subplot(211);
plot(F_domain,Amp);
title('频谱')
xlabel('频率(Hz)')
ylabel('频谱幅值')
傅里叶变换的Matlab代码与注释
最新推荐文章于 2024-03-01 18:14:26 发布