CR.rar_CR matlab_Monte Carlo_cr

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');