学士学位论文—-电子初级技能培训论文电子秤

1、郑州航空工业管理学院电子信息系统仿真课程设计2014级电子信息工程专业班级题目 基于 Hamming 窗的 FIR 数字滤波器设计姓名 学 号二O 一六年土二月二十五日一、 MATLAB 软件简介MATLAB 是矩阵实验室( Matrix Laboratory )之意。除具备卓越 的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可 视化建模仿真和实时控制等功能。 MATLAB 的基本数据单位是矩阵， 它的指令表达式与数学 ,工程中常用的形式十分相似 ,故用 MATLAB 来解算问题要比用 C,FORTRAN 等语言完相同的事情简捷得多。MATLAB 拥有 数百个内部函 数的主包和三十

2、几种工具包 (Toolbox).工具包又可以分为功能性工具包和学科工具包.功能工具包 用来扩充 MATLAB 的符号计算 ,可视化建模仿真 ,文字处理及实时控 制等功能 .学科工具包是专业性比较强的工具包,控制工具包 ,信号处理工具包 ,通信工具包等都属于此类 . 开放性使 MATLAB 广受用户 欢迎.除内部函数外 ,所有 MATLAB 主包文件和各种工具包都是可读 可修改的文件 ,用户通过对源程序的修改或加入自己编写程序构造新 的专用工具包。 语言特点一种语言之所以能如此迅速地普及，显示 出如此旺盛的生命力， 是由于它有着不同于其他语言的特点， 正如同 FORTRAN和C等高级语言使人们摆

3、脱了需要直接对计算机硬件资源 进行操作一样，被称作为第四代计算机语言的 MATLAB ，利用其丰 富的函数资源， 使编程人员从繁琐的程序代码中解放出来。 MATLAB 最突出的特点就是简洁。 MATLAB 用更直观的，符合人们思维习惯 的代码，代替了 C 和 FORTRAN 语言的冗长代码。 MATLAB 给用户 带来的是最直观，最简洁的程序开发环境。二、理论分析1. FIR数字滤波器结构特征FIR数字滤波器的冲激响应只能延续一定的时间，N阶因果有限冲激M响应滤波器的差分方程表达式为：y(k) h(i)x(k i)。其卷积和是一个有限和，可以直接计算y(k),所以其涉及的基本操作就是简单的乘和

4、 加，其计算仅需要所给出的初始样本值及所要求的相关样本值。通过 精确设计，有限冲激响应滤波器在整个频率范围内均能提供精确的线 性相位，由于系统的单位脉冲序列为有限序列，当输入有限时，输出也必然为有限，这样其稳定性总可以独立于滤波器系数之外。因此，在很多情况下，有限冲激响应滤波器成为首选，只要确定能满足要求的转移 序列或者脉冲响应的常数，就可以准确地设计出满足要的FIR数字滤波器。通常所采用的设计方法主要有窗函数、频率采样法和等波纹最 佳逼近法等，其中窗函数法是从时域进行设计的，其算法简单、物理意 义清晰，因此得到了较为广泛的应用。2. 算法研究为了建立一个具有线性相位和稳定的非递归特性的有限脉

5、冲响应 滤波器即FIR滤波器,要考虑两个方面：一是使用有限长的单位取样响 应来逼近理想低通，二是单位取样响应对(N -1)/2对称,保证线性相位 根据FIR数字滤波器输入输出关系的差分方程描述，有n 4y(x) * h(i)x(n -1)i上式反映了有限列长的单位取样响应N _1 滤波网络传输函数为：H(z)八hd(n)Zn =0利用傅立叶变换得单位取样响应 hd( n)为1hd(n)N j jn血(Wc(n 一 a)-H d(e)e d = 一 (n _a)W是截止频率错误！未找到引用源。，a=(N_1)/2从而保证线性相位。 在窗函数法设计中使用有限列长的 错误！未找到引用源。h(n)逼近

6、 错误！未找到引用源。hd(n)，用Hamming窗函数错误！未找到引用 源。w(n)将错误！未找到引用源。hd(n)截断，进行加权处理后可得 错 误！未找到引用源。h(n)=hd(n)w(n)。N 4此时频率响应函数为：Hd(ej )八h( n)e-j nn=0Hamming窗函数wh(n)二0.54 -0.46遇# Rn(n)。其中错误！未找到引用源。rn (n)是幅度为1，长度为N的矩形序列。3. 窗函数法的原理可以从时域或频域出发来设计 FIR滤波器，从频域出发的方法称为频率采样设计法，从时域出发的方法称为窗函数法，也称为傅立叶级 数法,这种方法应用较广泛。窗函数法的设计思想是按照所要

7、求的理 想滤波器频率响应 错误！未找到引用源。Hd(ej ),设计一个FIR滤波N 4器,使之频率响应H Q)八h(n)n来逼近错误！未找到引用源。n=0Hd(e)。因为设计是在时域中进行的，需先由错误！未找到引用源Hd(e)的傅立叶反变换导出序列 错误！未找到引用源。hd( n)，即1 严jjccnhd(n) =2? JH d(e )ed 由于错误！未找到引用源。Hd(ej )是矩形频率特性，有 错误！未找到 引用源。hd(n)是一无限长的序列，且是非因果的,而要设计的FIR滤波 器的冲激响应序列是有限长的，所以要用有限长的序列 错误！未找到 引用源。h(n)来逼近无限长的序列hd(n)，最

