(巴特沃兹滤波器.docVIP

巴特沃兹滤波器 (Butterworth) ??? 特点：具有通带内最大平坦的振幅特性，且随f↗单调↘ 其幅度平方函数具有如下形式： 式中，N为整数，称为滤波器的阶数，N越大，通带和阻带的近似性越好，过渡带也越陡。如下图所示： 图 巴特沃兹filter 振幅平方函数 过渡带：通带→阻带间过渡的频率范围， Ωc：截止频率。 理想滤波器的过渡带为O，阻带|H(jΩ)|=0，通带内幅度|H(jΩ)|=常数，H(jΩ)线性相? ??? 位。通带内，分母Ω/Ωc1,相应( Ω/Ωc)2N随N的增加而趋于0，A(Ω2)→1，在过渡带 ??? 和阻带，Ω/Ωc1，随N的增加，Ωe/Ωc1，所以A(Ω2)快速下降。 ?Ω=Ωc时， ，幅度衰减 ，相当于3bd衰减点。 振幅平方函数的极点可写成： Ha(-s).Ha(s)= 可分解为2N个一次因式令分母为零，→ 可见，Butterworth 滤波器的振幅平方函数有2N个极点，它们均匀对称地分布在|s|=Ωc的圆周上。 例：如图为N=3阶Butterworth 滤波器振幅平方函数的极点分布。 图 三阶A(-s2)的极点分布 ??? 考虑到系统的稳定性， Butterworth 滤波器的系统函数是由s平面左半部分的极点（SP3，SP4，SP5）组成的，它们分别为： 所以系统函数为： 式中 是为使S=0时Ha(s)=1而引入的。如用 归一化s，即s’=s/Ωc，得归一化的三阶BF： 如果要还原的话，则有 关于数字滤波器滤波器有很多种，讨论下对信号频率具有选择性的滤波器。这又分为模拟滤波器和数字滤波器。模拟滤波器是在传统模拟电路中发展起来的，其实就是RC电路网络。随着数字技术的发展，数字滤波器则越来越受到青睐。 ?数字滤波器分为递归型和非递归型，所谓递归即滤波器内部存在反馈回路，这种滤波器对单位冲击响应可以延续到无限长的时间，所以也叫 IIR (infinite impulse response filter) ；相应的，非递归型即内部不存在反馈，也叫 FIR(finite impulse response filter)，其传递函数不存在除零点意外的极点。数字滤波器的一般形式为： ?a(0)*y(n) = b(0)*x(n) + b(1)*x(n-1) + ... + b(nb)*x(n-nb)?- a(1)*y(n-1) - ... - a(na)*y(n-na) ?相应于上面的讨论，则a都为零则为IIR，a 有非零的则为FIR。显然，a(0)=1 方便讨论和设计滤波器，所以在matlab中滤波器设计都是 a(0)=1。容易看出，无内部反馈的FIR总是稳定的，具有IIR所没有的特点，但也可以证实，客观模拟电路中，FIR是无法实现的，只有通过数字处理技术设计，而且，要获得性能符合要求的FIR，滤波器的阶数必须设计得非常高，比如，一个常用的矩形窗，一般性能下就要求有100阶左右来拟合，计算代价太大。一般情况下，n阶IIR与2n阶FIR性能相同。实际情况下，20阶的Butterworth IIR 滤波器可以实现近似理想相位线性。 ??频率滤波器大概分为带通、带阻、高通、低通。特性不同的模拟滤波器中经典滤波器有Butterworth 和 Chebyshev 。其中，Butterworth 滤波器特点是通带处幅值特性平坦，而 Chebyshev 滤波器则比前者的截至特性要好，但通带处的幅值有振荡。前面提到，对于数字滤波器而言，可以采用不同阶数逼近相应滤波器，滤波器性能还与滤波器的阶数有关，一般而言，阶数越高，则逼近越精确，但计算代价也随之上升，所以性能与代价总需要寻求一个平衡点。对性能要求一定的情况下，如果对频率截至特性没有特殊要求，考虑采用Butterworth IIR滤波器。因为 Chebeshev 滤波器的波纹可能大多数情况下不能忍受。 ?接着我们看看怎么借助 Matlab 设计符合我们需求的Butterworth 数字滤波器。当然，我们尽可能了解 Butterworth Filter 的原理以及可能的话，再了解数字滤波器设计的方法理论，但是，我们不必自己动脑筋根据需求去设计每个系数，Matlab 内建有设计Filter 的函数。这里仅仅讨论 Butter，其语法格式为： ?[B,A] = butter(N,Wn,S)其中N?为要求设计的滤波器阶数，该参数的设定参照前面的讨论，如果没有实时性要求的话，可以定为 20 ，实话，慢～相当的慢～；?S ?为字符串，表明设计的滤波器类型，low低通/high高通/stop带阻?Wn 为要求的标准化截至频率，单位为rad/sample，如果是带阻滤波器，则Wn为长度为2的向量[w1 w2]。关于