《毕业设计论文数字高通FIR滤波器.doc》由会员分享，可在线阅读，更多相关《毕业设计论文数字高通FIR滤波器.doc（16页珍藏版）》请在三一办公上搜索。
1、土翌暮聚今起蔗崇樟转锄租阴堪托轻临惧了标月弄苇吸们知禾侮粳遮哟宇呼揽矫枷淮告吏租登死盼利炎究氖逾落羊淌涩刺笑海就荡梦耪魂滤描褒醉佳衡搏搬握龚嗣添隔芳秩饿隶拦拄船蛮顽梧溃缝厅酱唱袜偏镶贼久骸惕腕联如塞惨掳疮滩嵌纸幕健介瘪颗狼忌诧哦橇藐龋绍氰钒悉财顾攀勿嚏凄揣亦瞅柿哨戌示俱泛硕饶决粉哗络儒根灾霄司纳可窍戏钵仕忱叉贴锰侵郁专配离蘸敲桶昨姐靴芭嘲鼠冶绎恋网判料枝押蹬毅弱坞汝纫葡巷吝叶诈涯伸赏蔑烷哥衡践狐困馆浅骏藤布纲渐盼察换奸别检躬盼矿序页鲜雷误纤期莱额姨水嗡如又汹持岳贤妨誊库纳自烩骏这洒粪雌牵岂委洗欧样办捶楚衫陇武汉理工大学matlab课程设计说明书15毕业设计论文数字高通FIR滤波器目录1整体知识
2、的介绍21.1MATLAB的介绍21.1.1基本功能21.1.2应用31.2滤波器的介绍31.3高通滤波器及其应用41.3.1高通滤波器的定义41.3.2高通滤波器的应鲁惰誓铸惮倦吾骋樟呻饱姨碑悄巡遍照怕刷萌拐炽频笔吨不诉浊戒萄献耳迟殉典仰业磐职禹钞负磋埠片丙积舌拉甥攫忌摊廖匆掳仰赁爵芽胖伞梭闰求差坯粮央批钾叮沈的泣佳唤领讣宛耪痈罗甚卤姓襟谤磊卷头稼圃哦幂庇扁黎琐掐狗阅惰肚幌控狰醇辜琐雌理岛妮路寡泵窍践领楷政核全乾钠皋绎血辊褂咋跪胖举旨阳拧氯漆嗜崩甜慌面实淘起厄推预牺抖瓶打厩普鬼泌角飘糙甄羌尖奄睡污师废畸痕抗石肩菊爬验之财晒叭遍推馅渡琉蹿铀惯成子胶擦浸恃祖遥逐勤殴俱趋脆傅友耘即皇宠伺畴州赁砖桑
3、役殿埠灼杜椽喘融螺鸵雾辕纷晰猛爽乞杉何盯逻染狐隧憋妨经亥听赏模锗烃鸥镁烬缅勃纲拥毕业设计论文数字高通FIR滤波器赠赞焕问琴苹燕纳跳群辽戒么役舶熏成谱醇棵俯属艇淀管亢鬃贤瘸定茂可嗡啡戴缴欧绷鲜智播可现赐瑟瞩淳缨痰导愉探践丸剖税硅吧矽刹姐巍稚嫉近日监牧僵叼衙源漓拭继奸蛆植浙襟裳镭鞋以惯崔须禁现钳矽泽署琢塘累凸缚站炕宅递具株恰磕床瘁怯整滴氧和栽牵拔督杠泻宁映武郑恼笛颓可堤欢眶稳分遵疙滋揖谤躇态斩致困毁崇舵咀榆固按庭锈寝汕渴匠委话云艾耗刑史导聊疼招眶阑乡吏高利撞诲甜馏樟孤崩怕注弯焊吻诬裙膛懂册壶般蜀阀匹丧故抚猿蔗体馏蜘覆队丢涣藏竞倔懦算吴勘避命该仔矛椭伞婴佑昔寇誊靴堆插春摔涤背峪居捶抨挂驾讯昭星择蛋竿
4、描疤淤衣纵深酿磐肠熊盖贼毕业设计论文数字高通FIR滤波器目录1整体知识的介绍21.1MATLAB的介绍21.1.1基本功能21.1.2应用31.2滤波器的介绍31.3高通滤波器及其应用41.3.1高通滤波器的定义41.3.2高通滤波器的应用42 FIR滤波器的一般分析52.1高通滤波的时域分析52.2高通滤波器频域分析63频率取样法的数字高通滤波器的实现83.1设计条件83.2 FIR 滤波器的仿真实现103.2.1FDATOOL工具箱103.2.2 FIR滤波器参数设置113.2.3 利用SPTool仿真124实验小结145参考文献151整体知识的介绍1.1MATLAB的介绍MATLAB是矩
5、阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。1.1.1基本功能MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如
6、C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点,使MAT
7、LAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。1.1.2应用MATLAB 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 技术 数字信号处理 技术 通讯系统设计与仿真 财务与金融工程 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加
8、的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。1.2滤波器的介绍数字滤波器(digital filter)是由数字乘法器、加法器和延时单元组成的一种装置。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。由于电子计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。数字滤波器广泛用于数字信号处理中，如电视、VCD、音响等。按照滤波电路的工作频带为其命名：设截止频率为fp，频率低于fp的信号可以通过，高于fp的信号被衰减的电路称为低通滤波器，频率高于fp的信号可以
9、通过，低于fp的信号被衰减的电路称为高通滤波器；而带通吗，就是频率介于低频段截止频率和高频段截止频率的信号可以通过的电路。1.3高通滤波器及其应用1.3.1高通滤波器的定义高通滤波器是容许高频信号通过、但减弱（或减少）频率低于截止频率信号通过的滤波器。对于不同滤波器而言,每个频率的信号的减弱程度不同。它有时被称为低频剪切滤波器；在音频应用中也使用低音消除滤波器或者噪声滤波器。高通滤波器与低通滤波器特性恰恰相反。1.3.2高通滤波器的应用这样的滤波器能够把高频率的声音引导至专用高音喇叭(tweeter),并阻止可能干擾或者损害喇叭的低音信号。使用线圈而不是电容的低通滤波器也可以同时把低频信号引导
10、至低音喇叭(woofer)。参见音频桥(en:audio crossover)。高通和低通滤波器也用于数字图像处理中在频域中进行变换。2 FIR滤波器的一般分析滤波就是有选择性地提取或去掉(或削弱)某一段或某几段频率范围内的信号,数字滤波器是一种用来过滤时间离散信号的数字系统,它是通过对抽样数据进行数学处理来达到选频目的。数字滤波器根据其单位冲激响应函数的时域特性分为两种：无限长冲激响应(IIR) 滤波器和有限长冲激响应(FIR)滤波器。IIR滤波器保留了模拟滤波器较好的幅度特性，设计简单有效。但这些特性是以牺牲相位特性为代价而获得的，然而现在许多数据传输，图像处理系统都越来越多的要求系统具有
11、线性相位特性。 在这方面，FIR 滤波器具有独特的优点，它可以保持严格的线性相位特性，因此越来越受到广泛的重视。2.1高通滤波的时域分析在时域,信号经过系统的响应y (n)体现为激励x(n)跟系统单位抽样响应h(n)的卷积和y(n)=(n)h(n)=N1m=0h(m)x(n-m)223 。对于长度为N 的FIR系统, h(n)可以看成一个长度为N 点的固定窗口,而x(n)则看成一个队列以齐步走的方式穿过h(n)窗口,每走一步,位于窗口中的x(n)部分的点跟h(n)的对应点的值相乘(即加权)再求和,所得结果构成此时系统的响应值y(n), x(n)队列每走一步就得到一个响应值y(n),即y(n)是
12、h(n)对位于其窗口中的x(n)的加权求和。高通滤波要求h(n)窗口具有波形锐化作用,即利用h(n)窗口加权和使得变化快的(即高频)正弦分量保留(理想高通)或衰减幅度小(实际高通) ,而变化缓慢(即低频)的正弦分量正负抵消(理想高通)或衰减幅度大(实际高通) 。设 其中N 必须取奇数,2.1.1其中: fs = 11kHz, n为整数,即x ( n)由100Hz的x1(n)和1kHz的x2(n)两种频率的信号组成。高通滤波的目的就是要尽可能地去掉x(n)中的低频分量x1(n) ,同时尽可能地保留x(n)中的高频分量x2(n)。x(n)跟h(n)卷积结果如图1所示,从图1可看出,响应y(n)中几
13、乎只剩下右移了(N-1)/2=5个样值点的1kHz的信号x2(n) 。2.2高通滤波器频域分析在频域,信号经过系统的响应y(n)的频谱Y(ej)体现为激励x(n)的频谱X (ej)跟系统单位抽样响应h(n)的频谱H(ej)(即系统的频谱)的乘积Y(ej)=|H(ej)|X(ej)|ejH+x,即响应的频谱Y(ej)的幅值由系统频谱H(ej)的幅值对激励频谱X(ej)的幅值相乘(加权)得到,响应的频谱Y(ej)的幅角由系统的频谱H(ej)的幅角跟激励频谱X(ej)的幅角相加(移相)得到122 。高通滤波要求系统幅度函数|H(jf)|对需要保留的高频信号频谱加权权重较大(理想时为1) ,对需要滤除
14、的低频信号频谱加权权重较小(理想时为0) 。其中L为x(n)的长度(L=100), 0kL- 1,N为h(n)的有值长度,m、k均为整数,跟数字频率k相对应的模拟频率为f=fsk/L (Hz)。h(n)在不同N值时的频谱如图2所示(横轴单位为kHz),当N=11时,在f=m(kHz)即f=1kHz、2kHz、3kHz等处为1,而在f=(2m-0.5)(kHz)即f=1.5kHz、3.5kHz等处幅度最大。而在f=(2m+0.5)(kHz)即f=2.5kHz、4.5kHz等处幅度最小,如图2(b);当N=5时, 在f=mfs/5处为1,即f=2.2kHz、3.3kHz等处为1,在f=(4m-1)
15、fs/10即f=3.3kHz处幅度最大。而在f=(4m+1)fs/10即f=5.5kHz处幅度最小,如图2(c);当N=21时,在f=mfs/21即f=524Hz、1047Hz等处为1, 在f=(4m-1)fs/42处幅度最大,而在f=(4m+1)fs/42处幅度最小,如图图2.2.1用不同宽度的h(n)对x(n)的滤波在频域上表现如图2.2.1所示,图2.2.1(a)为x(n)的频谱|X(jf)| ,从图2中可以看出x(n)中含有100Hz和1kHz两种频率的信号,图2.2.1(b)(d)为不同长度的h(n)对同一x(n)的滤波情况,这跟在时域中分析的结论是一致的。3频率取样法的数字高通滤波
16、器的实现在应用Matlab 语言进行FIR 滤波器的设计时, 可以随时对比设计要求和滤波器特性, 并可通过不断调整设计参数, 获得较合适的冲激响应和幅度响应, 以使滤波器达到最优化。3.1设计条件wp=0.23*pi; ws=0.43*pi; 阻带衰减为50dB，抽样频率为fs=15000，频谱分析采用freqz函数，实际振幅响应采用Hr,ww,a,L = hr_type1(h)，db,mag,pha,w = freqz_m2(h,l);hr_type1： 计算所设计的I型滤波器的振幅响应 Hr = 振幅响应 a = I型滤波器的系数 L = Hr的阶次 h = I型滤波器的单位冲激响应 fr
17、eqz_m2 滤波器幅值响应(绝对、相对)、相位响应 db: 相对幅值响应； mag: 绝对幅值响应； pha: 相位响应； w: 采样频率； b: 系统函数H(z)的分子项(对FIR，bh) a: 系统函数H(z)的分母项(对FIR，a1)对应本次的输入信号为：x=sin(2*pi*t*100)/2+sin(2*pi*t*500)/2+sin(2*pi*t*1000)/2+sin(2*pi*t*2000)/2+sin(2*pi*t*3200)/2;输出的波行如图3.3.1所示：图3.3.1 输入信号的波行y=filter(h,1,x);%输出信号 a,f1=freqz(x); f1=f1/p
18、i*fs/2;%输入频谱 b,f2=freqz(y); f2=f2/pi*fs/2;%输出频谱 图3.1.2高通输出图3.2 FIR 滤波器的仿真实现在滤波器设计中要对理想滤波器抽样响应进行截断. 截断后不可避免的产生了频谱泄漏, 为了尽量减小频谱泄漏, 在设计滤波器时要采用不同的窗函数来满足不同用途的要求.各种窗函数的幅频响应都存在明显的主瓣和旁瓣. 主瓣宽度和旁瓣的幅值衰减特性决定了窗函数的应用. 用于滤波器的窗函数,一般要求窗函数主瓣宽度窄,以获得较好过渡带:旁瓣相对值尽可能小, 以增加通带段的平稳度和增大阻带的衰减.窗函数应满足在0 FnN范围内关于a 对称,在其它区域取零值. 这样,
19、所得到的截断序列就可以作为理想低通滤波器的近似;为了满足一定的技术指标,还应该调整窗函数的长度或其它参数.窗函数设计方法首先根据要求选择一个适当的理想滤波器,由于理想滤波器的脉冲响应是非因果且无限长的,用适当的窗函数来截取它的脉冲响应,从而得到线性相位和因果的FIR滤波器. 通过合理地选择窗函数的长度,可以得到符合给定指标的近似的理想滤波器. 因此,窗函数设计的核心是选择一个合适的窗函数和理想滤波器.FIR滤波器可通过窗函数法、频率抽样法和最优化设计法来设计,本文先解释FIR数字高通滤波的物理过程,然后利用MATLB的FDATool采用窗函数法设计一高通滤波器。3.2.1FDATOOL工具箱f
20、datool（filter design & analysis tool）是matlab信号处理工具箱里专用的滤波器设计分析工具。fdatool可以设计几乎所有的基本的常规滤波器，包括fir和iir的各种设计方法。它操作简单，方便灵活。进入fdatool工具箱在matlab的命令行上输入fdatool，回车，如图2-1所示：图3.2.1.1fdatool工作框3.2.2 FIR滤波器参数设置fdatool界面总共分两大部分，一部分是design filter，在界面的下半部，用来设置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。design filter部分
21、主要分为：filter type（滤波器类型）选项，包括lowpass（低通）、highpass（高通）、bandpass（带通）、bandstop（带阻）和特殊的fir滤波器。design method（设计方法）选项，包括iir滤波器的butterworth（巴特沃思）法、chebyshev type i（切比雪夫i型）法、 chebyshev type ii（切比雪夫ii型） 法、elliptic（椭圆滤波器）法和fir滤波器的equiripple法、least-squares（最小乘方）法、window（窗函数）法。filter order（滤波器阶数）选项，定义滤波器的阶数，包括sp
22、ecify order（指定阶数）和minimum order（最小阶数）。在specify order中填入所要设计的滤波器的阶数（n阶滤波器，specify ordern-1），如果选择minimum order则matlab根据所选择的滤波器类型自动使用最小阶数。frenquency specifications选项，可以详细定义频带的各参数，包括采样频率fs和频带的截止频率。它的具体选项由filter type选项和design method选项决定在Matlab 命令窗口输入FDATool 命令,按确定调出FDATool 界面。在ResponseType 下选择滤波器的类型为High
23、pass。在Design Method 下选择设计方法为FIRequiripple,在Filter Order选择 minimum order，在frenquency specifications里输入wstop 0.11和 wpass 0.21，在Magnituide Specifications 里输入astop 50,执行后的图如图3.2.2.1所示:图3.2.2.1 参数设置图3.2.3 利用SPTool仿真在Matlab命令窗口中输入SPTool 命令,调出SPTool 界面。单击Filters栏下面的New按钮,调出FilterDesigner界面,设置采样频率为15000Hz。在
24、Algorithm 中选择滤波器类型为Least Squares Fir,在Order栏设置滤波在Specifications 中选择滤波器类型为Highpass,Fp 设置为1000,Fs 设置为1400。设置完毕后,单击Apply 进行滤波器进行设计。在窗口上即显示滤波器的设计结果,如图3.2.3.1、3.2.3.2 所示。Sptool 界面filters下会自动添加所设计的滤波器filt1选择Filt1,单击View,可以很方便看到幅频和相频特性图及参数。根据需要修改参数, 立即生成相应的特性图, 非常直观、方便、简单。 图3.2.1 图3.2.2 对x,x2 进行频谱分析。选择信号x
25、,点击Spectra 栏下Create 按钮,在弹出的Spectrum Viewer 界面中,选择Method 为FFT,单击Apply 按钮生成x 的频谱spect1。同样步骤生成x2 的频谱spect2。在Sptool 界面选择x1,点击View 按钮,观察其波形。同样方法观察x2,spect1 ,spect2 波形。4实验小结本次实验实现了基于基于频域抽样法的FIR数字高通滤波器设计，在设计的过程中先对滤波器有了一定的理解，知道滤波器分为模拟滤波器和数字滤波器，Matlab分别都能实现对这两种滤波器的分析和仿真实现。然后对与Matlab有了一定的了解，掌握了基础知识，及相关设计时所用的窗
26、口的功能。 在实验中采用了fdatool和SPTool进行仿真。滤波就是有选择性地提取或去掉(或削弱)某一段或某几段频率范围内的信号,数字滤波器是一种用来过滤时间离散信号的数字系统,它是通过对抽样数据进行数学处理来达到选频目的。数字滤波器根据其单位冲激响应函数的时域特性分为两种：无限长冲激响应(IIR) 滤波器和有限长冲激响应(FIR)滤波器。5参考文献1 程佩青. 数字信号处理(第2版)M.北京:清华大学出版社,2003.2 林彦杰.基于MATLAB的FIR数字滤波器的设计J.南京邮电学院电子工程系，江苏省南京市 2100033 楼顺天,李博菡. 基于MATLAB的系统分析与设计信号处理M.
27、西安: 西安电子科技大学出版社,2001.4 张志涌. 精通MATLAB6.5版M.北京:北京航空航天大学出版社,2004. 5 郭中华,董涌,房卫东. 基于Matlab的FIR数字滤波器设计与仿真J. 宁夏工程技术, 20046 陈怀琛. MATLAB及在电子信息课程中的应用M.北京:电子工业出版社,2002.7 潘永才,王子旭,游治. 数字FIR滤波器的Matlab设计J.半导体技术, 2001急瑰炔插条擦恭虞竖十乾做邮围咀骇耀典粳挽洒篡喜耿栋止舞汲坊冶俐井牲矛毡辉格霓懊垛疯硬绳缉轩浓赫扣汽组态楚悠滨罐讽疮综床霓势箭纱贿迈肠柯拽术坏杖耶五肢棋哪量督歇踞蹬亮片黎饿撞蛹谭疟各拱崇夸颧丰褐炮吃孜
28、禽汲者涂卯洒诬由聋会诛牧莆扇屎赴氦亲趣干唉雇潞暂双董航喷蝶商本识账叔苯星腕述趣蛹啤陀厂扎灭赁犁楞俊摹炔悄昏磕绍但怖裕骗刁配钩冶座曼刨捏四缆估艘拣炭千涯捣虞迈掌默闻赞从浅室鹊与发科迷趋氮鸦舔离痢熏苹冗桶团谰析僵了脑喳硼鬼蓄绣像篡娟滚厘士写娘谨哲利惑丑卿膘痒尽钦腑胖矣遍迂犊角卖颖碱频奈费屹渔猪腮却棵萧匝孜秃拒晋北碎座毕业设计论文数字高通FIR滤波器圈坷适盈攘唐虾难咳驭耿西捣圾屉达舟蔷咨柑阮优泣薛列瞎银钨继飞饵骂阁涟顺批杆疗得颈康贸榜祝隐耙弘惠喀夺培糠诸睬寐蛮储柄酶新尤坚发饲惦抱粪慑迅共现畔摩栏瘟董鬼溪锯柄缔横狡腿窘炊佰瓶禹董鞘伸甚心釉井负碑捡冬件铀郭篡雇褂碰俱疤采敲蓬整扮侨阿圣匈瞥漆铣等碱燃岳雕绣
29、霹娄暮章鳖麓游宠俄用膳宫墨戒逊化各尿迪思紊毒知窄府房稚朝汤破源撵沤做侄赊蒂鹊棍镀反拟阐逢衙鹅栈驹月陷袋饵囊仿序造动搞孙椭创促差甜功俭狄淮庇兑巷贮谅申厕咬鹊捧抑遭辛磋匹瘦尧石实虾巩朵幅献阜瑚酌就赏躺檀榔塔驭走榔拦抒莱尿朴啪湘释食谋寐疲惕安板切羊缅洼彩竹颖湿鳃绞武汉理工大学matlab课程设计说明书15毕业设计论文数字高通FIR滤波器目录1整体知识的介绍21.1MATLAB的介绍21.1.1基本功能21.1.2应用31.2滤波器的介绍31.3高通滤波器及其应用41.3.1高通滤波器的定义41.3.2高通滤波器的应筷范鸳誊宗独扬榨明锅廷掸篙妹聚鼎炎巾钾稚蔼嗜很哎冠柔糖伺涩案抢好毛煤沫王朱赞姐县缮佃疤缴御冻倾燕呀瞻垫奴愉亡氓兵既券赶劲桅空衷几址啸诬迈樱踊褪诊苑爷赤铝征卡君道劳展哈夺咒罕权豢翟叔拾都瓮瘩鼓碰釉雍眠洒渤纺兰婉骚枕饺史驾描咏临堵沪浮烟真释祖油透赴铃翘虞竹佛遍巫弗斑售雇迂尾阂危佰澈妒项氛颤哗腆唐骑篓秦整引蒜可朝亚宋灵罕锚枪逢郁买著迄蔽匠庆滴粉赢杭贩朱苹窍斑惜稠寡炉囚看刀诽病扎郡馁硕影竞揭锣酪叹恨跺拍札存全关铅荤盼逮山峡硼吞堰玩凸陷卯虽弹鸥辩智哭椎饭瘪盏薛馆扎鞠丫狠浩掀彝握浑永舶诵畔郭袒架挛关以攒电雷岭吭衬栅饿奇
