脉冲响应不变法设计巴特沃斯低通滤波器

1、脉冲响应不变法设计巴特沃斯低通滤波器姓 名学 号专 业文档第一章摘要滤波器是自动控制、信号处理和通信领域的重要组成部分，广泛地应用于各种系统中。MATLAB语言是一种简单、高效的高级语言,是一种容丰 富、功能强大的分析工具，其应用围几乎覆盖了所有的科学和工程计算领域。MATLAB中提供了丰富的用于模拟滤波器设计的函数,通过编程可以很容易实现低通、高通、带通、带阻滤波器，并能画出滤波器的幅频特性曲线，大大简化了模拟滤波器设计。本文将通过利用MATLAB滤波器设计函数直接实现巴特沃斯低通滤波器的设计，找到应用MATLAB来设计低通滤波器的方法。介绍了用脉冲响应不变法设计巴特沃斯低通滤波器的基本理

2、论和设计思想，给出了基于MATLAB计巴特沃斯低通滤波器的具体步骤， 文中还给出了利用 MATLAB产生的滤波器的幅频及相频特性曲线。详细介 绍了在基于MATLAB计低通滤波器过程中常用到的工具和命令。关键字：巴特沃斯低通滤波器脉冲响应不变法第二章引言滤波器设计在电子工程、应用数学和计算机科学领域都是非常重要 的，。而低通滤波器在信号处理中的作用等同于其它领域如金融领域中移 动平均数(moving average)所起的作用，低通滤波器有很多种，其中， 最通用的就是巴特沃斯滤波器和切比雪夫滤波器。巴特沃斯滤波器的特点是通频带的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。在振幅

3、的对数对角频率的波得图上，从某一边界角频率开始，振幅随着角频率的增加而逐步减少，趋向负 无穷大。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论 阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数 越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和 地结束的振幅对角频率有不同的形状。第三章基本原理3.1 滤波器的基本理论3.1.1 滤波器的分类滤波器分两大类：经典滤波器和现代滤波器。经典滤波器是假定输入信号x（n）中的有用成分和希望取出的成分各自占有不同的频带。这样，当x（n）通过一个线性系统 （即滤波器）后可将欲去除的成分有效的去除。现代滤波器理论研究的主

4、要容是从含有噪声的数据记录（又称时间序列）中估计出信号的某些特征或信号本身。经典滤波器分为低通、高通、带通、带阻滤波器。每一种又有模拟滤波器（AF）和数字滤波器（DF）。对数字滤波器，又有IIR滤波器和FIR滤波器。IIR DF的转移函数是:H(z)丫X(z)Mbrz rr 0N1 akZk 1FIR DF的转移函数是:H(z)N 1h(n)z nn 0FIR滤波器可以对给定的频率特性直接进行设计，而IIR滤波器目前最通用的方法是利用已经很成熟的模拟滤波器的设计方法进行设计。滤波器的技术要求低通滤波器： p :通带截止频率S :阻带下限截止频率(p, s的单位dB)p :通带上限角频率s ：阻

5、带下限角频率ps3.2模拟滤波器的设计巴特沃思(Butterworth) 滤波器的介绍：因为我们设计的滤波器的冲击响应一般都为实数，所以有* 2G(s)G (s) G(s)G(s)sj G(j s)2这样，如果我们能由 p， p， s， s求出G(j )，那么就容2易得到所需要的G(s)。不同类型的 G(j )的表达式，代表了几种不同类型的滤波器。巴特沃思 (Butterworth) 滤波器为：G(j )2c2(2)C为待定常数，N为待定的滤波器阶次。3.2.2 巴特沃思模拟低通滤波器的设计由于每一个滤波器的频率围将直接取决于设计者的应用目的，因此必然是千差万别。为了使设计规化，我们需要将滤波

