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 , Ws1 2.9203rad/s, Ws2 2.9603rad/s, As ?10dB。 Wp [2.813,2.9880];Ws [2.9203,2.9603]; Ap 1;As 10; [N,Wc] buttord Wp/pi,Ws/pi,Ap,As [numd,dend] butter N,Wc,stop 运行结果 N 2 numd [0.9522, 3. 7327, 5.5624, 3.7327, 0.9522] dend [1.0000, 3.8242 , 5.5601, 3.6412, 0.9067] 利用IIR数字滤波器处理脉搏信号 ① 中医或西医领域的多种疾病诊断分析（心血管疾病等） ② 疲劳状态预测（司机疲劳驾驶等） ③ 情绪变化分析（学生、孕妇情绪监测） ④ 心理学分析（测谎、犯罪等心理研究） ① ② ③ ④ 采集脉搏信号，经数字滤波、特征提取等