rng (“默认”）;％重置随机数发生器simParameters = [];%清除simParameters变量simParameters。NumFramesSim = 20;%以10毫秒帧数计算的模拟时间simparameters.schedulingtype = 0;%设置值为0(基于插槽调度)或1(基于符号调度)。

simParameters。NumCells = 3;％细胞数量Simparameters.cellradius = 500;每种细胞的％半径（以米为单位）

simParameters。GNBPosition= [1700 600; 3000 600; 2500 2000];%单元内UEs的个数。每个单元包含相同数量的uesimParameters。NumUEsCell = 4; simParameters.GNBTxPower = 32;GNB的％TX功率（在DBM中）％将信道带宽设置为10 MHz和子载波间距（SCS）至30 kHz%定义见3GPP TS 38.104章节5.3.2。simparameters.numrbs = 24;simparameters.scs = 30;％kHz.simParameters。DLCarrierFreq = 2.635 e9;％Hz.simParameters。DLBandwidth = 10 e6;％Hz.％SINR到CQI映射表0.1 BLERSimparameters.sinr90pc = [-5.46 -0.46 4.54 9.05 11.54 14.04 15.54 18.04......20.04 22.43 24.93 25.43 27.43 30.43 33.43];对PDSCH分配的RB的最大限制％。传输限制是％适用于新的PDSCH分配，而不是用于重新定位simparameters.rballocationlimitdl = 15;PDSCH的百分比

％指定感兴趣的小区ID。所有可视化和度量标准都显示为此单元格simparameters.cellofInterest = 2;％将0从0设置为numcells-1% CQIVisualization和RBVisualization参数控制显示%这些可视化:(i) CQI可视化RBs (ii) RB赋值％可视化。默认情况下，“rbvisualization”绘图已禁用。你%可以通过设置为“true”来启用它simparameters.cqivisualization.= true; simParameters.RBVisualization = false;％输出指标绘图在周期性的NumMetricsSteps中更新％模拟持续时间simParameters。NumMetricsSteps = 20;要写入日志的% mat文件。它们用于后期仿真分析和可视化simparameters.ParametersLogfile =.'simparameters'；用于记录模拟参数的％simparameters.simulationlogfile.=“simulationLogs”；用于记录模拟日志的％

%设置终端定时DL应用流量模型。以下%配置适用于每个单元格simparameters.dlpacketperiodicitues = 10;%在gNB上终端生成DL包的周期性(毫秒)simparameters.dlpacketsize = 2e4;在GNB处为UE生成的DL数据包的％大小（以字节为单位）％验证模拟配置hNRIntercellInterferenceValidateConfig (simParameters);

simparameters.duplexmode = 0;％FDD.simparameters.ulbandwidth = 10e6;％Hz.simParameters。ULCarrierFreq = 2.515 e9;％Hz.simParameters。NumUEs = simParameters.NumUEsCell;%单元内UEs的个数simparameters.ncellidlist = 0：SimParameters.Numcells-1;%物理单元id列表％CSI-RS资源配置。假设所有UE都在相同的CSI-RS资源上测量信道质量simParameters。CSIRSRowNumber = 2;%单个发射天线的可能行号为1和2simparameters.subbandsize = 8;%以RBs数量表示的CQI报告子频带的大小％将bsrperiodicity设置为“inf”，因为没有UL流量simParameters。BSRPeriodicity =正;百分比所选SCS的％插槽持续时间和10 ms帧中的插槽数slotDuration = 1 / (simParameters.SCS / 15);百分比numslotsframe = 10 / slotduration;10 ms帧中的次数数numslotssim = simparames.numframessim * numslotsframe;%模拟中的槽数指标可视化更新的时间间隔％插槽。由于一个插槽是模拟的最佳时间粒度，使%确定MetricsStepSize是一个整数simparameters.metricsstepsize = ceil（numslotssim / simparameters.nummetricssteps）;如果mod(numSlotsSim, simParameters.NumMetricsSteps) ~= 0如果numslotssim不是，％更新nummetricssteps参数%完全被它整除simparameters.nummetricssteps =楼层（numslotssim / simparameterm.metricsstepsize）;结束numLogicalChannels = 1;在此示例中，每个UE中仅在每个UE中假设1个逻辑信道％逻辑信道配置适用于模拟中的所有节点（UE和GNB）simParameters.LCHConfig.LCID = 4;％RLC实体方向。值0仅表示DL，1%只代表UL, 2代表UL和DL%的方向。设置实体方向同时具有UL和DLsimParameters.RLCConfig.EntityDir = 0;％创建RLC通道配置结构rlcchannelconfigstruct.lcgid = 1;%逻辑通道与逻辑通道组ID的映射关系rlcChannelConfigStruct。优先级= 1;每个逻辑信道的％优先级rlcchannelconfigstruct.pbr = 8;%每个逻辑通道的优先比特率(PBR)，单位为千字节/秒rlcChannelConfigStruct。BSD = 10;％桶尺寸持续时间（BSD），在MS中，每个逻辑通道rlcChannelConfigStruct。EntityType = simParameters.RLCConfig.EntityDir;rlcChannelConfigStruct。LogicalChannelID = simParameters.LCHConfig.LCID;每3GPP TS 38.323，％最大RLC服务数据单元（SDU）长度（以字节为单位）simparameters.maxrlcsdulength = 9000;％在每个单元格中生成UE的位置simparameters.ueposition = generateuepositions（simparameters）;模拟中的UE总数simParameters。MaxReceivers = simParameters。NumCells* simParameters.NumUEsCell;如果~ isfield (simParameters'schedulingtype'）||simparameters.schedulingtype == 0.%不指定调度类型或指定基于槽位的调度rbAssignmentPlotPeriodicity = numSlotsFrame;%每帧更新RB分配可视化(10毫秒)Tickgranularity = 14;simparameters.puschmappeType ='一种'；simparameters.pdschmapptype =.'一种'；其他的基于％符号的调度rbAssignmentPlotPeriodicity = 1;%更新RB分配可视化每个插槽tickgranularity = 1;simparameters.puschmappeType ='B'；simparameters.pdschmapptype =.'B'；结束

gNB =细胞(simParameters。NumCells, 1);问题=细胞(simParameters。NumCells那simParameters。NumUEsCell);％创建DL数据包分发对象dlpacketdistributionobj = hnrpacketdistribution（simparameters，0）;％0对于DL％创建UL数据包分发对象ulPacketDistributionObj = hNRPacketDistribution(simParameters, 1);UL的％1为cellIdx = 1: simParameters。NumCells每个单元格%simparameters.ncellid = simparameters.ncellidlist（cellidx）;％单元IDsimparameters.position = [simparameters.gnbposition（cellidx，:) 0];% gNB在(x,y,z)坐标中的位置gnb {cellidx} = hnrgnb（simparameters）;%创建gNB节点Scheduler = HNRSchedulerProportionalFair（SimParameters）;%创建比例公平调度程序Addscheduler（GNB {CellidX}，调度程序）;％为GNB添加调度程序gnb {cellidx} .phyentity = hnrgnbphy（simparameters）;％创建PHY图层实例configurephy（gnb {cellidx}，simparameters）;%配置PHY层setphyInterface（gnb {cellidx}）;％设置接口到PHY层每个单元格的百分比，创建一组UE节点并随机将它们放在小区半径内为UEIDX = 1：SimParameters.Numuescell SimParameters.Position = [SimParameters.ueposition {CellidX}（UEIDX，:) 0];UE中的％位置（x，y，z）坐标UES {CellidX，UEIDX} = HNRUE（SIMPParameters，UEIDX）;UES {CellidX，UEIDX} .Phyentity = HNRuephy（SimParameters，UEIDX）;％创建PHY图层实例configurePhy(问题{cellIdx, ueIdx}, simParameters);%配置PHY层setphyInterface（UE {cellidx，UEIDX}）;%设置接口到PHY%为终端在gNB设置逻辑通道configurelogicalChannel（GNB {CellidX}，UEIDX，RLCChannelConfigstruct）;UE的％设置逻辑通道configurelogicalChannel（UE {Cellidx，UEIDX}，UEIDX，RLCChannelConfigstruct）;%向gNB节点添加数据流量模式生成器packetsize = simparameters.dlpacketsize;%使用命令计算开-关流量模式的数据速率(kbps)％数据包大小（以字节为单位）和数据包间隔（在MS中）dataRate = cell (1000/simParameters.DLPacketPeriodicityUEs) * packetSize * 8e-3;％将生成的应用程序包的大小限制为最大值％RLC SDU大小。最大支持的RLC SDU万博1manbetx大小为9000字节如果数据包> SimParameters.maxrlcsdulength packetsize = simparameters.maxrlcsdulength;结束%为指定的开关网络流量模式创建一个对象％UE并将其添加到GNB中。该对象生成下行链路数据为UE的GNB上的％流量app = networktrafficonoff（'包裹'packetSize,“GeneratePacket”， 真的，......“定时”, simParameters。NumFramesSim/100,'反转'0,“DataRate”， 数据速率）;gnb {cellidx} .addapplication（ueidx，simparameters.lchconfig.lcid，app）;结束%设置UL和DL分组分发机制hNRSetUpPacketDistribution(simParameters, gNB{cellIdx}， UEs(cellIdx，:)， dlPacketDistributionObj, ulPacketDistributionObj);结束

