脉冲响应不变法设计切比雪夫II型IIR数字低通滤波器汇总

1、河由7才比殍课程设计课程设计名称：基于脉冲响应不变法设计切比雪夫II型IIR数字低通滤波器专业班级：学生姓名： 墨蓝的星空学 号：指导教师：课程设计时间： 2013年6月ii数字信号处理 专业课程设计任务书学生姓名星空专业班级学号题目基于脉冲响应不义法设计Chebyshev-II型IIR数字低通滤波器课题性质工程技术研究课题来源自拟课题指导教师同组姓名主要内容设一个具有如下性能指标的曰liptic 型IIR数字低通滤波器： 3=0.25兀，初=0.4n, 5p =0.01 , 6s=0.001 ;分析最后设计结果性 能。任务要求1 .掌握脉冲响应不变法设计数字滤波器的基本流程；2 .掌握Ell

2、iptic型模拟低通原型滤波器；3 .求出所设计滤波器的传输函数；4 .绘制所设计滤波器的幅频特性图。奔f文献1 .胡广书.数字恬号处理理论、算法与实现M.北京：清华大学出版社，1997.2 . R. Lyons. Understanding Digital Signal Processing M. 2nd ed. Prentice Hall3 . A.V.奥本海姆，R.W.谢弗and J.R.巴克.离散时间佶号处理M.第二 版.西安交通大学出版社，2001.PTR., 2004.4 . S. K. Mitra. Digital Signal Processing: A Computer-Ba

3、sed ApproachM. 3rd ed. McGraw-Hill, 2005.审查意见指导教师签字：教研室主任签字：年月日说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页1需求分析切比雪夫数字(Chbyshev)滤波器的振幅特性具有等波纹特性，低通滤波器振幅特性在通带内是单调的、在阻带内是等波纹的。阻带内具有等波纹的欺负特性，而在通带内是单调的、平滑的，阶数越高，频率特性曲线越接近矩形，传输函数既有极点又有零点。本设计要求切比雪夫II型的数字滤波器所需的四个参数分别是归一化的通带截止频率 6 P =0.25 n ,阻带截止频率6 s =0.4 n ,通带误

4、差 容限6p =0.01,通带误差容限6s=0.001;由此得到对应的模拟原型低通滤波器的 p各个主要参数为Wp=0.25*pi/Ts;Ws=0.4*pi/Ts;Rp=20*log10(1/0.99);Rs=20*log10(1/0.001);2概要设计本设计采用经典设计法设计IIR数字低通滤波器，就是先根据技术指标设 计出来相应的模拟滤波器，然后把设计好的模拟滤波器通过脉冲响应不变法转换 成IIR数字滤波器，它能很好地重现的原型模拟滤波器频率特性。基本实现流程 如下图所示图2.1Chebyshev-II 型IIR数字低通滤波器设计流程图3运行环境操作系统： Windows 7 软件：MATL

5、AB4开发工具和编程语言MATLAB口 MATLA褊程语言5详细设计(完整代码见最后一页)数字滤波器采用经典低通滤波器作为连续域上的设计模型， 通过频域变换得 到IIR数字滤波器，最后还要进行离散化处理。用 MATLA提供的低通模拟滤波 器原型函数cheb2ap频域变换函数包括lp21P ；离散化处理函数impinvar。(1)确定数字低通滤波器的技术指标：通带截止频率、阻带截止频率、通带 最小衰减和阻带最小衰减。0P =0.25 n , 0s =0.4 n, 6P =0.01 , 1=0.001;(2)将数字滤波器的技术指标转换成模拟低通滤波器的技术指标(Ts=0.01)。Wp=0.25*p

6、i/Ts;%计算模拟通带截止频率Ws=0.4*pi/Ts;%计算模拟阻带截止频率Rp=20*log10(1/0.99);%设置峰值通带波纹Rs=20*log10(1/0.001);%设置最小阻带衰减(3)按照模拟低通滤波器的技术指标设计模拟低通滤波器N,Wn=cheb20rd(Wp,Ws,Rp,Rs,s); %确定模拟滤波器传输函数最低阶数 Nz,p,k =cheb2ap(N,Rs);%设计切比雪夫II型模拟滤波器b,a=zp2tf(z,p,k);%将零点极点增益形式转换为传输函数形式(4)按频率变换设计模拟滤波器bt,at=lp21P(b,a,Wn); %实现低通模拟原型滤波器至低通滤波器的

7、频率转换(5)用脉冲响应不变法模拟离散化，+将模拟滤波器传输函数转换为数字滤波 器传输函数。num,den=impinvar(bt,at,fs);%脉冲响应不变法得到数字滤波器的传输函数6调试分析第一次设计滤波器时直接由N,Wn=cheb20rd(Wp,Ws,Rp,Rs,s);再利用cheby2(N,Rs,Wn)生成数字低通滤波器，而且绘出的幅频相应也满足要求。但是，经过检查发现这种方法是完全设计法，就是利用MATLA暇计工具箱中提供的完全工具函数直接生成数字低通滤波器，并不符合课程设计的要求， 也没有用到脉冲响应不变法。然后我查阅了 MATLAB计低通滤波器的相关书籍 后，列出了低通滤波器的

