设信号的数据长度为L,其N点DFT定义为信号的DTFT在奈奎斯特区间 上N个等间距的频率点处的频谱密度值,即信号的N点DFT结果是N个频谱密度值。 在奈奎斯特区间上均匀分布的N个DFT频率为
根据DFT的定义,将这些频率值代入DTFT的频谱公式中,得
在算法上,DFT计算式的一种最简单的实现是:套用前一节所讲的任意区间上的DTFT计算方法,将区间的起点和终点分别给定为0和 即可,如下面的程序实现。
#include <comlx.h> void dtftr(); void dft (int L, double * x, int N, complex * X) { double pi=4*atan(1.0); for(k=0;k<N;k++) dtftr(L, x , N , X , 0.0 , 2*pi); }