新数字信号处理 教学课件 尹为民 21 6.2IIR DF设计(一).pptVIP

6.2 IIR数字滤波器设计 基本概念 模拟低通滤波器的设计 冲激响应不变法 双线性变换法 设计IIR滤波器的频率变换法 一、 基本概念 1 滤波器的技术指标 1 滤波器的技术指标 2 滤波器的主要参数 3 IIR滤波器设计方法 3 IIR滤波器设计方法 3 IIR滤波器设计方法 二、模拟低通滤波器的设计 1 实际低通滤波器的技术指标 1 实际低通滤波器的技术指标 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 4 巴特沃思低通滤波器的实际设计步骤 5 MATLAB实现 [例]试设计巴特沃思低通滤波器，要求通带截止频率fp=2kHz，通带衰减 ? p ? 3dB ，阻带截止频率fs=4KHz ，阻带衰减?s ? 15dB。 [解] 代入 s = s / ?c 得到所需的巴特沃思低通滤波器 其中 ?c = 2? f c = 4000? (rad/s) [例]试设计巴特沃思低通滤波器 Ha(s) ，使其满足下列技术指标： [解] 取N = 6 [num,den]=butter(N,wc,s) 确定阶数为N，3-dB截频为wc(radian/s)的Butterworth filter分子和分母多项式。s 表示模拟域。 [N,wc]=buttord(wp,ws,Ap,As,s) 确定模拟Butterworth filter的阶数N和3-dB截频wc。wc是由阻带参数确定的。 s 表示模拟域。 [z,p,k]=buttap(N) 确定N阶归一化的Butterworth filter 的零点、极点和增益(gain) 例：设计满足下列条件的模拟Butterworth低通滤波器 fp=10kHz, fs=15kHz, ?p=1dB, ?s=15dB Wp=2*pi*10000;Ws=3*pi*10000;Ap=1;As=15; [N,Wc]=buttord(Wp,Ws,Ap,As,s); fprintf(Order of the filter=%.0f\n,N) [num,den] = butter(N,Wc,s); disp(Numerator polynomial); fprintf(%.4e\n,num); disp(Denominator polynomial); fprintf(%.4e\n,den); omega=[Wp Ws];h = freqs(num,den,omega); fprintf(Ap= %.4f\n,-20*log10(abs(h(1)))); fprintf(As= %.4f\n,-20*log10(abs(h(2)))); omega = [0: 200: 40000*pi]; h = freqs(num,den,omega); gain=20*log10(abs(h));plot(omega/(2*pi),gain); xlabel(Frequency in Hz);ylabel(Gain in dB); 例：设计满足下列条件的模拟Butterworth低通滤波器 fp=10kHz, fs=15kHz, ?p=1dB, ?s=15dB ?p=0.9203dB, ?s=15dB 一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。以低通滤波器为例，如图6-2（称容限图）所示， 频率响应有通带、 过渡带及阻带三个范围（而不是理想的陡截止的通带、阻带两个范围）。图中?1为通带的容限，?2为阻带的容限。 式中，ωp, ωs分别为通带截止频率和阻带截止频率，它们都是数字域频率。幅度响应在过渡带（ωs-ωp）中从通带平滑地下降到阻带，过渡带的频率响应不作规定。 当只需要考虑幅度逼近而不考虑相位时，用幅度平方响应来设计时很方便的！ 例如，选择最小均方误差准则。它是指在一组离散的频率｛ωi｝(i=1, 2, …, M)上，所设计出的实际频率响应幅度|H(ejω)|与所要求的理想频率响应幅度|Hd(ejω)|的均方误差ε最小。 此外还可以有其他许多种误差最小的准则，如最大误差最小准则等。 第二步，求在此最佳准则下滤波器系统函数的系数ak, bk。一般是通过不断改变滤波器系数ak、bk，分别计算ε; 最后，找到使ε为最小时的一组系数ak, bk，从而完成设计。这种设计需