FIR数字滤波器课程设计报告要点.docx

1、FIR数字滤波器课程设计报告要点吉林建筑大学电气与电子信息工程学院数字信号处理课程设计报告设计题目： 专业班级： 学生姓名：FIR 数字滤波器的设计学号：指导教师： 设计时间：教师评语：成绩评阅教师日期1目录1、设计目的 3 2、设计内容 3 3、设计原理 33.1数字低通滤波器的设计原理 33.1.1 数字滤波器的定义和分类 3 3.1.2 数字滤波器的优点 3 3.1.3 FIR 滤波器基本原理 43.2 变换方法的原理 7 4、设计步骤 8 5、数字低通滤波器 MATLAB 编程及幅频特性曲线 95.1 MATLAB 语言编程 95.2幅频特性曲线 106、总结 11 7、参考文献 13

2、2-3一、设计目的课程设计是理论学习的延伸，是掌握所学知识的一种重要手段，对于贯彻 理论联系实际、提高学习质量、塑造自身能力等于有特殊作用。本次课程设计 一方面通过 MATLAB 仿真设计内容，使我们加深对理论知识的理解，同时增强其 逻辑思维能力，另一方面对课堂所学理论知识作一个总结和补充二、设计内容(1)设计一线性相位 FIR 数字低通滤波器，截止频率 c 0.2 ，过渡带宽度 0.4 ,阻带衰减 As 30 dB 。(2)设计一线性相位 FIR 数字低通滤波器，截止频率 c 0.2， 过渡带宽度，阻 带 衰0.4减 As 50 dB 。三、设计原理3.1 数字低通滤波器的设计原理3.1.1

3、 数字滤波器的定义和分类数字滤波器是指完成信号滤波处理功能的，用有限精度算法实现的离散时 间线性非时变系统，其输入是一组数字量，其输出是经过变换的另一组数字量。 因此，数字滤波器本身既可以是用数字硬件装配成的一台完成给定运算的专用 的数字计算机，也可以将所需要的运算编成程序，让通用计算机来执行。从数字滤波器的单位冲击响应来看，可以分为两大类:有限冲击响应(FIR) 数字滤波器和无限冲击响应(IIR)数字滤波器。滤波器按功能上分可以分为低通滤波器(LPF)、高通滤波器(HPF)、带通滤波器(BPF)、带阻滤波器(BSF)4。3.1.2 数字滤波器的优点相对于模拟滤波器，数字滤波器没有漂移，能够处

4、理低频信号，频率响应 特性可做成非常接近于理想的特性，且精度可以达到很高，容易集成等，这些 优势决定了数字滤波器的应用将会越来越广泛。同时 DSP 处理器(Digital Signal Processor)的出现和 FPGA(FieldProgrammable Gate Array)的迅速发 展也促进了数字滤波器的发展，并为数字滤波器的硬件实现提供了更多的选择。数字滤波器具有以下显著优点:精度高:模拟电路中元件精度很难达到 10 ，以上，而数字系统 17 位字长3-5就可以达到 10 精度。因此在一些精度要求很高的滤波系统中，就必须采用数 字滤波器来实现。灵活性大:数字滤波器的性能主要取决于乘

5、法器的各系数，而这些系数是 存放在系数存储器中的，只要改变存储器中存放的系数，就可以得到不同的系 统，这些都比改变模拟滤波器系统的特性要容易和方便的多，因而具有很大的 灵活性。可靠性高:因为数字系统只有两个电平信号:1”和“0，受噪声及环境条 件的影响小，而模拟滤波器各个参数都有一定的温度系数，易受温度、振动、 电磁感应等影响。并且数字滤波器多采用大规模集成电路，如用 CPLD 或 FPGA 来实现，也可以用专用的 DSP 处理器来实现，这些大规模集成电路的故障率远 比众多分立元件构成的模拟系统的故障率低。易于大规模集成:因为数字部件具有高度的规范性，便于大规模集成，大 规模生产，且数字滤波电

