%首先，創建包含了期望的頻率響應的矩陣Hd
[f1,f2] = freqspace(21,'meshgrid');
Hd = zeros(21); 
Hd(7:15,7:15) = 1;
Axis([-1 1 –1 1 0 1.2]),
colormap(jet(64))
mesh(f1,f2,Hd)

%期望的濾波器的頻率響應見圖7-16
%然后，設計通過此響應的濾波器
h = fsamp2(Hd);
figure,freqz2(h)