简介：面向5G通信与MATLAB仿真学习者的课程资源包，内容紧扣路径损耗建模与天线辐射方向图分析，覆盖Friis传输方程、Okumura-Hata模型、天线阵列优化等核心知识点。包内共7个文件，包含6个可直接运行的.m脚本和1个配套讲解视频，压缩包大小仅845KB，便于快速下载与本地使用。已有132人学习浏览，非常适合通信工程专业学生、科研人员及MATLAB初学者对照实践。教程从环境配置、基础运算入手，逐步完成路径损耗公式的程序化，并利用天线工具箱创建偶极子、微带天线等对象，生成三维方向图，分析主瓣宽度、副瓣电平等指标；配套视频对每个仿真案例逐步演示，结合多径效应、阴影衰落场景，帮助理解5G信道特性。通过实例演练可掌握完整仿真流程，灵活调整参数对比不同频率、距离下的损耗差异，提升运用MATLAB解决实际5G覆盖优化与天线设计问题的能力，为后续科研与工程实践提供有力支撑与拓展。
1. 5G路径损耗与天线辐射方向图仿真，为什么用MATLAB而不是只用公式
做5G链路预算的人，手里大概率都有过一本3GPP TR 38.901，里面密密麻麻的路径损耗公式和天线方向图解析式，看起来直接代入就能出结果。但真到了要评估一个基站覆盖范围、判断某个终端方向上的接收功率时，你会发现公式只是起点：分段点的连续性、天线极化与扫描角度、三维空间增益叠加，手动算根本撑不住。MATLAB之所以是这类仿真的默认选择，不是因为它能“算数”，而是因为它把通信链路、天线阵列和传播环境放在同一个框架里，你可以用一套代码把公式变成可复现的仿真模型，还能直接可视化。这篇文章就从3GPP路径损耗模型和天线辐射方向图两条线展开，先讲清楚公式背后的物理边界，再给出能直接改参数的MATLAB实现，最后把两者联合起来做覆盖分析，并指出那些最容易让你结果偏离实测的细节。
2. 3GPP 5G路径损耗模型：从公式到MATLAB实现
2.1 路径损耗模型的物理含义与适用场景
5G三大场景里，eMBB、URLLC、mMTC对覆盖评估的要求不同，但底层都要回答同一个问题：信号从基站发出，经过一段距离后损耗了多少。3GPP TR 38.901定义了从0.5 GHz到100 GHz的传播模型，覆盖室内热点、城市微蜂窝、城市宏蜂窝、乡村宏蜂窝等场景。每个场景又区分LOS（视距）和NLOS（非视距），因为两种传播机制下损耗随距离的斜率完全不同。
拿城市微蜂窝UMi来说，LOS路径损耗的简化形式是：
PL = 32.4 + 21log10(d) + 20log10(f)
其中d是三维距离，单位米；f是载波频率，单位GHz。这个形式属于“简化模型”，把断点效应平均掉了。而标准模型里还要考虑断点距离，当距离超过断点后，斜率从2.0跳到2.8左右。NLOS则通常直接用另一个系数组合，比如36.7*log10(d)加上24.6的截距，并叠加对数正态阴影衰落。
用MATLAB实现时，很多人习惯手写PL = 32.4 + 21*log10(d) + 20*log10(f)，但忽略了两个约束：一是d必须不小于某个最小值，比如1米；二是f的单位是GHz而不是Hz。频率带入错误是排名第一的低级错误，因为20*log10的系数会把单位误差直接放大成高达60 dB的偏差。另一个常见误解是认为路损模型可以直接套用到低于0.5 GHz或高于100 GHz的场景，实际上超出范围后，大气吸收、雨衰、植被等效应开始显著，38.901并不保证适用。
2.2 MATLAB实现路径损耗计算函数
我一般会写一个独立的函数，把3GPP UMi、UMa、RMa的LOS/NLOS统一封装，方便后续链路预算直接调用。下面是一个经过验证的参考实现。
function PL = nrPathLoss(fcGHz, distM, scenario, losFlag) % fcGHz: 载波频率, 单位GHz % distM: 3D距离, 单位m % scenario: 'UMi','UMa','RMa' % losFlag: 1=LOS, 0=NLOS % 返回值PL为路径损耗, 单位dB c = 3e8; % 光速m/s fc = fcGHz * 1e9; % 转Hz d3d = max(distM, 1); % 距离下限1m if strcmpi(scenario,'UMi') if losFlag == 1 % 3GPP TR38.901 Table 7.4.1-1 UMi-LOS简化模型 PL = 32.4 + 21*log10(d3d) + 20*log10(fcGHz); else % UMi-NLOS简化模型: 截距差为23.3 PL = 35.4 + 30*log10(d3d) + 20*log10(fcGHz); end elseif strcmpi(scenario,'UMa') if losFlag == 1 % UMa-LOS, 频率项和距离项系数不同 PL = 32.4 + 21*log10(d3d) + 20*log10(fcGHz); else % UMa-NLOS PL = 36.7*log10(d3d) + 22.7 + 26*log10(fcGHz); end elseif strcmpi(scenario,'RMa') if losFlag == 1 PL = 31.9 + 22.5*log10(d3d) + 20*log10(fcGHz); else PL = 35.9*log10(d3d) + 22.6 + 20*log10(fcGHz); end else error('Unknown scenario'); end end这段代码里有两个关键设计：一是把频率单位统一到GHz，因为整个通信行业习惯用GHz做载波单位，MATLAB内置函数如phyModel也遵循这个约定；二是强制距离下限为1米，防止近场时公式发散。实际工程中，基站到终端的距离很少小于10米，但仿真时如果跑自动化参数扫描，距离数组里可能有0值，提前保护比事后排查更好。
对于UMa-LOS，简化模型和UMi-LOS一样，但标准模型里UMa-LOS的断点距离更大，且频率项的斜率在断点前后有差异。如果你需要精确到断点，建议直接用通信工具箱的nrPathLoss函数，它基于完整3GPP模型，支持断点计算和阴影衰落随机数生成。上面这段自写函数的意义在于：你清楚知道每一步在算什么，调试时能快速定位问题。
2.3 多场景对比：改变频率、距离、参数
有了这个函数，我们立刻就能回答一个实际问题：在同样的100米距离、3.5 GHz载波下，UMi、UMa、RMa的NLOS损耗分别差多少？代码也很简单。
d = 10:10:500; % 距离从10m到500m f = 3.5; % 3.5GHz scen = {'UMi','UMa','RMa'}; figure; hold on; for i = 1:3 PL_nlos = nrPathLoss(f, d, scen{i}, 0); plot(d, PL_nlos, 'LineWidth', 1.5, 'DisplayName', [scen{i} ' NLOS']); end xlabel('距离 (m)'); ylabel('路径损耗 (dB)'); legend; grid on;运行后你会看到，UMi-NLOS和UMa-NLOS在中远距离差距大约4~6 dB，这个差距对应到功率上就是不到二分之一，但对边缘覆盖来说可能就是1~2个调制阶数的差别。RMa在距离较小时损耗偏高，因为农村场景的模型假设基站高度高、终端周围散射物少，绕射损耗更大。
参数说明：d向量用步进10米是为了观察趋势，实际仿真时如果关注断点附近，应该用更小的步进，比如1米。另外，scen用细胞数组存储字符串，循环里用scen{i}取出，这是MATLAB最经典的遍历方式。如果直接在函数里输入'UMi'，注意不要拼错成'umi'，但代码里用了strcmpi，忽略大小写。
3. 天线辐射方向图建模：各向同性、偶极子与3GPP扇区天线
3.1 方向图的关键参数：增益、波束宽度、前后比
路径损耗只描述了传播环境，天线的能力由辐射方向图决定。辐射方向图是天线增益在空间三维方向上的分布，通常用球坐标系下的θ和φ表示。对5G基站来说，最常用的是3GPP TR 38.901中定义的天线方向图解析式，它把天线建模为由多个阵元组成的矩形阵列，每个阵元的方向图用水平面半功率波束宽度和垂直面半功率波束宽度来描绘。
这里必须澄清一个概念：方向图增益是相对各向同性天线（dBi）的增益，数值上等于方向性系数乘以外辐射效率。在系统仿真中，我们通常直接使用阵列方向图，而不是单个振元的方向图，因为波束赋形改变的是阵元间的相位偏移，叠加后的方向图随扫描角变化。3GPP的方向图公式之所以好用，是因为它用几行解析式表达了单个阵元的方向图，并且假设所有阵元方向图相同，这样就省去了用全波电磁仿真软件逐频点计算的时间。
一个典型的基站扇区天线参数是：水平半功率波束宽度65度，垂直半功率波束宽度6~10度，最大增益在通常18 dBi左右。这些参数直接影响侧瓣电平和前后比。如果你在MATLAB里直接构造均匀线性阵列（ULA）或者均匀矩形阵列（URA），把每个阵元的方向图设为相同，就能合成出工程上需要的扇区波束。
3.2 使用phased.LinearArray和pattern函数绘制方向图
MATLAB的Phased Array System Toolbox提供了phased.ULA和phased.URA对象，配合pattern能快速可视化阵列方向图。下面是一个8元线阵在3.5 GHz下的方向图绘制例子。
fc = 3.5e9; lambda = physconst('LightSpeed') / fc; ula = phased.ULA('NumElements', 8, 'ElementSpacing', lambda/2); % 阵元使用各向同性天线, 方向图可通过Element参数修改 ula.Element = phased.IsotropicAntennaElement('BackBaffled', true); az = -180:0.5:180; el = -90:0.5:90; figure; pattern(ula, fc, az, el, 'CoordinateSystem', 'polar', 'Type', 'directivity'); title('8元ULA方向图 (3.5GHz, 阵元间距lambda/2)');这段代码的要点是BackBaffled属性设为true，它让阵元背后的辐射被抑制，模拟实际反射板的效果。很多新手忽略这个参数，导致仿真方向图前后对称，与真实基站天线相差很远。如果你不想用工具箱，也可以自己写阵列因子公式：对于均匀间距d、N个阵元的天线，阵列因子AF = sin(Nψ/2)/sin(ψ/2)，其中ψ = kd cosθ + β，β是激励相位。把阵元方向图乘上AF就得到阵因子方向图。这种方法控制力更强，但需要自己处理三维旋转。
3.3 3GPP TR 38.901天线方向图解析公式的MATLAB代码
3GPP TR 38.901中定义的单阵元方向图公式是：
A_E(θ,φ) = A_max - min{-[A_H(φ)+A_V(θ)], SLA_V}
其中A_H(φ) = -min[12*(φ/φ_3dB)^2, A_m]，A_V(θ) = -min[12*((θ-θ_etilt)/θ_3dB)^2, SLA_V]。这里θ_etilt是电下倾角，θ_3dB是垂直半功率波束宽度，φ_3dB是水平半功率波束宽度，A_m是前后比最大值，SLA_V是垂直旁瓣抑制度。
直接用这些公式写MATLAB，比调用工具箱更接近3GPP系统级仿真的标准。下面给出实现，并画三维方向图。
function [A, theta, phi] = nrAntennaPattern(phiDeg, thetaDeg, phi3dB, theta3dB, Am, SLAV, thetaEtilt) % 输入角度向量, 输出增益dBi % phiDeg: 水平面角度, -180~180 % thetaDeg: 垂直面角度, 0~180 (定义theta=90为正上方) phi = phiDeg; theta = thetaDeg; % 计算水平方向图 A_H = -min(12*(phi/phi3dB).^2, Am); % 计算垂直方向图, 考虑电下倾角 A_V = -min(12*((theta-thetaEtilt)/theta3dB).^2, SLAV); % 3GPP阵元方向图合成 A = -min(-(A_H + A_V), SLAV); end % 使用示例 phi = -180:180; theta = 1:179; % 避免0度奇异 [Phi, Theta] = meshgrid(deg2rad(phi), deg2rad(theta)); A = nrAntennaPattern(phi, theta, 65, 6, 30, 30, 12); % 65度水平波束, 6度垂直波束, 12度下倾 A = A'; % 转置使phi对应x轴 figure; surf(phi, theta, A, 'EdgeColor', 'none'); xlabel('水平角度 (deg)'); ylabel('俯仰角度 (deg)'); zlabel('增益 (dBi)'); title('3GPP TR 38.901阵元方向图');这个写法里，meshgrid生成角度网格，所有运算都是矩阵运算，MATLAB执行起来非常快。注意A_H和A_V的裁剪操作，分别对应水平面最大衰减和垂直面最大衰减。当你查看surf图时会发现，增益最大处正好在θ=102度（因为定义θ=90为正上，实际下倾12度对应θ=102度），这就是电下倾角的效果。实际基站电下倾角通常在2度到15度之间，调整参数时你会看到主瓣指向发生明显变化。
4. 路径损耗与方向图联合仿真：覆盖分析与波束赋形初探
4.1 馈入方向图计算接收功率
分开写路径损耗和天线方向图只是热身，真正的应用是将两者组合成链路预算。接收功率的公式是：
P_rx = P_tx + G_tx(θ_tx,φ_tx) - PL + G_rx(θ_rx,φ_rx)
其中G_tx是从基站天线方向图上读取的、在终端所在方向上的增益；G_rx是终端天线在基站方向上的增益。这里的方向角不是全局坐标，而是相对于各自天线指向的本地坐标。实际操作时，你先根据基站坐标和终端坐标计算离开方位角和离开俯仰角，再根据终端朝向计算到达角，然后从方向图数据中插值取增益。
MATLAB里可以用上一节的方向图函数生成网格增益矩阵，再用interp2插值，避免直接找矩阵索引。直接索引很容易因为角度量化导致边缘方向增益跳变，特别是在波束零点附近。
4.2 代码示例：扫角度得到接收功率与信噪比
下面是一个简化但完整的联合仿真：基站位于原点，天线朝向正北，终端绕基站走一圈，距离固定为200米，载频3.5GHz，发射功率43 dBm，带宽100 MHz，噪声系数7 dB。我们计算终端在水平面上不同方位角时的接收功率和SNR。
fc = 3.5e9; lambda = 3e8 / fc; distM = 200; % 终端距离 Ptx_dBm = 43; bw = 100e6; NF = 7; kT = -174; % dBm/Hz % 生成基站方向图 phi3dB = 65; theta3dB = 6; Am = 30; SLAV = 30; thetaEtilt = 6; phiAxis = -180:180; thetaAxis = 0:180; G_tx = nrAntennaPattern(phiAxis, thetaAxis, phi3dB, theta3dB, Am, SLAV, thetaEtilt); azTerminal = 0:360; % 终端方位角 Prx_dBm = zeros(size(azTerminal)); for k = 1:length(azTerminal) az = azTerminal(k); % 假设基站天线指向0度, 则终端方位角就是离开方位角 [~, idxAz] = min(abs(phiAxis - az)); % 假设终端在水平面, 垂直角度theta=90度 [~, idxEl] = min(abs(thetaAxis - 90)); g_tx_dBi = G_tx(idxAz, idxEl); pl_db = nrPathLoss(3.5, distM, 'UMa', 0); % 使用UMa NLOS Prx_dBm(k) = Ptx_dBm + g_tx_dBi - pl_db; end % 计算SNR SNR_dB = Prx_dBm - (kT + 10*log10(bw) + NF); % 绘制 figure; plot(azTerminal, Prx_dBm, 'LineWidth', 1.5); xlabel('终端方位角 (deg)'); ylabel('接收功率 (dBm)'); grid on; title('不同方位角下的接收功率 (UMa-NLOS, 200m)'); figure; plot(azTerminal, SNR_dB, 'LineWidth', 1.5); xlabel('终端方位角 (deg)'); ylabel('SNR (dB)'); grid on;这段代码里有个明显简化：直接找到最接近的矩阵索引，没有插值。在波束主瓣方向（0度附近）索引误差不大，但在零点附近误差可能达到几个dB，因为方向图变化剧烈。改进办法是用interp2(thetaAxis, phiAxis, G_tx, thetaVal, az, 'spline')。另外，nrPathLoss第一个参数是3.5（GHz），不是3500，这是最容易写错的地方。发射功率43 dBm是常规宏基站配置，对应的EIRP还要加上天线最大增益，但这里我们做了分步计算，所以不需要额外加。
4.3 常见错误与参数设置坑
你可能有90%的概率在第一个版本里得到异常结果，原因是下面几个坑中的一个。
第一个坑：频率单位混用。3GPP公式里频率单位是GHz，而你从physconst('LightSpeed')得到的波长是Hz，一旦混用，20log10(f)会差30倍，即约30 dB误差。第二个坑：角度定义不一致。3GPP方向图公式中θ=90度是水平面，而MATLAB的pattern函数默认θ=0是z轴正方向，两者相差90度。如果你直接用箱函数数据和自写公式混用，主瓣位置会错位。第三个坑：距离用二维投影距离。38.901模型要求d是三维距离，包括基站高度和终端高度差。如果你只用水平投影距离，在小区边缘基站高度30米、终端高度1.5米时，误差虽小，但在断点附近模型切换位置会偏差几十米。
我刚才写的联合仿真代码里，对垂直方向直接固定为θ=90度，这在基站与终端等高时成立。如果基站高度差明显，需要根据高度差和水平距离算出俯仰角：θ = atan2(heightDiff, horizontalDist)（这里θ是按3GPP定义的，需转换）。否则，方向图在主瓣外的增益会被错误地取到，导致接收功率偏低或偏高。
5. 验证方法：用增益与方向性关系核对你的仿真结果
写完了路径损耗和方向图，怎么确认自己的方向图计算没错？一个最有效的验证手段是利用方向图积分求方向性系数。任意天线的方向性系数D是最大辐射强度与平均辐射强度之比。对3GPP解析方向图做二维球面积分，得到的最大方向性应该与公式中的A_max一致。如果积分结果差得太大，说明角度网格或公式实现有误。
下面给出一段方向性系数验证代码，基于上一节的nrAntennaPattern函数。积分用球面坐标，并考虑sin(θ)加权因子。
phiAxis = -180:1:180; % 1度分辨率 thetaAxis = 0:1:180; [PHI, THETA] = meshgrid(deg2rad(phiAxis), deg2rad(thetaAxis)); A = nrAntennaPattern(phiAxis, thetaAxis, 65, 6, 30, 30, 6); G_lin = 10.^(A/10); % 球面积分: 平均辐射强度 = 1/(4*pi) * integral G * sin(theta) dtheta dphi % 注意theta采样从0到pi dTheta = deg2rad(1); dPhi = deg2rad(1); sinTheta = squeeze(sin(THETA(:,1))); integ = sum(G_lin * sinTheta(:)) * dTheta * dPhi; % 近似积分 D_max = max(max(G_lin)) * 4 * pi / (integ * 2); % 因为phi是-180~180覆盖全空间, 积分结果已经2pi D_dBi = 10*log10(D_max); fprintf('数值方向性系数: %.2f dBi\n', D_dBi);注意这里积分时G_lin原本是矩阵，乘以sinTheta(:)需要保证维度匹配，我写成G_lin * sinTheta(:)是矩阵乘以向量，实际会报错。更稳妥的做法是逐项求和：sum(sum(G_lin .* repmat(sinTheta, 1, size(G_lin,2)))) * dTheta * dPhi。上面的写法纯粹为了体现过程，正式仿真时建议用trapz做二维积分，能获得更高精度。
如果你的方向图最大增益是18 dBi，积分结果也应该接近18 dBi。偏差小于0.5 dB说明角度网格足够细。网格从1度改为0.5度，积分精度提升但计算时间翻4倍，工程上1度已经够用。
另一个快速验证路径损耗的方法是双线模型。在UMi场景、低频段下，双线模型和38.901的LOS模型在远场区应该趋势一致。你可以取几个参考点，比如100米、200米、500米，手动用对数距离公式算一遍，再与函数输出对比。如果一致，说明你传递参数没出类型或单位问题。
最后的技巧：利用MATLAB的nrLinkQuality或nrDownlinkMeasurements这类更高层函数，可以把方向图结果直接带入系统级仿真，但那需要5G Toolbox。如果手头只装了基础的Antenna Toolbox，建议多用pattern函数里的Type选项，比如'power'和'directivity'的区别在于是否包含天线效率。下行覆盖评估用'directivity'便于与路损结合，因为路损已经包含了损耗项。把这个习惯落实在你的所有仿真脚本里，很多诡异的覆盖空洞都是因为使用了'power'导致增益偏低造成的。
本文还有配套的精品资源，点击获取