8、有效的方法是截断 错误! 未找到引用源。hd(n)，或者说用一个有限长度的窗口函数 错误！未 找到引用源。w(n)序列来截取错误！未找到引用源。hd(n)，即h(n) = hd(n)w( n)按照复卷积公式，在时域中的乘积关系可表示成在频域中的周期性卷积关系，即可得所设计的FIR滤波器的频率响应i八1 兀H(ej )亍(诃(-心其中,错误！未找到引用源。WC)为截断窗函数的频率特性。由此可见,实际的FIR数字滤波器的频率响应 错误！未找到引用源。H (ej )逼近理想滤波器频率响应 错误！未找到引用源。Hd(ej )的好坏， 完全取决于窗函数的频率特性 错误！未找到引用源。W(.)。4. 窗函

9、数法的设计步骤(1) 给定所要求的频率响应 错误！未找到引用源。Hd (ej );(2) 由错误！未找到引用源。Hd(ej )的傅立叶反变换求出 错误！未找 到引用源。hd(n);(3) 根据过渡带宽及阻带最小衰减，选择窗的形状并估计窗口长度 N;(4) 计算滤波器的单位脉冲响应h(n)二hd (n)w(n)(5) 验算技术指标是否满足要求。(1)窗谱主辨尽可能窄，以获得较陡的过渡带；(2)尽量减少窗谱的最大旁瓣的相对幅度，使能量尽量集中于主辨，进而增加阻带的衰减。给定的滤波器指标参数一般为通带截止频率 错误！未找到引用源。、 阻截止频率错误！未找到引用源。.s、实际通带波动错误！未找到引 用

10、源。Rp和最小阻带衰减 错误！未找到引用源。As。窗函数设计的 经验公式为归一化过渡带，二一 错误！未找到引用源。滤波器阶数2 二M二乓-7.95错误！未找到引用源。14.36M当A错误！未找到引用 源。 50时，错误！未找到引用 源。匕=0.1102(As -8.7)；当 21 W错误！未找到引用源。 50 时， =0.5842(As -21)0.4 0.07886(As -21)三、设计内容1. 低通滤波器的设计指标。通带截止频率 错误！未找到引用源。=02：，阻带截止频率 错误！未 找到引用源。叫=0.财,最小阻带衰减错误！未找到引用源。ASOdB, 抽样频率为2000Hz2. 设计分析

11、由以上参数可得过度带宽 trwidth = ws-wp , wc = (ws + wp)/2。考虑选择 Hamming 窗函数，则列长 N 二 6.64*pi/trwidth ,hd(n)=sin(wc*(n-(N-1)/2)/(pi*(n-(N-1)/2), h(n)=hd (n)*wh(n)。利用 MATLAB 编程时，先计算 N、trwidth、hd(n)和h(n),并利用 MATLAB 提供的相关工具计算该低通滤波器的幅度特性值 , 最终显 示 Hamming 窗频谱、实际单位取样响应频谱、理想单位取样响应频 谱和该 FIR 滤波器幅度特性曲线。3编写 Matlab 程序(1) 实现

12、FIR 滤波器的主程序 :wp=0.2*pi;ws=0.4*pi;%给出通带和阻带的截止频率trwidth=ws-wp;%得出过渡带宽度N=ceil(6.64*pi/trwidth)+1;%滤波器的长度n=0:1:N-1;wc=(ws+wp)/2;% 计算理想低通截止频率hd=ideallp(wc,N);% 计算理想低通单位冲激响应w_ham=(hamming(N);%求海明函数h=hd.*w_ham;%应用窗函数加权后的有限长单位冲激响应序列db,mag,pha,grd,w=freqz_m(h,1); %幅检查值响应、相位响应和群延迟响应 delta_w=2*pi/1000;Ap=-(min

13、(db(1:1:wp/delta_w+1);%检查通带最大衰减Ar=-round(max(db(ws/delta_w+1:1:501);%检查阻带最小衰减subplot(221);stem( n,hd);title(理想单位脉冲响应 hd( n)subplot(222);plot(n,w_ham);title(海明窗 w(n) subplot(223);stem( n,h);title(实际单位脉冲响应h(n)subplot(224);plot(w/pi,db);title( 幅度响应 (dB) n=0:1023;Fs=2000;%抽样频率dt=1/Fs;t=n*dt;Tp=1024*dt;f

14、1=450;f2=600;x=sin(2*pi*f1*t)+sin(2*pi*f2*t); y=filter(w_ham,1,x);figure(2);subplot(2,1,1); plot(t,x);title(输入信号 x(t); xlabel(t/s); ylabel(x(t);axis(0,Tp/6,min(x),max(x);subplot(2,1,2); plot(t,y);title(输出信号 y(t); xlabel(t/s);ylabel(y(t);axis(0,Tp/6,mi n(y),max(y);(2)计算理想低通的单位冲激响应的程序：fun cti on hd=id

15、eallp(wc,N)%理想线性相位低通滤波器计算%hd=ideallp(wc ,n)% hd=0N1之间的理想脉冲响应%可。=截止频率(弧度)%N=理想线性相位滤波器的长度%tao=(N-1)/2;n=0:1:(N-1);m=n-tao+eps;%加一个极小数eps,以避免hd的分母为零hd=si n( wc*m)./(pi*m);检查幅值响应、相位响应和群延迟响应fun cti on db,mag,pha,grd,w = freqz_m(b,a);H,w = freqz(b,a,1000,whole);H = (H(1:1:501); w = (w(1:1:501);mag = abs(H

16、);db = 20*log10(mag+eps)/max(mag);pha = an gle(H);grd = grpdelay(b,a,w);输出结果如下图所示：理想单位胁冲响应hdCn)J0.20.10*0 1 _实际单位脉冲响应h(n)海明窗0110203040 W 203040NameValue 0+0301 55536735ApArAsMinMax0.03010.03015555535367673535由上图可以看到：通带纹波为0.0301dB,最小阻带衰减为55dB。满足设计指标的要求把滤波器的阶数人为的设置成75 (将 N二ceil(6.64*pi/trwidth)+1;改为N=

17、75)理想单位脉冲响应hd(n)海明窗w(n)0.5Name ValueMinMax-Ap0.02840.02840.0284-Ar585858As535353-M676767-N757575由上图可以看到：通带纹波为0.0284dB,最小阻带衰减为 58dB把滤波器的阶数设置为15:hCTG1由話/胪0 30 20.1O51O15理想单位脉冲响应hd(n)0.30.20.1O-0.1头琢单位脉冲响应h(n)O51O1500.010 020 030.040 050 060 070,08谕入信弓刈9t/sValue1.S322MinMax1.63221.63221616535367671515由

18、上图可以看到：通带纹波为 1.6322dB,最小阻带衰减为16dB。不能满足指标要求 综上所述：当滤波器的阶数增加时，可以减小通带纹波，增加最小阻 带衰减。但同时计算量也随之增大，成本也要同样增加。用 MATLAB中的FilterDesignandAnalysis工具进行同样的设计，设置阶数为35,低通，窗口选择 hamming，归一化频率-0.3。幅频特性曲线(un_DCU - QGnELId00,10,20.30+0.50.6D.7Q.&0 9Normalized Frequency g rad/sampler相频特性曲线0510152025M35Samples4JP3亡auv冲击响应错误！未找到引用源。*1-2*1.5-1-0.50 OS 1Real Part1s O45七足严6匚一Be E 一零极点图四、总结通过本次课程设计，使我更加扎实的掌握了有关数字滤波器方面的知 识，在这次课程设计中虽然遇到了很多的困难，但是经过一次又一次的查阅资料，一次一次的认真思考，终于找出了原因，通过这次课程 设计，了解知识还是需要通过实践