欢迎来到沃文网! | 帮助中心 分享知识,传播智慧!
沃文网
换一换
首页 沃文网 > 资源分类 > DOC文档下载
 

利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带通IIR滤波器。

  • 资源ID:20015       资源大小:234.36KB        全文页数:25页
  • 资源格式: DOC        下载权限:游客/注册会员/VIP会员    下载费用:10积分 【人民币10元】
快捷注册下载 游客一键下载
会员登录下载
三方登录下载: QQ登录   微博登录  
下载资源需要10积分 【人民币10元】
邮箱/手机:
温馨提示:
支付成功后,系统会自动生成账号(用户名和密码都是您填写的邮箱或者手机号),方便下次登录下载和查询订单;
支付方式: 微信支付    支付宝   
验证码:   换一换

加入VIP,免费下载资源
 
友情提示
2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,既可以正常下载了。
3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带通IIR滤波器。

课程设计任务书学生姓名专业班级指导教师工作单位信息工程学院题目利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带通IIR滤波器。初始条件①MATLAB软件②数字信号处理与图像处理基础知识要求完成的主要任务利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带通IIR滤波器。参考书目时间安排指导教师签名2011年月日系主任(或责任教师)签名2011年月日武汉理工大学MATLAB应用课程设计说明书2目录目录................................................................2Abstract............................................................41绪论..............................................................52IIR数字滤波器设计的原理与方法....................................62.1IIR数字滤波器设计的原理.....................................62.2IIR数字滤波器设计的基本方法................................83IIR带通滤波器的MATLAB设计......................................103.1IIR带通滤波器的设计流程....................................113.2IIR带通滤波器的设计步骤....................................12心得与体会.........................................................23参考文献...........................................................24武汉理工大学MATLAB应用课程设计说明书3摘要数字滤波是数字信号分析中最重要的组成部分之一,数字滤波与模拟滤波相比,具有精度和稳定性高、系统函数容易改变、灵活性高、不存在阻抗匹配问题、便于大规模集成、可实现多维滤波等优点。数字滤波器的作用是利用离散时间系统的特性对输入信号波形(或频谱)进行加工处理,或者说利用数字方法按预定的要求对信号进行变换。从广义讲,数字滤波是由计算机程序来实现的,是具有某种算法的数字处理过程。MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。关键词IIR数字滤波器,MATLAB,仿真武汉理工大学MATLAB应用课程设计说明书4AbstractDigitalfilterdigitalsignalanalysisisthemostimportantpartofdigitalfilterandsimulationfilter,itischaracterizedbyhighprecisionandstability,systemfunctiontochange,flexibility,highimpedancematchingproblemdoesnotexist,easyforlargescaleintegrated,mayrealizethemultidimensionalfiltering,etc.Theroleofdigitalfilterisusingdiscretetimethecharacteristicsofthesystemofthesignalwaveorspectrumprocess,orusingadigitalaccordingtotherequirementsofthescheduledtosignaltrans.Broadlyspeaking,digitalfilterisbythecomputerprogramtorealize,ishassomekindofalgorithmdigitalprocessingprocess.MATLABisreleasedbytheAmericanmathworkscompanymainlyfaceofscientificcalculation,visualizationandinteractiveprogramdesignofthehigh-techcomputingenvironment.Itwillnumericalanalysis,calculation,scientificdatavisualizationandnonlineardynamicsystemandsimulation,andmanyotherstrongfunctionintegrationinaneasytouseWindowsenvironment,forscientificresearch,engineeringdesign,andtoeffectivelythenumericalcalculationmanyfieldsofscienceprovidesacomprehensivesolution,andtoalargedegreefromthetraditionaltheinteractiveprogramminglanguagesuchasC,Fortraneditmode,whichrepresentsthecurrentinternationalscientificcomputingsoftwareadvancedlevel.KeywordsIIRdigitalfilters,MATLAB,thesimulation武汉理工大学MATLAB应用课程设计说明书51绪论在现代通信系统中,由于信号中经常混有各种复杂成分,所以很多信号的处理和分析都是基于滤波器而进行的。数字滤波器可以理解为是一个计算程序或算法,将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中,使信号按预定的形式变化。数字滤波器有多种分类,根据数字滤波器冲激响应的时域特征,可将数字滤波器分为两种,即无限长冲激响应IIR滤波器和有限长冲激响应FIR滤波器。切比雪夫滤波器,又名“车比雪夫滤波器”,是在通带或阻带上频率响应幅度等波纹波动的滤波器。切比雪夫滤波器来自切比雪夫分布,以“切比雪夫”命名,是用以纪念俄罗斯数学家巴夫尼提列波维其切比雪夫。根据频率响应曲线波动位置的不同,切比雪夫滤波器可以分为以下两种○1I型切比雪夫滤波器在通带(或称“通频带”)上频率响应幅度等波纹波动的滤波器○2II型切比雪夫滤波器在阻带(或称“阻频带”)上频率响应幅度等波纹波动的滤波器称为“II型切比雪夫滤波器”。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。本文设计一个IIR数字带通滤波器。IIR数字滤波器具有无限宽的冲激响应,与模拟滤波器相匹配,所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤武汉理工大学MATLAB应用课程设计说明书6波器设计法。在对滤波器实际设计时,整个过程的运算量是很大的。设计完成后对已设计的滤波器的频率响应要进行校核,要得到幅频、相频响应特性,运算量也是很大的。在这种情况下,滤波器设计就要进行大量复杂的运算,单纯的靠公式计算和编制简单的程序很难在短时间内完成。利用MATLAB强大的计算功能进行计算机辅助设计,可以快速有效地设计数字滤波器,大大地简化了计算量。2IIR数字滤波器设计的原理与方法2.1IIR数字滤波器设计的原理由信号与系统的知识可知滤波器的实现实际上就是频域上的信号与系统冲激响应函数的卷积运算。一个理想的低通或高通滤波器其冲激响应函数是时域无限长的,显然其卷积运算是不可能直接实现的。相比于FIR数字滤波器利用窗函数将冲击函数截短,使用有限长度的冲激响应函数与信号样值序列进行卷积运算产生吉不斯现象。IIR滤波器没有对冲激响应函数进行截短来实现滤波器。由于冲激响应函数无限长,所以IIR滤波器不可能由频域卷积方法来设计,通常由频域直接设计。设IIR数字滤波器的输入序列为nx,则IIR滤波器的输入序列nx和输出序列ny之间的关系可用下列方程式表示NjjMiijnyainxbny10(公式2-1)其中ib和ja滤波器的系数,且中至少有一个非零。具有无限持续时间冲激响应,需要用递归模型来实现,其系统函数为011MkkkNkkkbzYzHzXzaz(公式2-2)由传递函数可以发现无限长单位冲激响应滤波器应有如下特点(1)单位冲激响应是无限长。武汉理工大学MATLAB应用课程设计说明书7(2)系统传递函数在有限平面上有极点存在(3)结构上存在着输出到输入的反馈,也就是结构上是递归型。设计IIR滤波器的任务就是寻求一个因果、物理上可实现的系统函Hz,使其频率响应jwHe满足所希望得到的频域指标,即符合给定的通带截止频率、阻带截止、通带衰减和阻带衰减。利用冲激响应不变法设计数字滤波器时可能会导致频域混叠现象,为了克服这一问题,需要找到由s平面到z平面的另外的映射关系,这种关系应保证1s平面的整个jΩ轴仅映射为z平面单位圆上的一周;2若Gs是稳定的,由Gs映射得到的Hz也应该是稳定的;3这种映射是可逆的,既能由Gs得到Hz,也能由Hz得到Gs;4如果Gj01,那么0jHe=1。双线性Z变换满足以上4个条件的映射关系,其变换公式为211szsTz(公式2-3)双线性变换法的特点(1)由于映射是单值变换关系,避免了冲激响应不变法的频率响应混迭现象。(2)当频率增加时,变换关系存在严重的非线性。因此,若采用双线性变换法进行设计时,需要对模拟频率进行预畸变之后再设计模拟滤波器。可将模拟系统函数分解成并联或级联的低阶子系统后再对每个子系统函数分别进行双线性变换。即分解为低阶的方法可在模拟系统函数上进行,而模拟系统函数的分解有大量的图表可以利用,分解起来较方便。双线性Z变换的基本思路首先将整个S平面压缩到一条从-π/Ts变换到2π/Ts的横带里,然后通过标准的变换关系sTze将横带变换到整个Z平面上去,这样就得到了S平面与Z平面间的一一对应的单值关系。武汉理工大学MATLAB应用课程设计说明书8图2-1双线性变换法S平面到z平面的映射关系2.2IIR数字滤波器设计的基本方法IIR数字滤波器的设计方法有多种方法,如频率变换法、数字域直接设计以及计算机辅助设计等。可以借助于模拟滤波器的设计方法设计出模拟滤波器,利用冲激响应不变法或双线性变换法转换成数字滤波器,再用硬件或软件实现;另可以直接在频域或时域中进行设计,设计时需要计算机作辅助工具。频率变换法也分为三种(1)先把一个归一化原型模拟低通滤波器经模拟频带变换成所需要的模拟滤波器,然后再通过冲激响应不变法或双线性变换法数字化为所需类型的数字滤波器。(2)把两步合成一步来实现,即把归一化原型模拟低通滤波器经模拟频带变换成所需要类型的模拟滤波器的公式与用双线性变换法得到响应数字滤波器的公式合并,就可直接从模拟低通滤波器归一化模型通过一定的频率变换关系,一步完成各类型数字滤波器的设计。(3由模拟低通原型先变换成数字低通滤波器,然后利用数字频带变换法变换为各型数字滤波器。基于冲激响应有频率混迭效应,故采用双线性变换法,则IIR数字滤波器设计的基本步骤如下(1)根据任务,确定性能指标。在设计一个滤波器之前,首先根据工程实际Re[z]σσjΩs平面jΩ1s平面jImtz平面武汉理工大学MATLAB应用课程设计说明书9的需要确定滤波器的技术指标如边界频率wp通带截止频率,ws阻带截止频率,wc;阻带最小衰减As和通带最大衰减RP;(2)对截止频率进行预畸变,以减轻双线性变换的压缩效应。3将数字滤波器的技术指标转换成模拟滤波器指标。利用冲激响应不变法与双线性变换法进行频率间的转换,主要是边界频率Wp与Ws的转换。4用模拟滤波器设计方法得到模拟滤波器的传输函数Has;可借助巴特沃斯Butterworth滤波器、切比雪夫Chebyshev滤波器、椭圆Cauer滤波器、贝塞尔Bessel滤波器等,这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。5进行适当的模拟频率变换,将得到的模拟低通滤波器变换为满足第一步要求的模拟低通,高通,带通或带阻滤波器。(6)映射实现。利用双线性变换法将模拟滤波器Has转换成数字滤波Hz。7用有限精度算法实现这个系统函数Hz包括选择运算结构、选择合适的字长、有效数字处理方法。8用适当的软、硬件技术实现。包括采用通用计算机软件、数字滤波器硬件或者软硬件结合,确定DF采用的结构将会影响其精度、稳定性、经济性及运算速度等很多重要性质。前面介绍第i级二阶基本节的差分方程为2121nwanwianxzwiiiii(公式2-4)2121nwbnwibnwzyiiiii(公式2-5)结合前面讲过的第i级二级基本节的递归函数,只需存储中间变量及其上两个时刻值,则由当前时刻的输入及本级二阶节的滤波器系数即可求出本级二级节的输出。整个滤波器程序可以二阶节乘虚的多次循环来实现。武汉理工大学MATLAB应用课程设计说明书103IIR带通滤波器的MATLAB设计IIR带通滤波器的设计框图如下图3-1IIR带通滤波器的设计框图模拟低通滤波器的设计模拟频带变换低通变高通、带通等滤波器变换模拟变换为数字IIR滤波器武汉理工大学MATLAB应用课程设计说明书113.1IIR带通滤波器的设计流程图3-2IIR带通滤波器的设计流程本文设计的IIR带通滤波器是从低通变换过来的,利用的是双线性变换以及切比雪夫II滤波器的原型,其具体的设计流程为上图所示。首先根据题目要求确定带通滤波器的技术指标,先要进行频率的预畸变,并且归一化频率,再设计出切比雪夫II模拟低通滤波器,并求出其阶数等相关参数。其次利用双线性变换法设计数字带通滤波器,再调用函数进行双线性变换,并求出分子、分母的系数向量。然后通过画图求出其幅频响应、相频响应、幅度特性曲线与零极点,并画出波形图。最后进行验证,看所设计的滤波器能否达到要求的指标,若能达到,则说明该滤波器设计符合要求。开始确定带通滤波器的技术指标设计切比雪夫低通滤波器并求出相关参数求相应的幅频响应与相频响应双线性变换模拟低通滤波器转化为带通滤波器结束武汉理工大学MATLAB应用课程设计说明书123.2IIR带通滤波器的设计步骤(1)根据设计流程,首先确定所要设计的数字带通滤波器的相关指标①通带截止频率wp10.4π,wp20.6π,通带最大衰减Rp2dB;②阻带截止频率ws10.2π和ws20.8π,阻带最小衰减Rs30dB;③取样间隔T0.1s。其实现程序如下(程序中pi代表π)Ts0.1;Fs1/Ts;%取样周期或频率Rp2;%通带最大衰减Rs30;%阻带最小衰减wp10.4*pi;%通带、阻带上、下限截止频率wp20.6*pi;ws10.2*pi;ws20.8*pi;(2)频率的预畸变。双线性变换中无法避免的一个问题即是频率的非线性偏移,因为数字频率的最大值为π,而模拟频率可以向无穷延伸,两者之间又要保持一一对应的映射关系。双线性变换中的模拟角频率与数字角频率之间的关系为2tan2T公式3-1表明S平面与Z平面是单值的一一对应关系,即频率轴是单值变换关系。虽然避免了脉冲响应不变法的频率响应的混叠现象,但是经过变换后,得到的幅频响应特性各分段边缘频率不能保持原来的比例关系,必须通过预修正加以校正。做法是将数字频率按=2/T*tanw/2的关系,变成模拟频率,利用这组做过修正的模拟频率来设计模拟带通滤波器作为模拟原型。Wp12/T*tanwp1/2;Wp22/T*tanwp2/2;Wp[Wp1,Wp2];%模拟滤波器的通带截止频率Ws12/T*tanws1/2;武汉理工大学MATLAB应用课程设计说明书13Ws22/T*tanws2/2;Ws[Ws1,Ws2];%模拟滤波器的阻带截止频率BWWs2-Ws1;%模拟滤波器的带宽Omegaw0sqrtWs1*Ws2;%模拟滤波器的中心频率(2)设计切比雪夫模拟低通滤波器。%求模拟低通滤波器的阶数与边缘频率[N,OmegaC]cheb2ordWp,Ws,Rp,Rs,s;%求切比雪夫II型模拟低通滤波器的零、极点与增益[z0,p0,k0]cheb2apN,Rs;利用函数[N,OmegaC]cheb2ordWp,Ws,Rp,Rs,’s’,通过给定滤波器的技术指标Wp、Ws、Rp、Rs,求得滤波器的阶数N与边缘频率OmegaC。Wp、Ws、与OmegaC均在[0,1]区间归一化,以π弧度为单位。利用函数[z,p,k]cheb2apN,Rs,来设计一个阶数为N,阻带波动为Rs的归一化切比雪夫II型原型滤波器,得到左半平面零极点。数组Z中返回零点,数组P中返回极点,并且返回增益K。利用上面知识可设计一个ChebyshevI型带通模拟滤波器,带通滤波器的阶数为十阶,其通带波值为3dB,低带边缘为/5,高带边缘为[z,p,k]cheb1ap5,3;[A,B,C,D]zp2ssz,p,k;将模拟滤波器转换为状态空间形式u10.1*2*pi;u20.5*2*pi;单位是rad/sBwu2-u1;Wosqrtu1*u2;[At,Bt,Ct,Dt]lp2bpA,B,C,D,Wo,Bw;[b,a]ss2tfAt,Bt,Ct,Dt;从状态空间形式转换为转换函数所需的参数形式wlinspace0.01,1,500*2*pi;计算频率向量hfreqsb,a,w;计算频率响应semilogyw/2/pi,absh,grid;绘制图像武汉理工大学MATLAB应用课程设计说明书14xlabelFreqyencyHz;[Ad,Bd,Cd,Dd]bilinearAt,Bt,Ct,Dt,2,0.1;利用双线性变换法设计一个采样频率为2Hz,低带边缘频率为0.1Hz的切比雪夫Ⅰ模拟滤波器[bz,az]ss2tfAd,Bd,Cd,Dd;将其参数转换为状态空间形式fvtoolbz,az图3-3ChebyshevI型带通模拟滤波器图3-4Matlab实现滤波器的离散化00.10.20.30.40.50.60.70.80.9-150-100-500NormalizedFrequencyrad/sampleMagnitudedBMagnitudeResponsedB00.10.20.30.40.50.60.70.80.9110-710-610-510-410-310-210-1100FreqyencyHz武汉理工大学MATLAB应用课程设计说明书15图3-4Matlab实现滤波器的离散化图3-5双线性变换后的滤波器幅度与相位响应(3)设计归一化的模拟原型带通滤波器求原型滤波器的分子系数AnalogBk0*realpolyz0;求原型滤波器的分母系数AnalogArealpolyp0;模拟低通到模拟带通的分子、分母系数的变换[BandB,BandA]lp2bpAnalogB,AnalogA,Omegaw0,BW;双线性变换模拟带通与数字带通的分子分母系数的变换[DigitalB,DigitalA]bilinearBandB,BandA,Fs;变为二阶节级联结构[sos,G]tf2sosDigitalB,DigitalA;利用函数ppolyA来计算模拟滤波器的分子、分母系数向量,因其为实数,因此用real函数取其实部,即可得模拟滤波器的分子、分母系数向量。这两个函数实现的功能可以用函数[Bs,As]zp2tfz,p,k直接求得传递函数的分子、分母系数向量。利用函数[BandB,BandA]lp2bpAnalogB,AnalogA,Omegaw0,BW,将模拟域的低通变为带通,并且得到模拟带通滤波器的分子、分母系数向量,Omegaw0取为中心频率,BW为带宽。利用函数[DigitalB,DigitalA]00.10.20.30.40.50.60.70.80.9-140-120-100-80-60-40-200NormalizedFrequencyrad/sampleMagnitudedBMagnitudeResponsedB武汉理工大学MATLAB应用课程设计说明书16bilinearBandB,BandA,Fs,双线性变换为数字带通滤波器的指标,如分子、分母的系数向量。函数[sos,G]tf2sosDigitalB,DigitalA,即把z变换传递函数的直接形式转换成级联形式。需要注意的是,这个函数是针对以z的负幂排列的多项式开发的,虽然可以推广到s域,但连续系统传递函数是按s的正幂排列的,要使两者一致,关键是使分子、分母系数向量同长,两序列中各元素的幂次排列一致。(4)求数字带通滤波器的幅频、相频特性、及其群延迟%求数字带通滤波器的幅频特性[Hz,Wz]freqzDigitalB,DigitalA,1024,whole;%将数字带通滤波器的幅频特性转化为分贝表示dbHz20*log10absHzeps/maxabsHz;%求数字带通滤波器的相频特性φ=angleHz%求数字带通滤波器的群延迟特性grdgrpdelayDigitalB,DigitalA,Wz;函数[Hz,Wz]freqzDigitalB,DigitalA,1024,whole可以求数字带通滤波器的幅频特性,而其幅度(即模值)的最大值可以归一化为1,则其模值(单位为dB)即可以用公式表示为dbHz20*log10absHzeps/maxabsHz。函数φ=angleHz可求得其相频特性,而对于一个滤波器来说,要满足其线性相位,则其群延迟要为一常数,则其相位特性必须为一直线,即满足t。因此,利用群延迟函数grdgrpdelayDigitalB,DigitalA,Wz可以判断所设计的滤波器是否是线性相位,如果不符合,可以更改参数加以较正或者用其他方法重新设计,从而方便了设计。(5)作图武汉理工大学MATLAB应用课程设计说明书174IIR带通滤波器程序及波形4.1参考程序Ts0.1;Fs1/Ts;Rp2;Rs30;Wp10.4*pi;Wp20.6*pi;Ws10.2*pi;Ws20.8*pi;武汉理工大学MATLAB应用课程设计说明书18频率的预畸变Wp12/Ts*tanwp1/2;Wp22/Ts*tanwp2/2;Wp[Wp1,Wp2];模拟滤波器的通带截止频率Ws12/Ts*tanws1/2;Ws22/Ts*tanws2/2;Ws[Ws1,Ws2];模拟滤波器的带宽BWWp2-Wp1;Omegaw0sqrtWs1*Ws2;模拟滤波器的中心频率Omegaw0sqrtWp1*Wp2;求模拟低通滤波器的阶数与边缘频率[N,OmegaC]cheb2ordWp,Ws,Rp,Rs,s;[N,OmegaC]cheb1ordWp,Ws,Rp,Rs,s求切比雪夫模拟低通滤波器的零、极点与增益[z0,p0,k0]cheb2apN,Rs;[z0,p0,k0]cheb1apN,Rp;设计归一化的模拟原型带通滤波器求原型滤波器的分子系数AnalogBk0*realpolyz0;求原型滤波器的分母系数AnalogArealpolyp0;模拟低通到模拟带通的分子、分母系数的变换[BandB,BandA]lp2bpAnalogB,AnalogA,Omegaw0,BW;双线性变换模拟带通与数字带通的分子分母系数的变换[DigitalB,DigitalA]bilinearBandB,BandA,Fs;变为二阶节级联结构[sos,G]tf2sosDigitalB,DigitalA;求数字带通滤波器的幅频、相频特性、及其群延迟求数字带通滤波器的幅频特性武汉理工大学MATLAB应用课程设计说明书19[Hz,Wz]freqzDigitalB,DigitalA,1024,whole;将数字带通滤波器的幅频特性转化为分贝表示dbHz20*log10absHzeps/maxabsHz;求数字带通滤波器的相频特性求数字带通滤波器的群延迟特性grdgrpdelayDigitalB,DigitalA,Wz;作图subplot2,3,1;plotWz/pi,absHz;title幅频响应;xlabel,ylabel幅度|Hz|;axis[0,1,0,1.1];setgca,XTickMode,manual,XTick,[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1];grid;subplot2,3,4;plotWz/pi,dbHz;title模值dB;xlabel频率(单位\pi);ylabel分贝(dB);axis[0,1,-40,5];setgca,XTickMode,manual,XTick,[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1];setgca,YTickmode,manual,YTick,[-50,-30,-2,0];setgca,YTickLabelMode,manual,YTickLabels,[50;30;2;0];grid;subplot2,3,2;plotWz/pi,angleHz/pi;title相频响应;xlabel;ylabel单位\pi;axis[0,1,-1,1];setgca,XTickMode,manual,XTick,[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1];grid;subplot2,3,5;ylabel单位\dB;xlabel单位\pi;zplaneDigitalB,DigitalA;axis[-1.1,1.1,-1.1,1.1];subplot2,3,3;plotWz/pi,grd;title群延迟xlabel频率(单位\pi);ylabel样本;axis[0,1,0,8]setgca,XTickMode,manual,XTick,[0,ws1/pi,wp1/pi,wp2/pi,ws2/pi,1];画高刻度线武汉理工大学MATLAB应用课程设计说明书20setgca,YTickmode,manual,YTick,[00.510];grid;setgcf,color,w00.20.40.60.8100.51幅频响应幅度|Hz|00.20.40.60.813020模值dB频率(单位)分贝(dB)00.20.40.60.81-1-0.500.51相频响应单位-101-101RealPartImaginaryPart00.20.40.60.8100.511.522.533.544.555.566.577.58群延迟频率(单位)样本图4-1以上是根据切比雪夫II型滤波器设计的带通滤波器,如果以切比雪夫I型滤波器设计,则只需改变几个参数则可,其程序如下所示(只附改变的参数与程序)。%模拟滤波器的带宽BWWp2-Wp1;%模拟滤波器的中心频率Omegaw0sqrtWp1*Wp2;%求切比雪夫I型模拟低通滤波器的阶数与边缘频率[N,OmegaC]cheb1ordWp,Ws,Rp,Rs,s;%求切比雪夫I型模拟低通滤波器的零、极点与增益武汉理工大学MATLAB应用课程设计说明书21[z0,p0,k0]cheb1apN,Rp;00.20.40.60.8100.51幅频响应幅度|Hz|00.20.40.60.813020模值dB频率(单位)分贝(dB)00.20.40.60.81-1-0.500.51相频响应单位-101-101RealPartImaginaryPart00.20.40.60.8100.511.522.533.544.555.566.577.58群延迟频率(单位)样本图4-2根据上述的仿真波形,可以看出(1)由ChebyshevI型、ChebyshevII型低通原型变换成带通模型的幅频特性可以用分贝形式表示。对于ChebyshevI型的幅频响应在通带内为等纹波衰减、阻带内为单调减小的,且通带(归一化)截止频率在波形上显示为0.4与0.6,阻带(归一化)截止频率在波形上显示为0.2与0.8,与设计要求基本一致。对于ChebyshevII型的幅频响应来说,在通带内单调减小、阻带内为等纹波衰减的,而通带、阻带(归一化)截止频率为0.4与0.6,0.2与0.8,均满足设计要求。(2)此次设计的IIR带通滤波器的阶数可由函数[N,OmegaC]cheb1ordWp,Ws,Rp,Rs,s或[N,OmegaC]cheb2ordWp,Ws,Rp,Rs,s求得,N3。则从ChebyshevI型带通武汉理工大学MATLAB应用课程设计说明书22滤波器的幅频响应曲线上看,当N=3时,在0(归一化后为0--1)范围内通带波动只有两个波谷,阻带则单调减小,这与理论内容一致。对于ChebyshevII型带通滤波器来说,其幅频响应曲线的通带波动在0(归一化后为0--1)范围内阻带波动只有两个波峰,通带则单调减小,这与理论内容一致。因此,对于滤波器阶数这一指标来说,满足要求。(3)对于分贝化后的幅度特性,ChebyshevI型带通滤波器的通带最大衰减大约为2dB,阻带衰减最小为30dB,而ChebyshevII型带通滤波器的阻带最小衰减大约为30dB,通带最小衰减大约为2dB,这与带通滤波器的技术指标相同。因此,在衰减性能上,所设计的滤波器达到了要求。(4)对于相频响应来说,ChebyshevI型带通滤波器在通带内的相频响应曲线接近为一条平滑曲线,在阻带内则存在衰减、畸变。ChebyshevII型带通滤波器在通带截止频率以内的相频响应曲线几乎为一条直线,但在通带截止频率处有很大角度的转折,使相频特性在通带内产生了的畸变,而在阻带内存在很大的畸变。因此,对于双线性变换法,其相位特性得不到满足,必须用其它方法加以校正或采用其它方法重新设计滤波器,使之满足线性相位。武汉理工大学MATLAB应用课程设计说明书23心得与体会此次Matlab课程设计是在一定的理论基础之上进行的,在先修课程信号与系统中,有Matlab设计方面的知识积累,但是到了自己真正设计的时候,遇到了很多问题。而且经过了自己的亲身实践,学到了许多实践方面的知识。首先,因为并没有学过与IIR滤波器设计相关的知识,设计程序时没有切入点。后来,通过上网查询相关文档和去图书馆查阅有关Matlab设计的书籍,仔细思考,逐步摸索,利用课余时间,读了大量关于滤波器设计的知识,慢慢理解,让我豁然开朗柳暗花明。通过自学,我了解到IIR数字滤波器的设计和模拟滤波器的设计有着紧密的关系。通常要先设计出适当的模拟滤波器,再通过一定的频带变换把它转换成为所需的数字IIR滤波器。在比较设计滤波器的方法上应该明确其技术指标以及某些参数的实际意义。比如本文用双线性变换法设计数字带通滤波器时,必须先将频率归一化,并且进行频率预畸变,然后设计模拟滤波器,再利用频率变换法将模拟低通变为模拟带通,最后经过双线性变换法将模拟带通变换为数字带通。如果不进行频率预畸变,那么设计出来的带通滤波器的幅频特性与相频特性将会产生很严重的畸变,使设计的结果不满足给定的要求,在实际中会造成很严重的危害。在做本次课程设计的过程中,我深深地感受到了自己所学到知识的有限,明白了只学好课本上的知识是不够的,要通过图书馆和互联网等各种渠道来扩充自己的知识。在实验过程中我们曾经遇到过问题。但是从中我们学习到了如何对待遇到的困难,进一步培养了我们一丝不苟的科学态度和不厌其烦的耐心,所有的这些心得会对我以后的学习和工作有帮助作用。武汉理工大学MATLAB应用课程设计说明书24参考文献[1]董长虹,Matlab信号处理与应用.国防工业出版社,2005.1[2]陈怀琛编著,MATLAB及在理工课程中的应用指南.西安电子科技大学出版社,2000[3]刘泉,阙大顺编著,数字信号处理原理与实现.北京电子工业出版社,2005.8[4]胡广书,数字信号处理理论、算法与实现.北京清华大学出版社,1997[5]VinayK.Ingle编著,数字信号处理及其MATLAB实现.北京电子工业出版社,1998[6]陈怀琛编著,数字信号处理教程MATLAB释义与实现(第二版).北京电子工业出版社,2008.10武汉理工大学MATLAB应用课程设计说明书25本科生课程设计成绩评定表姓名尹嘉笛性别女专业、班级通信0906题目MATLAB课程设计利用matlab软件结合双线性变换法设计切比雪夫IIR带通滤波器答辩或质疑记录成绩评定依据最终评定成绩(以优、良、中、及格、不及格评定)指导教师签字年月日........忽略此处.......

注意事项

本文(利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带通IIR滤波器。)为本站会员(星星008)主动上传,沃文网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知沃文网(发送邮件至2622162128@qq.com或直接QQ联系客服),我们立即给予删除!

温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服点击这里,给沃文网发消息,QQ:2622162128 - 联系我们

版权声明:以上文章中所选用的图片及文字来源于网络以及用户投稿,由于未联系到知识产权人或未发现有关知识产权的登记,如有知识产权人并不愿意我们使用,如有侵权请立即联系:2622162128@qq.com ,我们立即下架或删除。

Copyright© 2017-2019 www.wodocx.com ,All Rights Reserved |陕ICP备19002583号  

陕公网安备 61072602000132号     违法和不良信息举报:0916-4228922