数字信号处理及Matlab实现第4章.pptVIP

因此,N点的FFT的运算量为 复乘: mF =（N/2）L=（N/2） log2 N 复加: aF =N L=N log2 N 由于计算机的乘法运算比加法运算 所 需的时间多得多，故以乘法作为比较 基准。 (0)=X0(0) X1(0) X2(0) X3(0)=X(0) (4)=X0(1) X1(1) X2(1) X3(1)=X(1) (2)=X0(2) X3(2)=X(2) (6)=X0(3) X3(3)=X(3) (1)=X0(4) X1(4) X2(4) X3(4)=X(4) (5)=X0(5) X3(5)=X(5) (3)=X0(6) X3(6)=X(6) (7)=X0(7) X1(7) X2(7) X3(7)=X(7) W W W W N 0 N 0 N 0 N 0 -1 -1 -1 -1 W W W W N 0 N 2 N 0 N 2 -1 -1 -1 -1 W W W W N N N N 0 1 2 3 . . . . . . . . . . . 三.DIT的FFT算法的特点 1.原位运算 输入数据、中间运算结果和最后输出均用同一存储器。 设用m(m=1,2,… ,L)表示第m列;用k,j表示蝶形 输入数据所在的（上/下）行数(0,1,2,… ,N-1);这时任何一个蝶形运算可用下面通用式表示， 即 由运算流图可知,一共有N个输入/出行,一共有log2 N=L列（级）蝶形运算(基本迭代运算). 所以,当m=1时,则有（前两个蝶形） 当m=2时,则有（前两个蝶形） 当m=3时,则有（前两个蝶形） 可见,在某列进行蝶形运算的任意两个节点(行)k和j的节点变量 就完全可以确定蝶形运算的结果 ,与其它行(节点)无关。 这样,蝶形运算的两个输出值仍可放回蝶形运算的两个输入所在的存储器中,即实现所谓原位运算。每一组(列)有N/2个蝶形运算,所以只需N个存储单元，可以节 省存储单元。 2 倒位序规律 由图可知,输出X(k)按正常顺序排 列在存储单元,而输入是按顺序: 这种顺序称作倒位序,即二进制数 倒位。 n =0 0 n =1 0 n =0 1 n =1 1 n =0 1 n =1 1 0 1 0 1 0 1 0 1 (n2) x(000) 0 乾 x(100) 4 兑 x(010) 2 离 x(110) 6 震 x(001) 1巽 x(101) 5 坎 x(011) 3 艮 x(111) 7 坤 （偶） （奇） 这是由奇偶分组造成的,以N=8为例 说明如下: 3.倒位序实现 输入序列先按自然顺序存入存储单 元,然后经变址运算来实现倒位序排列 设输入序列的序号为n,二进制为 (n2 n1 n0 )2 ,倒位序顺序用 表示,其倒位序二进制为(n0 n1 n2 )2 ,例如 ，N=8时如下表： 0 0 0 0 0 0 0 0 1 0 0 1 1 0 0 4 2 0 1 0 0 1 0 2 3 0 1 1 1 1 0 6 4 1 0 0 0 0 1 1 5 1 0 1 1 0 1 5 6 1 1 0 0 1 1 3 7 1 1 1 1 1 1 7 自然顺序n 二进制n n n 倒位序二进制n n n 倒位顺序n ^ 2 1 0