实验三无限冲激响应滤波器(IIR)算法实验.docxVIP

实验三 无限冲激响应滤波器(IIR)算法实验一．实验目的1．掌握设计IIR 数字滤波器的原理和方法。2．熟悉IIR 数字滤波器特性。3．了解IIR 数字滤波器的设计方法。二．实验设备PC 兼容机一台操作系统为Windows2000(或Windows98，WindowsXP，以下默认为Windows2000)安装Code Composer Studio 3.1 软件三．实验原理1．无限冲激响应数字滤波器的基础理论。2．模拟滤波器原理（巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器）。3．数字滤波器系数的确定方法。4．根据要求设计低通IIR 滤波器：要求：低通巴特沃斯滤波器在其通带边缘1kHz 处的增益为-3dB，12kHz 处的阻带衰30dB，采样频率25kHz。设计：- 确定待求通带边缘频率fp1Hz、待求阻带边缘频率fs1Hz 和待求阻带衰减-20logδsdB。模拟边缘频率为：fp1=1000Hz，fs1=12000Hz阻带边缘衰减为：-20logδs=30dB-用Ω=2πf/fs 把由Hz 表示的待求边缘频率转换成弧度表示的数字频率，得到Ωp1 和Ωs1。Ωp1=2πfp1/fs=2π1000/25000=0.08π弧度Ωs1=2πfs1/fs=2π12000/25000=0.96π弧度- 计算预扭曲模拟频率以避免双线性变换带来的失真。由w=2fs tan(Ω/2)求得wp1 和ws1，单位为弧度/秒。wp1=2fs tan(Ωp1/2)=6316.5 弧度/秒ws1=2fs tan(Ωs1/2)=794727.2 弧度/秒- 由已给定的阻带衰减-20logδs 确定阻带边缘增益δs。因为-20logδs=30，所以logδs=-30/20，δs=0.03162- 计算所需滤波器的阶数：因此，一阶巴特沃斯滤波器就足以满足要求。- 一阶模拟巴特沃斯滤波器的传输函数为：H(s)=wp1/(s+wp1)=6316.5/(s+6316.5)由双线性变换定义s=2fs(z-1)/(z+1)得到数字滤波器的传输函数为：因此，差分方程为：y[n]=0.7757y[n-1]+0.1122x[n]+0.1122x[n-1]。5．程序流程图：四．实验步骤1．实验准备：- 设置软件仿真模式。- 启动CCS。2．打开工程，浏览程序，工程目录为C:\ICETEK\VC5509Ae\VC5509Ae\Lab0502-IIR\IIR.pjt。3．编译并下载程序。4．打开观察窗口：* 选择菜单View-Graph-Time/Frequency…进行如下图所示设置。* 选择菜单View-Graph-Time/Frequency…，进行如下设置：5．清除显示：在以上打开的窗口中单击鼠标右键，选择弹出式菜单中“Clear Display”功能。6．设置断点：在程序iir.c 中有注释“break point”的语句上设置软件断点。7．运行并观察结果：⑴ 选择“Debug”菜单的“Animate”项，或按Alt+F5 键运行程序。⑵ 观察“IIR”窗口中时域图形；观察滤波效果。8．退出CCS。五．实验结果输入波形为一个低频率的正弦波与一个高频的余弦波叠加而成。如图：输入波形图输入的FFT波形图输出波形图滤波后FFT波形图通过观察频域和时域图，得知：输入波形中的低频波形通过了滤波器，而高频部分则被衰减。六．问题与思考1、试微调(±0.0001)改变程序中 fU 的取值，观察步长因子μ在自适应算法中所起的作用。答、如图可以得知，试微调(±0.0001)改变程序中 fU 的取值，步长因子μ在自适应算法中所起的作用是为自适应滤波提供了更好的收敛特性，在平滑随机扰动和加快快速收敛速度方面起到了一定的作用。