IIR切比雪夫低通数字滤波器

./数字信号处理课程设计报告课题名称：系别：学号：姓名：班级：指导教师：.目录1.1题目41.2目的41.3容及要求41.4论文格式42.1数字滤波器的概念62.2数字滤波器的分类62.3数字滤波器的设计要求82.4IIR数字滤波器的设计9IIR数字滤波器的设计步骤9双线性变换法设计IIR数字滤波器103.1设计步骤153.2MATLAB程序15.摘要[][]在现代通信系统中,由于信号中经常混有各种复杂成分,所以很多信号分析都是基于滤波器而进行的,而数字滤波器是通过数值运算实现滤波,具有处理精度高、稳定、灵活、不存在阻抗匹配问题,可以实现模拟滤波器无法实现的特殊滤波功能。数字滤波器根据其冲激响应函数的时域特性,可分为两种,即无限长冲激响应<IIR>数字滤波器和有限长冲激响应<FIR>数字滤波器。实现IIR滤波器的阶次较低,所用的存储单元较少,效率高,精度高,而且能够保留一些模拟滤波器的优良特性,因此使用很广。Matlab软件以矩阵运算为基础,把计算、可视化及程序设计有机融合到交互式工作环境中,并且为数字滤波的研究和使用提供了一个直观、高效、便捷的利器。尤其是Matlab中的信号处理工具箱使各个领域的研究人员可以直观方便地进行科学研究和工程使用。本文首先介绍了数字滤波器的概念,分类以及设计要求。接着利用MATLAB函数语言编程,用信号处理图形界面FDATool来设计滤波器以及Sptool界面设计的方法,并用FDATool模拟IIR数字滤波器处理信号。重点设计ChebyshevI型和ChebyshevII型数字低通滤波器,并介绍最优化设计。关键字:IIR;滤波器;FDATool;Sptool;SimulinkABSTRACTInmoderncommunicationsystems,Becauseoftenmixedwithvarioussignalcomplexcomponents,Somanysignalanalysisisbasedonfilters,andthedigitalfilterisrealizedthroughnumericalcomputation,digitalfiltersfilterwithhighprecision,stabilityandflexibility,don'texist,canrealizetheimpedancematchingsimulatingthespecialfiltercannotachievefilterfunction.Digitalfilteraccordingtoitsimpulseresponsefunctionandcharacteristicsofthetimecanbedividedintotwokinds,namelytheinfiniteimpulseresponse<IIR>digitalfilterandfiniteimpulseresponse<FIRdigitalfilters>.TheorderofrealizingIIRfilterisused,lowandhighefficiencylessstorageunit,highprecision,andcankeepsomesimulationcharacteristicsoffilter,soitiswidelyused.Matlabsoftwarebasedonmatrixcomputation,thecalculation,visualizationandprogramdesignoforganicintegrationtointeractiveenvironmentfordigitalfilter,andtheresearchandapplicationofprovidesanintuitive,efficientandconvenienttool.EspeciallyintheMatlabsignalprocessingtoallareasofresearchtoolboxpersonnelcaneasilyforscientificresearchandengineeringapplication.Thispaperintroducestheconceptofdigitalfilter,classificationanddesignrequirements.ThenusingMATLABlanguageprogramming,withfunctionsofsignalprocessingFDAToolgraphicalinterfacedesignofinterfacedesignandSptoolfilter,andFDAToolanalogsignalprocessingIIRdigitalfilter.KeydesignChebyshevtypeIandIIdigitalChebyshevlowpassfilter,andintroducesoptimizationdesign.Keywords:IIR;Filter;FDATool;Sptool;Simulink任务书题目设计原型为切比雪夫型的数字IIR低通滤波器。目的学会MATLAB的使用,掌握MATLAB的程序设计方法。掌握数字信号处理的基本概念、基本理论和基本方法。掌握MATLAB设计IIR数字滤波器的方法。学会用MATLAB对信号进行分析和处理。容及要求设计原型为切比雪夫型的数字IIR低通滤波器〔本文用双线性变换法,原型为切比雪夫I型。低通数字滤波器的技术指标：通带截止频率,通带衰减,阻带截止频率,阻带衰减。求出所设计滤波器的Z变换。用MATLAB画出幅频特性图。论文格式任务书及摘要。目录,中间打点号。章节,如：第1章,每一章都要另起一页,居中；二级目录要打头写,如：1.1；每个图都要有图标和名称,如：图1-1++++++。课程设计总结,主要包括课程设计意义、感受、遇到的问题及解决办法。参考文献,格式如：[1]阎石.数字电子技术基础.:高等教育,1998.数字滤波器的设计数字滤波器的概念滤波器是指用来对输入信号进行滤波的硬件和软件。数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波器可以理解为是一个计算程序或算法,将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中,使信号按预定的形式变化。数字滤波实质上是一种运算过程,实现对信号的运算处理。数字滤波器和模拟滤波器相比,因为信号的形式和实现滤波的方法不同,数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。输入数字信号〔数字序列通过特定的运算转变为输出的数字序列,因此,数字滤波器本质上是一个完成特定运算的数字计算过程,也可以理解为是一台计算机。描述离散系统输出和输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则,使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性：其中、分别是数字滤波器的输出序列和输入序列的频域特性〔或称为频谱特性,是数字滤波器的单位取样响应的频谱,又称为数字滤波器的频域响应。输入序列的频谱经过滤波后,因此,只要按照输入信号频谱的特点和处理信号的目的,适当选择,使得滤波后的满足设计的要求,这就是数字滤波器的滤波原理。数字滤波器的分类按照不同的分类方法,数字滤波器有许多种类,但总起来可以分成两大类：经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的频率成分和希望滤除的频率成分占有不同的频带,通过一个合适的选频滤波器滤除干扰,得到纯净信号,达到滤波的目的。但是,如果信号和干扰的频谱相互重叠,则经典滤波器不能有效地滤除干扰,最大限度地恢复信号,这时就需要现代滤波器,例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。现代滤波器是根据随机信号的一些统计特性,在某种最佳准则下,最大限度地抑制干扰,同时最大限度地回复信号,从而达到最佳滤波的目的。经典数字滤波器从滤波特性上分类,可以分为：低通滤波器、高通滤波器、带通滤波器、带阻滤波器。图STYLEREF2\s2.2.SEQ图\*ARABIC\s21各种理想滤波器的幅频特性数字滤波器根据其冲激响应函数的时域特性,可分为两种,即无限长冲激响应<IIR>数字滤波器和有限长冲激响应<FIR>数字滤波器。IIR数字滤波器的特征是,具有无限持续时间冲激响应,需要用递归模型来实现。其差分方程为：〔式2-2-1系统函数为：〔式2-2-2设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H<z>,使其频率响应H<z>满足所希望得到的频域指标,即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。数字滤波器的设计要求滤波器的指标常常在频域给出。数字滤波器的频响特性函数一般为复函数,所以通常表示为：〔式2-3-1其中,称为幅频特性函数,Φ<w>称为相频特性函数。幅频特性表示信号通过该滤波器后各频率成分的衰减情况,而相频特性反映各频率通过滤波器后在时间上的延时情况。一般IIR数字滤波器,通常只用幅频响应函数来描述设计指标,相频特性一般不作要求。IIR滤波器指标参数如下图所示。图中,和分别为通带边界频率和阻带边界频率；δ1和δ2分别为通带波纹和阻带波纹；允许的衰减一般用dB数表示,通带所允许的最大衰减〔dB和阻带允许的最小衰减〔dB分别为和表示：〔式2-3-2〔式2-3-3一般要求：当时,；图STYLEREF2\s2.3.图STYLEREF2\s2.3.SEQ图\*ARABIC\s21低通滤波器的技术要求IIR数字滤波器的设计[][]丁玉美,高西全.数字信号处理[M].XX：XX电子科技大学出版社,2011.IIR数字滤波器是一种离散时间系统,其系统函数为假设M≤N,当M＞N时,系统函数可以看作一个IIR的子系统和一个<M-N>的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和,它是数学上的一种逼近问题,即在规定意义上〔通常采用最小均方误差准则去逼近系统的特性。如果在S平面上去逼近,就得到模拟滤波器；如果在z平面上去逼近,就得到数字滤波器。IIR数字滤波器的设计步骤IIR数字滤波器的设计一般有两种方法：一个是借助模拟滤波器的设计方法进行。其设计步骤是,先设计模拟滤波器,再按照某种方法转换成数字滤波器。这种方法比较容易一些,因为模拟滤波器的设计方法已经非常成熟,不仅有完整的设计公式,还有完善的图表供查阅；另外一种直接在频率或者时域进行,由于需要解联立方程,设计时需要计算机做辅助设计。为了保证转换后的稳定且满足技术指标要求,对转换关系提出两点要求:因果稳定的模拟滤波器转换成数字滤波器,仍是因果稳定的。数字滤波器的频率相应模仿模拟滤波器的频响特性,s平面的虚轴映射为z平面的单位圆,相应的频率之间呈线性关系。利用模拟滤波器成熟的理论设计IIR数字滤波器的过程是：确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减、阻带截止频率、阻带最小衰减。将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。用所选的转换方法,将模拟滤波器转换成数字低通滤波器系统函数。IIR数字滤波器的设计流程图如下：图STYLEREF2\s2.4.SEQ图_\*ARABIC\s21IIR数字滤波器的设计步骤流程图双线性变换法设计IIR数字滤波器如前所述,IIR数字滤波器的设计方法有间接法和直接法。间接法有两种：脉冲响应不变法和双线性变换法。这两种方法已经是很成熟的模拟滤波器设计方法,此处仅讨论双线性变换法。设计原理[][]脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点,可以采用非线性频率压缩方法,将整个频率轴上的频率围压缩到之间,再用转换到Z平面上。也就是说,第一步先将整个S平面压缩映射到S1平面的一条横带里；第二步再通过标准变换关系将此横带变换到整个Z平面上去。这样就使S平面和Z平面建立了一一对应的单值关系,消除了多值变换性,也就消除了频谱混叠现象,映射关系如REF_Ref344331489\h图2.4.2所示。图STYLEREF2\s2.4.SEQ图\*ARABIC\s22双线性变换的映射关系为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的段上,可以通过以下的正切变换实现〔式2-6式中,T仍是采样间隔。当Ω1由经过0变化到时,Ω由-∞经过0变化到+∞,也即映射了整个jΩ轴。将式〔2-6写成〔式2-7将此关系分析延拓到整个S平面和S1平面,令jΩ=s,jΩ1=s1,则得〔式2-8再将S1平面通过以下标准变换关系映射到Z平面〔式2-9从而得到S平面和Z平面的单值映射关系为：〔式2-10〔式2-11式〔2-10和式〔2-11是S平面和Z平面之间的单值映射关系,这种变换都是两个线性函数之比,因此称为双线性变换式〔2-6和式〔2-10的双线性变换符合映射变换应满足的两点要求。首先,把,可得〔式2-12即S平面的虚轴映射到Z平面的单位圆。其次,将代入式〔2-12,得〔式2-13因此〔式2-14由此看出,当时,;当时,.也就是说,S平面的左半平面映射到Z平面的单位圆,S平面的右半平面映射到Z平面的单位圆外,S平面的虚轴映射到Z平面的单位圆上。因此,稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。双线性变换法优缺点双线性变换法和脉冲响应不变法相比,其主要的优点是避免了频率响应的混叠现象。这是因为S平面和Z平面是单值的一一对应关系。S平面整个jΩ轴单值地对应于Z平面单位圆一周,即频率轴是单值变换关系。这个关系如式〔2-10所示,重写如下：〔式2-15上式表明,S平面上Ω和Z平面的ω成非线性的正切关系,如REF_Ref344334040\h图2.4.2。由REF_Ref344334040\h图2.4.2看出,在零频率附近,模拟角频率Ω和数字频率ω之间的变换关系接近于线性关系；但当Ω进一步增加时,ω增长得越来越慢,最后当Ω→∞时,ω终止在折叠频率ω=π处,因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象,从而消除了频率混叠现象。图STYLEREF2\s2.4.SEQ图\*ARABIC\s23双线性变换法的频率变换关系但是双线性变换的这个特点是靠频率的严重非线性关系而得到的,如式〔2-12及REF_Ref344334040\h图2.4.3所示。由于这种频率之间的非线性变换关系,就产生了新的问题。首先,一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器,不再保持原有的线性相位了；其次,这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的,即某一频率段的幅频响应近似等于某一常数〔这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性,不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变,如REF_Ref344334166\h图2.4.4所示。图STYLEREF2\s2.4.SEQ图\*ARABIC\s24双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器,双线性变换后,仍得到幅频特性为分段常数的滤波器,但是各个分段边缘的临界频率点产生了畸变,这种频率的畸变,可以通过频率的预畸变来加以校正。也就是将临界模拟频率事先加以畸变,然后经变换后正好映射到所需要的数字频率上。切比雪夫低通数字滤波器的MATLAB设计设计步骤[][]刘舒帆,费诺,陆辉.数字信号处理实验<MATLAB版>[M].XX：XX电子科技大学出版社,2008.176-183.这里采用MTALAB双线性变换法设计IIR切比雪夫低通滤波器,其步骤为：输入给定的数字滤波器设计指标根据公式进行修正,将数字滤波器设计指标转换成模拟滤波器设计指标；确定模拟滤波器的最小结束和截止频率；计算模拟低通切比雪夫原型滤波器的系统传递函数；用双线性变换法将模拟滤波器转换为数字滤波器。MATLAB程序低通数字滤波器的技术指标：通带截止频率,通带衰减,阻带截止频率,阻带衰减。程序如下：wp=0.2*pi;ws=0.6*pi;rp=3;rs=40;%数字滤波器指标Omgp=tan<wp/2>Omgs=tan<ws/2>%转化为模拟滤波器指标〔K=1Fs=0.5;%模拟原型滤波器计算[N,wo]=cheb1ord<Omgp,Omgs,rp,rs,'s'>%计算阶数N和截止频率wo[b,a]=cheby1<N,rp,wo,'s'>%求原型滤波器系数b,a[Ha,Omg]=freqs<b,a>;%求模拟系统频率特性dbHa=20*log<<abs<Ha>+eps>/max<abs<Ha>>>;%求分贝值,加eps以避开0点%计算数字滤波器[bd,ad]=bilinear<b,a,0.5>;%用双线性变换法求数字滤波器系数bd,ad[H,w]=freqz<bd,ad>;%求数字系统频率特性dbH=20*log10<abs<H>/max<abs<H>>>;%化为分贝值%画图subplot<3,2,1>,plot<Omg*Omgp/2/pi,dbHa>,gridtitle<'a模拟幅度响应〔dB>'>;ylabel<'幅度/dB'>;xlabel<'f/Hz'>;subplot<3,2,2>,plot<w/pi,abs<H>>title<'b幅度响应'>;axis<[0,1,0,1]>;subplot<3,2,3>,plot<w/pi,dbH>,gridtitle<'c幅度响应<dB'>;axis<[0,1,-300,100]>;subplot<3,2,4>,plot<w/pi,dbH>,gridtitle<'d幅度响应<dB'>;axis<[0,0.7,-60,5]>;subplot<3,2,5>,plot<w/pi,angle<H>/pi>title<'e相位响应'>;axis<[0,1,-1,1]>;subplot<3,2,6>,zplane<bd,ad>title<'f零极图'>