一级项目滤波器

1、信息与电气工程学院电子信息工程CDIO一级项目（2013/2014学年第二学期）题 目 ： 数字滤波器设计 专业班级 ： 信息1101 学生姓名 ： 学 号： 指导教师 ： 设计周数 ： 第 周 设计成绩 ： 2014年月日 数字滤波器设计1. 项目设计目的1.1 掌握数字滤波器原理1.2 掌握数字滤波器设计思路1.3 掌握数字滤波器设计流程1.4 熟练使用MATLAB软件1.5 熟练掌握运用MATLAB设计数字滤波器2. 项目设计正文2.1 项目设计内容 运用MATLAB 软件设计出数字滤波器，利用窗函数法、频率取样法及优化设计方法使之能正常处理信号。2.2 软件介绍 MATLAB 

2、;(Matrix Laboratory)为美国Mathworks公司1983年首次推出的一套高性能的数值分析和计算软件，其功能不断扩充，版本不断升级，1992年推出划时代的4.0版，1993年推出了可以配合Microsoft Windous使用的微机版，95年4.2版，97年5.0版，99年5.3版，5.X版无论是界面还是内容都有长足的进展，其帮助信息采用超文本格式和PDF格式，可以方便的浏览。至2001年6月推出6.1版，2002年6月推出6.5版，继而推出6.5.1版, 2004年7月MATLAB7和Simulink6.0被推出，目前的最新版本为7.1版。&#

3、160;      MATLAB将矩阵运算、数值分析、图形处理、编程技术结合在一起，为用户提供了一个强有力的科学及工程问题的分析计算和程序设计工具，它还提供了专业水平的符号计算、文字处理、可视化建模仿真和实时控制等功能，是具有全部语言功能和特征的新一代软件开发平台。   MATLAB 已发展成为适合众多学科，多种工作平台、功能强大的大型软件。在欧美等国家的高校，MATLAB已成为线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的基本教学工具。成为攻读学位的本科、硕士、博士生必须掌

4、握的基本技能。在设计研究单位和工业开发部门，MATLAB被广泛的应用于研究和解决各种具体问题。在中国，MATLAB也已日益受到重视，短时间内就将盛行起来，因为无论哪个学科或工程领域都可以从MATLAB中找到合适的功能。2.3FIR数字滤波器设计原理  FIR 滤波器具有严格的相位特性，对于信号处理和数据传输是很重要的。目前 FIR滤波器的设计方法主要有三种：窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。常用的是窗函数法和切比雪夫等波纹逼近的最优化设计方法。本实验中的窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求高的时候是比较灵活方便的。

5、60;  如果 FIR 滤波器的 h(n)为实数, 而且满足以下任意条件,滤波器就具有准确的线性相位：  第一种：偶对称，h(n)=h(N-1-n)，   ()=-(N-1)/2    第二种：奇对称，h(n)=-h(N-1-n),  ()=-(N-1)/2+pi/2  对称中心在n=(N-1)/2处 。  FIR滤波器的基本结构可以理解为一个分节的延时线，把每一节的输出加权累加，可得到滤波器的输出。FIR

6、滤波器的冲激响应h(n)是有限长的，数学上M阶FIR滤波器可以表示为：FIR滤波器的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数的问题，设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。2.4 用窗函数设FIR 滤波器的基本方法 窗函数设计法是一种通过截短和计权的方法使无限长非因果序列成为有限长脉冲响应序列的设计方法。通常在设计滤波器之前，应该先根据具体的工程应用确定滤波器的技术指标。在大多数实际应用中，数字滤波器常常被用来实现选频操作，所以指标的形式一般为在频域中以分贝值给出的相对幅度响应和相位响应。用窗函数法设计FIR滤波器的步骤如下：（1）根据过渡带宽及阻带衰减要求，

