% Tdscdmamapini.m: 設(shè)置基站位置和小區(qū)結(jié)構(gòu)，確定計(jì)算區(qū)域
% 
% 畢海
% 2005.1.29
% Version 1
% 
if ((pathlossModel == 7) | (pathlossModel == 8)) % The map generated from nps/x 
   load(linklossFile);                           % loads the linkloss-data
   resolution = lossData(1).resolution;          % meters
   area = lossData(1).area;
   if (pathlossModel == 8)
      load(vectFile);                            % loads the building vector info
   else
      vectMap = [];
   end
else
   resolution = 50;
   area = (WraparoundCircle*2+1)/2*BSspace*[-1 1 -1 1]; % meters
   lossData = [];
   vectMap = [];
end

xmin = -((WraparoundCircle - 1)*1.5 + 1)*sqrt(3)/3*BSspace;
xmax = ((WraparoundCircle - 1)*1.5 + 1)*sqrt(3)/3*BSspace;
ymin = -1*(WraparoundCircle - 0.5)*BSspace - rem((WraparoundCircle-1),2)*BSspace/2;
ymax = WraparoundCircle*BSspace - rem((WraparoundCircle-1),2)*BSspace/2;

numBSs = 1;
for i = 1:(WraparoundCircle*2 - 1)
    for j = 1:(WraparoundCircle*2 - 1)
        BSposition(numBSs,1) = (i - 1)*sqrt(3)*BSspace/2 - (WraparoundCircle - 1)*sqrt(3)*BSspace/2;
        BSposition(numBSs,2) = ((j - 1)*2 + rem((i-1),2))*BSspace/2 - (WraparoundCircle - 1)*BSspace - rem((WraparoundCircle-1),2)*BSspace/2;
        numBSs = numBSs + 1;
    end;
end;

scatter(BSposition(:,1),BSposition(:,2));
axis([xmin xmax ymin ymax]);

yPixels = ceil((abs(ymax) - resolution/2)/resolution) + ceil((abs(ymin) - resolution/2)/resolution) + 1; 
xPixels = ceil((abs(xmax) - resolution/2)/resolution) + ceil((abs(xmin) - resolution/2)/resolution) + 1; 

[xx yy] = meshgrid(1:xPixels, 1:yPixels);
xx = xx*resolution + floor(xmin/resolution)*resolution;
yy = (yy-1)*resolution + floor(ymin/resolution)*resolution;