[f1,f2] = freqspace(25,'meshgrid');
Hd = zeros(25,25); d = sqrt(f1.^2 + f2.^2) < 0.5;
Hd(d) = 1;
mesh(f1,f2,Hd)
%用ftrans2函数设计一个通带在0.1和0.6间(量化频率,这里1.0对应采样频率的一半或π弧度)的近似圆对称二维通带滤波器。
%ftrans2函数由一FIR滤波器创建二维滤波器
%首先用信号处理工具箱函数remez设计一个一维FIR通带滤波器
subplot(121)
colormap(jet(64))
b = remez(10,[0 0.05 0.15 0.55 0.65 1],[0 0 1 1 0 0]);
[H,w] = freqz(b,1,128,'whole');
plot(w/pi-1,fftshift(abs(H)));
%应用函数ftrans2和缺省 McClellan变换矩阵创建近似圆对称滤波器
subplot(122)
h = ftrans2(b);
freqz2(h)
MATLAB图像处理—孙兆林编著—清华大学出版社