数字信号处理-数字滤波器设计方法论.ppt

* 例： 设计模拟带通滤波器，通带带宽Bw=2π×200rad/s，中心频率Ω0=2π×1000rad/s，通带内最大衰减ap=3dB，阻带Ωsl=2π×830rad/s, Ωsu=2π×1200rad/s，阻带最小衰减as=15dB。 解 (1) 模拟带通的技术要求： Ω0=2π×1000rad/s, ap=3dB Ωsl=2π×830rad/s, Ωsu=2π×1200rad/s as=15dB, Bw=2π×200rad/s; (2) 频率转换： Ωp=1, * （3）设计模拟归一化低通滤波器G(p)，采用巴特沃斯型，有 （4）求模拟带通Hd(s)： * Matlab实现 %Butt_BP wp=2*pi*[905,1105];%905*1105 approximates 1000*1000 ws=2*pi*[830,1200];Ap=3;As=15; [N,wc]=buttord(wp,ws,Ap,As,s); [B,A]=butter(N,wc,s); fk=0:2000/512:2000;wk=2*pi*fk; Hk=freqs(B,A,wk); plot(fk/1000,20*log10(abs(Hk))); xlabel(Frequency(kHz)); ylabel(Magnitude(dB)); axis([0,2,-20,5]); * (3) 低通到带阻的频率变换公式 (5.2.58) 复频域： 频域： (5.2.59) 归一化低通滤波器G(p)频率Ω 高通滤波器Hd(s)频率Ω 通带 [0,1] [-,0] [0,Ωpl] [Ωpu,∞] 阻带 [Ωs,+∞] [-∞,-Ωs] [Ωsl,Ω0] [ Ω0,Ωsl] * 低通与带阻滤波器的幅频特性 * 例： 设计模拟带阻滤波器，其技术要求为： Ωpl=2π×905rad/s, Ωsl=2π×980rad/s, Ωsu= 2π×1020rad/s, Ωpu=2π×1105rad/s, ap=3dB, as=25dB。试设计巴特沃斯带阻滤波器。 解 （1）模拟带阻滤波器的技术要求： Ω20=ΩplΩpu=4π2×1000025, Bw=Ωsu-Ωsl=2π×200; （2）频率转换： Ωp=1, * （3）设计归一化低通滤波器G(p)： （4）带阻滤波器的Hd(s)为： * Matlab实现 %Butt_BS wp=2*pi*[905,1105]; ws=2*pi*[980,1020];Ap=3;As=25; [N,wc]=buttord(wp,ws,Ap,As,s); [B,A]=butter(N,wc,stop,s); fk=0:2000/512:2000;wk=2*pi*fk; Hk=freqs(B,A,wk); plot(fk/1000,20*log10(abs(Hk))); xlabel(Frequency(kHz)); ylabel(Magnitude(dB)); axis([0,2,-30,5]); * 离散时间信号—序列 * 当： N=0时，C0(x)=1 N=1时，C1(x)=x N=2时，C2(x)=2x2-1 N=3时，C3(x)=4x3-3x 可归纳出高阶切比雪夫多项式的递推公式为 CN+1(x)=2xCN(x) - CN-1(x) (5.2.25) CN(x)特性： (1)切比雪夫多项式的过零点在|x|≤1的范围内； (2)当|x|1时，|CN(x)|≤1,在|x|1范围内具有等波纹性； (3)当|x|1时，CN(x)是双曲线函数，随x单调上升。 * N=0,4,5切比雪夫多项式曲线 Observations: CN(x)有N个过零点 N为偶数时，在x正半轴有N/2个过零点 N为奇数时，在x正半轴有(N+1)/2个过零点 CN(x)在x正半轴的过零点对应着|Ha(jΩ)|=1 * 切比雪夫Ⅰ型滤波器幅频特性 N=3 N=4 * |Ha(jΩ)|的性质： x = 0(或Ω = 0) x = 1(或Ω = Ωp) 0 ≤ x ≤1(或0≤Ω≤Ωp)