华科机械工程测试信息信号分析课件ch6-03数字信号分析讲解.ppt

此FFT算法，是在时间上对输入序列的次序是属于偶数还是属于奇数来进行分解的，所以称作按时间抽取的算法(DIT)。 二、运算量 由上述分析可知，N=8需三级蝶形运算N=2 =8，由此可知，N=2L 共需L级蝶形运算,而且每级都由N/2个蝶形运算组成，每个蝶形运算有一次复乘，两次复加。 3 因此，N点的FFT的运算量为： 复乘: mF=(N/2)L=(N/2)log2N 复加: aF=N L=Nlog2N 由于计算机的乘法运算比加法运算所需的时间多得多，故以乘法运算作为比较基准。 三、DIT的FFT算法的特点 -1 -1 -1 -1 -1 -1 -1 -1 . . . . . 1.原位运算 输入数据、中间运算结果和最后输出均用同一存储器。 x(0)=X0(0) x(4)=X0(1) x(2)=X0(2) x(6)=X0(3) x(1)=X0(4) x(5)=X0(5) x(3)=X0(6) x(7)=X0(7) X2(0) X2(1) X2(2) X2(3) X2(4) X2(5) X2(6) X2(7) X3(0)=X(0) X3(1)=X(1) X3(2)=X(2) X3(3)=X(3) X3(4)=X(4) X3(5)=X(5) X3(6)=X(6) X3(7)=X(7) X1(0) X1(1) X1(2) X1(3) X1(4) X1(5) X1(6) X1(7) 由运算流图可知，一共有N个输入/出行，一共有log2N=L列(级)蝶形运算(基本迭代运算)。 设用m(m=1,2, … ,L)表示第m列；用k,j表示蝶形输入数据所在的(上/下)行数(0,1,2,… ,N-1)；这时任何一个蝶形运算可用下面通用式表示，即： 所以，当m=1时，则有（前两个蝶形） 当m=2时，则有(前两个蝶形) 当m=3时，则有(前两个蝶形) 可见，在某列进行蝶形运算的任意两个节点(行)k和j的节点变量Xm-1(k)，Xm-1(j)就完全可以确定蝶形运算的结果Xm(k)，Xm(j)，与其它行(节点)无关。 这样，蝶形运算的两个输出值仍可放回蝶形运算的两个输入所在的存储器中，即实现所谓原位运算。每一组(列)有N/2个蝶形运算，所以只需N个存储单元，可以节 省存储单元。 2. 倒位序规律 由图可知，输出X(k)按正常顺序排列在存储单元，而输入是按顺序： 这种顺序称作倒位序，即二进制数倒位。 n0=0，(偶) n1=0 n1 =1 n1 =0 n1 =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为例说明如下： n0=1，(偶) 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 二进制n2 n1 n0 倒位序二进制n0 n1 n2 倒位顺序 A(1) A(2) A(3) A(4) A(5) A(6) A(7) A(8) x(0) x(1) x(2) x(3) x(4) x(5) x(6) x(7) x(0) x(4) x(2) x(6) x(1) x(5) x(3) x(7) 变址处理方法 存储单元 自然顺序 变址 倒位序 4. 蝶形运算两节点的距离：2m-1 其中，m表示第m列，且m =1,… ,L 例如N=8=23，第一级(列)距离为21-1=1， 第二级(列)距离为22-1=2，