priceutility.asv

来自「随着蜂窝移动通信市场的巨大发展」· ASV 代码 · 共 41 行

clear;
clc;

L=64;                                                                   %每个包所含的信息比特数
M=80;                                                                   %每个包所含的总比特数
R=10^4;                                                                 %码片数率
W=10^6;                                                                 %比特数率
const=7.75*10^(-3);                                                     %常数
a=5*10^(-21);                                                           %接收端噪声功率密度
r1=fun1(M,W,R);                                                         %信干比（小数点后两位）
r=floor(r1*10)/10;                                                      %信干比（小数点后一位）
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;
U=L*R*((1-exp(-0.5*r))^M)/M./Ptransmit;                                 %效用函数

loglog(d,U,'r+-');                                                      %效用函数随发射端到基站的距离变化曲线
hold on
axis([0.1 1 10^3 10^7]);

%含代价因子
r=11.7;                                                                 %含价格因子的信干比
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;
U=L*R*((1-exp(-0.5*r))^M)/M./Ptransmit;                                 %含价格效用函数

loglog(d,U,'ro-');
axis([0.1 1 10^3 10^7]);

xlabel('distance between terminal and base station(km)');
ylabel('Utility (b/J)');
text(10^-0.4,10^4.8,'Price=0');
text(10^-0.3,10^6.2,'Price=50p');