基于MATLAB的IIR数字滤波器的设计说明

1、 PAGE43 / NUMPAGES48课程设计说明书题目：基于MATLAB的IIR数字滤波器的设计姓 名：院 （系）：专业班级：学 号：指导教师：成 绩：时间：年月日至 年月日课 程 设 计 任 务 书题目基于MATLAB的IIR数字滤波器的设计专业、班级学号主要容、基本要求、主要参考资料等：主要容：利用四种模拟原型滤波器（巴特沃斯、切比雪夫I型、切比雪夫II型、椭圆型）和两种模/数转换方法（脉冲响应不变法、双线性变换法）分别进行IIR数字滤波器的设计。基本要求：根据给定的各类滤波器的技术指标，分别设计实现数字高通滤波器、数字带通滤波器和数字带阻滤波器，并据此进行分析总结：1、在一样的技术指

2、标要求下，用不同的模拟原型滤波器实现有何异同。2、在一样的技术指标要求下，用不同的模/数转换方法实现有何异同。主要参考资料：1、数字信号处理教程（第三版），程佩青著，清华大学，2007。2、数字信号处理教程MATLAB释义与实现（第2版），怀琛著，电子工业，2008。完 成 期 限：指导教师签名：课程负责人签名：年 月日基于MATLAB的IIR数字滤波器的设计摘 要利用MATLAB 设计滤波器,可以按照设计要求非常方便地调整设计参数,极减轻了设计的工作量,有利于滤波器设计的最优化。Matlab因其强大的数据处理功能被广泛应用于工程计算,其丰富的工具箱为工程计算提供了便利,利用Matlab信号处

3、理工具箱可以快速有效地设计各种数字滤波器,设计简单方便。本文介绍了在MATLAB R2009a 环境下滤波器设计的方法和步骤。关键词MATLAB IIR数字滤波器模拟滤波器目 录 TOC o 1-3 h z u HYPERLINK l _Toc360313998摘要IHYPERLINK l _Toc3603139991数字滤波器 PAGEREF _Toc360313999 h 1HYPERLINK l _Toc3603140001.1数字滤波器的概念 PAGEREF _Toc360314000 h 1HYPERLINK l _Toc3603140011.2数字滤波器的分类 PAGEREF _T

4、oc360314001 h 1HYPERLINK l _Toc3603140021.3数字滤波器的设计要求 PAGEREF _Toc360314002 h 3HYPERLINK l _Toc3603140032 IIR数字滤波器的设计 PAGEREF _Toc360314003 h 4HYPERLINK l _Toc3603140042.1 IIR数字滤波器的设计步骤 PAGEREF _Toc360314004 h 4HYPERLINK l _Toc3603140052.2 用脉冲响应不变法设计IIR数字滤波器 PAGEREF _Toc360314005 h 5HYPERLINK l _Toc

5、3603140062.3 双线性变换法设计IIR数字滤波器 PAGEREF _Toc360314006 h 7HYPERLINK l _Toc3603140073 IIR滤波器的MATLAB设计 PAGEREF _Toc360314007 h 10HYPERLINK l _Toc3603140083.1巴特沃斯数字滤波器的设计 PAGEREF _Toc360314008 h 11HYPERLINK l _Toc3603140093.1.1数字高通滤波器的设计 PAGEREF _Toc360314009 h 12HYPERLINK l _Toc3603140103.1.2 数字带通滤波器的设计

6、PAGEREF _Toc360314010 h 14HYPERLINK l _Toc3603140113.1.3数字带阻滤波器的设计 PAGEREF _Toc360314011 h 17HYPERLINK l _Toc3603140123.2切比雪夫数字滤波器的设计 PAGEREF _Toc360314012 h 19HYPERLINK l _Toc3603140133.2.1数字高通滤波器的设计 PAGEREF _Toc360314013 h 20HYPERLINK l _Toc3603140143.2.2数字带通滤波器的设计 PAGEREF _Toc360314014 h 24HYPERL

7、INK l _Toc3603140153.2.3数字带阻滤波器的设计 PAGEREF _Toc360314015 h 29HYPERLINK l _Toc3603140163.3 椭圆形滤波器的设计 PAGEREF _Toc360314016 h 33HYPERLINK l _Toc3603140173.3.1数字高通滤波器的设计 PAGEREF _Toc360314017 h 34HYPERLINK l _Toc3603140183.3.2数字带通滤波器的设计 PAGEREF _Toc360314018 h 36HYPERLINK l _Toc3603140193.3.3数字带阻滤波器的设计

