%---------------------------------------------------------------------------
% exa070902_fir2.m,  for example 7.9.2 and to test fir2,
%---------------------------------------------------------------------------
clear all;

f=[0 0.19 0.2 0.3 0.31 0.59 0.6 0.8 0.81 1]; 
% 給定頻率軸分點；
m=[0 0 1 1 0 0 1 1 0 0];
% 給定在這些頻率分點上理想的幅頻響應
N1=30;
N2=90;
% 取兩種不同的濾波器長度；
b1=fir2(N1,f,m);
b2=fir2(N2,f,m);
% 得到兩個濾波器；

subplot(311);
stem(b1,'.');grid;
subplot(312);
stem(b2,'.');grid;

M=128;
[h1,w]=freqz(b1,1,M,1);
[h2,w]=freqz(b2,1,M,1);
subplot(313);
plot(w,abs(h1),'b-',w,abs(h2),'g-');grid;