matlab iir带通滤波器,利用MATLAB辅助设计IIR数字带通滤波器

1引言随着计算机和信息科学的极大发展,信号处理已经逐步发展成为一门独立的学科,是信息科学的重要组成部分,在语音处理、图像处理、雷达、航空航天、地质勘探、通信、生物医学等众多领域得到了广泛的应用。信号是数字信号处理领域最基本最重要的概念,信号是信息的载体是信息的物理体现。而数字滤波器作为信号处理一项关键技术是数字信号处理的重要基础,在对信号的过滤、检测、与参数估算等处理过程中,它是使用最为广泛的一种线性系统。数字滤波器按照其冲激响应函数的时域特性,可分为无限长冲激响应(infinitimpulseresponse,IIR)滤波器和有限长冲激响应(finiteim-pulseResponse,FIR)滤波器。在满足相同指标下,IIR滤波器的阶数明显小于FIR,硬件实现容易且大大减少了运算量,在不要求严格线性相位的情况下,IIR滤波器的应用相当广泛。Matlab是一种交互式的以矩阵为基础的软件,它用于科学与工程项目的计算与可视化,它只需要其它编程语言的几分之一的时间即可以解决复杂的数值计算问题。它的强有力也在于那些相对简单的编程功能和提供的非常方便简单的不同学科的工具箱。由于只有低通滤波器的AFD表是现成的,于是为了设计其它滤波器(高通,带通,带阻)我们往往需要一些频率的变换。而本文旨在通过介绍一些常用的方法来解决IIR带通滤波器的设计问题。2IIR滤波器基本概念一个N阶递归型数字滤波器的差分方程为My(n)=i=0bix(n-r)-Ni=1aiy(n-i)IIR滤波器的系统函数H(Z)=Mr=0brz-r1+NK=1akz-k从以上的系统函数可知,设计IIR滤波器的任务就是通过计算寻求一个因果、物理上可实现的系统函数H(Z),使其频率响应H(ejw)满足所希望得到的频域指标,即符合给定的通带截止频率、阻带截止、通带衰减和阻带衰减。3数字滤波器的设计过程数字滤波器的一般设计过程可分为以下四个步骤:1)按照实际的需要,确定滤波器的性能要求或指标;2)用一个因果稳定的离散线性时不变系统,去逼近这一性能指标;3)用有限精度的运算实现所设计的系统;4)通过模拟,验证所设计的系统是否符合性能要求4应用实例1)用切比雪夫函数直接设计IIR数字带通滤波器在工程中模拟滤波器可以设计出相应的数字低通滤波器,此方法也能推广到设计其它类型的选频滤波器。通过对低通滤波器的频率轴进行变换可以得到其它类型的选频器。而matlab提供了能在s平面上实现频带转换的函数,我们下面首先说明下z平面映射的用法,接着再通过一个切比雪夫滤波器的应用实例来印证这种方法。通过给定的低通滤波器原型,并设H(z)是期望的选频数字滤波器,通过一系列的变换得到:nZ-1=G(z-1)=k=1z-1-k1-kz-1限于篇幅下表只给出了原型低通到带通的转换,如需对高通和带阻有所研究的读者可以自行查资料。表1数字滤波器的频带变换变换类型变换参数带通z-1z-2-1z-1+22z-2-1z-1+1w1=低截止频率wu高截止频率1=-2k/(k+1)2=(k-1)/(k+1)k=cos[(wu+wl)/2]cos[(wu-wl)/2]k=coswu-wetanwcCheby1函数设计IIR带通数字滤波器的调用方式[b,a]=cheby1(N,r,Wn)其中Wn为二元向量即Wn=[w1,w2],则此函数返回的是一个2n阶的带通切比雪夫I型带通滤波器的设计结果,其通带为[w1,w2]例1设计一个数字带通滤波器,要求如下:设计一个通带为0.6*pi到0.8*pi的8阶滤波器,其采样频率为1000Hz,其中那个滤波器在通带