8、 PAGEREF _Toc360314019 h 38HYPERLINK l _Toc360314020总结 PAGEREF _Toc360314020 h 41HYPERLINK l _Toc360314021参考文献 PAGEREF _Toc360314021 h 421数字滤波器1.1数字滤波器的概念滤波器是指用来对输入信号进行滤波的硬件和软件。数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波实质上是一种运算过程，实现对信号的运算处理。数

9、字滤波器和模拟滤波器相比，因为信号的形式和实现滤波的方法不同，数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性:（式1-1）其中、分别是数字滤波器的输出序列和输入序列的频域特性（或称为频谱特性）,是数字滤波器的单位取样响应的频谱，又称为数字滤波器的频域响应。输入序列的频谱经过滤波后,因

10、此，只要按照输入信号频谱的特点和处理信号的目的，适当选择，使得滤波后的满足设计的要求，这就是数字滤波器的滤波原理。1.2 数字滤波器的分类按照不同的分类方法，数字滤波器有许多种类，但总起来可以分成两大类：经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的频率成分和希望滤除的频率成分占有不同的频带，通过一个合适的选频滤波器滤除干扰，得到纯净信号，达到滤波的目的。但是，如果信号和干扰的频谱相互重叠，则经典滤波器不能有效地滤除干扰，最大限度地恢复信号，这时就需要现代滤波器，例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。现代滤波器是根据随机信号的一些统计特性，在某种最佳准则下，最大

11、限度地抑制干扰，同时最大限度地回复信号，从而达到最佳滤波的目的。经典数字滤波器从滤波特性上分类，可以分为：低通滤波器、高通滤波器、带通滤波器、带阻滤波器。图1-1 各种理想滤波器的幅频特性数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。IIR 数字滤波器的特征是，具有无限持续时间冲激响应，需要用递归模型来实现，其差分方程为：（式1-2）系统函数为： （式1-3）设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z)，使其频率响应H(z)满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止频率、通带衰

12、减系数和阻带衰减系数。1.3数字滤波器的设计要求滤波器的指标常常在频域给出。数字滤波器的频响特性函数一般为复函数，所以通常表示为：（1-4）其中，|称为幅频特性函数，(w)称为相频特性函数。幅频特性表示信号通过该滤波器后各频率成分的衰减情况，而相频特性反映各频率通过滤波器后在时间上的延时情况。一般IIR数字滤波器，通常只用幅频响应函数|来描述设计指标，相频特性一般不作要求。IIR滤波器指标参数如下图所示。图中，p和s分别为通带边界频率和阻带边界频率；1和2分别为通带波纹和阻带波纹；允许的衰减一般用dB数表示，通带所允许的最大衰减（dB）和阻带允许的最小衰减（dB）分别为p和s表示：（式1-5）

13、（式1-6）一般要求：当时，；当时，。图1-2 低通滤波器的技术要求2IIR数字滤波器的设计2.1IIR数字滤波器的设计步骤IIR数字滤波器的设计一般有两种方法：一个是借助模拟滤波器的设计方法进行。其设计步骤是，先设计模拟滤波器，再按照某种方法转换成数字滤波器。这种方法比较容易一些，因为模拟滤波器的设计方法已经非常成熟，不仅有完整的设计公式，还有完善的图表供查阅；另外一种直接在频率或者时域进行，由于需要解联立方程，设计时需要计算机做辅助设计。其设计步骤是：先设计过渡模拟滤波器得到系统函数，然后将按某种方法转换成数字滤波器的系统函数。这是因为模拟滤波器的设计方法已经很成熟，不仅有完整设计公式，还

14、有完善的图表和曲线供查阅；另外，还有一些典型的优良滤波器类型可供我们使用。为了保证转换后的稳定且满足技术指标要求，对转换关系提出两点要求:因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。数字滤波器的频率相应模仿模拟滤波器的频响特性，s平面的虚轴映射为z平面的单位圆，相应的频率之间呈线性关系。利用模拟滤波器成熟的理论设计IIR数字滤波器的过程是：(1)确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减、阻带截止频率、阻带最小衰减。(2)将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。(3)按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。(4)用所选的转换方法，将模

