课程设计（论文）-基于巴特沃斯法数字IIR带阻滤波器的设计精选.doc

基于巴特沃斯法数字IIR带阻滤波器的设计 数字带阻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滤波器的设计步骤 2、带阻滤波器的设计规范： 本设计中分别用h(t)、H(s)、Ha(j)表示模拟滤波器的单位脉冲相应、系统函数、频率响应函数，三者的关系如下： Ha(s) = LT[ha(t)] =edt Ha(j)=FT[ha(t)] =edt 可以用h(t)、H(s)、Ha(j)中任一个描述模拟滤波器，也可以用线性常系数微分方程描述模拟滤波器。但是设计模拟滤波器时，设计指标一般由幅频相应函数|Ha(j)|给出，而模拟滤波器设计就是根据设计指标，求系统函数H(s)。 工程实际中通常用所谓的损耗函数即衰减函数A（）来描述滤波器的幅频响应特性，对归一化幅频响应函数，A（）定义如下： A（）=-20lg|Ha(j)|=-10lg|Ha(j)| dB 模拟带阻滤波器的设计指标参数有、、、、。和分别表示带阻滤波器的通带下边界频率和通带上边界频率；和分别表示带阻滤波器的阻带下边界频率和阻带上边界频率。1和2分别称为通带和阻带波纹幅度 图2、带阻滤波器 带阻滤波器的技术指标要求必须是几何对称的。如果带阻滤波器不是几何对称的，为了使对称的带阻滤波器满足最小阻带衰减要求，在由非对称的带阻滤波器变换为对称的带阻滤波器时，保留非对称带阻滤波器的阻带截止频率调整两个通带截止频率中的一个，调整过程步骤如下： 计算= 计算=，如果,用代替； 如果，计算=，并用代替； 如果A A，选择A=min{ A, A}。 所有类型的变换，包括各自与低通原型的频率变换如表1所示。用频率变换法设计一个模拟滤波器的步骤如下： 1、确定低通、高通、带通和带阻模拟滤波器的技术要求。当要求的带通或带阻滤波器是非几何对称时，根据上面描述的过程使这些技术要求成几何对称的； 2、根据表1中的变换关系，确定归一化低通滤波器的技术要求：通带截止频率为，阻带截止频率为，阻带衰减为A(dB),阻带衰减为A (dB)； 3、根据、、A 和A，用巴特沃斯设计归一化低通滤波器； 4、根据表1最右栏的变换关系，得到要求的非归一化模拟滤波器。 表1 模拟滤波器的频率变换 滤波器类型 归一化低通滤波器H(p)的技术指标要求 要求的滤波器H(s), H(s)= H(p) 低通H(p)低通H(s) P= 低通H(p)高通H(s) P= 低通H(p)带通H(s) P= 低通H(p)带阻H(s) P= 下面我们用巴特沃斯滤波法来设计一个数字IIR带阻型滤波器： 滤波器规范：通带下边界频率： 50HZ 通带上边界频率： 450HZ 阻带下边界频率： 200HZ 阻带上边界频率： 300HZ 通带最大衰减： 3dB 阻带最小衰减： 20dB 抽样频率： 1KHZ 3、带阻滤波器的系数计算： 首先确定所要求的带阻滤波器是否为几何对称的： =2*50*2*450=4*2.25*10 =2*200*2*300=4*6*10 因为，而且设计的是带阻滤波器，所以需要调整这个带阻滤波器的通带上下边界频率。 因为： =，=，则 = 用值代替值，即令==2*133.3HZ。 因此，所要求的带阻滤波器的指标调整为： =2*133.3/1000 rad/s =2*450/1000 rad/s =2*200 /1