贝塞尔滤波器 编程题.docVIP

编程题 --贝塞尔滤波器 1.研究热点及原理 电子学和信号处理中，贝赛尔(Bessel)滤波器是具有最大平坦的群延迟（线性相位响应）的线性过滤器。贝赛尔滤波器常用在音频天桥系统中。模拟贝赛尔滤波器描绘为几乎横跨整个通频带的恒定的群延迟，因而在通频带上保持了被过滤的信号波形。滤波器的名字来自于Friedrich贝赛尔，一位德国数学家（1784–1846），他发展了滤波器的数学理论基础贝塞尔(Bessel)滤波器具有最平坦的幅度和相位相应。带通（通常为用户关注区域）的相位响应近乎呈线性。Bessel滤波器可用于减少所有IIR滤波器固有的非线性相位失真。 贝塞尔(Bessel)线性相位滤波器正是由于具有向其截止频率以下的所有频率提供等量延时的特性，才被用于音频设备中，在音频设备中，必须在不损害频带内多信号的相位关系前提下，消除带外噪声。另外，贝塞尔滤波器的阶跃响应很快，并且没有过冲或振铃，这使它在作为音频DAC输出端的平滑滤波器，或音频ADC输入端的抗混叠滤波器方面，是一种出色的选择。贝塞尔滤波器还可用于分析D类放大器的输出，以及消除其它应用中的开关噪声，来提高失真测量和示波器波形测量的精确度。 虽然贝塞尔滤波器在它的通频带内提供平坦的幅度和线性相位（即一致的群延时）响应，但它的选择性比同阶（或极数）的巴特沃斯(Butterworth)滤波器或切比雪夫(Chebyshev)滤波器要差。因此，为了达到特定的阻带衰减水平，需要设计更高阶的贝塞尔滤波器，从而它又需要仔细选择放大器和元件来达到最低的噪声和失真度。620kHz的贝塞尔低通滤波器为例，从软件方面来了解贝塞尔低通滤波器的特性及其实现方法。Matlab程序流程如图1所示 3.贝塞尔低通滤波器的Matlab仿真实现程序代码及主要注释 [b,a]=besself(5,620000*2*pi); t= 6e-6; impulse(b,a,t); w=[O:lO000*2*pi:1000000*2*pi]; h1=freqs(b,a,w1); mag=abs(h1); mag=20*log10(mag); f=w/(2*pi); subplot(2,1,1); plot(f,mag);grid; xlabel(Frequency(Hz)); ylabel(Magnitude(dB)); phase=angle(h1); phase=phase*180/pi; subplot(2,1,2);grid; plot(f,phase); xlabel(Frequency(Hz)); ylabel(Phase(degrees));[z,p,k]=Besself(5,620000*2*pi0); format long delay=O; for w=0:6.2e5/1O0:6.2e5*2*pi; for i=1:5 rp=real(p(i)); ip=imag(p(i)); tempp=rp/(rp^2+(ip-w)^2); delay=delay-tempp; end delay delay=0; end 程序中大量使用了matlab信号处理工具箱中的函数，现将几个主要的工具箱函数简要介绍如下： （1）Bsself 功能：贝塞尔模拟滤波器的设计，包括低通、高通、带通、带阻滤波器。 格式：[b,a]=besself（n,Wn） [b,a]=besself(n,Wn,’ftype’) [z,p,k]=besself(…) [A,B,C,D]= besself(…) 说明：besself设计低通、高通、带通、带阻贝塞尔模拟滤波器。模拟贝塞尔滤波器最主要特征是通带内固定的群时延，这样能够保持通带内波形的形状。数字贝塞尔滤波器没有这个特性，besself不支持数字贝塞尔滤波器的设计。 [b,a]=besself（n,Wn）可设计出一个n阶截止频率为wn弧度／秒的贝塞尔低通滤波器。依s幂次的增长，滤波器传递函数的系数分别放在n+l维向量a和b中。其传递函数为 贝塞尔滤波器的幅度响应在截止频率处小于3dh，且随着阶数n的增大，衰减也越大。 如果Wn是一个二维向量，wn：[wl w2]，且wlw2，则besself(n，wn)设计出的是一个2n阶带宽为wlww2的带通模拟滤波器。 [b，a]=besself(n，Wn，‘ftype’)可用于设计高通和带阻滤波器。其中‘ftype’可为： ‘high’,n阶高通模拟滤波器，截止频率为Wn； ‘stop’,2n阶带阻模拟滤波器如果Wn为二维向量wn=[wl w2]，则阻带为wlww2。 [z，P，k]=besself(…)返回的是零、极点、增益形式转移函数，括号中的参数意义同上。 [A，B，C，D]= besself(…)返回的是滤波器的状态方程形式，如下： 其