文章目录
1、 fft中的幅度
直流分量:abs(X[0])/N
其他频点的分量: abs(X[f])/(N/2)
幅值转db时注意pow2db()和db()的用法,前者是功率,后者是幅度。
下面为具体的fft程序,程序中更通用的是对于非直流成分的fft的幅度值
function fft_plot(y,fs,s_name)
% L_i = length(y)*100;
L_i = 2^nextpow2(length(y)*100);
s_i_fft = fft(y,L_i);
s_i_fftshfit = fftshift(s_i_fft);
P = abs(s_i_fftshfit/length(y))*2; %对于非直流分量或者0来说是除以N乘以2,对于直流分量来说是直接除以N
% fshift = (-L_i/2:L_i/2-1)*(fs/L_i);
fshift = linspace(-fs/2,fs/2,L_i);
figure;
plot(fshift,P);
title([s_name,'的双边谱 ']);
xlabel('f (Hz)');
ylabel('|P(f)|');
end