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