Question

求Matlab大神，问下fft做完以后如何写代码求得对应某个频率的幅值？

这是小弟的代码，谁能帮我改一下，输出52k频率信号的幅值，谢谢了。function dc = myfft(signal, fs, len, caption)if len == 0N = length(signal);elseN = len;endN;ts = 1/fs;data = signal(1:N);S = fft(data,N);CS = fftshift(S);freq = ([-N/2: N/2-1] + 1/2)/(N*ts);mag = abs(CS) /(N/2);magS = abs(S) / (N/2);dc = magS(2);figure;plot(freq/1e3,mag);set(gca, 'LineWidth', 1.5 ); set(gca, 'TickLength',[0.005 0.0025]);set(gca, 'FontSize', 10);set(gca, 'FontWeight', 'bold');box ontitle(caption, 'FontSize', 14, 'FontWeight', 'bold');xlabel('Freq. /kHz', 'FontSize', 12);ylabel('Mag. /fd ', 'FontSize',12);end

Answer 1

你fft得到频谱时
得到了频率横坐标 freq
幅度纵坐标 mag

如果你要知道某个频率 f0 的幅值
假如你的 f0 是你得到freq 坐标中的其中一个值
mag0=mag(freq==f0);
mag0就是对应f0的幅度值

如果你的f0 不在freq里 （例如 freq=[ 1 2 3 4 5 ...],你的f0是2.5），那么就需要插值
mag0=interp1(freq,mag,f0);