15、拟滤波器转换成数字低通滤波器系统函数。IIR数字滤波器的设计流程图如下：数字滤波器技术指标 模拟滤波器 技术指标 数字滤波器 模拟滤波器 模拟滤波器 设计方法变换变换图2-1 IIR数字滤波器的设计步骤流程图成熟的模拟滤波器设计方法主要有脉冲响应不变法和双线性变换法。2.2用脉冲响应不变法设计IIR数字滤波器一、设计原理利用模拟滤波器来设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可以从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应ha(t)，即将ha(t)进行等间隔采样，使h(n)正好等于ha(t)的

16、采样值，满足h(n)=ha(nT) 式中,T是采样周期。如果令Ha(s)是ha(t)的拉普拉斯变换，H(z)为h(n)的Z变换，利用采样序列的Z变换与模拟信号的拉普拉斯变换的关系得(式2-1)则可看出，脉冲响应不变法将模拟滤波器的S平面变换成数字滤波器的Z平面，图2-2 脉冲响应不变法的映射关系这个从s到z的变换z=esT是从S平面变换到Z平面的标准变换关系式。由（2-1）式，数字滤波器的频率响应和模拟滤波器的频率响应间的关系为(式2-2)这就是说，数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。正如采样定理所讨论的，只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以时，即(式2-3)

17、才能使数字滤波器的频率响应在折叠频率以重现模拟滤波器的频率响应，而不产生混叠失真，即 (式2-4)但是，任何一个实际的模拟滤波器频率响应都不是严格限带的，变换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。图2-3脉冲响应不变法中的频响混叠现象对某一模拟滤波器的单位冲激响应ha(t)进行采样，采样频率为fs，若使fs增加，即令采样时间间隔（T=1/fs）减小，则系统频率响应

18、各周期延拓分量之间相距更远，因而可减小频率响应的混叠效应。二、脉冲响应不变法优缺点从以上讨论可以看出，脉冲响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好，而且模拟频率和数字频率之间呈线性关系=T。因而，一个线性相位的模拟滤波器（例如贝塞尔滤波器）通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器。脉冲响应不变法的最大缺点是有频率响应的混叠效应。所以，脉冲响应不变法只适用于限带的模拟滤波器(例如，衰减特性很好的低通或带通滤波器)，而且高频衰减越快，混叠效应越小。至于高通和带阻滤波器，由于它们在高频部分不衰减，因此将完全混淆在低频响应中。如果要对高通

19、和带阻滤波器采用脉冲响应不变法，就必须先对高通和带阻滤波器加一保护滤波器，滤掉高于折叠频率以上的频率，然后再使用脉冲响应不变法转换为数字滤波器。当然这样会进一步增加设计复杂性和滤波器的阶数。2.3 双线性变换法设计IIR数字滤波器脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率围压缩到-/T/T之间，再用z=esT转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-/T/T一条横带里；第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。这样就使S平面

20、与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图2-4所示。图2-4双线性变换的映射关系为了将S平面的整个虚轴j压缩到S1平面j1轴上的-/T到/T段上，可以通过以下的正切变换实现（式2-5）式中,T仍是采样间隔。当1由-/T经过0变化到/T时，由-经过0变化到+，也即映射了整个j轴。将式（2-5）写成（式2-6）将此关系解析延拓到整个S平面和S1平面，令j=s，j1=s1，则得 （式2-7）再将S1平面通过以下标准变换关系映射到Z平面 （式2-8）从而得到S平面和Z平面的单值映射关系为：（式2-9）（式2-10）式（2-9）与式（2-10）是S平面与Z

21、平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换式（2-5）与式（2-9）的双线性变换符合映射变换应满足的两点要求。首先,把z=ej，可得（式2-11）即S平面的虚轴映射到Z平面的单位圆。其次，将s=+j代入式（2-11），得（式2-12）因此 （式2-13）由此看出，当0时，|z|0时，|z|1。也就是说，S平面的左半平面映射到Z平面的单位圆，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。二、双线性变换法优缺点双线性变换法与脉冲响应不变法相比，其主要的优点是避免了频率响应的混叠现象。这是因为S平面与Z平面是单值的一一对应关系。S平面整个j

