数字信号处理22-6.2IIR DF设计(二).pptVIP

1 冲激响应不变法的变换原理 为了能保持AF的频率响应特性，s平面的虚轴要映射到z平面的单位圆上 为了保持系统的稳定性，s左半平面上的点要映射到单位圆内部 （1）映射要求 （2）变换原理 DF的 h(n)模仿AF的 ha(t)。 T为抽样周期 * 1 冲激响应不变法的变换原理 s 平面 ? ? z = r e j? r ? Re[z] Re[s] j Im[z] j Im[s] ? /T －? /T －3? /T 3? /T 虚轴(? = 0) ←→ 单位圆(r = 1) 右半平面(? 0) ←→ 单位圆外(r 1) z 平面 多值映射 左半平面(? 0) ←→ 单位圆内(r 1) r e j? = e ?T e j ?T s = ? + j ? z = e sT * 1 冲激响应不变法的变换原理 DF的频率响应是AF频率响应的周期延拓 （3）混叠失真 * 1 冲激响应不变法的变换原理 若AF是限带的，且最高频率满足莱奎斯特采样定理，则不产生混叠失真; 任何实际AF都不是限带信号，故都产生混叠失真; 冲激响应不变法只适于带通和低通滤波器。 （3）混叠失真 * 1 冲激响应不变法的变换原理 （4）优缺点 缺点：存在频谱混叠，故不能用冲激响应不变法设计高通、带阻等滤波器。 优点：数字滤波器和模拟滤波器的频率关系为线性 w =WT * 2 模拟滤波器的数字化方法 思路 Ha(s) ha(t) h(n) H(z) 若H(s)具有单阶级点，则该方法较为简便 ILT t=nT ZT * 2 模拟滤波器的数字化方法 S平面极点与z平面极点存在映射z = esT；零点没有这种代数对应关系 为使增益稳定，作如下修正 h(n) = T ha(nT) 于是有 * [例]若模拟滤波器系统函数为 利用冲激响应不变法，设计IIR数字滤波器。 * * 直接代入公式，有 设 T=1，有 * * Ha(j?)和 H(e j?)的幅度谱 、 如下： * * 3 冲激响应不变法的设计步骤 1. 将数字滤波器的频率指标{?}转换为 模拟滤波器的频率指标{?} 2. 由模拟滤波器的指标设计模拟滤波器的H(s)。 3. 利用冲激响应不变法，将H(s)转换H(z)。 ?p,?s ?p,?s H(s) H(z) w=WT 设计模拟 滤波器 冲激响应不变 * 例：利用AF-BW filter及冲激响应不变法设计一DF，满足 ?p=0.2p, ?s=0.6p, ?p?2dB, ?s?15dB。 解： (1) 将数字低通指标转换成模拟低通指标w=W T. ?p=0.2p/T, ?s=0.6p/T, ?p?2dB, ?s?15dB (2) 设计模拟低通滤波器 （BW型） =2 =0.801 3 /T * * 例：利用AF-BW filter及冲激响应不变法设计一DF，满足 ?p=0.2p, ?s=0.6p, ?p?2dB, ?s?15dB。 解： (3) 将模拟低通滤波器转换成数字低通滤波器 极点为s1=(-0.567 8 + 0.565 4j)/T, s2=(-0.567 8 - 0.565 4j)/T 利用 可得DF的系统函数为 * * 4 MATLAB实现 [numd,dend] = impinvar(num,den,Fs) num,den：AF分子、分母多项式的系数向量 Fs=1/T：抽样频率 numd,dend：DF分子、分母多项式的系数向量 * 例：利用AF-BW filter及冲激响应不变法设计一DF，满足 ?p=0.2p, ?s=0.6p, ?p?2dB, ?s?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; Fs=1; %Sampling frequency(Hz) %Analog Butterworth specfication Wp=wp*Fs; Ws=ws*Fs; %determine the order of AF filter [N,Wc]=buttord(Wp,Ws,Ap,As,s); %determine the AF-BW filter [numa,dena]=butter(N,Wc,s); %determine the D