资源描述:
《ch4_4双线性变换法》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数字信号处理(DigitalSignalProcessing)信号与系统系列课程组国家电工电子教学基地IIR数字滤波器设计的基本思想模拟低通滤波器设计模拟域频率变换脉冲响应不变法双线性变换法IIR数字滤波器的基本结构利用MATLAB设计IIRDFIIR数字滤波器的设计问题的提出双线性变换法的基本原理双线性变换法设计DF的步骤双线性变换法问题的提出如何将模拟滤波器转变为数字滤波器?1.脉冲响应不变法2.双线性变换法Wp,Wswp,wsH(s)H(z)频率变换设计模拟滤波器AF到DF的转换问题的提出
2、采用脉冲响应不变法DF的频谱有混叠解决方法采用双线性变换法上节例题利用AF-BWfilter及脉冲响应不变法设计一DF,满足Wp=0.2p,Ws=0.6p,Ap2dB,As15dB。00.10.20.30.40.50.60.70.80.91-18-16-14-12-10-8-6-4-20NormalizedfrequencyGain,dBAs=14.2dB双线性变换法的基本原理基本思想:将非带限的模拟滤波器映射为最高频率为的带限模拟滤波器模拟频率与数字频率的关系为s域到z域的映射关系双线性变
3、换法的基本原理双线性变换双线性变换法的基本原理稳定性分析令s=+j,则有双线性变换法的基本原理稳定性分析1)s<0,
4、z
5、<1S域左半平面映射到z域单位元内2)s=0,
6、z
7、=13)s>0,
8、z
9、>1S域虚轴映射到z域单位圆上S域右半平面映射到z域单位圆外因果、稳定的AF系统映射为因果、稳定的DF系统双线性变换法的基本原理W和w的关系)2/tan(2W=TwWp)(WjeH)(wjHpWsWWpwsww双线性变换法的基本原理双线性变换法的优缺点缺点:幅度响应不是常数时会产生幅度失真优点:无混
10、叠双线性变换法设计DF的步骤1.将数字滤波器的频率指标{Wk}转换为模拟滤波器的频率指标{wk}2.由模拟滤波器的指标设计模拟滤波器的H(s)。3.利用双线性变换法,将H(s)转换H(z)。双线性变换法设计DF的步骤Wp,Wswp,wsH(s)H(z)设计模拟滤波器双线性变换[numd,dend]=bilinear(num,den,Fs)num,den:AF分子、分母多项式的系数向量Fs:抽样频率numd,dend:DF分子、分母多项式的系数向量利用MATLAB例:利用BW型模拟低通滤波器和双线
11、性变换法设计满足指标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的取值和最终的设计结
12、果无关。为简单起见,一般取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不满足给
13、定指标。而双线性变换法不存在频谱混叠,所设计的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)*[11];a=[1-exp(-Wp)];a1=[1+tan(Wp/2)tan(Wp/2)-1];w=linspace(0,pi,512);h=freqz(b,a,w
14、);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/pi1]);grid;例:利用AF-BWfilter及双线性变换法设计一DF,满足Wp=0.2p,Ws=0.6p,Ap2dB,As15dB解:(1)将数字低通指标转换成模拟低通指标,取T=2Ap