iir模拟低通滤波器设计.docx

信息科学与工程学院 综合性实验报告 姓 名： 学 号 班 级： 实验项目名称：IIR模拟低通滤波器设计 实验项目性质： 设计性实验 实验所属课程： 数字信号处理 实验室（中心）： 指导教师： 实验完成时间： 年 月 日 教师评阅意见： 签名： 年 月曰 实验成绩： 、实验目的 1、 加深对无限冲激响应IIR滤波器的常用指标和设计过程的理解。 2、 学会用冲激响应不变法把模拟滤波器设计成数字滤波器。 3、 进一步掌握matlab在数字信号处理中的应用，以便以后的学习。 、实验内容及要求 实验内容： 要求按照设计指标设计无限冲激响应 IIR巴特沃什模拟低通滤波器。 实验要求： 必须掌握IIR巴特沃什模拟低通滤波器的各个指标代表的含义，搞清 楚次实验的原理，有可能的话，用冲激响应不变法把模拟滤波器设计 成数字滤波器。 三、实验原理 1、低通滤波器的技术要求用图形表示如下: 和 H(ej )1 0.707- 0.707 - ■ [ Pc s 通带过渡带阻带 通带 过渡带 阻带 如上图表示了一个频域设计、一维低通滤波器的技术要求图。 p和 分别称为通 p s 带截止频率和阻带截止频率。通带频率范围为 ° p ，阻带频率范围 为 s n。从p到s称为过渡带，用 表示，在过渡带里，幅频 特性单调下降。在通带和阻带内的衰减一般用数 dB表示。 通带内允许最大衰减是 p ，阻带内允许最小衰减是 通带内允许最大衰减是 p ，阻带内允许最小衰减是 ，定义分别为: s 4 H(eJ p) s 2°严 H(eJs) 2、Matlab信号处理工具箱中提供了设计巴特沃思模拟滤波器的函数 2、 buttap和butter ，格式如下： .[N ,WC] buttord （WP,WS, RP, RS, s 用于计算巴特沃思模拟低通滤波器的阶 N和3dB截止频率Wc（即本书中的符号 c）。 其中，Wp和Ws分别是滤波器的通带截止频率 p和阻止截止频率 s，单位为rad/s ; Rp和Rs分别是通带最大衰减系数 p和阻带最小衰减系数 s，单位为dB。 一 [z, p,G] buttap（N） 用于计算N阶巴特沃思归一化（c=1）模拟低通滤波器系统函数的零、 极点和增益 因子，返回长度为N的向量z和p分别给出N个零点和极点，G是滤波器增益。得 到的滤波器系统函数形式如下： Ha(s)Qa(S)? (S ZQ(S Z2)L (S Z Ha(s) Pa(s) (s pi)(s P2)L (s Pn) 其中，Zk和Pk分别是向量z和p的第k个元素。如果要从零、极点得到系统函数的 分子和分母多项式系数向量 B和A,可以调用结构转换函数。（实验中没有用到） 亠[B, A] zp2tf （z, p,G），结构转换后系统函数的形式为 Ha(s)B(s)A(s) Ha(s) B(s) A(s) M 1 | ds L N 1 ■ a〔s L 其中，M是向量B的长度，N是向量A的长度，bk和ak分别是向量B和A的第k个元 素。 （3）[B, A] butter （N ,Wc, ftype , s） 用于计算巴特沃什模拟滤波器系统函数中分子和分母多项式系数向量 B和A,其中N 和Wc分别是滤波器的阶和3dB截止频率c,返回向量B和A中的元素ak和bk分别是上 面的Ha(s)表示式中的分母和分子系数。ftype缺省时，设计低通滤波器，S缺省时，设 计数字滤波器。 四、 实验仪器、材料 PC机一台、MATLAB^件实验平台 五、 实验过程及原始记录 1、巴特沃什模拟低通 滤波器设计指标如下： wp=0.4*pi ，ws=0.45*pi， Rp=5, Rs=10; 实验所用到的程序如下： wp=0.4*pi ;ws=0.45*pi ;Rp=5;Rs=10; [N,wc]=buttord(wp,ws,Rp,Rs, s); [B,A]=butter(N,wc, s); w=0:0.001:6; [H,W]=freqs(B ,A, w); H=20*log10(abs(H)); plot(w,H),grid on; xlabel(频率 /(rad/s) );ylabel(幅度 /dB) 经MATLA仿真以后得到的滤波器如下图所示： 巴特沃什模拟低通滤波器  运行后在matlab控制窗口上显示如下结果: N=7 2、为了方便观察故设计了不同阶次 N的滤波器，即Butterworth模拟原型低通滤波器, 其程序代码如下： n=0:0.01:2; for ii=1:4 switch ii case 1,N=2; case 2,N=5; case 3,N=10; case 4,N=20; end [z,p,k]=buttap(N); [b,a]=zp2tf( z,p,k); [H,