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

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

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

3、.6壬默玉.基于MATLAB勺图形用户界面的构造方式与应用J.现代电力，2002.19(1):7682.-1-信息工程学院课程设计论文7吴正茂.IIR数字滤波器的MATLA段计与实现J.长江工程职业技术学院学报，2004,21(2):2931.8陈婷，赵曜，吴刚亮。关于内模控制结构中二阶反馈滤波器的研究J.陕西工学院学报(自然科学版)，2004,20(2):22-24.9景军.一种新的双输入双输出电流模式多功能滤波器J.陕西，学院学报(自然科学版),2002,18(1):1-3.指导教师签字教研室主任签字信息工程学院课程设计论文信息工程系课程设计成绩评定表学生姓名学号：专业（班级）：电子信息科

4、学与技术课程设计题目：基于频率抽样法的FIR高通数字滤波器设计指导教师评语：成绩：指导教师：年月日信息工程学院课程设计论文摘要数字滤波是语音和图象处理、模式识别、频谱分析等应用中的一个基本处理算法，数字滤波技术是信号消噪的基本方法。根据噪声频率分量的不同，可选用具有不同滤波特性的数字滤波器。当噪声的频率高低于信号的频率时，就应该选用高通滤波器。本设计采用了频率抽样法设计的FIR高通数字滤波器，其目的是为了让中高频率的信号通过，而且利用频率抽样法的优点是可以在频域直接设计，并且适合最优化设计。FIR滤波器为有限长冲激响应滤波器，因其在线性相位特性方面具有独特的优点，因此也越来越受到广泛的重视。关

5、键词：FIR数字滤波器频率抽样高通AbstractKeyword:FIRdigitalfilterThedigitalfilteristhevoiceandimageprocessing,patternrecognition,spectrumanalysisintheapplicationofabasicalgorithm,digitalfilteringtechnologyisthebasicmethodofsignaldenoising.Accordingtodifferentfrequencycomponentsofnoise,canbeusedwithdifferentfilterin

6、gcharacteristicsofdigitalfilter.Whenthenoiseofhighfrequencybelowthesignalfrequency,shouldchooseahigh-passfilter.ThisdesignusesthefrequencysamplingdesignmethodofFIRhighpassdigitalfilter,itspurposeistomakethehighfrequencysignalthrough,andbyusingthefrequencysamplingmethodcanbedirectlyinthefrequencydoma

7、indesign,andissuitableforoptimizationdesign.FIRfilterforfiniteimpulseresponsefilters,duetoitslinearphasecharacteristichasuniqueadvantages,soitattractsmoreandmoreattention.FrequencysamplingHighpass信息工程学院课程设计论文目录1数字滤波器62频率抽样法设计的FIR高通数字滤波器72.1 频率抽样法白基本原理72.2 基于频率抽样法的FIR高通数字滤波器的设计83任务提出与方案论证103.1 设计任务与要

