基于频率抽样法的FIR高通数字滤波器的设计

1、信息工程系课程设计报告书题目:基于频率抽样法的 FIR高通数字滤波器设计专 业： 电子信息科学与技术班 级：学 号：学生姓名：指导教师：年月 日信息工程学院课程设计论文信息工程系课程设计任务书学号学生姓名专业(班级)电子信息科学与技术设计题目其于频率抽样法的FIR高通数字滤波器设计设计技术参数所需频率采样点个数M=32;通带截止频率 Wp=0.6*pi; 阻频带上的采样点 m=0:M/2;阻带截止频率 Wm=2*pi*m.(M+1);设 计 要 求要求数字咼通滤波器的通带截止频率为0.6 rad通带衰减不大于 IdB阻带截止频率为 0.96 rad阻带衰减不小于 25dB工作量共15页工 作

2、计 划2011.8.29-2011.8.31选题，并查阅相关资料；2011.9.1 -2011.9.5在MatIab环境下进行仿真并给出仿真结果；2011.9.6 -2011.9.9制作用户界面图形(GUI);2011.9.10-2011.9.11完成初步设计报告。参 考 资 料1 丁玉美，高西全.数字信号处理M .西安：西安电子科技大学出版社，2002.2 John G Proakis.数字信号处理M .方艳梅，刘永清，译.4版.北京：电子工 业出版社，2007.3 程佩青.数字信号处理教程M .北京：清华大学出版社，2007.4 吴镇扬.数字信号处理M .北京:高等教育出版社，2004.5

3、 胡广书.数字信号处理导论M .北 京:清华大学出版社，2006. 壬默玉.基于 MATLAB勺图形用户界面的构造方式与应用J.现代电力，2002 . 19(1) : 76 82.7 吴正茂. R数字滤波器的MATLA设计与实现J.长江工程职业技术学院学 报，2004, 21(2) : 29 31 .8 陈婷，赵曜，吴刚亮。关于内模控制结构中二阶反馈滤波器的研究J.陕西工学院学报(自然科学版)，2004, 20(2) : 22 24.9 景军一种新的双输入双输出电流模式多功能滤波器J.陕西工学院学报(自然科学版),2002 , 18(1) : I 3.指导教师签字教研室主任签字-5 -信息工程

4、系课程设计成绩评定表学生姓名学号：专业（班级）：电子信息科学与技术课程设计题目：基于频率抽样法的FIR高通数字滤波器设计指导教师评语：成绩： 指导教师： 年 月 日摘要数字滤波是语音和图象处理、模式识别、频谱分析等应用中的一个基本处理算法，数字滤波技术是信号消噪的基本方法。根据噪声频率分量的不同，可选用具有不同滤波特性的数字滤波器。当噪声的频率高低于信号的频率时，就应该选用高通滤波器。 本设计采用了频率抽样法设计的FIR高通数字滤波器，其目的是为了让中高频率的信号通过，而且利用频率抽样法的优点是可以在频域直接设计，并且适合最优化设计。FIR滤波器为有限长冲激响应滤波器，因其在线性相位特性方面具

5、有独特的优点，因此也越来越受到广泛的重视。关键词：FIR数字滤波器频率抽样高通AbStraCtThe digitalfilter isthe voice andimage PrOCeSS ing,Patter nrecognition, SPeCtrUman alysisin the applicationof a basic algorithm,digitalfilteri ng tech no IOgy is the basic method Of Sig nal denoising. ACCOrd ing to different frequency compOnentS of nois

6、e, Can be USed With different filtering CharaCteriStiCS of digital filter. When the noise of high frequency below the Signal frequency, should choose a high-pass filter. ThiS design USeS the frequency SamPIing design method of FIR high PaSS digital filter, its PUrPOSe is to make the high frequency S

7、ignal through, and by USing the frequency SamPIing method Can be directly in the frequency domain design, and is SUitabIe for OPtimiZation design. FIR filter for finite impulse response filters, due to its Iinear PhaSe CharaCteriStiC has UniqUe advantages, so it attracts more and more atte ntio n.Ke

