設(shè)計帶通濾波器
[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);

h1=fwind1(Hd,hamming(21));%一維漢明窗
figure,freqz2(h1)
h2=fwind1(Hd,bartlett(21));
figure,freqz2(h2)
h3=fwind1(Hd,hanning(21));
figure,freqz2(h3)
h4=fwind1(Hd,blackman(21));
figure,freqz2(h4)


二維高斯窗設(shè)計近似圓對稱的帶通濾波器
[f1,f2]=freqspace(21,'meshgrid');
Hd=ones(21);
r=sqrt(f1.^2+f2.^2);
Hd((r<0.1)|(r>0.5))=0;%通頻帶為0.5－0.1
colormap(jet(64));
mesh(f1,f2,Hd);
%生成旋轉(zhuǎn)對稱的二維高斯窗函數(shù)
win=fspecial('gaussian',21,2);
win=win./max(win(:));%將窗口的最大值歸1
mesh(win)
h1=fwind2(Hd,win);%用前一步生成的窗函數(shù)設(shè)計濾波器
freqz2(h1)


二維拉普拉斯－高斯窗設(shè)計近似圓對稱的高通濾波器
[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);
win=fspecial('log',21,2);
win=win./max(win(:));
mesh(win);
h2=fwind2(Hd,win);
freqz2(h2)