7、选择窗函数的类型并估计窗口长度N（或阶数M=N-1）。窗函数类型可根据最小阻带衰减AS独立选择，因为窗口长度N对最小阻带衰减AS没有影响。在确定窗函数类型以后，可根据过渡带宽小于给定指标确定所拟用的窗函数的窗口长度N。设待求滤波器的过渡带宽为，它与窗口长度N近似成反比。窗函数类型确定后，其计算公式也确定了，不过这些公式是近似的，得出的窗口长度还要在计算中逐步修正。原则是在保证阻带衰减满足要求的情况下，尽量选择较小的N。在N和窗函数类型确定后，即可调用MATLAB中的窗函数求出窗函数wd(n)。（2）根据待求滤波器的理想频率响应求出理想单位脉冲响应hd(n)。如果给出待求滤波器的频率响应为Hd(

8、ej)，则理想的单位脉冲响应可以用下面的傅里叶反变换式求出：在一般情况下，hd(n)是不能用封闭公式表示的，需要采用数值方法表示。从=0到=2采样N点，采用离散傅里叶反变换(IDFT)即可求出。（3）计算滤波器的单位脉冲响应h(n)。它是理想单位脉冲响应和窗函数的乘积，即h(n)=hd(n)·wd(n)，在MATLAB中用点乘命令表示为h=hd·wd。（4）验算技术指标是否满足要求。为了计算数字滤波器在频域中的特性，可调用freqz子程序，如果不满足要求，可根据具体情况，调整窗函数类型或长度，直到满足要求为止。使用窗函数法设计时要满足以下两个条件：（1）窗谱主瓣尽可能地窄，

9、以获得较陡的过渡带；（2）尽量减少窗谱的最大旁瓣的相对幅度，也就是使能量尽量集中于主瓣，减小峰肩和纹波，进而增加阻带的衰减。根据工程经验，给定的滤波器指标参数一般为通带截止频率p、阻带截止频率s、实际通带波动Rp和最小阻带衰减As。窗函数设计的经验公式为：在实际工程中常用的窗函数有五种，即矩形窗、三角窗、汉宁窗、海明窗和凯泽窗。这些窗函数在MATLAB中分别用boxcar、triang、hanning、hamming、kaiser实现，它们之间的性能比较如表1所示。表1 5种窗函数性能比较窗类型旁瓣峰值主瓣峰值最小阻带衰减矩形窗13dB4/M21dB三角窗25dB8/M25dB汉宁窗31dB8

10、/M44dB海明窗41dB8/M53dB凯泽窗57dB12/M74dB2.5窗函数的MATLAB实现在窗函数法的Matlab实现中，程序中经常使用的函数有fir1和kaiserord。 程序中fir1函数的用法：b=fir1(n,Wn,ftype,window)n为滤波器的阶数Wn为滤波器的截止频率，它是一个0到1的数。如果Wn是一个含有两个数的向量，则函数返回一个带通滤波器ftype为滤波器的类型，ftype=high时，设计的是高通滤波器；ftype=stop时，设计的是带阻滤波器；没有此参数时，设计的是低通滤波器window为指定的窗函数，矩形窗为boxcar(n)，汉宁窗为hannin

11、g(n)，海明窗为hamming(n)，布莱克曼窗为blackman(n)，凯撒窗为kaiser(n,beta)，没有此参数时，默认为hamming窗程序中kaiserord函数的用法：n,Wn,beta,ftype=kaiserord(f,a,dev,Fs)f是一个向量，为设计滤波器过渡带的起始点和结束点a是一个向量，指定频率段的幅度值dev是一个向量，长度和a相同，为各个通带和阻带内容许的幅度最大误差n为能够满足要求的滤波器的最小阶数Wn为滤波器的截止频率ftype为根据待设计滤波器的要求得到的滤波器的类型高通滤波器是容许高频信号通过、但减弱（或减少）频率低于截止频率信号通过的滤波器。对于

12、不同滤波器而言,每个频率的信号的减弱程度不同。它有时被称为低频剪切滤波器；在音频应用中也使用低音消除滤波器或者噪声滤波器。低通滤波器与高通滤波器特性恰恰相反。利用窗函数法设计低通滤波器设计要求：使用hamming窗，采样频率2000Hz通带截频0.1，阻带截频0.17通带衰减小于等于0.1dB，阻带衰减大于等于50dB程序参见附录二中的1-（1）利用窗函数法设计低通滤波器图2-1 窗函数法设计低通滤波器的增益响应从参考程序及图2-1可以得到所设计出滤波器的参数如下：滤波器的采样频率为2000Hz，滤波器的阶数为266滤波器的通带截频0.1 ，阻带截频0.17 ，过渡带宽0.07 通带衰减为0.

