数字信号处理_12.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 双线性变换法设计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.5851 例：利用AF-BW filter及双线性变换法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 解： (3) 用双线性变换法将模拟低通滤波器转换成数字低通滤波器 例：利用AF-BW filter及双线性变换法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB %Design DF BW low-pass filter using impulse invariance %DF BW LP specfication Wp=0.2*pi; Ws=0.6*pi; Ap=2; As=15; T=2;Fs=1/T; %Sampling frequency(Hz) %Analog Butterworth specfication wp=2*tan(Wp/2)/T;ws=2*tan(Ws/2)/T; %determine the order of AF filter and the 3-db cutoff frequency [N,wc]=buttord(wp,ws,Ap,As,s) %determine the AF-BW filter [numa,dena]=butt