Matlab问题。这个fft代码在做什么？

Question

我是 Matlab 的新手，正在执行信号处理。我想了解这段代码在做什么？我们如何以及为什么要确定指数奈奎斯特和频谱？

spectrum = fft(Signal,k);              %generate spetrum of signal with FFT to k points
indexNyquist = round(k/2+1);        %vicinity of nyquist frequency        
spectrum = spectrum(1:indexNyquist);       %truncate spectrum to Nyquist frequency
spectrum = spectrum/(length(Signal));      %scale spectrum by number of points  
spectrum(2:end) = 2 * spectrum3(2:end);  %compensate for truncating negative frequencies, but not DC component

Answer 1

对于纯实输入信号，对应的 FFT 将是关于奈奎斯特频率的复共轭对称的，因此顶部 N/2 个 bin 中没有有用的附加信息。因此，我们可以只取底部的 N/2 个 bin 并将它们的幅度乘以 2 以获得没有冗余的（复杂）频谱。此频谱表示从 0 到 Nyquist 的频率（及其混叠等效频率）。

请注意，bin 0（0 Hz aka DC）是纯真实的，不需要加倍，因此您的 Matlab 代码中的注释。