matlab巴特沃斯滤波器设计.docVIP

数字IIR带阻滤波器的设计 （基于巴特沃斯法） 1、 数字带阻IIR滤波器设计 IIR数字滤波器在很多领域中有着广阔的应用。与FIR数字滤波器相比,它可以用较低的阶数获得高选择性,所用存储单元少,经济而效率高,在相同门级规模和相同时钟速度下可以提供更好的带外衰减特性。下面介绍一种设计实现IIR数字滤波器的方法。 设计步骤如下： 步骤1：将数字滤波器H(z)的技术指标ωp和ωs，通过Ω＝tan(ω/ 2)转变为模拟滤波器G(s)的技术指标Ωp和Ωs，作归一化处理后，得到ηp＝1，ηs＝Ωs/Ωp； 步骤2：化解为模拟原型滤波器G(s)的技术指标； 步骤3：设计模拟原型滤波器G(p)； 步骤4：将G(p)转换为模拟滤波器的转移函数G(s)； 步骤5：将G(s)转换成数字滤波器的转移函数H(z)，s=(z–1)(z+1)。 所谓原型滤波器是指归一化的低通滤波器。本节主要讨论通过IIR数字滤波器的原型转换设计法和IIR数字滤波器的直接设计方法来设计数字高通、带通及带阻滤波器其转换方法主要有3种：一是直接由模拟低通滤波器转换成数字高通、带通或带阻滤波器；二是先由模拟低通滤波器转换成模拟高通、带通或带阻滤波器，然后再把它转换成 1 相应的数字滤波器；三是将模拟低通滤波器先转换成数字低通滤波器，再通过变量代换变换成高通、带通或带阻滤波器。 数字IIR滤波器的设计结构如图： 图1、IIR滤波器的设计步骤 2、带阻滤波器的设计规范： 本设计中分别用ha(t)、Ha(s)、Ha(jΩ)表示模拟滤波器的单位脉冲相应、系统函数、频率响应函数，三者的关系如下： ∞ Ha(s) = LT[ha(t)] =?ha(t)e-stdt -∞ ∞ Ha(jΩ)=FT[ha(t)] =?ha(t)e-jΩtdt -∞ 可以用ha(t)、Ha(s)、Ha(jΩ)中任一个描述模拟滤波器，也可以 用线性常系数微分方程描述模拟滤波器。但是设计模拟滤波器时，设计指标一般由幅频相应函数|Ha(jΩ)|给出，而模拟滤波器设计就是根据设计指标，求系统函数Ha(s)。 工程实际中通常用所谓的损耗函数即衰减函数A（Ω）来描述滤波器的幅频响应特性，对归一化幅频响应函数，A（Ω）定义如下： A（Ω）=-20lg|Ha(jΩ)|=-10lg|Ha(jΩ)|2 dB 模拟带阻滤波器的设计指标参数有Ωp1、Ωs1、Ω0、Ωsu、Ωpu。Ωp1和Ωpu分别表示带阻滤波器的通带下边界频率和通带上边界频率；Ωs1和Ωsu分别表示带阻滤波器的阻带下边界频率和阻带上边界频率。δ1和δ2分别称为通带和阻带波纹幅度 图2、带阻滤波器 带阻滤波器的技术指标要求必须是几何对称的。如果带阻滤波器不是几何对称的，为了使对称的带阻滤波器满足最小阻带衰减要求，在由非对称的带阻滤波器变换为对称的带阻滤波器时，保留非对称带阻滤波器的阻带截止频率调整两个通带截止频率中的一个，调整过程 3 步骤如下： 1、 2、 3、 4、 计算Ω02=Ωs12Ωsu Ω0计算 ，如果 Ωp1 =Ωp1 ΩΩpu2p1,用 Ωp1 代替Ωp1； 如果 Ωp1 ≤Ωp1，计算ΩpuΩ0=，并用Ωpu代替Ωpu； Ωp1如果Ap1≠ Ap2，选择AP=min{ Ap1, Ap2}。 所有类型的变换，包括各自与低通原型的频率变换如表1所示。用频率变换法设计一个模拟滤波器的步骤如下： 1、确定低通、高通、带通和带阻模拟滤波器的技术要求。当要求的带通或带阻滤波器是非几何对称时，根据上面描述的过程使这些技术要求成几何对称的； 2、根据表1中的变换关系，确定归一化低通滤波器的技术要求：通带截止频率为λp，阻带截止频率为λs，阻带衰减为 AP(dB),阻带衰减为As (dB)； 3、根据λp、λs、AP 和As，用巴特沃斯设计归一化低通滤波器； 4、根据表1最右栏的变换关系，得到要求的非归一化模拟滤波器。 表1 模拟滤波器的频率变换 下面我们用巴特沃斯滤波法来设计一个数字IIR带阻型滤波器： 滤波器规范：通带下边界频率： 50HZ 通带上边界频率： 450HZ 阻带下边界频率： 200HZ 阻带上边界频率： 300HZ 通带最大衰减： 3dB 阻带最小衰减： 20dB 抽样频率： 1KHZ 3、带阻滤波器的系数计算： 首先确定所要求的带阻滤波器是否为几何对称的： Ω Ωp1Ωpu=2π*50*2π*450=4πΩsu=2π*200*2π*300=4π p12*2.25*104 *6*104 2s1因为ΩΩpu≠Ωs1Ωsu，而且设计的是带阻滤波器，所以需要调整这个带阻滤波器的通带上下边界频率。 因为： Ω0 ，Ω02=ΩΩp1 =Ωpu2s1Ωsu，则 Ωp1 = 2 π*200 *2 π*300 Ω2 π*450p1 用 Ωp1 值