MATLAB用fwind1函数去设计一个近似圆对称的二维带通滤波器

%用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窗')

MATLAB图像处理—孙兆林编著—清华大学出版社