无限长单位脉冲响应(IIR)滤波器设计 (1)

实验四 IIR 数字滤波器的设计及网络结构一、实验目的1．了解IIR 数字滤波器的网络结构。

2.掌握模拟滤波器、IIR 数字滤波器的设计原理和步骤。

3．学习编写数字滤波器的设计程序的方法。

二、实验内容数字滤波器：是数字信号处理技术的重要内容。

它的主要功能是对数字信号进行处理，保留数字信号中的有用成分，去除信号中的无用成分。

1．数字滤波器的分类滤波器的种类很多，分类方法也不同。

（1）按处理的信号划分：模拟滤波器、数字滤波器 （2）按频域特性划分；低通、高通、带通、带阻。

（3）按时域特性划分：FIR 、IIR2．IIR 数字滤波器的传递函数及特点数字滤波器是具有一定传输特性的数字信号处理装置。

它的输入和输出均为离散的数字信号，借助数字器件或一定的数值计算方法，对输入信号进行处理，改变输入信号的波形或频谱，达到保留信号中有用成分去除无用成分的目的。

如果加上A/D 、D/A 转换，则可以用于处理模拟信号。

设IIR 滤波器的输入序列为x(n)，则IIR 滤波器的输入序列x(n)与输出序列y(n)之间的关系可以用下面的方程式表示：1()()()M Ni j i j y n b x n i a y n j ===-+-∑∑（5-1）其中，j a 和i b 是滤波器的系数，其中j a 中至少有一个非零。

与之相对应的差分方程为：10111....()()()1....MM NN b b z b z Y z H Z X z a z a z ----++==++ （5-2）由传递函数可以发现无限长单位冲激响应滤波器有如下特点： （1） 单位冲激响应h(n)是无限长的。

（2） 系统传递函数H(z)在有限z 平面上有极点存在。

（3） 结构上存在着输出到输入的反馈，也就是结构上是递归型的。

3．IIR 滤波器的结构IIR 滤波器包括直接型、级联型和并联型三种结构：① 直接型：优点是简单、直观。

但由于系数bm 、a k 与零、极点对应关系不明显，一个bm 或a k 的改变会影响H(z)所有零点或极点的分布，所以一方面，bm 、a k 对滤波器性能的控制关系不直接，调整困难；另一方面，零、极点分布对系数变化的灵敏度高，对有限字长效应敏感，易引起不稳定现象和较大误差。

FIR和IIR滤波器设计滤波器是信号处理中常用的工具，用于去除信号中的噪声、增强或抑制特定频率成分等。

FIR（有限脉冲响应）和IIR（无限脉冲响应）是两种常见的滤波器设计方法。

FIR滤波器是一种线性相位的滤波器，其脉冲响应是有限长度的，因此被称为有限脉冲响应。

它的频率响应是通过一个线性组合的单位样本响应来实现的。

在设计FIR滤波器时，可以通过窗函数法或频率采样法来选择滤波器的系数。

窗函数法适用于要求较为简单的滤波器，而频率采样法适用于要求较高的滤波器。

窗函数法是一种基于原始滤波器响应的方法。

它通过将滤波器响应乘以一个窗函数，从而使得脉冲响应在时间上截断。

常用的窗函数有矩形窗、汉明窗、布莱克曼窗等。

通过选择不同窗函数可以得到不同的滤波器特性，如频带宽度、峰值纹波等。

频率采样法是一种通过等间隔采样得到频率响应的方法。

首先确定滤波器的截止频率和带宽，然后选择一组频率点进行采样。

根据采样得到的频率响应，可以通过逆傅里叶变换得到滤波器的脉冲响应，进而得到滤波器的系数。

频率采样法可以灵活地选择频率点，从而得到更精确的滤波器特性。

与FIR滤波器不同，IIR滤波器的脉冲响应是无限长度的，因此被称为无限脉冲响应。

IIR滤波器的频率响应是通过递归方式的单位样本响应来实现的。

