?function?PFS1_meng08aug18（）
clear?all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%　　　　　　　　　　　　　記錄仿真系統(tǒng)文檔?????????????????????????????????%
%比例公平調(diào)度
%neraspace?revised?on?Aug?182008
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fpname?=?strcat（mfilename?‘.txt‘）;
fp?=?fopen（fpname‘a(chǎn)+‘）;
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‘*-------------------------設(shè)置系統(tǒng)參數(shù)-------------------------------*\n‘）;
Tx?=?1;

Rx?=?1;

fprintf（fp‘\n‘）;

fprintf（fp‘*-------------------------設(shè)置仿真參數(shù)-------------------------------*\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（