%用fwind1函數去設計一個近似圓對稱的二維帶通濾波器
%通帶頻率范圍為[0.1  0.5]
%生成的矩陣Hd包含了期望帶通頻率響應
%使用freqspace可以創建頻率向量f1和f2 
[f1,f2] = freqspace(21,'meshgrid');
Hd = ones(21); 
r = sqrt(f1.^2 + f2.^2);
Hd((r<0.1)|(r>0.5)) = 0;
colormap(jet(64))
mesh(f1,f2,Hd)
%用一維bartlett窗、blackman窗、hamming窗、kaiser窗等設計濾波器
figure;
subplot(221)
h = fwind1(Hd,bartlett(21));
freqz2(h)
title('bartlett窗')
subplot(222)
h = fwind1(Hd,blackman(21));
freqz2(h)
title('blackman窗')
subplot(223)
h = fwind1(Hd,hamming(21));
freqz2(h)
title('hamming窗')
subplot(224)
h = fwind1(Hd,kaiser(21,2.5));
freqz2(h)
title('kaiser窗')