非相干全息记录和再现算法

clear all
clc

row=256;col=256;

%%%%%%%%%% 定义坐标 %%%%%%%%%%
xx=linspace（-5.12/25.12/2row）;
yy=linspace（-5.12/25.12/2col）;
[xs1ys1]=meshgrid（xxyy）; % the location of the source point
[xfyf]=meshgrid（xxyy）; % the lens plane
[xmym]=meshgrid（xxyy）; % the SLM plane
[xy]=meshgrid（xxyy）; % the location of the CCD plane
[uv]=meshgrid（xxyy）;   % the reconstruction plane


%%%%%%%%%% 定义各参数 %%%%%%%%%%
zs_1=250;xs=0;ys=0;
f=250;a=520;d1=132;d2=260;
theta1=0; theta2=2/3*pi; theta3=4/3*pi;
lamda=532*（10^（-6））:2*（10^（-6））:542*（10^（-6））;



%%%%%%%%%%% 计算全息图 %%%%%%%%%%
    H1_1=zeros（rowcol）;
    H2_1=zeros（rowcol）;
    H3_1=zeros（rowcol）; 
    

   
for ii=1:row
    for jj=1:col 
         
         if g_1（iijj）~=0
            xs=xs1（iijj）;
            ys=ys1（iij