资源描述:
《ch4_4双线性变换法.ppt》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、数字信号处理(DigitalSignalProcessing)信号与系统系列课程组国家电工电子教学基地IIR数字滤波器设计的基本思想模拟低通滤波器设计模拟域频率变换脉冲响应不变法双线性变换法利用MATLAB设计IIRDFIIR数字滤波器的设计变换推导W和w的关系双线性变换法设计DF的步骤双线性变换法问题的提出采用脉冲响应不变法00.20.40.60.81-9-6-30NormalizedfrequencyGain,dbWp=0.2pWp=0.6p例如上节用脉冲响应不变法和一阶巴特沃思低通滤波器,设计的3d
2、B截频为Wp的数字滤波器为DF的频谱有混叠解决方法采用双线性变换法双线性变换法的基本原理原理:利用数值积分将模拟系统变换为数字系统用梯形面积近似积分原理:利用数值积分将模拟系统变换为数字系统双线性变换法的基本原理将H(z)和H(s)比较可得双线性变换法的基本原理稳定性分析令s=+j,则有双线性变换法的基本原理稳定性分析1)s<0,
3、z
4、<1s域左半平面映射到z域单位元内2)s=0,
5、z
6、=13)s>0,
7、z
8、>1s域虚轴映射到z域单位圆上s域右半平面映射到z域单位圆外因果、稳定的AF系统映射为因果、稳
9、定的DF系统双线性变换法的基本原理W和w的关系令s=jw,z=ejW,则有双线性变换法的基本原理W和w的关系)2/tan(2W=TwWp)(WjeH)(wjHpWsWWpwsww双线性变换法的基本原理双线性变换法的优缺点缺点:幅度响应不是常数时会产生幅度失真优点:无混叠双线性变换法设计DF的步骤1.将数字滤波器的频率指标{Wk}转换为模拟滤波器的频率指标{wk}2.由模拟滤波器的指标设计模拟滤波器的H(s)。3.利用双线性变换法,将H(s)转换H(z)。双线性变换法设计DF的步骤Wp,Wswp,wsH(s
10、)H(z)设计模拟滤波器双线性变换[numd,dend]=bilinear(num,den,Fs)num,den:AF分子、分母多项式的系数向量Fs:抽样频率numd,dend:DF分子、分母多项式的系数向量利用MATLAB解:设双线性变换中的参数为T(1)将DF的频率指标转换为AF的频率指标(2)设计3dB截频为wp的一阶BW型模拟低通滤波器,即N=1,wc=wp故例:用双线性变换法和一阶巴特沃思低通滤波器,设计一个3dB截频为Wp的数字低通滤波器,并与脉冲响应不变法设计的DF比较。例:用双线性变换法和
11、一阶巴特沃思低通滤波器,设计一个3dB截频为Wp的数字低通滤波器,并与脉冲响应不变法设计的DF比较。解:(3)用双线性变换法将模拟滤波器转换为数字滤波器结论:参数T的取值和最终的设计结果无关。为简单起见一般取T=2例:用双线性变换法和一阶巴特沃思低通滤波器,设计一个3dB截频为Wp的数字低通滤波器,并与脉冲响应不变法设计的DF比较。解:双线性变换法设计的DF的系统函数为脉冲响应不变法设计的DF的系统函数为取Wp=0.6p,令z=ejW,可分别获得两者的幅度响应例:用双线性变换法和一阶巴特沃思低通滤波器,设
12、计一个3dB截频为Wp的数字滤波器,并与脉冲响应不变法设计的DF比较。00.6100.71NormalizedfrequencyAmplitude脉冲响应不变法双线性变换法Wp=0.6p脉冲响应不变法存在频谱混叠,所设计的DF不满足给定指标。而双线性变换法不存在频谱混叠,所设计的DF满足给定指标。3dB例:用双线性变换法和一阶巴特沃思低通滤波器,设计一个3dB截频为Wp的数字低通滤波器,并与脉冲响应不变法设计的DF比较。H双(z)和H脉(z)幅度响应比较的MATLAB实现Wp=0.6*pi;b=[1-e
13、xp(-Wp)];b1=tan(Wp/2)*[11];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('Normalizedfrequency');ylabel('Amplitude');set(gca,'ytick',[00.71]);set(gca,'xtick',[0Wp/
14、pi1]);grid;例:利用AF-BWfilter及双线性变换法设计一DF,满足Wp=0.2p,Ws=0.6p,Ap2dB,As15dB解:(1)将数字低通指标转换成模拟低通指标,取T=2Ap2db,As15dB(2)设计模拟低通滤波器(BW型)=2=0.5851例:利用AF-BWfilter及双线性变换法设计一DF,满足Wp=0.2p,Ws=0.6p,Ap2dB,As15dB解:(3)用双线性变换法将模拟低通滤