m=-log（2）/（log（cosd（70）））;?????????????????????????????????????????????????????????????????
lx=5;?ly=5;?lz=3;????????????????
h=2.15;?????????
%[XTYT]=meshgrid（[-lx/4?lx/4][-ly/4?ly/4]）;?
Nx=lx*14;?Ny=ly*14;???????????????
x=linspace（-lx/2lx/2Nx）;
y=linspace（-ly/2ly/2Ny）;
[XRYR]=meshgrid（xy）;
I_0?=?0.73*3600;

XT1=1.25;?YT1=1.25;????
D1=sqrt（（XR-XT1（11））.^2+（YR-YT1（11））.^2+h.^2）
theta=acos（h./D1）????%θ
theta_d=（theta*180）/pi
q=cosd（theta_d）.^m
lux1=（（I_0）*q.*cosd（theta_d））./（D1.^2）;
%?E=（I（0）（m+1）*cos?（AoD）^m）./（2*pi.*D.^2）

XT2=-1.25;?YT2=1.25;?????
D2=sqrt（（XR-XT2（11））.^2+（YR-YT2（11））.^2+h.^2）
theta=ac