无限冲击响应滤波器算法.docVIP

DSP 实 验 报 告 无 限 冲 击 滤 波 器 算 法 无限冲击响应滤波器(IIR)算法 实验目的 掌握设计IIR数字滤波器的原理和方法； 熟悉IIR数字滤波器的特性； 了解IIR数字示波器的设计方法。 实验设备 PC兼容机一台，操作系统为Windows2000以上，安装Code Coposer Studio 2.0软件 实验原理 无限冲击响应数字滤波器的基础理论； 模拟滤波器原理（巴特沃斯、切比雪夫、椭圆滤波器、贝塞尔滤波器等）。 数字滤波器的确定方法； 根据要求设计低通IIR滤波器 要求：低通巴特沃斯滤波器在其通带边缘1KHz处的增益为-3dB，12KHz处的阻带衰 减为30dB,采样频率为25KHz。设计： —确定待求通带边缘频率fp1Hz、待求阻带边缘频率fs1Hz和待求阻带衰减-20logδsdB。 模拟边缘频率为：fp1 =1000Hz，fs1=12000Hz 阻带衰减为：-20logδdB=30dB —用Ω=2πf/ fs把由HZ表示的待求边缘频率转换成弧度表示，得到Ωp1和Ωs1。 Ωp1=2πfp1/ fs=2π1000/ 25000=0.08πrad Ωs1=2πfs1/ fs=2π12000/ 25000=0.96πrad —计算预扭曲模拟频率以避免双线性变化带来的失真。 由w=2fs tan(Ω/2)求得wp1和ws1，单位是 rad/s Wp1=2fs tan(Ωp1/2)= 6316.5 rad/s Ws1=2fs tan(Ωs1/2)= 794727.2 rad/s —由已给定的阻带衰减-20logδsdB确定阻带边缘增益δs 因为-20logδs=30，所以δs=0.03162 —计算滤波器的阶数 因此，一阶巴特沃斯滤波器就足以满足要求。 —一阶模拟巴特沃斯滤波器的传输函数为： H(s)=wp1/(s+wp1)=6316.5/(s+6316.5) 由双线性变换定义： 因此差分方程：y[n]=0.7757y[n-1]+0.1122x[n]+0.1122x[n-1] 程序流程图 开始 开始 初始化工作变量 调用波形发生子程序产生混叠的波形(高频+低频) 调用IIR滤波子程序计算当前输出 波形发生 计算步长 用标准C的sin函数和cos函数计算当前波形值 返回波形值 IIR滤波 用滤波器系数乘以保存的N-1个输入输出值和当前输入值并求和 返回计算结果 实验步骤 实验准备 —设置软件仿真模式 —启动CCS 打开工程，浏览程序，工程目录为D：\dsp\t7\iir\iir.pjt 编译并下载程序 打开窗口：选择菜单 View-Graph-Time/Frquency…，进行如下设置 5.清楚显示：在以上的窗口单机鼠标右键，选择弹出菜单中“Clear Display” 6.设置断点：在程序iir.c中加入软件断点 7.运行观察结果 (1)选择“Debug”菜单中的“RUN”项运行程序 ?(2)观察“IIR”窗口中时域波形；观察滤波效果。 实验结果 输出的波形为一个低频率的正弦波与一个高频的余弦波叠加而成，如图： Matlab仿真IIR滤波器的单位脉冲响应 源程序 B=[0.1122 0.1122]; A=[1 -0.7757]; 源程序 B=[0.1122 0.1122]; A=[1 -0.7757]; [H w]=freqz(B,A); plot(w/2/pi,abs(H)) 绘制显示出中幅频特性曲线，如图所示： 通过观察时域和频域图，得知：输入波形中的低频波形通过了滤波器，而高频部分则被衰减。Matlab的仿真也可以看出来，IIR滤波器在高频处的衰减使得信号中的高频分量被衰减，从而实现低通滤波器的功能。 实验分析： 输入波形为一个低频率的正弦波与一个高频的正弦波叠加而成。通过观察频域和时域图，得知：输入波形中的低频波形通过了滤波器，而高频部分则大部分被滤除。通过对滤波器参数的设计，可以改变该滤波器的滤波性能。 实验心得：通过这个实验把以前学过的数字信号处理的知识有了更深的巩固了了解，尤其是IIR滤波器的设计和应用方面，从模拟滤波器的设计，到模拟数字角频率的转换，再到最后数字滤波器的滤波有了很好的认识，这个对我们加深数字信号处理的认识有很好的帮助，更加直观的感受到了信号处理的作用和意义，对我们以后的学习有很大的帮助。 程序代码 #include DSP281x_Device.h // DSP281x Headerfile Include File #include DSP281x_Examples.h // DSP281x Examples Include Fi