22、轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。这个关系如式（2-9）所示，重写如下：（式2-14）上式表明，S平面上与Z平面的成非线性的正切关系，如图2-4所示。由图2-4看出，在零频率附近，模拟角频率与数字频率之间的变换关系接近于线性关系；但当进一步增加时，增长得越来越慢，最后当时，终止在折叠频率=处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，从而消除了频率混叠现象。图2-5双线性变换法的频率变换关系但是双线性变换的这个特点是靠频率的严重非线性关系而得到的，如式（2-12）与图2-4所示。由于这种频率之间的非线性变换关系，就产生了新的问题。首先，一个

23、线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数（这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如图2-5所示。图2-6双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变，这种频率的畸变，可以通过频率的预畸变来加以校正。也就是将临界模拟频率事先加以畸变，然后经变换后正好映射到

24、所需要的数字频率上。3 IIR滤波器的MATLAB设计用MATLAB进行模拟原型的数字滤波器的设计，一般步骤如下:（1）按一定规则将给出的数字滤波器的技术指标转换成模拟低通滤波器的技术指标；（2）根据转换后的技术指标使用滤波器阶数选择函数，确定最小阶数N和固有频率Wn，根据选用的模拟低通滤波器的类型可分别用: buttord,cheblord，cheb2ord，ellipord等函数；（3）运用最小阶数N产生模拟滤波器原型，模拟低通滤波器的创建函数有:buttap,cheblap,cheb2ap,ellipap,besselap等；（4）运用固有频率Wn把模拟低通滤波器原型转换成模拟低通、高通

25、、带通、带阻滤波器，可分别用函数lp2lp,lp2hp,lp2bp,lp2bs；（5）运用冲激响应不变法或双线性变换法把模拟滤波器转换成数字滤波器，分别用函数impinva和bilinear来实现。低通Chebyshevl型数字滤波器的设计:设计中需要限定其通带上限临界频率Wp,阻带临界滤波频率Ws，在通带的最大衰减Rp，阻带的最小衰减Rs。其步骤如图3-1所示。数字滤波器的参数指标参数变化模拟滤波器的参数数字滤波器模拟滤波器离散化标模拟滤波器的设计图3-1 数字滤波器设计步骤3.1巴特沃斯数字滤波器的设计 （4.1）式中，N为正整数，称为滤波器的阶数。N值越大，通带和阻带的近似特性就越好。在

26、截止频率处，幅度平方响应为=0处的1/2，相当于幅度响应的3dB衰减点。其系统函数为 （4.2）式中，为归一化常数，一般；为s平面左半平面的极点。低通巴特沃斯滤波器设计步骤如下：确定阶数N。 取N=4求极点， 归一化低通原型系统函数为由N=4直接查表得到：极点：归一化低通滤波器系统函数为式中，0.0000,0.0999,0.1914,0.0252将去归一化最终得到通过计算可以总结出过程太麻烦，而且容易出错，结果不直观。3.1.1数字高通滤波器的设计，通带最大衰减为1dB，阻带最小衰减为100dB。用脉冲响应不变法设计巴特沃斯高通滤波器，程序如下：Wp=0.6*pi/T; Ws=0.4*pi/T

27、; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wn=buttord(Wp,Ws,Ap,As,s);%调用butter函数确定巴特沃斯滤波器阶数B,A=butter(N,Wn,high,s); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,1000*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W

28、/pi,abs(Hz); %绘出巴特沃斯数字高通滤波器的幅频特性曲线grid on;xlabel(Frequency/Hz);ylabel(Magnitude);图3-2脉冲响应不变法设计的巴特沃斯高通滤波器运行波形用双线性变换法设计巴特沃斯高通滤波器，程序如下： Wp=0.6*pi/T; Ws=0.4*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wn=buttord(Wp,Ws,Ap,As,s);%调用butter函数确定巴特沃斯滤波器阶数B,A=butter(N,Wn,hi

29、gh,s); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出巴特沃斯数字高通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-3双线性变换法设计的巴特沃斯高通滤波器运行波形3.1.2数字带通滤波器的设计，通带最大衰减为1dB，阻带最小衰减为100dB。用

