雍鹏亮基于MatlabSimulink的线性相位FIR滤波器设计副本

1、目录1.绪论31.1数字滤波器定义31.2数字滤波器的分类31.3数字滤波器的技术指标32.FIR数字滤波器的条件和特点52.1 线性相位FIR数字滤波器52.2第一类线性相位对h(n)的约束条件62.3线性相位FIR滤波器的幅度特性Hg(w)的特点72.4线性相位FIR滤波器的零极点分布特点83.FIR数字滤波器的MATLAB设计93.1 通过窗函数法设计FIR滤波器93.1.1 窗函数设计原理93.1.2 典型窗函数介绍11用窗函数法设计FIR录波器的步骤13通过 M文件实FIR数字语音低通滤波器的设计133.2通过等波纹最佳逼近法设计FIR数字滤波器173.2.1 remez和remez

2、ord函数滤波器指标17通remez和remezord函数实现FIR数字语音低通滤波器的设计183.3通过FDATool实现FIR数字滤波器的设计及分析194.利用SIMULINK对FIR滤波器进行仿真234.1Simulink模块库简介234.2搭建语音信号FIR低通滤波仿真模型23总结27致谢281.绪论当前，通信技术高速发展，业务范围不断扩大，人们对产品的需求迅速增长，滤波器在这些产品电路中就扮演着重要的角色，当然数字滤波器更影响着人们生活的方方面面。数字滤波器既可以是用数字硬件装配成德一台完成给定运算的专用数字计算机，也可将所需的运算编成程序，让通用计算机来执行。MATlAB是“矩阵实

3、验室”（MATrix LABoratoy）的缩写，它是由美国Mathworks公司于1984年正式推出的，是一种以矩阵运算为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求。MATLAB是功能强大的科学及工程计算软件，它不但具有以矩阵计算为基础的强大数学计算和分析功能，而且还具有丰富的可视化图形表现功能和方便的程序计算能力。MATLAB的应用领域极为广泛，除数学计算和分析外，还被广泛地应用于自动控制、系统仿真、数字信号处理、图形图像分析、数理统计、人工智能、虚拟现实技术、通信工程、金融系统等领域，因此，MATLAB是面向21世纪的计算机程序设计及科学计算语言。1.1数字滤波器定义所谓数

4、字滤波器，是指输入、输出均为数字信号，通过数字运算处理改变入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。因此数字滤波的概念和模拟滤波相同，只是信号的形式和实现滤波的方法不同。正因为数字滤波器通过数值运算实现滤波，所以数字滤波器的处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器不能实现的特殊滤波功能。1.2数字滤波器的分类数字滤波器按照不同的分类方法，有许多种类，总体上可分为两大类。一类是经典滤波器，即一般的线性系统滤波器。另一类即所谓的现代滤波器，它的理论建立在随机信号处理的理论基础上，利用了随机信号内部的统计特性对信号进行滤波。从功能上可以

5、分为低通、高通、带通和带阻滤波器，根据数字滤波器冲激响应的时域特性，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器1.3数字滤波器的技术指标1)通带边界频率2)阻带截止频率3)通带内允许最大衰减4)阻带内允许最小衰减2.FIR数字滤波器的条件和特点 与IIR滤波器相比，FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到具有严格的线性相位特性。本章用N表示FIR滤波器单位脉冲响应h(n)的长度，其系统函数H(z)为：H(z)是的N-1次多项式，所以它在z平面内有N-1个零点，在原点z=0处有N-1个重极点，因此H(z)永远稳定。稳定性和相位线性是F

6、IR滤波器最为特出的优点。2.1 线性相位FIR数字滤波器 对于长度为N的h(n)，频率响应函数为： （2-1） （2-2）式中为幅度特性，为相位特性，值得注意的是，不同于，为实函数，其值可正可负，总是正。线性相位FIR滤波器是指是w的线性函数，即 为常数 (2-3)如果 满足下式 是起始相位 (2-4)严格的说此时不具有线性相位特性，但以上两种情况都满足群延时是一常数，即：也称这种情况为线性相位。一般称满足(2-3)式是第一类线性相位，满足(2-4)式是第二类线性相位。本论文的FIR滤波器设计仅采用第一类线性相位。2.2第一类线性相位对h(n)的约束条件 线性相位FIR滤波器的时域约束条件是

