FFT的C语言编程归类.pdfVIP

FFT 的 C 语言编程 1.程序： #include stdio.h #include math.h #define re 0 /* re=0 ，用 re 表示实部 */ #define im 1 /* im=1 ，用 im 表示虚部 */ main() { float x[128][2],w[2],temp[2]; /* x[128][2] ： 复数变量； x[i][re] ： 第 i 个复数变量的实部； x[i][im] ： 第 i 个复数变量的虚部； w[2] ： 存储旋转因子 WN P，w[re] 、w[im] 分别代表旋转因子的实部和虚部； temp[2] ： 蝶形计算中的临时变量， temp[re] 、temp[im] 分别代表其实部 和虚部； */ float arg,wreal,wimag; /* arg 存储旋转因子指数 p( 数值上相差 - 2 π/N) 。wreal 存储 cos(arg) ， wimag 存储 -sin(arg) */ float tem,tr,ti; int L,M,B,j,i,k, N,N2; char c=i; scanf(%d %d,N,M); /* 输入复数信号长度 N，蝶形运算级数 M */ N2=N1; for(j=0;jN;j++) /* 输入复数信号实部和虚部 */ { scanf(%f %f,x[j][re],x[j][im]); } printf(\n); /* 输入倒序 */ for(j=0,i=1;iN-1;i++) { k=N2; while(k=j) { j=j-k; k=k1; } j=j+k; if(ij) { tr=x[j][re]; ti=x[j][im]; x[j][re]=x[i][re]; x[j][im]=x[i][im]; x[i][re]=tr; x[i][im]=ti; } } /*FFT 三重循环模块 */ for(L=1; L=M; L++) /* 逐级进行计算共 M 级*/ { B=1L-1; /* 第 L 级共有 B=2L-1 个不同的旋转因子 */ arg=-acos(-1)/B; /* 旋转因子初始化注释见结尾处 */ w[re]=cos(arg); w[im]=-sin(arg); for(j=0; jB; j++) /* j 代表第 L 级不同旋转因子的个数 */ { /* 旋转因子 */ arg=acos(-1)/B; /* arg= π/ Β */ wreal=cos(arg); wimag= -sin(arg); tem=w[re]*wreal-w[im]*wimag; w[im]=w[re]*wimag+w[im]*wreal; w[re]=tem; for(k=j; kN; k+=2*B) /* 第 L 级具有相同旋转因子蝶形计算，每个蝶形相距 2L=2B 个点 */ { temp[re]=x[k+B][re]*w[re]-x[k+B][im]*w[im]; temp[im]=x[k+B][im]*w[re]+x[k+B][re]*w[im]; x[k+B][re]=x[k][re]-temp[re]; x[k+B][im]=x[k][im]-temp[im]; x[k][re]=x[k][re]+temp[re]; x[k][im]=x[k][im]+temp[im]; /* 编写蝶形运算程序 */ /* 第 L