8、y word: FIR digital filterFreqUe ncy SamPIi ngHigh PaSS信息工程学院课程设计论文目录1数字滤波器 62频率抽样法设计的FlR高通数字滤波器 72.1频率抽样法的基本原理 72.2基于频率抽样法的FIR高通数字滤波器的设计 83任务提出与方案论证 103.1设计任务与要求 103.2方案设计论证与选择 104程序设计与调试 114.1程序设计与说明 114.2仿真结果与分析 115总结 14参考文献 15181数字滤波器数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号(数字序列)

9、通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提出运算规 则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性：YejW =X ejw H ejw1 -1其中YejW、X ejw分别是数字滤波器的输出序列和输入序列的频域特性(或称为频谱特性)，HejW是数字滤波器的单位取样相应的频谱，又称为数字滤波器的频域响应。输入序列的频谱X ejw经过滤波后XejWHejW，因此，只要按照输入信号频谱的特点和处 理信号的目的，适当选择HejW，使得滤波后的XejW

10、H ejW满足设计的要求，这就是数 字滤波器的滤波原理。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。IIR数字滤波器的特征是，具有无限续时间冲级响应，需要用递归模型来实现，其差分方程为：NNy(n)7 ajX(n i) ' 0y(n i)z(12)2频率抽样法设计的FIR高通数字滤波器2.1频率抽样法的基本原理频率采样法是从频域出发，对给定的理想滤波器的频响 HdejW 进 行N等间 隔 采 样,即Hdk=HdejW w±*,然后以此H d ( k)作为实际F IR滤波器的频率特性采样值H (

11、k),即令：H k =Hdk =HdejW, k = 0,1,2, -,N-1(2.1.1)由DFT定义，可以用这个N个频域的采样值 H k来唯一确定FIR的单位脉冲响应hn ,即:h n =IDFT IH k IyHk eNjknN k =O(2.1.2)F面对设计出的滤波器频率响应特性进行分析。由频域采样定理中的内插公式可以知道，利用这N个频域采样值H k同样可以求得FlR滤波器的频率响应HejW ,这个HejW将逼近理想滤波器的频响HdejW。H ejw的内插公式为:N JH (ejw )= H (kk=0WkN(2.1.3)式中：W是内插函数。W=N1 Sin Nw/2 -weSin

12、w/2(2.1.4)把式4代入3 ,化简后可得:HejW =e 2、丄ek Z0 N怡HkJWrsin.r7:二 k(2.1.5)从式5可以看到，在各频率采样点上，设计的滤波器，实际的频率响应严格地与理想滤波J je Nk)之间的频率响应是由各采样点的加权内插函数叠加而形成的,器的频率响应数值相等，即H=H k = Hdk =Hdej .：|2-。但是在采样点WkN因而有一定的逼近误差。 该误内插值越接近理想值。逼近差大小取决于理想频率响应的开关，理想频响特性变化越平缓，则内插值与理想值之间的误差误差越小；反之，如果采样点之间的理想频响特性变化越陡， 越大，因而在理想滤波器不连续点的两边，就会

13、产生尖峰，而在通带和阴带就会产生波纹。2.2基于频率抽样法的FIR高通数字滤波器的设计频率抽样法设计的FIR高通数字滤波器，其目的是为了让中高频率的信号通过，而通过频率采样法的优点是可以在频域直接设计，并且适合最优化设计；缺点是采样频率只能等于2 /N的整数倍，因而不能确保截止频率 C的自由取值，要想实现自由地选择截止频率， 必须增加采样点数 N ,但这又使计算量加大。 下图是通过本实例对 FIR滤波器滤波前后进行 的MATLAB仿真对比。时间楸图2.2.1信号滤波前时域和频域图图222 信号滤波后时域和频域图3任务提出与方案论证3.1设计任务与要求利用MATLAB仿真软件系统结合双线性变换法

