数字滤波器设计(3)

4.3. 布莱克曼窗 4.3.1. 程序 Fs=20000; fp=4000; fs=4500; m=[1 1 0 0]; wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; Rp=0.5; Rs=50; Bt=ws-wp;

N0 = ceil (6.6* pi /Bt); N=N0+mod(N0+1,2); wc=(wp+ws)/2/pi;

hn=fir1(N-1,wc,blackman(N)); freqz(hn,1,512);

plot (w,20*log(abs(hn))); grid;

axis ([0 ,1, -1000 , 100]); xlabel ('频率/kHz'); ylabel ('幅值/dB' );

4.3.2. 滤波器h(n)的波形及其幅频响应损耗函数曲线

5. 用频率采样法设计FIR低通滤波器 5.1. 程序

T=input('T=') fp=4000; fs=4500; Fs=20000;

wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; Rp=0.5; As=50; Bt=ws-wp; m=1;

N=ceil((m+1)*2*pi/Bt); N=N+mod(N+1,2); Np=fix(wp/(2*pi/N)); Ns=N-2*Np-1;

Hk=[ones(1,Np+1),zeros(1,Ns),ones(1,Np)]; Hk(Np+2)=T; Ak(N-Np)=T;

thetak=-pi*(N-1)*(0:N-1)/N; hdk=Hk.*exp(j*thetak); hn=real(ifft(hdk)); hw=fft(hn,1024);

wk=2*pi*[0:1023]/1024; hgw=hw.*exp(j*wk*(N-1)/2); Rp=max(20*log10(abs(hgw))); hgmin=min(real(hgw)); As=20*log10(abs(hgmin));

[N,wc]=buttord(wp,ws,Rp,As,'s'); [B,A]=butter(N,wc,'s');

k=0:511;fk=0:20000/512:20000;wk=2*pi*fk; Hk=freqs(B,A,wk);

plot(fk/1000,20*log10(abs(Hk))); grid on;

xlabel('频率/kHz'); ylabel('幅度/dB'); T=1 T =

1

5.2. 滤波器h(n)幅频响应损耗函数曲线

参考文献

[1]高西全，丁玉美.数字信号处理.西安：西安电子科技大学出版社，2008. [2]陈后金，数字信号处理.北京：高等教育出版社，2004.