function [H_Channel,d_lp,h_vertical,LED_coord,testot ] = VLC_Ch_diffuse_4single(LED_ID,mode,h_receiver)
% VLC信道传输
% 信道模型，3*3*3室内模型，16*4LED
% 只考虑一次反射
%% 房间模型参数
% Radi_half = 30/180*pi;
% m = -1/log2(cos(Radi_half)); % 朗伯辐射系数2
m = 1; % LED朗伯系数，对应半功率角为60°
n = 1; % 室内墙面的漫反射朗伯系数，一般为1
p_wall = 0.6; % 墙壁平均反射系数
A_wall = 1e-2; % 每小块反射面积0.01m2
reflect_x = 0.2 : 0.10 : 2.8; % 反射点高度范围
reflect_y = 0.5 : 0.10 : 2.6; % 每堵墙反射点的坐标
H_led = 3.0;
%% 接收机位置
test_x = 0.45:0.15:2.55; % train 0.4:0.1:2.6 test 0.45:0.15:2.55
test_y = 0.45:0.15:2.55;
% [X,Y] = meshgrid(test_x,test_y);
la = length(test_x);
testot = zeros(la*la,2);
count = 1;
for i = 1:la
for j = 1:la
testot(count,1) = test_x(i);
testot(count,2) = test_y(j);
count = count + 1;
end
end
% h_receiver = 0.4;
h_vertical = H_led - h_receiver;
%% 发射机模型参数
LED_coord = [0.7,0.7; 0.7,2.2; 2.2,2.2; 2.2,0.7]; % LED中心坐标
Xl_1 = 0.65 : 0.05 : 0.85; % 每个大LED灯中，led间隔0.05m
Yl_1 = 0.65 : 0.05 : 0.85; % 10cm的筒灯
Xl_2 = 0.65 : 0.05 : 0.85;
Yl_2 = 2.15 : 0.05 : 2.35;
Xl_3 = 2.15 : 0.05 : 2.35;
Yl_3 = 2.15 : 0.05 : 2.35;
Xl_4 = 2.15 : 0.05 : 2.35;
Yl_4 = 0.65 : 0.05 : 0.85;
d_lp = zeros(la*la,1);
if LED_ID == 1
Xl = Xl_1; % 第一个大灯
Yl = Yl_1; % 每个包含25个灯珠
for i = 1:size(testot,1)
d_lp(i) = sqrt((testot(i,1) - LED_coord(LED_ID,1))^2 + (testot(i,2) - LED_coord(LED_ID,2))^2 + h_vertical^2); % 对应LED中心坐标到pd的距离
end
end
% Led_num = length(Xl)*length(Xl);
if LED_ID == 2
Xl = Xl_2; % 第er个大灯
Yl = Yl_2; % 每个包含25个灯珠
for i = 1:size(testot,1)
d_lp(i) = sqrt((testot(i,1) - LED_coord(LED_ID,1))^2 + (testot(i,2) - LED_coord(LED_ID,2))^2 + h_vertical^2);
end
% Led_num = length(Xl)*length(Xl);
end
if LED_ID == 3
Xl = Xl_3; % 第3个大灯
Yl = Yl_3; % 每个包含25个灯珠
for i = 1:size(testot,1)
d_lp(i) = sqrt((testot(i,1) - LED_coord(LED_ID,1))^2 + (testot(i,2) - LED_coord(LED_ID,2))^2 + h_vertical^2);
end
% Led_num = length(Xl)*length(Xl);
end
if LED_ID == 4
Xl = Xl_4; % 第4个大灯
Yl = Yl_4; % 每个包含25个灯珠
for i = 1:size(testot,1)
d_lp(i) = sqrt((testot(i,1) - LED_coord(LED_ID,1))^2 + (testot(i,2) - LED_coord(LED_ID,2))^2 + h_vertical^2);
end
end
% Led_num = length(Xl)*length(Xl);
% case LED_ID == 5
% Xl = [Xl_1 Xl_2 Xl_3 Xl_4]; % 4个大灯
% Yl = [Yl_1 Yl_2 Yl_3 Yl_4]; % 每个包含25*4个灯珠
% Led_num = (length(Xl)/4)^2*4;
Zl = H_led*ones(length(Xl));
% pt = 0.5; % LED功率W,总功率就是0.5*484W
Aeff = 1e-4; % 接收器有效面积cm
r = 1.5; % 接收器表面折射指数
Ts = 1; % 接收机光滤波器增益
FOV = 1*pi/4; % 接收机视角，45（半角）°;led辐射角60°
g = r^2/sin(FOV)^2; % 接收机光学集中器增益，当入射角小于FOV时为定值，否则为0
T = 0:0.1:25; %响应时间（ns）,采样精度为0.1ns
c = 0.3; % 光速
% k = 0.95; %Nlos链路导致的功率损失
%% mode = 0,只考虑直射；mode = 1，考虑一次反射
if mode == 0
H_Channel = zeros(la^2,length(T));
%% 不考虑漫反射
for loop = 1:size(testot,1)
coordinate = [testot(loop,1) testot(loop,2) h_receiver]; % 接收机坐标
h0 = zeros(size(T));
f = zeros(size(T));
for i = 1 : length(Xl)
for j = 1 : length(Xl)
[d,cos_out,cos_in ] = DirectParameter(Xl(i),Yl(j),Zl(1),coordinate(1),coordinate(2),coordinate(3)); % 计算两点之间的距离及角度
for t = 1:length(T)
f(t) = Dirac(T(t) - round(10*d/c)/10);%单位冲激函数，其中对信号量化处理
end
temp = Aeff*(m+1)/(2*pi)*g*Ts*(cos_out).^m.*cos_in./d.^2*f;% 计算功率
if cos_in < cos(FOV)
temp = 0;%角度大于FOV，光学集中器增益为0
end
h0 = h0 + temp;
% [aaa I] = max(h0)
% figure
% plot(T,h0);
H_Channel(loop,:) = h0;
end
end
end
% figure(1)
% plot(T,h0);
% title('(任一点信道冲激响应(不考虑漫反射)')
% xlabel('时间');
% ylabel('幅度');
end
if mode ==1
docount = 0;
for loop = 1:size(testot,1)
coordinate = [testot(loop,1) testot(loop,2) h_receiver]; % 接收机坐标
%% 不考虑漫反射
h0 = zeros(size(T));
f = zeros(size(T));
for i = 1 : length(Xl)
for j = 1 : length(Xl)
[d,cos_out,cos_in ] = DirectParameter(Xl(i),Yl(j),Zl(1),coordinate(1),coordinate(2),coordinate(3)); % 计算两点之间的距离及角度
for t = 1:length(T)
f(t) = Dirac(T(t) - round(10*d/c)/10);%单位冲激函数，其中对信号量化处理
end
temp = Aeff*(m+1)/(2*pi)*g*Ts*(cos_out).^m.*cos_in./d.^2*f;% 计算功率
if cos_in < cos(FOV)
temp = 0;%角度大于FOV，光学集中器增益为0
end
h0 = h0 + temp;
end
end
% figure(1)
% plot(T,h0);
% title('(任一点信道冲激响应(不考虑漫反射)')
% xlabel('时间');
% ylabel('幅度');
%% 考虑漫反射
h2_e = zeros(length(t),length(t));
for i = 1 : length(Xl)
for j = 1 : length(Xl)
% 一次反射相关参数，发射机Xl(i),Yl(j),Zl(i)，反射点（xp,yp,zp），接收机t(x),t(y),h_receiver
% d1为发射机到反射点距离
xp_e = 6; yp_e = reflect_x; zp_e = reflect_y;
for y = 1 : length(reflect_y)
for x = 1 : length(reflect_x)
[d1_e(x,y),angle_out1_e(x,y),angle_in1_e(x,y) ] = OneReflectParameter_in(Xl(i),Yl(j),Zl(i),xp_e,yp_e(x),zp_e(y)); % 计算两点之间的距离及角度
[d2_e(x,y),angle_out2_e(x,y),angle_in2_e(x,y) ] = OneReflectParameter_out(xp_e,yp_e(x),zp_e(y),coordinate(1),coordinate(2),coordinate(3)); % 计算两点之间的距离及角度
d_temp = d1_e(x,y) + d2_e(x,y);
for t = 1:length(T)
f(t) = Dirac(T(t) - round(10*d_temp/c)/10);%单位冲激函数，其中对信号量化处理
end
h1_e = A_wall*(m+1)/(2*pi)*(cos(angle_out1_e(x,y)/180*pi)).^m.*cos(angle_in1_e(x,y)/180*pi)./d1_e(x,y).^2;% 计算反射点的照度
h2_e_temp = f*Aeff*p_wall*(n+1)/(2*pi)*h1_e*g*Ts*(cos(angle_out2_e(x,y)/180*pi)).^n.*cos(angle_in2_e(x,y)/180*pi)./d2_e(x,y).^2;% 计算照度
if angle_in2_e(x,y)/180*pi > FOV
h2_e_temp = 0;
end
h2_e = h2_e + h2_e_temp;
end
end
end
end
% figure(2)
% plot(T,h2_e);
% title('(任一点信道冲激响应(考虑漫反射)')
% xlabel('时间');
% ylabel('幅度');
h2_w = zeros(length(t),length(t));
for i = 1 : length(Xl)
for j = 1 : length(Xl)
% 一次反射相关参数，发射机Xl(i),Yl(j),Zl(i)，反射点（xp,yp,zp），接收机t(x),t(y),h_receiver
% d1为发射机到反射点距离
xp_w = 0; yp_w = reflect_x; zp_w = reflect_y;
for y = 1 : length(reflect_y)
for x = 1 : length(reflect_x)
[d1_w(x,y),angle_out1_w(x,y),angle_in1_w(x,y) ] = OneReflectParameter_in(Xl(i),Yl(j),Zl(i),xp_w,yp_w(x),zp_w(y)); % 计算两点之间的距离及角度
[d2_w(x,y),angle_out2_w(x,y),angle_in2_w(x,y) ] = OneReflectParameter_out(xp_w,yp_w(x),zp_w(y),coordinate(1),coordinate(2),coordinate(3)); % 计算两点之间的距离及�
- 1
- 2
前往页