第5讲-DSP应用系统设计.pptVIP

N=8的FFT运算举例 （2）快速傅里叶变换FFT 例：时间抽取的FFT算法DSP C语言实现实例。 FFT运算函数与主函数为 #includemath.h //数学函数头文件 #define PI 3.1415926 #define N 128 //采样次数N void InitForFFT( ); //FFT初始化函数 void MakeWave( ); //波形发生函数 void finv(int N1, float *xr, float *xi); //倒序运算函数f(N1,Xr,Xi)，对输入序列倒序 int INPUT[N], DATA[N]; float fWaveR[N], fWaveI[N], w[N]; float sin_tab[N],cos_tab[N]; //正余弦函数表 int Mum; //Mum为蝶形运算的级数 （3） FFT的C语言例程 * void FFT(float Xr[N],float Xi[N]) //时间抽取法FFT程序，要求采样点数N为2的整数幂次方 { //Xr[ ], Xi[ ]分别为输入序列的实部和虚部 int S,B; //S为旋转因子的幂数, B为蝶形运算输入数据的距离， //也即各级旋转因子的个数 int m, j, k; float X,Y; finv(N, Xr, Xi); //倒序运算函数，对输入序列倒序 for (m=1; m=Mum; m++) { B=(int)(pow(2,m-1)+0.5); //B=2^(m-1) for(j=0; jB; j++) //每级需要进行B种蝶形运算 { S=j*(int)(pow(2,Mum-m)+0.5); //S=2^(Mum-1) for(k=j; k=N-1; k+=(int)(pow(2,m)+0.5)) （3） FFT的C语言例程 //每种蝶形运算在某一级中需要进行N/pow(2,m)次 { //蝶形运算展开，结果的实部和虚部 //分别存储在原实部和虚部位置 X=Xr[k+B]*cos_tab[S]+Xi[k+B]*sin_tab[S]; Y=Xi[k+B]*cos_tab[S]-Xr[k+B]*sin_tab[S]; Xr[k+B]=Xr[k]-X; Xi[k+B]=Xi[k]-Y; Xr[k]=Xr[k]+X; Xi[k]=Xi[k]+Y; } } } for ( m=0;mN/2;m++ ) { w[m]=sqrt(Xr[m]*Xr[m]+Xi[m]*Xi[m]); //计算功率普 } } （3） FFT的C语言例程 main() { int i; InitForFFT( ); //FFT初始化函数 MakeWave( ); //波形发生函数 for ( i=0;iN;i++ ) { fWaveR[i]=INPUT[i]; fWaveI[i]=0.0; w[i]=0.0; } Mum=(int)(0.5+log(N)/log(2)); //Mum为蝶形运算的级数, N=2^Mum FFT(fWaveR,fWaveI); for ( i=0;iN;i++ ) DATA[i]=w[i]; while ( 1 ); } （3） FFT的C语言例程 void InitForFFT() //FFT初始化函数，建立正余弦函数表 { int i; for ( i=0;iN;i++ ) { sin_tab[i]=sin(PI*2*i/N); cos_tab[i]=cos(PI*2*i/N); } } void MakeWave() //波形发生函数 { int i; for ( i=0;iN;i++ ) { INPUT[i]=sin(PI*2*i/N*3)*1024; //f=3Hz, 正弦函数 } } （3） FFT的C语言例程 FFT函数包含的函数finv(N,Xr,Xi)为倒序运算，函数代码如下。 //倒序运算函数finv(N1,Xr,Xi)，对输入序列倒序 //N1为序列长度; Xr[ ], Xi[ ]分别为输入序列的实部和虚部 //倒序原理：倒序数的加1是在最高位加1， //满2向次高位进1，最高位变0，依次往下 //从当前倒序值可求下一倒序值 void finv(int N1, float *xr, float *xi) //倒序运算函数f(N1,Xr,Xi)，对输入序列倒序 { int m, n,N2,k; // m为正序数；n为到序数；k为各个权值；N2为最高位的权值 float T; //临时变量T N2=N1/2; //最高位加1相当于十进制加上最高位的权N1/2