武汉理工大学通信工程基于MATLAB和双线性变换法的数字巴特沃斯高通IIR滤波器 19页

武汉理工大学《MATLAB原理与应用》课程设计说明书课程设计任务书学生姓名：专业班级：指导教师：工作单位：信息工程学院题目：基于MATLAB和双线性变换法的数字巴特沃斯高通IIR滤波器初始条件：1）MATLAB软件2）数字信号处理与图像处理基础知识3）双线性变换法的原理和算法4）巴特沃斯高通滤波器的性能指标要求完成的主要任务：（包括课程设计工作量及技术要求，以及说明书撰写等具体要求）1、课程设计工作量：1周。2、技术要求：1）利用MATLAB仿真软件系统结合双线性变换法设计一个数字巴特沃斯高通IIR滤波器2）在数字信号处理平台上（PC机、MATLAB仿真软件）进行软件仿真设计并进行调试和数据分析。3、查阅至少5篇参考文献。按《武汉理工大学课程设计工作规范》要求撰写设计报告书。全文用A4纸打印，图纸应符合绘图规范。时间安排:1、_一年—_月_—日，布置课设具体实施计划与课程设计报告格式的要求说明。2、_一年—_月_—日至年—M—日，方案选择和电路设计。3、——年_—月—一日至年一M—日，电路调试和设计说明书撰写。4、——年_—月—一日，上交课程设计成果及报告，同时进行答辩。指导教师签名：年月日系主任（或责任教师）签名: 武汉理工大学《MATLAB原理与应用》课程设计说明书目录摘要IAbstractII1数字滤波器11.1数字滤波器的基本概念11.2数字滤波器的分类11.3数字滤波器的MATLAB现12双线性变换法32.1双线性变换法知识简介32.2双线性变换法设计数字滤波器原理33设计任务及方案选择53.1设计任务53.2方案设计论证与选择54数字滤波器的设计74.1设计方法74.2设计步骤84.3设计实例及仿真85总结与体会146参考文献15 武汉理工大学《MATLAB原理与应用》课程设计说明书摘要数字滤波器是数字信号处理的基础，用来对信号进行过滤、检测与参数估计等处理，在通信、图像、语音、雷达等许多领域都有着十分广泛的应用。尤其在图像处理、数据压缩等方面取得了令人瞩目的进展和成就。鉴于此，数字滤波器的设计就显得尤为重要。为了得到较纯净的真实信号，对巴特沃斯数字滤波器的幅频特性，设计方法以及设计步骤进行了研究，利用MATLA的真软件结合双线性变换法设计出了数字巴特沃斯高通IIR滤波器，并且分析了巴特沃斯数字滤波器的相关特性。关键词：数字滤波器；巴特沃斯；MATLAB双线性变换法I 武汉理工大学《MATLAB原理与应用》课程设计说明书AbstractDigitalfilteristhebasisofdigitalsignalprocessing,usedforsignalfiltering,detectionandparameterestimation,incommunications,image,voice,radar,andinmanyotherfieldshaveverywiderangeofapplications.Especiallyinimageprocessing,datacompressionhasmaderemarkableprogressandachievements.Inviewofthis,digitalfilterdesignisparticularlyimportant.Inordertoobtainmorepuretruesignal,theamplitudefrequencycharacteristicofButterworthdigitalfilter,designmethodanddesignstepswerestudiedbyusingMATLABsimulationsoftware,combinedwiththebilineartransformationmethodforthedesignofButterworthhighpassdigitalIIRfilter,andananalysisoftherelatedcharacteristicsofthedigitalfilterwithBartVos.Keywords:digitalfilter;Butterworth;MATLAB;bilineartransformationmethodii 武汉理工大学《MATLAB原理与应用》课程设计说明书1数字滤波器1.1数字滤波器的基本概念数字滤波器是数字信号处理的重要基础，是对信号都是过滤检测与参数估计等处理过程中，它是使用最为广泛的一种线性系统。数字滤波器处理的对象是经由采样期间将模拟信号转换而得到的数字信号。数字滤波器是指完成信号滤波处理功能的，用有限精度算法实现的离散时间线性非时变系统。数字滤波器的输入是一组数字量。它本身既可以是用数字硬件装配而成的一台用于完成给定运算的专用数字计算机，也可以是将所需的运算编写的程序通过计算机来执行。数字滤波器具有稳定性高、精度高、灵活性大等突出优点。随着数字技术的发展，用数字技术实现滤波器的功能愈来愈受到人们的重视，并得到了广泛的应用。数字滤波器的数学运算通常有两种实现方式。一种是频域法，即利用FFT快速运算方法对输入信号进行离散傅里叶变换，分析其频谱，然后再根据所希望的频率特性进行滤波，再利用傅里叶反变换得到时域信号。这种方法具有较好的频域选择特性和灵活性，并且由于信号频率与所希望的频率特性是简单的相乘关系，所以它比计算等价的时域卷积要快得多。另一种方法是时域法，这种方法是通过离散的抽样数据做差分数学运算来达到滤波目的的。1.2数字滤波器的分类数字滤波器按照不同的分类方法，可分为许多种，但总体来讲可以分成两大类。一类称为经典滤波器，即一般滤波器，特点是输入信号中有用的频率成分和希望滤除的频率成分各占有不同的频带，通过一个合适的选频滤波器达到滤波的目的。但对于一般滤波器，如果信号和干扰的频带互不重叠，则不能完成对干扰的有效滤除，这是需要采用另一类所谓的现代滤波器，例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器，这些滤波器可按照随机信号内部的一些统计分布规律，从干扰中最佳地提取信号。从功能上,一般数字滤波器可以分为低通、高通、带通、带阻和全通等，此种分类方法是和模拟滤波器一样的。1.3数字滤波器的MATLA取现MATLAB软件为滤波器的设计应用提供了丰富而简便的方法，如函数方法和图形工具方法等，使原来非常繁琐复杂的程序设计变成了简单的函数调用，为滤波器的设计和实现开辟了广阔的天地。数字滤波器的一般设计过程为：按照实际需要，确定滤波器的性能要15 武汉理工大学《MATLAB原理与应用》课程设计说明书求；用一个因果的、稳定的离散线性时不变系统，去逼近这一性能指标；用有限精度的运算实现所设计的系统；通过模拟，验证所设计的系统是否符合给定性能要求。15 武汉理工大学《MATLAB原理与应用》课程设计说明书2双线性变换法2.1双线性变换法知识简介脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Z平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法。双线性变换法是一种从模拟滤波器变换数字滤波器的方法，其主要思想是利用非线性压扩的方法将整个s平面压缩成宽为2n/T的带状区，即将整个s平面频率轴上的频率范围压缩到-九/T〜sT兀/T之间，再用z=e转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到sTS1平面的-冗/T〜兀/T一条横带里；第二步再通过标准变换关系Z=e将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象。由图2-1看出，在零频率附近，模拟角频率。与数字频率⑴之间的变换关系接近于线性关系；但当。进一步增加时，⑴增长得越来越慢，最后当。一8时，⑴终止在折叠频率⑴=冗处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，从而消除了频率混叠现象。图2.1双线性变换频率对应关系2.2双线性变换法设计数字滤波器原理将s平面j。轴压缩变换到s1平面j。轴上的-n/T到兀/T一段，可以采用以下变换关系:二ctan15 武汉理工大学《MATLAB原理与应用》课程设计说明书其中C为常数，一般来说，为了使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系，就引入该待定常数Co这样Q=±g变换到^Y/，—。变换到仪肛可将上式写成j11e=c-M」皿2-e2☆s=jC,s=jCi,解析延拓到整个s平面和si平面，可得S1TS1T~^2SiTe2-e21-es=c-S1TS1T=c.tST22~21eee将si平面为2兀/T的带状区域与整个z平面之间的标准映射关系z=es1T带入上式得s平面和z平面的单值映射关系为1-z」二cJ1z在MATLABK双线性z变换可以通过bilinear函数实现，其调用格式为:lBZ1AZ1=bilinearB,A,FS;其中B,A为模拟滤波器传递函数G(s)的分子分母多项式的系数向量，而BZ,AZ为数字滤波器的传递函数H(z)的分子分母多项式的系数向量3设计任务及方案选择3.1设计任务课题要求设计一个IIR数字滤波器，高通，采用双线性变换法，用巴特沃斯实现，用MATLA歆件对其进行仿真与调试。本设计将先说明用双线性法设计IIR数字滤波器的原理，然后写出基于MATLAB勺软件设计流程。在对设计进行调试，分析实验数据。3.2方案设计论证与选择IIR滤波器设计的主要方法是先设计低通模拟滤波器，然后转换为高通、带通或带阻数字滤波器。对于其他如高通，带通，则通过频率变换转换为设计相应的高通，带通等。在设计的全过程的各个步骤，MATLABfB15 武汉理工大学《MATLAB原理与应用》课程设计说明书提供相应的工具箱函数，使得IIR数字滤波器设计变得非常简单。总的来说，我的设计方案主要有以下两种：方案一：从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。图3.1先频率变换再离散变换方案二：先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在Z域内经数字频率变换为所需类型的数字滤波器。换法图3.2先离散再频率变换以上两种方案都是可行的，我选择了第一种方案进行设计，即先在模拟域内经频率变换成为所需类型的模拟滤波器，然后进行双线性变换，由S域变到Z域，从而的到所需类型的数字滤波器。因为高通滤波器的设计方法是先将要设计的滤波器的技术特性指标通过频率转换成模拟低通滤波器的技术指标，再根据这些性能指标设计出低通滤波器的传递函数，为了设计程序的简洁，故选择了方案一。15 武汉理工大学《MATLAB原理与应用》课程设计说明书4数字滤波器的设计4.1设计方法方案一：冲激响应不变法冲激相应不变法是从时域出发，要求数字滤波器的激响应h(n)对应于模才以滤波器ha(t)的等间隔抽样，h(n)=ha(nT),其中T是抽样周期，因此时域逼近良好。优点：1)h(n)完全模仿模拟滤波器的单位抽样响应时域逼近良好2)线性相位模拟滤波器转变为线性相位数字滤波器缺点：1)对时域的采样会造成频域的“混叠效应”，故有可能使所设计数字滤波器的频率响应与原来模拟滤波器的频率响应相差很大。2)不能用来设计高通和带阻滤波器。只适用于限带的低通、带通滤波器方案二：双线性变换法双线性变换法是从频域出发，使DF的频率响应与AF的频率响应相似的一种变换法。直接使数字滤波器的频率响应，逼近模拟滤波器的频率响应，进而求得H(z)。优点：1)避免了频率响应的混迭现象。3)在特定AF和特定DF处，频率响应是严格相等的，它可以较准确地控制截止频率的位置。4)它是一种简单的代数关系，设计十分方便。缺点：1)除了零频率附近，。与C之间严重非线性，即线性相位模拟滤波器变为非线性相位数字滤波器。5)要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变6)对于分段常数型AF滤波器，经双线性变换后，仍得到幅频特性为分段常数的DF.但在各个分段边缘的临界频率点产生畸变，这种频率的畸变，可通过频率预防变加以校正。方案三：频率变换法设计思想：1)从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。2)先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在Z域内经数字频率变换为所需类型的数字滤波器。综上所述，频率变换法是基于双线性变换法或冲激响应不变法(用于频带变换)，此处选择双线性变换法主要是基于要获得严格的频率响应，以及较准确地控制截止频率的位置。而由于此种方法是一种简单的代数关系，设计也十分方便。15 武汉理工大学《MATLAB原理与应用》课程设计说明书4.2设计步骤⑴用公式Q=-tanQ；对数字滤波器H(z)边界频率与进行预畸变，得到同类型的模拟滤T2波器Hd(S附边界频率C(即使给定的边界频率以模拟频率的形式给出，也要先对这些频率进行数字化0=2炳/Fs,然后再进行预防变。)(2)用相应的频率变换公式将Hd(S)的频率指标变换成原型低通滤波器Hlp(S)的频率指标。⑶用相应的方法设计巴特沃斯模拟低通滤波器HLP(S)。(4)用相应的频率逆变换将模拟低通滤波器系统函数Hlp(S)变换为Hd(S)。21-Z」(5)用s=y的双线性变换式将HD(S)变换为H(z)oT1Z4.3设计实例及仿真设计一个数字巴特沃斯高通IIR滤波器，其技术指标为：通带临界频率fp=400HZ,通带p内衰减小于Rp=1dB,阻带临界频率fs=300Hz,阻带内衰减大于Rs=20dB,采样频率为Fs=1000HZ。(1)将以模拟频率的形式给出的边界频率数字化fp=400;fs=300;Rp=1;%通带最大衰减Rp=1dBRs=20;%阻带最小衰减Rs=20dBwp=fp*2*pi;ws=fs*2*pi;Fs=1000;T=1/Fs;%采样频率与采样周期Wp=wp/(Fs);%＞模拟频率数字化Ws=ws/(Fs);(2)频率预畸变，得到同类型模拟滤波器的边界频率wp2,ws2wp2=2*tan(Wp/2)/T;%预畸变求滤波器通带临界频率(3)模拟巴特沃斯低通滤波器的设计[N,Wn]=buttord(wp2,ws2,Rp,Rs,"s");[z,p,k]=buttap(N);%^"J建Buttord低通滤波器原型[Bap,Aap]=zp2tf(z,p,k);%由零极点转换为传递函数的形式15 武汉理工大学《MATLAB原理与应用》课程设计说明书figure⑴freqs(Bap,Aap);%莫拟低通滤波器的频率响应title("模拟滤波器(低通原型)的频率响应‘)下图是模拟巴特沃斯低通原型滤波器的仿真结果模拟滤波器(低通原型)的频率响应20O001210，：10Frequencytra±s)1010Frequency：rads)001图4.1模拟巴特沃斯低通原型滤波器分析：1.要设计的是巴特沃斯滤波器，通过巴特沃斯最小阶数选择函数[N,Wn]=buttord(wp,ws,Rp,Rs,"s");其中wp为通带边界频率，ws为阻带边界频率，单位是rad/s。Rp,Rs分别是通带波纹和阻带衰减，单位为dB.函数返回值N为模拟滤波器的最小阶数，Wn为模拟滤波器的截止频率，单位是rad/s。MATLAB空制台返回N=4,说明该滤波器是4阶。2.buttap是模拟低通滤波器原型设计函数。用[z,p,k]=buttap(N)；得到[z,p,k]后，可求出滤波器系数B,Ao其中z表示零点，p表示极点，k表示增益。3.freqs求模拟滤波器HLP(S)的频率响应函数。(4)模拟巴特沃斯高通滤波器的设计[Bbs,Abs]=lp2hp(Bap,Aap,Wn);%莫拟低通变高通figure(2)freqs(Bbs,Abs);15 武汉理工大学《MATLAB原理与应用》课程设计说明书title("模拟滤波器的频率响应‘)通过转换后将模拟低通原型滤波器转换为了模拟高通滤波器，其波形如下图所示6W1-O模拟滤•波器的频率响应100行Frequency(rad/s)1口”(sio>sl』5wip)