% =========================================================================% 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);endfor m=1:N3     aveSNR_set3(m) = rand(1);endsum_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);endfor 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);endfor 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);endfigurexlabel('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;