《数字信号处理》实验5 IIR数字滤波器设计.pptxVIP

《数字信号处理实验与课程设计教程--面向工程教育》 实验目的：本实验包含4个实验子项目，通过matlab编程仿真，掌握IIR数字滤波器设计方法，包括：1.掌握模拟滤波器设计，能利用频带变换实现从低通到高通/带通/带阻滤波器的转换。2.掌握用脉冲响应不变法设计IIR数字滤波器。3.掌握用双线性变换法设计IIR数字滤波器。4.熟悉IIR数字滤波器的应用实例。实验五 IIR数字滤波器设计 实验五 IIR数字滤波器设计5.1 模拟滤波器设计5.2 脉冲响应不变法/双线性变换法设计IIR数字滤波器5.3 IIR数字滤波器的应用实例 5.1 模拟滤波器设计一、实验原理1.数字滤波器及其分类 选频型数字滤波器的功能是对输入信号进行数值运算，保留信号中的有用频率成分而无用频率成分，即：能实现选频功能的数字硬件系统或程序模块。数字滤波器通常分为：无限长单位单位脉冲响应(IIR)滤波器与有限长单位单位脉冲响应(FIR)滤波器，其中：IIR数字滤波器可借助模拟滤波器设计方法：先根据技术指标设计模拟滤波器，再将其转换成数字滤波器。模拟滤波器设计理论非常成熟，且有多种性能优良的滤波器可供选择，如巴特沃斯、切比雪夫和椭圆滤波器等，设计公式和图表十分完善。2.模拟滤波器的技术指标 模拟滤波器的技术指标有：通带截止频率Ωp,阻带截止频率Ωs,通带最大衰减αp，阻带最小衰减αs。3.模拟滤波器设计步骤(1)根据信号处理的要求确定设计指标并选定滤波器类型。(2)设计模拟原型低通滤波器：计算滤波器阶数N和3dB截止频率Ωc；通过查表计算滤波器系统函数Ha(s)。(3)通过频带变换将模拟原型低通滤波器转换为所需类型的滤波器H(s)。 4.本实验相关Matlab函数 (1)由技术指标求模拟滤波器阶数N和3dB截止频率Ωc。1）巴特沃斯低通滤波器MATLAB函数：[N,wc]=buttord(wp,ws,ap,as,’s’)其中：N---滤波器阶数；wc—模拟滤波器3dB截止频率(rad/s)；wp--通带截止频率Ωp；ws---阻带截止频率Ωs；ap----通带最大衰减αp；as---阻带最小衰减αs。2）切比雪夫I型低通滤波器MATLAB函数：[N,wpo]=cheb1ord(wp,ws,ap,as,’s’)其中：N---滤波器阶数；wpo—模拟滤波器通带截止频率(rad/s)，其他参数含义与1)相同。3）切比雪夫II型低通滤波器MATLAB函数：[N,wso]=cheb2ord(wp,ws,ap,as,’s’)其中：N---滤波器阶数；wso—模拟滤波器阻带截止频率(rad/s)，其他参数含义与1)相同。4）椭圆低通滤波器MATLAB函数：[N,wpo]=ellipord(wp,ws,ap,as,’s’)其中：N---滤波器阶数；wpo—模拟滤波器通带截止频率(rad/s)，其他参数含义与1)相同。 (2) 利用频带变换实现从低通到高通/带通/带阻滤波器的转换1）巴特沃斯滤波器MATLAB函数：[B,A]=butter(N,wc,’ftype’,’s’)其中：其中：N---滤波器阶数；wc—模拟滤波器3dB截止频率(rad/s)；ftype---滤波器类型，ftype=high----高通滤波器；ftype=stop----带阻滤波器；缺省为低通滤波器；当w=[wc1 wc2]时，为带通滤波器。B/A---模拟滤波器系统函数H(s)的分子/分母多项式系数。H(s)= 2）切比雪夫I型和II型滤波器I型：MATLAB函数：[B,A]=cheby1(N,ap,wpo,’ftype’,’s’)II型：MATLAB函数：[B,A]=cheby2(N,as,wso,’ftype’,’s’)3)椭圆滤波器MATLAB函数：[B,A]=ellip(N,ap,as,wpo,’ftype’,’s’)(3)模拟滤波器的频响特性MATLAB函数：[H,Omeg]=freqs(B,A)其中：H----模拟滤波器H(s)的频率响应；B,A----分别是H(s)分子/分母多项式系数; Omeg ---模拟滤波器的角频率(rad/s)。二、实验环境1.计算机1台。2. Windows 7 以上操作系统。3. Matlab 7.0 以上版本软件。 三、实验参考与实验内容1. 模拟低通滤波器设计 (1)巴特沃斯滤波器设计 设计一个巴特沃斯模拟低通滤波器，通带截止频率fp=1kHz,阻带截止频率fs=5kHz,通带 最大衰减αp=1dB,阻带最小衰减αs=40dB。 实验参考程序：参见: 实验项目参考源码\ch5\ch5prog1.m实验内容：(1