信号分析与处理实验报告程序+报告内容非常全.docVIP

实验四 离散傅里叶变换及其快速算法 一、实验目的 1、掌握计算序列的离散傅里叶变换（DFT）的方法。 2、掌握实现时间抽取快速傅里叶变换（FFT)的编程方法。 3、通过编写程序，深入理解快速傅里叶变换算法（FFT）的含义，比较DFT与FFT的运算效率。 4、加深对DFT与序列的傅里叶变换和Z变换之间关系的理解。 二、实验原理 2.1 离散傅里叶变换（DFT）算法 有限列长为的序列的变换为 其逆变换为 2.2 快速傅里叶变换基本思路 当很大的时候，求一个点的DFT要完成次复数乘法和次复数加法，其计算量相当大。快速傅里叶变换巧妙地利用因子的周期性和对称性，构造了一个DFT快速算法。其基本思路如下： 1、利用的对称性使运算中有些项合并 2、利用的周期性和对称性使长序列的分解为更小点数的 为了讨论方便，设，其中为整数。如果不满足这个条件，可以认为得加上若干零点来达到。由的定义知： 其中是列长为的输入序列，把它按的奇偶分成两个序列： 又由于， 则 上式表明了一个N点的DFT可以被分解为两个N/2点的DFT。同时，这两个N/2点的DFT按照上式又可以合成为一个N点的DFT。 为了要用点数为N/2点的X(k)、X(k)来表达N点的X(k)值还必须要用W系数的周期性，即 这样可得 即 同理可得 另外再加上W的对称性 就可以将X(k)的表达式分为前后两个部分： 前半部分 后半部分 由以上分析可见，只要求出区间内各个整数k值所对应的X(k)、X(k)的值，即可求出区间内的全部X（k）值，这一点恰恰是FFT能大量节省计算的关键所在。 2.3 时间抽取的快速傅里叶变换计算机算法 1、程序输入序列的元素数目必须为2的整数次幂，即,整个运算需要M级蝶形运算； 2、输入序列应该按二进制的码位倒置排列，输出序列按自然序列排列； 每个蝶形运算的输出数据占用其他输入数据的存储单元，实现“即位运算”； 每级包括N/2个基本蝶形运算，共有M*N/2个基本蝶形运算； 3、第L级中有 个群，群与群的间隔为 ； 4、同一级的各个群的系数Wn分布相同，第L级的群中有 个系数； 5、处于第L级的群的系数是 (p=1,2,3，……., )； 6、对于第L级的蝶形运算，两个输入数据的间隔为。 三、实验程序及注释 3.1 用C语言实现FFT与IFFT #include stdio.h #include math.h #include stdlib.h /*清屏命令clrscr()需加载的预处理函数*/ #define swap(a,b) temp=(a);(a)=(b);(b)=temp /*宏定义的交换函数*/ /*快速傅利叶变换程序，数组A、B分别是带变换序列的实部和虚部ap是变换的类别，ap=1是做FFT; ap=-1是做IFFT */ void fft(float A[],float B[],unsigned M,int ap) {unsigned long N,I,J,K,L,LE,LE1,P,Q,R; int i; float Wr,Wi,W1r,W1i,WTr,WTi,theta,Tr,Ti,temp; N=pow(2,M); /*N=2M 是序列的总长度*/ J=0; for(I=0;IN-1;I++) /*码位倒置*/ {if(JI) {swap(A[I],A[J]); swap(B[I],B[J]);} K=N1; /*即k=k/2 */ while(K=2J=K) {J-=K; K=1;} J+=K;} for(L=1;L=M;L++) /*外层循环由级数L控制，执行M次*/ {LE=1L; /* LE=2L是群间隔*/ LE1=LE/2; /* LE1=2L-1是每个群的系数W数目*/ Wr=1.0; Wi=0.0; theta=(-1)*ap*3.1415926536/LE1; /*ap用于在IFFT时改变 极性*/ W1r=cos(theta); W1i=sin(theta); for(R=0;RLE1;R++) /*中层循环由群系数 控制，执行 次*/ {for(P=R;PN-1;P+=LE) /*R是群系数的编号，P、Q是基本蝶形运 算两个输入数据在数组中的编号，循环每 次完成同一个系数 的蝶形运算*/ {Q=P+LE