example7_6.m

来自「用于图像处理与分析」· M 代码 · 共 20 行

%首先，创建包含了期望的带通响应的矩阵Hd。代码如下：
[f1,f2] = freqspace(21,'meshgrid');
Hd = ones(21); 
r = sqrt(f1.^2 + f2.^2);
Hd(r<0.5) = 0;
colormap(jet(64))
mesh(f1,f2,Hd)
 
%期望频率响应矩阵如图7-13所示
%然后，生成旋转对称的二维拉普拉斯—高斯窗函数
win = fspecial('log',21,2);
% 将窗口的最大值归1
win = win ./ max(win(:)); 
mesh(win)

%二维拉普拉斯—高斯窗函数如图7-14所示
%最后，用前一步生成的窗函数设计滤波器
h2= fwind2(Hd,win);
freqz2(h2)