14、设计一个数字巴特沃斯高通、低通 和带通FIR滤波器。说明：自己确定滤波器的不同参数得到不同的结果，最后自己分 析原因。3.2方案设计论证与选择本次滤波器设计的主要方法是先设计低通模拟滤波器，然后转换为高通、带通或带阻数字滤波器。对于其他如高通，带通，则通过频率变换转换为设计相应的高通，带通等。在设 计的全过程的各个步骤，matlab都提供相应的工具箱函数，使得IIR数字滤波器设计变得非常简单。总的来说，我的设计思路主要有以下两种：方案一：从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟 滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。模拟频归一化模拟

15、低通原型模拟高、带通或带阻域变换频率抽样法：数字高、带通或带阻图3-1先频率变换再离散方案二：先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在Z域内经数字频率变换为所需类型的数字滤波器。图3-2先离散再频率变换以上两种方案都可以，我最后选择了第一种方案进行设计，即先在模拟域内经频率变 换成为所需类型的模拟滤波器；然后进行频率抽样变换，由S域变换到Z域，而得到所需类型的数字滤波器。因为高通，带通滤波器的设计方法是先将要设计的滤波器的技术特性指标 通过频率转换成模拟低通滤波器的技术指标，再根据这些性能指标设计出低通滤波器的传递函数，为了设计程序的简洁，故选择了方案一。4程序设计与

16、调试4.1程序设计与说明具体设计步骤如下：(1) 确定数字高通的技术指标：WP，Wm, RP，RS。(2) 将高通数字滤波器的技术指标转换成高通模拟滤波器的设计指标，预畸变校 正得到模拟边界频率：4.1.1214.1.2tan ST2(3) 通过频率变换将模拟低通转换成相应的过渡模拟高通滤波器。(4) 采用频率抽样法将相应类型的过渡模拟高通滤波器转换成数字高通滤波器。4.2仿真结果与分析将设计好的程序在 matlab软件下进行仿真，得到仿真波形。4.2.1仿真结果记录：当选择“运行”之后得到如图所示：图421 GUl 操作界面图图422幅频特性图6.5.O3.2O.F i.gur R 2EiI

17、e Edi tew Ijisert IClOlS Besktop Jtindow I* ¾¾n® < Bo图423相频特性图附部分主程序如下：M=32;%所需频率采样点个数Wp=0.6*pi;%通带截止频率m=0:M/2;%阻频带上的采样点Wm=2*pi*m./(M+1);%阻带截止频率mtr=ceil(Wp*(M+1)(2*pi)% Ad=Wm>=Wp;Ad(mtr)=0.28;Hd=Ad.*exp(-j*0.5*M*Wm);%Hd=Hd conj(fliplr(Hd(2:M/2+1);%fliplr函数实现矩阵的左右翻转h=real(ifft(Hd

18、);%hw=li nspace(O,pi,1OOO);%H=freqz(h,1,w);% figure(1)Plot(Wpi,20*log10(abs(H)% XIabeI('wpi');ylabel(' 幅度(dB)'); axis(0 1 -50 0);向正方向舍入 ceil(3.5)=4;ceil(-3.2)=-3;构造频域采样向量H(k)conj是求复数的共轭(n) =IDFTH(k)用于产生0,pi之间的1000点行矢量 滤波器的幅频特性图参数分别是归一化频率与幅值figure(2)plot ( w/pi,angle(H);仿真结果分析：由仿真结果可知，设计出的数字高通的幅频响应均符合自身的特性，从滤波器特性曲线可以看出，所设计的数字高通滤波器的通带截止频率为0.2 rad,通带衰减不大于带截止频率为0.3 rad,阻带衰减不小于25dB ,符合本次设计的要求。1dB,阻5总结MATLB仿真结果验证了用频率采样法设计F IR高通数字滤波器这一数字信号处理中的基本理论，有助于我们深入理解并掌握这一重要的FIR滤波器设计方法。需要强调的是，频率采样法是从频域出发直接设计滤波器的，而窗函数法是从时域出发设计滤波器的，两种设计方法各有优缺点。窗函数法设计 FIR数字滤波器是傅里叶变换的典 型运用，