8、求103.2 方案设计论证与选择104程序设计与调试114.1 程序设计与说明114.2 仿真结果与分析115总结14参考文献15信息工程学院课程设计论文1数字滤波器数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号(数字序列)通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提出运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性：Yejw=XejwHejw1-1其中YbjwX(ejw历

9、别是数字滤波器的输出序列和输入序列的频域特性(或称为频谱特性)，H(ejw是数字滤波器的单位取样相应的频谱，又称为数字滤波器的频域响应。输入序列白频谱X(ejw羟过滤波后X(ejwH(ejw),因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择H(ejw)，使得滤波后的X(ejwH(ejw瞒足设计的要求，这就是数字滤波器的滤波原理。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。IIR数字滤波器的特征是，具有无限续时间冲级响应，需要用递归模型来实现，其差分方程为：NNy(n)=、aix(n-i)xbiy(n-

10、i)yy1-2信息工程学院课程设计论文2频率抽样法设计的FIR高通数字滤波器2.1频率抽样法的基本原理频率采样法是从频域出发，对给定的理想滤波器的频响Hdejw)进行N等间隔采样,即Hd(k)=Hd(ejw(为，然后以此Hd(k)作为实际FIR滤波器的频率特性采样值H(k),即令：H(k)=Hd(k)=Hd(ejw，k=0,1,2,N-1(2.1.1)由DFT定义，可以用这个N个频域的采样值H(k)来唯一确定FIR的单位脉冲响应h(n),即:.11Njknhni-IDFTHkI“HkeNNkf(2.1.(2)卜面对设计出的滤波器频率响应特性进行分析。由频域采样定理中的内插公式可以知道，利用这N

11、个频域采样值H(k刑样可以求得FIR滤波器的频率响应H(ejw),这个H(ejw腾逼近理想滤波器的频响HdfejwH(ejw和内插公式为:N1H(ejw)=ZH(k%k=0w-kN(2.1.(3)式中：中(w诞内插函数。1sinNw/2tTwesinw/2(2.1.(4)把式(4M弋入(3),化简后可得:_jN-wN11Hejw=e2'-ek=0NsinI熏咻sin-_2N(2.1.5)从式(5间以看到，在各频率采样点上，设计的滤波器，实际的频率响应严格地与理想滤波信息工程学院课程设计论文rj包、器的频率响应数值相等，即HeNk=H(k)=Hd(k)=Hd(ejW)1卬一五。但是在采样

12、点一旷<J之间的频率响应是由各采样点的加权内插函数叠加而形成的，因而有一定的逼近误差。该误差大小取决于理想频率响应的开关，理想频响特性变化越平缓，内插值越接近理想值。逼近误差越小；反之，如果采样点之间的理想频响特性变化越陡，则内插值与理想值之间的误差越大，因而在理想滤波器不连续点的两边，就会产生尖峰，而在通带和阴带就会产生波纹。2.2基于频率抽样法的FIR高通数字滤波器的设计频率抽样法设计的FIR高通数字滤波器，其目的是为了让中高频率的信号通过，而通过频率采样法的优点是可以在频域直接设计，并且适合最优化设计；缺点是采样频率只能等于2兀/N的整数倍，因而不能确保截止频率coc的自由取值，要

13、想实现自由地选择截止频率，必须增加采样点数N，但这又使计算量加大。下图是通过本实例对FIR滤波器滤波前后进行的MATLAB仿真对比。信岩波波前时域图时间Z秒图2.2.1信号滤波前时域和频域图信息工程学院课程设计论文信居更感后时地图时间，,同频率，林兹图2.2.2信号滤波后时域和频域图信息工程学院课程设计论文3任务提出与方案论证3.1 设计任务与要求利用MATLAB仿真软件系统结合双线性变换法设计一个数字巴特沃斯高通、低通和带通FIR滤波器。说明：自己确定滤波器的不同参数得到不同的结果，最后自己分析原因。3.2 方案设计论证与选择本次滤波器设计的主要方法是先设计低通模拟滤波器，然后转换为高通、带

14、通或带阻数字滤波器。对于其他如高通，带通，则通过频率变换转换为设计相应的高通，带通等。在设计的全过程的各个步骤，matlab都提供相应的工具箱函数，使得IIR数字滤波器设计变得非常简单。总的来说，我的设计思路主要有以下两种：方案一：从归一化模拟低通原型出发，先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行双线性变换，由S域变换到Z域，而得到所需类型的数字滤波器。归一化模拟低通原型模拟频_域变痪.模拟高、带通或带阻1频率抽样法二数字高、带通或带阻图3-1先频率变换再离散方案二：先进行双线性变换，将模拟低通原型滤波器变换成数字低通滤波器；然后在Z域内经数字频率变换为所需类型的数字滤波器。图

15、3-2先离散再频率变换以上两种方案都可以，我最后选择了第一种方案进行设计，即先在模拟域内经频率变换成为所需类型的模拟滤波器；然后进行频率抽样变换，由S域变到Z域，而得到所需类型的数字滤波器。因为高通，带通滤波器的设计方法是先将要设计的滤波器的技术特性指标通过频率转换成模拟低通滤波器的技术指标，再根据这些性能指标设计出低通滤波器的传递函数，为了设计程序的简洁，故选择了方案一。10信息工程学院课程设计论文4程序设计与调试4.1 程序设计与说明具体设计步骤如下：(1)确定数字高通的技术指标：Wp,Wm,Rp,Rs。(2)将高通数字滤波器的技术指标转换成高通模拟滤波器的设计指标，预畸变校正得到模拟边界

16、频率：一2,111phtan-p4.1.121zsh=-tan''''s4.1.2T2(3)通过频率变换将模拟低通转换成相应的过渡模拟高通滤波器。(4)采用频率抽样法将相应类型的过渡模拟高通滤波器转换成数字高通滤波器。4.2 仿真结果与分析将设计好的程序在matlab软件下进行仿真，得到仿真波形。4.2.1 仿真结果记录：当选择“运行”之后得到如图所示：图4.2.1GUI操作界面图11信息工程学院课程设计论文匚金宿更图4.2.2幅频特性图)Figure2HelpFileEditEiewInsertIoolsdesktopWindow图4.2.3相频特性图附部分主

17、程序如下：M=32;%所需频率采样点个数Wp=0.6*pi;%通带截止频率m=0:M/2;%阻频带上的采样点Wm=2*pi*m./(M+1);%阻带截止频率12信息工程学院课程设计论文mtr=ceil(Wp*(M+1)/(2*pi);%Ad=Wm>=Wp;Ad(mtr)=0.28;Hd=Ad.*exp(-j*0.5*M*Wm);%Hd=Hdconj(fliplr(Hd(2:M/2+1);%fliplr函数实现矩阵的左右翻转向正方向舍入ceil(3.5)=4;ceil(-3.2)=-3;构造频域采样向量H(k)conj是求复数的共轲h=real(ifft(Hd);%hw=linspace(

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