clc;?clear;
%讀圖
%-----------------------------
A?=?imread（‘3.jpg‘）;
TA?=?rgb2gray（A）;?%轉灰度
%?TA?=?imread（‘pout.tif‘）;
%-----------------------------

%FFT
%-----------------------------
FA?=?fft2（TA）;????%FFT2
MFA?=?abs（FA）;????%幅度取絕對值
AFA?=?angle（FA）;??%相位
%-----------------------------

%IFFT
%-----------------------------
FC?=?MFA.*exp（j*AFA）;%反變換
C?=?ifft2（FC）;
NFC?=?1.*exp（j*AFA）;?%無幅度譜，純相位譜??subplot（337）----灰度信息丟失
NC?=?ifft2（NFC）;
NNFC?=?MFA.*exp（j*1）;%無相位譜，純幅度譜??subplot（338）----位置信息丟失
NNC?=?ifft2（NNFC）;
%-----------------------------

%繪圖
%-----------------------------
figure;
set?（gcf‘Position‘[800120800800]?‘color‘‘w‘）????????%xy???w?h
?%變換
subplot（331）;imshow（TA）;?????????????????title（‘源圖像‘）;???%原圖
subplot（332）;imshow（fftshift（MFA/10000））;title（‘幅度譜‘）;???%幅度譜
subplot（333）;imshow（（AFA+pi）/（2*pi））;????title（‘相位譜‘）;???%相位譜
?%反變換
subplot（334）;imshow（abs（C）/250）;??title（‘反變換源圖像‘）;???%反變換圖像
subplot（337）;imshow（abs（NC）*100）;?title（‘純相位譜重建‘）;???%純相位譜圖像重建
subplot（338）;imshow（abs（NNC）/290）;title（‘純幅度譜重建‘）;???%純幅度譜圖像重建

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

?????文件???????1189??2017-12-01?19:16??根據相位譜和幅度譜重建圖像.m

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

?????????????????1189????????????????????1