optimum_search_6_users.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 = 6;

% 计算用户位置
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);
% % % % %     [opW_jain(i,:), negopTh_jain(i)] = OT_jain(N, aimFair_jain(i), mean, sigma)
% % % % %     opTh_jain(i) = - negopTh_jain(i);
% % % % % 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);
    [opW_log(i,:), negopTh_log(i)] = OT_log(N, aimFair_log(i), mean, sigma)
    opTh_log(i) = - negopTh_log(i);
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);
    [opW_gini(i,:), negopTh_gini(i)] = OT_gini(N, aimFair_gini(i), mean, sigma)
    opTh_gini(i) = - negopTh_gini(i);
end