5G路径损耗与天线方向图MATLAB仿真：从模型选择到联合链路预算-北京尧图网络

简介围绕5G路径损耗与天线辐射方向图的MATLAB仿真教程面向通信工程学生、无线网络规划人员及MATLAB初学者帮助掌握高频段信号传播损耗计算与天线增益、方向图分析方法。资源包共7个文件含6个可直接运行的.m脚本和1份教学讲解mp4压缩后仅845KB轻量易得。已有132人学习适合课程设计或自学辅助。教程从MATLAB基础操作、Friis传输方程与Okumura-Hata等路径损耗模型到天线辐射机理与三维方向图生成均配有具体仿真示例同时结合5G多径效应、阴影衰落场景演示不同传播条件下的链路预算评估并利用天线工具箱模拟偶极子、微带、相控阵等天线分析主瓣宽度、副瓣电平等指标。通过脚本注释与视频讲解可快速上手基于MATLAB的无线通信仿真为5G网络规划与天线优化提供实用技能。1. 5G路径损耗与天线辐射方向图先搞懂这两件事在仿真里是什么做5G链路仿真时见过太多人拿着频段和功率就开始叠链路预算算出来的接收功率却和实际路测对不上。问题多半不在收发机而在两处路径损耗模型没按场景选对天线方向图被当成全向点源处理。路径损耗决定信号随距离衰减的“骨架”天线方向图决定能量在空间里的“分布”两者叠加才是真实的覆盖结果。这篇博文围绕“matlab 5G路径损耗和天线辐射方向图的matlab仿真”展开给出从模型选择、方向图构建到联合仿真、参数排错的完整套路并用MATLAB把每一步写成可运行代码。无论你是做覆盖仿真的无线工程师、做物理层算法验证的研究生还是刚接手5G基站天线指标的测试开发都能直接照着改参数跑出自己的结果。2. 5G路径损耗模型怎么选从自由空间到TR 38.901的MATLAB实现2.1 路径损耗模型的适用场景与频段边界3GPP TR 38.901是现在5G NR链路预算和覆盖仿真的默认参考模型里面定义了UMi、UMa、RMa三类场景分别对应城区微站、城区宏站、郊区宏站。和4G时代常用的Hata/COST-231模型比38.901最大的变化是显式区分了LOS/NLOS并且把断点距离breakpoint distance作为分段函数的切换点这对5G高频段特别重要。选择模型时要先问三个问题基站高度和用户高度差是多少用户距离集中在10米到500米还是更长工作在FR1的3.5GHz频段还是FR2的28GHz频段自由空间模型只适合无遮挡的视距短距场景比如室内热点或回传链路真正的蜂窝覆盖仿真建议直接用nrPathLoss函数或按38.901分段公式手写。38.901模型的另一个特点是频率适用范围扩展到0.5GHz到100GHz这正好覆盖5G的n77/n78/n79频段。若用老模型外推到毫米波算出来的路径损耗会明显偏高导致基站数量被高估。下面给出一个可直接粘贴到MATLAB里的38.901 UMi手写函数不使用任何工具箱也能工作。2.2 手写TR 38.901 UMi路径损耗函数function pl pl_tr38901_umi(fc, d3d, los, hBS, hUT) % fc 单位 GHzd3d 单位 mlos 为 true/false % hBS、hUT 单位 m用于计算断点距离 dBP c 2.99792458e8; dBP 4 * (hBS - 1) * (hUT - 1) * (fc * 1e9) / c; if los if d3d dBP pl 32.4 20 * log10(fc) 21 * log10(d3d); else pl 32.4 20 * log10(fc) 40 * log10(d3d) - 19 * log10(dBP); end return; end % NLOS 时取 LOS 与 NLOS 表达式的较大者这是 38.901 的常用简化形式 pl_los pl_tr38901_umi(fc, d3d, true, hBS, hUT); pl_nlos 35.3 22 * log10(d3d) 12.5 * log10(fc); pl max(pl_los, pl_nlos); end核心逻辑是先算出断点距离dBPLOS时在断点前用21的斜率断点后用40的斜率NLOS时用22的斜率但最终取LOS和NLOS中较大的那个因为实际传播在NLOS场景下不可能优于LOS。调用时注意fc和d3d的单位不能搞混fc是GHzd3d是米unit的换算都在公式内部处理。例如仿真3.5GHz、基站高25米、用户高1.5米、距离200米的情况pl pl_tr38901_umi(3.5, 200, false, 25, 1.5); fprintf(UMi NLOS path loss at 200m: %.2f dB\n, pl);这段代码在LOS与NLOS交界处的过渡比官方查表要平滑适合快速链路预算。若使用5G Toolbox直接用nrPathLoss函数代替上面的手写实现输出一致且支持向量距离输入。2.3 用nrPathLoss替代手写5G Toolbox的官方实现5G Toolbox中nrPathLoss函数的语法是pl nrPathLoss(model, fc, d)model支持free-space、UMi、UMa、RMa四种。下面的代码画一条0到500米的路径损耗曲线顺便对比自由空间模型直观看出遮挡损耗的量级d 10:5:500; fc 3.5; pl_fs nrPathLoss(free-space, fc, d); pl_umi_nlos nrPathLoss(UMi, fc, d); plot(d, pl_fs, --, d, pl_umi_nlos, -, LineWidth, 1.5); grid on; legend(Free Space, UMi NLOS, Location, northwest); xlabel(Distance (m)); ylabel(Path Loss (dB)); title(3.5GHz Path Loss Comparison);nrPathLoss内部的UMi默认NLOS概率与距离相关在宏站场景下直接用UMa更合适。参数fc以GHz为单位d可以传向量函数内部自动做逐点计算。用free-space和UMi对照的好处是两者差距就是阴影衰落和绕射损耗的近似值能帮助判断链路是否还有余量。2.4 用实测CSV校正模型从路测数据到仿真参数从现场路测拿到的RSRP数据通常以CSV形式存在可以先导入MATLAB做预处理再与仿真值对比。常见做法是用readmatrix读取CSV再按距离分箱并计算每个bin的平均RSRP差从而修正仿真模型的偏置data readmatrix(drive_test.csv); % 假设第1列距离第2列RSRP d_meas data(:, 1); rsrp_meas data(:, 2); [tf, idx] ismember(round(d_meas / 10) * 10, round(d / 10) * 10); if ~all(tf) warning(部分距离超出仿真范围请检查CSV数据); end这段脚本把实测距离舍入到10米分箱去匹配仿真距离向量方便后续做差。注意readmatrix在R2021a之后的中文路径下可能报错建议把CSV放在纯英文目录再导入。数据校正不是一次完成的事通常要跑多轮将模型的损耗偏置作为常数补偿量加进后续链路预算。3. 天线辐射方向图建模从各向同性到5G阵列单元方向图3.1 方向图的数据结构与phased.CustomAntennaElement路径损耗告诉我们“信号走多远衰减多少”方向图则回答“每个方向上辐射多少能量”。5G Massive MIMO的方向图不是简单一个喇叭或贴片而是由阵列因子与单元因子相乘得到。MATLAB的Phased Array System Toolbox把这两种因子分开建模单元方向图定义在方位角和俯仰角的网格上阵列用phased.URA、phased.UCA等方式描述最终由pattern函数合成方向图。先用各向同性单元演示最简单的全向方向图再到贴片、喇叭和阵列。CustomAntennaElement可以指定方位角、俯仰角网格和对应幅度响应适合导入天线厂商提供的实测文件。需要注意的是MagPattern矩阵的行数必须等于俯仰角数量列数等于方位角数量否则报错DimensionMismatch。3.2 构建一个5G基站贴片单元方向图并做3D可视化5G宏站通常使用贴片天线作为阵列单元水平半功率波束宽度约65度垂直面约6到12度。为方便教学用高斯函数近似贴片单元的水平切面az -180:2:180; % 方位角2度步进 el -90:2:90; % 俯仰角2度步进 [az_mesh, el_mesh] meshgrid(az, el); % 高斯束近似水平3dB宽度约65度垂直面用余弦压缩 response exp(-(az_mesh.^2) / (2 * 35^2)) .* cosd(el_mesh); patch_elem phased.CustomAntennaElement(... AzimuthAngles, az, ... ElevationAngles, el, ... MagPattern, response);这段代码中35度的标准差对应约65度的半功率波束宽度。cosd(el_mesh)让单元在水平面上辐射最强向下和向上随俯仰角衰减。建好单元后直接绘制单个单元的空间方向图pattern(patch_elem, 3.5e9, -180:1:180, -90:1:90, ... CoordinateSystem, polar, Type, powerdb); colorbar;pattern函数中间的坐标参数fc是载波频率单位为Hz。CoordinateSystem设为polar更直观看到波瓣指向z轴方向。要导出场强值做后处理可以用pattern函数带返回值的形态[pat, az_out, el_out] pattern(patch_elem, 3.5e9, az, el);pat是一个len(el_out)行len(az_out)列的矩阵单位默认是dBi。拿到这个矩阵后可以自己画等值线图或乘以发射功率后作为方向性增益叠加到路径损耗上。3.3 8x8 UPA阵列方向图单元因子与阵列因子如何叠加单个贴片方向图只能覆盖一个波束5G基站实际用的是8x8或更大规模的均匀平面阵列。构建UPA阵列关键是设置阵元间距常见为半波长。工作在3.5GHz时波长是0.0857米半波长约0.0428米fc 3.5e9; lambda physconst(LightSpeed) / fc; array phased.URA(Size, [8 8], ... ElementSpacing, [lambda/2 lambda/2], ... Element, patch_elem);phased.URA的ElementSpacing用的是[m n]形式分别对应行方向和列方向的间距。如果不指定Element默认使用各向同性单元。绘制方向图pattern(array, fc, -180:1:180, -90:1:90, ... CoordinateSystem, rectangular, Type, powerdb);对比单元方向图和阵列方向图能发现阵列方向图的主瓣明显变窄旁瓣出现在固定角度上而且栅瓣出现的角度由阵元间距决定。当阵元间距超过一个波长时可见空间内会出现与主瓣增益相当的栅瓣这就是为什么5G Massive MIMO的标准间距都取半波长。3.4 方向图参数对覆盖的影响半功率波束宽度与阵元数UPA的主瓣宽度由阵元数决定不取决于单元方向图。8元线阵在法向的半功率波束宽度约为0.886lambda/(Nd)弧度代入N8、dlambda/2得到约12.7度。把8乘8阵列的主瓣增益近似为10*log10(64)单元增益通常在24到27dBi之间。这些数值会影响5G波束管理的波束数量设计波束越窄覆盖同样扇区需要越多波束。4. 5G基站天线与路径损耗联合仿真一个完整的链路预算例子4.1 波束赋形后的EIRP计算有了路径损耗和方向图就可以从发射端一路算到接收端。先把“功率、增益、损耗”统一到dB域计算EIRPptx_w 10^(30 / 10) / 1000; % 30 dBm - 1W array_gain_dbi 24; % 8x8 UPA主瓣方向近似增益 eirp_dbm 30 array_gain_dbi; % EIRP 发射功率阵列增益 fprintf(EIRP %.2f dBm\n, eirp_dbm);EIRP是网络部署中受法规限制的指标3GPP和各国监管机构都规定了不同频段的最大EIRP。实际工程中阵列增益需在主瓣方向减去馈线损耗、插入损耗等通常按2到3dB预留。4.2 接收功率公式的MATLAB表达路径损耗与方向图合流用户位于阵列方向图的某个偏角时不能直接使用最大增益。要同时查天线在真实方向上的增益值和路径损耗。完整的接收功率表达式为tx_power_dbm 30; rx_lat 31.2304; rx_lon 121.4737; % 示例坐标 dist 120; % 单位m pl nrPathLoss(UMa, 3.5, dist); % 宏站场景 % 计算用户相对基站的方位角和俯仰角简化示例 az_to_user 35; el_to_user -8; % 利用pattern返回单元方向图增益矩阵 [pat, az_out, el_out] pattern(array, 3.5e9, -180:1:180, -90:1:90); [~, az_idx] min(abs(az_out - az_to_user)); [~, el_idx] min(abs(el_out - el_to_user)); g_user_dbi pat(el_idx, az_idx); rx_power_dbm tx_power_dbm g_user_dbi - pl 14; % 接收天线增益约14dBi fprintf(RSRP %.2f dBm\n, rx_power_dbm);注意这行代码的天线增益g_user_dbi取自双线性插值直接用min(abs(...))会损失零点几dB但对链路预算足够。接收天线增益14dBi是终端或CPE的典型值。真实系统中如果用户处于阵列旁瓣方向图增益可能比主瓣低20dB以上只看EIRP会高估覆盖这也是为什么LTE时代“路损万金油”的做法在5G行不通。4.3 在3D地图上画覆盖栅格算法逻辑与数据出口覆盖仿真最终要落在栅格上常见做法是把基站周围区域划分为10米或20米分辨率网格逐点计算路径损耗与天线增益叠加结果。下面的代码演示核心循环并给出输出矩阵x -500:20:500; % 东西方向单位m y -500:20:500; % 南北方向单位m [x_grid, y_grid] meshgrid(x, y); d_grid sqrt(x_grid.^2 y_grid.^2); az_grid atan2d(x_grid, y_grid); % 方位角 el_grid atan2d(3, d_grid) * 180 / pi; % 假设基站相对地面高3m得到俯仰角 gain_grid zeros(size(d_grid)); for i 1:numel(d_grid) [~, az_idx] min(abs(az_out - az_grid(i))); [~, el_idx] min(abs(el_out - el_grid(i))); gain_grid(i) pat(el_idx, az_idx); end pl_grid nrPathLoss(UMa, 3.5, d_grid); rsrp_grid tx_power_dbm gain_grid - pl_grid 14;这段循环用了向量化前的直观写法数据量大时可用interp2替代两次min查找。栅格矩阵rsrp_grid可直接用imagesc显示成热力图也可以写成CSV给GIS工具或PostGIS做渲染。栅格物理位置的经纬度转换需要引入UTM投影或球面坐标换算仿真中常用等距近似。4.4 结合PSS/SSS同步信号做覆盖验证5G NR的PSS/SSS在时频结构上有固定位置SSB波束扫描顺序可配置。做覆盖仿真时经常需要把RSRP计算结果对应到具体SSB index上。5G Toolbox可用nrSSBurst配置SSB的频域位置和波束数nrb 273; % 30kHz SCS下100MHz带宽的RB数 ssb nrSSBurst; ssb.BlockPattern Case B; ssb.FrequencySSB nrb * 12 / 2; % 把SSB放在载波中心 info nrSSBurstInfo(ssb);但这一步只决定SSB在时频格上的位置真正的RSRP测量值是信道估计后在SSB时频资源上取平均得到的。在纯链路预算层面用4.2节计算出的RSRP作为SSB参考信号功率即可不需要展开OFDM波形级仿真。若要做波形级验证则要将路径损耗和方向图结果放入nrPerfectChannelEstimate中构建信道再对PSS/SSS解调做互相关峰值检测。5. 参数设置与仿真坑点5G频段、单位、阵列互耦5.1 n77/n78/n79频段对应的模型参数5G FR1的n773.3-4.2GHz、n783.3-3.8GHz、n794.4-5.0GHz在仿真中除了频率不同大部分链路预算结构一致。差别主要体现在三个地方一是路径损耗模型中的频率项二是阵元间距要按各自波长算三是不同频段的穿透损耗实测值不同。下面列一个常用参数表参数n78典型值n79典型值n77典型值中心频率3.5 GHz4.8 GHz3.7 GHz波长0.0857 m0.0624 m0.0810 m子载波间隔30 kHz30 kHz30 kHz每小区SSB波束数888UMI NLOS损耗100m约105 dB约109 dB约106 dB注意n79的路径损耗比n78高约3到4dB这意味着同站址要覆盖同样边界需要更高的发射功率或更小的站间距。若有穿透损耗需求可以在公式中额外增加5到15dB的building penetration loss具体按ITU-R P.2108取。5.2 dB与线性值换算陷阱log10、pow2db与db2pow链路预算最容易翻车的点是单位的混用。在MATLAB中dBm转W很常见但经常有人直接写ptx_w 10^(30/10)忘了除以1000。推荐的做法是使用内置函数pow_dbm 30; pow_w db2pow(pow_dbm - 30); % dBm - W pow_w2 (10^(pow_dbm / 10)) / 1000; % 等价写法 % 验证两者是否一致 assert(abs(pow_w - pow_w2) 1e-12);db2pow输入的是dB若输入dBm必须先减30转到dBW。反向转换用pow2db(pow_w * 1000)。在4.2节这种全公式都在dB域时不存在换算问题但若把路径损耗与发射功率直接相乘计算结果会错得离谱。这类毛病在接手别人的MATLAB脚本时尤其常见强烈建议在脚本开头用assert做单位自检。5.3 UPA的Taper与单元方向图的互耦近似实际5G天线阵列的每个单元之间存在互耦方向图会发生畸变波束指向也会偏移。全波仿真软件可以算精确结果但MATLAB的phased.URA默认不包含互耦效应。工程上常用两种近似第一种是在阵列方向图乘一个单元因子后再做幅度加权第二种通过Taper属性做泰勒加权来压低旁瓣taper hamming(8) * hamming(8).; array.Taper taper(:); pattern(array, fc, -180:1:180, -90:1:90, Type, powerdb);使用hamming窗可以将旁瓣压低到-40dB左右代价是主瓣展宽约1.3倍、增益损失约1dB。实际5G设备的波束配置中广播波束常采用较宽的垂直面波束业务波束才用窄波束高增益这个切换就是通过Taper和相位加权一起实现的。互耦的更高精度修正可以参考天线厂商提供的S参数但不属于纯MATLAB链路预算的范畴。5.4 MATLAB版本与工具箱缺失时的替代方案运行上述代码需要Communication Toolbox、Phased Array System Toolbox和5G Toolbox。如果机器上没有装5G Toolbox路径损耗部分用2.2节的手写函数即可完全不影响结果。方向图部分若没有Phased Array Toolbox可以用自编阵列因子公式叠加单元方向图公式为AF(theta, phi)sum(w_n * exp(jkdnsin(theta)*cos(phi)))。从工程角度看手写方向图的好处是方便做扫描角与幅相加权的调试缺点是没有内置的3D绘图接口。版本方面从R2021a开始5G Toolbox的nrPathLoss已经支持全部四种模型R2023b之后pattern绘制3D方向图的速度有大幅提升。代码迁移到老版本的主要坑是某些新函数名不可用比如nrSSBurstInfo在R2022a才加入用之前先执行exist检查更稳妥。即使考虑最近几版MATLAB安装包问题只要工具箱完整上述代码在新版本里都能直接运行。6. 把仿真做成可复用脚本一行命令生成“距离-损耗-方向图”组合图6.1 封装simulateCoverage.m函数真正好用的仿真不是一堆脚本到处粘贴而是把链路预算封装成一个函数。输入基站位置、频段、发射功率和波束指向输出RSRP栅格和方向图数据function [rsrp_grid, gain_grid] simulateCoverage(fc, tx_power_dbm, az_center, el_center) % fc: GHz; tx_power_dbm: 发射功率; az_center/el_center: 波束指向 lambda physconst(LightSpeed) / (fc * 1e9); array phased.URA(Size, [8 8], ... ElementSpacing, [lambda/2 lambda/2], ... Element, patch_elem); stv phased.SteeringVector(SensorArray, array, ... PropagationSpeed, physconst(LightSpeed)); sv stv(fc * 1e9, [az_center; el_center]); pattern(array, fc * 1e9, -180:1:180, -90:1:90, ... Weights, sv, Type, powerdb); % 随后复用第4.3节网格计算逻辑 end函数输出rsrp_grid后可以用subplot把三个信息画到一张图上子图1是“距离-路径损耗”曲线子图2是阵列方向图切面子图3是RSRP栅格热力图。这样做最大好处是当需要给基站或天线厂商汇报时同一行命令出的图可以直接贴进报告。将函数与脚本分离之后调参也不需要翻代码。6.2 波束扫描动画技巧在5G基站的波束管理测试中经常要观察波束从-60度扫描到60度时覆盖栅格的变化。在MATLAB中用循环配合drawnow即可生成动画for az_scan -60:5:60 sv stv(fc * 1e9, [az_scan; 0]); pattern(array, fc * 1e9, -180:1:180, -90:1:90, ... Weights, sv, Type, powerdb); hold on; plot3([0 500], [0 0], [10 10], r--, LineWidth, 2); drawnow; pause(0.1); end这段循环每步修改波束指向并重绘方向图红色虚线标出需要重点观察的覆盖方向方便快速判断主瓣是否对准目标区域。若屏幕上同时显示RSRP栅格需要把方向图与栅格计算结果都放入循环内复杂度会随栅格点数上升。数据量大时建议先把波束扫描到栅格的增益表预计算好再用矩阵乘法快速得到每个扫描角下的覆盖图。将方向图与路径损耗结合之后这套仿真还能继续扩展到28GHz频段、RIS反射面辅助覆盖、以及网络优化中的自动波束对齐场景只需要替换路径损耗模型和天线阵列参数即可。用这个思路5G链路预算从一个Excel表变成了可交互、可复现、可扩展的MATLAB工具链。本文还有配套的精品资源点击获取

