% =========================================================================
% Copyright (c)2006， 東南大學移動通信國家重點實驗室
% All rights reserved.
%
% 文件名稱：HyNAFair.m
% function描述：基于Hybrid N/A調度算法的公平性計算，Jain Index

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

function hyNAFair = HyNAFair(K, vecN, aveSNR)
for i = 1:K
    for j= 1:vecN(i)
        Th(i,j) = quadl(@HyNAIntegralGene, 10^(-3), 10^3, 10^(-12), [], K, vecN, aveSNR, i, j);
    end
end
hyNAFair = (sum(sum(Th))) .^ 2 ./ (sum(vecN) .* sum(sum(Th .^ 2)));