在设计IIR滤波器时，可以通过模拟滤波器的方法来选择滤波器的结构和参数。

常用的模拟滤波器有巴特沃斯滤波器、切比雪夫滤波器等。

巴特沃斯滤波器是一种最优近似设计的滤波器，其特点是在通带和阻带中都具有等级衰减。

切比雪夫滤波器是一种在通带和阻带中都具有等级衰减，同时具有较窄过渡带的滤波器。

这两种滤波器的设计方法都是基于频率变换的思想，首先将模拟滤波器的频率响应映射到数字滤波器上，然后利用一定的优化算法来得到滤波器的参数。

FIR和IIR滤波器在滤波器设计中有不同的特点和适用范围。

FIR滤波器具有线性相位特性，因此适用于对信号的相位要求较高的应用，如音频处理、图像处理等。

实验四无限长单位脉冲响应滤波器设计一、实验目的1.掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。

2.观察双线性变换及脉冲响应不变法设计的滤波器的频率特性，了解双线性变换法及脉冲响应不变法的特点。

3.熟悉巴特沃思滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。

二、实验原理(1)实验中有关变量的定义：fc通带边界频率，fr阻带边界频率，tao通带波动，at 最小阻带衰减，fs采样频率，t采样周期。

(2)设计一个数字滤波器一般包括以下两步：a.按照任务要求，确定滤波器性能指标b.用一个因果稳定的离散时不变系统的系统函数去逼近这一性能要求(3)数字滤波器的实现：对于IIR滤波器，其逼近问题就是寻找滤波器的各项系数，使其系统函数逼近一个所要求的特性。

先设计一个合适的模拟滤波器，然后变换成满足约定指标的数字滤波器。

用双线形变换法设计IIR数字滤波器的过程：a.将设计性能指标中的关键频率点进行“预畸”b.利用“预畸”得到的频率点设计一个模拟滤波器。

c.双线形变换，确定系统函数三、实验内容1、设计一切比雪夫高通滤波器，性能指标如下：通带边界频率f c=0.4kHz，通带波动δ=0.5dB，阻带边界频率f r=0.3kHz，阻带最小衰减At=20dB，采样频率f s=1000Hz，观察其通带波动和阻带衰减是否满足要求。

（绘制对数幅度谱）2、设计一巴特沃思低通滤波器，性能指标如下：通带边界频率f c=0.4kHz，通带波动δ=1dB，阻带边界频率f r=0.6kHz，阻带最小衰减At=40dB，采样频率f s=2000Hz，分别用脉冲响应不变法和双线性变换法进行设计，比较两种方法的优缺点。

（绘制线性幅度谱）3、用双线性变换法设计巴特沃思、切比雪夫和椭圆低通滤波器，其性能指标如下：通带边界频率f c=1.8kHz，通带波动δ≤1dB，阻带边界频率f r=2.6kHz，阻带最小衰减A t≥50dB，采样频率f s=8kHz。

无限长单位脉冲响应IIR 滤波器的设计一．设计目的1．掌握数字滤波器的设计过程；2．了解IIR 的原理和特性；3．熟悉设计IIR 数字滤波器的原理和方法；4．学习IIR 滤波器的DSP 实现原理；5．学习使用CCS 的波形观察窗口观察输入/输出信号波形和频谱变化情况。

二．设计内容1．通过MATLAB 来设计一个低通滤波器，对它进行模拟仿真确定IIR 滤波器系数。

2．用DSP 汇编语言及C 语言进行编程，实现IIR 运算，对产生的合成信号，滤除信号中高频成分，观察滤波前后的波形变化。

三．设计原理IIR 滤波器与FIR 滤波器相比具有相位特性差的特点，但它的结构简单，运算量小，具有经济高效的特点，并且可以用较少的阶数获得很高的选择性，因此也得到了广泛的应用。

