（全新）数字信号处理(王震宇 张培珍编 ) 第六章.pptVIP

(3) 对采样，得到数字滤波器的单位脉冲响应序列为 例6.2 求系统函 在采样间隔分别为0.5，0.25，0.02时数字滤波器的幅频响应。 实际系统不可能严格限带，都会存在混叠失真，在 ?处衰减越快，失真越小，如图所示。 (2) 双线性变换法设计的数字滤波器幅频响应曲线形状与原模拟滤波器的幅频特性曲线的形状偏离较大，这是由于变换算法的非线性造成的，ω小时，非线性的影响少一些，所以适合于片断常数滤波器的设计。故双线性变换只能用于设计低通、高通、带通、带阻等选频滤波器。 要使稳定的系统函数Hp(u)变换成稳定的H(z)，u的单位圆内部必须对应于z的单位圆内部， u的单位圆必须对应于z的单位圆。 解 首先利用式(6.40)将数字滤波器的指标转化成模拟指标 根据给定指标可以求解模拟滤波器的系统函数后(前一章已介绍其求解方法)，利用双线性变换法 将其数字化，具体过程可以通过程序求得。 例6.6 试分别用脉冲响应不变法和双线性变化法将下图所示的RC低通滤波器转换成数字滤波器，其中1/RC=1，T=0.5 解： 模拟RC滤波器的频率响应函数为 可以得出对应的系统函数为 例 6 利用脉冲响应不变法转换，数字滤波器的系统函数H1(z)为 利用双线性变换法转换，数字滤波器的系统函数H2(z)为 例 6 clear a=[0 1]; b=[1 1]; freqs(a,b); figure [az,bz]=bilinear(a,b,2); freqz(az,bz); figure [az1,bz1]=impinvar(a,b,2) freqz(az1,bz1); 模拟滤波器幅度和相位特性 例 6 脉冲响应不变法设计数字滤波器 幅频和相频特性 双线性变换法设计数字滤波器 幅频和相频特性 例 6 数字滤波器最优化设计法的特点是适合设计任意幅度特性的滤波器，几种常用的方法： IIR数字滤波器最优化设计法 6 1．零极点累积法 2．频域最小均方误差法 3．时域逼近法 设计数字滤波器的频率变换 (a) 方法一 (b) 方法二 设计数字滤波器的频率变换 6 1．模拟低通 6.3.1 低通变换 6 数字低通 （1）将归一化的Ha1(s)变换成截止频率为Ωc的低通滤波器，低通变换模拟角频率 （2）模拟角频率与数字角频率转化 （3）脉冲响应不变法或双线性变换法， 将模拟低通滤波器化为所需的数字低通滤波器 利用脉冲响应不变法设计一个低通IIR数字滤波器H(z)，T=2性能指标如下： 解 : N=4 Ωc P135 例 6 Matlab程序 wp=0.25*pi; ws=0.55*pi; rp=0.5; rs=15; T=2; Omegap=wp/T; Omegas=ws/T; [N, Omegac]=buttord(Omegap,Omegas,rp ,rs,s); [z,p,k]=buttap(N); [num,den]=zp2tf(z,p,k); [numt,dent] = lp2lp(num,den, Omegac); w=linspace(0,4*Omegac,512); h1=freqs(num,den,w); h2=freqs(numt,dent,w); 例 6 L1=plot(w,20*log(abs(h1)),-); hold on L2=plot(w,20*log(abs(h2)),-.); [azt,bzt]=impinvar(numt,dent,1/T); [h3,w]=freqz(azt,bzt); L3=plot(w,20*log(abs(h3)),:); legend([L1 L2 L3],模拟低通原型, 截止频率Wc模拟低通,截止频率Wc数字低通); 程序运行结果 N=N1=4；Omegac=0.5633 所以滤波器的阶数为4，Ωc=0.5633 例 6 例 6 （1） 确定数字低通滤波器的技术指标： ωp、Ap、ωs、As。 （2） 将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。边界频率ωp和ωs的转换， Ap和As指标不变。 （3） 求得滤波器的阶数，查表归一化模拟滤波器系统函数。 （4）设计过渡模拟低通滤波器 （5） 将模拟滤波器Ha(s)转换成数字低通滤波器系统函数H(z)[脉冲相应不变法或双线性变换法）。 设计IIR数字低通滤波器的步骤 6 设变换前z平面定义为u平面，变换后z平面仍为z平面。 平面间变换关系为 低通变换映射函数(全通函数N=1) 2．数字低通---数字低通 6 数字滤波器低通原型的系统函数Hp(z