直流分量: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
本文发布于:2024-01-31 19:56:33,感谢您对本站的认可!
本文链接:https://www.4u4v.net/it/170670219330984.html
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。
留言与评论(共有 0 条评论) |