16点基4dit-fft程序c++

//2011-11-04
//general
//points 16 radix4 dit-fft




#include 
using namespace std;

int main（）
{
float dataR[16]dataI[16];

int x0x1x2x3;
int Ljkbpi;
float TRTITBRTBIT2BRT2BIT3BRT3BI;

float COS[10]SIN[10];

int xx;
cout<<“16 points fft “<int a[16];


for（i=0;i<16;i++）
{
	a[i]=0;
}

for（i=0;i<16;i++）
{
	a[i]=i;
	cout<<“a[“<}



/////////////////////////////



cout<<“\n“<<“invert sequence: “<
///////////////////////////
for（i=0;i<16;i++）
{
	dataR[i]=a[i];
}


////////////////////////////////////////////wnp table

{
COS[0]=1;SIN[0]=0;
COS[1]=0.92388f;SIN[1]=-0.382683f;
COS[2]=0.707107f;SIN[2]=-0.707107f;
COS[3]=0.382683f;SIN[3]=-0.92388f;
COS[4]=0;SIN[4]=-1;
COS[5]=-0.382683f;SIN[5]=-0.92388f;
COS[6]=-0.707107f;SIN[6]=-0.707107f;
COS[7]=-0.92388f;SIN[7]=-0.382683f;
COS[8]=-1;SIN[8]=0;
COS[9]=-0.923879f;SIN[9]=0.382684f;


}



/********** following code invert sequence ************/
for（i=0;i<16;i++）
{ x0=x1=x2=x3=0;
x0=i&0x01; x1=（i/2）&0x01; x2=（i/4）&0x01; x3=（i/8）&0x01;
xx=x1*8+x0*4+x3*2+x2;
dataI[xx]=dataR[i];

}


for（i=0;i<16;i++）
{ dataR[i]=dataI[i]; dataI[i]=0; 


}
for（i=0;i<16;i++）
{
	cout<<“iv[“<}
cout</************** following code FFT *******************/
for（L=1;L<=2;L++） 
{ /* for（1） */
	b=1; i=L-1;
	cout<<“\n“<<“L=“<	w