clear;
i=imread（‘Tulips.jpg‘）;
grayim=rgb2gray（i）;%轉為灰度圖%
[width?height]=size（grayim）;
figure（1）;imshow（i）;title（‘原圖‘）;%原圖%
figure（2）;imshow（grayim）;title（‘原灰度圖‘）;%原圖%
f=fft2（grayim）;%dft%
f_shift=log（abs（fftshift（f）））;%把直流分量移到中間%
figure（3）;imshow（f_shift[]）;title（‘DFT結果‘）;colormap（jet（64））?colorbar;%彩色圖直流分量移至中間dft結果%
f1=f;f1（log（abs（f1））>11）=0;%f1為f去低頻分量%
f2=f;count=0;pixels=width*height;%去高頻分量%
for?j=1:width;
????for?k=1:height;
????????if?log（abs（f2（jk）））<8
????????????f2（jk）=0;
????????????count=count+1;
????????end;
????end;
end;
f3=f;f3（11）=0;%去直流分量%

iff=ifft2（f）;%idft%
iff1=ifft2（f1）;%idft%
iff2=ifft2（f2）;%idft%
iff3=ifft2（f3）;%idft%
?iff4=abs（f）;%求幅值%
pyy=angle（f）;%求相位%
p=180/pi*pyy;
figure（4）;imshow（iff[]）;title（‘沒有壓縮重構的圖像‘）;%圖像復原%
figure（5）;i

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

?????文件?????595284??2011-04-26?21:39??dft\Hydrangeas.jpg

?????文件??????72660??2011-05-01?20:01??dft\H傅里葉變換幅度.jpg

?????文件?????117190??2011-05-01?20:01??dft\H傅里葉變換相位.jpg

?????文件??????54885??2011-05-01?20:01??dft\H幅度+T相位.jpg

?????文件?????620888??2009-07-14?12:52??dft\Tulips.jpg

?????文件??????68860??2011-05-01?20:02??dft\T傅里葉變換幅度.jpg

?????文件?????107379??2011-05-01?20:02??dft\T傅里葉變換相位.jpg

?????文件??????51593??2011-05-01?20:01??dft\T幅度+H相位.jpg

?????文件??????90737??2011-03-20?17:06??dft\傅里葉變換.jpg

?????文件???????1255??2011-04-26?21:27??dft\傅里葉變換.m

?????文件??????47600??2011-03-20?17:05??dft\原圖灰度圖.jpg

?????文件??????53453??2011-03-20?17:06??dft\去低頻輪廓.jpg

?????文件??????49399??2011-03-20?17:07??dft\去直流分量.jpg

?????文件????????964??2011-05-01?20:00??dft\幅相變換1.m

?????文件??????48598??2011-03-20?17:06??dft\無壓縮重構.jpg

?????文件??????49922??2011-03-20?17:07??dft\進行圖像壓縮.jpg

?????目錄??????????0??2011-05-01?20:02??dft

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

??????????????2030667????????????????????17