%用fwind2函数设计近似圆对称二维带通滤波器
%通带频率为[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)
%由fspecial函数创建二维高斯窗函数
win = fspecial('gaussian',21,2);
win = win ./ max(win(:));
% 窗函数的最大值为1
mesh(win)
%用上面生成的窗函数设计滤波器
h = fwind2(Hd,win);
freqz2(h)
MATLAB图像处理—孙兆林编著—清华大学出版社