华科信号课设三份 FFT 卷积 滤波器.docxVIP

《FFT的C语言实现》设计报告姓 名：江海啸学 号：U2010119511班 级： 电气1006专 业：电气工程及其自动化指导教师：李开成时 间：2012年八月至九月华中科技大学 电气与电子工程学院姓名：江海啸 班级：电气1006班时间：2012年八月至九月 编程软件：microsoft visual studio2010简介：FFT（Fast Fourier Transformation），即为快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。它对傅氏变换的理论并没有新的发现，但是对于在计算机系统或者说数字系统中应用离散傅立叶变换，可以说是进了一大步。本程序可以实现简单的FFT变换。它具有较简单的人机交互界面，可以选择FFT的变换点数（最低为4点，最高为512点），对输入不足的点自动补零，并返回输入的点数，计算结果较为精确。编程思想及原理：本程序主要由两部分组成，一部分是FFT变换函数，另一部分是实现人机交互的主函数，其中FFT主函数由复数加法模块，复数减法模块，复数乘法模块，倒序模块，蝶形计算模块组成。主要程序框图见图 1，FFT算法的原理介绍将不再在本篇报告里面赘述，详细原理见科学出版社出版的《信号与系统基础——应用Web和MATLAB（第三版）》第128面。图 1 系统整体框图模块介绍：主程序模块主程序程序框图见图 2，它的主要功能是人机交互，可以让用户选择FFT变换点数和输入需要变换的点数，并对未输入的点进行自动补零，并提示用户输入的点数。图 2 主函数程序框图FFT模块本模块实现了FFT的运算，先对输入数组进行倒序排列，倒序排列使用了雷德算法，然后再对重排以后的数列进行蝶形算法。下面先对雷德算法进行简要介绍。表 1倒序前后对比倒序排列前倒序排列后000（0）000（0）001（1）100（4）010（2）010（2）011（3）110（6）100（4）001（1）101（5）101（5）110（6）011（3）111（7）111（7）由 可以看出，排序以前的数最低位增1，排序以后的数最高位就增1，并且排序以后的数是向低位进位的，所以根据这个原理，雷德算法可以实现倒序排列。雷德算法的程序框图见图 3图 3 雷德算法程序框图本程序的蝶形算法采用迭代变换系数的方法，蝶形图见图 4图 4蝶形图 编程思想是：先计算蝶形级数，然后计算运算系数和一级里面蝶形结的宽度，再计算参与运算的两个数之间的距离（蝶形结宽度的一半），具有相同系数的运算先计算，然后改变系数重复运算，最终得到结果。实验结果：用户界面（图片经过反色处理），选择变换点数：提醒用户输入点数：输入点1,2,3,4后得出结果：改变运算点数后得出结果：输入点数不足时(用户可以在最后一个输入的数后加入“，0”选择退出输入)，会自动提示用户输入的点数:最后为了验证计算结果的准确性，我选择使用Matlab进行验证，验证结果见图 5，对比Matlab的结果可以知道该程序的结果十分精确。图 5MATLAB验证小结：经过这次课程设计，我对C语言的了解又更深了一步，也知道了FFT在C语言中是如何实现的，对FFT的理解也更加深入。参考书籍Edward W.Kamen Bonnie S.Heck著.信号与系统基础-应用Web和MATLAB（第三版）信号系统-英文.北京：科学出版社，2011丁毓峰 编著.MATLAB函数查询及应用案例.MATLAB软件.北京：电子工业出版社.2011年谭浩强 编著.C程序设计.C语言-程序设计.北京：清华大学出版社.2004年附录程序代码：#include stdio.h#include stdlib.h#include math.h#define PI 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068struct compx{float real,imag;};//定义复数结构体struct compx* s;//复数指针struct compx EE(struct compx a,struct compx b)//复数乘法{struct compx c;c.real=a.real*b.real-a.imag*b.imag;c.imag=a.real*b.imag+a.imag*b.real;return(c);}struct compx PP(struct compx a,struct compx b)//复数加法{struct comp