％显示网络拓扑plotNetwork (simParameters);

simSchedulingLogger =细胞(simParameters。NumCells, 1);simPhyLogger =细胞(simParameters。NumCells, 1);为cellIdx = 1: simParameters。NumCells%创建用于MAC DL调度信息可视化和日志记录的对象simparameters.ncellid = simparameters.ncellidlist（cellidx）;simschedulinglogger {cellidx} = hnrschedulinglogger（simparameters，0）;％0对于DL为PHY层度量日志创建一个对象simPhyLogger{cellIdx} = hNRPhyLogger(simParameters, 0);％0对于DL结束存储感兴趣的单元格ID的索引cellOfInterestIdx =找到(simParameters。CellOfInterest = = simParameters.NCellIDList);%为MAC和PHY指标创建指标可视化对象(DL方向表示为0)视觉型的人= hNRMetricsVisualizer (simParameters,“MACLogger”, simSchedulingLogger {cellOfInterestIdx},“PhyLogger”, simPhyLogger {cellOfInterestIdx},'VisualizationFlag',0);％运行处理循环slotnum = 0;numsymbolssim = numslotssim * 14;以符号持续时间为单位的仿真时间%执行模拟中的所有符号为symbolnum = 1：划分：numsymbolssim如果mod(symbolNum - 1,14) == 0 slotNum = slotNum + 1;结束％SCS在相同SC上操作的所有单元都相同为cellIdx = 1: simParameters。NumCells每个单元格%%运行gNB的MAC和PHY层运行（gnb {cellidx}）;%运行ue的MAC层和PHY层为UEIDX = 1：SimParameters.Numuescell Run（UE {CellidX，UEIDX}）;结束% MAC日志logcellschedulingstats（simschedulinglogger {cellidx}，symbornum，gnb {cellidx}，ues（cellidx，:)，0）;％0对于DL％phy logging.logcellphystats（simphylogger {cellidx}，symbornum，gnb {cellidx}，UE（cellidx，:)）;结束%的可视化％检查插槽边界如果symbolnum> 1 &&（（simparameters.schedulingtype == 1 && mod（symbolnum，14）== 0）||（simparameters.schedulingtype == 0 && mod（symbolnum-1,14）== 0））% RB赋值可视化(如果启用)如果simParameters。RBVisualization如果mod(slotNum, rbAssignmentPlotPeriodicity) == 0%如果达到更新周期，在插槽边界绘制plotrbgrids（simschedulinglogger {cellofInterestidx}）;结束结束% CQI网格可视化(如果启用)如果simparameters.cqivisualization.如果mod（slotnum，numslotsframe）== 0帧边界的％绘图plotcqirbrids（simschedulinglogger {cellofInterestidx}）;结束结束%如果达到更新周期，绘制调度器指标和PHY指标可视化槽边%如果mod（slotnum，simparameters.metricsstepsize）== 0 plotmetrics（可视化器，slotnum）;结束结束为cellIdx = 1: simParameters。NumCells％autain定时器为gnb和ues的刻度按每个插槽的符号数advanceTimer (gNB {cellIdx}, tickGranularity);为ueIdx = 1: simParameters。NumUEsCell advanceTimer(UEs{cellIdx, ueIdx}, tickGranularity);结束结束结束

