% =========================================================================% Copyright (c)2006， 東南大學移動通信國家重點實驗室% All rights reserved.%% 文件名稱：evalution_Gini_12.m% function描述：調度算法性能評估，12用戶，Gini Index% 當前版本： 1.0% 作    者：梁浩% 完成日期： 2006年10月%==========================================================================clear;clc;%隨機集合的再現、歸一化、從小到大排序rand('state',0);N1 = 12;for m=1:560    test_rand = rand(1);end for m=1:N1    aveSNR(m) = rand(1);endsum_set = sum(aveSNR);aveSNR = aveSNR / sum_set * 120 ;%從小到大排列aveSNR_sorted = sort(aveSNR);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 分組#1(折中分組，多取一個轉置)，Hybird N/A K=2N = 12;clear aveSNR_grouped;K = 2;uses_per_group = N ./ (12/K);vecN = uses_per_group .* ones(1,(12/K));for k = 1:12/K    aveSNR_grouped(k,:) = aveSNR_sorted( [( (k-1) * uses_per_group + 1) : (k * uses_per_group)] );endaveSNR_grouped = aveSNR_grouped';clear uses_per_group;clear vecN;uses_per_group = N ./ K;vecN = uses_per_group .* ones(1,K);%計算Hybrid N/A方法的Throughput、FairnesshyNAFair12 = HyNAFair_Gini(K, vecN, aveSNR_grouped)hyNATh12 = HyNATh(K, vecN, aveSNR_grouped)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 分組#2(折中分組，多取一個轉置)，Hybird N/A K=3N = 12;clear aveSNR_grouped;K = 3;uses_per_group = N ./ (12/K);vecN = uses_per_group .* ones(1,(12/K));for k = 1:12/K    aveSNR_grouped(k,:) = aveSNR_sorted( [( (k-1) * uses_per_group + 1) : (k * uses_per_group)] );endaveSNR_grouped = aveSNR_grouped';clear uses_per_group;clear vecN;uses_per_group = N ./ K;vecN = uses_per_group .* ones(1,K);%計算Hybrid N/A方法的Throughput、FairnesshyNAFair2 = HyNAFair_Gini(K, vecN, aveSNR_grouped)hyNATh2 = HyNATh(K, vecN, aveSNR_grouped)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 分組#3(折中分組，多取一個轉置)，Hybird A/N K=4N = 12;clear aveSNR_grouped;K = 4;uses_per_group = N ./ (12/K);vecN = uses_per_group .* ones(1,(12/K));for k = 1:12/K    aveSNR_grouped(k,:) = aveSNR_sorted( [( (k-1) * uses_per_group + 1) : (k * uses_per_group)] );endaveSNR_grouped = aveSNR_grouped';clear uses_per_group;clear vecN;uses_per_group = N ./ K;vecN = uses_per_group .* ones(1,K);%計算Hybrid A/N方法中用到的中間參數aveSNRnfor n = 1:K    aveSNRn(n) = quadl(@HyANmodpdf, 10^(-3), 10^3, 10^(-12), [], vecN, aveSNR_grouped, n);end%計算Hybrid A/N方法的Throughput、Fairnessfor i = 1:K    for j= 1:vecN(i)        Th12(i,j) = double_int(10^(-4),'x_hi',10^(-3) ,10^3, K, vecN, aveSNR_grouped, aveSNRn, i, j)    endendhyANFair12 = Fairness_Gini(N, Th12)hyANTh12 = sum(sum(Th12))%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 分組#4(折中分組，多取一個轉置)，Hybird A/N K=6N = 12;clear aveSNR_grouped;K = 6;uses_per_group = N ./ (12/K);vecN = uses_per_group .* ones(1,(12/K));for k = 1:12/K    aveSNR_grouped(k,:) = aveSNR_sorted( [( (k-1) * uses_per_group + 1) : (k * uses_per_group)] );endaveSNR_grouped = aveSNR_grouped';clear uses_per_group;clear vecN;uses_per_group = N ./ K;vecN = uses_per_group .* ones(1,K);%計算Hybrid A/N方法中用到的中間參數aveSNRnfor n = 1:K    aveSNRn(n) = quadl(@HyANmodpdf, 10^(-3), 10^3, 10^(-12), [], vecN, aveSNR_grouped, n);end%計算Hybrid A/N方法的Throughput、Fairnessfor i = 1:K    for j= 1:vecN(i)        Th2(i,j) = double_int(10^(-4),'x_hi',10^(-3) ,10^3, K, vecN, aveSNR_grouped, aveSNRn, i, j)    endendhyANFair2 = Fairness_Gini(N, Th2)hyANTh2 = sum(sum(Th2))%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 其它幾種調度算法：max C/I、PF、M-LWDFmaxCIFair= MaxCIFair_Gini(N, aveSNR)maxCITh = MaxCITh(N, aveSNR)pfFair = PFFair_Gini(N,aveSNR)pfTh = PFTh(N,aveSNR)rrFair = RRFair_Gini(N,aveSNR)rrTh = RRTh(N,aveSNR)[PBFS_Th_100, PBFS_Fair_100] = PBFS_Gini(N, aveSNR, ones(1, N), 1000, 50)[PBFS_Th_500, PBFS_Fair_500] = PBFS_Gini(N, aveSNR, ones(1, N), 1000, 250)[PBFS_Th_3000, PBFS_Fair_3000] = PBFS_Gini(N, aveSNR, ones(1, N), 1000, 1500)% [S_PBFS_Th, S_PBFS_Fair] = S_PBFS_Gini(N, aveSNR, ones(1, N), 1000)% load mat_precise_12_24_36.matfigureplot(aimFair_set1_gini,opTh_set1_gini,'k.-')hold on;plot(rrFair,rrTh,'kh')hold on;plot(maxCIFair,maxCITh , 'ks');hold on;plot(pfFair,pfTh , 'kd');hold on;plot(hyNAFair12, hyNATh12, 'k<')hold on;plot(hyNAFair2, hyNATh2, 'k>')hold on;plot(hyANFair12, hyANTh12, 'k^')hold on;plot(hyANFair2, hyANTh2, 'kv')hold on;plot(PBFS_Fair_100,PBFS_Th_100, 'k*')hold on;plot(PBFS_Fair_500,PBFS_Th_500, 'k+')hold on;plot(PBFS_Fair_3000,PBFS_Th_3000, 'kp')hold on;% % % plot(S_PBFS_Fair,S_PBFS_Th, 'k*')% % % hold on;grid on;legend('Optimal Curve','RR','Max C/I','PF','Hybrid N/A K=2','Hybrid N/A K=3','Hybrid A/N K=4','Hybrid A/N K=6', 'PBFS D=100', 'PBFS D=500', 'PBFS D=3000', 4)xlabel('Fairness');ylabel('Throughput (bps/Hz)');