fft负频率部分.docVIP

FFT运算中的负频率成分 1. fft 函数的负频率 在使用FFT运算进行离散时域信号的频谱分析时，我们知道这样一个事实，那就是：如果时域信号是一组实值序列x，那么由FFT运算得到的频域矢量X则具有共轭对称的属性，而相应的幅度谱矢量S abs X 则具有偶对称的属性。因此，我们在考察时域信号的频谱时，只需要取X或者S的前半部分就可以了。 以一组余弦信号为例： x t A/2*[cos 2πf1t + cos 2πf2t ] 其中：f1 1 Hz，f2 58 Hz，信号振幅A 58 Hz，离散时间t为0到1秒的共128个采样点。因此，采样间隔为Ts 1/ 128-1 0.007874 s，采样频率fs 127 Hz。 使用Matlab的fft 函数进行频谱计算： X fft x ; f 0:127; figure stem f,abs X ; 最终得到以下的频谱： 对于常规的频谱分析，目的已经达到了，因为我们可以看到在频谱图的左半，f1 1 Hz 和 f2 58 Hz 两个频点很清晰地显示了出来，而零频率点则是整条频谱的第一个点，以绿色标出。 确实，我们也看到了，除去零频率点，整个频谱图的确是偶对称的，那么频谱图的右对半又有什么物理意义呢？ 其实，这右对半就是频谱的负频率成分。 2. 调整fft 的频谱图 在Matlab中，为了把这部分负频率成分还原到频谱图的左半，以符合常规的数学坐标表达方式，我们需要使用fftshift 函数来进行还原。可参考以下代码： X1 fftshift X ; f -64:63; figure stem f,abs X1 ; 最终得到以下的频谱： 可以看到，绿色的零频点已经从频谱图的最左端搬移到了正中。这样就得到了常规的频谱坐标图象。 备注： 如果要使用ifft 函数将经过调整的频谱矢量X1还原为原始的离散时间信号x t ，则需要使用ifftshift 函数将X1还原为左正右负的方 式（其实就是未经过调整的频谱矢量X），有一个等式精辟概括了fft，ifft，fftshift和ifftshift这四个函数的关系： ifft ifftshift fftshift fft anything anything 3. 确定频谱图的横坐标 在第二节里，我们看到了，对于128个值点的频谱矢量X，它的奈奎斯特采样频率是fs 127 Hz，于是混叠频率fd 63.5 Hz，按照对称的关系，横坐标f矢量的取值应该是 -63.5，-62.5 …… -1.5，-0.5，+1.5 …… 62.5，63.5 。 可是，实际的f矢量的取值却是 -64，-63 …… -1，0，+1 …… 63 ，为什么会这样呢？ 这是因为，在常规情况下，fft 得到的频谱矢量，第一个点都总是零频率点，即便是经过fftshift 调整后，零频率点也是存在的。所以，采用 -63.5，-62.5 …… -1.5，-0.5，+1.5 …… 62.5，63.5 这种横坐标范围是错的。 这种错误在频谱矢量值点数为偶数个（如N 128）时容易出现，因为此时的频域间隔是奇数个 N-1 127 。其实这是的频谱图并不是绝对对称的，因为代表负频率的值点有N/2个，而代表正频率的值点只有N/2 – 1个，剩下的一个零频率点。 当然，如果频谱矢量值点数为奇数个（如N 129）时，此时的频谱图是绝对对称的，就不存在上述确定横坐标的困惑。 不过对奇数个离散序列使用fft ，实际上做的还是常规的DFT运算，这样运行效率得不到保证，所以还是应该尽量保证值点个数为2的幂。