毕业论文GSM移动通信系统中低通滤波器的设计与仿真16791.doc

《毕业论文GSM移动通信系统中低通滤波器的设计与仿真16791.doc》由会员分享，可在线阅读，更多相关《毕业论文GSM移动通信系统中低通滤波器的设计与仿真16791.doc（53页珍藏版）》请在三一办公上搜索。

1、西安欧亚学院 本科毕业论文（设计） 题 目： GSM移动通信系统中低通滤 波器的设计与仿真 学生姓名： xxx 指导教师： xxx 所在分院： 信息工程学院 专 业： 通信工程 班 级： 统本通工0901 二O一三 年 五月GSM移动通信系统中低通滤波器的设计与仿真摘要：本文应用MATLAB软件中的Simulink平台分别研究了FIR数字低通滤波器和IIR数字低通滤波器的设计方法和实现方法，并进行模拟仿真设计，以显示本文所介绍数字低通滤波器的多种特性。在无限脉冲响应（IIR）数字低通滤波器设计中，应用了冲击响应不变法和双线性Z变换法，先进行模拟滤波器的设计，然后把模拟滤波器转换为数字滤波器，最

2、后进行滤波器的频带转换，使用MATLAB软件按两种方法编辑程序并通过图形展示IIF低通滤波器的各种参数和性能指标。在有限脉冲响应（FIR）低通数字滤波器设计中，讨论了FIR线性相位滤波器的特点和用窗函数法设计FIR滤波器两个问题。这两类低通滤波器整个设计过程都是按照理论分析、编程设计、模拟仿真、具体实现的步骤进行的。关键词：数字低通滤波器、MATLAB、SimulinkGSM mobile communication system design and simulation of the low-pass filterAbstract：In this paper, the applicatio

3、n of MATLAB software Simulink platform respectively studied the FIR digital low-pass filter and IIR low-pass filter design methods and implementation method, and carries on the simulation design,to display a variety of characteristics of digital low-pass filter, presented in this paper. In infinite

4、impulse response (IIR) digital low-pass filter design, the application of shock response is not political reform and bilinear z-transform method, first for the design of analog filter, then converted to a digital filter, the analog filter to filter the frequency conversion, using MATLAB software acc

5、ording to the two methods of editing program and through the graphic display of IIF low-pass filter parameters and performance indicators. In finite impulse response (FIR) low-pass digital filter design, discussed the characteristics of the FIR linear phase filter and FIR filter with window function

6、 method with two problems. These two types of low-pass filter the whole design process is in accordance with the theoretical analysis, programming design, simulation and implementation steps.Keywords: digital low-pass filter, the MATLAB, Simulink,目 录1 绪论11.1 课题研究的意义11.2 国内外研究现状11.3 毕业论文的内容11.4 MATLA

7、B软件及Simulink平台的概述22 IIR数字滤波器设计及其MATLAB实现52.1 IIR数字滤波器的原理及特点52.2 IIR滤波器的结构52.2.1 直接型52.2.2 级联型72.2.3 并联型82.3 常用的滤波器原型82.3.1 巴特沃斯低通滤波器92.3.2 切比雪夫低通滤波器92.4 模拟数字滤波器变换及其MATLAB实现112.4.1 脉冲响应不变法112.4.2 双线性Z变换法142.4.3 脉冲不变法和双线性Z变换法的比较172.5 滤波器的阶数与频率关系183 FIR滤波器203.1 FIR滤波器的结构203.1.1 直接型结构203.1.2 级联型结构203.1.

8、3 频率抽样型结构213.1.4 快速卷积型结构213.2 常用窗函数介绍213.2.1 矩形窗(Rectangular window)223.2.2 三角窗(Triangular window)223.2.3 汉宁窗(Hanning window)223.2.4 海明窗(Hamming window)223.2.5 布拉克曼(Blackman window)233.3 各种窗函数的实现与比较233.4 海明窗设计低通滤波器244 图形用户界面设Simulink设计数字滤波器284.1 Simulink滤波器设计与分析工具（FDATOOL）的设计技术284.2 建立滤波器系统的Simlink模

