比例公平调度matlab源程序

比例公平调度matlab源程序，按照Goldsmith的无线通信上的算法编程。也可以嵌入到OFDMA资源管理的算法中。 该程序有详细注释。程序函数名称为PFS1_08aug18.m

 function PFS1_meng08aug18（）
clear all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%　　　　　　　　　　　　　记录仿真系统文档                                 %
%比例公平调度
%neraspace revised on Aug 182008
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fpname = strcat（mfilename ‘.txt‘）;
fp = fopen（fpname‘a+‘）;
fprintf（fp‘\n\n*****************************************************************************\n‘）;
fprintf（fp‘+ +Proportional Fair SchedulingThe Simulated Date: %12s  + + + +\n‘datestr（now））;
fprintf（fp‘+ +             Note: SISO  Rayleigh fading channels      　　　         + + + +\n‘）;
fprintf（fp‘*****************************************************************************\n‘）;
fprintf（fp‘*-------------------------设置系统参数-------------------------------*\n‘）;
Tx = 1;

Rx = 1;

fprintf（fp‘\n‘）;

fprintf（fp‘*-------------------------设置仿真参数-------------------------------*\n‘）;
Simtime = 10000;
fprintf（fp ‘The Simulated time slot Number = %d \n‘Simtime）;
NumUser=10;  

EsN0dB_Aveage=0;                        
EsN0_Aveage=10^（EsN0dB_Aveage/10）;     
Num_group=2;                            
Num_eachgroup=ceil（NumUser/Num_group）; 
EsN0dB_unsymmetrical=[-5 +5]*1;
SNRgap=ceil（max（EsN0dB_unsymmetrical）-min（EsN0dB_unsymmetrical））;
EsN0_unsymmetrical = 10.^（EsN0dB_unsymmetrical/10）;
if SNRgap==0
    fprintf（