双线性变换法.ppt
《双线性变换法.ppt》由会员分享，可在线阅读，更多相关《双线性变换法.ppt（33页珍藏版）》请在冰豆网上搜索。
广州大学物理与电子工程学院广州大学物理与电子工程学院第四章第四章IIRIIR数字滤波器的设计数字滤波器的设计4.44.4双线性变换法双线性变换法主要内容主要内容一、一、双线性变换法的基本原理双线性变换法的基本原理二、二、双线性变换法的设计方法双线性变换法的设计方法重点与难点重点与难点重点重点1、双线性变换法的基本原理双线性变换法的基本原理2、双线性变换法、双线性变换法的的设计方法设计方法难点难点1、双线性变换法的基本原理双线性变换法的基本原理问题的提出问题的提出如何将模拟滤波器转变为数字滤波器如何将模拟滤波器转变为数字滤波器?
1.1.脉冲响应不变法脉冲响应不变法脉冲响应不变法脉冲响应不变法2.2.双线性变换法双线性变换法双线性变换法双线性变换法WWp,WWsp,sH（s）H（z）频率频率变换变换设计模拟设计模拟滤波器滤波器AF到到DF的转换的转换问题的提出问题的提出采用脉冲响应不变法采用脉冲响应不变法采用脉冲响应不变法采用脉冲响应不变法DFDF的频谱有混叠的频谱有混叠的频谱有混叠的频谱有混叠解决方法解决方法解决方法解决方法采用双线性变换法采用双线性变换法采用双线性变换法采用双线性变换法上节例题上节例题上节例题上节例题利用模拟利用模拟BW滤波器及脉冲响应不变法设计一滤波器及脉冲响应不变法设计一数数字滤波器字滤波器，满足，满足WWp=0.2pp,WWs=0.6pp,Ap2dB,As15dB。
00.10.20.30.40.50.60.70.80.91-18-16-14-12-10-8-6-4-20NormalizedfrequencyGain,dBAs=14.2dB1、双线性变换法的基本原理、双线性变换法的基本原理基本思想：
基本思想：
基本思想：
基本思想：
将将将将非带限非带限非带限非带限的模拟滤波器映射为最高频的模拟滤波器映射为最高频的模拟滤波器映射为最高频的模拟滤波器映射为最高频率为率为率为率为/T/T的的的的带限带限带限带限模拟滤波器。
模拟滤波器。
模拟滤波器。
模拟滤波器。
模拟频率与数字模拟频率与数字频率的关系为：
频率的关系为：
再线性变换到数字频率。
再线性变换到数字频率。
再线性变换到数字频率。
再线性变换到数字频率。
ss域到域到域到域到zz域的映射关系域的映射关系域的映射关系域的映射关系1、双线性变换法的基本原理、双线性变换法的基本原理双线性变换双线性变换1、双线性变换法的基本原理、双线性变换法的基本原理稳定性分析稳定性分析稳定性分析稳定性分析令令s=+j，则有：
，则有：
1、双线性变换法的基本原理、双线性变换法的基本原理稳定性分析稳定性分析稳定性分析稳定性分析1）1）0,0,|z|0,0,|z|1S域域虚轴虚轴映射到映射到z域域单位圆上单位圆上S域域右半平面右半平面映射到映射到z域域单位圆外单位圆外因果因果因果因果、稳定的稳定的稳定的稳定的AFAF系统映射为因果系统映射为因果系统映射为因果系统映射为因果、稳定的稳定的稳定的稳定的DFDF系统系统系统系统1、双线性变换法的基本原理、双线性变换法的基本原理数字频率数字频率数字频率数字频率WWWW变换到变换到变换到变换到模拟频率模拟频率模拟频率模拟频率）2/tan（2W=TwWppWsWpwsw1、双线性变换法的基本原理、双线性变换法的基本原理双线性变换法的优缺点双线性变换法的优缺点双线性变换法的优缺点双线性变换法的优缺点缺点：
缺点：
缺点：
缺点：
幅度响应不是幅度响应不是常数时会产生幅度失真常数时会产生幅度失真优点优点优点优点：
无混叠无混叠）2/tan（2W=TwWp）（WjeH）（wjHpWsWWpwsww2、双线性变换法的设计方法、双线性变换法的设计方法第一步，将数字滤波器的频率指标第一步，将数字滤波器的频率指标WWk转换为模拟滤转换为模拟滤波器的频率指标波器的频率指标k第二步，由模拟滤波器的指标设计模拟滤波器的第二步，由模拟滤波器的指标设计模拟滤波器的H（s）第三步，利用双线性变换法，将第三步，利用双线性变换法，将H（s）转换转换H（z）。
双线性变换法的设计步骤：
双线性变换法的设计步骤：
2、双线性变换法的设计方法、双线性变换法的设计方法WWp,WWsp,sH（s）H（z）设计模拟设计模拟滤波器滤波器双线性变换双线性变换例例例例1:
1:
利用利用BW型模拟低通滤波器和双线性变换法设计满型模拟低通滤波器和双线性变换法设计满足指标足指标WWp=pp/3，Ap=3dB，N=1的数字低通滤波器，并与的数字低通滤波器，并与脉冲响应不变法设计的脉冲响应不变法设计的DF比较。
比较。
解：
解：
解：
解：
设双线性变换中的参数为设双线性变换中的参数为T
（1）将将DF的频率指标转换为的频率指标转换为AF的频率指标的频率指标
（2）设计设计3dB截频为截频为p的一阶的一阶BW型模拟低通滤波器，即型模拟低通滤波器，即N=1,c=p故故例例例例1:
1:
利用利用BW型模拟低通滤波器和双线性变换法设计型模拟低通滤波器和双线性变换法设计满足指标满足指标WWp=pp/3，Ap=3dB，N=1的数字低通滤波器，的数字低通滤波器，并与脉冲响应不变法设计的并与脉冲响应不变法设计的DF比较。
比较。
解：
解：
解：
解：
（3）用双线性变换法将模拟滤波器转换为数字滤波器用双线性变换法将模拟滤波器转换为数字滤波器结论：
结论：
结论：
结论：
参数参数T的取值和最终的设计结果无关。
的取值和最终的设计结果无关。
为简单起见，一般取为简单起见，一般取T=2例例例例1:
1:
利用利用BW型模拟低通滤波器和双线性变换法设计型模拟低通滤波器和双线性变换法设计满足指标满足指标WWp=pp/3，Ap=3dB，N=1的数字低通滤波器，的数字低通滤波器，并与脉冲响应不变法设计的并与脉冲响应不变法设计的DF比较。
比较。
解：
解：
解：
解：
双线性变换法设计的双线性变换法设计的DF的系统函数为：
的系统函数为：
脉冲响应不变法设计的脉冲响应不变法设计的DF的系统函数为：
的系统函数为：
令令z=ejWW，可分别获得两者的幅度响应。
，可分别获得两者的幅度响应。
例例例例1:
1:
利用利用BW型模拟低通滤波器和双线性变换法设计型模拟低通滤波器和双线性变换法设计满足指标满足指标WWp=pp/3，Ap=3dB，N=1的数字低通滤波器，的数字低通滤波器，并与脉冲响应不变法设计的并与脉冲响应不变法设计的DF比较。
比较。
脉冲响应不变法脉冲响应不变法双线性变换法双线性变换法脉冲响应不变法脉冲响应不变法存在频谱混叠，所设计的存在频谱混叠，所设计的DF不满足给定指标。
不满足给定指标。
双线性变换法双线性变换法不存在频不存在频谱混叠，所设计的谱混叠，所设计的DF满足满足给定指标。
给定指标。
3dB3dB例例例例1:
1:
利用利用BW型模拟低通滤波器和双线性变换法设计型模拟低通滤波器和双线性变换法设计满足指标满足指标WWp=pp/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）;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;例例例例22：
利用利用BW模拟滤波器及模拟滤波器及双线性变换法双线性变换法设计一低通设计一低通DF，满足，满足WWp=0.2pp,WWs=0.6pp,Ap2dB,As15dB解：
解：
解：
解：
（1）将数字低通指标转换成模拟低通指标，取将数字低通指标转换成模拟低通指标，取T=2Ap2db,As15db
（2）设计模拟低通滤波器设计模拟低通滤波器（BW型）型）=2=0.5851例例例例22：
利用利用BW模拟滤波器及模拟滤波器及双线性变换法双线性变换法设计一低通设计一低通DF，满足，满足WWp=0.2pp,WWs=0.6pp,Ap2dB,As15dB解：
解：
解：
解：
（3）用用双线性变换法将模拟低通滤波器转换成数字低通滤波器双线性变换法将模拟低通滤波器转换成数字低通滤波器例例例例22：
利用利用BW模拟滤波器及模拟滤波器及双线性变换法双线性变换法设计一低通设计一低通DF，满足，满足WWp=0.2pp,WWs=0.6pp,Ap2dB,As15dB%DesignDFBWlow-passfilterusingimpulseinvariance%DFBWLPspecficationWp=0.2*pi;Ws=0.6*pi;Ap=2;As=15;T=2;Fs=1/T;%Samplingfrequency（Hz）%AnalogButterworthspecficationwp=2*tan（Wp/2）/T;ws=2*tan（Ws/2）/T;%determinetheorderofAFfilterandthe3-dBcutofffrequencyN,wc=buttord（wp,ws,Ap,As,s）%determinetheAF-BWfilternuma,dena=butter（N,wc,s）例例例例22：
利用利用BW模拟滤波器及模拟滤波器及双线性变换法双线性变换法设计一低通设计一低通DF，满足，满足WWp=0.2pp,WWs=0.6pp,Ap2dB,As15dB%determinetheDFfilternumd,dend=bilinear（numa,dena,Fs）%plotthefrequencyresponsew=linspace（0,pi,1024）;h=freqz（numd,dend,w）;plot（w/pi,20*log10（abs（h）;axis（01-500）;grid;xlabel（Normalizedfrequency）;ylabel（Gain,dB）;%computerApAsofthedesignedfilterw=WpWs;h=freqz（numd,dend,w）;fprintf（Ap=%.4fn,-20*log10（abs（h
（1）;fprintf（As=%.4fn,-20*log10（abs（h
（2）;例例例例22：
利用利用BW模拟滤波器及模拟滤波器及双线性变换法双线性变换法设计一低通设计一低通DF，满足，满足WWp=0.2pp,WWs=0.6pp,Ap2dB,As15dBAp=0.3945dBAs=15.0000dB例例例例22：
利用利用BW模拟滤波器及模拟滤波器及双线性变换法双线性变换法设计一低通设计一低通DF，满足，满足WWp=0.2pp,WWs=0.6pp,Ap2dB,As15dB将将双线性变换法与脉冲响应不变法所双线性变换法与脉冲响应不变法所设计设计DF的结果比较。
的结果比较。
双线性变换法双线性变换法Ap=0.3945dBAs=15.0000dB脉冲响应不变法脉冲响应不变法双线性变换法双线性变换法脉冲响应不变法脉冲响应不变法Ap=1.1187dBAs=12.3628dB33、非低通、非低通IIRIIR数字滤波器的设计数字滤波器的设计WWp,WWsp,sH（z）w=W/T脉冲响应不变法脉冲响应不变法双线性变换法双线性变换法H（s）模拟频模拟频率变换率变换设计原型设计原型低通滤波器低通滤波器复复频频率率变变换换注意：
注意：
注意：
注意：
脉冲响应不脉冲响应不变法不能设计变法不能设计
升级会员 