快速傅里变换FFT实验.docVIP

快速傅里叶变换FFT实验. 快速傅里叶变换FFT源程序： #includeiostream.h #includemath.h //using namespace std #define pi 3.1415926536 //定义pai void swap(double a,double b) //定义交换函数swap {double t=a; a=b; b=t; } void fft(double A[],double B[],int M) //定义fft快速傅里叶函数 { long int J,I,L,LE,LE1,P,Q,R; double N,K,Wr,Wi,W1r,W1i,WTr,WTi,theta,Tr,Ti; J=0; N=pow(2,M); //雷达算法码位倒置程序 for(I=0;IN-1;I++) { if(JI) { swap(A[I],A[J]); swap(B[I],B[J]); } K=N/2; //向低位进位 while(K=2J=K) { J-=K; K=K/2; } J+=K; } //蝶形运算程序 //最外层循环由蝶形图的级数确定 for(L=1;L=M;L++) { LE=pow(2,L); //各级群间隔 LE1=LE/2; //输入间隔 Wr=1.0; //Wn的实部 Wi=0.0; //Wn的虚部 theta=(-1)*pi/LE1; W1r=cos(theta); W1i=sin(theta); //中层循环由Wn的个数确定 for(R=0;RLE1;R++) { //内层循环由每一级群的个数确定 for(P=R;PN-1;P+=LE) { Q=P+LE1; Tr=Wr*A[Q]-Wi*B[Q]; Ti=Wr*B[Q]+Wi*A[Q]; A[Q]=A[P]-Tr; B[Q]=B[P]-Ti; A[P]+=Tr; B[P]+=Ti; } //计算每一级的Wn的值 WTr=Wr; WTi=Wi; Wr=WTr*W1r-WTi*W1i; Wi=WTr*W1i+WTi*W1r; } } } void main() { int M; double N; cout实验人的信息： 姓名：张超华 学号：200801000728 专业班级：电力电子0801endl; cout输入序列的个数N(N必须为2^M)，N=; cinN; //判断输入N的个数是否符合2^M个 while(1) { double b=log10(N)/log10(2); int a=b/1; if(a-b==0) { M=a; break; } else { cout输入序列的个数须为2^M个,重新输入N=; cinN; } } double *A=new double[N]; //为序列的实部分配动态空间 double *B=new double[N]; //为序列的虚部分配动态空间 for(int i=0;iN;i++) { A[i]=0.0; B[i]=0.0; } // 输入序列 cout请在下方输入序列的实部：endl; for(int j=0;jN;j++) { cinA[j]; } cout请在下方输入序列的虚部：endl; for(j=0;jN;j++) { cinB[j]; } fft(A,B,M); //调用fft函数 cout经快速傅里叶变换后得到的结果为：endl; //输出变换结果 for(int m=0;mN;m++) { cout A[m]+j(B[m])endl;