6、器的频率参数作归一化处理。设所给的实际频率为(或f )，归一化后的频率为，对低通模 拟滤波器，令1 /显然，p ', s s/ p。又令归一化复数变量为 p, p j ,显然P j j / p s/ p所以巴特沃思模拟低通滤波器的设计可按以下三个步骤来进行。(1)将实际频率 归一化C,N可按下式求解:C210 P/101lg10 s/10 1/1010 p1lg s若令p = 3dB,则C= 1,这样巴特沃思滤波器的设计就只剩一个参数N,这时2 1G(j )1(/ p)2N(3)确定 G(s)因为p j，根据上面公式有G(p)G(P)2N1(P/J)1N 2N1(1) pN 2N由 1

7、( 1) P0解得k N 1pk exp(j 2N 2 ),k = 1 , 2,，2N这样可得Gk(P)(P Pk)( P Pn 1 k )2 p cos(2k N2N求得G(P)后，用s/P代替变量P，即得实际需要得G(s)。3.3数字滤波器的设计按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。根据转换后的技术指标设计模拟低通滤波器G(s)；再按一定的规则将G(s)转换成H (z)。若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标， 然后按上述步骤2)设计出低通G(s)，再将G(s)转换为所需的H (z)。第四章设计过程4.1 设计低通数

8、字滤波器4.1.1 数字低通滤波器设计结构图用脉冲响应不变法转换原理如下图：Ha (s)拉氏逆变换 ha (t)等间隔采样ha(nT)-h(n) Z变换H(z) 图4.1结构框图4.1.2设计步骤设计要求：脉冲响应不变法设计数字低通滤波器，给定技术指标为：fp 100Hz， fs 300Hz， p 3dB， s 20dB，采样频率 Fs 1000 Hzo首先应该得到角频率，然后再按如下步骤求解：(1 )将数字滤波器的技术要求转换为模拟滤波器的技术要求。(2 )设计低通滤波器 G(S)由 =/ p依次求出 p,s，再求出N,可得G(p)G(s) G(p) p 亠然后由p转换成G(s)(3)由 G

9、(s)求 H(z)4.2 滤波器的MATLAB实现4.1.1 在MATLABF，设计巴特沃斯IIR滤波器可使用butter函数。IIR滤波Butter函数可设计低通、高通、带通和带阻的数字和模拟 器，其特性为使通带的幅度响应最大限度地平坦，但同时损失截止频率处的下降斜度。在期望通带平滑的情况下，可使用butter函数。butter函数的用法为：b,a=butter(n,Wn,/ftype/)其中n代表滤波器阶数，Wn代表滤波器的截止频率，这两个参数可使用buttord函数来确定。buttord函数可在给定滤波器性能的情况下，求出巴特沃斯滤波器的最小阶数n,同时给出对应的截止频率 Wn butt

10、ord 函数的用法为：n,Wn= buttord(Wp,Ws,Rp,Rs)其中Wp和 Ws分别是通带和阻带的拐角频率(截止 频率)，其取值围为0至1之间。当其值为1时代表采样频率的一半。Rp和Rs分别是通带和阻带区的波纹系数。不同类型(高通、低通、带通和带阻)滤波器对应的Wp和Ws值遵循以下规则：a .高通滤波器：Wp和 Ws为一元矢量且Wp>Wsb .低通滤波器：Wp和 Ws为一元矢量且Wp<Wsc .带通滤波器：Wp和Ws为二元矢量且Wp<Ws 如Wp=0.2,0.7,Ws=0.1,0.8;d .带阻滤波器：Wp和 Ws为二元矢量且 Wp>Ws 如Wp=0.1,0.

11、8,Ws=0.2,0.7。4.1.2 Matlab的信号处理工具箱提供了频谱分析函数：freqz 函数。freqz函数功能：离散时间系统的频率响应。格式：h,w=freqz(b,a, n)h,f=freqz(b,a ,n, Fs)h=freqz(b,a,w) h=freqz(b,a,f,Fs) freqz(b,a ,n)说明：freqz用于计算数字滤波器H(Z)的频率响应函数 H(ej w )。h,w=freqz(b,a,n)可得到数字滤波器的n点复频响应值，这 n个点均匀地分布在0, n 上，并将这n个频点的频率记录在w中，相应的频响值记录在h中。要求n为大于零的整数，最好为2的整数次幕，以

