clc;
clear all;
theta=60; %半功率角
m=-log10(2)/log10(cosd(theta));
P_LED=0.02;
nLED=60;
P_total=nLED*nLED*P_LED;
Ar=1e-4;
Ts=1;
index=1.5;
FOV=70;
R=0.4; %光电转换系数
Crf=7.8e-8; %??
G_Con=(index.^2)./(sind(FOV).^2);
rho=0.8; %反射系数
lx=5;
ly=5;
h=2.5
h1=0.5;
h2=1;
h3=1.5;
h4=2;
T1=[1,1,h];
T2=[4,1,h];
T3=[1,4,h];
x=linspace(0,lx,6);
y=linspace(0,ly,6);
[XR,YR]=meshgrid(x,y);
%h=0.5时 来自LED1的接收光功率分布
D11=sqrt((XR-T1(1,1)).^2+(YR-T1(1,2)).^2+h1.^2);
cosphi_A11=h1./D11;
receiver_angle=acosd(cosphi_A11);
H11=(m+1).*Ar.*cosphi_A11.^(m+1).