7、指满足相位线性时，对h(n)的约束条件。第一类线性相位的FIR滤波器的相位函数为，联立(2-1)式和 (2-2)式可得到： （2-5）由式（2.1.5）可得： （2-6）将（7.2.2）中两式相除可得到：即：移相得： （2-7） 函数h(n)sinw(n-)若能关于(N-1)/2奇对称，则为（2-7）的一组解，显然sinw(n-)关于n=奇对称，若取=(N-1)/2，则可取h(n)关于n=(N-1)/2成偶对称，即满足如下条件： （2-8）由以上推论可知：如果要求单位脉冲响应为h(n),长度为N的FIR滤波器具有第一类相位特性，只需满足即可。2.3线性相位FIR滤波器的幅度特性Hg(w)的特点

8、实质上，幅度特性特点就是FIR滤波器线性相位的频率域约束条件，将时域约束条件代入（2-1）式中，即可推导出FIR线性相位的Hg(w)的约束条件，当N取奇偶不同时，Hg(w)的约束条件不同，本论文仅针对第一类相位特性进行讨论，所以只有两种情况。为了书写方面，引入一下两个参数： ， 表示取不大于的最大整数，显然，当为奇数时=，当为偶数是-=。情况1：，为奇数。所以 （2-8）由于关于三点偶对称，因此情况1可以实现低通、高通、带通、带阻滤波器。情况2：，为偶数仿照情况1可推导： （2-9）此时关于偶对称，关于奇对称，此时不能实现高通和带阻滤波器。2.4线性相位FIR滤波器的零极点分布特点 将代入上式

9、 （2-10）由（2-10）可以分析出，若是的零点，则也是其零点，又因为是实序列，所以的零点必成共轭对出现，所以和也是其零点。(z)是的N-1次多项式，所以它在z平面内有N-1个零点，在原点z=0处有N-1个重极点，因此H(z)永远稳定。3.FIR数字滤波器的matlab设计3.1 通过窗函数法设计FIR滤波器3.1.1 窗函数设计原理设希望逼近的滤波器频率响应函数为,，其单位脉冲响应为。 （3-1） （3-2）如果由已知的得出，经过Z变换可得到滤波器的系统函数，但理想器的幅频响应在边界频率处往往是不连续的，因而是无限时宽的，且是非因果序列。例如，线性相位理想低通滤波器的频率响应函数为： (3

10、-3)其单位脉冲响应为： （3-4）由上式看到，理想低通滤波器的单位脉冲响应是无限长，且是非因果序列。的波形如图3-1（a）所示。为了构造一个长度为N的第一类线性相位的FIR滤波器，只有截取一段，并保证截取的一段关于偶对称。设截取的一段用表示，即： (3-5)式中为一矩形序列，长度为N，波形如图3-1（b）所示。由该图可知，当取值为时，截取的一段关于偶对称，保证所设计的滤波器具有线性相位。我们所设计的滤波器单位脉冲响应为，长度为N，其系统函数为，这样用一个有限序列代替，肯定会引起误差，表现在时域就是吉布斯效应。 在(3-5)中，就是起对无限时宽非因果序列截断作用，所以形象地将称为窗函数，将称为

11、对的加窗处理。图3-1 窗函数设计法的时域波形图加窗后的幅频特性如图3-2 图3-2 幅频响应函数曲线3.1.2 典型窗函数介绍 1）矩形窗其频率响应为：其幅度响应为：2）汉宁窗升余弦窗其频率响应为：其幅度响应为：3）布莱克曼窗其频率响应为：其幅度响应为：表3-1 5种窗函数基本参数类型窗函数的旁瓣峰值过渡带宽度加窗后滤波器的阻带最小衰减rectwin-134p/N-21bartlet三角-258p/N-25hanning-318p/N-44hamming-418p/N-53blackman-5712p/N-743.1.3用窗函数法设计FIR录波器的步骤 根据窗函数设计FIR滤波器的步骤如下：

