利用MATLAB仿真软件系统结合双线性变换法设计一个数字巴特沃斯高通IIR滤波器

Matlab课程设计任务书学生姓名：管行专业班级：通信0606指导教师：王晟工作单位：信息工程学院题目:利用MATLAB仿真软件系统结合双线性变换法设计一个数字巴特沃斯高通IIR滤波器基础强化训练目的1．理论目的《数字信号处理》课程设计的目的之一是为了巩固课堂理论学习，并能用所学理论知识正确分析数字信号处理的基本问题和解释数字信号处理的基本现象。2．实践目的《数字信号处理》课程设计的目的之二是通过设计具体的各种滤波器掌握滤波器设计方法和步骤。训练内容和要求利用MATLAB仿真软件系统结合双线性变换法设计一个数字巴特沃斯高通IIR滤波器在数字信号处理平台上（PC机﹑MATLAB仿真软件系统和TC++编程环境）进行软件仿真设计，并进行调试和数据分析。课程设计题目由指导教师提供，每人隶属一组完成任务，每组完成的内容不能雷同（按学号分组）初始条件①MATLAB软件②数字信号处理与图像处理基础知识时间安排：第21周，安排任务（鉴3-204，7月14日）第21周，仿真设计（鉴主13楼计算机实验室）第21周，完成（答辩，提交报告，演示）指导教师签名：2008年7月14日系主任（或责任教师）签名：年月日目录摘要........................................................................1Abstract....................................................................21设计项目要求与说明........................................................32系统设计..................................................................32.1设计思路............................................................32.2设计方法对比........................................................42.3典型模拟滤波器比较......................................................52.4设计步骤............................................................53仿真程序的设计与调试......................................................63.1数字域指标变换成模拟域指标..........................................63.2数字域频率进行预畸变................................................63.3模拟滤波器的设计....................................................73.4模拟滤波器变成数字滤波器............................................93.5理论计算数字滤波器的仿真...........................................114.程序调试中出现的问题.....................................................125.总结与体会..............................................................13参考文献...................................................................14附录一总程序如下.........................................................15附录二设计数字滤波器函数总结.............................................17摘要此报告重点介绍了用双线性不变法设计IIR数字滤波器的基本流程，比较了各种设计方法的优缺点，总结了模拟滤波器的性能特征。最后以双线性不变法设计了一个高通巴特沃斯IIR数字滤波器，介绍了设计步骤，然后在Matlab环境下进行了仿真与调试，实现了设计目标。关键词：Matlab双线性不变法IIR数字滤波器巴特沃斯高通AbstractThisreportintroducedwithemphasisofthebasicflowofdesigningtheIIRdigitfilterbythebilinearpoliticalreform,comparedwitheachkindofdesignmethod'sgoodandbadpoints,summarizedanalogfilter'sperformancecharacteristic.FinallydesignonebythebilinearpoliticalreformtopassButterworthhighIIRdigitfilter,introducedthedesignprocedure,thenhascarriedonthesimulationandthedebuggingundertheMatlabenvironment,hasachievedtheprojectobjective.Keyword:MatlabbilinearitypoliticalreformtheIIRdigitalfilterpasshighbutterworthMatlab课程设计——设计一个数字巴特沃斯高通IIR滤波器1设计项目要求与说明课题要求设计一个IIR数字滤波器，高通，采用双线性变换法，用巴特沃斯实现，用matlab软件对其进行仿真与调试。本设计将先说明用双线性法设计IIR数字滤波器的原理，然后写出基于matlab的软件设计流程。在对设计进行调试，分析实验数据。2系统设计此部分将详细介绍IIR数字滤波器的设计流程，比较各种设计方案的优劣。2.1设计思路IIR滤波器设计的主要方法是先设计低通模拟滤波器，然后转换为高通、带通或带阻数字滤波器。对于其他如高通，带通，则通过频率变换转换为设计相应的高通，带通等。在设计的全过程的各个步骤，matlab都提供相应的工具箱函数，使得IIR数字滤波器设计变得非常简单。总的来说，我的设计思路主要有以下两种：思路一：从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。模拟域冲激响应不变法频率变换双线性变换法图2-1先频率变换再离散思路二：先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在Z域内经数字频率变换为所需类型的数字滤波器。数字域双线性变换法频率变换图2-2先离散再频率变换归一化模拟低通原型模拟高,带通或带阻数字高,带通或带阻归一化模拟低通原型数字原型低通数字高,带通或带阻以上两种思路都可以，我最后选择了第一种思路进行设计，即先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。2.2设计方法对比方案一：冲激响应不变法冲激相应不变法是从时域出发，要求数字滤波器的激响应h(n)对应于模拟滤波器ha(t)的等间隔抽样，h(n)=ha(nT)，其中T是抽样周期，因此时域逼近良好。优点：○1h(n)完全模仿模拟滤波器的单位抽样响应时域逼近良好○2线性相位模拟滤波器转变为线性相位数字滤波器缺点：○1对时域的采样会造成频域的“混叠效应”，故有可能使所设计数字滤波器的频率响应与原来模拟滤波器的频率响应相差很大，○2不能用来设计高通和带阻滤波器。只适用于限带的低通、带通滤波器方案二：双线性变换法双线性变换法是从频域出发，使DF的频率响应与AF的频率响应相似的一种变换法。直接使数字滤波器的频率响应，逼近模拟滤波器的频率响应,进而求得H(z)。优点：○1避免了频率响应的混迭现象○2在特定AF和特定DF处，频率响应是严格相等的，它可以较准确地控制截止频率的位置。○3它是一种简单的代数关系，设计十分方便。缺点：○1除了零频率附近，与之间严重非线性，即线性相位模拟滤波器变为非线性相位数字滤波器○2要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变○3对于分段常数型AF滤波器，经双线性变换后，仍得到幅频特性为分段常数的DF.但在各个分段边缘的临界频率点产生畸变，这种频率的畸变，可通过频率预畸变加以校正。方案三：频率变换法设计思想：○1从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。○2先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在Z域内经数字频率变换为所需类型的数字滤波器。综上所述，频率变换法是基于双线性变换法或冲激响应不变法（用于频带变换），此处选择双线性变换法主要是基于要获得严格的频率响应，以及较准确地控制截止频率的位置。而由于此种方法是一种简单的代数关系，设计也十分方便。2.3典型模拟滤波器比较1、Butterworth巴特沃斯滤波器：它具有单调下降的幅频特性；即最平幅度。2、Chebyshev切比雪夫滤波器：在通带或阻带等波纹，可提高选择性。3.Bessel贝塞尔滤波器：在通带内有较好的线性相位特性。4.Ellipse椭圆滤波器：其选择性相对前三种是最好的。此处选择巴特沃斯主要是想获得最平稳的幅频响应。而不计较相位特性，而用双线性变换法也会将线性相位变为非线性相位。2.4设计步骤如设计一个数字低通滤波器，其技术指标为：通带临界频率fp，通带内衰减小于rp；阻带临界频率fs，阻带内衰减大于αs；采样频率为FS○1将指标变为角频率wp=fp*2*pi;ws=fs*2*pi;（式2-1）○2将数字滤波器的频率指标{Wk}由wk=(2/T)tan(Wk/2)转换为模拟滤波器的频率指标{wk},由于是用双线性不变法设计，故先采取预畸变。22wpptgT;22wsstgT;（式2-2）○3将高通指标转换为低通指标，进而设计高通的s域模型○4归一化处理1ppp；11pp；（式2-3）ssp；1ss；（式2-4）/10/101011lg/lg()2101SpSN；（式2-5）由式2-3，2-4，2-5计算出N，查表可得模拟低通滤波器的阶数，从而由下式确定模拟高通滤波器的参数。11111111H(z)=H()|H()|caLzzspzzsp（式2-6）3仿真程序的设计与调试3.1数字域指标变换成模拟域指标其程序为：fp=400;fs=300;Rp=1;Rs=20;wp=fp*2*pi;ws=fs*2*pi;FS=1000;T=1/FS；程序执行结果为：wp=2.5133e+003ws=1.8850e+003与实际计算结果相符。3.2数字域频率进行预畸变其程序为：wp2=2*tan(Wp/2)/T;ws2=2*tan(Ws/2)/T;经过预畸变，可以发现频率变为：wp2=6.1554e+003ws2=2.7528e+0033.3模拟滤波器的设计其程序为%设计模拟滤波器[N,Wn]=buttord(wp2,ws2,Rp,Rs,'s')[z,p,k]=buttap(N);%创建Buttord低通滤波器原型[Bap,Aap]=zp2tf(z,p,k);%由零极点转换为传递函数的形式figure(1)freqs(Bap,Aap);%模拟低通滤波器的频率响应title('模拟滤波器（低通原型）的频率响应')[Bbs,Abs]=lp2hp(Bap,Aa