数字信号处理第三四章 DFT及FFT.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 时间偶奇分 频率前后分 计算量：一个复数乘法，两个复数加法 1.本级有4个（N/2）蝶形 2.N点DFT可由两个N/2DFT计算 8点FFT示意图1 N/2点DFT可由两个N/4DFT计算 本级共有有4个（N/2）蝶形 8点FFT示意图2（基本原则同示意图1） 时间偶奇分 频率前后分 8点DFT的FFT示意图 共有3（M）级，每级有4（N/2）个蝶形，共M* N/2个蝶形 两点DFT，蝶形 不能再分解了 基本运算：蝶形运算环节 FFT的运算量为： M* N/2个复数乘法， M* N个复数加法 3 1）同址（原位）运算 每一级的蝶形的输入与输出在运算前后可以存储在同一地址的存储单元中。 同址运算的优点：节省存储单元，降低了硬件实现成本 2）蝶形运算两节点间的“距离” 第一级：1 第二级 2 第三级：4 蝶形两节点的“距离”： 其中：L:第几级蝶形运算 3）旋转因子的变化规律 第一级 第二级 第三级 4)蝶形运算规律 第L级 5)变址运算 混 序 正序输入 正序输出 变址运算 变址运算的目的：实现FFT算法的正序输入、正序输出 变址计算：IJ，x(I)与x(J)交换位置 6）编程思想 三、频率抽选基2-FFT算法Decimation-In-Frequency FFT---DIF FFT 对于任何流程图，只要保持各节点所连支路及其传输系数不变，则不论节点怎样排列，所得到的流程图总是等效的。 对于任何流程图，只要保持各节点所连支路及其传输系数不变，则不论节点怎样排列，所得到的流程图总是等效的。 四、快速逆傅里叶变换---IFFT 对上式两端取共轭 X*(k)的DFT 五、实序列的FFT 进一步减小计算量 1.构建两个新序列 2.用FFT计算出DFT[y(n)] 3.由共轭对称性 4.由蝶形运算和实序列性质 3.6 快速傅里叶变换的应用The Application of FFT 一、用FFT对连续信号进行谱分析 3.3 频率采样Frequency Sampling 一、序列的几种变换之间的关系 二、由N个频率取样值X(k)恢复序列x(n) 三、从N个取样值恢复X(Z)或X(ejw) 3.4 利用循环卷积计算线性卷积 L=5点循环卷积 L=6点循环卷积 L=7点循环卷积 L=8点循环卷积 分段卷积 3.5 快速傅里叶变换Fast Fourier Transform——FFT 一、DFT的计算量 直接利用DFT的定义式计算DFT，如下式所示 二、时间抽选的基2-FFT( )Decimation-In-Time—DIT FFT 二、离散傅里叶变换的性质 B、序列分解 C、DFT的共轭对称性 有限长序列 有限长共轭对称序列 有限长共轭反对称序列 Re[DFT[x(n)]] Im[DFT[x(n)]] Re[DFT[xep(n)]] Im[DFT[xep(n)]] Re[DFT[xop(n)]] Im[DFT[xop(n)]] Re[DFT[x(n)]] Im[DFT[x(n)]] Re[DFT[xep(n)]] Im[DFT[xep(n)]] Re[DFT[xop(n)]] Im[DFT[xop(n)]] Re[DFT[x(n)]] Im[DFT[x(n)]] Re[Xep] Im[Xop] M=16点循环卷积 M=8点循环卷积 X1(n)，N1=4 X2(n)，N2=4 循环卷积的圆形示意图 第三章 DFT 主要内容： 引言 频谱连续,非周期 频谱离散,非周期 ? 非周期离散时间（或有限长序列x[n]）的Fourier 变换 若 0≤n ≤N－1 ?离散Fourier变换 (DFT) ? 周期离散时间序列的Fourier 变换—DFS 频谱连续,周期 频谱离散, 周期 频谱离散 一、离散傅里叶级数（DFS） 3.1 周期序列的ZT不收敛 可用离散傅里叶级数表示离散周期序列，即用周期为N的复制数序列来表示 令WN = e－j 2π/N旋转因子 可见 -2N/M -N/M 0 N/M 2N/M e.g. M=4,N=12 …… …… X(k) 二、离散傅里叶级数的性质 3.1 离散傅里叶变换及其性质 一、离散傅里叶变换定义 8、DFT与ZT、序列Fourier的关系 x(n)的离散傅里叶变换的结果与变换区间N