1 数字滤波器.docVIP

摘要 数字滤波技术是数字信号处理的一个重要组成部分，滤波器的设计是信号处理的核心问题之一。本文根据IIR滤波器的设计原理提出了MATLAB环境下基于双线性变换法和基于脉冲响应不变法来实现IIR数字滤波器。这些方法需要根据数字滤波器的性能指标设计相应的滤波器。MATLAB软件是一个很好的应用软件，MATLAB环境为设计IIR滤波器提供了一个可靠而有效的工作平台。 关键词: IIR数字滤波器；系统函数；脉冲相应不变法；双线性变换法；MATLAB实现 目录 1 数字滤波器 2 1.1 简述数字滤波器 2 1.1.1 数字滤波器的基本知识 2 1.1.2 数字滤波器设计的基本流程 3 1.2 IIR数字滤波器 3 1.2.1 IIR数字滤波器的设计原理 3 1.2.2 脉冲响应不变法 4 1.2.3 双线性变换法 4 1.3 IIR数字滤波器的基本特点 5 1.4 IIR数字滤波器的设计 6 1.4.1 IIR 数字滤波器设计的基本步骤 6 1.4.2 IIR数字滤波器的两种设计方法 7 2 MATLAB 8 2.1 MATLAB的语言特点 8 2.2 MTALAB的功能 9 3 利用MATLAB来实现IIR滤波器 12 3.1 线性相位变换法 12 3.2 脉冲响应不变法的MATLAB实现 16 结束语 23 心得体会…………………………………………………………………………………………..24 参考文献…………………………………………………………………………………………..25 1 数字滤波器 1.1简述数字滤波器 1.1.1 数字滤波器的基本知识 数字信号处理（Digital Signal Processing，简称DSP）是一门涉及许多科学而又广泛应用于许多领域的新兴科学，DSP有两种含义：数字信号处理，数字信号处理器。数字信号处理是将信号以数字方式表示并处理的理论和技术。数字信号处理与模拟信号处理是信号处理的子集。其目的是对真实世界的连续模拟信号进行测量或滤波。 说到数字信号处理，我们就不得不说到滤波器。滤波器是指用来对输入信号进行滤波的硬件和软件。它是一种用来消除干扰杂讯的器件，将输入或输出经过过滤而得到纯净的交流电具有无限持续时间冲激响应,需要用递归模型来实现,其系统函数为: (1.2.1) 设计IIR 滤波器的任务就是寻求一个因果、物理上可实现的系统函数H (z) 使其频率响应满足所希望得到的频域指标即符合给定的通带截止频率、阻带截止、通带衰减和阻带衰减.利用响应不变法设计数字滤波器时可能会导致频域混叠现象为了克服这一问题，需要找到由s平面到z平面的另外的映射关系，这种关系应保证： 1) s平面的整个jΩ轴仅映射为z平面单位圆上的一周； 2) 若G(s)是稳定的，由G(s)映射得到的H(z)也应该是稳定的； 3) 这种映射是可逆的，既能由G(s)得到H(z)，也能由H(z)得到G(s)；G(j0)=1，那么＝1。 如果不考虑频谱混叠现象，用这种方法设计的数字滤波器会很好的重现原模拟滤波器的频响特性，另外数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲击响应波形，时域特性逼近好。但是由于脉冲响应不变法会产生不同程度的频谱混叠失真，只适用于低通、带通滤波器的设计，而不适用于高通、带阻滤波器的设计。 1.2.3 双线性变换法 为了克服脉冲响应不变法会产生频谱混叠失真的这个缺点，人们又找到了双线性变换法来设计IIR数字滤波器。 双线性Z变换的基本思路：首先将整个S平面压缩到一条从－π/Ts变换到2π/Ts 的横带里，然后通过标准的变换关系将横带变换到整个Z平面上去，这样就得到了S平面与Z平面间的一一对应的单值关系。 双线性Z变换满足4个条件的映射关系，其变换公式为 (1.2.3) 图 1.2.3就是表现双线性变换法S平面到Z平面的映射关系的。 图1.2.3 双线性变换法S 平面到Z 平面的映射关系 自然双线性变换法也是有其缺点的，但在这里我们就不详细介绍了。我们要知道只要掌握了这些基本原理，在工程实际中的设计就十分容易了。 1.3 IIR数字滤波器的基本特点 要想知道如何设计IIR数字滤波器，我们就必须了解它的基本特点。而IIR数字滤波器的基本特点有以下几点： ①数字滤波器的系统函数可以写成封闭函数的形式 ②IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。 ③I