clc;
clear;
L=64; %每个包所含的信息比特数
M=80; %每个包所含的总比特数
R=10^4; %比特速率
W=10^6; %频谱带宽
a=5*10^(-21); %接收端噪声功率密度
const=7.75*10^(-3); %常数
d=[0.32,0.46,0.57,0.66,0.74,0.81,0.88,0.94,1]; %终端到基站的距离
h=const./(d.^3.6); %路径增益
r=fun2(h,W,R);
r=floor(r.*10)./10;
N=floor(1+(W/R)./r); %终端数
Preceive=r*a*W/((W/R)-(N-1).*r); %接收端功率
Ptransmit=Preceive./h; %发射端功率
Ptransmit=Ptransmit.*10^11; %转换为标准单位(km->m)
U=L*R*((1-exp(-0.5.*r)).^M)./M./Ptransmit; %效用函数
loglog(d,U,'o'); %效用函数随发射端到基站的距离变化曲线
hold on;
axis([0.1 1 10^3 10^7]); %限制坐标范围
title('效用函数随用户与基站之间距离(km)变化曲线'); %标题
xlabel('distance between terminal and base station(km)'); %X坐标自变量
ylabel('utility function'); %Y坐标函数
hold on;
%figure;
r1=fun1(M,W,R); %信干比 r=12.4205
r=floor(r1*10)/10; %信干比 r=12.4000
N=floor(1+(W/R)/r); %终端数
Preceive=r*a*W/((W/R)-(N-1)*r); %接收端功率
d=[0.32,0.46,0.57,0.66,0.74,0.81,0.88,0.94,1]; %终端到基站的距离
h=const./(d.^3.6); %路径增益
Ptransmit=Preceive./h; %发射端功率
Ptransmit=Ptransmit.*10^11; %转换为标准单位(km->m)
U=L*R*((1-exp(-0.5*r))^M)/M./Ptransmit; %效用函数
loglog(d,U,'+'); %效用函数随发射端到基站的距离变化曲线
hold on;
axis([0.1 1 10^3 10^7]); %限制坐标范围
legend('价格因子与路径增益有关的曲线','不含价格因子的曲线','location','northwest');
