%例程4－16    利用漢寧窗設計Ⅰ型數字高通濾波器

clear all;
Wp=0.6*pi;
Ws=0.4*pi;
tr_width=Wp-Ws;                        %過渡帶寬度
N=ceil(6.2*pi/tr_width)              %濾波器長度
n=0:1:N-1;
Wc=(Ws+Wp)/2;                          %理想低通濾波器的截止頻率
hd=ideal_hp1(Wc,N);                     %理想低通濾波器的單位沖激響應
w_han=(hanning(N))';                    %漢寧窗
h=hd.*w_han;                           %截取得到實際的單位脈沖響應
[db,mag,pha,w]=freqz_m2(h,[1]);        %計算實際濾波器的幅度響應
delta_w=2*pi/1000;
Ap=-(min(db(Wp/delta_w+1:1:501)))        %實際通帶紋波
As=-round(max(db(1:1:Ws/delta_w+1))) %實際阻帶紋波
subplot(221)
stem(n,hd)
title('理想單位脈沖響應hd(n)')
subplot(222)
stem(n,w_han)
title('漢寧窗w(n)')
subplot(223)
stem(n,h)
title('實際單位脈沖響應hd(n)')
subplot(224)
plot(w/pi,db)
title('幅度響應(dB)')
axis([0,1,-100,10])