8、设计流程图，才使自己设计过程清晰明了起来。7测试结果h,W=freqs(b,a);H,w=freqz(num,den);(1)plot(W/pi,abs(h),r);%绘出切比雪夫II型模拟低通滤波器的幅频特性曲10!IIljr11 . aJrr -:：i ： ：1!：ii1-111/1切比雪夫理模拟低通滤波器幅频响应0.2蚂0 60.40.30 90 00.10.20.30,40 50.60 7频率以(2)plot(w/pi,abs(H),r); %绘出切比雪夫II型数字低通滤波器的幅频特性曲线(3)plot(w/pi,angle(H),r); %绘出切比雪夫II型数字低通滤波器的相频特性曲

9、线切比雪天理数字低通滤波器幅频响应(粤频率改(4)zplane(num,den);% 绘出零极图零极图参考文献1万永革，数字信号处理的 MATLA殴现，中国林业出版社，2006.6.252陈亚勇，MATLABF号处理详解.人民邮电出版社，2001.43周辉，数字信号处理的MATLA映现，科学出版社，2006.24 程佩青，数字信号处理教程， 清华大学出版社： 2012.06.5 胡广书，数字信号处理，清华大学出版社， 2012.1020076 刘泉 阙大顺 , 数字信号处理原理与实现. 北京：电子工业出版社，心得体会本次数字信号处理课程设计过程中，我掌握了设计数字滤波器的思路和方法， 了解了模

10、拟滤波器的特性， 同时锻炼了自己把课本上的理论知识运用到时间中的能力培养了科学严谨的态度。 当然， 在设计的过程中， 不可避免的遇到了很多问题，如刚开始思路比较混乱，没有明确的方向。经过参考书籍和网上查阅，总结出了模拟原型法的基本步骤。 设计主要是如何将理论计算的模型转换为仿真模型。 因为在理论上， 将低通转换成高通， 一般是变换将高通频率特征转换成低通原型频率特征； 而在软件设计中， 是对其传递函数进行修改， 即变换其 z 域的 表达式， 设计初期一直不知道如何将其联系起来。 后来发现， 其实变换传递函数， 也就是变换频率特征， 是将变换后的频率代入原低通模型， 而后得到高通模型的。 数字滤

11、波器的设计流程有了大致的了解， 掌握了一些设计滤波器的基本方法， 提 高了将理论用于实际， 并通过软件实现相关的函数功能能力， 最终设计出按自己 要求的滤波器，掌握了更多专业相关的使用知识与技能。其次课程设计使我进一步熟悉了计算工具软件-MATLAB,并进一步掌握了 MATLAB勺使用方法.对MATLAB言的发展和特点有了更深的了解，熟悉其工 作环境，并掌握了 MATLAB勺基本语法，在其应用方面也有了更深的了解，了解 了 MATLA毓序书写的一般思路及一般步骤，学会了用MATLAB单决复杂信号处理 等问题的方法，使我对 MATLA蛤号处理工具箱内的模型函数有了深刻的体会和 应用。 但更为重要

12、的是， 对于解决一个问题的思路更加清晰， 我找到了属于自己 的方法，那就是要善于合理利用资源，积极查阅查阅相关书籍，独立自主，耐心 细致，虚心向别人请教。同时，我也发现了自己在动手实践中很多的不足，在以 后的学习中， 我一定会扎实学好基本技能， 努力实践， 不断提高自己的综合能力。脉冲响应不变法设计切比雪夫II 型 IIR 数字低通滤波器完整代Ts=0.01;% 计算模拟通带截止频率% 计算模拟阻带截止频率%设置峰值通带波纹%设置最小阻带衰减fs=1/Ts;Wp=0.25*pi/Ts;Ws=0.4*pi/Ts;Rp=20*log10(1/0.99);Rs=20*log10(1/0.001);N

13、,Wn=cheb2ord(Wp,Ws,Rp,Rs,s); %确定模拟滤波器传输函数最低阶数N 和 3dB 截止频率z,p,k =cheb2ap(N,Rs);b,a=zp2tf(z,p,k);bt,at=lp2lp(b,a,Wn);h,W=freqs(b,a);num,den=impinvar(bt,at,fs);H,w=freqz(num,den);subplot(2,2, 1);plot(W/pi,abs(h),r);线%设计切比雪夫II 型滤波器%将零点极点增益形式转换为传输函数形式% 调用脉冲响应不变法得到数字滤波器的传输函数% 绘出切比雪夫II 型模拟低通滤波器的幅频特性曲xlabel( 频率/pi);ylabel( 幅度);title( 切比雪夫 II 型模拟低通滤波器幅频响应);axis(0,1,0,1.1);grid on;subplot(2,2, 2);plot(w/pi,abs(H),r);%绘出切比雪夫II 型数字低通滤波器的幅频特性曲线xlabel(频率八pi);ylabel( 幅度 );title( 切比雪夫 II 型数字低通滤波器幅频响应);axis(0,1,0,1.