s109.m

来自「简单的数字图像处理程序」· M 代码 · 共 21 行

% 产生一维汉明窗并计算其频率响应
f=hamming(20);  % 用来产生汉明窗
figure(10),plot(f);     %如图示  汉明窗 
g=fft(f,512*512);
g=reshape(g,512*256,2);
g=fliplr(g);
g=reshape(g,512*512,1);
x=-(512*256-1)/(512*512):1/(512*512):1/2;
figure(1),plot(x,abs(g));   % 如图示 带通滤波器
hd=zeros(20,20);
hd(4:8,4:8)=1;
[f1,f2]=freqspace(20,'meshgrid');
figure(2),mesh(f1,f2,hd);    %如图 带通滤波器
h=fwind1(hd,hamming(20));
figure(3),freqz2(h);      % 如图示 带通滤波器
%  利用所设计滤波器进行图像滤波
I=imread('3.jpg');
I1=rgb2gray(I);
ff=filter2(h,I1);
figure(5),imshow(ff)      %如图示  波滤效果