9、型并进行仿真分析294.2 .1 建立用数字滤波器设计模块实现滤波器系统的仿真模型304.2.2 系统模块参数设置304.2.3 系统仿真参数的设置314.3 系统仿真与分析324.4 仿真结果的演示324.4.1 Simulink下滤波器设计技术应用324.4.2 IIR和FIR滤波器比较375 利用MATLAB辅助DSP实现滤波器的方法385.1 数字信号处理概述385.2 以巴特沃斯滤波器为例用MATLAB辅助DSP设计滤波器385.2.1 设计巴特沃斯滤波器程序滤波器385.2.2 DSP的系数对称滤波器的实现方法405.2.3 DSP滤波算法实现405.2.4 MATLAB辅助设计4

10、06 总结42致谢43参考文献441 绪论1.1 课题研究的意义随着科学的突飞猛进的发展，一些数字滤波器被广泛使用在多准数字信号处理领域中。比如，在通信工程、图像的编码、语音编码、雷达、HDTV(高清晰度电视)等。在实际学习中对滤波器比较注重理论，而其本身比较抽象、复杂、公式繁多、不易理解，为了更好的增强我们的认识理解，本论文通过对FIR与IIR两种滤波器的对比分析，将所涉及到的多种类型的滤波器逐个系统的仿真学习和分析，提高学习的效率。应用软件对GSM中的低通滤波器的设计与仿真，掌握滤波器的设计方法，对输入滤波器的模拟信号或数字信号进行处理，对模拟信号抽样量化后进行适当的限定，分析滤波器的参数

11、对滤波器性能的影响，并掌握其结构功能等。最终使用MATLAB平台实现滤波器的仿真，并完成论文。1.2 国内外研究现状数字滤波器在数字信号处理中有广泛的应用，将对滤波器作详细的介绍。如果滤波器的输入信号、输出信号都是离散的时间信号，那么，该滤波器的冲击响应也必然是离散的，我们称这样的滤波器为数字滤波器。当只用硬件实现一个数字滤波器时，所需的元件是延迟器、乘法器和加法器等。当用计算机软件实现时，它就是一段线性褶积（或卷积）的程序。然而，模拟滤波器只能用硬件的方法来实现，其元器件是电阻、电感、电容及运放或开关等。因此，数字滤波器的实现要比模拟滤波器实现更加简单，而且易获得更加理想的效果。滤波器的种类

12、有非常多，分类方法也不同，例如可以从实现方法、设计方法、功能上来分等等。但从总体来看的话，滤波器可以现代滤波器和分经典滤波器两大类。经典滤波器是假定输入信号x(n)中的有效信号和噪声（或干扰）信号成分各自在不同的频带内，当x(n)通过滤波系统后，可以将噪声信号和有用信号有效的分开，同时抑制噪声信号。但是，有效信号和噪声信号的频率带一般会有相互重叠，这样的话经典的滤波器将无法有效区分。有时滤波效果较好，有时效果较差。现代滤波理论研究的主要内容是从含有噪声的数据记录中估测出信号的某些特征或信号本身。如果信号被估计出，那么估计出的信号会有高的信噪比。现代滤波器把信号和噪声都视为随机信号，利用统计学的

13、方法推导出来一套最佳的估值算法，然后结合硬件和软件实现滤波。目前现代滤波器主要有：线性预测器、卡尔曼滤波器、维纳滤波器、自适应滤波器等。而在国内将智能控制引入到传统控制方法中，如无差拍控制，神经网络控制等，建立统一能用于实际工程中的滤波器模型是智能控制待解决的问题。1.3 毕业论文的内容 GSM移动通信系统的数字滤波器的设计，GSM移动通信系统中常见的是低通滤波器，其中低通滤波器适用于数字基带信号的滤波， GSM系统的一个信道带宽为200KHz。利用MATLAB模拟GSM移动通信的信道和信号，并将其仿真。低通滤波器可分为FIR与IIR两种滤波器，其中FIR滤波器可用窗函数、巴特沃斯、契比雪夫等

