数字信号处理—基于计算机的方法2014cph9IIR的设计.ppt

73 确定传输函数H(z)的过程称为数字滤波器设计。 num = [0.008 -0.033 0.05 -0.033 0.008]; den = [1.0 2.37 2.7 1.6 0.41]; k = 256; w = 0:pi/(k-1):pi; h = freqz(num, den, w); subplot(2,1,1) plot(w/pi,abs(h)); grid; title(幅度谱); xlabel(\omega/\pi); ylabel(幅度); subplot(2,1,2) plot(w/pi,angle(h)); grid; title(相位谱); xlabel(\omega/\pi); ylabel(相位, 弧度); IIR滤波器设计的基本思想 基于模拟滤波器设计的思想： 1）将数字滤波器的设计转为模拟滤波器的设计。 2）设计满足技术指标的模拟滤波器。 3）将模拟滤波器转换为数字滤波器。 9.1.0 模拟滤波器的设计 (P611, 附录) Butterworth(巴特沃兹)模拟低通滤波器 切比雪夫I型模拟低通滤波器 切比雪夫II型模拟低通滤波器 椭圆低通滤波器 与Butterworth和Chebyshev两种滤波器相比，在同样误差指标下，椭圆滤波器阶数最小。或者说，同样阶数N下，通带到阻带变化最陡峭。 2、由模拟滤波器映射到数字滤波器的条件 9.3 低通IIR滤波器设计 （4）解归一化 Wp = 0.4142; Ws = 1.1708; %模拟指标 [n,Wn] = buttord(Wp, Ws, 0.5, 15, s); %P157 [b,a] = butter(n, Wn, s); %P157 [b,a] = bilinear(b, a, 0.5); %0.5=1/T (T=2) [h,omega] = freqz(b, a, 512); figure(1); plot (omega/pi, abs(h)); figure(2); plot (omega/pi, 20*log10(abs(h))); %第2种实现方法： Wp = 0.25; Ws = 0.55; %数字指标 [n,Wn] = buttord(Wp, Ws, 0.5, 15); [b,a] = butter(n, Wn); [h, omega] = freqz(b, a, 512); figure(1); plot (omega/pi, abs(h)); figure(2); plot (omega/pi, 20*log10(abs(h))); 例: 设计一个二阶IIR陷波滤波器 模拟 数字 例：设计一个Butterworth型的低通IIR数字滤波器G(z)，指标如下: （1）数字截止频率转换为模拟截止频率: （2）求阶数 N 和Ωc: (3)三阶归一化低通Butterworth传输函数: (3dB截止频率Ωc=1) （1）通过查表得到上式；(表格见下页) (3)三阶归一化低通Butterworth传输函数: (3dB截止频率Ωc=1) （1）通过查表得到上式； （2）通过Matlab计算得到上式： N=3; [z, p, k] = buttap(N); [num, den] = zp2tf(z,p,k); 结果：num = [0 0 0 1] den = [1 2 2 1] (denominator) (numerator) Butterworth analog lowpass filter prototype (返回zeros, poles, gain) N=3; [z, p, k] = buttap(N); [num, den] = zp2tf(z,p,k); 计算3dB截止频率为1、阶数为N的归一化模拟Butterworth低通滤波器的传输函数的零点、极点和增益因子。传输函数形式如下： （Zero-pole to transfer function）根据零点、极点和增益因子，计算分子、分母的系数（按s的降幂排列）。 （5）双线性变换，得到数字滤波器的传输函数： buttord返回阶数和3dB截止频率(rad/s) 增益(dB) 幅度 * —— 数字信号处理 —— 李 永 杰 Liyj@uestc.edu.cn 电子科技大学生命科学与技术学院 第9章 IIR数字滤波器的设计 Design of IIR Digital Filter 预备知识 模拟滤波器及其设计方法 数字滤波器的指标 滤波器类型的选择 IIR数字滤波器设计的基本方法 IIR数字滤波器阶数的估计 双线性变换 低通IIR滤波器设计 主要内容 在大多数应用中，设计数字滤波器的关键