第二章数字信号处理.ppt

表码位倒置顺序自然顺序二进码表示码位倒置倒置顺序0000000010011004201001023011110641000011510110156110011371111117在实际运算中，一般直接将输入数据x（n）按码位倒置的顺序排好输入很不方便，总是先按自然顺序输入存储单元，然后再通过变址运算将自然顺序的存储转换成码位倒置顺序的存储，然后进行FFT的原位计算。目前有许多通用DSP芯片支持这种码位倒置的寻址功能。（4）蝶形类型随迭代次数成倍增加观察8点FFT的三次迭代运算:第一级迭代，有一种类型的蝶形运算系数W08，两个数据点间隔为1第二级迭代，有二种类型的蝶形运算系数W08、W28，参加运算的两个数据点间隔为2。第三级迭代，有四类蝶形运算系数W08、W18、W28、W38，参加运算的两个数据点间隔为4。结论：每迭代一次，蝶形类型增加一倍，数据点间隔也增大一倍。每一级的取数间隔和蝶形类型种类均为2i-1，i=1,2,…M。回节首3、按频率抽取的FFT对于N=2M情况下的另外一种普遍使用的FFT结构是频率抽取法。对于频率抽取法，将输入序列是按前后对半分开：将X（k）按k分解为偶数组和奇数组：偶数组奇数组令a(n)=x(n)+x(n+N/2)b(n)=[x(n)-x(n+N/2]WnN都是N/2点的序列两个N/2点的DFT运算将一个N点的DFT分解为两个N/2点的DFT上式的运算关系可用下图表示：以N=8的频率抽取为例按频率抽取将8点DFT分解成四个2点DFTN=8的按频率抽取FFT运算流图一个N=2M点的DFT通过M次分解后，最后只剩下全部是2点的DFT。2点DFT实际上只有加减运算。频率抽取法的流图是时域抽取法流图的左右翻转。频率抽取法FFT的运算特点：（1）蝶形运算对于任何一个2的整数幂N=2M，总是可以通过M次分解最后完全成为2点的DFT运算。这样的M次分解，就构成从x(n)到X(k)的M级运算过程。将频率抽取法的流图反转，并将输入变输出，输出变输入，得到的便是时间抽取法流图（反映了时域，频域的对称性）。频率抽取法也共有M级运算（N=2M），其运算量与时间抽取法相同。（2）原位计算类似于时间抽取法，当数据输入到存储器中以后，每一级运算的结果仍然储存在同一组存储器中，直到最后输出，中间无需其它存储器，所以频域抽取法也可进行原位计算。（4）蝶形类型随迭代次数成倍减少(与时间抽取法相反)第一级迭代中有N/2种蝶形运算系数，参加蝶形运算的两个数据相隔N/2，随后每次迭代，蝶形类型比前一级减少一倍，间距也减少一倍，最后一级迭代，蝶形类形只有一种W0N，数据间隔为1。由这几点规律可以看出:频率抽取法与时间抽到法是两种等价的FFT运算。回节首（3）序数重排它的输入正好是自然顺序。但它的输出却是码位倒置的顺序。因此运算完毕后，要通过变址运算将码位倒置的顺序转换为自然顺序，然后输出，变址方法同时间抽取法。4.N为组合数的FFT（任意基数的FFT算法）FFT的基本思想是将DFT的运算尽量分小,因此,在N=PQ情况下,将N点的DFT分解为P个Q点DFT(或Q个P点DFT)以减少计算量。步骤：分别为0,1,…，Q-1;分别为0,1,…，P-1。N点DFT写为如N为复合数，可分解为两个整数P与Q的乘积。令n0为参变量Q个P点DFT令恢复X(k)的正常顺序分别为0,1,…，Q-1;分别为0,1,…，P-1。例P=3，Q=4,N=12(1)先将x(n)通过x(n1Q+n0)改写成x(n1,n0)。因为Q=4,n1=0,1,2,n0=0,1,2,3,故输入是按自然顺序的，即x(0,0)=x(0)x(0,1)=x(1)x(0,2)=x(2)x(0,3)=x(3)x(1,0)=x(4)x(1,1)=x(5)x(1,2)=x(6)x(1,3)=x(7)x(2,0)=x(8)x(2,1)=x(9)x(2,2)=x(10)x(2,3)=x(11)(2)求Ｑ个Ｐ点的DFTn0为参变量Q个P点DFTX1(0,0)X1(0,1)X1(0,2)X1(0,3)X1(1,0)X1(1,1)