第6章 快速傅立叶变换.ppt

第6章 快速傅立叶变换 电气学院 自动化系 6.1 直接计算DFT的问题及改进途径 累加式中一项的运算量：4次实数乘法，两次实数加法 计算X(k)的一个数值需要4N次实数乘法和（4N-2）次实数加法 完成X(k)的所有数值计算需要：4N2次实数乘法和N（4N-2）次实数加法; 或N2次复数乘，N(N-1)次复数加 6.2 按时间抽取（DIT）的基2FFT算法 基本思想： 计算时间～N2 缩短序列长度可以大大节省计算时间 将一个长序列分成若干个短序列，分别计算DFT，再对它们进行组合得到原序列的DFT 序列分解——将x(n)按n的奇偶性分成两序列分别计算DFT 前提：序列的长度N=2L 组合两个子序列的DFT得到原序列的DFT 前N/2个值 后N/2个值 FFT的 蝶形运算图： 前N/2个值的组合公式 X1(k)和 X2(k)的周期均为N/2 后N/2个值的组合公式 FFT算法中迭代运算的思想 将子序列继续分解计算，直至所有的子序列均为2点序列 FFT算法只对2点序列计算DFT：(也是蝶形运算） -1 x(0)+WN0x(4) x(0)-WN0x(4) 例：利用FFT算法计算16点序列x(n),n=0,…,15的DFT，请在算法左侧填入序列x(n)的输入顺序 8点序列DFT 8点序列DFT 16点序列DFT 4点序列DFT 4点序列DFT 4点序列DFT 4点序列DFT 2点序列DFT 2点序列DFT 2点序列DFT 2点序列DFT 2点序列DFT 2点序列DFT 2点序列DFT 2点序列DFT X(1) X(2) . . . . . . . . . X(15) FFT算法中迭代运算的思想 FFT的运算时间分析 N=2L时，有L级蝶形运算，每级蝶形运算的个数都为N/2个，每个蝶形运算有一次复数乘和二次复数加 例：对于一个32点序列，若使用FFT算法计算其离散傅立叶变换，则一共需要________次蝶形运算。 DFT与FFT的运算时间比 【例6-1】若使用某DSP处理芯片来计算2048点DFT，已知一次复数乘法运算需要400ns，一次复数加法运算需要100ns，则直接计算DFT需要多少时间？采用基-2 DIT FFT需要多少时间？ DFT需要2.0969s FFT需要0.0068s IFFT 6.7 DFT的应用 1.利用FFT对连续信号进行频谱分析 对连续信号x(t)进行抽样和截取 “截取”——x(n)RN(n)相当于频域信号的卷积 “截取”过程会带来“频谱泄漏”问题，导致对原连续信号的频谱特性无法做准确的分析 (π) (π) Ω/π 频谱泄漏 Ω/π n w(n) n N=16 xN(n) Ω/π n w(n) Ω/π n xN(n) N=32 解决频谱泄漏的方法：加宽窗函数 6.7 DFT的应用 2 离散傅立叶变换参数选择 频率分辨力——两个离散谱峰的频率间隔（单位：Hz），表示DFT对两个最近频谱峰的分辨能力 对连续信号的时域记录点数 提高频率分辨力的方法：增加记录点数或降低抽样频率 两个离散点之间的频谱分布不明——“栅栏效应” 例：已知对某连续信号f(t)的截止频率为500Hz，若对该信号通过FFT分析其频谱特性，要求频率分辨力为0.1Hz，则对f(t)的抽样点数N至少为____________。 6.7 DFT的应用 3.信号消噪 信噪比：用于衡量信道的通信质量，单位（dB） 噪声的可消除条件：噪声频率与信号频率范围不同