数字信号处理---第六章 无限脉冲响应数字滤波器的设计.pptVIP

第6章 无限脉冲响应数字滤波器的设计 6.1 数字滤波器的基本概念 6.2 模拟滤波器的设计 6.3 用脉冲响应不变法设计IIR数字低通滤波器 6.4 用双线性变换法设计IIR数字低通滤波器 6.5 数字高通、带通和带阻滤波器的设计 6.1 数字滤波器的基本概念 经典滤波器从功能上分又可分为： 低通滤波器（LP):Low pass filter 带通滤波器(BP):Bandpass filter 高通滤波器(HP):High pass filter 带阻滤波器(BS):Bandstop filter 我们通常用的数字滤波器一般属于选频滤波器。假设数字滤波器的传输函数H(e jω)用下式表示： 6.2 模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。 巴特沃斯低通滤波器的幅度平方函数|Ha(jΩ)|2用下式表示： 将幅度平方函数|Ha(jΩ)|2写成s的函数： 为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，而右半平面的N个极点构成Ha(-s)。 Ha(s)的表示式为 四.低通巴特沃斯滤波器的设计步骤 (1)根据技术指标Ωp,αp,Ωs和αs，求出滤波器的阶数N。 (2) 求出归一化极点pk，得到归一化传输函数Ha(p)。 (3)将Ha(p)去归一化。将p=s/Ωc代入Ha(p)，得到实际的滤波器传输函数Ha(s)。 例6.2.1 已知通带截止频率fp=5kHz，通带最大衰减αp=2dB，阻带截止频率fs=12kHz，阻带最小衰减αs=30dB，按照以上技术指标设计巴特沃斯低通滤波器。 解 (1) 确定阶数N。 (2) 求极点 上式分母可以展开成为五阶多项式，或者将共轭极点放在一起，形成因式分解形式。 这里不如直接查表6.2.1简单，由N=5，直接查表得到： 极点：-0.3090±j0.9511,-0.8090±j0.5878; -1.0000 (3) 为将Ha(p)去归一化，先求3dB截止频率Ωc。 按照(6.2.19)式，得到： 6.3 用脉冲响应不变法设计IIR数字低通滤波器 为了保证转换后的H(z)稳定且满足技术要求，对转换关系提出两点要求： (1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。 (2)数字滤波器的频率响应模仿模拟滤波器的频响，s平面的虚轴映射z平面的单位圆，相应的频率之间成线性关系。 设模拟滤波器的传输函数为Ha(s),其单位冲激响应是ha(t) 对ha(t)进行等间隔采样，采样间隔为T，得到： 5 脉冲响应不变法优缺点 缺点 会产生频率混叠，适合低通、带通滤波器的设计，不适合高通、带阻滤波器的设计。 6. 修正公式 假设 没有频率混叠现象，即满足 7.数字滤波器二阶基本节 一般Ha(s)的极点si是一个复数，且以共轭成对的形式出现，一对复数共轭极点形成一个二阶基本节。如果模拟滤波器的二阶基本节的形式为 Bz, Az］=impinvar(B, A) 实现用脉冲响应不变法将分子和分母多项式系数向量为B和A的模拟滤波器系统函数Ha(s)转换成数字滤波器的系统函数H(z)，H(z)的分子和分母多项式系数向量为Bz和Az。 例6.3.1 已知模拟滤波器的传输函数Ha(s)为 用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。 解 首先将Ha(s)写成部分分式： 例6.4.1试分别用脉冲响应不变法和双线性不变法将图6.4.4所示的RC低通滤波器转换成数字滤波器。 解 首先按照图6.4.4写出该滤波器的传输函数Ha(s)为 利用双线性变换法转换，数字滤波器的系统函数H2(z)为 例6.4.2 设计低通数字滤波器，要求在通带内频率低于0.2πrad