30、脉冲响应不变法设计巴特沃斯带通滤波器，程序如下：Wp1=0.6*pi/T; Wp2=0.7*pi/T; Ws1=0.3*pi/T;Ws2=0.8*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wn=buttord(Wp,Ws,Ap,As,s);%调用butter函数确定巴特沃斯滤波器阶数B,A=butter(N,Wn,bandpass,s); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,1000*pi); %指

31、定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出巴特沃斯数字带通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-4脉冲响应不变法设计的巴特沃斯带通滤波器运行波形用双线性变换法设计巴特沃斯带通滤波器，程序如下：Wp1=0.6*pi/T; Wp2=0.7*pi/T; Ws1=0.3*pi/T;Ws2=0.8*pi/T; %设置归一化通带和阻带截止频

32、率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wn=buttord(Wp,Ws,Ap,As,s);%调用butter函数确定巴特沃斯滤波器阶数B,A=butter(N,Wn,bandpass,s); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,1000*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性转换法 Hz=freqz(D,C,W); %返回频率响

33、应 plot(W/pi,abs(Hz); %绘出巴特沃斯数字带通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-5双线性变换法设计的巴特沃斯带通滤波器运行波形 3.1.3数字带阻滤波器的设计，通带最大衰减为1dB，阻带最小衰减为100dB。用脉冲响应不变法设计巴特沃斯带阻滤波器，程序如下：Wp1=0.2*pi/T; Wp2=0.8*pi/T; Ws1=0.3*pi/T;Ws2=0.7*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大

34、和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wn=buttord(Wp,Ws,Ap,As,s);%调用butter函数确定巴特沃斯滤波器阶数B,A=butter(N,Wn,stop,s); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,1000*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出巴特沃斯数字带阻滤波器的幅频特性曲线grid on;xlabel

35、(Frequency);ylabel(Magnitude);图3-6脉冲响应不变法设计的巴特沃斯带阻滤波器运行波形用双线性变换法设计巴特沃斯带阻滤波器，程序如下：Wp1=0.2*pi/T; Wp2=0.8*pi/T; Ws1=0.3*pi/T;Ws2=0.7*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wn=buttord(Wp,Ws,Ap,As,s);%调用butter函数确定巴特沃斯滤波器阶数B,A=butter(N,Wn,stop

36、,s); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出巴特沃斯数字带阻滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-7双线性变换法设计的巴特沃斯带阻滤波器运行波形3.2切比雪夫数字滤波器的设计切比雪夫型滤波器在通带幅度特性是等波纹的，在阻带是单调的

37、。切比雪夫滤波器在通带是单调的，在阻带是等波纹的。切比雪夫型低通滤波器的幅度平方函数为：上式结果是一个正实数，式中的是一个小于1的正数，它与通带的波纹有关，越大，波纹越大，为通带的截止频率，时N阶切比雪夫多项式，定义为：切比雪夫滤波器的传递函数参数确定：，N。当N大于或等于1时，切比雪夫多项式递推公式为：切比雪夫滤波器的阶数N等于通带等幅波动的次数，即最大值和最小值的总数，可由下式确定， (4.6)其中均为与衰减有关的参数，Ap为通带的最大衰减，As为阻带的容许衰减，它们的单位为。因为滤波器的幅度平方函数为：进而化简得到：，又因为，由式上式得到：故滤波器的阶数N为： 若要求阻带截止频率上的衰减

38、越大，或过渡带幅度特性越陡，则所需的阶数N越高。从式中可知切比雪夫滤波器的幅度平方函数也只有极点没有零点，且只需求出左半平面的极点，所以一旦求出，N，就能确定。3.2.1数字高通滤波器的设计，通带最大衰减为1dB，阻带最小衰减为100dB。用脉冲响应不变法设计切比雪夫1型高通滤波器，程序如下：Wp=0.6*pi/T; Ws=0.4*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wc= cheb1ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫滤波器阶数B,A=ch

39、eby1(N,Ap,Wc,high,s); %调用 函数设计切比雪夫滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫数字高通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-8脉冲响应不变法设计的切比雪夫1型滤波器运行波形用双线性变换法设计切播学府1型高通滤波器，程序如下：Wp=

40、0.6*pi/T; Ws=0.4*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wc= cheb1ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫滤波器阶数B,A=cheby1(N,Ap,Wc,high,s); %调用 函数设计切比雪夫滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=freqz(D,C,W);

41、 %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫数字高通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-9双线性变换法设计的切比雪夫1型滤波器运行波形用脉冲响应不变法设计切比雪夫2型高通滤波器，程序如下：Wp=0.6*pi/T; Ws=0.4*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wc= cheb2ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫2滤波器阶数B,A=ch

42、eby2(N,Ap,Wp,s); %调用 函数设计切比雪夫2滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫2数字高通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-10脉冲响应不变法设计的切比雪夫2型滤波器运行波形用双线性变换法设计切比雪夫2型高通滤波器，程序如下：Wp=0.

43、6*pi/T; Ws=0.4*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wc= cheb2ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫2滤波器阶数B,A=cheby2(N,Ap,Wp,s); %调用 函数设计切比雪夫2滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=freqz(D,C,W); %返回频

44、率响应 plot(W/pi,abs(Hz); %绘出切比雪夫2数字高通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-11双线性变换法设计的切比雪夫2型滤波器运行波形3.2.2数字带通滤波器的设计，通带最大衰减为1dB，阻带最小衰减为100dB。用脉冲响应不变法设计切比雪夫1型带通滤波器，程序如下:Wp1=0.5*pi/T; Wp2=0.6*pi/T; Ws1=0.3*pi/T;Ws2=0.7*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设

45、置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wc= cheb1ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫滤波器阶数B,A=cheby1(N,Ap,Wc,bandpass,s); %调用 函数设计切比雪夫滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫数字带通滤波器的幅频特性曲线grid on;xl

46、abel(Frequency);ylabel(Magnitude);图3-12脉冲响应不变法设计的切比雪夫1型带通滤波器滤波器运行波形用双线性变换法设计切比雪夫1型带通滤波器，程序如下:Wp1=0.5*pi/T; Wp2=0.6*pi/T; Ws1=0.3*pi/T;Ws2=0.7*pi/T; %设置归一化通带和阻带截止频率T=1; fs=1/T; Wp=Wp1,Wp2;Ws=Ws1,Ws2;Ap=1;As=100; N,Wc= cheb1ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫滤波器阶数B,A=cheby1(N,Ap,Wc,bandpass,s); %调用 函数设计切

47、比雪夫滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫数字带通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude); %设置通带最大和最小衰减图3-13双线性变换法设计的切比雪夫1型带通滤波器滤波器运行波形用脉冲响应不变法设计切比雪夫2型带通滤波器，程序如下：Wp1=0.5*pi/T; Wp

48、2=0.6*pi/T; Ws1=0.3*pi/T;Ws2=0.7*pi/T; %设置归一化通带和阻带截止频率T=1; fs=1/T; Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wc= cheb2ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫2滤波器阶数B,A=cheby2(N,Ap,Wc,bandpass,s); %调用 函数设计切比雪夫2滤波器W=linspace(0,pi,1000*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %

49、调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); grid on;xlabel(Frequency);ylabel(Magnitude);图3-14脉冲响应不变法设计的切比雪夫2型带通滤波器滤波器运行波形用双线性变换法设计的切比雪夫2型带通滤波器，程序如下：Wp1=0.6*pi/T; Wp2=0.7*pi/T; Ws1=0.3*pi/T;Ws2=0.8*pi/T; %设置归一化通带和阻带截止频率T=1; fs=1/T; Wp=Wp1,Wp2;Ws=Ws1,Ws2;Ap=1;As=100; N,Wc= cheb2ord(Wp,Ws,Ap,

50、 As,s);%调用 函数确定切比雪夫2滤波器阶数B,A=cheby2(N,Ap,Wp,bandpass,s); %调用 函数设计切比雪夫2滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); grid on;xlabel(Frequency);ylabel(Magnitude);图3-15双线性变换法设计的切比雪夫2型带通滤波器滤波器运行波形3.2.3数字带阻滤

