tdscdmabestcell.m

TD-SCDMA移动通信系统的系统级仿真平台

% tdscdmaBestCell.m; 计算导频覆盖功率和导频覆盖概率，确定每一栅格点的最佳服务小区
% % 毕海
% 2005.1.29
% Version 1
% 
bar1 = waitbar(0,'计算每个栅格点的导频覆盖功率');
PilotPower = zeros(size(linkloss));
maxBSOneAntennaPerCodePower = maxBSOneAntennaPower - 10*log10(16);
maxPccPCHPowerOneAntenna = maxBSOneAntennaPerCodePower + 3;
maxPccPCHPowerTotalAntenna = maxPccPCHPowerOneAntenna + 10*log10(numAntenna);
PilotPower = (maxPccPCHPowerTotalAntenna + AntennaGain) - linkloss;

for i = 1:xPixels
    waitbar(i/xPixels,bar1);
    for j = 1:yPixels
        BestCellPower(j,i) = max(PilotPower(:,j,i));
        if BestCellPower(j,i) < -110
            BestCellPower(j,i) = -200;
        end;
        for k = 1:numBSs
            %temp = squeeze(PilotPower(k,j,i));
            if BestCellPower(j,i) == PilotPower(k,j,i)
                BestCellid(j,i) = k;
            elseif BestCellPower(j,i) == -200
                BestCellid(j,i) = -200;    
            end;
        end;
    end;
end;
close(bar1);
CellColor = randperm(numBSs);
bar1 = waitbar(0,'计算每个栅格点的最佳服务小区');
for i = 1:xPixels
    waitbar(i/xPixels,bar1);
    for j = 1:yPixels
        if BestCellid(j,i) == -200
            BestCellColor(j,i) = 0;
        else
            BestCellColor(j,i) = CellColor(BestCellid(j,i));
        end;
    end;
end;
close(bar1);
fig1 = surf(xx1,yy1,BestCellPower);title('导频覆盖'),view([0 0 1]),shading interp,axis([xmin xmax ymin ymax]);
hold,scatter3(BSposition(:,1),BSposition(:,2),BSposition(:,3)-200);
figure;
fig2 = surf(xx1,yy1,BestCellColor);title('最佳服务小区'),view([0 0 1]),shading interp,axis([xmin xmax ymin ymax]);
hold,scatter3(BSposition(:,1),BSposition(:,2),BSposition(:,3)-100);

%pause;
clear BestCellColor BestCellPower BestCellid PilotPower;