%%?初始化

clc;
clear;

%%?生成一個基準蜂窩

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%總體思路：首先確定一個基準中心蜂窩小區，之后從區群最左一列開始逐列繪制整個區群。

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

r?=?1;?%?半徑
n?=?6;?%?六邊形蜂窩6個角
a?=?0?:?2?*?pi?/?n?:?2?*?pi;?%?基準蜂窩6個角定位參數
x0?=?r?*?cos（a）;?%?基準蜂窩橫坐標
y0?=?r?*?sin（a）;?%?基準蜂窩縱坐標

%%?生成一個區群

layer?=?6;
cell_num?=?1?+?3?*?layer?*?（layer?-?1）;?%?整個網絡蜂窩總數每層增加6（n-1）個蜂窩小區
x1?=?zeros（n?+?1?cell_num）;?%?存儲小區橫坐標
y1?=?zeros（n?+?1?cell_num）;?%?存儲小區縱坐標
m?=?layer?-?1;?%?區群列數計數
layer_count?=?layer;?%?區群列內小區個數計數
s?=?1;?%?區群換列調整參數
bs_x?=?zeros（1?cell_num）;?%?小區中心橫坐標
bs_y?=?zeros（1?cell_num）;?%?小區中心縱坐標
c?=?1;?%?小區中心坐標計數

for?i?=?-m?:?m
????if?mod（layer_count?2）?==?1?%?區群層數為奇數，最外列為奇數個小區
????????p?=?0;?%?區群列內小區編號計數
????????for?j?=?-?fix（layer_count?/?2）?:?fix（layer_count?/?2）
????????????bs_x（c）?=?i?*?1.5?*?r;?%?小區中心橫坐標，即橫坐標偏移量
????????????bs_y（c）?=?j?*?sqrt（3）?*?r;?%?小區中心縱坐標，即縱坐標偏移量
????????????c?=?c?+?1;
????????????x1（:?p?+?s）?=?x0?+?i?*?1.5?*?r;