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

重叠相加法例子 (2)重叠保留法 ? 这种方法和第一种方法稍有不同，即将上面分段序列中补零的部分不是补零，而是保留原来的输入序列值，这时，如利用DFT实现h(n)和xi(n)的循环卷积，则每段卷积结果中有N1-1个点不等于线性卷积值需舍去。 重叠保留法与重叠相加法的计算量差不多，但省去了重叠相加法最后的相加运算。 A、先将x(n)分成几个短序列xi(n),每段长度为N点，M为短序列的长度，这样分段后，应在每一段的前边再补上前一段保留下来的M-1个输入序列值，组成L点短序列：L=N+M-1，进行圆周卷积，而对第一段，由于没有前一段保留信号，则需要在它前面填充M-1个零值点 B、计算N点FFT，H(K)=DFT[h(n)], Xi(K)=DFT[xi(n)] C、相乘 Yi(K)=H(K)Xi(K) D、计算N点IFFT yi(n)=IDFT[Yi(K)] E、由于每段圆周卷积结果的前M-1个点不等于线性卷 积的值，必须舍去，再把各相邻输出段留下来的序 列衔接起来，就构成了最后的正确输出。 （2）重叠保留法例子 * * * （3）序数重排 它的输入正好是自然顺序。但它的输出却是码位倒置的顺序。因此运算完毕后，要通过变址运算将码位倒置的顺序转换为自然顺序，然后输出，变址方法同时间抽取法。 （4）蝶形类型随迭代次数成倍减少(与时间抽取相反) 第一级迭代中有N/2种蝶形运算系数，参加蝶形运算的两个数据相隔N/2，随后每次迭代，蝶形类形比前一级减少一倍，间距也减少一倍，最后一级迭代，蝶形类形只有一种W0N，数据间隔为1。 由这几点规律可以看出，频率抽取法与时间抽取法是两种等价的FFT运算。 1.相同点 (1)进行原位运算； (2)运算量相同,均为（N/2) Log2N次复乘,N Log2N次复加。 2.不同点 (1)DIT输入为倒位序,输出为自然顺序； DIF正好与此相反。但DIT也有输入为自然顺序,输出为倒位序的情况。 (二).DIF法与DIT法的异同 (2)蝶形运算不同 a.(DIT) 用矩阵表示 = 1 1 b.(DIF) 用矩阵表示 = 1 1 (3)两种蝶形运算的关系 互为转置（矩阵）； 将流图的所有支路方向都反向,交换输入和输出，即可得到另一种蝶形。 a.(DIT) b.(DFT) 1 1 1 1 DIF与DIT根本区别：在于蝶形结不同。 DIT的复数相乘出现在减法之前。 DIF的复数相乘出现在减法之后。 §2-6 IFFT算法 一.稍微变动FFT程序和参数可实现IFFT 比较两式可知,只要DFT的每个系数 换 成 ,最后再乘以常数1/N就可以得到IDFT 的快速算法--IFFT。 另外,可以将常数1/N分配到每级运算中, ∵ ,也就是每级蝶形运算均乘 以1/2。 利用FFT计算IFFT时在命名上应注意： (1)把FFT的时间抽取法，用于IDFT运算时，由于输入变量由时间序列x(n)改成频率序列X(k),原来按x(n)的奇、偶次序分组的时间抽取法FFT，现在就变成了按X(k)的奇偶次序抽取了。 （2）同样，频率抽取的FFT运算用于IDFT运算时，也应改变为时间抽取的IFFT。 二、改变FFT流图系数的方法 1、在IFFT的运算中，常常把1/N分解为(1/2)m，并且在M级运算中每一级运算都分别乘以1/2因子，就可得到IFFT的两种基本蝶形运算结构。 A B A B (b)时间抽取IFFT的蝶形运算 (a)频率抽取IFFT的蝶形运算 2、IFFT的基本蝶形运算 三.不改(FFT)的程序直接实现IFFT 这就是说,先将X(k)取共轭,即将X(k)的 虚部乘-1，直接利用FFT程序计算DFT；然后 再取一次共轭；最后再乘1/N,即得 (n)。所 以FFT,IFFT可用一个子程序。 四、直接利用FFT流图方法的注意点 （1）FFT与IFFT连接应用时，注意输入输出序列的 排列顺序，是自然顺序还是倒序。 （2）FFT和IFFT共用同一个程序时，也应注意利用 FFT算法输入输出的排列顺序，即应注意自然 顺序还是例位序 （3）当把频率抽取FFT流图用于IDFT时，应改称时间 抽取IFFT流图。 （4）当把时间抽取FFT流图用于IDFT时，应改称频率 抽取IFFT流图。 §2-7 任意基数的FFT算法 以上讨论的都是以2为基数的FFT算法，即N=2M。 优点：程序简单，