张军《数字信号处理》ch4-3脉冲响应不变法-end.pptVIP

（4）虽然脉冲响应不变法能保证S平面极点与Z平面极点有这种代数对应关系，但是并不等于整个S平面与Z平面有这种代数对应关系，特别是数字滤波器的零点位置就与模拟滤波器零点位置没有这种代数对应关系，而是随H(s)的极点pl以及系数Al两者而变化。 数字滤波器频率响应幅度还与采样间隔T成反比： 如果采样频率很高，即T很小，数字滤波器可能具有太高的增益，这是不希望的。为了使数字滤波器增益不随采样频率而变化，可以作以下简单的修正，令 h(n)=Tha(nT) 则有: 模拟滤波器的频率响应H(jw)以及数字滤波器的频率响应H(ejW)分别为: 脉冲响应不变法的基本原理 脉冲响应不变法的优缺点 缺点：存在频谱混叠，故不能用脉冲响应不变法设 计高通、带阻等滤波器。 优点：数字滤波器和模拟滤波器的频率关系为线性 W =wT 即一个线性相位的模拟滤波器通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器。 脉冲响应不变法设计DF的步骤 1. 将数字滤波器的频率指标{Wk}转换为 模拟滤波器的频率指标{wk} 2. 由模拟滤波器的指标设计模拟滤波器的H(s)。 3. 利用脉冲响应不变法，将H(s)转换H(z)。 Wp,Ws wp,ws H(s) H(z) w=W/T 设计模拟 滤波器 脉冲响应不变 脉冲响应不变法设计DF的步骤 [numd,dend] = impinvar(num,den,Fs) num,den：AF分子、分母多项式的系数向量 Fs=1/T：抽样频率 numd,dend：DF分子、分母多项式的系数向量 脉冲响应不变法的MATLAB实现 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB。 解： (1) 将数字低通指标转换成模拟低通指标w=W / T. wp=0.2p/T, ws=0.6p/T, Ap?2dB, As?15dB (2) 设计模拟低通滤波器 （BW型） =2 =0.801 3 /T 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 。 解： (3) 将模拟低通滤波器转换成数字低通滤波器 极点为s1=(-0.567 8 + 0.565 4j)/T, s2=(-0.567 8 - 0.565 4j)/T 利用 可得DF的系统函数为 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?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=buttord(wp,ws,Ap,As,s); %determine the 3-db cutoff frequency of BW filter from pass-band specfication wc=wp/(10^(0.1*Ap)-1)^(1/N/2); %determine the AF-BW filter [numa,dena]=butter(N,wc,s); 例：利用AF-BW filter及脉冲响应不变法设计一DF，满足 Wp=0.2p, Ws=0.6p, Ap?2dB, As?15dB 。 %determine the DF filter [numd,dend]=impinvar(numa,dena,Fs); %plot the frequency response w=linspace(0,pi,1024); h=freqz(numd,dend,w); norm=max(abs(h)); numd=numd/norm; plot(w/pi,20*log10(abs(h/norm))); xlabel(Normalized frequency); ylabel(Ga