51、波器的设计，通带最大衰减为1dB，阻带最小衰减为100dB。用脉冲响应不变法设计切比雪夫1型带阻滤波器：程序如下：Wp1=0.3*pi/T; Wp2=0.7*pi/T; Ws1=0.4*pi/T;Ws2=0.6*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wc= cheb1ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫滤波器阶数B,A=cheby1(N,Ap,Wc,stop,s); %调用 函数设计切比雪夫滤波器W=l

52、inspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫数字高通滤波器的幅频特性曲线grid on;xlabel(Frequency);ylabel(Magnitude);图3-16冲激响应比变法设计的切比雪夫1型带阻滤波器运行波形用双线性变换法设计切比雪夫1型带阻滤波器，程序如下：Wp1=0.3*pi/T; Wp2=0.7*pi/T; Ws1=0.4*pi/T

53、;Ws2=0.6*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wc= cheb1ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫滤波器阶数B,A=cheby1(N,Ap,Wc,stop,s); %调用 函数设计切比雪夫滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=

54、freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫数字带阻滤波器的幅频特性曲grid on;xlabel(Frequency);ylabel(Magnitude);图3-17双线性变换法设计的切比雪夫1型带阻滤波器运行波形用脉冲响应不变法设计切比雪夫2型带阻滤波器，程序如下：Wp1=0.3*pi/T; Wp2=0.7*pi/T; Ws1=0.4*pi/T;Ws2=0.6*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;W

