关闭

关闭

关闭

封号提示

内容

首页 基于MATLAB的滤波器的设计

基于MATLAB的滤波器的设计.doc

基于MATLAB的滤波器的设计

蘑菇头仿不来
2019-01-26 0人阅读 举报 0 0 暂无简介

简介:本文档为《基于MATLAB的滤波器的设计doc》,可适用于IT/计算机领域

关于基于MATLAB的滤波器的设计.doc文档,大发快三大小单双口诀拥有内容丰富的相关文档,站内每天千位行业名人共享最新资料。

基于MATLAB高通滤波器的设计摘要:滤波器可以分为模拟滤波器和数字滤波器两大类其中模拟滤波器的设计已经很成熟不仅有完整的设计公式还有完善的图表和曲线供查阅另外还有一些典型的优良滤波器类型共我们使用。数字滤波器可以分为IIR滤波器和FIR滤波器其中IIR滤波器的设计方法有间接法和直接法两种FIR滤波器不能采用间接法设计其常用的方法有窗函数法、频率采样法和切比雪夫等纹逼近法。关键词:模拟滤波器、数字滤波器、FIR滤波器、IIR滤波器。一.IIR模拟高通滤波器的设计设计方法选频型模拟滤波器按幅频特性可分成低通、高通、带通和带阻滤波器但设计滤波器时总是先设计低通滤波器再通过频率变换将低通滤波器转换成希望类型的滤波器。低通滤波器的设计方法有技术指标和逼近法、巴特沃斯滤波器和切比雪夫滤波器设计法等。设计步骤)通过高通与低通滤波器的映射关系式将希望设计的高通滤波器的指标转换成相应的低通滤波器的指标。)设计相应的低通系统函数)再用低通到高通滤波器的映射关系将设计好的模拟低通滤波器的函数式转换成希望设计的高通滤波器的系统函数。实例设计本例将设计一个巴特沃斯模拟高通滤波器要求通带边界频率为kHz,阻带边界频率为kHz通带最大衰减为dB阻带最小衰减为dB。 )低通到高通滤波器的映射关系式公式公式将希望设计的高通滤波器的指标转换成相应低通滤波器的指标。为了计算简单令=。则由公式可求的归一化阻带边界频率为转换得到低通滤波器的指标为:通带边界频率=阻带边界频率=通带最大衰减阻带最小衰减。)为了设计方便该处将调用MATLAB函数buttord和butter来设计。)用公式将转换成希望设计的高通滤波器的系统函数。该处调用MATLAB函数lphp实现低通到高通的变换。基于MATLAB的模拟高通滤波器设计wp=ws=Rp=As=N,wc=buttord(wp,ws,Rp,As,'s')B,A=butter(N,wc,'s')wph=*pi*BH,AH=lphp(B,A,wph)fk=::wk=*pi*fkQ=freqs(B,A,wk)subplot(,,)plot(wk,*log(abs(Q)))gridonxlabel('归一化频率')ylabel('幅度(dB)')axis(,,,)title('归一化低通滤波器Q(p)')fkh=:wkh=*pi*fkhH=freqs(BH,AH,wkh)subplot(,,)plot(fkh,*log(abs(H)))gridonxlabel('频率(Hz)')ylabel('幅度(dB)')axis(,,,)title('高通滤波器H(s)')仿真如下:图二.IIR数字高通滤波器的设计设计方法IIR数字滤波器的设计方法有两种一种是脉冲响应不变法另一种是双线性变换法。本例设计的是高通滤波器如果要用脉冲响应不变法设计为避免频谱混叠则要在高通滤波器之前加保护滤波器滤除高于折叠频率以上的频带颇为不便所以此处将用双线性变换法进行设计。在利用MATLAB设计的时候可以更简便可以直接调用数字高通滤波器函数进行设计设计步骤)确定所需类型数字滤波器的技术指标)将所需类型数字滤波器的边界频率转换成相应类型模拟滤波器的边界频率转换公式为公式。)将相应类型模拟滤波器计数指标转换成模拟低通滤波器技术指标。)设计模拟低通滤波器。)通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。)采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。数字高通滤波器的传统方法设计设计一个数字高通滤波器要求通带介质频率通带衰减不大于dB阻带截止频率阻带衰减不小于dB。希望采用巴特沃斯型滤波器。下边是传统设计方法)确定数字高通技术指标:)将高通数字滤波器的及数字表转换成高通模拟了滤波器的及数字表:令T=s预畸变校正得到模拟边界频率:)模拟低通滤波器的技术指标计算如下:对带通边界频率归一化即将和带入公式求出归一化低通滤波器的阻带截止频率=)设计归一化模拟滤波器取查表可得归一化模拟低通原型系统函数为)利用频率变换公式即公式将转换成模拟高通:)利用双线性变换法将模拟高通数字高通:基于MATLAB的数字高通滤波器的设计传统的设计步骤用MATLAB实现起来比较麻烦对于新手而言程序过长容易出错所以该处用比较简单的另一种方法设计直接调用MATLAB中的数字高通滤波器函数进行设计该方法程序短容易理解和应用。设计程序如下:wpz=wps=rp=rs=N,wc=buttord(wps,wpz,rp,rs)Bz,Az=butter(N,wc,'high')f=::w=f*piH=freqz(Bz,Az,w)plot(f,*log(abs(H)))gridonxlabel('wpi')ylabel('幅度dB')axis(,,,)title('高通数字滤波器耗损函数')仿真图如下图三.FIR数字滤波器的设计设计方法FIR滤波器的设计和IIR滤波器的设计方法有很大的不同FIR滤波器不能用间接法进行设计其设计方法主要有以下几种:窗函数法、频率采样法和切比雪夫波纹逼近法等在本设计中将采用较为简单的窗函数法设计。其中窗函数设计法里面又分为矩形窗、三角窗、汉宁窗、哈明窗、布莱克曼窗以及凯瑟窗等。设计步骤)根据对阻带衰减及过渡带的指标要求选择窗函数的类型并估计窗口长度N。选择原则为:在保证阻带衰减满足要求的情况下尽量选择主瓣窄的窗函数。)构造希望逼近的频率响应函数。)计算。)加窗得到设计结果:。实例设计用窗函数法设计线性相位高通FIR滤波器要求通带介质频率阻带截止频率通带最大衰减阻带最小衰减)选择窗函数计算窗函数长度。已知阻带最小衰减查表可知汉宁窗和哈明窗均满足要求我们选择汉宁窗。由于其中解的,对于高通滤波器必须取奇数取由汉宁窗公式有。)构造:式中,)求出)加窗基于MATLAB的窗函数设计程序:wp=piws=piBt=wpwsN=ceil(*piBt)N=Nmod(N,)wc=(wpws)pihn=fir(N,wc,'high',hanning(N))n=:N  plot(n,hn)stem(n,hn,'')xlabel('n')ylabel('h(n)')axis(,,,)title('h(n)波形')真图如下仿:

类似资料

该用户的其他资料

学前教育系顶岗实习手册.doc

工业设计项目申请与评审制度.doc

矿山区队长的剖析材料.doc

幼儿园礼仪.doc

山东省农业旅游示范点检查标准.doc

职业精品

精彩专题

建筑施工合同范本大全

施工合同亦称“工程合同”或“包工合同”。指发包方 (建设单位) 和承包方 (施工单位) 为完成商定的建筑安装工程施工任务,明确相互之间权利、义务关系的书面协议。对于我们日常接触比较多,且需求量比较大的施工合同,这些模板也许能够帮到你。

用户评论

0/200
上传我的资料

热门资料排行换一换

  • 林中乐队3

  • 北京体育大学本科生毕业论文开题报告

  • 攀成钢1、2 转炉检修投标

  • 3电力工程施工安全操作规程

  • 张艺谋:电影导演工作与电影制作

  • Haworthia瓦苇属学名—中文…

  • 药品包装材料与药物相容性试验指导原…

  • 财务会计报告练习题

  • 四年级上册语文复习计划

  • 资料评价:

    / 12
    所需积分:0 立即下载

    VIP

    在线
    客服

    免费
    邮箱

    大发快三大小单双口诀服务号

    扫描关注领取更多福利