张军《数字信号处理》ch4-4双线性变换法-end.pptVIP

由上图看出，在零频率附近，模拟角频率ω与数字频率Ω之间的变换关系接近于线性关系；但当ω进一步增加时，Ω增长得越来越慢，最后当ω→∞时，Ω终止在折叠频率= π处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象， 从而消除了频率混叠现象。 但是双线性变换的这个特点是靠频率的严重非线性关系而得到的。由于这种频率之间的非线性变换关系，就产生了新的问题。 首先，一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了； 其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数（这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如下图 所示。 图 双线性变换法幅度和相位特性的非线性映射 例 设计一个一阶数字低通滤波器，3dB截止频率为Ωc =0.25π，将双线性变换应用于模拟巴特沃思滤波器。 解： 数字低通滤波器的截止频率为Ωc=0.25π，相应的巴特沃思模拟滤波器的 3 dB截止频率是ωc，就有 模拟滤波器的系统函数为 由上题可知，T不参与设计。 将双线性变换应用于模拟滤波器，有 双线性变换法设计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 例：用双线性变换法和一阶巴特沃思低通滤波器，设计一个3dB截频为Wp的数字滤波器，并与脉冲响应不变法设计的DF比较。 解：设双线性变换中的参数为T (1) 将DF的频率指标转换为AF的频率指标 (2) 设计3dB截频为wp的一阶BW型模拟低通滤波器，即 N=1, wc = wp 故 例：用双线性变换法和一阶巴特沃思低通滤波器，设计一个3dB截频为Wp的数字滤波器，并与脉冲响应不变法设计的DF比较。 解：设双线性变换中的参数为T (3) 用双线性变换法将模拟滤波器转换为数字滤波器 结论：参数T的取值和最终的设计结果无关。 为简单起见，一般取T=2 例：用双线性变换法和一阶巴特沃思低通滤波器，设计一个3dB截频为Wp的数字滤波器，并与脉冲响应不变法设计的DF比较。 解： 双线性变换法设计的DF的系统函数为 脉冲响应不变法设计的DF的系统函数为 取Wp =0.6p，令z=ejW ，可分别获得两者的幅度响应。 例：用双线性变换法和一阶巴特沃思低通滤波器，设计一个3dB截频为Wp的数字滤波器，并与脉冲响应不变法设计的DF比较。 0 0.6 1 0 0.7 1 Normalized frequency Amplitude 脉冲响应不变法 双线性变换法 Wp =0.6p 脉冲响应不变法存在频谱混叠，所设计的DF不满足给定指标。而双线性变换法不存在频谱混叠，所设计的DF满足给定指标。 ?3dB 例：用双线性变换法和一阶巴特沃思低通滤波器，设计一个3dB截频为Wp的数字滤波器，并与脉冲响应不变法设计的DF比较。 H双(z)和H脉(z)幅度响应比较的MATLAB实现 Wp=0.6*pi; 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.2p, Ws=0.6p, Ap?2dB, As?15dB 解： (1) 将数字低通指标转换成模拟低通指标，取T=2 Ap?2db, As?15db (2) 设计模拟低通滤波器 （BW型） =2 =0.585 1 例：利用AF-BW fil