optimum_search_sqp_12_users_jain.m

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

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

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

clear
clc



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 初始化
% 设置系统参数
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);

% 设置用户数目
N = 10;

% 计算用户位置
for i = 1:N
    d(i) = 2 ./ (N+1) .* i;
end

% 计算用户接收SNR的均值和标准差
mean = - 10 .* alpha .* log10(d) - L0 + mean_m + 10 .* log10(Pt./N0);
sigma = ones(1,N) .* sqrt(sigma_original.^2 + var_m);

% '8888888888888888888888'
mean
sigma

% 计算信噪比平均值的dB值
aveSNR_dB = 10 .* log10(exp(mean./beta + sigma.^2./2./beta.^2))

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Jain最优曲线搜索
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_jain
end



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Log最优曲线搜索
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_log
end



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Gini最优曲线搜索
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_gini
end