close;
clc;
clear?all;
A?=?zeros（64）;
A（15:2020:40）?=?1;
A（15:5020:25）?=?1;
A（45:5020:40）?=?1;
A（30:3420:35）?=?1;

figure;
imshow（abs（A）[]）;
Fa?=?fft2（fftshift（A））;
Fs?=?fftshift（Fa）;
Am?=?abs（Fs）;%幅度
Ph?=?angle（Fs）;%相位
s?=?11;
cgh?=?zeros（64*s）;
th?=?max（max（abs（Fs）））;
qq?=?th/1.2;
Am（Am>qq）?=?qq;
q?=?1:s;
w?=?（s+1）/2;

for?m?=?1:64
????for?n?=?1:64
?????