g_yi.m

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

% =========================================================================
% Copyright (c)2006， 东南大学移动通信国家重点实验室
% All rights reserved.
%
% 文件名称：G_yi.m
% function描述：供double_int.m调用的子程序，用于计算可变上限多重积分

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

function f = G_yi(SNR, x1, x2, K, vecN, aveSNR, aveSNRn, i, j)
SNR = SNR(:);
length_SNR = length(SNR);
if isstr(x1)==1;xx1=feval(x1,SNR);else xx1=x1*ones(size(SNR)); end
if isstr(x2)==1;xx2=feval(x2,SNR);else xx2=x2*ones(size(SNR)); end

for p=1:length_SNR
    f(p) = log2 (1 + SNR(p)) .* exp(- SNR(p) ./ aveSNR(i,j)) ./ aveSNR(i,j);
    for k= 1:vecN(i)
      if k == j
         f(p) = f(p);
      else
         f(p) = f(p) .* (1 - exp(- SNR(p) ./ aveSNR(i,k)));
      end
    end
    
%%%自己加的测试语句
%f= log2 (1 + SNR) .* exp(- SNR ./ aveSNR(i,j)) ./ aveSNR(i,j);
%    for k= 1:vecN(i)
%      if k == j
%         f = f;
%      else
%         f = f .* (1 - exp(- SNR ./ aveSNR(i,k)));
%      end
     for n = 1:K
         if n == i
           f(p) = f(p);
         else
            f(p)=f(p) .* quadl(@sub_gene,xx1(p),xx2(p),10^(-12),[],SNR(p),n, vecN, aveSNR, aveSNRn, i, j);
         end
     end
end
f=f(:);