%用fsamp2函數設計一個近似對稱的二維帶通濾波器，通帶范圍[0.1 0.5]
%創建一個包含期望帶通響應的矩陣Hd
%首先用函數freqspace創建頻率范圍向量f1和f2
subplot(121)
[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)
%設計一個通過這個響應的濾波器
subplot(122)
h = fsamp2(Hd);
freqz2(h)