//----------------------------------------------------------------------
//?按頻率抽取的基4-FFT
//----------------------------------------------------------------------
#include?
#include?

#include?“fft_4_float_DIF_ok.h“
//----------------------------------------------------------------------
//尋找對應當前下標的反轉下標，即倒位序
int?FindIndex12bit（int?k）{
	int?j;
	k=k&（0xfff）;
	j=（（k&0x003）<<10）+（（k&0x00c）<<6）+（（k&0x030）<<2）+（（k&0x0c0）>>2）+（（k&0x300）>>6）+（（k&0xc00）>>10）;
	return?j;
}

int?FindIndex10bit（int?k）{
	int?j;
	k=k&（0x3ff）;
	j=（（k&0x003）<<8）+（（k&0x00c）<<4）+（k&0x030）+（（k&0x0c0）>>4）+（（k&0x300）>>8）;
	return?j;
}

uchar?FindIndex8bit（uchar?k）{
	uchar?j;
	k=k&（0xff）;
	j=（（k&0x003）<<6）+（（k&0x00c）<<2）+（（k&0x030）>>2）+（（k&0x0c0）>>6）;
	return?j;
}

uchar?FindIndex6bit（uchar?k）{//64點
	uchar?j;
	k=k&（0x3f）;//去除高位
	j=（（k&0x003）<<4）+（k&0x00c）+（（k&0x030）>>4）;
	return?j;
}

uchar?FindIndex4bit（uchar?k）{//16點
	uchar?j;
	k=k&（0x3f）;//去除高位
	j=（（k&0x003）<<2）+（（k&0x00c）>>2）;
	return?j;
}

//反轉下標
void?ReverseIndex（complex?*x）{
int?ij;
complex?tmp;
for?（i=0;i	j=FindIndex（i）;
	if（j>i）{
		tmp=x[i];
		x[i]=x[j];
		x[j]=tmp;
		}
	}
}
//----------------------------------------------------------------------
void?fft4（complex?*xint?M）
{
	int?ijktLi1i2i3nn1B;
	double?pe;
	complex?y[4]temp[3]wn[3];
	n=1<<（2*M）;
	B=n;
	for（L=1;L<=M;L++）	//第一層循環k控制級數
	{
		n1=B;			//n1=n相當于4^M-L（0……N-1）
		B=B/4;			//B=4^m-k，相當于4^M-L-1（0……N/4-1）
		e=PI_2/n1;		//e=2PI/N
		p=0;			//p為旋轉因子的指數，即p
		for（j=0;j		{				//第二層循環，更新j值，j表示各個子塊的第j個蝶形。因為每個子塊的同地位蝶形具有相同的wn，所以用第二層循環控制wn
			wn[0].real=cos（p）;	//旋轉因子的計算，第一次為0，所以e不必參與進來
			wn[1].real=cos（2*p）;
			wn[2].real=cos（3*p）;
			wn[0].imag=sin（p）;
			wn[1].imag=sin（2*p）;
			wn[2].imag=sin（3*p）;
			p=（j+1）*e;	//下一組的旋轉因子指數
			for（i=j;iN）。
			{
				i1=i+B;	//每個蝶形的4個輸入數據間距B
				i2=i1+B;
				i3=i2+B;
				
				y[0].real=x[i].real+x[i1].real+x[i2].real+x[i3].real;//計算一個DFT的第一個分式的實部
				y[0].imag=x[i].imag+x[i2].imag+x[i1].imag+x[i3].imag;//計算一個DFT的第一個分式的虛部

				temp[0].real=x[i].real-x[i1].imag-x[i2].real+x[i3].imag;//計算一個DFT的第二個分式的實部
				temp[1].real=x[i].real-x[i1].real+x[i2].real-x[i3].real;//計算一個DFT的第三個分式的實部
				temp[2].real=x[i].real+x[i1].imag-x[i2].real-x[i3].imag;//計算一個DFT的第四個分式的實部
				
				temp[0].imag=x[i].imag+x[i1].real-x[i2].imag-x[i3].real;//計算一個DFT的第二個分式的虛部
				temp[1].imag=x[i].imag-x[i1].imag+x[i2].imag-x[i3].imag;//計算一個DFT的第三個分式的虛部
				temp[2].imag=x[i].imag-x[i1].real-x[i2].imag+x[i3].real;//計算一個DFT的第四個分式的虛部
				
				y[1].real=wn[0].real*temp[2].real+wn[0].imag*temp[2].imag;//第二個分式的實部進行加權計算
				y[1].imag=wn[0].real*temp[2].imag-wn[0].imag*temp[2].real;//第二個分式的虛部進行加權計算
				y[2].real=wn[1].real*temp[1].real+wn[1].imag*temp[1].imag;//第三個分式的實部進行加權計算
				y[2].imag=wn[1].real*temp[1].imag-wn[1].imag*tem

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????35328??2010-12-14?13:55??FFT_4_float_DIF_（OK）\Debug\FFT_4_float_DIF_（OK）.exe

?????文件?????382420??2010-12-14?13:55??FFT_4_float_DIF_（OK）\Debug\FFT_4_float_DIF_（OK）.ilk

?????文件?????437248??2010-12-14?13:55??FFT_4_float_DIF_（OK）\Debug\FFT_4_float_DIF_（OK）.pdb

?????文件???????7392??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\BuildLog.htm

?????文件????????663??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\FFT_4_float_DIF_（OK）.exe.embed.manifest

?????文件????????728??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\FFT_4_float_DIF_（OK）.exe.embed.manifest.res

?????文件????????621??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\FFT_4_float_DIF_（OK）.exe.intermediate.manifest

?????文件??????16108??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\FFT_4_float_DIF_OK.obj

?????文件?????????60??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\mt.dep

?????文件??????60416??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\vc90.idb

?????文件??????61440??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug\vc90.pdb

?????文件???????4026??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）.vcproj

?????文件???????1413??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）.vcproj.AOTIAN-L.傲天.user

?????文件???????4624??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\FFT_4_float_DIF_OK.cpp

?????文件????????735??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\fft_4_float_DIF_ok.h

?????文件?????478208??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）.ncb

?????文件????????931??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）.sln

?????文件??????11776??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）.suo

?????目錄??????????0??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）\Debug

?????目錄??????????0??2010-12-14?13:55??FFT_4_float_DIF_（OK）\Debug

?????目錄??????????0??2010-12-14?13:55??FFT_4_float_DIF_（OK）\FFT_4_float_DIF_（OK）

?????目錄??????????0??2010-12-14?13:55??FFT_4_float_DIF_（OK）

-----------?---------??----------?-----??----

??????????????1504137????????????????????22