%例程4－13   利用海明窗設(shè)計Ⅱ型數(shù)字低通濾波器
clear all;
Wp=0.2*pi;
Ws=0.4*pi;
tr_width=Ws-Wp;                        %過渡帶寬度
N=ceil(6.6*pi/tr_width)+1              %濾波器長度
n=0:1:N-1;
Wc=(Ws+Wp)/2;                          %理想低通濾波器的截止頻率
hd=ideal_lp1(Wc,N);                     %理想低通濾波器的單位沖激響應(yīng)
w_ham=(hamming(N))';                    %海明窗
h=hd.*w_ham;                           %截取得到實際的單位脈沖響應(yīng)
[db,mag,pha,w]=freqz_m2(h,[1]);        %計算實際濾波器的幅度響應(yīng)
delta_w=2*pi/1000;
Ap=-(min(db(1:1:Wp/delta_w+1)))        %實際通帶紋波
As=-round(max(db(Ws/delta_w+1:1:501))) %實際阻帶紋波
subplot(221)
stem(n,hd)
title('理想單位脈沖響應(yīng)hd(n)')
subplot(222)
stem(n,w_ham)
title('海明窗w(n)')
subplot(223)
stem(n,h)
title('實際單位脈沖響應(yīng)hd(n)')
subplot(224)
plot(w/pi,db)
title('幅度響應(yīng)(dB)')
axis([0,1,-100,10])

%--------------------------------------------------------
function[db,mag,pha,w]=freqz_m2(b,a)
%濾波器的幅值響應(yīng)（相對、絕對）、相位響應(yīng)
%db:相對幅值響應(yīng)
%mag：絕對幅值響應(yīng)
%pha: 相位響應(yīng)
%w 采樣頻率； 
%b 系統(tǒng)函數(shù)H(z)的分子項（對FIR，b=h）
%a 系統(tǒng)函數(shù)H(z)的分母項（對FIR，a=1）

[H,w]=freqz(b,a,1000,'whole');
H=(H(1:1:501))';
w=(w(1:1:501))';
mag=abs(H);                                 %絕對幅植響應(yīng)
db=20*log10((mag+eps)/max(mag));            %相對幅植響應(yīng)
pha=angle(H);                               %相位響應(yīng)