4-20.m

matlab数字信号与应用源代码4

%例程4－20   利用布拉克曼窗设计Ⅱ型数字带通滤波器

clear all;                     
Wpl=0.4*pi;
Wph=0.6*pi;
Wsl=0.2*pi;
Wsh=0.8*pi;
tr_width=min((Wpl-Wsl),(Wsh-Wph));      %过渡带宽度
N=ceil(11*pi/tr_width)+1                %滤波器长度
n=0:1:N-1;
Wcl=(Wsl+Wpl)/2;                        %理想低通滤波器的截止频率
Wch=(Wsh+Wph)/2;
hd=ideal_bp1(Wcl,Wch,N);               %理想低通滤波器的单位冲激响应
w_bman=(blackman(N))';                 %布拉克曼窗
h=hd.*w_bman;                          %截取得到实际的单位脉冲响应
[db,mag,pha,w]=freqz_m2(h,[1]);        %计算实际滤波器的幅度响应
delta_w=2*pi/1000;
Ap=-(min(db(Wpl/delta_w+1:1:Wph/delta_w+1)))   %实际通带纹波
As=-round(max(db(Wsh/delta_w+1:1:501)))   %实际阻带纹波
subplot(221)
stem(n,hd)
title('理想单位脉冲响应hd(n)')
subplot(222)
stem(n,w_bman)
title('布拉克曼窗w(n)')
subplot(223)
stem(n,h)
title('实际单位脉冲响应hd(n)')
subplot(224)
plot(w/pi,db)
title('幅度响应(dB)')
axis([0,1,-100,10])