function?[image_out]?=?Bfilter（image_in?D0?N）
%?Butterworth濾波器，在頻率域進行濾波
%?輸入為需要進行濾波的灰度圖像，Butterworth濾波器的截止頻率D0，階數(shù)N
%?輸出為濾波之后的灰度圖像
[row?colu]?=?size（image_in）;
P?=?2?*?row;
Q?=?2?*?colu;
fp?=?zeros（P?Q）;
%對圖像填充0并且乘以（-1）^（x+y）?以移到變換中心
for?i?=?1?:?row
????for?j?=?1?:?colu
????????fp（i?j）?=?double（image_in（i?j））?*?（-1）^（i+j）;
????end
end
%?對填充后的圖像進行傅里葉變換?到頻域
F1?=fft2（fp）;?%不要用?fftshift
%?生成Butterworth濾波函數(shù)，中心在（m+1n+1）
Bw?=?zeros（P?Q）;
a?=?D0^（2?*?N）;%D0?截止頻率，?N?階數(shù)
for?u?=?1?:?P
????for?v?=?1?:?Q
????????temp?=?（u-（row+1.0））^2?+?（v-（colu+1.0））^2;
????????Bw（u?v）?=?1?/?（1?+?（temp^N）?/?a）;
????end
end
%進行濾波
G?=?F1?.*?Bw;
%?反傅里葉變換
gp?=?ifft2（G）;
%?處理得到的圖像
image_out?=?zeros（row?colu?‘uint8‘）;
gp?=?real（gp）;%?取實部
g?=?zeros（row?colu）;
for?i?=?1?:?row
????for?j?=?1?:?colu
????????g（i?j）?=?gp（i?j）?*?（-1）^（i+j）;
????????%取消圖像的乘數(shù)
????end
end
mmax?=?max（g（:））;
mmin?=?min（g（:））;
range?=?mmax-mmin;
for?i?=?1?:?row
????for?j?=?1?:?colu
????????image_out（ij）?=?uint8（255?*?（g（i?j）-mmin）?/?range）;?%重新量化到min?max
????end
end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????158711??2015-11-13?12:32??輸入圖像\lena.jpg
?????文件??????572892??2013-07-14?20:43??輸入圖像\moon.tif
?????目錄???????????0??2018-12-03?23:34??輸入圖像\
?????文件?????5103414??2018-12-03?09:43??輸出圖像\高斯高通.bmp
?????文件?????5103414??2018-12-03?10:00??輸出圖像\頻域拉普拉斯增強.bmp
?????文件?????5103414??2018-12-03?09:44??輸出圖像\理想高通.bmp
?????文件?????5103414??2018-12-03?15:35??輸出圖像\理想低通.bmp
?????文件?????5103414??2018-12-03?09:44??輸出圖像\布特沃斯高通.bmp
?????文件?????5103414??2018-12-03?15:35??輸出圖像\布特沃斯低通.bmp
?????文件?????5103414??2018-12-03?09:09??輸出圖像\實驗2二維傅里葉變換.bmp
?????文件?????5103414??2018-12-03?09:18??輸出圖像\實驗2三維.bmp
?????文件?????5103414??2018-12-03?09:00??輸出圖像\實驗1一維傅里葉變換.bmp
?????目錄???????????0??2018-12-03?23:31??輸出圖像\
?????文件?????5103414??2018-12-03?15:34??輸出圖像\高斯低通.bmp
?????文件?????5103414??2018-12-03?09:56??輸出圖像\頻域高斯低通與空域高斯低通對比.bmp
?????文件????????1379??2018-12-03?23:21??源程序\Exp3_3.m
?????目錄???????????0??2018-12-03?23:27??源程序\
?????文件????????1218??2018-12-03?23:17??源程序\Bfilter.m
?????文件????????1204??2018-12-03?23:19??源程序\Bhighfilter.m
?????文件????????1668??2018-12-03?23:15??源程序\Exp1.m
?????文件????????1315??2018-12-03?23:16??源程序\Exp2.m
?????文件????????1671??2018-12-03?23:26??源程序\Exp2_3D.m
?????文件????????1969??2018-12-03?23:20??源程序\Exp3_1.m
?????文件????????2012??2018-12-03?23:20??源程序\Exp3_2.m
?????文件??????572892??2013-07-14?20:43??源程序\moon.tif
?????文件????????1152??2018-12-03?23:22??源程序\Exp4.m
?????文件?????????886??2018-12-03?23:18??源程序\gauss_highpass.m
?????文件?????????878??2018-12-03?23:17??源程序\gauss_lowpass.m
?????文件?????????302??2018-12-03?23:18??源程序\gaussfiltbox.m
?????文件????????1088??2018-12-03?23:19??源程序\ideal_highpass.m
?????文件?????????943??2018-12-03?23:16??源程序\ideal_lowpass1.m
............此處省略1個文件信息