双线性变换法精要.pptVIP

例2：利用BW模拟滤波器及双线性变换法设计一低通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]=butter(N,wc,s) 例2：利用BW模拟滤波器及双线性变换法设计一低通DF，满足Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB %determine the DF filter [numd,dend]=bilinear(numa,dena,Fs) %plot the frequency response w=linspace(0,pi,1024); h=freqz(numd,dend,w); plot(w/pi,20*log10(abs(h))); axis([0 1 -50 0]);grid; xlabel(Normalized frequency); ylabel(Gain,dB); %computer Ap As of the designed filter w=[Wp Ws];h=freqz(numd,dend,w); fprintf(Ap= %.4f\n,-20*log10( abs(h(1)))); fprintf(As= %.4f\n,-20*log10( abs(h(2)))); 例2：利用BW模拟滤波器及双线性变换法设计一低通DF，满足Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB Ap= 0.3945dB As= 15.0000dB 例2：利用BW模拟滤波器及双线性变换法设计一低通DF，满足Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 将双线性变换法与脉冲响应不变法所设计DF的结果比较。 双线性变换法 Ap= 0.3945dB As= 15.0000dB 脉冲响应不变法 双线性变换法 脉冲响应不变法 Ap= 1.1187dB As= 12.3628dB 3、非低通IIR数字滤波器的设计 Wp,Ws wp,ws H(z) w =W/T 脉冲响应不变法 双线性变换法 H(s) 模拟频 率变换 设计原型 低通滤波器 复频率变换 注意： 脉冲响应不变法不能设计高通和带阻数字滤波器 方法一 3、非低通IIR数字滤波器的设计 方法二 Wp,Ws H(z) w =W/T 脉冲响应不变法 双线性变换法 数字频 率变换 设计原型 低通滤波器 z域 变换 例4：试设计满足下列指标的BW型数字带阻滤波器 Wp1=2.8113rad/s, Wp2=2.9880rad/s, Ap?1dB ,  Ws1=2.9203rad/s, Ws2=2.9603rad/s, As ?10dB 。 解: 脉冲响应不变法不适合设计数字带阻滤波器，因此采用双线性变换法设计。 (1) 将数字带阻滤波器指标转换成模拟带阻滤波器指标 取T=2，利用 得模拟带阻指标为： wp1=6rad, wp2=13rad, ws1=9rad, ws2=1rad1,Ap?1dB, As ?10dB 例4：试设计满足下列指标的BW型数字带阻滤波器 Wp1=2.8113rad/s, Wp2=2.9880rad/s, Ap?1dB ,  Ws1=2.9203rad/s, Ws2=2.9603rad/s, As ?10dB 。 解: (2) 将模拟带阻滤波器指标转换成模拟低通滤波器指标 Ap?1dB, As ?10dB wp1=6rad, wp2=13rad, ws1=9rad, ws2=1rad1,Ap?1dB, As ?10dB 模拟带阻指标 例4：试设计满足下列指标的BW型数字带阻滤波器 Wp1=2.8113rad/s, Wp2=2.9880rad/s, Ap?1dB ,  Ws1=2.9203rad/s, Ws2=2.9603rad/s, As ?10dB 。