13、019dB，阻带衰减为53dB 对比设计要求与所设计出滤波器的参数可知，其各项参数均满足设计指标，所设计出的滤波器即为设计所要求的滤波器。图2-2 信号滤波前的时域图和频域图图2-3 信号滤波后的时域图和频域图从图2-2和图2-3的图像中可以看到：输入信号是由两个不同频率的正弦信号叠加而成，信号频域图中位于滤波器通带内的频率分量保留了下来，位于滤波器阻带内的频率分量被滤除，滤波器的效果符合设计要求。（2） 利用窗函数法设计带通滤波器设计要求：使用Kaiser窗，采样频率8000Hz通带截频0.325与0.5525，阻带截频0.25与0.6025阻带衰减大于等于40dB，通带和阻带波纹0.01程

14、序参见附录二中的1-（2）利用窗函数法设计带通滤波器图2-4 窗函数法设计带通滤波器的增益响应从参考程序及图2-4可以得到所设计出滤波器的参数如下：滤波器的采样频率为8000Hz，滤波器的阶数为90滤波器的通带截频0.325与0.5525，阻带截频0.25与0.6025，过渡带宽0.075与0.05阻带衰减为40dB，通带和阻带的波纹均为0.01对比设计要求与所设计出滤波器的参数可知，其各项参数均满足设计指标，所设计出的滤波器即为设计所要求的滤波器。图2-5 信号滤波前的时域图和频域图图2-6 信号滤波后的时域图和频域图从图2-5和图2-6的图像中可以看到：输入信号是由四个不同频率的正弦信号叠

15、加而成，信号频域图中位于滤波器通带内的频率分量保留了下来，位于滤波器阻带内的频率分量被滤除，滤波器的效果符合设计要求。 (3) 利用窗函数法设计多通带滤波器设计要求：使用Kaiser窗，采样频率200Hz通带截频0.2、0.4、0.7、0.8阻带截频0.1、0.5、0.6、0.9阻带衰减大于等于30dB，通带和阻带波纹0.01程序参见附录二中的1-（3）利用窗函数法设计多通带滤波器图2-7 窗函数法设计多通带滤波器的增益响应从参考程序及图2-7可以得到所设计出滤波器的参数如下：滤波器的采样频率为200Hz，滤波器的阶数为46滤波器的通带截频0.2、0.4、0.7、0.8，阻带截频0.1 、0.

16、5、0.6、0.9，过渡带宽均为0.1阻带衰减为38dB，通带和阻带的波纹均为0.01对比设计要求与所设计出滤波器的参数可知，其各项参数均满足设计指标，所设计出的滤波器即为设计所要求的滤波器。图2-8 信号滤波前的时域图和频域图图2-9 信号滤波后的时域图和频域图从图2-8和图2-9的图像中可以看到：输入信号是由六个不同频率的正弦信号叠加而成，信号频域图中位于滤波器通带内的频率分量保留了下来，位于滤波器阻带内的频率分量被滤除，滤波器的效果符合设计要求。3. 程序代码 FIR低通 clear;close all;fs=4000;x1=wavread('xinhao.wav');w

17、p=2*pi*1000/fs;ws=2*pi*1200/fs;Rp=1;Rs=100;wdelta=ws-wp;n=ceil(100-7.95)/(2.286*wdelta);bta=0.1102*(100-8.7);wn=(wp+ws)/2; b,a=fir1(n-1,wn/pi,kaiser(n,bta); figure(1);subplot(311);freqz(b,a,512,fs);grid onxlabel('频率/Hz')ylabel('频率响应幅度')title('FIR低通滤波器频率响应'); f2=filtfilt(b,a,x

18、1); %f2为滤波后输出figure(2);subplot(211);plot(x1);title('FIR低通滤波器滤波前的时域波形');subplot(212);plot(f2);title('FIR低通滤波器滤波后的时域波形');sound(f2,20000)F0=fft(f2,1024); %F0为滤波后fftf=fs*(0:511)/1024; figure(3);y2=fft(x1,1024); %y2为原始信号的fftsubplot(2,1,1);plot(f,abs(y2(1:512);axis(0,700,0,100);title('

19、FIR低通滤波器滤波前的频谱');xlabel('频率/HZ');ylabel('幅值');subplot(2,1,2);plot(f,abs(F0(1:512);axis(0,700,0,100);title('FIR低通滤波器滤波后的频谱');xlabel('频率/Hz');ylabel('幅值');4. 结果分析 图4-1由图4-1可以看出在f=1200Hz时满足阻带衰减要求，即As=100dB,且通带截止频率超过所给指标，即在f=1000Hz时，Ap<dB。滤波器的相位也满足线性特性。 图4-2 图4-3由上两图可以看出，信号通过FIR低通滤波器后，信号基本不发生变化，原始信号在低通滤波器通