6、路主要工作在截止或饱和状态，对电路参数要求不严 格。因此产品的成品率高，价格也日趋降低。相对于模拟滤波器，数字滤波器 在体积、重量和性能方面的优势己越来越明显。比如在用一些用模拟网络做的 低频滤波器中，网络的电感和电容的数值会大到惊人的程度，甚至不能很好地 实现，这时候若采用数字滤波器则方便的多。并行处理：数字滤波器的另外一个最大优点就是可以实现并行处理，比如 数字滤波器可采用 DSP 处理器来实现并行处理。TI 公司的 TMS320C5000 系列的 DSP 芯片采用 8 条指令并行处理的结构，时钟频率为 100MHZ 的 DSP 芯片，可高 达 100MIPs(即每秒执行百万条指令)。3.

7、1.3 FIR 滤波器基本原理 FIR 数字滤波器的特点及结构在数字信号处理应用中往往需要设计线性相位的滤波器，FIR 滤波器在保 证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性。FIR 滤波 器不断地对输入样本 x(n)延时后，再作乘法累加算法，将滤波结果 y(n)输出， 因此，FIR 实际上是一种乘法累加运算。在数字滤波器中，FIR 滤波器的最主要的特点是没有反馈回路，故不存在 不稳定的问题，同时，可以在幅度特性是随意设置的同时，保证精确的线性相4位。稳定和线性相位特性是 FIR 滤波器的突出优点。另外，它还有以下特点： 设计方式是线性的;硬件容易实现;滤波器过渡过程具有有限区

8、间;相对 IIR 滤波器而言，阶次较高，其延迟也要比同样性能的 IIR 滤波器大得多。3FIR 数字滤波器系统的传递函数为:(3.1)通过反 z 变换,数字滤波器的差分方程为:(3.2)由此得到系统的差分方程:(3.3)由上式可以得出如下图 3.1 所示的直接型结构，这种结构又可以称为卷积型结 构。将转置理论应用于图 3.1 可以得到转置直接型结构。将式中的系统函数 H(z)分解成若干一阶和二阶多项式的连乘积: (1.4)(3.4)则可构成如图 1.1 所示的级联型结构。其中 一阶节;为为二阶节。每个一阶节、二阶节可用图 3.2 所示的直接型结构实现。当 M1= M2时，即得到图 3.3 所示

9、的具体结构。这种结构的每一节都便于控制零点，在需要控制传输零点时可以采用。 但是它所需要的系数 a 比直接型的 h(n)多，所需要的乘法运算也比直接型多。 在对滤波器计算时间没有特殊要求的时候可以采用这种形式。若需要严格考虑 滤波器的计算时间则需要折衷它们的优点和缺点来设计。这在算法设计时候要 使用软件编辑环境来计算运行的时间问题。通常 FIR 的计算时间都较长。很多 时候我们需要牺牲时间来获得想要得到的滤波器功能。5图 3.1 FIR 滤波器直接型结构图图 3.3级联型具体结构2 FIR 滤波器的优点：可以在幅度特性随意设计的同时，保证精确、严格 的线性相位；由于 FIR 滤波器的单位脉冲

10、h(n)是有限长序列，因此 FIR 滤波器 没有不稳定的问题；由于 FIR 滤波器一般为非递归结构，因此，在有限运算下 不会出现递归型结构中的极限振荡等不稳定现象误差较小；FIR 滤波器可以采 用 FFT 算法实现，从而提高了运算效率。3 FIR 和 IIR 滤波器的比较在很多实际应用中如语音和音频信号处理中，数字滤波器来实现选频功能。 因此，指标的形式应为频域中的幅度和相位响应。在通带中，通常希望具有线 性相位响应。在 FIR 滤波器中可以得到精确的线性相位。在 IIR 滤波器中通带 的相位是不可能得到的,因此主要考虑幅度指标。IIR 数字滤波器的设计和模拟 滤波器的设计有着紧密的联系，通常

