信号处理四大算法教程第6讲FFT算法设计.pptVIP

方法二软件分析：已一个字节(N=256)的倒序为例 A[0], A[1], ……., A[255] (下标从0000_0000变化到1111_1111) /* 定义两个标志位F0，F1 */ for(i=0;i≤(255/2);i++) //除2是因为只要判断前半部 { j=0; F0=i0x01; //取最低位 F1=i0x80; //取最高位 if(F0) j=j|0x80; //最低位换到最高位 if(F1) j=j|0x01; //最高位换到最低位 F0=i0x02; //取次低位 F1=i0x40; //取次高位 if(F0) j=j|0x40; //最次位换到次高位 if(F1) j=j|0x02; //最次位换到次低位 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] ; } 只需单层循环即可，共需要循环（128-2）次 算法改进一： 前面的算法可以进一步优化为： for(i=0;i≤(255/2);i++) { j=0; for(k=0;k4;k++) { F0=i(0x01k); //取最低位 F1=i(0x80k); //取最高位 if(F0) j=j|(0x80k); //最低位换到最高位 if(F1) j=j|(0x01k); //最高位换到最低位 } if(ij) //前半部与后半部交换，相等时无需交换 A[i] ? A[j] ; } 这个算法只是针对8位的，如果是任意N位的应该如何做？ 这里的N必须满足N=2M 算法改进二：针对任意N=2M的情况 for(i=0;i(N/2);i++) //或者for(i=0;i≤((N-1)/2);i++) { j=0; for(k=0;k(M/2);k++) //当N=256时，M=8 { m= 0x01k; n= 0x80k; F0=i m; F1=i n; if(F0) j=j | n; if(F1) j=j | m; } if(ij) //前半部与后半部交换，相等时无需交换 A[i] ? A[j] ; } FFT软件示例 #include math.h #define PI 3.1415926 #define N 128 #define M 7 #define A0 255.0 //定义输入波形的幅值 void main(void) { int i,j,k; int p,J,L,B; float SinIn[N]; float dataR[N],dataI[N]; float dataA[N]; float Tr,Ti,temp; //输入波形 for(i = 0; i N; i++) { SinIn[i] = A0 * (sin(2*PI*i/25)+sin(2*PI* i * 0.4 )); dataR[i] = SinIn[i]; //输入波形的实数部分 dataI[i] = 0; //输入波形的虚数部分 dataA[i] = 0; //输出波形的幅度谱为0 } //输入序列倒序 j = N / 2; //第0个数(二进制数都为0)和最后一个第N-1个数(二进制数都为1)不需倒序 for(i = 1; i N - 2; i++) { if(i j) { temp = dataR[i]; data