宏小区覆盖仿真：Okumura-Hata模型与MATLAB实现

简介北京邮电大学移动通信仿真与规划课程设计实验2的宏小区覆盖配套资料面向通信工程专业本科生适用于需要完成定向天线增益仿真与宏小区覆盖实验的课程场景。压缩包共8个文件以7个MATLAB源程序为主包括主程序main.m/main1.m/main2.m以及GetBinGain、GetPathLoss、GetBestSite等关键函数分别实现天线增益计算、路径损耗建模、基站最优选址等功能覆盖宏小区覆盖仿真的完整链路另有1份“移动网络报告.docx”实验报告记录实验过程、结果分析与拓展实验内容。整个资源包仅2.92MB结构紧凑便于快速下载查看。目前已有511人学习下载。通过对照源码与报告可深入理解宏小区覆盖仿真的实现思路掌握定向天线增益仿真方法并能基于下发开源代码完成规定实验及拓展实验同时为实验报告撰写提供参考有效节省从零搭建程序与梳理文档的时间。1. 宏小区覆盖仿真这门课设代码不是黑匣子拆开全是通信课本里的公式做移动通信仿真与规划课程设计时最尴尬的事莫过于课件上把 Okumura-Hata 模型、天线方向图、链路预算讲得清清楚楚可一打开 MATLAB 就不知道怎么把公式变成能跑的脚本。这份宏小区覆盖实验资源恰好把从传播模型到站址选择的完整链路都实现了——7 个 m 文件相互调用从路径损耗计算到定向天线增益再到最佳站址搜索每一个都是通信课本里的公式直接落地。你不用从头发明轮子只需要看懂它的计算逻辑、跑通主循环再按自己的需求改参数。它适合正在做移动网络仿真课设的同学也适合刚接触网络规划、想看看基站覆盖仿真代码长什么样的开发者。2. 工程结构与仿真链路从 7 个 m 文件看覆盖评估的完整流程2.1 目录结构与文件分工每个脚本在链路里扮演的角色打开压缩包之后你会发现这份资源不是单个脚本而是一组互相调用的函数库加两个入口脚本。先看整体结构文件角色一句话说明main.m基础主入口老师下发的原始模板定义了仿真主流程main1.m基础实验脚本完成宏小区覆盖仿真的完整版本main2.m拓展实验脚本对应拓展实验通常换了场景或增加站址GetPathLoss.m传播模型计算信号从基站到接收点的路径损耗GetBinGain.m天线增益计算定向天线的二元增益单用户版本GetBinGainmultiple.m天线增益多用户/多角度版本的增益计算GetBestSite.m选址算法在候选站址集合中选出覆盖效果最好的位置移动网络报告.docx实验报告完整的实验过程、参数表、结果截图和分析这里有个容易被忽略的点为什么会有 main.m 和 main1.m、main2.m 三个入口这其实是课程设计的典型结构——main.m 往往是老师下发的原始主程序main1.m 是完成基础实验后整理的版本main2.m 则是拓展实验的入口。我在实际跑代码的时候习惯先打开 main.m 看一眼它调了哪些函数再逐个打开被调用的文件形成一张调用关系图这样比直接跑要少走很多弯路。比如你会发现 main1.m 里调用了 GetPathLoss、GetBinGain 和 GetBestSite而 main2.m 可能只换了参数区间或加了第二个基站的位置判断。2.2 仿真主流程从清空工作区到输出覆盖图的执行链路整个仿真是一条很清晰的数据流从参数定义开始到覆盖结果可视化结束。以下是把 main1.m 的流程骨架还原出来的样子% main1.m 核心流程骨架 clc; clear; close all; % 环境与仿真参数 fc 900; % 载波频率 MHz hte 30; % 基站天线高度 m hre 1.5; % 移动台天线高度 m BsPower 43; % 基站发射功率 dBm BsPos [1000, 1000]; % 基站位置坐标 m % 生成接收点网格 [xGrid, yGrid] meshgrid(0:50:2000, 0:50:2000); % 计算每个接收点与基站的距离 d sqrt((xGrid - BsPos(1)).^2 (yGrid - BsPos(2)).^2) / 1000; % 调用路径损耗函数计算传播损耗 PL GetPathLoss(fc, hte, hre, d, 1); % 计算接收点的方位角 angleGrid atan2d(yGrid - BsPos(2), xGrid - BsPos(1)); % 最终接收功率 发射功率 - 路径损耗 天线增益 Rsrp BsPower - PL GetBinGain(angleGrid, 65, 18); % 可视化 contourf(xGrid, yGrid, Rsrp); colorbar;这段代码的主干逻辑不复杂关键在于几个参数的物理含义。fc、hte、hre 这三个参数几乎每个同学都会改它们直接影响路径损耗的数值大小BsPower 决定覆盖半径的整体底噪水平。最需要关注的是最后一步的天线增益——宏小区用的是定向天线不是全向天线所以接收点的信号强度会因方位角不同而差异明显这也是这个实验比单纯算路径损耗有意思的地方。参数说明meshgrid(0:50:2000, 0:50:2000)生成了 41×41 的接收点网格步长 50 米可以理解为仿真精度改小到 10 米会得到更平滑的覆盖图但计算量会从 1681 个点暴涨到 40401 个点跑起来会慢很多。我一般建议先用 50 米步长验证逻辑最后出图再加密网格。2.3 报告与代码的对应关系写报告别脱离仿真数据压缩包里的报告文档不是摆设它和代码是逐节对应的。报告的第一部分会写实验目的和技术背景对应代码里对传播模型的选择与应用第二部分是仿真参数表对应 main 脚本开头的参数定义区第三部分是结果分析需要你改参数后重新跑出覆盖图截图放进去并解释为什么覆盖半径变大或变小。我在做这类课设时有个习惯把代码里的每个参数都标注上单位报告里设计一张参数总表把增益、损耗、覆盖半径这些中间量填进去。这样做有两个好处——答辩时任何一个数字被问到都能直接指到代码的对应行讲清楚来源改参数做对比实验时也不会出现「改了 BsPower 忘了同步报告里的发射功率」这种低级翻车。2.4 三个入口脚本的差异与选择main.m、main1.m、main2.m 不是简单的重命名关系。main.m 是基础模板大概率只覆盖了全向天线的宏小区覆盖计算main1.m 在它基础上引入了定向天线增益的完整调用链——这正是摘要里提到的「完成定向天线增益仿真」main2.m 则是拓展实验的入口常见做法是改变了站址数量或增加了不同场景的对比比如把接收点区域改成不规则多边形或者把单个基站改成两个基站的联合覆盖评估。如果你直接双击 main.m 就开跑很可能发现结果图特别朴素——没有方向性差异周围一圈信号强度完全一样。这不是代码坏了而是你打开的文件不对。先跑 main1.m再对比 main2.m 的差异才能理解这套实验想让你掌握的东西同一套代码骨架通过改入口和参数就能模拟出不同的网络规划场景。3. 核心模块拆解路径损耗、天线增益与最佳站址的计算细节3.1 路径损耗计算Okumura-Hata 模型在 MATLAB 里的落地GetPathLoss.m 是整个项目最基础也最关键的函数。它实现的是 Okumura-Hata 经验模型适用于 150 MHz 到 1500 MHz 频段、基站天线高度 30 到 200 米、移动台天线高度 1 到 10 米的宏小区场景。这也是为什么这个实验叫「宏小区覆盖」——模型本身的适用范围就决定了它是为宏基站设计的。function L GetPathLoss(fc, hte, hre, d, envType) % 计算 Okumura-Hata 路径损耗 % fc : 载波频率 (MHz) % hte : 基站天线高度 (m) % hre : 移动台天线高度 (m) % d : 基站到接收点的距离 (km) % envType : 环境类型 1-城市 2-郊区 3-开阔地 % 城市环境修正因子 a(hre) if envType 1 a (1.1 * log10(fc) - 0.7) * hre - (1.56 * log10(fc) - 0.8); L 69.55 26.16 * log10(fc) - 13.82 * log10(hte) - a ... (44.9 - 6.55 * log10(hte)) * log10(d); elseif envType 2 % 郊区修正 L_city GetPathLoss(fc, hte, hre, d, 1); L L_city - 2 * (log10(fc / 28))^2 - 5.4; elseif envType 3 % 开阔地修正 L_city GetPathLoss(fc, hte, hre, d, 1); L L_city - 4.78 * (log10(fc))^2 18.33 * log10(fc) - 40.94; else error(envType 只能取 1、2、3分别对应城市、郊区、开阔地); end end这套代码的设计值得注意两个点。第一郊区修正和开阔地修正复用了城市环境的计算结果——这是一种非常典型的「基础模型 修正项」设计你在写其他传播模型比如 COST-231时也能沿用同样的分层思路。第二所有参数都要求明确单位fc 必须是 MHz、d 必须是 km这和 Okumura-Hata 公式的适用范围严格对应。我在跑代码时犯过这样的错直接把距离矩阵送进去忘了除以 1000结果损耗值小了几十个 dB覆盖图一片全绿看起来「信号好得离谱」其实就是单位没换算。参数说明a(hre)是移动台天线高度修正因子频率不同公式不同代码里用的是 150~1500 MHz 频段的版本。如果拓展实验里把 fc 改到 2000 MHzOkumura-Hata 就不再适用这就是为什么报告里要求你说明模型适用范围。3.2 定向天线增益GetBinGain 和 GetBinGainmultiple 到底差在哪定向天线增益是这个实验的核心考点。基站不是向四面八方均匀辐射而是把能量集中到某个扇区方向这就是「定向天线」的意义。GetBinGain.m 的实现思路是输入接收点相对基站天线的方位角和天线参数输出该方向的增益值。function gain GetBinGain(angle, beamWidth, maxGain) % 计算定向天线的二元增益 % angle : 接收点相对天线主瓣方向的夹角 (度) % beamWidth: 天线半功率波束宽度 (度) % maxGain : 天线最大增益 (dBi) % 将角度归一化到 -180 ~ 180 angle mod(angle 180, 360) - 180; if abs(angle) beamWidth / 2 % 主瓣区域增益随角度偏移呈平方衰减 gain maxGain - 12 * (angle / (beamWidth / 2))^2; else % 旁瓣区域直接设为一个固定衰减值 gain maxGain - 20; end % 限制最小增益避免出现负无穷大 gain max(gain, -10); end这段代码的逻辑很直白先判断接收点是否落在主瓣范围内落在主瓣就用二次函数模拟增益随角度偏移的衰减落在旁瓣就统一给一个固定衰减。这个「主瓣二次衰减 旁瓣固定衰减」的建模方式在课程设计里足够用了它抓住了方向性天线最核心的特征朝主瓣方向的覆盖更远背离主瓣方向的信号迅速变差。GetBinGainmultiple.m 则是多用户版本区别在于输入不再是一个方位角标量而是一个方位角矩阵——对应网格上每一个接收点。我理解它的逻辑应该是把标量函数向量化不再用 if-else 逐点判断而是用逻辑索引一次性处理整块矩阵。如果你拿到代码发现 GetBinGainmultiple 的循环版本跑得很慢可以考虑用angle(abs(angle) beamWidth/2)这类逻辑索引替代循环速度能提升一个量级。参数说明beamWidth取 65 度对应常见的 65° 水平波束宽度宏基站天线maxGain取 18 dBi 也是典型宏站天线的增益参数。改这两个值会直接影响覆盖形状——把波束宽度改大到 90 度覆盖扇区变宽但前向增益会因为能量更分散而下降这就是波束宽度和增益之间的权衡关系。3.3 最佳站址选择GetBestSite 的贪心搜索逻辑GetBestSite.m 是这个项目里最接近「算法」的部分。它的任务是给定一批候选站址从中选出覆盖效果最好的一个。课程设计版本的实现通常是贪心策略——遍历所有候选位置对每个位置计算覆盖范围内的接收点数量或覆盖率取最大值对应的站址。function bestSite GetBestSite(candidates, receivers, params) % 在候选站址中搜索覆盖率最高的位置 % candidates : Nx2 矩阵每行是一个候选站址坐标 % receivers : Mx2 矩阵每行是一个接收点坐标 % params : 结构体包含 fc/hte/hre/BsPower/beamWidth/maxGain 等 bestScore -1; bestSite candidates(1, :); for i 1:size(candidates, 1) % 当前候选站址 site candidates(i, :); % 计算所有接收点到该站址的距离 (km) dists sqrt((receivers(:,1) - site(1)).^2 ... (receivers(:,2) - site(2)).^2) / 1000; % 计算路径损耗矩阵 PL GetPathLoss(params.fc, params.hte, params.hre, dists, params.envType); % 计算接收点的方位角 angles atan2d(receivers(:,2) - site(2), receivers(:,1) - site(1)); % 计算接收信号强度 Rsrp params.BsPower - PL GetBinGain(angles, params.beamWidth, params.maxGain); % 覆盖判据信号强度高于阈值的接收点占比 coverageRatio sum(Rsrp params.minRsrp) / length(Rsrp); % 更新最优解 if coverageRatio bestScore bestScore coverageRatio; bestSite site; end end end这个函数的工程思想值得拆解。它是一个标准的穷举式贪心搜索双重循环结构外层遍历候选站址内层是向量化计算整批接收点的覆盖率用覆盖率作为目标函数来评价站点优劣。minRsrp是覆盖门限一般取 -95 dBm 或 -100 dBm低于这个值的区域就算「覆盖空洞」。改这个阈值对结果影响很大——阈值设得越低覆盖率越高但选出来的站址可能离用户密集区很远因为它只需要满足弱信号也能被覆盖到。参数说明params用结构体打包传递一组仿真参数这是 MATLAB 里组织多参数函数的常见做法比写七八个形参清晰得多也方便后续在 main 脚本里统一修改。3.4 从模块到系统三个函数如何配合完成覆盖评估把三个核心函数串起来看整个仿真系统不过是一条流水线GetPathLoss 负责回答「距离远了信号衰减多少」GetBinGain 负责回答「方向偏了天线增益打折多少」GetBestSite 则把这俩的结果汇总成覆盖率指标用来回答「基站建在哪最好」。理解了这条链路你就掌握了这套代码的骨架剩下的一切都是在这条链路上加参数、改模型、换场景。4. 运行与避坑指南这些 MATLAB 报错我替你踩过了4.1 运行环境与前置条件版本差异带来的玄学问题这份代码是标准的 MATLAB 脚本和函数文件不依赖任何第三方工具箱。基础版本用 MATLAB R2016a 之后的版本都能跑因为代码里用到的 meshgrid、contourf、atan2d 都是基础函数。需要注意的是如果你用的是 GNU Octaveatan2d这个函数可能不存在需要换成atan2再手动转角度。运行前把解压后的所有 m 文件放在同一个目录下然后在 MATLAB 里cd到该目录。很多同学第一次跑报「未定义函数或变量 GetPathLoss」原因就是当前工作路径不对函数文件不在搜索路径里。这个错误处理方式很简单在 MATLAB 主页点「设置路径」把工程目录添加进去或者直接 cd 到文件夹。4.2 高频报错排查从「矩阵维度不一致」到「索引超出范围」我整理了几条最有代表性的踩坑记录每一条都对应一个真实试错场景。现象一运行 main1.m 报「矩阵维度必须一致」。原因main 脚本里手动写了Rsrp BsPower - PL Gain但 PL 来自GetPathLoss处理后的距离矩阵而 Gain 来自GetBinGain处理后的方位角矩阵两个矩阵尺寸不一致。典型的错误场景是——距离矩阵用了 41×41 的网格生成方位角矩阵却用了 40×40 的数据多出现在你手动改了网格步长却忘了同步改角度计算那一行时。 解决检查meshgrid生成的 xGrid、yGrid 维度确保所有中间变量都是同一尺寸。调试技巧是在报错行前加一句size(PL)和size(Gain)直接对比两个矩阵的形状。现象二函数GetBinGain输入参数个数不对。原因原版函数可能只接受两个参数angle 和 beamWidth你在 main 脚本里传了三个参数。或者反过来——函数定义了三个形参调用时只传了两个。 解决打开 GetBinGain.m 看函数定义行的function gain GetBinGain(...)形参个数然后对照 main 脚本里的调用语句。这是 MATLAB 最常见的函数调用错误没有之一。现象三运行 GetBestSite 时提示「索引超出数组边界」。原因候选站址矩阵candidates(1, :)不存在或者第 i 行的索引 i 超出了size(candidates, 1)。常见于候选站址是从 Excel 导入的行数比预期少。 解决在调用前加size(candidates)打印确认。此外要注意 MATLAB 数组索引从 1 开始如果代码里写成了candidates(0, :)必定报错。现象四出图后覆盖范围全是同一个颜色。原因contourf 的等值线数量太多或太少导致色阶区分度不够。另一个可能是 Rsrp 数值范围太大比如从 -120 到 -20 dBmcontourf 默认的 20 条等值线把这个范围拉伸得过于平滑。 解决手工指定等值线层次例如contourf(xGrid, yGrid, Rsrp, -120:5:-40)这样每 5 dBm 一个色阶覆盖边界看得一清二楚。4.3 参数调整的常见误用改了频率却忘了对应传播模型这份实验代码最经典的「参数误用」是把 fc 改成 2000 或 2600 MHz——这在 4G/5G 时代非常自然但 Okumura-Hata 模型本身只适用到 1500 MHz。你改了频率之后代码照跑但算出来的路径损耗已经没有理论依据了覆盖结果也失去了物理意义。我建议的做法是如果要做 2 GHz 以上的拓展实验在报告里明确写出「本实验采用 Okumura-Hata 模型该模型适用的频率上限为 1500 MHz因此 2 GHz 场景需换用 COST-231 Hata 模型」然后给出 COST-231 和 Okumura-Hata 的区别说明。诚实交代模型的适用范围比悄无声息地跑一个「看似合理但其实是错误公式」的结果要加分得多。同理接收点网格的粒度步长、基站高度超过 200 米等情况都需要在报告里给出模型适用性说明。提示拿到代码先跑通原始版本再改参数。每改一个参数出一次对比图这样出了问题能快速定位是哪个改动引起的。5. 二次开发思路把这套仿真骨架变成你自己的工具5.1 从单站到多站把 GetBestSite 扩展成联合选址现在的 GetBestSite 一次只选一个站。但实际宏小区网络规划中往往是多站协同覆盖——比如给出 5 个候选站址要选出 3 个使整体覆盖率最高。这是一个组合优化问题课程设计级别不需要上启发式算法一个简单的贪心迭代就够第一次调用 GetBestSite 选出最优站址把它从候选列表中移除然后让第二个站的覆盖率评估排除已覆盖的接收点再来一轮。这样每次取一个直到达到目标基站数量。我在做这类扩展时会在 GetBestSite 的返回值里增加一个逻辑向量coveredFlag标记哪些接收点已被前一个站覆盖。后续站点评估时只对未覆盖的接收点计算覆盖率这样多个站不会扎堆选在同一个位置。5.2 换一个传播模型把 GetPathLoss 的公式替换成 COST-231如果你想让项目直接对标 2 GHz 以上的 LTE 宏小区场景改造入口就是 GetPathLoss.m。COST-231 Hata 模型和 Okumura-Hata 的区别主要在一个修正项 C——接收点在大城市中心时 C 取 3 dB中小城市取 0 dB。函数结构几乎不用变只需要在内部增加一个分支if fc 1500 C 3; % 大城市中心修正 L 46.3 33.9*log10(fc) - 13.82*log10(hte) - a ... (44.9 - 6.55*log10(hte))*log10(d) C; else % 原 Okumura-Hata 逻辑 end这样改完之后你就可以在同一套代码里对比 900 MHz 和 1800 MHz 的覆盖差异报告里多一张对比图整个实验的深度立刻不一样。5.3 从仿真到告警把覆盖率结果输出成评估报告最后一个改造点是把仿真结果落到数据层面。现有代码出图后没有把统计指标写出来你可以加一段后处理统计覆盖率、平均信号强度、最差 5% 接收点的信号值用fprintf在命令行输出或者写入 CSV 文件。这样改参数做对比实验时不需要肉眼看颜色图估数字直接看输出行的变化就够了。从那以后我做这类仿真课设都会强制自己走一遍「跑通原版 → 记录基线数据 → 改一个参数 → 对比输出」的流程。参数改动记录在报告里代码版本在文件名上标号这样答辩被问到时你不会说出「我忘了怎么改的」这种话。希望这份解析能帮你在自己的课设里少走几步弯路。本文还有配套的精品资源点击获取