冲激响应不变法设计数字巴特沃思低通滤波器.docxVIP

题目要求：抽样频率fs=10kHz,在频率小于1kHz的通带内，幅度特性下降小于1dB;在频率大于fst=1.5kHz的阻带内,衰减大于15dB. 用冲激响应不变法设计数字巴特沃思低通滤波器程序 %性能指标 clear all fc=1000;omgc=2*pi*fc; fst=1500;omgst=2*pi*fst; fs=10000;Ap=1;As=15; [n,omgc3]=buttord(omgc,omgst,Ap,As,s); [z,p,k]=buttap(n); b0=k*real(poly(z)); a0=real(poly(p)); [H1,Omg0]=freqs(b0,a0); [b1,a1]=lp2lp(b0,a0,omgc3); [b2,a2]=impinvar(b1,a1,fs); [H2,W]=freqz(b2,a2); dbH=20*log10((abs(H2)+eps)/max(abs(H2))); %输出幅度特性曲线 subplot(3,1,1),plot(W/pi,abs(H2)) xlabel(w(/π)); ylabel(幅度); title(幅度特性曲线) set(gca,XTickMode,manual,XTick,[0:0.2:1]); set(gca,YTickMode,manual,YTick,[0:0.2:1]); grid on