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

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，成像物體

%%
%傳輸過(guò)程。
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））‘