14、方法实现，而IIR滤波器可用巴特沃斯、契比雪夫和椭圆滤波器等仿真，具体内容可分为以下几章：1 介绍课题的意义及制作背景；2 MATLAB软件的相关介绍，以及在使用中的要点和出现的问题分析；3 介绍GSM移动通信系统的原理，根据其特点推导分析其性能指标；4 用MATLAB仿真其信道，再输入一个信号，做出输出信号，对整个过程分析（图像、公式等）；5 通过仿真描述在实际应用中的所需的器件，所注意的问题，应该避免的错误，以及解决方法；6 对整个流程进行分析总结，以及个人通过论文写作提出对滤波器的展望等。 设计原理图1-1所示：输出信号低通滤波器输入信号图 1-1 滤波器设计流程图1.4 MATLAB软

15、件及Simulink平台的概述 MATLAB名字是由Matrix和Laboratory两词中的前三个字母组合而成，它是由美国MathWorks公司所推出，用于数值计算和图形处理的数学计算环境。MATLAB环境下，用户可以集成的进行数值计算、程序设计、输入输出、图形绘制、文件管理等多项操作。它优秀的数值计算能力和强大的数据可视化能力使其在同类软件中独占鳌头。MATLAB系统最初由Cleve Moler用FORTRAN语言设计，现在的MATLAB程序是由MathWorks公司用C语言开发的。它的第一版发行于1984年，经过多年的改进，版本不断更新和升级，其所包含的工具箱功能也越来越丰富，应用越广泛

16、，使用越来越方便。MATLAB语言之所以能如此迅速的普及推广，是由于它有着不同于其他语言的特点，归纳如下：1、 简单易学：MATLAB是一个开发软件，也是一门编程语言。其语法简单易学，就像在演算纸上列出公式与求解问题，所以被称为演算纸式科学算法语言。 2、计算功能强大：MATLAB拥有庞大的统计、数学及工程函数，用户可通过MATLAB平台很快的实现所需的强大数学计算功能。通过长期的研究和开发的数值计算程序，使用了可靠、成熟、安全的算法。此外，MATLAB平台中还有很多个工具箱，可用来处理应用中的数学、工程等方面的问题。3、具有先进的可视化工具：MATLAB提供功能强大的和交互式的平面和立体的绘

17、图的功能，用户可以创直观的和动态的彩色图形。可视化的工具还包括了图像显示、光源、曲面渲染、伪彩图线框图、动画等等。4、开放性和可扩展性强：M文件是可见的MATLAB上的程序，所有用户可以查看源代码。用户能够测试各种算法的准确，修改原有的函数，或者加入自行设计的新建的部件和子模块。5、能够特殊应用的工具箱：MATLAB的工具箱还强大了对工程及科学中的很多特殊的应用的支持。此中的工具箱还和MATLAB相同是完全用户化，可扩展性是非常强的。将许多个工具箱和MATLAB二者联合起来使用，还能够得到一个具有强大功能的计算组合包，能够满足每个用户的特殊要求8。MATLAB数字信号处理工具箱和滤波器设计工具

18、箱专门应用于数字信号处理领域。工具箱提供了丰富而简便的设计，使原本很繁琐的程序设计简化成为函数的调用。只要以正确的指标参数调用相应的滤波器的设计好的程序或者工具箱中的函数，就将重新得到正确的结果，使用也非常的方便。Simulink是MATLAB最重要的组件之一，它提供一个动态的系统的仿真、建立模型和综合性能分析的集成环境。MATALB的环境中用户可以不用书写大量的具体程序，用户可以通过既简单而又直观的操作，就能构创建出繁杂的系统。Simulink还具有结构清晰和流程明了、适应范围广及仿真结构精细、使用效率高、与实际贴近、灵活多变等诸多的优点，通过以上优点Simulink被大量应用于数字信号处理

19、和控制理论中的相繁杂多变的系统仿真和设计中去。同时它还有很多的第三方的辅助硬件和软件可应用于或被要求应用于Simulink。Simulink编辑的特点：（1） 强大的可扩充的预定义模块库；（2） 予以设计功能的层次性来分割模型，实现对复杂的设计的管理；（3） 应用Model Explorer 搜索、配置、创建模型、导航中的随机信号、任意参数、属性，产生出模型的代码；（4） 具有图形化调试器，判断系统设计的性能和系统的异常行为，应用剖析器来检查仿真结果；（5） 其中的API工具用于和别的仿真程序的链接或者是和集成手写代码功能；（6） Embedded模块在Simulink平台与嵌入式的系统运行中

