vmean2tpandsir.m

fading 在无线中的应用

function [pul, sirul] = vmean2tpandsir(vmean, user, t)
par = setpar;       % default setting

par.km = 1;        % km^2+lm^2+km*lm => the number of sites
par.lm = 1;        % related to km above
par.sps = 1;       % Number of sectors per site
par.ncluster = 3;  % number of clusters equals ncluster^2
par.kpc = 1;       % number of allocated channels per cel
par.raa = 0;       % lognormal correlation down link (typical 0.5)
par.offtraf = 10;  % average number of offered calls to a cell     
par.sirmin = -inf; % C/I level under which a call is dropped  
par.maxiter = 1;   % number of iterations in the main loo
par.mht = 0;       % mean holding time. In this case snapshot
par.vmean = vmean; % varege speed
par.dt = diff(t(1:2)); % time step between iterations
[res1,par,sys]  = runef(par); % make an initial traffic situation

par.mht = inf;     % keep all ongoing calls alive forever.
pmax = 3;         % max power per user  [dBm]
sirtarget = 10;     % target C/I for each user
for iter = 1:length(t)
 [res1, par]  = runef(par, res1{1}, sys);
 % power control 
 res1{1}.pul = min(pmax,res1{1}.pul+(sirtarget-res1{1}.sirul));
 % collect result
 sirul(iter,1,:) = flatten(res1{end}.sirul,3); 
 pul(iter,1,:)   = flatten(res1{end}.pul  ,3); 
end