% =========================================================================
% Copyright (c)2006， 東南大學移動通信國家重點實驗室
% All rights reserved.
%
% 文件名稱：PFFair_log.m
% function描述：PF調度算法，在Log-of-Proportion Definition下的公平性計算

% 當前版本： 1.0
% 作    者：梁浩
% 完成日期： 2006年10月
%==========================================================================

function pfFair = PFFair_log(N, aveSNR)
for i = 1:N
    Th(i) = quadl(@PFIntegralGene_log, 10^(-3), 10^3, 10^(-12), [], N, aveSNR, i);
end

sum_all = 0;
for k = 1:N
    sum_all = sum_all + Th(k);
end

for i = 1:N
     p(i) = Th(i) / sum_all;
end

entropy = 0;
for i = 1:N
    entropy = entropy - p(i) * log2( p(i) );
end

pfFair = entropy / log2( N );