仅供参考 locx1和locy1为循环后对应的矩阵
clc;
clear all;
theta=40;%半功率角
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=60;
R=0.4;%??
Crf=7.8*1e-6;%??
G_Con=(index.^2)/(sind(FOV).^2);
rho=0.8; %反射系数
lx=5;
ly=5;
h=0;
T1=[1,1];
T2=[4,1];
T3=[3,4];
%Nx=lx*2;
%Ny=ly*2;
x=linspace(0,lx,26);
y=linspace(0,ly,26);
[XR,YR]=meshgrid(x,y);
D1=sqrt((XR-T1(1,1)).^2+(YR-T1(1,2)).^2);%
H1=(m+1)*Ar./(2*pi.*D1.^2);
Pr1=P_total.*H1.*Ts.*G_Con;
D2=sqrt((XR-T2(1,1)).^2+(YR-T2(1,2)).^2);
H2=(m+1)*Ar./(2*pi.*D2.^2);