室内可见光通信光照分布MATLAB代码

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