12、（1） 根据阻带衰减指标，选择合适的窗函数类型，再根据过渡带指标，估算出窗口长度N。具体参数参考表3.1.2.1.（2） 构造希望逼近的频率响应函数。很据所设计滤波器（低通、高通、带通、带阻）及通带边界频率和阻带截止频率构造希望逼近的频率响应函数。一般想滤波器的截止频率取: ,为通带边界频率，为阻带截止频率。（3） 计算。（4） 加窗得到设计结果：。3.1.4通过 M文件实FIR数字语音低通滤波器的设计 实际设计时一般使用matlab工具箱函数。可调用工具箱函数fir1实现窗函数设计。Fir1是利用窗函数法设计线性相位fir数字滤波器的工具箱函数，实现线性相位FIR数字滤波器的标准窗函数设计。

13、这里所指的“标准”，是指设计低通、高通、带通和带阻FIR滤波器时，分别表示相应的线性相位理想低通、高通、带通和带阻滤波器的频率响 应函数。因而将所设计的滤波器的频率响应函数称为标准频率响应。Fir1的调用格式及功能： hn=fir1(m,wc),返回6dB截止频率为的M阶滤波器系数向量hn，默认选用哈明窗。滤波器单位脉冲响应与的关系为：其中为对归一化的数字角频率，。当=wcl,wcu得到的是带通滤波器，其-6dB通带为。hn=fir(M,wc,ftype)，可设计高通和带阻滤波器。当ftype=high时，设计高通FIR滤波器；当ftype=stop，且=wcl,wcu时，设计带阻FIR滤波器

