双线性变换法北京交通大学陈后金教授信号处理课件PPT课件.ppt

近代数字信号处理(Advanced Digital Signal Processing) 电子信息工程学院 信号与图像处理研究室 IIR数字滤波器设计的基本思想 模拟低通滤波器设计 模拟域频率变换 脉冲响应不变法 双线性变换法 IIR数字滤波器的基本结构 利用MATLAB设计IIR DF IIR数字滤波器的设计 问题的提出 双线性变换法的基本原理 双线性变换法设计DF的步骤 双线性变换法 问题的提出 如何将模拟滤波器转变为数字滤波器? 1. 脉冲响应不变法 2. 双线性变换法 问题的提出 采用脉冲响应不变法 DF的频谱有混叠 解决方法 采用双线性变换法 上节例题利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap2dB, As15dB 。 As = 14.2dB 双线性变换法的基本原理 基本思想：将非带限的模拟滤波器映射为最高频率为的带限模拟滤波器 模拟频率与数字频率的关系为 双线性变换法的基本原理 稳定性分析 1) s0, |z|1 S域左半平面映射到z域单位元内 2) s=0, |z|=1 3) s0, |z|1 S域虚轴映射到z域单位圆上 S域右半平面映射到z域单位圆外 因果、稳定的AF系统映射为因果、稳定的DF系统 双线性变换法的基本原理 W和w 的关系 双线性变换法的基本原理 双线性变换法的优缺点 缺点：幅度响应不是常数时会产生幅度失真 优点：无混叠 双线性变换法设计DF的步骤 1. 将数字滤波器的频率指标{Wk}转换为 模拟滤波器的频率指标{wk} 2. 由模拟滤波器的指标设计模拟滤波器的H(s)。 3. 利用双线性变换法，将H(s)转换H(z)。 双线性变换法设计DF的步骤 Wp,Ws wp,ws H(s) H(z) 设计模拟 滤波器 双线性变换 [numd,dend] = bilinear(num,den,Fs) num,den：AF分子、分母多项式的系数向量 Fs：抽样频率 numd,dend：DF分子、分母多项式的系数向量 利用MATLAB 例: 利用BW型模拟低通滤波器和双线性变换法设计满足指标Wp=p/3，Ap=3dB，N=1的数字低通滤波器，并与脉冲响应不变法设计的DF比较。 解：设双线性变换中的参数为T (1) 将DF的频率指标转换为AF的频率指标 (2) 设计3dB截频为wp的一阶BW型模拟低通滤波器，即 N=1, wc = wp 故 例: 利用BW型模拟低通滤波器和双线性变换法设计满足指标Wp=p/3，Ap=3dB，N=1的数字低通滤波器，并与脉冲响应不变法设计的DF比较。 解：设双线性变换中的参数为T (3) 用双线性变换法将模拟滤波器转换为数字滤波器 结论：参数T的取值和最终的设计结果无关。 为简单起见，一般取T=2 例: 利用BW型模拟低通滤波器和双线性变换法设计满足指标Wp=p/3，Ap=3dB，N=1的数字低通滤波器，并与脉冲响应不变法设计的DF比较。 解： 双线性变换法设计的DF的系统函数为 脉冲响应不变法设计的DF的系统函数为 令z=ejW ，可分别获得两者的幅度响应。 例: 利用BW型模拟低通滤波器和双线性变换法设计满足指标Wp=p/3，Ap=3dB，N=1的数字低通滤波器，并与脉冲响应不变法设计的DF比较。 脉冲响应不变法 双线性变换法 脉冲响应不变法存在频谱混叠，所设计的DF不满足给定指标。而双线性变换法不存在频谱混叠，所设计的DF满足给定指标。 3dB 例: 利用BW型模拟低通滤波器和双线性变换法设计满足指标Wp=p/3，Ap=3dB，N=1的数字低通滤波器，并与脉冲响应不变法设计的DF比较。 H双(z)和H脉(z)幅度响应比较的MATLAB实现 Wp=pi/3; b=[1-exp(-Wp)];b1=tan(Wp/2)*[1 1]; a=[1 -exp(-Wp)];a1=[1+tan(Wp/2) tan(Wp/2)-1]; w=linspace(0,pi,512); h=freqz(b,a,w);h1=freqz(b1,a1,w); plot(w/pi,(abs(h)),w/pi,(abs(h1)) ); xlabel(Normalized frequency); ylabel(Amplitude); set(gca,ytick,[0 0.7 1]); set(gca,xtick,[0 Wp/pi 1]); grid; 例：利用AF-BW filter及双线性变换法设计一DF，满足 Wp=0.2