g_a_w_log.m

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

% =========================================================================
% Copyright (c)2006， 东南大学移动通信国家重点实验室
% All rights reserved.
%
% 文件名称：g_A_w_Log.m
% function描述：Log-of-Proportion Definition下的目标函数

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

function [g A] = g_A_w_Log(w, N, aimFair, aveSNR, precision)

%%% 为求梯度，将w再各个方向上按照precision散开
for gradient_num = 1:length(w)
    w_for_gradient_low(gradient_num,:) = w;
    w_for_gradient_high(gradient_num,:) = w;
    
    w_for_gradient_low(gradient_num,gradient_num) = w_for_gradient_low(gradient_num,gradient_num) - precision;     % precision为梯度求解的精度，即横坐标间隔
    w_for_gradient_high(gradient_num,gradient_num) = w_for_gradient_high(gradient_num,gradient_num) + precision;
end

%%% 求中心点处的取值
ww = [1, w];                % ww是将w向量的第一个元素固定为1之后的结果
for i = 1:N
    Th(i) = quadl(@IntegralGene, 10^(-1), 10^2, 10^(-9), [], N, ww, aveSNR, i);
end
T = 0 - sum(Th);
F = Fairness_Log(N, Th) - aimFair;

%%% 求中心点两侧点处的取值
for gradient_num = 1:length(w)
    
%     gradient_num
    
    ww_for_gradient_low = [1, w_for_gradient_low(gradient_num,:)];                % ww是将w向量的第一个元素固定为1之后的结果
    ww_for_gradient_high = [1, w_for_gradient_high(gradient_num,:)];                % ww是将w向量的第一个元素固定为1之后的结果
    
    for i = 1:N
        Th_low(i) = quadl(@IntegralGene, 10^(-1), 10^2, 10^(-6), [], N, ww_for_gradient_low, aveSNR, i);
        Th_high(i) = quadl(@IntegralGene, 10^(-1), 10^2, 10^(-6), [], N, ww_for_gradient_high, aveSNR, i);
    end
    T_low = 0 - sum(Th_low);
    F_low = Fairness_Log(N, Th_low) - aimFair;

    T_high = 0 - sum(Th_high);
    F_high = Fairness_Log(N, Th_high) - aimFair;
    
    g_gradient_temp = gradient([T_low T T_high], precision);
    A_gradient_temp = gradient([F_low F F_high], precision);
    
    g(gradient_num) = g_gradient_temp(2);
    A(gradient_num) = A_gradient_temp(2);
    
    clear ww
    clear ww_for_gradient_low
    clear ww_for_gradient_high
end