% =========================================================================% Copyright (c)2006， 東南大學(xué)移動(dòng)通信國(guó)家重點(diǎn)實(shí)驗(yàn)室% All rights reserved.%% 文件名稱：Optimum_Search_6_users.m% function描述：針對(duì)6用戶場(chǎng)景的最優(yōu)折中曲線搜索% 當(dāng)前版本： 1.0% 作    者：梁浩% 完成日期： 2007年6月%==========================================================================clearclc%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 初始化% 設(shè)置系統(tǒng)參數(shù)alpha = 3.52249;L0 = 140.4193;sigma_original = 8;mean_m = -2.5;var_m = 31;Pt = 10.^(40./10)./1000;N0 = 10.^(-100./10)./1000;beta = 10 ./ log(10);% 設(shè)置用戶數(shù)目N = 18;% 計(jì)算用戶位置for i = 1:N    d(i) = 2 ./ (N+1) .* i;end% 計(jì)算用戶接收SNR的均值和標(biāo)準(zhǔn)差mean = - 10 .* alpha .* log10(d) - L0 + mean_m + 10 .* log10(Pt./N0);sigma = ones(1,N) .* sqrt(sigma_original.^2 + var_m);% '8888888888888888888888'meansigma% 計(jì)算信噪比平均值的dB值aveSNR_dB = 10 .* log10(exp(mean./beta + sigma.^2./2./beta.^2))%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Jain最優(yōu)曲線搜索maxCIFair_jain = MaxCIFair_jain(N, mean, sigma)maxCITh_jain = MaxCITh_jain(N, mean, sigma)pfFair_jain = PFFair_jain(N, mean, sigma, aveSNR_dB)pfTh_jain = PFTh_jain(N, mean, sigma, aveSNR_dB)for i=1:11    aimFair_jain(i)=maxCIFair_jain + (pfFair_jain - maxCIFair_jain) ./ 10 .* (i-1);    [op_th_jain(i) op_fair_jain(i) op_w_jain(i,:)] = SQP_Optimization_Jain(N, aimFair_jain(i), mean, sigma);    'Jain 8888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888'    op_th_jain     op_fair_jain    op_w_jainend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Log最優(yōu)曲線搜索maxCIFair_log = MaxCIFair_log(N, mean, sigma)maxCITh_log = MaxCITh_log(N, mean, sigma)pfFair_log = PFFair_log(N, mean, sigma, aveSNR_dB)pfTh_log = PFTh_log(N, mean, sigma, aveSNR_dB)for i=1:11    aimFair_log(i)=maxCIFair_log + (pfFair_log - maxCIFair_log) ./ 10 .* (i-1);    [op_th_log(i) op_fair_log(i) op_w_log(i,:)] = SQP_Optimization_Log(N, aimFair_log(i), mean, sigma);    'Log 8888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888'    op_th_log     op_fair_log    op_w_logend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Gini最優(yōu)曲線搜索maxCIFair_gini = MaxCIFair_gini(N, mean, sigma)maxCITh_gini = MaxCITh_gini(N, mean, sigma)pfFair_gini = PFFair_gini(N, mean, sigma, aveSNR_dB)pfTh_gini = PFTh_gini(N, mean, sigma, aveSNR_dB)for i=1:11    aimFair_gini(i)=maxCIFair_gini + (pfFair_gini - maxCIFair_gini) ./ 10 .* (i-1);    [op_th_gini(i) op_fair_gini(i) op_w_gini(i,:)] = SQP_Optimization_Gini(N, aimFair_gini(i), mean, sigma);    'Gini 8888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888888'    op_th_gini     op_fair_gini    op_w_giniend