search_set3.m

TDM系统中的最优折中曲线搜索及调度算法性能评估

% =========================================================================
% Copyright (c)2006， 东南大学移动通信国家重点实验室
% All rights reserved.
%
% 文件名称：search_set3.m
% function描述：针对第3组用户的最优折中曲线搜索

% 当前版本： 1.0
% 作    者：梁浩
% 完成日期： 2006年10月
%==========================================================================

rand('state',0);
N3 = 36;
for m=1:1400
     test_rand = rand(1);
end
for m=1:N3
     aveSNR_set3(m) = rand(1);
end
sum_set3 = sum(aveSNR_set3);
aveSNR_set3 = aveSNR_set3 / sum_set3 * 360;

 PFFair_set3_jain = PFFair_jain(N3,aveSNR_set3);
 PFTh_set3_jain = PFTh_jain(N3,aveSNR_set3);
 maxCIFair_set3_jain = MaxCIFair_jain(N3, aveSNR_set3);
 maxCITh_set3_jain = MaxCITh_jain(N3, aveSNR_set3);
 
  PFFair_set3_gini = PFFair_gini(N3,aveSNR_set3);
 PFTh_set3_gini = PFTh_gini(N3,aveSNR_set3);
 maxCIFair_set3_gini= MaxCIFair_gini(N3, aveSNR_set3);
 maxCITh_set3_gini = MaxCITh_gini(N3, aveSNR_set3);
 
 PFFair_set3_log = PFFair_log(N3,aveSNR_set3);
 PFTh_set3_log = PFTh_log(N3,aveSNR_set3);
 maxCIFair_set3_log= MaxCIFair_log(N3, aveSNR_set3);
 maxCITh_set3_log = MaxCITh_log(N3, aveSNR_set3);
 
for i=1:11
    aimFair_set3_jain(i)=maxCIFair_set3_jain + (PFFair_set3_jain - maxCIFair_set3_jain) ./ 10 .* (i-1);
    [opW_set3_jain(i,:), negopTh_set3_jain(i)] = OT_jain(N3, aimFair_set3_jain(i), aveSNR_set3)
    opTh_set3_jain(i) = - negopTh_set3_jain(i);
end

for i=1:11
    aimFair_set3_log(i)=maxCIFair_set3_log + (PFFair_set3_log - maxCIFair_set3_log) ./ 10 .* (i-1);
    [opW_set3_log(i,:), negopTh_set3_log(i)] = OT_log(N3, aimFair_set3_log(i), aveSNR_set3)
     opTh_set3_log(i) = - negopTh_set3_log(i);
end

for i=1:11
    aimFair_set3_gini(i)=PFFair_set3_gini + (maxCIFair_set3_gini - PFFair_set3_gini) ./ 10 .* (i-1);
    [opW_set3_gini(i,:), negopTh_set3_gini(i)] = OT_gini(N3, aimFair_set3_gini(i), aveSNR_set3)
     opTh_set3_gini(i) = - negopTh_set3_gini(i);
end

figure
xlabel('Fairness');
ylabel('Throughput (bps/Hz)');
plot(aimFair_set3_jain, opTh_set3_jain, 'k^-');
hold on;
plot(aimFair_set3_gini, opTh_set3_gini, 'ko-');
hold on;
plot(aimFair_set3_log, opTh_set3_log, 'ks-');
hold on;