✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。
🍎个人主页:Matlab科研工作室
🍊个人信条:格物致知。
⛄ 内容介绍
利用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;
end
end
% 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]);
⛄ 运行结果
⛄ 参考文献
[1] 陈岚,万国春,冯志彪.基于Matlab的移动通信中多径衰落信道的仿真[J].江西科技师范学院学报, 2004(5):78-81.DOI:10.3969/j.issn.1007-3558.2004.05.023.
[2] 张建昌.在蜂窝网中应用D2D通信技术研究[D].电子科技大学[2023-07-09].DOI:10.7666/d.D499416.