《用脉冲响应不变法设计数字的滤波器》由会员分享，可在线阅读，更多相关《用脉冲响应不变法设计数字的滤波器（22页珍藏版）》请在金锄头文库上搜索。
1、word皖西学院数字信号处理课程设计报告题 目 用脉冲响应不变法设计数字滤波器学 院 信息工程学院 专 业 通信工程专业 班 级 * 班 学生某某 陈* 孙* 指导教师 吴* 二 0一二年十二月 / 前言 数字信号处理课程设计是在学生完成数字信号处理和MATLAB结合后的根本实验后开设的，本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的根本实验有本质性的区别。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。 I
2、IR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采用在模拟滤波器设计的根底上进一步变换的方法。其设计方法主要有间接设计法、直接设计法和最大平滑滤波器设计方法。间接法是借助于模拟滤波器的设计方法进展的。其设计步骤是：先设计过度模拟滤波器得到系统函数，然后将其按某种方法转换成数字滤波器的系统函数。这是因为模拟滤波器的设计方法已经成熟，不仅有完整的设计公式，还有完善的图表和曲线供查阅；另外还有一些优良的滤波器可供我们使用。直接法直接在频域或者时域中设计数字滤波器，由于要解联立方程，设计时需要计算机辅助设计。FIR数字滤波器的单位脉冲响应应是有限长序列。它的设计问题实
3、质上是确定能满足要求的转移序列或脉冲响应的常数问题，它不能采用间接法，设计方法主要有窗函数法、频率采样法和等波纹最优逼近法等。目 录第1章绪论22课程设计题目描述与要求33345第2章 MATLAB介绍52.1 MATLAB的简介52.2 MATLAB的优势和特点6第3章 IIR数字滤波器的设计63.1 IIR数字滤波器的设计概述63.2 IIR数字滤波器的设计思想：773.4 巴特沃斯滤波器的设计原理11第4章利用脉冲响应不变法设计数字滤波器的过程16课程设计的解题思路与过程1617第5章总结20参考文献21第1章 绪论电子信息工程专业的培养目标是具备电子技术的根本理论和应用技术，能从事电子
4、、信息、通信、电信等领域的工作，具有高素质、宽口径、创新晋升的专业人才。对本专业学生的培养要进展工程素质培养、拓宽专业口径、注重根底和开展潜力。特别是培养学生的创新能力，以实现技术为主线多进展实验技能的培养。学生通过数字信号处理课程设计这一重要环节，可以将本专业的主干课程-数字信号处理从理论学习到实践应用，对数字信号处理技术有较深的了解，进一步增强学生动手能力和适应实际工作的能力。课程设计题目描述与要求设计题目：脉冲响应不变法设计数字滤波器。设计技术指标：wpp，ap2dB；阻带截止频率wsp ，阻带最小衰减as15dB。设计要求：设计模拟巴特沃斯滤波器并脉冲响应不变法转化成数字滤波器。数字滤
5、波器一词出现在60年代中期。由于电子计算机技术和大规模集成电路的开展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。所谓数字滤波器，是指输入、输出均为数字信号，通过数值运算处理改变输入信号所含频率成分的相比照例，或者滤除某些频率成分的数字器件或程序。因此，数字滤波器的概念和模拟滤波器一样，只是信号的形式和实现滤波方式不同。应用数字滤波器处理模拟信号时，首先须对输入模拟信号进展限带、抽样和模数转换。数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即1/2抽样频率点呈镜像对称。为得到模拟信号，数字滤波器处理的输出数
6、字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、体积小、重量轻、不存在阻抗匹配问题、可程控改变特性或复用、便于集成等优点，可以实现模拟滤波器无法实现的特殊滤波功能。数字滤波器在语言信号处理、图像信号处理、医学生物信号处理以与其他应用领域都得到了广泛应用。 按照不同的分类方法，数字滤波器有许多种类，按照不同的特性具有不同的分类，具体分类如下：按功能分：低通、高通、带通、带阻、全通滤波器。如果我们设计这些理想滤波器是不可能实现的，因为它们的单位脉冲响应是非因果且无限长的，我们只能按照某些准如此设计滤波器，使之在误差容限内逼近理想滤波器，理想滤波器可作为逼近的标准。按实现的网络结构或单位抽
7、样响应长度分类：无限脉冲响应滤波器IIR滤波器、有限脉冲响应滤波器FIR滤波器。它们的系统函数为： 按滤波器对信号的处理作用可分为：选频滤波器和其他滤波器。上述低通、高通、带通、带阻滤波器均属于选频滤波器，其他滤波器有微分器、希尔伯特变换器、频谱校正等滤波器。另外，它还可以被分为线性与非线性、因果与非因果等。其中，线性时不变的数字滤波器是最根本的类型；而由于数字系统可以对延时器加以利用，因此可以引入一定程度的非因果性，获得比传统的因果滤波器更灵活强大的特性；相对于IIR滤波器，FIR滤波器有着易于实现和系统绝对稳定的优势，因此得到广泛的应用；对于时变系统滤波器的研究如此导致了以卡尔曼滤波为代表
8、的自适应滤波理论。常用的数字滤波器一般属于选频滤波器。假设数字滤波器的频率响应函数H(ej)用下式表示： 式中，|H(e)|称为幅频特性函数; ()称为相频特性函数。图1.2表示低通滤波器的幅频特性，p和s分别称为通带边界频率和阻带截止频率。通带频率X围为0|p，在通带中要求(11)|H(ej)|1，阻带频率X围为s|，在阻带中要求|H(ej)|2。从p到s称为过渡带，过渡带上的频响一般是单调下降的。通常，通带内和阻带内允许的衰减一般用分贝数表示，通带内允许的最大衰减用ap表示，阻带内允许的最小衰减用as表示。ap 越小,通带波纹越小，通带逼近误差就越小；as越大,阻带波纹越小，阻带逼近误差就
9、越小；p与s间距越小, 过渡带就越窄。所以低通滤波器的设计指标完全由通带边界频率p、通带最大衰减ap阻带边界频率s和阻带最小衰减as确定。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。这种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法等等。随着MATLAB软件尤其是MATLAB的信号处理工作箱
10、的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。第2章 MATLAB介绍2.1 MATLAB的简介 MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以与交互式程序设计的高科技计算环境。MATLAB的根本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，并且吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进展下载就可以用。2.2 MATLAB的优势和特点1友好的工作平台
11、和编程环境MATLAB由一系列工具组成。2简单易用的程序语言。3出色的图形处理功能4应用广泛的模块集合工具箱。5应用软件开发包括用户界面在开发环境中，使用户更方便地控制多个文件和图形窗口。第3章 IIR数字滤波器的设计3.1 IIR数字滤波器的设计概述IIR滤波器系统函数的极点可以在单位圆内的任何位置，实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，又由于IIR数字滤波器能够保存一些模拟滤波器的优良特性，因此应用很广。设计IIR数字滤波器的方法主要有基于脉冲响应不变法的IIR数字滤波器设计，基于双线性Z变换法的IIR数字滤波器设计，数字高通、带通与带阻IIR滤波器设计，基于MATLAB
12、函数直接设计IIR数字滤波器。3.2 IIR数字滤波器的设计思想：目前，IIR数字滤波器设计最通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表供查询，因此，充分利用这些已有的资源将会给数字滤波器的设计带来很大方便，IIR数字滤波器的设计步骤是：(1)按一定规如此将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标。(2)根据转换后的技术指标设计模拟低通滤波器H(s)。(3)在按一定规如此将H(s)转换为H(z)。假如所设计的数字滤波器是低通的，那么上述设计工作可以完毕。假如所设计的是高通、带通或者带阻滤波器，那么
13、还有步骤(4)将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后按上述步骤(2)设计出模拟低通滤波器H(s)，再由冲击响应不变法或双线性变换将H(s)转换为所需的H(z)设计的关键问题：找出从到的转换方法，主要有冲激响应不变法、阶跃响应不变法、双线性变换法等。下面讨论脉冲响应不变法。设模拟滤波器的系统函数为，相应的单位冲激响应是，=LT。LT代表拉氏变换，对进展等间隔采样，采样间隔为T，得到，将作为数字滤波器的单位脉冲响应，那么数字滤波器的系统函数H(z)便是h(n)的Z变换。因此脉冲响应不变法是一种时域逼近方法，它使h(n)在采样点上等于h(t)。但是，模拟滤波的设计
14、结果是H(s)，所以下面基于脉冲响应不变法的思想，导出直接从H(s)到的转换公式。设模拟滤波器H(s)只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将H(s)用局部分式表示：(3.1) 式中si为H(s)的单阶极点。将Ha(s)进展逆拉氏变换, 得到:(3.2)式中, u(t)是单位阶跃函数。对ha(t)进展等间隔采样，采样间隔为T，得到：(3.3) 对上式进展Z变换，得到数字滤波器的系统函数H(z)，即(3.4) 比照(3.1)和(3.4)式，H(s)的极点映射到z平面的极点为，系数A不变。下面我们分析从模拟滤波器转换到数字滤波器，s平面和z平面之间的映射关系，从而找到这种转换方法的优缺点。这里以理想采样信号作为桥梁，推导其映射关系。设的理想采样信号用表示，即=对进展拉氏变换，得到：3.5 上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式表示： 3.6我们知道模拟信号h(t)的傅里叶变换Ha(j)和其采样信号的傅里叶变换之间的关系满足：将s=j代入上式，得设s=,z=re
