带通滤波器的模拟.docx

《带通滤波器的模拟.docx》由会员分享，可在线阅读，更多相关《带通滤波器的模拟.docx（15页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 带通滤波器的模拟一、摘要本文介绍低通滤波器的设计方法，设计一个简单的带通滤波器，同时用Matlab 仿真,Audition 模拟。还观察比较了贝塞尔滤波器，巴特沃斯滤波器、切比雪夫滤波器，椭圆滤波器的特点。Abstract This article mainly introduced the designing of a filter, to design a bandpass filter.Meanwhile , simulating it by using Matlab and Audition ,for observing when a signal passes 4different

2、 types of BPF. Then made a comparison about the effect of these filters of Bessel，Butterworth, Chebychev, and elliptic二、关键词带通滤波器设计 Matlab 仿真 Audition 模拟三、前言在通信（如信号传递）过程中，会被外界信号干扰，信号中会叠加一部分噪声等无用信号，有时候在产生某一种信号时，一定会产生噪声，当中也不可避免的会混入一部分干扰信号，不利于信号的接收，这就给现实生活带来很大的麻烦，因此，在信号传递过程中，消除或者减弱干扰信号是必须的。这就需要提取信号，也就是滤

3、波。滤波器的种类很多，有低通，带通，高通，带阻，还有以人名命名的滤波器，本文将以带通滤波器为主题，进行探讨。 四、正文(一)实现带通滤波器的方法主要分为两种：【1】1、直接法，通过电阻、电容、电感和运算放大器等元器件直接设计的频率选择性带通滤波器（如高通滤波器和低通滤波器串联），但是这样的滤波器的中心频率由许多元件决定，如果要直接调节中心频率，全部元件都要按照某一种正确的方式同时变化。在实际应用中会很麻烦，而且仅制作一个固定特性的带通滤波器是比较麻烦的。2、间接法，利用信号的相乘性质,用恰当地将信号移动的方法改变滤波器的中心频率,这样就实现了中心频率可变的带通滤波器,中心频率可以通过调谐旋钮来

4、实现.本文将采用该方法实现带通滤波器的模拟.(二)原理 现在考虑一个系统,如图所示, 输入一个信号 x(t)，然后乘以一个复指数信号图 1然后信号通过一个截止频率为而傅里叶变换的性质可以得到。的低通滤波器，出来的信号再乘以()( ) ( ) ( ) ( )w d w w * ww wY j = X j = X j -cc图 1 中的y(t)，w(t)，f(t)的频谱图像如图 2 所示。这样就等于向右移 ， 在中靠近附近的频谱就经过频移移到了低通滤波器的通带内。同样，可以得到，F 向左移动 就可以得到 。由图2 可以看到，图 1 所示的F整个系统可以等效为一个中心频率是，图 3通带为的带通滤波器

5、，如图3 所示，这样只要改变复指数信号的频率 ，就可以轻松的改变带通滤波器的中心频率。（三）设计带通滤波器 【2】模拟滤波器技术指标的转换, f模拟带通滤波器的设计指标有中心频率 f ，通带截止频率 f，阻带截止频0p1p2, f率 f，通带最大衰耗 ，阻带最小衰耗 。s1s2（1）把模拟频率转换成相应的数字频率 通过公式（式中，T 为采样周期）可以把给定的模拟频率转换成数字频率。（2）通过频率预畸变把数字频率转换成相应的模拟角频率：2 W 2 W = tan= tanwwwp1wp2 2 2 p1s10Tp2T2W2 W = tanT= tanw1s2s 2 2 s2T2Ww wB=tan-

6、=02p2p1T（3）将指标转换为归一化模拟低通指标：经查表我们可以得到从归一化模拟低通滤波器到实际模拟带通滤波器的变换s2+ w2s =B0公式s BBLw w表示通带带宽。ww w 表示通带的中心频率，B =-式中：=0p1 p2p2p1w w-22w=将 s 代以 w ，则有0jBwBLB因此频率作如下变换w-ww -w22222ww=p10p20wwBc1Bc2p1p2w ww w-222ww=s10s20wwBr1Br2s1s2而模拟低通滤波器通带最大衰耗与阻带最大衰耗和模拟带通滤波器的相同。归一化模拟低通滤波器的设计（1） 确定归一化截频 归一化的模拟低通滤波器取w=1，而阻带截频

7、w 取w ，w 中绝对值较crr1r2小的。（2） 确定归一化模拟低通滤波器阶数低通滤波器的阶数 N 满足下面的不等式10 a -1lgr10 a -10.1N p（1.3.1）lg(w w)rc（通常 N 取满足条件的最小的正整数）因此，可以通过理论计算来确定 N 的值。此外，利用 Matlab 中的命令语句也可以确定出 N 的数值。把归一化的模拟低通滤波器转换成模拟带通滤波器s +w22=将 s分别代入系统函数可得带通滤波器的频率响应函数。0sB滤波器特点简介贝塞尔(Bessel filt滤e r波)器是具有最大平坦的群延迟（线性相位响应）的线性过滤器。 贝赛尔滤波器常用在音频天桥系统中。