11、要设计出适当地模拟滤波器，再通过一定 的频带变换把它转换成为所需要的数字 IIR 滤波器。此外，任何数字信号处理 系统中也还不可避免地用到模拟滤波器，比如 A/D 变换器前的抗混叠滤波器及 D/A 转换后的平缓滤波器，因此模拟滤波器设计也是数字信号处理中应当掌握 的技术。6从性能上来说，IIR数字滤波器传递函数包括零点和极点两组可调因素，对 极点的唯一限制是在单位圆内。因此可用较低的阶数获得高的选择性，所用的 存储单元少、计算量小、效率高。但是这个高效率是以相位的非线性为代价的。 选择性越好，则相位非线性越严重。FIR滤波器传递函数的极点是固定在原点， 是不能动的，它只能靠改变零点位置来改变它

12、的性能，所以要达到高的选择性， 必须用高的阶数，对于同样的滤波器设计指标，FIR滤波器所要求的阶数可能比 IIR滤波器高5-10倍，结果成本高信号延时也较大，如果按线性相位要求来说， 则IIR滤波器就必须加全通网络进行相位校正，同样大大增加了滤波器的阶数和 复杂性。而FIR滤波器却可以得到严格的线性相位。从结构上看，IIR 滤波器必须采用递归结构来配置极点，并保证极点位置 在单位圆内。由于有限字长效应，运算过程中将对系数进行舍入处理，引起极 点的偏移，这种情况有时会造成稳定性问题，甚至造成寄生振荡。相反，FIR 滤波器只要采用非递归结构，不论在理论上还是实际的有限精度运算中都不存 在稳定性问题

13、，因此造成的频率特性误差也较小。此外 FIR 滤波器可以采用快 速傅立叶变换算法，在相同的阶数条件下运算速度可以快的多。3.2 变换方法的原理FIR 滤波器目前常用的设计方法有窗函数法和频率采样法，窗函数法是从 时域进行设计，而频率采样法是从频域进行设计。窗函数法由于简单、物理意 义清晰，因而得到了较为广泛的应用。窗函数法设计的基本思想是：首先根据 技术指标要求，选取合适的阶数 N 和窗函数的类型 w(n)，使其幅频特性逼近 理想滤波器幅频特性。其次，因为理想滤波器的 hd(n)是无限长的，所以需要 对 hd(n) 进行截断，数学上称这种方法为窗函数法。 简而言之，用窗函数法 设计 FIR 滤

14、波器是在时域进行的，先用傅里叶变换求出理想滤波器单位抽样相 应 hd(n)，然后加时间窗 w(n)对其进行截断，以求得 FIR 滤波器的单位抽样响 应 h(n)。四、设计步骤窗函数旁瓣峰值幅度/dB过渡带宽阻带最小衰减/dB矩形窗-134/-12三角形窗-258/-257汉宁窗-318/-44哈明窗-418/-53布莱克曼窗-5712/-74凯塞窗-5712/-80表 4.1 6 种窗函数的基本参数（1）根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口 长度 N。先按照阻带衰减选择窗函数的类型。原则是在保证阻带衰减满足要求 的情况下，尽量选择主瓣窄的窗函数。然后根据过渡带宽度估计