IIR 数字滤波器系统的传递函数为：H(z)= NN N N z a z a z b z b b z X z Y ----+⋯⋯+++⋯⋯++=111101)()( 它具有N 个极点和N 个零点，如果任何一个极点在单位圆外，则系统不稳定。

如果系数a j （j=1，… ，N ）全部为0，则滤波器变成非递归的FIR滤波器，系统总是稳定的。

对于IIR滤波器，有系数量化敏感的缺点。

由于系统对序列施加的算法是由加法、乘法和延时的基本运算的组合，所以可以用不同结构的数字滤波器来实现而不影响系统总的传输函数。

四．MATLAB设计IIR滤波器的方法我们所用滤波器设计方法为巴特沃夫Butterworth滤波器设计和切比雪夫Chebyshev滤波器设计。

MATLAB的butter函数可以设计低通、带通、高通和带阻数字滤波器，其特征可以使通带内的幅度响应最大限度的平坦，但会损失截止频率处的下降斜度，使幅度响应衰减较慢，因此butter函数主要用于设计通带平坦的数字滤波器。

如果期望幅度响应下降斜度大，衰减快，可以使用Elliptic(椭圆)或Chebyshev(切比雪夫)滤波器。

基于MATLAB的IIR数字滤波器设计与仿真一、概述在现代数字信号处理领域中，数字滤波器扮演着至关重要的角色。

其通过对输入信号的特定频率成分进行增强或抑制，实现对信号的有效处理。

无限脉冲响应（IIR）数字滤波器因其设计灵活、实现简单且性能优良等特点，得到了广泛的应用。

本文旨在基于MATLAB平台，对IIR数字滤波器的设计与仿真进行深入研究，以期为相关领域的研究与应用提供有益的参考。

IIR数字滤波器具有无限长的单位脉冲响应，这使得其在处理信号时能够展现出优秀的性能。

与有限脉冲响应（FIR）滤波器相比，IIR滤波器在实现相同性能时所需的阶数更低，从而减少了计算复杂度和存储空间。

在需要对信号进行高效处理的场合，IIR滤波器具有显著的优势。

MATLAB作为一款功能强大的数学软件，提供了丰富的函数和工具箱，使得数字滤波器的设计与仿真变得简单而高效。

通过MATLAB，我们可以方便地实现IIR滤波器的设计、分析和优化，从而满足不同应用场景的需求。

本文将首先介绍IIR数字滤波器的基本原理和特性，然后详细阐述基于MATLAB的IIR数字滤波器的设计方法和步骤。

接着，我们将通过仿真实验验证所设计滤波器的性能，并对其结果进行分析和讨论。

本文将总结IIR数字滤波器设计与仿真的关键技术和注意事项，为相关领域的研究人员和工程师提供有益的参考和启示。

1. IIR数字滤波器概述IIR（Infinite Impulse Response）数字滤波器是数字信号处理中常用的一类滤波器，它基于差分方程实现信号的滤波处理。

与FIR （Finite Impulse Response）滤波器不同，IIR滤波器具有无限长的单位脉冲响应，这意味着其输出不仅与当前和过去的输入信号有关，还与过去的输出信号有关。

这种特性使得IIR滤波器在实现相同的滤波效果时，通常具有更低的计算复杂度，从而提高了处理效率。

IIR滤波器的设计灵活多样，可以根据不同的需求实现低通、高通、带通和带阻等多种滤波功能。

实验四 无限冲激响应（IIR ）数字滤波器设计一、实验目的1．熟悉用双线性变换法设计IIR 数字滤波器的原理和方法；2．了解用脉冲响应不变法设计IIR 数字滤波器的原理和方法；3．掌握双线性变换及脉冲响应不变法设计的滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点；4．掌握数字滤波器的计算机仿真方法；二、实验原理介绍IIR 数字滤波器的系统函数为1z -的有理分式： 1011()1N kk Nk k b z H z a z -=-==+∑∑ 设计IIR 滤波器的系统函数，就是要确定()H z 的阶数N 及分子分母多项式的系数k a 和k b ，使其()()j j z e H e H z ωω==满足指定的频率特性。

由于模拟滤波器的设计有许多简单而严谨的设计公式和大量的图表可以利用，因此IIR 滤波器设计的方法之一是：先设计一个合适的模拟滤波器，然后将模拟滤波器通过适当的变换转换成满足给定指标的数字滤波器。

1、Butterworth 模拟低通滤波器221()1a N c H j Ω=⎛⎫Ω+ ⎪Ω⎝⎭幅度平方函数：其中，N 为滤波器的阶数，c Ω为通带截止频率。

2．Chebyshev 模拟低通滤波器 2221()1()a N c H j C εΩ=Ω+Ω幅度平方函数：3、脉冲响应不变法原理 用数字滤波器的单位脉冲响应序列h(n)逼近模拟滤波器的冲激响应()a h t ，让h(n)正好等于()a h t 的采样值，即：()()a h n h nT =其中，T 为采样间隔。

如果以()a H s 和H(z)分别表示()a h t 的拉氏变换及h(n)的Z 变换，则：12ˆ()()sT a a z e k H z H s H s j k T T π∞==-∞⎛⎫==- ⎪⎝⎭∑4、双线性变换法原理双线性变换法是通过两次映射采用非线性频率压缩的方法，将整个频率轴上的频率范围压缩到±π/T 之间，再用sTz e =转换到z 平面上，从而使数字滤波器的频率响应与模拟滤波器的频率响应相似。

脉冲响应不变法设计iir数字滤波器以脉冲响应不变法设计IIR数字滤波器引言：数字滤波器是数字信号处理中的重要组成部分，用于处理和改变数字信号的频率特性。

脉冲响应不变法（Impulse Invariance Method）是一种常用的IIR数字滤波器设计方法，其基本原理是通过将连续时间域中的模拟滤波器的脉冲响应与采样脉冲响应进行匹配，从而实现滤波器的设计。

一、脉冲响应不变法基本原理脉冲响应不变法的基本原理是将连续时间域中的模拟滤波器的脉冲响应与离散时间域中的数字滤波器的脉冲响应进行匹配。

在这种方法中，首先需要确定模拟滤波器的脉冲响应，然后通过采样得到数字滤波器的脉冲响应，最后将其离散化得到数字滤波器的差分方程。

二、脉冲响应不变法的设计步骤1. 确定模拟滤波器的脉冲响应：选择适当的模拟滤波器类型，并设计其频率响应。

根据滤波器的阶数和截止频率，确定模拟滤波器的差分方程。

2. 采样得到数字滤波器的脉冲响应：通过将连续时间域中的模拟滤波器的脉冲响应与采样脉冲进行卷积，得到数字滤波器的脉冲响应。

3. 离散化得到数字滤波器的差分方程：将数字滤波器的脉冲响应离散化，得到数字滤波器的差分方程。

根据差分方程，可以计算数字滤波器的各个系数。

三、脉冲响应不变法的优缺点脉冲响应不变法具有以下优点：1. 设计方法简单：通过匹配模拟滤波器和数字滤波器的脉冲响应，可以直接得到数字滤波器的差分方程，设计方法相对简单。

2. 精度较高：脉冲响应不变法可以保持模拟滤波器的频率响应特性，因此可以实现较高的滤波器精度。

3. 适用范围广：脉冲响应不变法适用于各种模拟滤波器类型和滤波器规格的设计。

然而，脉冲响应不变法也存在一些缺点：1. 频率响应失真：由于采样过程中的截断和抽样误差，脉冲响应不变法可能导致数字滤波器的频率响应失真。

2. 高阶滤波器设计困难：对于高阶滤波器的设计，脉冲响应不变法可能会导致数字滤波器的稳定性问题和数值计算问题。

四、脉冲响应不变法的应用领域脉冲响应不变法广泛应用于数字信号处理领域，特别是在音频信号处理、图像处理和通信系统中的滤波器设计中。

第五章 数字滤波器5-3 无限长单位脉冲响应数字滤波器设计一、概述数字滤波器从实现的网络结构或者从单位脉冲响应分类，可以分成无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。

它们的系统函数分别为：理想低通、高通、带通、带阻滤波器幅度特性2 数字滤波器的技术要求我们通常用的数字滤波器一般属于选频滤波器。

假设数字滤波器的传输函数H(e j ω)用下式表示：通带内和阻带内允许的衰减一般用dB 数表示，通带内允许的最大衰减用αp 表示，阻带内允许的最小衰减用αs 表示，αp 和αs 分别定义为：011()1()()Mrr r Nk kk N nn b z H z a z H z h n z -=-=--==+=∑∑∑)(e j ωH )(e j ωH )(e j ωH )(e j ωH 0低通0高通0带通0带阻ωωωωπ-π2-π2-π2-π-π-π-π-ππππ2π2π2π2()()()j j j H e H e e ωωωΩ=00()20lg ()()20lg()psj p j j s j H e dBH e H e dBH e ωωαα==如将|H(ej0)|归一化为1，则表示成：3. 数字滤波器设计方法概述IIR 滤波器和FIR 滤波器的设计方法是很不相同的。

IIR 滤波器设计方法有两类，经常用的一类设计方法是借助于模拟滤波器的设计方法进行的。

其设计步骤是：先设计模拟滤波器得到传输函数Ha(s)，然后将Ha(s)按某种方法转换成数字滤波器的系统函数H(z)。

二、模拟滤波器的设计模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。

图5-21理想滤波器的幅频特性1.模拟低通滤波器的设计指标及逼近方法模拟低通滤波器的设计指标有αp, Ωp,αs 和Ωs 。

第5章 无限长单位脉冲响应（IIR ） 数字滤波器的设计方法5.1 基本概念5.1.1 选频滤波器的分类可分为低通、高通、带通、带阻和全通满足奈奎斯特采样定理时，信号的频率特性只能限带于|ω|<π的范围。

5.1.2 滤波器的技术指标在通带内，幅度响应以最大误差±δ1逼近于1，即1111 δδωω+<≤-≤）（jw P e H在阻带内，幅度响应以误差小于δ2而逼近于零，即2 ,δπωω≤≤≤）（jw S e H1)()ωj e H (π1 δ阻带过渡带ωπ2 δ通带11δ+11δ-sw p wN 为奇数时实轴上有极点，N 为偶数时实轴上没有极点。

要称为稳定的滤波器)s (H a 表示为：∏=-Ω=Nk k N Ca )s s ()s (H 1方法2:,s s p A A P 、、由模拟ΩΩc ,Ω→N c ,Ω→N⎥⎥⎦⎤⎢⎢⎣⎡ΩΩ+=Ω-Ω=ΩNa P j H A 2c p 2p )(11lg 10)(lg 10 －＝在 ⎥⎦⎤⎢⎣⎡ΩΩ+=Ω-Ω=ΩN a P j H A 2c S 2S )(11lg 10)(lg 10 －＝在 ⎥⎥⎦⎤⎢⎢⎣⎡ΩΩ--=)(lg 2])110()110(lg[s p 1010sp A A N NA sc s 210110-Ω=Ω（通带指标改善）或 NA Pc P 210110-Ω=Ω（阻带指标改善）方法3: )()(s H s H a aN →NN N aN ss a s a s a s H ++++=--1122111)( caN a s s s H s H Ω'=→)()( 例 1 导出三阶巴特沃思模拟低通滤波器的系统函数， 设Ωc ＝2 rad/s 。

解 幅度平方函数是62)2/(11|)(|Ω+=Ωj H 令Ω2=-s 2即s =j Ω，则有)2/(11)()(66s s H s H a a -=-π⎥⎦⎤⎢⎣⎡-+=612212k j k es k =1, 2, …, 6会产生频率混叠，适合低通、带通滤波器的设计，不适合高通、带阻滤波器的设计。