theta?=?30;
ml=-log10（2）/log10（cosd（theta））;
P_LED=20;
nLED=60;
P_total=nLED*nLED*P_LED;
Adet=1e-4;
Ts=1;
index=1.5;
FOV=70;
G_Con=（index^2）/（sind（FOV）.^2）;
lx=5;?ly=5;?lz=3;
h=2.15;
i=2000;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
xt1=-0.5;
yt1=-0.85;

Nx=lx*5;?Ny=ly*5;
x=linspace（-lx/2lx/2Nx）;
y=linspace（-ly/2ly/2Ny）;
[XRYR]=meshgrid（xy）;
D1=sqrt（（XR-xt1）.^2+（YR-yt1）.^2+h^2）;
cosphi_A1=h./D1;
receiver_angle=acosd（cosphi_A1）;
H_A1=（i.*cosphi_A1.^（ml-1））./（D1.^2）;
P_rec_A1=H_A1;
P_rec_A2=fliplr（P_rec_A1）;

xt3=0;
yt3=1;
Nx=lx*5;?Ny=ly*5;
x=linspace（-lx/2lx/2Nx）;
y=linspace（-ly/2ly/2Ny）;
[XRYR]=meshgrid（xy）;
D3=sqrt（（XR-xt3）.^2+（YR-yt3）.^2+h^2）;
cosphi_A3=h./D3;
receiver_angle3=acosd（cosphi_A3）;
H_A3=（i.*cosphi_A3.^（ml-1））./（D3.^