数字信号处理(DSP)课程设计—利用Matlab设计IIR数字滤波器.docVIP

数字信号处理(DSP)课程设计—利用Matlab设计IIR数字滤波器 数字信号处理 课 程 设 计 报 告 姓名：一.蒲钇霖 二.舒琳 学号：一.201021030619 二.201021030414 学院：微固 利用Matlab设计IIR数字滤波器 一、前言 数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，这也就促进了数字信号处理这门技术的飞速发展，使其不仅自成一门学科，更是以不同形式影响和渗透到其它学科。如今，数字信号处理技术与国民经济息息相关，与国家建设紧密相连，影响着并改变着我们的生产、生活方式，因此它受到了社会的普遍关注。 而在数字信号处理当中，数字滤波器又占有相当重要的位置。因为任何需要我们处理的信号都包含有噪声，而滤波是去除噪声的基本手段。数字滤波器简单的说就是一种用来过滤时间离散信号的线性时不变系统，通过对采样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可以分为两类：无限冲激响应(IIR)滤波器和有限冲激响应（FIR）滤波器。数字滤波本质上是一种运算过程，输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器的本质上是一个完成特定运算的数字计算过程。描述离散系统输入和输出关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。 二、数字滤波器设计原理 在数字滤波过程中，我们主要讨论离散时间序列。设输入序列为x(n)，滤波器对单位采样序列δ(n)的响应为h(n)。因δ(n)在时域离散信号和系统中所起的作用相当于单位冲激函数在时域连续信号和系统中所起的作用。所以数字滤波器的序列y(n)就是以上两个序 列的卷积，即： ?? y(n)? k????h(k)x(n?k) (1) 同样，两个序列卷积的z变换等于各自z变换的乘积，即： Y(z)?H(z)X(z) (2) 用z?ej?代入上式，则可得到时域离散系统的频谱特性： Y(ej?)?H(ej?)X(ej?) (3) 式中Y(ej?)和X(ej?)分别是数字滤波器的输出序列和输入序列的频谱特性，而H(ej?)则是单位采样序列h(n)的频特性谱??又称为数字滤波器的频域响应。输入序列X(ej?)经过滤波后变为X(ej?)H(ej?)。因此，只需要按照输入信号频谱X(ej?)的特点和我们处理信号的目的，选择适当的频域响应H(ej?)，使得滤波后的输出信号频谱Y(ej?)满足我们设计的要求。 三、课程设计目的与内容 通过完成本次课程设计，加深我们对于在数字信号处理这门学科的印象，巩固我们在课堂上学习到的理论知识和实验方法，增强我们将理论知识化为实践技巧的能力，并且让我们进一步熟悉数字信号处理中常用软件Matlab的使用方法。在这个过程中同时培养我们发现问题、分析问题以及解决问题的能力。 在这次课程设计当中，我们需要分析和研究采用不同设计方法时，滤波器最低阶数与设计指标之间的关系，通过理论分析方式推导建立起相应的简单数学模型。在此基础上，自选设计指标，分别采用双线性变换法和冲激响应不变法进行设计Butterworth低通滤波器，然后通过仿真检验所建模型的正确性和局限性。 四、具体设计步骤 1. 推导滤波器最小阶数与设计指标的关系 滤波器阶数就是其系统函数的极点个数。为了避免滤波器的复杂程度与我们的要求不匹配，造成不必要的成本浪费，我们在滤波器设计前先需要确定其合适的阶数，即满足设计要求的最小阶数。下面就分别以Butterworth滤波器和Chebyshev1型滤波器为例，推导其阶数的数学模型。 ①对于Butterworth滤波器，若给出滤波器的设计指标为：通带截止频率?p，阻带截止频率?s，通带最大纹波Rp(dB)，阻带最小纹波Rs(dB)。因为滤波器幅频特性为： H(j?? （4） 其中?c为3dB截止频率，N为滤波器阶数。所以当???p以及 ???s时，可得到： H(j?p? H(j?s? (5) 然后由Rp(dB)与Rs(dB)的关系式可得到： 2?11 R?10lg?H(j?)?p?pRp2 H(j?p??1010 ? 211?R?10lg?H(j?)?sRs2?s H(j?s?101