20、直接调用MATLAB的相应算法；（7） 通过交互式的图像处理编辑器来重新组合、方便管理直观的各种模块图；（8） 也可以应用固定步长、可变步长来仿真，也可以通过仿真模式来具体解释性的方式来运行、也可以通过编译出C语言代码来运行整个系统仿真的模型；（9） 还可以通过访问MATLAB程序来对系统的结果进行具体的分析和可视化的效果，来具体的定制建模的环境，配置信号参数和检测数据模型，判断和处理工具来确保模型的统一性，找出模型中的错误之处7。（10） 在学习过程中应用Simulink的内容, 不仅可以帮助使用者对新问题深入思考解决, 并对问题进行建模直观的分析, 同时使所学到的知识能够学以致用。2 II

21、R数字滤波器设计及其MATLAB实现2.1 IIR数字滤波器的原理及特点 IIR滤波器，即为无限脉冲响应滤波器，其单位脉冲响应是无限长的，其系统的传递函数如公式（2-1）所示：（）（2-1）系统的差分方程可以写成如公式（2-2）所示：（2-2）通过系统函数可以看出，IIR滤波器有几个显著的特点：（1） IIR型的滤波器具有不为零的极点和零点。为保证滤波器成为稳定的系统，必须使系统的极点只能在单位圆的内部，系统的稳定性也就取决于系统函数的极点。（2） 因为线性相位滤波器中的各个零点、极点是关于单位圆而对称的，所以极点只能被允许在单位圆的原点上。因为IIR型的滤波器存在不为零的极点，所以只可能实现

22、近似的线性相位特性。也正因为IIR型的滤波器的非线性相位特性限制了起应用范围1。2.2 IIR滤波器的结构 数字滤波器中采用的数字器件有：加法器、延迟器、乘法器等，选择网络结构就是考虑整个系统的技术要求、使用环境和系统的稳定性、精度、误差、经济性及运算速度等。由于乘法运算较费时间，因此通常希望网络中乘法器尽可能少。另外，在用硬件实现时，寄存器和存储器的字长通常是十分有限，有限字长对数字滤波器实现的影响和结构有关，有时乘法器和延迟器的数目可能不一定最少，但它对有限字长的影响的敏感程度较低。因此实现时需根据实际要求来选用。2.2.1 直接型 IIR型的数字滤波器系统函数的一般形式和差分方程的一般形

23、式是 （2-3） （2-4） 由以上公式可以得出IIR型的数字滤波器的直接结构型式，如图 2-1。 图 2-1 直接型 图 2-2 直接型 由公式（2-3）和图 2-1可见，H(z)也可视为 与 两个子系统的串联形式，它的结构和个部件的串联的次序没有关系。所以，我们把它还可改写成图 2-2所示的结构，如若再合并其中中间的两条延时链的话，则得到图 2-2所示的另一种直接型结构了，这个结构与图 2-1 的结构来做对比，它的延时单元可以节省一半，这种系统类存容量可以节省系统的结构型式，这种结构也被称为正准型结构。 以上的两种结构均有一些缺点，首先缺点是系数、对滤波器性能的控制关系不直接，所以调整极为

24、不方便。更加严重的是，当滤波器的阶数较高时，其极点的位置灵敏度非常的大，对字长效应过于敏感。因而容易出现不稳定现象并产生较大的误差。极点对于系数的变化过于灵敏，从而使系频率响应对系数变化过于灵敏了，即对有限精度运算也过于灵敏，一般不稳定或产生过大的系统误差。因此，相对来说以下的两种结构比较好点1。2.2.2 级联型 N阶的系统函数也可用它的零点和极点来表示，如公式（2-5）所示： (2-5) 因为的系数、是实系数，因此此函数的零极点当且有两种情况，、为实零点和实极点，、 表示复数的共轭的零点，、表示复数的共轭的极点。把其中的每对共轭因子合并到一起就能构成具有一个实系数的二阶因子了。 (2-6)

