可见光通信计算室内光强分布

%%
theta = 70;
% semi-angle at half power
ml=-log10（2）/log10（cosd（theta））;
%Lambertian order of emission
P_LED=20;
%transmitted optical power by individual LED
nLED=60;
% number of LED array nLED*nLED
P_total=nLED*nLED*P_LED;
%Total transmitted power
Adet=1e-4;
%detector physical area of a PD
Ts=1;
%gain of an optical filter; ignore if no filter is used
index=1.5;
%refractive index of a lens at a PD; ignore if no lens is used
FOV=70;
%FOV of a receiver
G_Con=（index^2）/（sind（FOV）.^2）;
%gain of an optical concentrator; ignore if no lens is used
%%
lx=5; ly=5; lz=3;
% room dimension in meter
h=2.15;
%the distance between source and receiver plane
[XTYT]=meshgrid（[-lx/4 lx/4][-ly/4 ly/4]）;
% position of LED; it is assumed all LEDs are located at same point for
% faster simulation
% for one LED simulation located at the central of the room use XT=0 and YT=0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Nx=lx*5; Ny=ly*5;
% number of grid in the receiver plane
x=linspace（-lx/2lx/2Nx）;
y=linspace（-ly/2ly/2Ny）;
[XRYR]=meshgrid（xy）;
D1=sqrt（（XR-XT（11））.^2+（YR-YT（1