电工电子课程设报告最新.doc

电子系统课程设计报告 FIR数字滤波器设计 专 业： 07级电子信息工程 班 级： 电信一班 学 号： XXXXXXX 姓 名： XXX 指导老师： XXXXXXXXX 同组成员： XXX 时间：2010年7月1日 一 、设计目的 1、熟悉FIR滤波器设计的基本方法。 （4） 去逼近Hd(ejw)。一般情况下，一个理想的频率响应Hd(ejw)的傅立叶反变换： （5） 所得到的理想单位脉冲响应hd(n)是一个无限长序列。对hd(n)经过适当的加权、截断处理才能得到一个所需要的有限长脉冲响应序列。对应不同的加权、截断，就有不同的窗函数。而所要寻找的滤波器脉冲响应就等于理想脉冲响应和窗函数的乘积，即： （6） 由此可见，窗函数的性质就决定了滤波器的品质。由于窗函数有一个等于M的有限长度，所以它的响应有一个主瓣，其宽度与1/M成正比，还有高度比较小的旁瓣。主瓣在Hd(ejw)中产生一个过渡带，主瓣的宽度决定过渡带的宽度，即与1/M成正比，主瓣越宽，则过渡带越宽，旁瓣产生波动。 以下是几种常用的窗函数： 1、矩形窗 （7） 2、汉宁(Hanning)窗 （8） 3、哈明(Hamming)窗 （9） 4、布莱克曼(Blackman)窗 （10） 5、凯泽(Kaiser)窗 （11） 其中I0[.]是修正过的零阶贝塞尔(Bessel)函数，β为依赖于M的参数，选择M可以产生各种过渡带宽和接近最优的阻带衰减。对于相同的M，凯泽窗可以提供不同的过渡带宽，这是其它窗所做不到的。 窗函数法设计线性相位FIR滤波器可以按如下步骤进行： （1）确定数字滤波器的性能要求。确定各临界频率{ωk}和滤波器单位脉冲长度M。 （2）根据性能要求和M值，合理地选择单位脉冲响应h(n)的奇偶对称性，从而确定理想频率响应Hd(ejw)的幅频特性和相位特性。 （3）用傅立叶反变换公式（5）求得理想单位脉冲响应hd(n)。 （4）选择适当的窗函数w(n)，根据式（6），求得所设计的FIR滤波器单位脉冲响应。 （5）用傅立叶变换求得其频率响应H(ejw)，分析它的幅频特性，若不满足要求，可适当改变窗函数形式或长度M，重复上述过程，直到得到满意的结果。 三、MATLAB实现的内容及步骤 （一）编制实验用主程序及相应子程序 1、在实验编程之前，认真复习有关FIR滤波器设计的有关知识，尤其是窗函数法的有关内容，阅读本实验原理与方法，熟悉窗函数及四种线形相位FIR滤波器的特性，掌握窗函数设计滤波器的具体步骤。 2、编制窗函数设计FIR滤波器的主程序及相应子程序。 （1）根据Hd(ejw)，利用傅立叶变换，编制产生理想单位脉冲响应hd(n)的简单子程序。 （2）根据公式（7）到公式（10）编制相应的窗函数序列的子程序。 1矩形窗程序： function w=jvxing(x); M=x; for n=1:M; w(n)=1; end 2汉宁窗程序： function w=hn(x) M=x for n=1:M w(n)=0.5*[1-cos(2*pi*n/(M+1))] end 3汉明窗程序： function w=hm(x); M=x for n=1:M w(n)=0.54-0.46*cos(2*pi*(n-1)/(M-1)) end 4布莱克曼床程序： function w=blkm(x); M=x for n=1:M w(n)=0.42-0.5*cos(2*pi*(n-1)/(M-1))+0.08*cos((4*pi*(n-1)/(M-1))) end 5凯泽窗程序： function w=kz(x); M=x; beta=4（6，8）; for n=1:M; w(n)=besseli(0,(beta*sqrt(1-((2*n/(M-1))-1)^2)))/besseli(0,beta); end 窗函数图形： （3）在上述子程序的基础上设计主程序完成线性相位 FIR滤波器的窗函数法设计。 （二）上机实验内容 1、用Hanning窗设计一个线形相位带通滤波器，其长度M=30，技术指标要求如下：