8、模拟贝赛尔滤波器描绘为几乎横跨整个通频带的恒定的群延迟，因而在通频带上保持了被过滤的信号波形.虽然贝塞尔滤波器在它的通频带内提供平坦的幅度和线性相位（即一致的群延时）响应，但它的选择性比同阶（或极数）的巴特沃斯滤波器或切比雪夫滤波器要差。巴特沃斯滤波器(Butterworth filter)是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。切比雪夫滤波器(Chebychev filter),在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波

9、器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动I 型切比雪夫滤波器在通带（或称“通频带”）上频率响应幅度等波纹波动的滤波器称为“I 型切比雪夫滤波器 ”；II 型切比雪夫滤波器在阻带（或称“阻频带”）上频率响应幅度等波纹波动的滤波器称为“II 型切比雪夫滤波器”。椭圆滤波器（Elliptic filter），是在通带和阻带等波纹的一种滤波器。椭圆滤波器相比其他类型的滤波器，在阶数相同的条件下有着最小的通带和阻带波动。它在通带和阻带的波动相同，这一点区别于在通带和阻带都平坦的巴特沃斯滤波器，以及通带平坦、阻带等波纹或是阻带平坦、通带等波纹的切比雪夫滤波器。四种带通滤波器的

10、 Matlab 仿真Matlab 仿真模拟图 频率响应增益名称1贝 最23型滤波器 4型滤波器5椭 通圆 带滤 阻波 带器 都有波纹5（当阶数与巴特沃斯滤波器相同时，贝塞尔滤波器的增益小于0，对信号起衰减效应）(五)四种带通滤波器对粉红噪音信号的过滤模拟-基于 Audition 12（所设Transition Bandwith较窄，接近于理想带通滤波器）3.4.阶数增大，波纹条数增多信号通过滤波器后得到的频谱和波形图如下 （从左至右依次为贝塞尔，巴特沃斯，切比雪夫I 型，切比雪夫 II 型滤波器）频谱波形可以看出，当过渡带较窄时，所滤出的信号具有更加明显的通带-阻带区分。当切比雪夫滤波器阶数增

11、大时，其纹波条数增多。但当滤波器阶数越高时，滤波效果越明显。如何根据实际选择带通滤波器若要最大限度的保持信号不失真，一般选择贝塞尔带通滤波器，但是贝塞尔滤波器的的一个弊端是，在阶数相同的情况下，其选择性不如巴特沃斯等带通滤波器。只满足相频特性而不关心幅频特性。贝塞尔滤波器又称最平时延或恒时延滤波器。其相移和频率成正比，即为一线性关系。但是由于它的幅频特性欠佳，而往往限制了它的应用。从幅频特性提出要求，而不考虑相频特性。巴特沃斯滤波器具有最大平坦幅度特性，切贝雪夫滤波器也是从幅频特性方面提出逼近要求的，与巴特沃斯逼近特性相比较，这种特性虽然在通带内有起伏，但对 同样的 n 值在进入阻带以后衰减更

12、陡峭，更接近理想情况。切贝雪夫滤波器与巴特沃斯滤波器进行比较，切贝雪夫滤波器的通带有波纹，过渡带轻陡直，因此，在不允许通带内有纹波的情况下，巴特沃斯型更可取；从相频响应来看，巴特沃斯型要优于切贝雪夫型，只满足相频特性而不关心幅频特性。巴特沃斯型相频响应更接近于直线。至于椭圆滤波器，则是对于通带和阻带的波纹都有限制条件，而且还要求过渡带较小。五、结论本文给出了带通滤波器的模拟原理，并且给出了常用的带通滤波器的Matlab 程序以及彼此之间的异同，巴特沃斯滤波器的衰减因数小于椭圆滤波器，但是其阻带波纹小，最终趋向于 0，切比雪夫滤波器只是阻带或者通带有波纹。并通过Audition 进行了模拟滤波测

