a=imread（‘1.jpg‘）;????????????????????%imread（‘eee.jpg‘）;??????????
I=rgb2gray（a）;????????????%轉為灰度圖像
[m0n0w0]=size（I）;????????????????????????????????????????????????
figure（1）;
subplot（121）;
imshow（I）;
title（‘灰度圖像‘）
f1=double（I）;?????%?數據類型轉換，MATLAB不支持圖像的無符號整型的計算
g1=fft2（f1）;????????%?傅立葉變換
g2=fftshift（g1）;?????%?轉換數據矩陣
[MN]=size（g2）;
nn=2;???????????%?二階巴特沃斯（Butterworth）低通濾波器
d0=100;
m=fix（M/2）;?n=fix（N/2）;
for?i=1:M
???????for?j=1:N
???????????d=sqrt（（i-m）^2+（j-n）^2）;
???????????h=1/（1+0.414*（d/d0）^（2*nn））;??%?計算低通濾波器傳遞函數
???????????result（ij）=h*g2（ij）;
???????end
end
result=ifftshift（result）;
J1=ifft2（result）;
J2=uint8（real（J1））;
J3=imadjust（J2[0.04?0.7][0?1]1）;?
f2=double（J3）;?????%?數據類型轉換，MATLAB不支持圖像的無符號整型的計算
g3=fft2（f2）;????????%?傅立葉變換
g4=fftshi

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-04-26?10:32??源程序\
?????文件???????33047??2010-06-24?12:42??源程序\Form1.frm
?????文件??????695638??2010-06-24?12:42??源程序\Form1.frx
?????文件????????2180??2010-05-24?15:13??源程序\Module1.bas
?????文件?????????190??2010-05-24?15:00??源程序\MSSCCPRJ.SCC
?????文件?????????843??2010-06-24?12:42??源程序\工程1.vbp
?????文件??????????84??2010-06-24?12:45??源程序\工程1.vbw
?????文件?????????629??2010-04-27?09:13??工程2.vbp
?????文件??????????51??2010-04-29?10:27??工程2.vbw
?????文件?????2267134??2012-04-19?15:15??1.bmp
?????文件??????585720??2010-06-15?11:02??1.jpg
?????文件????????1694??2010-04-25?15:00??biaozhun.m
?????文件???????10318??2010-05-09?11:38??Form1.frm
?????文件????????1903??2010-04-27?09:13??Form2.frm
?????文件???????21814??2010-04-27?09:13??Form2.frx
?????文件????????1361??2010-05-06?11:13??Module1.bas
?????文件????????1536??2010-04-22?11:16??Untitled.m
?????文件?????????254??2010-04-25?15:34??Untitled3.m
?????文件?????????880??2010-05-09?11:38??工程1.vbp