基于matlab的低通滤波器.pdf

《基于matlab的低通滤波器.pdf》由会员分享，可在线阅读，更多相关《基于matlab的低通滤波器.pdf（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 基于 matlab 的低通滤波器 The Standardization Office was revised on the afternoon of December 13,2020 基于 matlab的低通滤波器 摘要:调用 MATLAB 信号处理工具箱中滤波通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。应用最广泛的是双线性变换法。基本设计过程是：先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标；设计过渡模拟滤波器；将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。MATLAB 信号处理工具箱中的各种 IIR 数字滤波器设计函数都是采用双线性变换法。关键词:滤波

2、器，matlab，c语言，声音 Abstract:call MATLAB signal processing toolbox filtering through the observation filter of input and output signals time domain waveform and spectrum,establish the concept of digital filter.One of the most widely applied is double linear transformation method.The basic design process

3、 is:(1)to a given digital filters index converted into analog filter transition index;(2)the design transition simulation filter;(3)transition simulation filter system function will be converted into digital filter system function.MATLAB signal processing toolbox digital filter function design of II

4、R is bilinear transformation method.Keywords:filter,matlab,the c language,the voice 一.任务：用 matlab 软件设计 IIR 模拟、数字以及各种窗函数的 FIR 低通滤波器 二.设计目的：（1）了解 matlab 软件的用途以及用法；（2）了解用冲激响应不变法设计模拟低通滤波器；（3）了解用脉冲响应不变法设计的巴特沃思数字低通滤波器；（4）了解基于汉明窗函数的 FIR 低通滤波器的设计。三.设计内容：用冲激响应不变法设计模拟低通滤波器 设计内容：要求按照设计指标设计无限冲激响应IIR 巴特沃什模拟低通滤波器