55、s=Ws1,Ws2;N,Wc= cheb2ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫2滤波器阶数B,A=cheby2(N,Ap,Wp,stop,s); %调用 函数设计切比雪夫2滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫2数字带阻滤波器的幅频特性曲线grid on;xlabel(Frequency/pi);ylab

56、el(Magnitude);图3-18脉冲响应不变法设计的切比雪夫2型带阻滤波器运行波形用双线性变换法设计设计切比雪夫2型带阻滤波器，程序如下：Wp1=0.3*pi/T; Wp2=0.7*pi/T; Ws1=0.4*pi/T;Ws2=0.6*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减Wp=Wp1,Wp2;Ws=Ws1,Ws2;N,Wc= cheb2ord(Wp,Ws,Ap, As,s);%调用 函数确定切比雪夫2滤波器阶数B,A=cheby2(N,Ap,Wp,stop,s); %调

57、用 函数设计切比雪夫2滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=bilinear(B,A,fs); %调用双线性变换法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出切比雪夫2数字带阻滤波器的幅频特性曲grid on;xlabel(Frequency/pi);ylabel(Magnitude);图3-19双线性变换法设计的切比雪夫2型带阻滤波器运行波形3.3 椭圆形滤波器的设计椭圆低通滤波器采用有限零点设计来逼近理想的低通滤波器的特性。因为零点的

58、位置与椭圆函数的许多特性有关，所以称之为椭圆低通滤波器。它的幅度平方函数为： (4.7)式中的是雅可比椭圆函数，是与通带衰减有关的参数，阶数N等于通带和阻带最大点和最小点的总和。阶数和系统函数由通带的截止频率，通带最大衰减，阻带截止频率与阻带的最小衰减等参数确定。双线性变换法设计步骤如下：（1）首先写出该滤波器的系统函数（2）利用双线性变换法转换，数字滤波器的系统函数为3.3.1数字高通滤波器的设计，通带最大衰减为1dB，阻带最小衰减为100dB。用脉冲响应不变法设计椭圆形高通滤波器，程序如下：Wp=0.6*pi/T; Ws=0.4*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采

59、样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wn=ellipord(Wp,Ws,Ap,As,s);%调用函数确定椭圆滤波器阶数B,A=ellip(N,Ap,As,Wn,high,s); %调用ellip函数设计椭圆滤波器W=linspace(0,pi,400*pi); %指定一段频率值 hf=freqs(B,A,W); %计算模拟滤波器的幅频响应 D,C=impinvar(B,A,fs); %调用脉冲响应不变法 Hz=freqz(D,C,W); %返回频率响应 plot(W/pi,abs(Hz); %绘出椭圆数字高通滤波器的幅频特性曲线

60、grid on;xlabel(Frequency/pi);ylabel(Magnitude);图3-20脉冲响应不变法设计的椭圆形高通滤波器运行波形用双线性变换法设计的椭圆形高通滤波器，程序如下：Wp=0.6*pi/T; Ws=0.4*pi/T; %设置归一化通带和阻带截止频率T=1; %设置采样周期为1fs=1/T; %采样频率为周期倒数Ap=1;As=100; %设置通带最大和最小衰减 N,Wn=ellipord(Wp,Ws,Ap,As,s);%调用函数确定椭圆滤波器阶数B,A=ellip(N,Ap,As,Wn,high,s); %调用ellip函数设计椭圆滤波器W=linspace(0,