快速傅里叶算法C语言完整实现

Invert_order(N,x,M);
printf("dao xu hou xu lie:\n");
for(i=0;i<DATA_LEN;i++)//倒序后再次输出
{
printf("%lf %lf ",x[i].Re,x[i].Im);
if(i%2==0)printf("\n");
x[k+Bank]=Sub(x[k],x_temp);
x[k].Re=temp.Re; x[k].Biblioteka m=temp.Im; }
}
//printf("zhu yi le:clk=%d\n",clk);
}
}
/******************************************************************************************
printf("%f\t%f\t\n",x[i].Re,x[i].Im);
// printf("%f\n",sqrt(x[i].Re*x[i].Re+x[i].Im*x[i].Im));
}
// printf("\n");
}
#endif
for(j=0;j<=Bank-1;j++)
{
Complex W_temp;
W_temp.Im=-1.0;
W_temp.Re=1.0;
r=Pow(2,M-L)*j ;
ps=2*pai/N*r ;
W.Re=cos(ps);
W.Im=-sin(ps);
// W.Re=ps;
#include<stdio.h>
#include<math.h>
#define uint unsigned int
#define uchar unsigned char
typedef struct Complex//复数结构体
{
double Re;
double Im;
}Complex;
log函数结束：
******************************************************************************************/
/******************************************************************************************
void Invert_order(uint N,Complex *x,uchar M)
{
uint power,i,j,m,tempp;//m是i的倒位序数
Complex temp;//是临时变量
for(i=0;i<N;i++)
{
power=1;
m=0;
for(j=0;j<M;j++)
函数说明：按时间抽取基二快速傅里叶算法
参数说明：N点数，等于2的M次方，*x序列的首地址
返回值：无
******************************************************************************************/
void Dit2Fft(uint N,Complex *x)
int i;
Complex x[DATA_LEN];
N=DATA_LEN;
M=Log(2,N);
printf("chu shi xu lie:\n");
for( i=0;i<DATA_LEN;i++)//序列初始化
{
x[i].Im=0;
x[i].Re=i;
}
printf("\n");
}
printf("\n");
Dit2Fft(N,x);
printf("FFT yun suan hou xu lie:\n");
for(i=0;i<DATA_LEN;i++)//FFT后再次输出
{
// printf("%lf+%lfi\t\t",x[i].Re,x[i].Im);
{
uint m=1;//指数函数值
uchar i;
for(i=0;i<IndexNumber;i++)
{
m=BaseNumber*m;
}
return m;
}
/******************************************************************************************
函数说明：指数运算
参数说明：BaseNumber底数,IndexNumber指数
返回值：无符号整数
******************************************************************************************/
uint Pow(uchar BaseNumber,uint IndexNumber)
//printf("%lf\t",x[i].Im);
// if(i%2==0)printf("\n");
}
// printf("\n");
printf("After...\nReal\t\tImag\n");
for(i=0;i<DATA_LEN;i++)//FFT后再次输出
{
}
/******************************************************************************************
函数说明：倒序运算
参数说明：N点数，等于2的M次方，*x序列的首地址
返回值：无
******************************************************************************************/
{
int clk=0;
Complex temp;//临时复数变量
uint L,M=Log(2,N);//M级蝶形图
uint k,j;
uint StepLength;//k值步长
uint Bank ;//两点间的距离
const double pai=3.1415926;
double ps;
函数说明：对数运算
参数说明：BaseNumber底数,AntiNumber真数
返回值：无符号整数
******************************************************************************************/
uint Log(uchar BaseNumber,uint AntiNumber)
{
uint m=0;//对数函数值
while(1)
{
AntiNumber=AntiNumber/BaseNumber;
if(AntiNumber)m++;
else break;
}
return m;
}
/******************************************************************************************
// W.Im=10.0*ps;
for(k=j;k<=N-1;k=k+StepLength)
{
Complex x_temp;
printf("%u\n",r);
x_temp=Mul(W,x[k+Bank]);
temp=Plus(x[k],x_temp);
Dit2Fft函数结束：
******************************************************************************************/
#if 1
#define DATA_LEN 16
void main()
{
uint N,M;
uint r;//旋转因子的幂
Complex W;//旋转因子
printf("M de zhi shi :M=%u\n",M);
for(L=1;L<=M;L++)
{
StepLength=Pow(2,L);
Bank=StepLength/2;
printf("L=%u\tStep=%u\tBank=%u\n",L,StepLength,Bank);
Pow函数结束：
******************************************************************************************/
/******************************************************************************************