利用MATLAB结合频率取样法设计数字高通FIR滤波器.docx
- 文档编号:1156835
- 上传时间:2022-10-18
- 格式:DOCX
- 页数:14
- 大小:167.78KB
利用MATLAB结合频率取样法设计数字高通FIR滤波器.docx
《利用MATLAB结合频率取样法设计数字高通FIR滤波器.docx》由会员分享,可在线阅读,更多相关《利用MATLAB结合频率取样法设计数字高通FIR滤波器.docx(14页珍藏版)》请在冰豆网上搜索。
利用MATLAB结合频率取样法设计数字高通FIR滤波器
课程设计任务书
学生:
专业班级:
指导教师:
工作单位:
题目:
利用MATLAB结合频率取样法设计数字高通FIR滤波器
要求完成的主要任务:
1.利用频率取样法设计一个数字高通FIR滤波器
2.画出高通滤波器的幅频响应
课程设计进度安排
序
号
阶段容
所需时间
1
方案设计
1天
2
软件设计
2天
3
系统调试
1天
4
撰写报告
1天
合计
5天
指导教师签名:
年月日
系主任(或责任教师)签名:
年月日
摘要
MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
本文介绍了如何利用MATLAB仿真软件系统及数字信号处理所学知识利用频率采样法设计一个数字高通滤波器。
以此来巩固课堂理论学习,并能用所学理论知识正确分析信号处理的基本问题和解释信号处理的基本现象。
关键字:
MATLAB;数字信号处理;数字滤波器;频率采样法
Abstract
MATLABisreleasedbytheUnitedStatesmathworksmainlyforscientificcomputing,visualizationandinteractiveprogramdesignedhigh-techcomputingenvironment.Itnumericalanalysis,matrixcomputation,scientificdatavisualizationaswellasnon-lineardynamicsystemsmodelingandsimulation,andmanyotherpowerfulintegratedinaneasy-to-useWindowsenvironment,scientificresearch,engineeringdesignandtheneedforeffectivenumericaltheeditmodemanyscientificfieldsprovidesacomprehensivesolution,andinlargeparttogetridofthetraditionalnon-interactiveprogramminglanguage(suchasC,Fortran),onbehalfoftheadvancedleveloftoday'sinternationalscientificcomputingsoftware.
ThisarticledescribeshowtouseMATLABsimulationsoftwaresystemsanddigitalsignalprocessinglearnedknowledgeusingfrequencysamplingmethodtodesignadigitalhigh-passfilter.Inordertoconsolidatethetheoreticalclassroomlearning,andbasicquestionsandexplainbasicsignalprocessingphenomenoncanbelearnedtheoreticalknowledgetocorrectlyanalyzethesignalprocessing.
Keyword:
MATLAB;digitalsignalprocessing;digitalfilter;frequencysamplingmethod
1FIR数字滤波器
1.1FIR滤波器的特点
FIR滤波器的脉冲响应h(n)是有限长的(0≤n≤N-1),其z变换为:
(式1)
它是z-1的(N-1)阶多项式,在有限z平面(0 FIR滤波器的基本结构可以理解为一个分节的延时线,把每一节的输出加权累加,可得到滤波器的输出,FIR滤波器的冲激响应h(n)是有限长的,数学上M阶FIR滤波器可以表示为: y(n)=(式2) 其系统函数为: H(z)=(式3)普通的直接型FIR滤波器结构如图1所示。 图1FIR滤波器的直接型结构 FIR滤波器最突出的优点有2个: 一是只要对h(n)附加一定的条件,很容易获得严格的线性相位特性; 二是由于H(z)的极点位于原点z=0处,始终满足稳定条件,所以FIR滤波器永远稳定。 三是FIR滤波器由于单位脉冲响应是有限长的,因而可以用快速傅里叶变换(FFT)算法来实现过滤信号,从而可大大提高运算效率。 但是,要取得很好的衰减特性,FIR滤波器H(z)的阶次比IIR滤波的要高。 1.2FIR数字滤波器设计方法 IIR滤波器设计中的各种变换法对FIR滤波器设计是不适用的,这是因为那里是利用有理分式的系统函数,而FIR滤波器的系统函数只是z-1的多项式。 FIR的设计任务是选择有限长度的脉冲响应h(n),得到系统函数H(z),使幅频特性满足技术指标要求,同时使相频特性达到线性相位。 常用设计方法: (1)窗函数法 (2)频率采样法 (3)切比雪夫等波纹逼近法。 人们最感兴趣的是FIR滤波器具有线性相位的相频特性。 对非线性相位的FIR滤波器,一般可以用IIR滤波器来代替,因为同样幅度特性,IIR滤波器所需阶数比FIR滤波器的阶数要少得多。 1.3线性相位FIR数字滤波器的条件和特点 1.3.1线性相位条件 对于长度为N的h(n),传输函数为 (式4) H(ejω)=Hg(ω)ejθ(ω)(式5) 式中,Hg(ω)称为幅度特性,θ(ω)称为相位特性。 注意,这里Hg(ω)不同于|H(ejω)|,Hg(ω)为ω的实函数,可能取负值,而|H(ejω)|总是正值。 H(ejω)线性相位是指θ(ω)是ω的线性函数,即 θ(ω)=-τω,τ为常数(式6) 如果θ(ω)满足θ(ω)=θ0-τω,θ0是起始相位 严格地说,此时θ(ω)不具有线性相位,但以上两种情况都满足群时延是一个常数,即 (式7) 也称这种情况为线性相位。 1.3.2线性相位FIR滤波器的幅度特性与相位特性 线性相位FIR滤波器的幅度特性与相位特性如下图: 图2线性相位FIR滤波器的幅度特性与相位特性一览表 在设计时,要注意选择合适的h(n)对称形式(奇或偶)和h(n)长度N(奇数或偶数)。 如要设计高通滤波器,只能选情况1和情况4;要设计低通滤波器,只能选情况1和情况2。 2利用频率采样法设计FIR滤波器 2.1用频率采样法设计滤波器的基本原理 待设计的滤波器的传输函数用Hd(ejω)表示,可按下列思路进行设计: 1它在ω=0到2π之间等间隔采样N点,得到Hd(k) (式8) 2N点Hd(k)进行IDFT,得到h(n) (式9) 式中,h(n)作为所设计的滤波器的单位取样响应。 3h(n)求系统函数H(z) (式10) 将插值公式重写如下 (式11) 此式就是直接利用频率采样值Hd(k)形成滤波器的系统函数。 用频率采样法设计线性相位滤波器的条件: FIR滤波器具有线性相位的条件是h(n)是实序列,且满足h(n)=h(N–1–n),其传输函数应满足的条件是 (式12) (式13) (式14) (式15) 且Hg(π)=0。 在ω=0~2π之间等间隔采样N点,将ω=ωk代入式(4~7)中,并写成k的函数: (式16) (式17) ,N为奇数(式18) ,N为偶数且(式19) (式20) 说明N等于奇数时Hg(k)对(N–1)/2偶对称,N等于偶数时,Hg(k)对N/2奇对称,且Hg(N/2)=0。 对于高通滤波器,这里N只能取奇数。 截止频率为ωc,采样点数N,Hg(k)和θ(k)用下面公式计算 (式21) 以上是用频率采样法设计滤波器的基本原理。 2.2线性相位的约束条件 以h(n)为偶对称,N为奇数的情况进行分析。 1)FIR的频响具有线性相位的一般表达式 当h(n)为偶对称,N为奇数时,则 (式22) 而且幅度函数H(w)应为偶对称,即 (式23) 2)采样值H(k)具有线性相位的约束 (式24) 其中,表示采样值的模(纯标量),表示其相角。 因此,在采样点上具有线性相位的条件应为: (式25) 而且,必须满足偶对称,即: (式26) 实际滤波器的传输函数,与理想的传输函数Hd(ejω)间存在误差,如图2 图3频率采样的响应 需要讨论逼近误差问题及其改进措施。 2.3逼近误差及其改进措施 2.3.1产生误差的原因 从图3可看出,实际的H(ejω)与理想的Hd(ejω)相比,误差主要体现在一是通带和阻带出现波动,二是过渡带加宽,与窗函数设计法情况类似,产生误差的原因可从时域和频域两方面进行分析。 从时域分析: 如果Hd(ejω)有间断点,那么相应单位取样响应hd(n)应是无限长的。 这样,由于时域混叠,引起所设计的h(n)和hd(n)有偏差。 为此,希望在频域的采样点数N加大。 N愈大,设计出的滤波器愈逼近待设计的滤波器Hd(ejω)。 从频域分析: 在采样点ω=2πk,k=0,1,2,…,N-1,Ф(ω-2πk/N)=1,因此,采样点处H(ejωk)(ωk=2πk/N)与H(k)相等,逼近误差为0。 在采样点之间,H(ejω)由有限项的H(k)Ф(ω-2πk/N)之和形成。 其误差和Hd(ejω)特性的平滑程度有关,特性愈平滑的区域,误差愈小;特性曲线间断点处,误差最大。 表现形式为间断点用倾斜线取代,且间断点附近形成振荡特性,使阻衰减减小,往往不能满足技术要求。 2.3.2减小误差的方法 最直观的想法是增加采样点数,即加大N值,由于过渡带就等于采样间隔(参看图3),即 (式27) 所以加大N,可使过渡带变窄,但增加要适当,否则会增加滤波器体积与成本。 但是,增加N并不会改善滤波器的阻带衰减特性,因为Hd(ejω)是理想矩形,无论怎样增多频率采样的点数,在通、阻带交界处,幅值总是从1突变到0,会引起较大的起伏振荡。 为使逼近误差更小,和窗口法的平滑截断一样,通过在理想频率响应的不连续点的边缘上加一些过渡的抽样点,减小频带边缘的突变,也就减小了起伏振荡,增大了阻带最小衰减。 一般过渡带取一、二、三点抽样值即可得到满意结果。 如在低通设计中,不加过渡点时,阻带最小衰减为-20dB,加三个过渡点(最优设计)则可达-80dB到-95dB左右。 加过渡点的示意如图4所示。 图4理想低通滤波器增加过渡点 增加过渡点,可使阻带衰减明显提高,但付出的代价是过渡带加宽,可通过下式加大N来调整。 m=0,1,2,3…(式28) 2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 利用 MATLAB 结合 频率 取样 设计 数字 FIR 滤波器