FFT之数组逆序C语言实现

前言：FFT算法是DFT的快速算法实现，FFT算法中有一个关键步骤就是将原始序列逆序。比如对于个点的FFT变换，首先对序列号进行编码。比如M=4阶，共N=16个点，编码为0000,0001,0010,0011，……，1111。至于逆序的过程，就是将当前编码从右向左书写重新编码，比如说，0001编码后为1000,对于1101编码后为1011，0011编码后为1100。逆序之后，需要把原来位置上序列号对应的值更换为更改序列号后对应的值，比如说，原始序列第2个位置即序列号0001，对应的值为50，第9个位置即序列号1000，对应的值为100，更换后的序列第2个位置值为100，第9个位置值为50。以下我们要探讨的程序，就是想要如何在C语言上实现上述过程。

具体代码如下：

#include<stdio.h>
#include<stdlib.h>
int main()//
{
	int X1[16]={0,50,2,3,4,5,6,7,8,100,10,11,12,13,14,15};//存储原始序列
	int FD[16];//存储逆序后的序列
	int p;
	for (int j=0;j<16;j++)
	{
		p=0;
		for (int i=0;i<4;i++)
		{
			if (j&(1<<i))
				p+=1<<(4-i-1);
		}
		FD[j]=X1[p];
	}
	for (int j=0;j<16;j++)
		printf("%d\n",FD[j]);
}

下面具体分析一下这段代码。0000的逆序还是本身，对于其它每一个原始序列号，由于本例中为4位编码，即每一个序列号编码的组成都可以由1000,0100,0010，和0001组合相加后形成，上面四个原码可以由1向左分别移3、2、1、0即1<<3,1<<2,1<<1,1<<0得到。举1101为例，1101=1000+0100+0001，逆序后1011=1000+0010+0001。所以本逆序算法的思路就是，假如原始序列号第i位上值为1，则逆序后的序列号上第4-i-1位上值为1。用原序列号和0001,0010,0100,1000相与后判断原序列号每一位的值是否为1，若为1，则逆序后第4-i-1位上值为1，将这些“1”组合后形成逆序后的序列号。