clc;
clear;
close all;
warning off;
addpath 'func\'
%B2的取值范围
B2set = [0.25:0.25:5];
lemdas = [0.25:0.25:5]*1e-5;
for bj = 1:length(B2set)
B2set(bj)
for lj = 1:length(lemdas)
%产生网络图
%User个数
Nu = 50;
Lemdau = 90;
%蒙特卡洛循环次数
MTKL = 10000;
%选择一个典型的用户
UserID = 15; %随便设置一个
for mt = 1:MTKL
RandStream.setDefaultStream(RandStream('mt19937ar','seed',mt));
L = 10e3;
N1 = 30;
N2 = 40;
Lemda1 = 5e-5;
Lemda2 = lemdas(lj);
%Macro层1的相关参数
%泊松分布
x1 = poissrnd(Lemda1*1e7,[1,N1]);x1=L*x1/max(x1);
y1 = poissrnd(Lemda1*1e7,[1,N1]);y1=L*y1/max(y1);
%Micro层2的相关参数
%泊松分布
x2 = poissrnd(Lemda2*1e7,[1,N2]);x2=L*x2/max(x2);
y2 = poissrnd(Lemda2*1e7,[1,N2]);y2=L*y2/max(y2);
%产生符合泊松随机分布的用户的区域位置
xu = poissrnd(Lemdau,[1,Nu]);xu=L*xu/max(xu);
yu = poissrnd(Lemdau,[1,Nu]);yu=L*yu/max(yu);
%参数的设定
epsl = 1e6;%数据包量
alpha1 = 4;
alpha2 = 3.6;
delta2 = -106;
B = 1e3;
Pt1 = 40;
Pt2 = 5;
h = 1;
DeltaB1= 1;
%**************************************************************************
%DeltaB2待求
DeltaB2= B2set(bj);
%计算每个用户的信号的强度
for i = 1:Nu
%针对Macro
%选择最近的一个基站,计算对应的距离
for j1 = 1:N1
dist_tmp1(j1) = sqrt((xu(i)-x1(j1))^2 + (yu(i)-y1(j1))^2);
end
dist1 = min(dist_tmp1);
P1(i) = Pt1*h*DeltaB1*dist1^(-alpha1);
%针对Micro
%选择最近的一个基站,计算对应的距离
for j2 = 1:N2
dist_tmp2(j2) = sqrt((xu(i)-x2(j2))^2 + (yu(i)-y2(j2))^2);
end
dist2 = min(dist_tmp2);
P2(i) = Pt2*h*DeltaB2*dist2^(-alpha2);
%选择较大的一个联结
[V,I] = max([P1(i),P2(i)]);
J(i) = I;
end
%计算得到的J为每个用户对应选择的基站标号
J;
%根据如下规则计算SINR
%定义与 Macro层BS连接的用户集合
U1 = find(J==1);
%定义与 Micro层BS连接的用户集合
U2 = find(J==2);
%计算SINR1和RATE1
%计算SINR2和RATE2
SINR1 = zeros(1,Nu);
SINR2 = zeros(1,Nu);
RATE1 = zeros(1,Nu);
RATE2 = zeros(1,Nu);
DeltaT1 = zeros(1,Nu);
DeltaT2 = zeros(1,Nu);
for i = 1:Nu
%计算SINR1和RATE1
if J(i) == 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for j1 = 1:N1
dist_tmp1(j1) = sqrt((xu(i)-x1(j1))^2 + (yu(i)-y1(j1))^2);
end
for j2 = 1:N2
dist_tmp2(j2) = sqrt((xu(i)-x2(j2))^2 + (yu(i)-y2(j2))^2);
end
[V1,I1]= min(dist_tmp1);
dist1 = V1;
FZ = Pt1*h*dist1^(-alpha1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ind1 = 0;
tmps = [];
for j1 = 1:N1
if (j1 < I1) | (j1 >I1)
ind1 = ind1 + 1;
tmps(ind1) = Pt1*h*dist_tmp1(j1)^(-alpha1);
end
end
FM1 = sum(tmps);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tmps = [];
for j1 = 1:N2
tmps(j1) = Pt2*h*dist_tmp2(j1)^(-alpha2);
end
FM2 = sum(tmps);
SINR1(i) = FZ/(FM1+FM2+10^(delta2/20)/1000);
RATE1(i) = B*log2(1+SINR1(i));
DeltaT1(i) = epsl/RATE1(i);
else
SINR1(i) = 0;
RATE1(i) = 0;
DeltaT1(i) = 0;
end
%计算SINR2和RATE2
if J(i) == 2
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for j1 = 1:N1
dist_tmp1(j1) = sqrt((xu(i)-x1(j1))^2 + (yu(i)-y1(j1))^2);
end
for j2 = 1:N2
dist_tmp2(j2) = sqrt((xu(i)-x2(j2))^2 + (yu(i)-y2(j2))^2);
end
[V2,I2]= min(dist_tmp2);
dist2 = V2;
FZ = Pt2*DeltaB2*h*dist2^(-alpha2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tmps = [];
for j1 = 1:N1
tmps(j1) = Pt1*h*dist_tmp1(j1)^(-alpha1);
end
FM1 = sum(tmps);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tmps = [];
ind2 = 0;
for j1 = 1:N2
if (j1 < I2) | (j1 >I2)
ind2 = ind2 + 1;
tmps(ind2) = Pt2*h*dist_tmp2(j1)^(-alpha2);
end
end
FM2 = sum(tmps);
SINR2(i) = FZ/(FM1+FM2+10^(delta2/20)/1000);
RATE2(i) = B*log2(1+SINR2(i));
DeltaT2(i) = epsl/RATE2(i);
else
SINR2(i) = 0;
RATE2(i) = 0;
DeltaT2(i) = 0;
end
end
%计算E
Pbs1 = zeros(1,Nu);
Pm1 = zeros(1,Nu);
Pbs2 = zeros(1,Nu);
Pm2 = zeros(1,Nu);
for i = 1:Nu
if J(i) == 1
Pbs1(i) = alpha1*Pt1;
else
Pbs1(i) = 0;
end
if J(i) == 2
Pbs2(i) = alpha2*Pt2;
else
Pbs2(i) = 0;
end
E(i) = DeltaT1(i)*Pbs1(i) + DeltaT2(i)*Pbs2(i);
end
E2(mt) = E(UserID);
end
E3(lj,bj) = mean(E2);
end
end
[X,Y] = meshgrid(lemdas,B2set);
figure
surf(X,Y,20*log10(E3'));
xlabel('lemda');
ylabel('B');
zlabel('E(db)');