%基本參數；
clc;
clear?all;
lamda=532e-9;
z=10;???????????????????????????????????????%遠場。
D=5e-3;?????????????????????????????????????%光源尺寸
deltax=lamda*z/D;???????????????????????????%散斑尺寸1.1mm，適宜成像的物體的細節尺寸3mm
%%
%構建物面（Dr面）尺寸50mm，步長1mm，步數50。Dt緊貼后物面
%構建點擴散函數矩陣
dx=1e-3;????????????????????????????????????%物面（Dr面）步長
M=50;???????????????????????????????????????%物面（Dr面）步數
ds=0.1e-3;??????????????????????????????????%光源面步長
M0=50;??????????????????????????????????????%光源面步數
N=M*M;??????????????????????????????????????%物面（Dr面）像素數???????????????????????????
beta=0.8;???????????????????????????????????%采樣率
K=N*beta;???????????????????????????????????%采樣數
source=zeros（M0M0K）;??????????????????????%方形光斑，尺寸5mm，步長0.1mm步數50每一個獨立子光源振幅在0-1之間隨機分布，相位在0-2pi之間隨機分布；
Lp=lamda*z/ds;??????????????????????????????%獨立子光源尺寸為ds，有效成像周期為53.2mm，成像物體應該小于該尺寸；如果不想看到偽像，Dr面尺寸也應該小于該尺寸。

D1=zeros（MM0）;
for?m=1:M?????
????n=1:M0?????????
????D1（mn）=?exp（1j*（pi/（lamda.*z））*（（m-M/2）*dx-（n-M0/2）.*ds）.^2）;??
end
D2=zeros（M0M）;
for?m=1:M0??????????
????n=1:M
????D2（mn）=exp（1j*（pi/（lamda.*z））*（（m-M0/2）*ds-（n-M/2）.*dx）.^2）;
end

%%
%成像物體：雙縫，縫距5mm
zero1=zeros（1550）;
zero2=zeros（2022）;
zero3=zeros（204）;
one=ones（201）;
obj=[zero1;zero2?one?zero3?one?zero2;zero1];
imwrite（obj‘obj.bmp‘）;?
x=reshape（obj25001）;????????????????????????????????%N*1，成像物體

%%
%傳輸過程。
Ir?=?zeros（MMK）;

%K次采樣
for?h=1:K
??source?=rand（M0M0）.*exp（1j*2*pi*rand（M0M0））;?????%贗熱光源振幅分布。
?%參考光路
??Er_k=D1*source*D2;
??Ir_k=abs（Er_k）.^2;??
??Ir（::h）=Ir_k;
end

I1=Ir（::1）
figure（1）;
imagesc（I1/max（max（I1）））;
title（‘第一幅散斑‘）;
imwrite（I1/max（max（I1））‘