源程序：蜂窝小区经典模型，7小区3扇区

clf;
d = 0:60:360;
ISD = 500;
side = ISD/sqrt(3);
side_sector = ISD/3;

offset_bs_7 = [
    0,0;
    ISD*sqrt(3)/2,ISD/2;
    0,ISD;
    -ISD*sqrt(3)/2,ISD/2;
    -ISD*sqrt(3)/2,-ISD/2;
    0,-ISD;
    ISD*sqrt(3)/2,-ISD/2;
];

ds = 30:60:390;

offset_sector = [
    
    sqrt(3)*side_sector/2, side_sector/2;
    -sqrt(3)*side_sector/2, side_sector/2;
    0, -side_sector;
    
];

for i = 1:7
    x = side * cosd(d) + offset_bs_7(i,1);
    y = side * sind(d) + offset_bs_7(i,2);
       
    plot(x,y);hold on;
    text(offset_bs_7(i,1),offset_bs_7(i,2),num2str(i-1));
    
    for j = 1:3
        xs = side_s * cosd(ds) + offset_sector(j,1) + offset_bs_7(i,1);
        ys = side_s * sind(ds) + offset_sector(j,2) + offset_bs_7(i,2);
        plot(xs,ys);
        text(xs(1)-sqrt(3)*side_sector/2,ys(1)-side_sector/2,num2str((i-1)*3+j-1))
    end
end

offset_bs_19 = [ 
    ISD*sqrt(3),ISD;
    ISD*sqrt(3)/2,1.5*ISD;
    0,2*ISD;
    -ISD*sqrt(3)/2,1.5*ISD;
    -ISD*sqrt(3),ISD;
    -ISD*sqrt(3),0;
    -ISD*sqrt(3),-ISD;
    -ISD*sqrt(3)/2,-1.5*ISD;
    0,-2*ISD;
    ISD*sqrt(3)/2,-1.5*ISD;
    ISD*sqrt(3),-ISD;
    ISD*sqrt(3),0;
];
for i = 8:19
    x = side * cosd(d) + offset_bs_19(i-7,1);
    y = side * sind(d) + offset_bs_19(i-7,2);
    plot(x,y);hold on;
    text(offset_bs_19(i-7,1),offset_bs_19(i-7,2),num2str(i-1));
end


axis equal