FFT算法设计(含程序设计).pptVIP

FFT算法设计;傅立叶变换将信号从时域转换为频域，可以进行模拟信号的频率分析 离散傅立叶变换(DFT)将信号从频域转换为数字(频)域，可以进行数字信号（模拟信号数字化）的频率分析 为了实现DFT在计算机上的快速实现，提出了快速离散傅立叶变换(FFT);如何有傅氏变换-DFT-FFT?;;推导分析;;;蝶形分解图示;N=8点FFT运算图示;N=16点FFT运算图示;蝶形运算规律;;;公式中的J就是流程图中公式的变量k， 流程图中： N表示阶数， M表示总级数， L表示当前级数， B表示每个蝶形的两个输入数据的间隔， P表示旋转因子指数;;看图推导软件编程规则：方法一;;;;看图推导方法二：;;外循环次数为级数L 中循环为根据当前L求出各个不同的p，循环次 数为p的个数2L-1 内循环为每级中每个p对应的蝶形运算个数(记为VTotal)，循环次数为2M-L ;每个蝶形的两个输入数据间隔（记为INd）：;输入序列倒序的算法设计;倒序规律;; //输入序列倒序软件程序 j = N / 2; //第0个数(二进制数都为0)和最后一个第N-1个数(二进制数都为1)不需倒序 for(i = 1; i N - 2; i++) { if(i j) { temp = dataR[i]; dataR[i] = dataR[j]; dataR[j] = temp; } k = N / 2; while(1) { if(j k) { j = j + k; break; } else { j = j - k; k = k / 2; } } };输入序列倒序的算法设计方法二;方法二软件分析：已一个字节(N=256)的倒序为例; F0=i0x04; F1=i0x20; if(F0) j=j|0x20; if(F1) j=j|0x04; F0=i0x08; F1=i0x10; if(F0) j=j|0x10; if(F1) j=j|0x08; if(ij) //前半部与后半部交换，相等时无需交换 A[i] ? A[j] ; } ;算法改进一：;算法改进二：针对任意N=2M的情况;FFT软件示例; //输入序列倒序 j = N / 2; //第0个数(二进制数都为0)和最后一个第N-1个数(二进制数都为1)不需倒序 for(i = 1; i N - 2; i++) { if(i j) { temp = dataR[i]; dataR[i] = dataR[j]; dataR[j] = temp; //因为波形虚数部分都为0,所以不用交换 //temp = dataI[i]; //dataI[i] = dataI[j]; //dataI[j] = temp; } k = N / 2; while(1) { if(j k) { j = j + k; break; } else { j = j - k; k = k / 2; } } }; //进行FFT //Xr[J] = Xr(J) + Tr //Xi[J] = Xi(J) + Ti //Xr[J+B] = Xr(J) - Tr //Xi[J+B] = Xi(J) - Ti //(其中 Xr为上一级的Xr, Xi为上一级的Xi) //其中Tr = Xr(J+B)cos(2.0*PI*p/N) + Xi(J+B)sin(2.0*PI*p/N) // Ti = Xi(J+B)cos(2.0*PI*p/N) - Xr(J+B)sin(2.0*PI*p/N) for(L=1; L=M;L++) //FFT蝶形级数L从1--M { //计算每个蝶形的两个输入数据相距 B = 2^(L-1); B = 1; i = L-1; while(i) { B = B * 2; i--; }//或采用运算,即B = 2^(L-1); B = (int)pow(2,L-1); ; //第L级蝶形有pow(2,L-1),即2的L-1次方个蝶形运算和pow(2,L-1)个旋转因子p for(J=0;J=B-1;J++)