5、。设计原理：低通滤波器的技术要求用图形表示如下：|)(|jeHpsc10.707112 通带 过渡带 阻带 如上图表示了一个频域设计、一维低通滤波器的技术要求图。ws 和 wp 分别称为通带截止频率和阻带截止频率。通带频率范围为 0wwp，阻带频率范围 wsw。从wp到 ws 称为过渡带，用 w表示，在过渡带里，幅频特性单调下降。在通带和阻带内的衰减一般用数 dB 表示，通带内允许最大衰减是 s，阻带内允许最小衰减是 p。定义分别为：Matlab 信号处理工具箱中提供了设计巴特沃思模拟滤波器的函数 buttord、buttap 和 butter，格式如下：,(,CPSPSN Wbuttord

6、W WRRs 用于计算巴特沃思模拟低通滤波器的阶 N和 3dB截止频率 Wc(即本书中的符号c)。其中，Wp 和 Ws 分别是滤波器的通带截止频率p和阻止截止频率s，单位为rad/s；Rp和 Rs 分别是通带最大衰减系数p和阻带最小衰减系数s，单位为 dB。,()z p Gbuttap N用于计算 N阶巴特沃思归一化（c=1）模拟低通滤波器系统函数的零、极点和增益因子，返回长度为 N的向量 z 和 p 分别给出 N个零点和极点，G是滤波器增益。得到的滤波器系统函数形式如下：1212()()()()()()()()()aNaaNQsszszszHsGP sspspsp 其中，kz和kp分别是向量

7、 z 和 p 的第 k 个元素。如果要从零、极点得到系统函数的分子和分母多项式系数向量 B和 A，可以调用结构转换函数：,2(,)B Azp tf z p G 结构转换后系统函数的形式为：111111()()()MMMaNNNbsbsbB sHsA sa sasa 其中，M 是向量 B的长度，N是向量 A的长度，kkba和分别是向量 B和 A的第 k个元素。,(,)CB Abutter N Wftype s 用于计算巴特沃什模拟滤波器系统函数中分子和分母多项式系数向量 B 和 A，其中 N和CW分别是滤波器的阶和 3dB截止频率c，返回向量 B和 A中的元素ka和kb分别是上面的()aHs表示

8、式中的分母和分子系数。ftype缺省时，设计低通滤波器，S 缺省时，设计数字滤波器。pj0pj(e)20lg(e)HHsj0sj(e)20lg(e)HH 设计程序如下：（1）巴特沃思模拟低通滤波器设计指标如下：wp=*pi,ws=*pi,Rp=5,Rs=10;实验所用到的程序如下 wp=*pi;ws=*pi;Rp=5;Rs=10;N,wc=buttord(wp,ws,Rp,Rs,s);B,A=butter(N,wc,s);w=0:6;H,W=freqs(B,A,w);H=20*log10(abs(H);plot(w,H),grid on;xlabel(频率/(rad/s);ylabel(幅度/

9、dB)经 matlab 仿真以后得到的滤波器如下图所示：图 巴特沃思模拟低通滤波器 运行后在 matlab控制窗口上显示如下结果：N=7 （2）为了方便观察故设计了不同阶次 N的滤波器，即 Butterworth模拟原型低通滤波器，其程序代码如下：n=0:2;for ii=1:4 switch ii case 1,N=2;case 2,N=5;case 3,N=10;case 4,N=20;endz,p,k=buttap(N);0123456-60-50-40-30-20-10010频 率/(rad/s)幅度/dB巴 特 沃 什 模 拟 低 通 滤 波 器b,a=zp2tf(z,p,k);H,

10、w=freqs(b,a,n);magH2=(abs(H).2;hold on;plot(w,magH2);end xlabel(w/wc);ylabel(|H(jw)|2);title(Butterworth模拟原型低通滤波器);text,n=2)text,n=5)text,n=10)text,n=20)grid on;仿真后结果图形如图：图 Butterworth模拟原型低通滤波器（3）利用冲激响应不变法，把系统函数为651)(2ssssHa的模拟滤波器变换成等价的数字滤波器，采样间隔 T=1S。先经过计算得到：2132651)(2ssssssHa，则有如下式子 21112130067.01

11、851.012209.00000.1112)(zzzzeTzeTzHTT Matlab 实验参考程序如下:B=1,1;00.20.40.60.811.21.41.61.8200.20.40.60.811.21.4w/wc|H(jw)|2Butterworth 模 拟 原 型 滤 波 器n=2n=5n=10n=20A=1,5,6;T=1;Fs=1/T;Bz,Az=impinvar(B,A,Fs);%用冲激响应不变法将模拟滤波器变换成数字滤波器 运行结果：Bz=Az=得到的数据和先前计算得到的数据相同，故利用冲激响应不变化法把模拟滤波器设计成数字滤波器正确。实验结果及分析（1）步骤 1 中产生的低

12、通模拟滤波器的系统函数，经计算为：经过仿真后的结果，与上式对比，所得结果进分析相同。此处为 N=7 的低通模拟滤波器（2）在 Butterworth模拟原型低通滤波器中，可以看到，滤波器的幅频特性随着滤波器阶次 N的增加而变得越来越好。（3）步骤 3 中，利用冲激响应不变法，把系统函数的模拟滤波器变换成等价的数字滤波器，在 matlab仿真软件中验证了冲激响应不变法的正确性。用脉冲响应不变法设计的巴特沃思数字低通滤波器 设计内容：设计巴特沃思数字低通滤波器和椭圆数字低通滤波器，要求通带边界频率 fp=，通带最大衰减 Rp=；阻带边界频率 fs=8kHZ，阻带最小衰减 Rs=30dB,采样频率为

13、Fs=20kHZ。设计步骤：用 MATLAB进行典型的数字滤波器的设计，一般步骤如下:（1）将设计指标归一化处理。如果采用双线性变换法，还需进行预畸变。（2）根据归一化频率，确定最小阶数 N和频率参数 Wn。可供选用的阶数择函数有:buttord，cheblord，cheb2ord，ellipord 等。（3）运用最小阶数 N设计模拟低通滤波器原型。模拟低通滤波器的创建函数有:buttap,cheblap,cheb2ap,ellipap 和 besselap，这些函数输出的是零极点式形式，还要用zp2tf函数转换成分子分母多项式形式。如果想根据最小阶数直接设计模拟低通滤波器原型，可用 butt

14、er,chebyl,cheby2,ellip,bessel等函数，只是注意要将函数中的 Wn 设为1。（4）根据第 2 步的频率参数 Wn，模拟低通滤波原型转换模拟低通、高通、带通、带阻滤波器，可用函数分别是：lp21p,lp2hp,lp2bp,lp2bs。（5）运用脉冲响应不变法或双线性变法把模滤波器转数字滤波器，调用的函数是impinvar和 bilinear。脉冲响应不变法适用于采样频率大于 4倍截止频率的锐截止低通带通滤波器，而双线性变换法适合于相位特性要求不高的各型滤波器。（6）根据输出的分子分母系数，用 tf函数生成 H（z）的表达式，再用 freqz 函数验证设计结果。设计程序如

15、下：fp=2100;fs=8000;）（）（1893.1s4144.152.024144.1483.12s1.4144)2.143Ss23625.3)(ssssHaFs=20000;Rp=;Rs=30;T=1/Fs;%设计指标 W1p=fp/Fs*2;W1s=fs/Fs*2;%求归一化频率 N,Wn=buttord(W1p,W1s,Rp,Rs,s);%确定 butterworth的最小介数 N 和频率参数Wn z,p,k=buttap(N);%设计模拟低通原型的零极点增益参数 bp,ap=zp2tf(z,p,k);%将零极点增益转换成分子分母参数 bs,as=lp2lp(bp,ap,Wn*pi

16、*Fs);%将低通原型转换为模拟低通 bz,az=impinvar(bs,as,Fs);%用脉冲响应不变法进行模数变换 sys=tf(bz,az,T);%给出传输函数 H(Z)H,W=freqz(bz,az,512,Fs);%生成频率响应参数 subplot(2,1,1);plot(W,20*log10(abs(H);%绘制幅频响应 grid on;%加坐标网格 xlabel(频率/Hz);ylabel(振幅/dB);subplot(2,1,2);plot(W,abs(H);grid on;xlabel(频率/Hz);ylabel(振幅/H);运行后的波形如下：图 巴特沃思和椭圆数字低通滤波器

17、 N=4 bz=az=基于汉明窗函数的 FIR 低通滤波器的设计 设计具体低通滤波器要求：设低通滤波器通带截止频率为0.2p，阻带截止频率为0.4s，通带纹波pA=，最小阻带衰减rA=50dB。设计内容及原理：对于线性相位FIR滤波器通常采用窗函数法设计。窗函数设计FIR滤波器基本思想，就是从时域出发,把理想的无限长的)(nhd用一定形状的窗函数截取成有限长的()h n，以此()h n来逼近()dhn，从而使所得到的频率响应()jwH e与所要求的理想频率响应()jwdHe相接近。理想的数字滤波器频率特性()jwdHe是无法实现的，FIR的设计就是要寻找一个可以得到的频率特性()jwH e来逼

18、近()jwdHe，这相当于用一个可实现的单位脉冲响应()h n去逼近一个理想单位脉冲响()dhn。要想用一个有限长的因果序列去逼近它，最简单的方法是截取n从0N-1的一段来表示它，即：()()dh nh n（0nN-1）；当n为其他的数时，()h n=0。同时，为了保证线性相位，还要满足偶对称：()(1)h nh Nn ，这就好像通过一个窗口观看到的一段()dhn，因此()h n就表示成()dhn和一个“窗口函数”的乘积。这样，对()h n的求解就变为：()()()dh nh nW n这里的()W n就称为窗口函数。窗口函数相当于对()dhn在时域上截取一段，势必造成频域的矩形窗口的失真，结果

19、就是截取出的信号也相应失真。为了补偿这种失真，只有改变原来窗口的形状，修正经过时域截取后的窗口失真。改善窗口函数的标准：（1）尽量减少窗口频谱的旁瓣，使能量集中在主瓣中，以减少波动的幅度和次数，提高阻带衰减。这样，主瓣就不能太窄；（2）主瓣的宽度应尽量窄，以得到尽量窄的过渡带。工程中常用的窗函数有以下六种：矩形窗，三角形窗，汉宁窗，汉明窗，布莱克曼窗，凯泽窗(可变窗函数)。表1为六种窗函数的性能指标。窗函数 窗普性能指标 加窗后滤波器性能指标 -旁瓣峰值 dB 主瓣宽度 2/n 过滤宽度 2/n 阻带最小衰减 矩形窗 -13 2 -21 三角窗 -25 4 -25 汉宁窗 -31 4 -44

20、汉明窗 -41 4 -53 布莱克曼窗 -57 6 -74 凯泽窗 -57 5 5 -80 表 常用传函数的性能指标 根据条件采用汉明窗设计程序如下：wp=*pi;wr=*pi;trwidth=wr-wp;%过渡带宽度 N=ceil*pi/trwidth)+1;%滤波器的长度 wc=(wr+wp)/2;hd=ideal_lp(wc,N);w_ham=(hamming(N);h=hd.*w_ham;db,mag,pha,grd,w=freqz_m(h,1);delta_w=2*pi/1000;Ap=-(min(db(1:1:wp/delta_w+1);Ar=-round(max(db(wr/de

21、lta_w+1:1:501);subplot(221);stem(n,hd);title(理想单位脉冲响应 hd(n)subplot(222);plot(n,w_ham);title(海明窗 w(n)subplot(223);stem(n,h);title(实际单位脉冲响应 h(n)subplot(224);plot(w/pi,db);title(幅度响应(dB)输出结果如下：图 由上图可以看到：通带纹波为，最小阻带衰减为 55dB。满足设计指标的要求。滤波器系数 N=35 具体系数如下：下面我把滤波器的阶数人为的设置成 15，我们得到如下的输出结果：图 由上图可以看到：通带纹波为，最小阻带衰

22、减为 16dB。显然不能满足指标要求。滤波器系数 N=15，具体系数如下：综上所述：当滤波器的阶数增加时，可以减小通带纹波，增加最小阻带衰减。但同时计算量也随之增大，成本也要同样增加。下面，我用 MATLAB 中的 Filter Design and Analysis 工具进行同样的设计，设置阶数为 35，低通，窗口选择 hamming，归一化频率0.3c。结果如下：图 幅频特性曲线 图相频特性曲线 图冲击响应 图零极点图 实验结论及分析 FIR滤波器以它优越的性能，在数字信号处理领域中占有很重要的地位，是现代电子技术中必须掌握的设计技术。本文通过讨论FIR数字滤波器的基本原理和设计步骤，阐述

23、了窗函数法设计FIR数字滤波器的思路。在MATLAB中,通过窗函数法设计了FIR数字滤波器，并很好的完成了设计要求的指标。参考文献 1张辉，曹丽娜，王勇，通信原理辅导.，西安：西安电子科技大学出版社，1992.2王兴亮，等.数字通信原理与技术.西安：西安电子科技大学出版社，2001.3宋祖顺，等.现代通信原理.北京：电子工业出版社，2001.4赵树杰等，数字信号处理，西安电子科技大学出版社，.5 陈怀琛等，MATLAB 及在电子信息课程中的应用，电子工业出版社出版，2002.6李丽 王振领，MATLAB 工程计算及应用，人民邮电出版社，2001.7 陈怀琛等译，数字信号处理及其 MATLAB 实现，电子工业出版社；2002.