#课程设计课程设计名称：基于脉冲响应不变法设计切比雪夫II型IIR数字低通滤波器 专业班级：学生姓名 ：墨蓝的星空学号：指导教师：课程设计时间:2013年6月#数字信号处理专业课程设计任务书学生姓名星空专业班级学号题目基于脉冲响应不变法设计Chebyshev-Il型IIR数字低通滤波器课题性质工程技术研究课题来源自拟课题指导教师同组姓名主要内容设计一个具有如下性能指标的 Elliptic 型IIR数字低通滤波器：%=0.25江，s=0.4兀，5p=0.01，6s=0.001 ;分析最后设计结果性 能。任务要求1. 掌握脉冲响应不变法设计数字滤波器的基本流程；2. 掌握Elliptic型模拟低通原型滤波器；3. 求出所设计滤波器的传输函数；4. 绘制所设计滤波器的幅频特性图。参考文献1 .胡广书数字信号处理理论、算法与实现M.北京：清华大学出版社，1997.2 .R. Lyons. Understanding Digital Signal Processing M. 2nd ed.Pren tice Hall3 . A.V.奥本海姆，R.W.谢弗and J.R.巴克.离散时间信号处理M.第二版.西安交通大学出版社，2001.PTR., 2004.4 . S. K. Mitra. Digital Signal Processing: A Computer-BasedApproachM. 3rd ed. McGraw-Hill, 2005.审查意见指导教师签字：教研室主任签字：年月日说明：本表由指导教师填写，由教研室主任审核后下达给选题学生，装订在设计（论文）首页1需求分析切比雪夫数字(Chbyshev)滤波器的振幅特性具有等波纹特性，低通滤波器 振幅特性在通带内是单调的、在阻带内是等波纹的。阻带内具有等波纹的欺负特 性，而在通带内是单调的、平滑的，阶数越高，频率特性曲线越接近矩形，传输 函数既有极点又有零点。本设计要求切比雪夫II型的数字滤波器所需的四个参数分别是归一化的通带截止频率 =0.25二，阻带截止频率=0.4二，通带误差 容限； p=0.01，通带误差容限；s=0.001 ;由此得到对应的模拟原型低通滤波器的 各个主要参数为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数字滤波器，它能很好地重现的原型模拟滤波器频率特性。基本实现流程如下图所示确定数字低通滤波器的各项性能指标j把器的技术指标转换成模拟低通滤波器的技术指标4 丄IF円按照模拟低通滤波器的技术指标设计模拟低通滤波器L. JVrr利用双线性变换法将模拟滤波器转变为数字低通滤波器利用模拟切比雪夫滤波器设计数字滤波器- 图2.1Chebyshev-ll 型IIR数字低通滤波器设计流程图3运行环境操作系统： Windows 7软件：MATLAB4开发工具和编程语言MATLAB口 MATLAffi程语言5详细设计(完整代码见最后一页)数字滤波器采用经典低通滤波器作为连续域上的设计模型， 通过频域变换得 到IIR数字滤波器，最后还要进行离散化处理。用 MATLA提供的低通模拟滤波 器原型函数cheb2 ap频域变换函数包括Ip2lp ;离散化处理函数imp in var。(1) 确定数字低通滤波器的技术指标：通带截止频率、阻带截止频率、通带 最小衰减和阻带最小衰减。p =0.25 二，,=0.4 二，、p =0.01， s =0.001 ；(2) 将数字滤波器的技术指标转换成模拟低通滤波器的技术指标(Ts=0.01)。Wp=0.25*pi/Ts;%计算模拟通带截止频率Ws=0.4*pi/Ts;%计算模拟阻带截止频率Rp=20*log10(1/0.99);%设置峰值通带波纹Rs=20*log10(1/0.001);%设置最小阻带衰减#(3) 按照模拟低通滤波器的技术指标设计模拟低通滤波器N,Wn=cheb2ord(Wp,Ws,Rp,Rs,s); %确定模拟滤波器传输函数最低阶数 Nz,p,k =cheb2ap(N,Rs);%设计切比雪夫II型模拟滤波器b,a=zp2tf(z,p,k);%将零点极点增益形式转换为传输函数形式(4) 按频率变换设计模拟滤波器bt,at=lp2lp(b,a,Wn); %实现低通模拟原型滤波器至低通滤波器的频率转换(5) 用脉冲响应不变法模拟离散化，+将模拟滤波器传输函数转换为数字滤波 器传输函数。nu m,de n=impi nv ar(bt,at,fs);%脉冲响应不变法得到数字滤波器的传输函数6调试分析第一次设计滤波器时直接由N,Wn=cheb2ord(Wp,Ws,Rp,Rs,s);再利用cheby2(N,Rs,Wn) 生成数字低通滤波器，而且绘出的幅频相应也满足要求。但是，经过检查发现这种方法是完全设计法，就是利用MATLAB设计工具箱中提供的完全工具函数直接生成数字低通滤波器，并不符合课程设计的要求， 也没有用到脉冲响应不变法。然后我查阅了 MATLAB设计低通滤波器的相关书籍 后，列出了低通滤波器的设计流程图，才使自己设计过程清晰明了起来。7测试结果h,W=freqs(b,a);H,w=freqz( nu m,de n);%绘出切比雪夫II型模拟低通滤波器的幅频特性曲 plot(W/pi,abs(h),r);线切比雪夭It型模拟低通滤波器幅频响应#(2)plot(w/pi,abs(H),r); %绘出切比雪夫II型数字低通滤波器的幅频特性曲线#(3)plot(w/pi,angle(H),r); %绘出切比雪夫II型数字低通滤波器的相频特性曲线(4)zpla ne(n um,de n);% 绘出零极图零极图#参考文献1 万永革，数字信号处理的 MATLABS现，中国林业出版社，20066252 陈亚勇，MATLA信号处理详解.人民邮电出版社，2001.43 周辉，数字信号处理的MATLA实现，科学出版社，2006.24 程佩青，数字信号处理教程， 清华大学出版社： 2012.06.5 胡广书，数字信号处理，清华大学出版社， 2012.102007刘泉 阙大顺 , 数字信号处理原理与实现 . 北京：电子工业出版社，#心得体会本次数字信号处理课程设计过程中，我掌握了设计数字滤波器的思路和方 法，了解了模拟滤波器的特性， 同时锻炼了自己把课本上的理论知识运用到时间 中的能力培养了科学严谨的态度。 当然，在设计的过程中， 不可避免的遇到了很 多问题，如刚开始思路比较混乱，没有明确的方向。经过参考书籍和网上查阅， 总结出了模拟原型法的基本步骤。 设计主要是如何将理论计算的模型转换为仿真 模型。因为在理论上， 将低通转换成高通， 一般是变换将高通频率特征转换成低 通原型频率特征； 而在软件设计中， 是对其传递函数进行修改， 即变换其 z 域的 表达式，设计初期一直不知道如何将其联系起来。 后来发现，其实变换传递函数， 也就是变换频率特征， 是将变换后的频率代入原低通模型， 而后得到高通模型的。 数字滤波器的设计流程有了大致的了解， 掌握了一些设计滤波器的基本方法， 提 高了将理论用于实际， 并通过软件实现相关的函数功能能力， 最终设计出按自己 要求的滤波器，掌握了更多专业相关的使用知识与技能。其次课程设计使我进一步熟悉了计算工具软件 -MATLAB,并进一步掌握 了 MATLAB勺使用方法.对MATLA语言的发展和特点有了更深的了解，熟悉其工 作环境，并掌握了 MATLAB勺基本语法，在其应用方面也有了更深的了解，了解 了 MATLAB?序书写的一般思路及一般步骤，学会了用MATLA解决复杂信号处理 等问题的方法，使我对 MATLA信号处理工具箱内的模型函数有了深刻的体会和 应用。但更为重要的是， 对于解决一个问题的思路更加清晰， 我找到了属于自己 的方法，那就是要善于合理利用资源，积极查阅查阅相关书籍，独立自主，耐心 细致，虚心向别人请教。同时，我也发现了自己在动手实践中很多的不足，在以 后的学习中， 我一定会扎实学好基本技能， 努力实践，不断提高自己的综合能力。#脉冲响应不变法设计切比雪夫 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,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.1);grid on;subplot(2,2,3);plot(w/pi,angle(H),r);%绘出切比雪夫 II 型数字低通滤波器的相频特性曲线xlabel(频率八pi);ylabel(相位);title( 切比雪夫 II 型数