clear?all;
close?all;
I=imread（‘C:\Documents?and?Settings\Administrator\桌面\lena.bmp‘）;
[heightwidth]=size（I）;
figureimshow（I）title（‘原始圖像‘）;

fftI=fft2（double（I））;
sfftI=fftshift（fftI）;
A=abs（sfftI）;
E=（A-min（min（A）））/（max（max（A））-min（min（A）））*255?;
figureimshow（E）title（‘圖像幅頻‘）;
B=angle（sfftI）;
F=（B-min（min（B）））/（max（max（B））-min（min（B）））*255;
figureimshow（F）title（‘圖像相頻‘）;

C=abs（ifft2（A））;
C=（C-min（min（C）））/（max（max（C））-min（min（C）））*255?;
figureimshow（C）title（‘幅頻反變換的圖像‘）;
for?i=1:height??????????????????????????????%B的范圍由（-pipi）調整至（02pi）
?????for?j=1:width
?????????if?B（ij）<0
?????????????B（ij）=B（ij）+2*pi;
?????????end
????end
end???????????????????
D=ifft2（ifftshift（exp（-j*B）））;

%D=ifft2（exp（-j*B））;
D=（D-min（min（D）））/（max（max（D））-min（min（D）））*255;

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

?????文件?????231322??2010-10-26?08:34??圖像處理\lena.bmp

?????文件????????878??2010-11-19?20:32??圖像處理\第一題\fourier.m

?????文件?????366258??2010-11-19?20:31??圖像處理\第一題\圖像幅頻.bmp

?????文件?????366258??2010-11-19?20:31??圖像處理\第一題\圖像相頻.bmp

?????文件?????366258??2010-11-19?20:30??圖像處理\第一題\幅頻反變換的圖像.bmp

?????文件?????366258??2010-11-19?20:32??圖像處理\第一題\相頻反變換的圖像.bmp

?????文件?????366258??2010-11-19?21:16??圖像處理\第二題\GLPF變換后的圖像.bmp

?????文件???????1137??2010-11-19?21:15??圖像處理\第二題\GPLF.m

?????目錄??????????0??2010-12-03?09:26??圖像處理\第一題

?????目錄??????????0??2010-12-03?09:26??圖像處理\第二題

?????目錄??????????0??2010-12-03?09:26??圖像處理

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

??????????????2064627????????????????????11