dlStats = getPerformanceIndicators (simSchedulingLogger {simParameters。CellOfInterest + 1});[~， avgBLERLogs] = getBLERLogs(simPhyLogger;CellOfInterest + 1});fprintf（'\nPeak DL吞吐量:%0.2f Mbps。达到的平均DL吞吐量:%0.2f Mbps'， dlStats(1,1)， dlStats(2,1);

fprintf（'\nAchieved average DL Goodput: %0.2f Mbps'，dlstats（5,1））;

fprintf（'\nPeak DL频谱效率:%0.2f bits/s/Hz。达到平均DL频谱效率:%0.2f bits/s/Hz'，dlstats（3,1），dlstats（4,1））;

DISP（['下行链路方向上每个UE的阻止错误率：['num2str（round（avgblerlogs（:, 1）'，2）））“]”]);

％获取日志如果(simParameters。DuplexMode = = 0)％FDD.loginfo = struct（“NCellID”[],“DLTimeStepLogs”[],'schedulingassignmentlogs'[],“BLERLogs”[],“AvgBLERLogs”，[]）;其他的loginfo = struct（“NCellID”[],'timesteplogs'[],'schedulingassignmentlogs'[],“BLERLogs”[],“AvgBLERLogs”，[]）;结束simulationlogs = cell（simparameters.numcells，1）;为cellidx = 1：simparameters.numcells loginfo.ncellid = simparameters.ncellidlist（cellidx）;如果(simParameters。DuplexMode = = 0)％FDD.logInfo。DLTimeStepLogs = getSchedulingLogs (simSchedulingLogger {cellIdx});其他的％TDD.loginfo.timesteplogs = getschedulinglogs（simschedulinglogger {cellidx}）;结束loginfo.schedulingassignmentlogs = getgrantlogs（simschedulinglogger {cellidx}）;计划分配日志loginfo.blerlogs = getBlerlogs（SimphyLogger {Cellidx}）;％块错误率日志simulationLogs{cellIdx, 1} = logInfo;结束保存(simParameters。ParametersLogFile,'simparameters'）;％在MAT文件中保存模拟参数保存（SimParameters.simulationLogfile，“simulationLogs”）;%将模拟日志保存到mat -文件中

函数Plotnetwork（SimParameters）%创建数字图('名称'那'网络拓扑可视化'那“单位”那'标准化'那'offormosition'，[0 0 1 1]，'可见的'那“在”）;标题('网络拓扑可视化'）;持有在；为cellIdx = 1: simParameters。NumCells％绘制圆圈th = 0:π/ 60:2 *π;xunit = simParameters。CellRadius * cos(th) + simParameters。GNBPosition (cellIdx, 1);yunit = simParameters。CellRadius * sin(th) + simParameters。GNBPosition (cellIdx 2);如果simparameters.cellofinterest == simparameters.ncellidlist（cellidx）h1 = plot（xunit，yunit，'颜色'那“绿色”）;感兴趣细胞%其他的H2 = plot(xunit, yunit，)'颜色'那'红色的'）;结束Xlabel（'X位置（米）') ylabel (的坐标(米)）%为gnb添加刀尖数据S1 =分散（SimParameters.gnboposition（cellidx，1），simparameters.gnboposition（cellidx，2），“^”那'markeredgecolor'那'品红'）;cellIdRow = dataTipTextRow (“细胞- - - - - -”，{num2str（simpparameters.ncellidlist（cellidx））}）;s1.datatiptemplate.datatiprows（1）= cellidrow;posRow = dataTipTextRow ('位置[x，y]：'，{[“(”num2str（simparameters.gnbposition（cellidx，:)）“]”]});s1.datatiptemplate.datatiprows（2）= posrow;%为终端添加工具提示数据Ueposition = simparameters.ueposition {cellidx};为UEIDX = 1：尺寸（UEPOSINE，1）S2 =分散（UEEDX，1），UEEPORION（UEIDX，2），'。'那'markeredgecolor'那'蓝色的'）;ueIdRow = dataTipTextRow (“问题- - - - - -”{num2str (ueIdx)});s2.DataTipTemplate.DataTipRows (1) = ueIdRow;posRow = dataTipTextRow ('位置[x，y]：'，{[“(”num2str (uePosition (ueIdx:))“]”]});s2.DataTipTemplate.DataTipRows (2) = posRow;结束结束％创造了传说如果Simparameters.numcells> 1个图例（[H1 H2 S1 S2]，“感兴趣的细胞”那干扰细胞的那“gNodeB”那“问题”那“位置”那“northeastoutside”）其他的图例（[H1 S1 S2]，“感兴趣的细胞”那“gNodeB”那“问题”那“位置”那“northeastoutside”）结束轴([0 4000 0 3000]);设置轴限持有从；daspect([1000、1000、1]);%设置数据宽高比结束函数uePositions = generateUEPositions (simParameters)％返回每个单元中UE的位置Uepositions = cell（simparameters.numcells，1）;为cellidx = 1：simparameters.numcells gnbxco = simparameters.gnbposition（cellidx，1）;% gNB坐标gnbyco = simparameters.gnbosition（cellidx，2）;% gNB坐标Theta = Rand（SimParameters.numuescell，1）*（2 * Pi）;％表达以计算细胞中UE的位置。默认情况下，%它将在单元格中随机放置uer =√兰德(simParameters。NumUEsCell, 1))*simParameters.CellRadius; x = round(gnbXCo + r.*cos(theta)); y = round(gnbYCo + r.*sin(theta)); uePositions{cellIdx} = [x y];结束结束