12、便采用 FFT计算，提高速度。缺省时 n =512。h,f=freqz(b,a,n,Fs)用于对 H(ej w)在0,Fs/2上等间隔采样n点，采样点频率及相应频响值分别记录在f和h中。由用户指定 FS(以 HZ为单位)值。h=freqz(b,a,w) 用于对H(ej w )在0,2 n 上进行采样，采样频率点 由矢量w指定。h=freqz(b,a,f,Fs) 用于对H(ej w )在0,FS上采样，采样频率点由 矢量f指定。freqz(b,a, n)用于在当前图形窗口中绘制幅频和相频特性曲线。4.1.3 采用脉冲响应不变法将模拟滤波器转换为数字滤波器可使用impinvar 函数。B Z,AZ

13、 = IMPINVAR(B,A,Fs)生成一个数字滤波器。BZ和AZ分别为脉冲值，结果与模拟滤波器的值相当。B和A为采样区间，采样周期为1/Fs。当Fs缺省时，默认采样频率为1Hz。第五章程序和仿真图5.1 设计程序fp=100;fs=300;Rp=3；Rs=20;Fs=1000;Wp=2*pi*fp;Ws=2*pi*fs;Nn=256;n=(0:100-1);% 采样点数N,Wn=buttord(Wp,Ws,Rp,Rs,'s');b,a=butter(N,Wn,'s');w=li nspace(1,400,100)*2*pi;H=freqs(b,a,w); %

14、在0,2 n 上进行采样，采样频率点由矢量W旨定figure(1);plot(w/(2*pi),20*log10(abs(H);title('巴特沃斯模拟滤波器幅频特性');xlabel('频率 /Hz');ylabel('幅度 /db');bz,az=imp in var(b,a,Fs);figure(2);zpla ne(bz,az);title('脉冲响应不变法零极点分布图');figure(3);Subplot(2,1,1);freqz(bz,az,Nn,Fs);对 H(ej w )在0,Fs/2上等间隔采样 Nn点gri

15、d on;title('巴特沃斯数字滤波器幅频特性');xlabel(' 幅度 /Hz');ylabel('频率 /db');Subplot(2,1,2);title('巴特沃斯数字滤波器相频特性');xlabel(' 幅度 /Hz');ylabel('频率 /db');grid on;figure(4);x=cos(2*pi*50* n/Fs)+cos(2*pi*500* n/Fs); plot( n,x);title('滤波前输入信号图像');xlabel('时间(s)&

16、#39;);ylabel('幅值');grid on;figure(5);y=filter(bz,az,x);%对原信号进行滤波plot (n /fs,y);title('滤波前输入信号图像');xlabel('时间(s)');ylabel('幅值');5.2 特性曲线仿真图File EditView lrsert Foos Deplete p Windew kelpbTh|o®cz*iainH|召特注斷模拟淤浪龍帳频特性-5*10S0 1001502002503003£0祖 0频率Nr-2U1o3sfEd F

17、.EERUJ-View In&ert T&cls Desktop Wind&w Helpi忆|j飞料跟朋 a t=i _ a脉;巾响应K亞法的霉樋点因Rfil Part图5.2脉冲响应不变法的零极点图-10q碰-70-3040100200-300File Edit View Insert Took De-sklup Wirduw Help为巴日 | 4 |鼻削辎归晟 |凰| 口討 口巴挣沃斯喲卞浇泯器悟频特性和 1001£0200250300350 4DD 4£0500巴昏沃斯数M連滾器相频特性50100150200250300350SCO图5.4滤波前输入信号图5.4滤波后输出信号图像第六章结语低通滤波器的应用十分广泛，利用Matlab语言，可以很容易地设计出巴特沃斯数字低通滤波器，与其他高级语言的程序设计相比,MA TLAB环境下更方便、快捷，节省大量的编程时间，提高编程效率，且参数的修改也十 分方便还可以进一步进行优化设计。巴特沃斯低通滤波器的振幅对角频率单调下降，无论阶数，振幅对角频率曲线都保持同样的形状的滤波器可以用较少的阶数获得很高的选择特性，在衰减特性要求不敏感的场合，适合用巴特沃斯型滤波器。心得体会通过一周的设计过程，让我感触最深的就是 M