clear
clc
A=imread（‘C:\wo.jpg‘）;??%從文件讀取圖象數據
B=imresize（A[512?512]）;???%重新設定圖像大小?
C=im2double（rgb2gray（B））;??%?把灰度圖像的數據類型轉換成雙精度浮點類型（有效數字可達15位，精確度更高）
a=imread（‘C:\gaosi.jpg‘）;
b=imresize（a[512512]）;
b=rgb2gray（b）;
c=im2double（b）;?
C0=abs（C）;?????????????????%?設定初始振幅
MB=abs（c）;?????????????????%?目標振幅
figuremesh（MB）;title（‘已知目標振幅分布‘）
for?p=1:600??????????????%?迭代循環100次
????C=fftshift（C）;???%將傅里葉變換的DC分量移到頻譜中心區
????D=ifft2（C）;?????%二維傅里葉逆變換????
????D1=fftshift（D

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

?????文件???????1132??2014-12-15?21:58??Gerchberg-??Saxton（GS）算法.m

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

?????????????????1132????????????????????1