dsp中的fir滤波器论文

《dsp中的fir滤波器论文》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、安徽工业大学研究生考试试卷评分考试科目：DSP技术阅卷人：陈德昱专业:检测技术与自动化装置学号：20100273姓名：王普注意事项1．考前研究生将上述项目填写清楚2．字迹要清楚、保持卷面清洁3．教师将试卷、答案一起送研究生学院归档二〇一一年六月24日13摘要随着信息与数字技术的发展，数字信号处理已经成为当今极其重要而学科与技术领域之一。它在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理的基本方法中，通常会涉及到变换、滤波、频谱分析、调制解调和编码解码等处理。其中滤波是应用非常广泛的一个环节，数字滤

2、波器的理论和相关设计也一直都是人们研究的重点之一。FIR滤波器的是非递归的，稳定性好，精度高；更重要的是，FIR滤波器在满足幅频响应要求的同时，可以获得严格的线性相位特征。因此，它在高保真的信号处理，如数字音频、图像处理、数据传输和生物医学等领域得到广泛应用。关键词：FIR滤波MATHLAB窗函数法CCS131．设计要求1）设计一FIR低通滤波器，实现对语音信号的滤波。2）用语音信号去检验该滤波器，查看滤波效果。根据老师的要求，我们自己拟定了滤波器的指标：Fs=8000hz，Wp=1500Hz，Ws=2000Hz，通带波纹为0.01，阻带波纹为0.1，N=

3、37。2.滤波器的设计原理说明2.1数字滤波器的设计原理数字滤波器的设计问题就是寻找一组系数ai和bi，使得其性能在某种意义上逼近所要求的特性。如果在s平面上去逼近，就得到模拟滤波器，如果在z平面上去逼近，则得到数字滤波器。数字滤波是将输入的信号序列，按规定的算法进行处理，从而得到所期望的输出序列。一个线性位移不变系统的输出序列y(n)和输入序列x(n)之间的关系，应满足常系数线性差分方程：（2.1.1）x(n)为输入序列，y(n)为输出序列，ai、bi为滤波器系数，N为滤波器的阶数。2.2FIR滤波器的基本结构在式2.1.1中，若所有的ai均为0，则得F

4、IR滤波器的差分方程：（2.2.1）对式（2.2.1）进行z变换，可得FIR滤波器的传递函数：（2.2.2）由此可得到FIR滤波器的结构如图1所示。13FIR滤波器的单位冲击响应h(n)是一个有限长序列。若h(n)为实数，且满足偶对称或奇对称的条件，即h(n)=h(N-1-n)或h(n)=-h(N-1-n)，则FIR滤波器具有线性相位特性。2.3FIR滤波器的主要特点1）单位冲击响应只有有限项；2）FIR滤波器无反馈回路，是一种无条件稳定系统；3）FIR滤波器可以设计成具有线性相位特性。3.FIR滤波器的设计方法3.1FIR滤波器设计FIR滤波器的设计问题

5、在于寻求一系统函数，使其响应逼近滤波器要求的理想频率响应。3.2窗函数设计的基本方法从时域出发，设计h(n)逼近理想hd(n)，设理想滤波器的单位脉冲响应为hd(n)，则有所求得的一般是无限长的，且是非因果的。要想得到一个因果的有限长的滤波器h(n)，最直接的方法是截断，或者说用一个窗口函数对进行加窗处理，即。因此，应选择合适的窗函数。3.3滤波器的算法实现FIR滤波器的输出表达式为式中，为滤波器系数；x(n)表示滤波器在n时刻的输入；y(n)为n时刻的输出。它的基本算法是一种乘法-累加运算，即不断地输入样本x(n)，经过延时后，再进行乘法-累加，最后输出

6、滤波结果y(n)。131）线性缓冲区法线性缓冲区法又称延迟线法。其特点：·对于N级的FIR滤波器，在数据存储器中开辟一个N单元的缓冲区（滑窗），用来存放最新的N个输入样本；·从最老样本开始取数，每取一个样本后，将此样本向下移位；·读完最后一个样本后，输入最新样本存入缓冲区的顶部。2）循环缓冲区法本次设计的FIR滤波器所采用的就是循环缓冲区法。循环缓冲区法的特点如下：·对于N级FIR滤波器，在数据存储器中开辟一个N单元的缓冲区（滑窗），用来存放最新的N个输入样本；·从最新样本开始取数；·读完最后一个样本（最老样本）后，输入最新样本来代替最老样本，而其他数据位

7、置不变；·用片内BK（循环缓冲区长度）寄存器对缓冲区进行间接寻址，使循环缓冲区地址首尾相邻。4.FIR滤波器的MATLAB实现MATLAB是一种功能强、效率高、便于进行科学和工程计算的交互式软件包，它集数值分析、矩阵运算、信号处理和图形显示于一体，为用户提供了方便、友好的界面环境。MATLAB中的工具箱（Toolbox）包含了许多实用程序。它提供了多种FIR滤波器设计方法。4.1用fir1函数设计FIR滤波器fir1函数用来设计标准频率响应的基于窗函数的FIR滤波器，可实现加窗线性相位FIR数字滤波器的设计。具体语法如下：b=fir1(n，Wn)b=fir

8、1(n，Wn，‘ftype’)b=fir1(n，Wn，Window