14、。hn=fir(M,wc,window)，可设计指定窗函数类型的FIR滤波器系数向量。缺省时默认为哈明窗。利用matlab设计FIR语音低通滤波器，技术指标要求：通带03000HZ，阻带3500HZ最小衰减50dB，采样频率为fs=22050HZ。实现程序：clear all;fs=22050Hz;f=0:fs-1;wp=3000*2*pi;ws=3500*2*pi; deltw=(ws-wp)/fs;N0=ceil(6.6*pi/deltw);N=N0+mod(N0+1,2);hn=fir1(N-1,(wp+ws)/(2*fs*pi);Hw=fft(hn,fs);abs_Hw=20*log1

15、0(abs(Hw);phase=angle(Hw);subplot(3,1,1);stem(0:N-1,hn);title(单位脉冲响应);xlabel(n);ylabel(h(n);subplot(3,1,2);plot(f,abs_Hw);title(幅频特性);xlabel(f);ylabel(20lg(|H(ejw)|);axis(0 11025 -150 50);grid onsubplot(3,1,3);plot(f,phase);title(相频特性);xlabel(f);ylabel(phase);axis(0 11025 -4 4);grid on图3-3 窗函数法实现的FI

16、R滤波器的响应特性有一个音乐文件，其信号频谱为：3-5 无噪音的音乐频谱 插入白噪声后的频谱为：3-5 加噪音的音乐频谱 利用所设计的FIR数字语音低通滤波器滤波后的频谱为：3-6 滤波后的音乐频谱实现程序：%-读取文件-y,fs,nbyte=wavread(c:因为爱情); M=length(y); f=(1:M)*(fs/M); %-加噪处理 - noise=0.01*randn(M,2); si=y+noise; Y=fft(y); abs_Yw=1/fs*abs(Y); Si=fft(si); abs_Siw=1/fs*abs(Si);%-滤波器设-wp=3000*2*pi;ws=35

17、00*2*pi; deltw=(ws-wp)/fs;N0=ceil(6.6*pi/deltw);N=N0+mod(N0+1,2);hn=fir1(N-1,(wp+ws)/(2*fs*pi);%-滤波处理- x=fftfilt(hn,si);Xw=fft(x);abs_Xw=1/fs*abs(Xw); %-画图部分- subplot(3,1,1); plot(f,abs_Yw);title(无噪声的信号频谱);xlabel(f/HZ);ylabel(幅度); subplot(3,1,2); plot(f, abs_Siw); title(无噪声的信号频谱);xlabel(f/HZ);ylabel

18、(幅度); subplot(3,1,3); plot(f,abs_Xw); title(处理后的信号频谱); xlabel(f/HZ);ylabel(幅度);%-结束-3.2通过等波纹最佳逼近法设计FIR数字滤波器等波纹最佳逼近法是一种优化设计方法，使最大误差最小化，并在整个逼近频段上均匀分布。用等文博最佳逼近法设计的FIR滤波器的幅频响应在整个通带和阻带都是等波纹的。而且可以分别控制通带和阻带的波纹幅度。最佳逼近是指在滤波器长度给定的条件下，使加权误差最小化。由于这种设计方法使最大误差均与分布，所以设计的滤波器性价比最高。阶数相同时，这种设计是最大逼近误差最小，即通带最大衰减最小，最带最小衰

19、减最大；指标相同时，滤波器的阶数最低。3.2.1 remez和remezord函数滤波器指标1）remez和remezord函数remez函数实现线性相位FIR数字滤波器的等文波最佳逼近设计。其调用格式为：hn=remez(M,f,m,w) 调用结果返回单位脉冲响应hn。M：滤波器的阶数，f和m给出希望逼近的特性。W为误差加权量。remezord函数可根据逼近指标估算出等文波最佳逼近FIR滤波器的最低阶数M，误差加权量w和归一化边界频率向量f。其返回参数供remez函数调用。其格式为：M,f0,m0,w=remezord(f,m,rip,Fs)2）滤波器指标低通滤波器设计指标 逼近通带：0，通

20、带最大衰减：dB；逼近阻带：， ，阻带最小衰减：dB，remezord调用参数：高通滤波器设计指标 逼近通带：，通带最大衰减：dB；逼近阻带：0， ，阻带最小衰减：dB，remezord调用参数：带通滤波器设计指标 逼近通带：，通带最大衰减：dB；逼近阻带：0， ，,，阻带最小衰减：dB，remezord调用参数：带阻滤波器设计指标 逼近通带：0， ，,，通带最大衰减：dB；逼近阻带：，阻带最小衰减：dB，remezord调用参数：3.2.2通remez和remezord函数实现FIR数字语音低通滤波器的设计利用matlab设计FIR语音低通滤波器，技术指标要求：通带03000HZ最大衰减1d

21、B，阻带3500HZ最小衰减50dB，采样频率为fs=22050HZ。Fs=22025;fps=3000,3500;m=1,0;rp=1;rs=50;dat1=(10(rp/20)-1)/(10(rp/20)+1);dat2=10(-rs/20);rip=dat1,dat2;M0,f0,m0,w=remezord(fps,m,rip,Fs); M=M0+mod(M0+1,2); hn=remez(M,f0,m0,w); n=0:M; Hw=fft(hn,Fs); abs_Hw=20*log10(abs(Hw); %-画图部分- subplot(2,1,1); stem(n,hn); title

22、(单位脉冲响应 );xlabel(n);ylabel(h(n); subplot(2,1,2); f=0:Fs-1; plot(f,abs_Hw);grid on; title(幅频响应);xlabel(f/HZ);ylabel(20*lg(|Hw|);图3-7 remez和remezord函数实现的滤波器响应特性分析图3-3和图3-7易得，在满足相同技术指标时，用等文波最佳逼近法设计fIR滤波器的阶数要低于用窗函数设计的阶数。3.3通过FDATool实现FIR数字滤波器的设计及分析FDATool（Filter Design and Analysis Tool）是一个强大的数字滤波器设计与分析

23、工具。利用它可方便地以设计出满足各种性能指标的滤波器，并可以查看滤波器各种分析图形。等待滤波器设计满意后，还可以把其系数直接导出为matlab变量，文本文件或C语言头文件。在命令窗中运行FDATool，启动滤波器分析工具，界面如下图3.3.1所示。图3-8 FDATool界面下面利用FDATool设计FIR语音低通滤波器，技术指标要求：通带03000HZ最大衰减1dB，阻带3500HZ最小衰减50dB，采样频率为fs=22050HZ。FDATool界面设计参数如图3-9图3-9 FDATool界面参数设置设置滤波器为低通滤波器 图3-10 滤波器类型设置用窗函数法设计FIR数字滤波器图3-11

24、设计方法设置采用哈明窗函数图3-12 窗函数类型设置根据过渡带带宽估算滤波器阶数M=N-1，图3-13 滤波器阶数设置采样频率，截止频率图3-14 滤波器采样频率设置图3-15、图3-16、图3-17是利用FDATool所设计的FIR数字语音低通滤波器的响应特性。图3-15 单位脉冲响应图3-16 幅频特性图3-17 相频特性由此可见利用FDATool工具箱可以很方便的设计我们所需的FIR数字滤波器，并可以生成.fda格式的工程文件，供simulink仿真使用。4.利用Simulink对FIR滤波器进行仿真 Simulink是用来对动态系统进行建模、仿真和分析的软件包。使用Simulink来建

25、模、分析和仿真各种动态系统，将是一种非常轻松的事情。它提供了一种图形化的交互环境，利用鼠标拖动大方法便能迅速地建立起系统框图模型。4.1Simulink模块库简介 Simulink模块库极其丰富，在本论文中不可能作全面介绍，下面只介绍在4.2章节中搭建语音信号FIR低通滤波仿真模型所需要的模块库。1） 信号源模块库脉冲输入（Discrete Impulse）读工作库模块（From Workspace）2）数学运算模块加法器（Add）3）信号处理模块数字滤波器设计（Digital Filter Design）4）输出模块示波器（Scope）频谱仪（Spectrum）扬声器（To Audio De

26、vice）4.2搭建语音信号FIR低通滤波仿真模型 通过4.1章节所介绍的模块库搭建如下仿真模型：图4-1 仿真模型图通过示波器Scope可观察语音信号的时域波形图4-2图4-2 音乐波形图通过频谱仪可观察到Spectrum Scope1、Spectrum Scope2、Spectrum Scope3、可观察信号加噪前后及处理后的频谱图。 图4-4 Spectrum Scope1所示信号频谱 图4-5 Spectrum Scope2所示信号频谱图4-6 Spectrum Scope3所示信号频谱总结滤波是信号处理中最基本有极为重要的技术，利用滤波器技术可以从复杂的信号中提取出所需要的信号，抑制

27、不需要的信号。FIR滤波器的优点FIR滤波器能严格做到线性相位或群延时为常数，而IIR滤波器只能逼近线性相位。FIR滤波器是全零点型滤波器，总是稳定的，不会因滤波运算的舍入误差而产生极限环振荡现象对同样幅度相应的滤波器，用FIR滤波器实现比用IIR滤波器实现需要较高的阶数，多达510倍。当滤波器的特性要求较高时，用FIR滤波器来实现，滤波过程需要较多的计算时间Matlab软件作为强大的计算软件， Matlab在滤波器设计方面还提供了编写图形用户界面的功能。所谓图形用户界面, 简称为GUI, 是由各种图形对象, 如图形窗口菜单按钮、文本框等构建的用户界面，它是用户和计算机之间交流的工具。Matlab 将所有GUl 支持的用户控件都集成在这个环境中并提供界面外观、属性和行为响应方式的设置方法,随着版本的提高,这种能力还会不断加强。而且具有强大的绘图功能,可以轻松