25、 如果把单实根因子也看成是二阶因子的特例，即就是二次项系数或等于零的二阶因子，那么整个函数就可完全分解成实系数的二阶因子的一种形式。滤波器也就可用几个二阶网络级联在一起来构成了，这些二阶网络也就被称为滤波器的二阶基本节，它的传递函数为 (2-7)假设其基本节用直接型来实现的话，整个滤波器就是它们的级联。如图 2-3所示。 图 2-3 级联型滤波器从整个级联型结构来看的话，它的每个基本节的都会关系到滤波器中的具体的某对极点或者是某对零点，调整函数的系数、只单独调整了滤波器第i对零点但是并不影响其他任何的零点。同样，调整系数、也只能是单独的调整了滤波器的第i对极点。因此，这种结构便于准确地实现单节

26、滤波器中的零极点，也便于性能的调整和控制。然而总体上来看，某频率成分在某一节阻断，也就是表现为整体上的阻断，必须在每一节畅通才能表现为整体畅通，因此，这种结构更加容易控制零点而很难调整极点。此外，级联结构也便于改变各子系统的运算的次序，优化滤波器的性能，减小运算中的误差。而其局限性在于，一般需要调整各级的质监电平的放大、缩小，使变量不能太大或者太小。在定点运算中，变量太大非常容易产生溢出，变量太小也会使信噪比太小1。也可以看出，这种结构可有许多不同的零极点搭配关系，基本节级联前后次序可不同，它们在实际工作时，由于精度限制，不同方案运算所带来的误差也是不同的，因此排列组合，将决定滤波器优劣。如何

27、才能得到最好的搭配，就是一个最优化问题。2.2.3 并联型将系统的传递函数展开成为部分分式之和的形式：（2-8）其中.此时H（z）可视为多个一阶和二阶滤波器的并联，若，则公式(_2-8)中不包括最后的一项 。其结构如图 2-4所示。图 2-4 并联型滤波器 从整体看，并联结构可以单独调整数字滤波器的极点，运算速度快，各二阶网路的误差均是互不影响的，总误差较小，对系统字长要求不低。但是通常是不能直接调整零点，所以多个二阶子系统的零点并不是整个系统函数的零点，如果需要较为准确的传输零点时，级联型是最合适不过的3。 2.3 常用的滤波器原型IIR滤波器设计的技术依靠现有的模拟滤波器来得到数字滤波器，

28、实际的工程当中把模拟滤波器称之为滤波器原型。实际工程中应用广泛的滤波器的有以下几种。2.3.1 巴特沃斯低通滤波器巴特沃斯低通滤波器的特征是通带和阻带都有平坦的幅度响应。N阶低通滤波器的平方幅度响应为 （2-9）试中：N为阶数：为截止频率（）4。 MATLAB提供了buttap函数设计N阶归一化巴特沃斯模拟滤波器。如下为5阶的低通滤波器的MATLAB程序：%巴特沃斯低通滤波器MATLAB实现clear allz,p,k=buttap(5); %阶数为5阶巴特沃斯低通滤波器num,den=zp2tf(z,p,k);freqs(num,den);plot(num,den)图 2-5 阶数为5巴特沃

