%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?高斯帶阻濾波器
%?濾波器的傳遞函數為H（uv）=1-exp（-1/2*（（D^2-D0^2）/（D*W））^2）
%?其中W是頻帶寬度D0是頻帶的中心半徑
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function?ImgOut=BandFilter（ImgInd0W）
s=fftshift（fft2（ImgIn））;???????????????%將灰度圖像的二維不連續Frourier變換的零頻率成分移到頻譜的中心
[MN]=size（s）;??????????????????????%分別返回s的行數到M中，列數到N中
n1=floor（M/2）+1;????????????????????%對M/2進行取整
n2=floor（N/2）+1;????????????????????%對N/2進行取整
for?i=1:M
????for?j=1:N
????????d=sqrt（（i-n1）^2+（j-n2）^2）;?????????%點（ij）到傅立葉變換中心的距離
????????h=1-exp（-1/2*（（d^2-d0^2）/（d*W））^2）;%GBEF濾波函數
????????s1（ij）=h*s（ij）;???????????????????%GBEF濾波后的頻域表示
????????H（ij）=h;
????end
end
figure;
subplot（131）imshow（log（abs（s）+1）[]）;title（‘濾波前圖像傅里葉頻譜取對數‘）;
subplot（132

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????17131??2011-09-08?18:54??帶通帶阻濾波器\1.jpg

?????文件???????3544??2011-09-08?19:19??帶通帶阻濾波器\BandFilter.asv

?????文件???????1200??2011-09-08?19:38??帶通帶阻濾波器\BandFilter.m

?????文件???????3810??2011-09-08?19:29??帶通帶阻濾波器\FilterTest.asv

?????文件????????585??2011-09-08?19:38??帶通帶阻濾波器\FilterTest.m

?????文件?????786486??2011-09-07?17:29??帶通帶阻濾波器\Lena.bmp

?????文件??????22538??2011-09-08?18:21??帶通帶阻濾波器\SinNoisy.jpg

?????目錄??????????0??2011-10-10?11:15??帶通帶阻濾波器

-----------?---------??----------?-----??----

???????????????835294????????????????????8