数字信号处理第4章 快速傅立叶变换(FFT).pptVIP

1.2.2. 序列的逆序排列 由于 x(n) 被反复地按奇、偶分组，所以流图输 入端的的排列不再是顺序的，但仍有规律可循： 对于任意 n 可以用M个二 进制码表示为： 当 n 反复按奇、偶分解时，即按二进制码的“0”， “1” 分解。 n0 0 1 n1 0 1 0 1 n2 0 1 0 1 0 1 0 1 … 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 0 1 逆序 （码位倒置） 比较 顺序 BIN DEC 000 100 010 110 001 101 011 111 0 4 2 6 1 5 3 7 BIN DEC 0 1 2 3 4 5 6 7 000 001 010 011 100 101 110 111 运算前 运算后 1.2.3. 同址运算 对于算法流图中的任意一个蝶形运算。 1.2.4. FFT的运算量 ∴可进行M次分解，即N 点FFT可经过M 次 迭代运算得到，每次迭代中的蝶形运算为 个, 每个蝶形 运算中有一次复乘，∴复乘次数 例如 FFT：复乘次数 DFT：复乘次数 相差 1260倍！ 1.3 基２ＦＦＴ按频率抽取算法 （sande-Tukey算法） 频域抽取基－2FFT算法（DIF-FFT） 算法的推导 频域抽取算法是把时间序列 前后对半分解为两个长为N/2点的序列，则： 当 k 取偶数时（ k = 2 r，r = 0 , 1 , ... , N / 2－ 1） ∴ 的N点DFT 按 k 的奇偶分组可分为两个N/2的DFT k 取奇数时（ k = 2 r＋1，r = 0 , 1 , ... , N / 2－ 1） 这一结论表明：求 的N点DFT 再次分解成 求两个N/2 点DFT DIF-FFT的蝶式运算流图 DIF-FFT的一次分解运算流图 先蝶式运算，后 DFT。例如：N=8时 DIF -FFT 1.3.1. 算法的推导 将x(n)按前后分为两组： 1.3.2. 流图的转置 频域中的第 l 点与时域中的第 r 点的关系为 频域中的第 r 点与时域中的第 l 点的关系为 从流图中也可看到： 流图转置： X(?)改写为x(?)，序号不变； x(?)改写为X(?)，序号不变； 所有箭头反向，系数不变。 DIT -FFT DIF -FFT 时间抽取 频率抽取 1.3.3. 两种蝶形运算流图的区别 1.4 FFT 算法的软件实现 (下页) （框图） 1.4.1. 算法的特点 1．基本算法为蝶形运算； 2．运算分为 次迭代进行； 3．运算可原址进行； 4．输入数据需进行逆序重排； 5 ．迭代过程中，基本蝶形运算的两个点之间 的间距及系数 的变化有规律可循。 在第L 轮迭代 中 （1）有 种蝶形运算系数 （2）系数分别为 （3）每种系数对应的蝶形运算有 个，相距 点 （4）参加同一蝶形运算的两个点之间的间距为 点 （返回） 1.4.2. 算法的程序框图 开始 逆序重排 蝶形运算 结束 L=1,M 输出 J=0,B-1 蝶形运算 三层循环： N1为相同系数的蝶形运算的间距（ N1= 2 L ） 2）中间层为蝶形运算组次J（J=1，2，…， 3）最内层为基本蝶形运算 1）最外层为迭代轮次L（L=1，2，…，M） 2 E- E1 -1） N2为参加同一蝶形运算的两点间的间距 （ N2= 2 L-1 = N1 /2 ） J为系数的序号（ J=0，1，…， N2 ） D为每组蝶形运算的第一个数组元素的序号 1.5 IFFT 算法 1.5.2. IFFT的一种算法 1.5.1. IFFT另一种算法流图 对照 可仿照FFT算法，且 ，而 的运算 即为二进制数的小数点左移一位，∴ 频率抽取IFFT 时间抽取IFFT * 第4章 快速傅立叶变换（FFT） 课程名称：数字信号处理 引言 四种不同傅立叶变换对 傅 里 叶 级 数(FS)：连 续 时 间 , 离 散 频 率 的 傅 里 叶 变 换 。 连 续 傅 里 叶 变 换(FT)：连 续 时 间 , 连 续 频 率 的 傅 里 叶 变 换 。 序 列 的 傅 里 叶 变 换(DTFT):离 散 时 间 , 连 续 频 率 的 傅 里 叶 变 换. 离 散 傅 里 叶 变 换(DFT):离 散 时 间 , 离 散 频 率 的 傅 里 叶 变 换 傅立叶变换 离散时间傅立叶变换 傅立叶级数 离散傅立叶变换 时间 ? 连续 离散