13、试，观察滤波器滤波效果。但是，对于滤波器的设计方法并没有给出过程，并没有用物理方法实现各种带通滤波器六、参考文献【1】奥本海姆.信号与系统（第二版）.刘树棠译【2】梁虹.信号与线性系统分析基于 MATLAB 的方法与实现.高等教育出版社 七、附录Matlab 代码基于 Matlab 设计的带通滤波器fs=20;fp1=2;fp2=4;fs1=1.5;fs2=4.5;Ap=0.0877;As=16.9897;wpa=2*pi*fp1/fs;wpb=2*pi*fp2/fs;wsa=2*pi*fs1/fs;wsb=2*pi*fs2/fs;c=sin(wpa+wpb)/(sin(wpa)+sin(wp

14、b);omegap=abs(c-cos(wpb)/sin(wpb);omegasa=(c-cos(wsa)/sin(wsa);omegasb=(c-cos(wsb)/sin(wsb);omegas=min(abs(omegasa),abs(omegasb);ep=sqrt(10(Ap/10)-1);es=sqrt(10(As/10)-1);N=ceil(log(es/ep)/log(omegas/omegap);omega0=omegap/ep(1/N);w=(0+eps):pi/300:pi;Hw2=1./(1+(c-cos(w)./(omega0*sin(w).(2*N);plot(w/p

15、i,Hw2);grid;巴特沃斯带通滤波器模拟程序wp1=0.45*pi; wp2=0.65*pi;ws1=0.3*pi;Rp=1;ws2=0.75*pi;Rs=40;Wp=wp1/pi,wp2/pi; Ws=ws1/pi,ws2/pi;N,Wn=buttord(Wp,Ws,Rp,Rs);B,A=butter(N, Wn);h,omega=freqz(B,A,64); subplot(1,2,1);plot(omega/pi,abs(h);grid;xlabel(omega/pi);ylabel(|H(s)|);subplot(1,2,2);plot(omega/pi,20*log10(abs

16、(h);grid;xlabel(omega/pi);ylabel(增益.dB);切比雪夫带通滤波器I 型模拟程序wp1=0.45*pi; wp2=0.65*pi;ws1=0.3*pi;Rp=1;ws2=0.75*pi;Rs=40;Wp1=tan(wp1/2); Wp2=tan(wp2/2);Ws1=tan(ws1/2); Ws2=tan(ws2/2);BW=Wp2-Wp1;W0=Wp1*Wp2; W00=sqrt(W0);WP=1;WS=WP*(W02-Ws12)/(Ws1*BW);N,Wn=cheb1ord(WP,WS,Rp,Rs,s);B,A=cheby1(N,Rp,Wn,s);BT,AT

17、=lp2bp(B,A,W00,BW);num,den=bilinear(BT,AT,0.5);h,omega=freqz(num,den,64);subplot(1,2,1);plot(omega/pi,abs(h);xlabel(omega/pi);ylabel(|H(s)|);subplot(1,2,2);plot(omega/pi,20*log10(abs(h);xlabel(omega/pi);ylabel(增益.dB);II 型模拟程序wp1=0.45*pi; wp2=0.65*pi;ws1=0.3*pi;Rp=1;ws2=0.75*pi;Rs=40;Wp1=tan(wp1/2);

18、Wp2=tan(wp2/2);Ws1=tan(ws1/2); Ws2=tan(ws2/2);BW=Wp2-Wp1;W0=Wp1*Wp2; W00=sqrt(W0);WP=1;WS=WP*(W02-Ws12)/(Ws1*BW);N,Wn=cheb2ord(WP,WS,Rp,Rs,s);B,A=cheby2(N,Rs,Wn,s);BT,AT=lp2bp(B,A,W00,BW);num,den=bilinear(BT,AT,0.5);h,omega=freqz(num,den,64);subplot(1,2,1);plot(omega/pi,abs(h);xlabel(omega/pi);ylabe

19、l(|H(s)|); subplot(1,2,2);plot(omega/pi,20*log10(abs(h);axis(0 1 -100 0);xlabel(omega/pi);ylabel(增益.dB);椭圆带通滤波器模拟程序wp1=0.45*pi; wp2=0.65*pi;ws1=0.3*pi;Rp=1;ws2=0.75*pi;Rs=40;Wp1=tan(wp1/2); Wp2=tan(wp2/2);Ws1=tan(ws1/2); Ws2=tan(ws2/2);BW=Wp2-Wp1;W0=Wp1*Wp2; W00=sqrt(W0);WP=1;WS=WP*(W02-Ws12)/(Ws1*B

20、W);N,Wn=ellipord(WP,WS,Rp,Rs,s);B,A=ellip(N,Rp,Rs,Wn,s);BT,AT=lp2bp(B,A,W00,BW);num,den=bilinear(BT,AT,0.5);h,omega=freqz(num,den,64);subplot(1,2,1);plot(omega/pi,abs(h);grid;xlabel(omega/pi);ylabel(|H(s)|);subplot(1,2,2);plot(omega/pi,20*log10(abs(h);grid;xlabel(omega/pi);ylabel(增益.dB);贝塞尔带通滤波器模拟程序

21、wp1=0.45*pi; wp2=0.65*pi;ws1=0.3*pi;Rp=1;ws2=0.75*pi;Rs=40;Wp=wp1/pi,wp2/pi; Ws=ws1/pi,ws2/pi;N,Wn=buttord(Wp,Ws,Rp,Rs);B,A=besself(N, Wn);h,omega=freqz(B,A,64);subplot(1,2,1);plot(omega/pi,abs(h);grid;xlabel(omega/pi);ylabel(|H(s)|);subplot(1,2,2);plot(omega/pi,20*log10(abs(h);grid;xlabel(omega/pi);ylabel(增益.dB);