DSP综合性设计实验(考试用).docVIP

湖南科技科学院 DSP综合性设计实验 设计课题： IIR滤波器设计 电子工程系 电信1403 班 设计者： 屈黎黎 张纪双 学号： 201406001127 201406002345 目录 一．设计目的 3 二．设计要求 3 三．总体概要设计 3 3.1 设计要求 3 3.2 IIR滤波器原理 3 四．差分方程及系数 5 五．实验步骤 6 六．源代码流程图 8 6.1 程序流程图： 8 6.2 程序源代码： 8 七．实验结果及分析 9 八. 总结 10 九．参考文献和辅助软件工具 11 9.1 参考文献 11 9.2 辅助软件工具：MATLAB 11 9.2.1 软件基本介绍 11 9.2.2 IIR滤波器的MATLAB实现 12 DSP综合性设计实验成绩评定 13 IIR滤波器设计 设计目的 运用之前所学的基础知识和基本理论，掌握设计 IIR 数字滤波器的原理和方法；熟悉 IIR 数字滤波器特性；了解 IIR 数字滤波器的设计方法。增强分析问题和解决问题的能力，加深对所学知识的理解，为今后的工作打下良好的基础。 二．设计要求 1．设计指标 2．设计要求 3．制作要求 自行装配和调试，并能发现问题和解决问题。 4．编写设计报告 写出设计与制作的全过程，附上有关资料和图纸，有心得体会。 三．总体概要设计 3.1 设计要求 本设计采用Butterworth型IIR低通滤波，其通带边缘 1kHz 处的增益为-3dB，10kHz 处的阻带衰减为60dB，采样频率 25kHz。 3.2 IIR滤波器原理 FIR滤波器的系统函数只有零点。而IIR滤波器除了具有极点以外，一般还存在零点。由于极点的存在，IIR滤波器用递归结构来实现较为简单。实现IIR滤波器的基本结构共有三中：直接型、级联型和并联型。下面简单介绍前两种IIR滤波器的结构。 直接型 利用公式 可以直接导出直接I型的IIR滤波器结构，可用下式来表示 可以看出，和反馈环节中的延时单元不能共用，需要个延时单元。为了减少延时单元，对于线性系统，上述公式也可以写成： 由此可以导出直接Ⅱ型的滤波器结构。图1表示了一个直接Ⅱ型IIR滤波器的结构，该滤波器的阶数是2阶，共存在两个反馈环节。 图1 直接Ⅱ型IIR滤波器结构 2、 级联型 对式(1-1)的分子、分母进行因式分解，由于中的系数都为实数，的极、零点只可能是实数或者复共轭对，对于复共轭对因子，可以复合成二阶因子： 式中系数都为实数。如果把实数因子(一阶因子)看成是二次项系数为0的二阶因子，则上式可以写成： 其中 由此可见，是一个2阶的IIR滤波器，N阶IIR滤波器可以看成是由多个2阶IIR滤波器级联而成的，即IIR滤波器可以采用下列结构来表示，也就是级联型结构(见图2)： 图2 级联型IIR滤波器结构 四．差分方程及系数 由双线性变换定义可得到IIR滤波器的系统传输函数表示为 IIR滤波器一般采用递归方式来实现。也就是说，滤波器当前输出是输入序列和以前各输出值、、…的函数，这可以用下列差分方程来表示： 运用辅助软件工具Matlab中的function模块将滤波器的系数a,b和阶数n计算出来。 分别输入： [n,Wn]=buttord(1/25,10/25,3,60) [a,b]=butter(n,Wn) freqz(a,b,128,1000) 可以得到阶数n=4，系数a={ 0.0003317,0.0009951,0.0009951,0.0003317 },系数b={ 1.0000,-2.7101,2.4608,-0.7480} 因此，差分方程为： y[n]=0.0003317y[n]+ 0.0009951y[n-1]+ 0.0009951y[n-2]+ 0.0003317y[n-3] + x[n]-2.7101x[n-1] +2.4608x[n-2] -0.7480x[n-3] 五．实验步骤 1．设置软件仿真模式并启动CCS： 打开软件，左击出现界面的左边第二行，选择Remove，删除此配置，重新选择配置C55xx Rev3.0 CPU FU…..，点击Add添加，再点击SaveQuit保存退出。 2．打开工程，浏览程序，工程目录为 C:\ICETEK\VC5509Ae\VC5509Ae\Lab0502-IIR\IIR.pjt。 3．打开main.c文件，将前面算出来的阶数n,系数a，b分别带入“#define IIRNUMBER n”、