第四章离散傅里叶变换和其快速算法.ppt

用 DFT 对连续时间信号进行频谱分析必然是近似的，其 近似程度与信号带宽，抽样频率和窗函数宽度等有关，近似 公式为： 2 1 0 6.27 6.28 ( ) ( ) [ ( )] 1 ( ) [ ( )] N j kn N N N n N N X k T x n e T DFT x n x n IDFT X k T ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? g g （ ） （ ） 频谱分析存在的问题 混叠现象 栅栏效应 频谱泄漏 混叠现象 时域离散 频域周期，混叠 要求： fs ≥2fmax 频域抽样： F=fs/N=1/(NTs)=1/tp 或 N=fs/F F 与信号的实际长度有关。 指在频域中能辨认的频率，即频域取样中两相邻点间的 频率间隔或基波频率。 F 愈小，频率分辨率愈高。 截断效应 x(n) 可能无限长 , 将 x(n) 截短，相当于 y(n)=x(n)R N (n). 则： ( ) 1 ( ) [ ( )] ( ) ( ) 2 1 ( ) ( ) 2 j j j N j j N Y e FT y n X e R e X e R e d ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 其中： 1 ( ) 2 sin( / 2) ( ) [ ( )] ( ) sin( / 2) N j j j N N N N R e FT R n e R e ? ? ? ? ? ? ? ? ? ? ? ? 例如， x(n)=cos(ω 0 n), ω 0 =π/4 其频谱为 ( ) [ ( 2 ) ( 2 )] 4 4 j l X e l l ? ? ? ? ? ? ? ? ? ? ? ??? ? ? ? ? ? ? ? ( ) 1 ( ) [ ( )] ( ) ( ) 2 1 ( ) ( ) 2 j j j N j j N Y e FT y n X e R e X e R e d ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 其频谱图如图所示： 所以，截断后序列的频谱与原序列的频谱有差别。 1 、泄漏现象：截断后，使原来的离散谱线向附近展宽， 通常称这种展宽为泄漏，它使频谱模糊，谱分辨率下降。 2 、谱间干扰：在主谱线两边形成很多旁瓣，引起不同频 率间的干扰，影响频谱分辨率。 上述两种现象又称为截断效应。 截断效应 例： 试利用 DFT 分析一连续信号，已知其最高频率 =1000Hz ，要求频率 分辨率 F ? 2Hz ， DFT 的点数必须为 2 的整数幂次，确定以下参数：最大的 抽样间隔，最少的信号记录时间，最少的 DFT 点数。 解： (1) 最大的抽样间隔 T max 为 ) ( 10 5 . 0 1000 2 1 2 1 3 max s f T m ? ? ? ? ? ? (2) 最少的信号持续时间 T pmin 为 ) ( 5 . 0 2 1 1 pmin s f T c ? ? ? ? (3) 由最大的抽样间隔 T max 与最少信号持续时间 T pmin ，可得最少 DFT 点数 N 为 1000 10 5 . 0 5 . 0 3 max min ? ? ? ? ? T T N p 选择 DFT 的 点数 N =1024 ，以满足其为 2 的整数幂次。 F 例：已知语音信号 x ( t ) 的最高频率为 f m =3.4kHz, 用 f sam =8kHz 对 x ( t ) 进 行抽样。如对抽样信号做 N =1600 点的 DFT ，试确定 X [ k ] 中 k =600 和 k =1200 点所分别对应原连续信号的连续频谱点 f 1 和 f 2 (kHz) 。 解：对连续信号 x ( t ) 按 f sam =8kHz 进行取样，得到对应的离散序列 x [ n ] ，在利用 离散序列 x [ n ] 的 DFT X [ k ] 分析连续信号 x ( t ) 的频谱时， X [ k ] 与 X (j ? ) 存在以下对 应关系： 当 k =600 时，由于 0 ? k ? N /2 ? 1 ，所以 kHz 3 600 1600 8 sam 1 ? ? ? ? m N f f 当 k =1200 时，由于 N /2 ? k ? N ，所以 kHz 2 ) 1600 1200 ( 1600 8 ) ( sam 2 ? ? ? ? ? ? ? N m N f f A ? ? m ?? m ? X (j ? ? A/T m X (e j ? ? ? ? 2 ? / N X [ m ] ?? ??? k k N 点 DFT 是在频率区间 [0 ，2π]上对信号频谱进行 N 点 等间隔采样，得到的是若干个离散的频谱点 X （ k ）