基于matlab模拟蜂窝网络

%% 初始化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;            y1(:, p + s) = y0 + j * sqrt(3) * r;            p = p + 1;        end    else % mod(layer_count, 2) == 0 区群层数为偶数，最外列为偶数个小区        p = 0; % 区群列内小区编号计数        for j = [- fix(layer_count / 2) : -1 , 1 : fix(layer_count / 2)]            if j > 0                t = 0.5;            else                t = -0.5;            end            bs_x(c) = i * 1.5 * r; % 小区中心横坐标，即横坐标偏移量            bs_y(c) = (j - t) * sqrt(3) * r; % 小区中心纵坐标，即纵坐标偏移量            c = c + 1;            x1(:, p + s) = x0 + i * 1.5 * r;            y1(:, p + s) = y0 + (j - t) * sqrt(3) * r;            p = p + 1;        end    end    s = s + layer_count;    if i < 0 % 列未达到中心列前，列内小区数逐列递增        layer_count = layer_count + 1;    else % i >= 0 列达到中心列后，列内小区数逐列递减        layer_count = layer_count - 1;    endend% bs_x = sum(x1, 1) / 6; % 各个蜂窝小区中心点横坐标（误）% bs_y = sum(y1, 1) / 6; % 各个蜂窝小区中心点纵坐标（误）figure;plot(x1, y1, 'k'); % 绘制蜂窝区群图axis equal;hold on;set(gcf,'Color', [1 1 1], 'Units', 'Normalized', 'Position', [0.25 0.25 0.5 0.5]);