15、窗口长度 N。 待求滤波器的过渡带宽度 B,近似等于窗函数的主瓣宽度，且近似与窗口长度 Ni 成反比。矩形窗的 A=4，哈明窗的 A=8 等。（2）构造希望逼近的频率响应函数 Hd ( ejw），即所谓的“标准窗函数法”，就是选择 Hd ( ejw）为线性相位理想滤波器（理想低通、理想高通、理想带通、理想带阻）。以低通滤波器为例， Hdg ( w)应满足：理想滤波器的截止频率 wc 近似于最终设计的 FIRDF 的过渡带的中心频率点，幅 度函数衰减一半。所以如果设计指标给定通带边界频率和阻带边界频率 wp 和 ws ， 一般取（3）计算 hd ( n ) 。如果给出待求滤波器的频响函数为 Hd

16、 ( e 响应用下式求出:jw) ，那么单位脉冲如果 Hd (ejw) 较复杂，或者不能用封闭公式表示，则不能用上式求出8hd ( n ) 。我们可以对 Hd (ejw) 从 w=0 到 w=2 采样 M 点,采样值为，k=0,1,2,，M-1，进行 M 点 IDFT(IFFT),得到：根据频域采样理论，与应满足如下关系：因此，如果 M 选的较大，可以保证在窗口内 hdM ( n ) 有效逼近 hd ( n ) 。对（公式 1）式给出的线性相位低通滤波器作为 Hd ( e求出单脉冲响应 hd ( n ) ：为保证线性相位特性，=(N-1)/2.（4）加窗得到设计结果： h ( n ) hd (

17、 n ) w( n ) 。jw) ，由（公式 2）式五、 数字低通滤波器 MATLAB 编程及幅频特性曲线 5.1 MATLAB 语言编程（1） delta=0.4*pi; Wc=0.2*pi; As=30; N=ceil(8*pi/delta)+1; win=bartlett(N); h=fir1(N-1,Wc/pi,win); omega=linspace(0,pi,512); mag=freqz(h,1,omega); magdb=20*log10(abs(mag); plot(omega/pi,magdb); axis(0 1 -100 0); grid; xlabel(归一化功率);

18、 ylabel(幅度/dB);9（2） delta=0.4*pi; Wc=0.2*pi; As=50; N=ceil(8*pi/delta)+1; win=hamming(N); h=fir1(N-1,Wc/pi,win); omega=linspace(0,pi,512); mag=freqz(h,1,omega); magdb=20*log10(abs(mag); plot(omega/pi,magdb); axis(0 1 -100 0); grid; xlabel(归一化功率); ylabel(幅度/dB);5.2幅频特性曲线图（1）10图（2）六、 总结本次 MATLAB 课程设计主

19、要任务是完成 FIR 数字滤波器的设计，对我来说这 个题目还是很有挑战性的，因为自己对 MATLAB 中有关滤波器的设计知识了解较 少，但正是这种有挑战性的题目才能提高自己的能力，才有研究价值，入手这 个题目后我查阅了相关的资料，也从网上获得了不少有关 MATLAB 设计滤波器的 资料，加上自己之间对 MATLAB 有一定的了解，因而设计思路渐渐明朗，经过自 己的不断尝试和探索，终于弄明白了 FIR 滤波器的工作原理，通过查阅相关资 料和研究 MATLAB 中提供的帮助信息，我也明白了 FIR 滤波器设计相关的一些函 数的使用方法，并用它们来设计 FIR 滤波器，最终完成了题目。在设计过程中，

20、我也遇到了很多不懂得地方，程序经常出现错误，尤其是 在利用所设计的滤波器对模拟出来的数字序列滤波时，出现很多错误，但经过 自己的不断努力和尝试，最终还是解决了问题。同时，在书写报告的过程中， 更加深入的学会了对 Word 的使用，其中遇到最大的困难就是如何在 Word 中添 加设计 FIR 数字滤波器过程中计算所需的公式，虽然书写起来很麻烦，但是最11终通过不断的探索和努力还是完成了。七、参考文献1数字信号处理原理与实践（第 2 版），方勇主编，清华大学出版社， 2010 年2MATLAB 信号处理，刘波、文忠、曾涯编著，电子工业出版社,2006.01 3数字信号处理实验指导书，吉林建筑大学电气与电子实验中心， 2013.094数字信号处理（第三版），高西全、丁玉美编著，西安电子科技大学出 版社，2008 年 8 月5 Matlab 基础及其应用教程 ，周开利、邓春辉编著，北京大学出版社, 200712