快速傅立叶变换 DSP.docVIP

DSP原理及应用 大作业 ——快速傅立叶变换 快速傅立叶变换 一、实验目的 1. 深入了解FFT 的原理及其算法； 2． 掌握用窗函数法设计FFT快速傅里叶的原理和方法； 3. 熟悉FFT快速傅里叶特性； 4. 研究采样点数对FFT的影响。 二、实验原理 快速傅里叶变换（FFT）是计算离散傅里叶变换（DFT）的快速算法。 DFT的定义式为 = FFT的基本思想：将大点数的DFT分解为若干个小点数DFT的组合，从而减少运算量。 因子具有以下两个特性，可使DFT运算量尽量分解为小点数的DFT运算： 周期性： 对称性： 利用这两个性质，可以使DFT运算中有些项合并，以减少乘法次数。 为了将大点数的DFT分解为小点数的DFT运算，要求序列的长度N为复合数，最常用的是的情况（M为正整数）。 一个N点DFT被分解为两个N/2点的DFT。 分析：和只有N/2个点（），则由只能求出的前N/2个点的DFT，要求出全部N点的，需要找出、和的关系，其中。由式子可得 ＝， 这样N点DFT可全部由下式确定出来： （＊） 上式可用一个专用的碟形符号来表示，这个符号对应一次复乘和两次复加运算。 图　蝶形运算符号 因为，N/2仍然是偶数，因此可以对两个N/2点的DFT再分别作进一步的分解，将两个N/2点的DFT分解成两个N/4点的DFT。 对于一个的DFT运算，其按时间抽取的分解过程及完整流图如下图所示。 三、实验内容 1、实验准备：设置软件仿真模式。 -启动CCS。 2、打开工程FFT.pjt，浏览程序。 3、编译并下载程序。 4、打开观察窗口：选择菜单View-Graph-Time/Frequency… 5、退出CCS。 四、实验程序 (1) FFT.c: #include myapp.h #include csedu.h #include scancode.h #include math.h #define PI 3.1415926 #define SAMPLENUMBER 128 void InitForFFT(); void MakeWave(); int INPUT[SAMPLENUMBER],DATA[SAMPLENUMBER]; float fWaveR[SAMPLENUMBER],fWaveI[SAMPLENUMBER],w[SAMPLENUMBER]; float sin_tab[SAMPLENUMBER],cos_tab[SAMPLENUMBER]; void FFT(float dataR[SAMPLENUMBER],float dataI[SAMPLENUMBER]); main() { int i; InitForFFT(); MakeWave(); for ( i=0;iSAMPLENUMBER;i++ ) { fWaveR[i]=INPUT[i]; fWaveI[i]=0.0f; w[i]=0.0f; } FFT(fWaveR,fWaveI); for ( i=0;iSAMPLENUMBER;i++ ) { DATA[i]=w[i]; } while ( 1 ); // break point } void FFT(float dataR[SAMPLENUMBER],float dataI[SAMPLENUMBER]) { int x0,x1,x2,x3,x4,x5,x6,x7,xx; int i,j,k,b,p,L; float TR,TI,temp; /********** following code invert sequence ************/ for ( i=0;iSAMPLENUMBER;i++ ) { x0=x1=x2=x3=x4=x5=x6=0; x0=i0x01; x1=(i/2)0x01; x2=(i/4)0x01; x3=(i/8)0x01;x4=(i/16)0x01; x5=(i/32)0x01; x6=(i/64)0x01; x7=(i/128)0x01,xx=x0*128+x1*64+x2*32+x3*16+x4*8+x5*4+x6*2+x7; dataI[xx]=dataR[i]; } for ( i=0;iSAMPLENUMBER;i++ ) { dataR[i]=dataI[i]; dataI[i]=0; } /************** following code FFT *******************/ for ( L=1;L=8;L++ ) { /* for(1) */ b=1; i=L-1