clc
clear
t=3.2e-006;%取樣間隔1
theta=2;?%入射角
W=632.8e-009;%為再現光波長
z=1;
I_1=double（imread（‘G:\quanxi1.bmp‘））;
I=I_1-mean（mean（I_1））;
[mn]=size（I_1）;
[XY]=meshgrid（1:m1:n）;
dx=W*z/（t*m）;%像面上的取樣間隔
C=exp（（i*2*pi/W）*（sin（theta*pi/180）*t.*Y））;%C為再現光
A=（exp（i*2*pi*z/W）/（i*W*z））*exp（（i*pi/（W*z））*（（dx）^2）.*（（X-m/2）.^2+（Y-n/2