5G路径损耗与天线方向图MATLAB仿真：从模型选择到联合链路预算

相关推荐

SpringBoot OAuth2授权登录实战与优化

Vue3个人网站模板源码：组合式API与数据驱动实践

Comsol中计算陈数：从Berry曲率到Wilson loop的完整仿真流程

更多资讯

基于土壤湿度控制继电器的 Wio Terminal 实战指南（IoT-For-Beginners 自动化植物浇水项目）

VisionMaster授权更新报错“本地LM通讯出错”排查与解决

Sanity React 性能规则详解:订阅派生布尔状态,降低组件重渲染频率

FPGA加法器实验全解析：从Verilog代码到上板调试

国产低噪声LDO LTP7792实测：从参数到ADC供电设计避坑指南

cvtColor内存泄漏排查：从Mat生命周期到输出缓冲区复用

今日更新

AWS无服务器应用开发指南：从Lambda到SAM的架构与实践

SpringBoot+Vue人事系统实战：解决HR数据断点与事务一致性

MATLAB手写紧束缚模型计算石墨烯能带结构

本周精选

AI驱动的自动化专利分析平台技术解析

LeetCode-Go 题解：1252. Cells with Odd Values in a Matrix（奇数值单元格计数）双解法剖析

在线答疑系统Java毕设实战：Spring Boot前后端分离与状态流转

本月热门

基于YOLO的焊缝缺陷检测：从数据标注到工业部署全流程实践

分布式服务的决策记录

创业产品的首版范围

读完文章，想聊聊您的网站？