clear?all
clear
N=256;
lamda=632.8e-9;%波長
k=2*pi/lamda;
%物
Uo=imread（‘guang.bmp‘）;
Uo=double（Uo（::1））;
[rc]=size（Uo）;
Lo=0.005;%物面尺寸
xo=linspace（-Lo/2Lo/2c）;
yo=linspace（-Lo/2Lo/2r）;
[xo?yo]=meshgrid（xoyo）;%物面坐標網格
Uo=ones（rc）*0.98-Uo/255*0.5;
figure（1）imshow（Uo[01]）;%顯示物（光學）分布
title（‘物（光學）分布‘）;
%===================================================================
%===================================================================
%主透鏡
fo=1.204;%主透鏡焦距
do=1.2041;%物面到主透鏡的距離
D=0.01;
L=r*lamda*do/Lo;
x1=linspace（-L/2L/2c）;
y1=linspace（-L/2L/2r）;
[x1?y1]=meshgrid（x1y1）;%透鏡前表面坐標網格
F0=exp（j*k*do）/（j*lamda*do）*exp（j*k/2/do*（x1.^2+y1.^2））;
F1=exp（j*k/2/do*（xo.^2+yo.^2））;
FU1=（Lo*Lo/r/r）.*fftshift（fft2（Uo.*F1））;
U1=F0.*FU1;%透鏡前表面上的復振幅分布
I1=U1.*conj（U1）;%