29、斯滤波器2.3.2 切比雪夫低通滤波器 切比雪夫滤波器则有两种，切比雪夫型滤波器在能带中具有等波动响应，而切比雪夫型滤波器在阻带具有等波动性。切比雪夫型滤波器的平方幅度响应为 （2-10）式中：N为滤波器阶数；为通带波动系数，它与有关切比雪夫型滤波器的平方幅度响应如公式（2-9）所示： （2-11）式（2-10）中参数与（2-11）参数均相同5。 以下为切比雪夫低通滤波器：%切比雪夫I型低通滤波器clear all; n=0:0.01:2;for i=1:2 switch i case 1 N=7; case 2 N=8; endRp=1; %通带纹波系数为1dBz,p,k=cheb1ap(N

30、,Rp); %函数cheb1-设计切比雪夫I型低通滤波器b,a=zp2tf(z,p,k); H,w=freqs(b,a,n); %函数freqs-求解模拟滤波器频率响应magH2=(abs(H).2;%绘制图形posplot=12 num2str(i); %输出图形subplot(posplot) %函数subplot-创建子图图区plot(w,magH2)axis(0 2 0 1);xlabel(w/wc);ylabel(|H(jw)|2);title(N= num2str(N);end图 2-6 不通阶数的切比雪夫滤波器2.4 模拟数字滤波器变换及其MATLAB实现设计完模拟低通滤波器后，

31、我们就可以以模拟滤波器为原型转变出数字滤波器。这些变换全都是复值映射，通过数字滤波器保持了模拟滤波器原有的不同特性，研究出不同的变换技术。其中，通用的方法有两种：脉冲响应不变法）和双线性Z变换法1。2.4.1 脉冲响应不变法脉冲响应不变法是让DL的单一抽样响应序列h(n)，仿照AL的脉冲响应而相应的改变g(t)10。假如系统的传递函数是G(s)的AL的单脉冲响应是g(t)的话，给整个系统的脉冲响应g(t)做等间隔的采样，使得DL单抽样响应h(n)和g(t)的采样值正好相等，即： (2-12)其中的为采样周期1。G(s)为AL的系统传递函数，又使H(z)为DL的系统传递函数。采样信号通过拉式变换

32、以后，相对应的采样序列通过Z变换以后被映射的关系式如公式（2-13）所示： (2-13)其中，系统函数G(s)和H(z)的关系如公式（2-14）所示： (2-14)式(2-13)首先将AL的系统函数G(s)做了周期的延迟以后，在通过如公式(2-14)的映射关系来变换，最后再对应的映射到Z平面上去，就可以直接得到DL的系统函数H(z)。且AL和DL之间2的频率满足=T的关系。经过式(2-13)的映射关系，s平面直接可以映射到Z平面上的单位圆内，因此，这个因果稳定的AL映射成为因果稳定的DL5。通过上面的比较和分析，应用脉冲响应不变法，通过AL的系统传递函数，可以直接求得DLD的系统函数，而它的具

33、体的设计方法归纳为：(1)利用=T，将DL指标，转换为AL指标，；(2)通过指标，,来设计AL的G(s)；(3)用部分分式展开的方法，把G(s)展成开，得出如公式（2-15）所示： (2-15)(4)最后再把AL的极点转换为数字极点，得到DL，如公式（2-16）所示： (2-16)通过以上所归纳的理论，在MATLAB的环境下用相应的函数来实现脉冲响应不变法设计一DL。 以下是脉冲响应不变法设计的AL：%脉冲响应不变法wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1; %采样周期%性能指标Rip=10(-Rp/20);Atn=10(-As/20);OmgP=wp*T;OmgS

34、=ws*T;N,OmgC=buttord(OmgP,OmgS,Rp,As,s); %模拟滤波器的阶数cs,ds=butter(N,OmgC,s); %所需的模拟低通滤波器 b,a=impinvar(cs,ds,T); %脉冲响应不变法进行转换 db,mag,pha,grd,w=freqz_m(b,a); %求得相对，绝对频响及相位、群延迟响应%绘制各条曲线subplot(2,2,1);plot(w/pi,mag);title(幅频特性);xlabel(w(/pi);ylabel(|H(jw)|);axis(0,1,0,1.1);set(gca,XTickMode,manual,XTick,0

35、0.2 0.3 0.5 1);set(gca,YTickMode,manual,YTick,0 Atn Rip 1);gridsubplot(2,2,2);plot(w/pi,db);title(幅频特性(dB);xlabel(w(/pi);ylabel(dB);axis(0,1,-40,5);set(gca,XTickMode,manual,XTick,0 0.2 0.3 0.5 1);set(gca,YTickMode,manual,YTick,-40 -As -Rp 0);gridsubplot(2,2,3);plot(w/pi,pha/pi);title(相频特性);xlabel(w(

36、/pi);ylabel(pha(/pi);axis(0,1,-1,1);set(gca,XTickMode,manual,XTick,0 0.2 0.3 0.5 1);gridsubplot(2,2,4);plot(w/pi,grd);title(群延迟);xlabel(w(/pi);ylabel(Sample);axis(0,1,0,12);set(gca,XTickMode,manual,XTick,0 0.2 0.3 0.5 1);grid图 2-7 脉冲不变法设计的低通滤波器2.4.2 双线性Z变换法利用脉冲响应不变法设计DL时，由于=T的频率关系是根据推导而得出的，所以是轴每隔2/T

37、便映射到单位圆上一周，则会引起了频域内的混叠。为克服这一现象，研究出来了另一种映射关系： (2-17) 这种关系称之为双线性Z变换法3。从中可以看出，双线性Z变换法是将整个s平面通过一定的关系压缩到s1平面的一条带宽为2/T的横带里去 ，之后再通过标准的变化关系将横带变换到整个Z平面上去，这样就可以得到S平面与Z平面间的对应关系，整个线性的对应的过程如图2-8所示：图 2-8 双Z变换法过程图由公式(2-17)变换得 (2-18)及 (2-19) (2-20)公式(2-17)和公式(2-18)给出了s和z之间的对应关系，而公式(2-19) 和公式(2-20)给出了和之间的映射关系，但这是一种非

38、映射关系，然而双线性Z变换法正是利用了正切函数的非线性这一特点，把整个j轴压缩到了单位圆的圆周上7。在MATLAB中，双线性Z变换可通过bilinear函数实现,其调用格式为：Bz,Az=bilinear(B,A,Fs)，其中B,A为AL的传递函数G(s)的分子分母多项式系数分量，而其中Bz,Az为数字滤波器的传递函数H(Z)的分子分母多项式的系数分量7。以下是以巴特沃斯滤波器为例，用双Z法设计的滤波器。%巴特沃斯数字低通滤波器的双Z变换%Fs为采样频率；fp为通带上限临界频率；fs为阻带下降临界频率；rp为通带允许最小衰减；rs为阻带内最小衰减fp=3000;fs=4000;Fs=20000

39、; %滤波器参数rp=3;rs=20; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; %计算数字滤波器的设计指标Fs=Fs/Fs; %使Fs为一wap=tan(wp/2);was=tan(ws/2);n,wn=buttord(wap,was,rp,rs,s); %计算数字滤波器的阶数和通带截止频率z,p,k=buttap(n); %计算数字滤波器系统函数bp,ap=zp2tf(z,p,k); %将零点极点增益形式转换为传递函数形式bs,as=lp2lp(bp,ap,wap); %将模拟原型低通滤波器的截止频率转换为4khzbz,az=bilinear(bs,as,Fs/2); %双

40、线性Z变换调用函数H,w=freqz(bz,az,256,Fs*1000); %计算数字滤波器的幅频响应plot(w,20*log10(abs(H);axis(0 500 -350 50)xlabel(omega/pi);ylabel(幅度/dB);title(模拟滤波器指标设计的数字滤波器)grid on;图 2-9 巴特沃斯为例设计双Z变换法2.4.3 脉冲不变法和双线性Z变换法的比较 以AL的理论为原型设计的DL，也就是说使DL能模仿AL的特性，这种模仿方法可从不同的角度出发。 脉冲不变法有如下重要特点：（3） 时域脉冲响应能模仿模拟滤波器时，通常是使用脉冲响应不变法。 （4） 而脉冲响

41、应不变法的特点是频率坐标的变换是线性的。（5） 假如Ha(s)是稳定的，映射到H(Z)也是稳定的。 （6） 脉冲响应不变法的缺点是有频谱周期延拓。 而双线性变换有如下特点：（1） 与脉冲响应不变法相比，双线性变换要靠频率的非线性关系得到S平面与Z平面的单值的一一对应关系（2） 双线性变换法的缺点是和的非线性关系，会导致DL的幅频响应相对于AL的幅频响应有一定的畸变。（3） 线性相位的AL经过双线性Z变换之后，滤波器也就不会再有线性相位特性。即便是双线性变换有这样的缺点，但它目前仍然还是使用得最为普遍、最有成效的一种设计工具。（4） 虽然发生了非线性变化，但其幅频特性还是保持了分段常数的特性。 （5） 双线性变换设计方法比脉冲响应法的设计方法计算更加直接、更加简单。 由于s与z之间的简单的代数关系，所以从AL传递函数可直接通过代数置换得到DL的传递函数6。 以下为脉冲响应不变法和双线性变换法比较：%脉冲响应不变法和双线性变换法比较