ch4_4双线性变换法-北京交通大学陈后金教授信号处理解说.ppt

例：试设计满足下列指标的BW型数字带阻滤波器 Wp1=2.8113rad/s, Wp2=2.9880rad/s, Ap?1dB ,  Ws1=2.9203rad/s, Ws2=2.9603rad/s, As ?10dB 解: (4) 由复频率变换将原型模拟低通转换为模拟带阻滤波器 (5) 由双线性变换模拟带阻滤波器转换成数字带阻滤波器 利用MATLAB实现IIR数字滤波器 确定数字滤波器的阶数及3dB截频Wc [N, Wc] = buttord(Wp, Ws, Ap, As) 其中Wp, Ws为归一化角频率。 例Wp=0.1p, 则Wp=0.1 若为带通或带阻滤波器，则 Wp=[Wp1, Wp2]; Ws=[Ws1, Ws2] BW型数字滤波器 利用MATLAB实现IIR数字滤波器 BW型数字滤波器 确定DF系统函数分子、分母多项式 低通 [num,den] = butter(N,Wc) 高通 [num,den] = butter(N,Wc,high) 带通 [num,den] = butter(N,Wc) 其中Wc=[W1, W2] 带阻 [num,den] = butter(N,Wc,stop) 其中Wc=[W1, W2] 利用MATLAB实现IIR数字滤波器 CB I型数字滤波器 确定数字滤波器的阶数及参数Wc [N, Wc] = cheb1ord(Wp, Ws, Ap, As) DF系统函数分子、分母多项式的确定 低通 [num,den] = cheby1(N,Ap,Wc) 高通 [num,den] = cheby1(N, Ap, Wc,high) 带通 [num,den] = cheby1(N, Ap, Wc) 带阻 [num,den] = cheby1(N, Ap, Wc,stop) * 带通、带阻中W=[W1, W2] 利用MATLAB实现IIR数字滤波器 CB II型数字滤波器 确定数字滤波器的阶数及参数Wc [N, Wc] = cheb2ord(Wp, Ws, Ap, As) DF系统函数分子、分母多项式的确定 低通 [num,den] = cheby2(N,As,Wc) 高通 [num,den] = cheby2(N, As, Wc,high) 带通 [num,den] = cheby2(N, As, Wc) 带阻 [num,den] = cheby2(N, As, Wc,stop) * 带通、带阻中W=[W1, W2] 利用MATLAB实现IIR数字滤波器 椭圆(C)型数字滤波器 确定数字滤波器的阶数及参数Wc [N, Wc] = ellipord(Wp, Ws, Ap, As) DF系统函数分子、分母多项式的确定 低通 [num,den] = ellip(N,Ap,As,Wc) 高通 [num,den] = ellip(N,Ap,As,Wc,high) 带通 [num,den] = ellip(N,Ap,As,Wc) 带阻 [num,den] = ellip(N,Ap,As,Wc,stop) * 带通、带阻中W=[W1, W2] 模拟低通原型 (prototype) buttap cheb1ap besselap ellipap cheb2ap 频率变换 lp2lp lp2hp lp2bp lp2bs 综合设计 butter cheby1 cheby2 ellip besself MATLAB中IIR设计的主要函数 IIR 滤波器最小阶数的计算 buttord cheb1ord cheb2ord ellipord 离散化 bilinear impinvar 例：利用MATLAB实现数字带阻滤波器 Wp1=2.8113rad/s, Wp2=2.9880rad/s, Ap?1dB ,