simParameters =结构（）;％Clear SimParameters变量包含所有密钥仿真参数simparameters.nframes = 2;% 10毫秒帧数simParameters。SNRIn = [-5 0 5];％SNR范围（DB）

simParameters。PerfectChannelEstimator = true;

simparameters.displaydiagnostics = false;

％SET波形类型和PDSCH数字（SCS和CP类型）simparameters.carrier = nrcarrierconfig;simparameters.carrier.nsizegrid = 51;资源块数的％带宽（51 rbs，30 kHz SCS为20 MHz BW）simparameters.carrier.subcarrierspacing = 30;％15,30,60,120,240（kHz）simparameters.carrier.CyclicPrefix =.'普通的';%“正常”或“扩展”(扩展CP只适用于60千赫SCS)simParameters.Carrier.NCellID = 1;%细胞身份％SS突发配置可以通过将ssbtransmited字段设置为所有零来禁用％突发simParameters.SSBurst =结构（）;simParameters.SSBurst.BlockPattern ='案例b';% 30千赫副载波间距simparameters.ssburst.ssbtransmited = [0 1 0 1];％位图指示在突发中传输的块simparameters.ssburst.ssbperiodicity = 20;％SS突发在MS（5,10,20,40,80,160）中的周期性％PDSCH / DL-SCH参数simParameters。PDSCH = nrPDSCHConfig;％此PDSCH定义是BLER仿真中所有PDSCH传输的基础simParameters。PDSCHExtension =结构();％该结构是为了保持DL-SCH和PDSCH的额外仿真参数％定义每个时隙PDSCH的时间 - 频率资源分配是满格（单满格BWP）simParameters.PDSCH.PRBSet = 0: simParameters.Carrier.NSizeGrid-1;% PDSCH PRB分配simParameters.PDSCH.SymbolAllocation = [0, simParameters.Carrier.SymbolsPerSlot];每个PDSCH分配的％启动符号和符号数simParameters.PDSCH.MappingType ='一种';％PDSCH映射类型（'a'（插槽），'b'（非插槽））％加扰标识符simparameters.pdsch.nid = simparameters.carier.ncellid;simparameters.pdsch.rnti = 1;％PDSCH资源块映射（TS 38.211第7.3.1.6节）simParameters.PDSCH.VRBToPRBInterleaving = 0;％禁用交错资源映射simparameters.pdsch.vrbbundleSize = 4;％定义将要使用的传输层的数量simparameters.pdsch.numlayers = 2;% PDSCH传输层数定义码字调制和目标编码率％代码字的数量直接依赖于图层的数量，以确保在获取码字编号之前首先设置％图层如果Simparameters.pdsch.numcodewords> 1％Multicodeword传输（当层的数目为> 4）simParameters.PDSCH.Modulation = {16 qam的那16 qam的};％'QPSK'，'16QAM'，'64QAM'，'256QAM'simParameters.PDSCHExtension.TargetCodeRate = [490490]/1024;%用于计算传输块大小的码率别的simparameters.pdsch.modulation =.16 qam的;％'QPSK'，'16QAM'，'64QAM'，'256QAM'simparameters.pdschextension.targetCoderate = 490/1024;%用于计算传输块大小的码率结束% DM-RS和天线端口配置(TS 38.211章节7.4.1.1)simparameters.pdsch.dmrs.dmrsportset = 0：simparameters.pdsch.numlayers-1;用于各层的DM-RS端口simparameters.pdsch.dmrs.dmrstypeaposition = 2;%仅映射类型A。第一个DM-RS符号位置(2,3)simparameters.pdsch.dmrs.dmrslength = 1;%前置DM-RS符号数(1个(单符号)，2个(双符号))simParameters.PDSCH.DMRS.DMRSAdditionalPosition = 0;％附加DM-RS的符号位置（最大范围0 ... 3）simparameters.pdsch.dmrs.dmrsconfiguringtype = 2;% DM-RS配置类型(1,2)simParameters.PDSCH.DMRS.NumCDMGroupsWithoutData = 1;%无数据的CDM组数simParameters.PDSCH.DMRS.NIDNSCID = 1;％扰频标识（0 ... 65535）simparameters.pdsch.dmrs.nscid = 0;％加扰初始化（0,1）％PT-RS配置（TS 38.211第7.4.1.2节）simParameters.PDSCH.EnablePTRS = 0;％启用或禁用PT-RS（1或0）simParameters.PDSCH.PTRS.TimeDensity = 1;％pt-rs时间密度（l_pt-rs）（1,2,4）simparameters.pdsch.ptrs.frequencyDensy = 2;PT-RS频率密度百分比(K_PT-RS)(2或4)simparameters.pdsch.ptrs.reoffset =.'00';％Pt-RS资源元素偏移（'00'，'01'，'10'，'11'）simparameters.pdsch.ptrs.ptrsportset = [];％PT-RS天线端口，DM-RS端口集的子集。空对应于较低的DM-RS的端口号％保留PRB图案，如果需要的话（对于CORESETs，前向兼容性等）simParameters.PDSCH.ReservedPRB{1}。SymbolSet = [];％保留的PDSCH符号simparameters.pdsch.reservedProb {1} .prset = [];%保留PDSCH PRBssimparameters.pdsch.reservedprb {1} .period = [];保留资源的百分比％额外的模拟和DL-SCH相关参数％％PDSCH PRB捆绑（TS 38.214第5.1.2.3节）simParameters.PDSCHExtension.PRGBundleSize = [];％2,4或[]表示“宽带”％HARQ过程和速率匹配/ TBS参数simParameters.PDSCHExtension.XOverhead = 6 * simParameters.PDSCH.EnablePTRS;当启用PT-RS时，％设置TBS（XOH）到6的顶部开销匹配的PDSCH速率匹配，否则为0simparameters.pdschextension.nharqprocesses = 16;%要使用的并行HARQ进程数simparameters.pdschextension.enableharq = true;％使用RV序列来实现每个进程的重传[0,2,3,1]LDPC解码器参数％可用算法：“信仰传播”，“分层信仰传播”，“归一化最小和”，“偏移最小和”simParameters.PDSCHExtension.LDPCDecodingAlgorithm =“分层信念传播”;simparameters.pdschextension.maximumLDPCiterationCount = 6;％在终点处定义整个传输天线几何％如果使用CDL传播信道，则整数的天线元件数量是%转换为天线面板配置时，通道模型对象被创建simparameters.ntxants = 8;PDSCH传输天线数量（1,2,4,8,16,32,64,128,256,512,1024）> = Numlayers如果Simparameters.pdsch.numcodewords> 1％多码字传输simParameters。NRxAnts = 8;％UE的数量接收天线（偶数> = numlayers）别的simParameters。NRxAnts = 2;% UE接收天线数量(1或偶数>= NumLayers)结束%定义通用CDL/TDL传播通道参数simparameters.delayprofile =.“CDL-C”;使用CDL-C模型(Urban macrocell model)simparameters.delayspread = 300e-9;Simparameter.maximumdopplershift = 5;％交叉检查PDSCH分层对沟道几何形状validateNumLayers (simParameters);

waveforminfo = nrofdminfo（simparameters.carrier）;%获取OFDM调制步骤后的基带波形信息

％构成的CDL或TDL信道模型对象如果包含（simparameters.delayprofile，'CDL'那'IGNORECASE'，true）channel = nrcdlchannel;％CDL信道对象％将天线的总数转换为特定的天线面板％阵列几何形状。配置的天线数量更新% nTxAnts不是(1、2、4、8、16、32、64、128、256、512、1024)或nRxAnts中的一个％不是1甚至1。[channel.TransmitAntennaArray.Size，channel.ReceiveAntennaArray.Size] =......harraygeometry（simparameters.ntxants，simparameters.nrxants）;ntxants = prod（channel.transmitantennaarray.size）;nrxants = prod（channel.receiveantenntennaray.size）;simparameters.ntxants = ntxants;simparameters.nrxants = nrxants;别的channel = nrtdlchannel;％TDL通道对象％设置频道几何形状通道。NumTransmitAntennas = simParameters.NTxAnts;通道。NumReceiveAntennas = simParameters.NRxAnts;结束％指定模拟信道参数和波形采样率与对象channel.delayprofile = simparameters.delayprofile;channel.delayspread = simparameters.delayspread;Channel.maximumdopplershift = SimParameters.maximumdopplershift;Channel.Samplerve = WaveformInfo.Sampleate;

Chinfo = Info（频道）;maxchdelay = ceil（max（chinfo.pathdelays * channel.samplerate））+ chinfo.channelfilterdelay;

％获取信息有关SS突发配置首先需要一些相关的参数赋值simparameters.ssburst.ncellid = simparameters.carrier.ncellid;simparameters.ssburst.samplerve = waveforminfo.sampleate;ssbinfo =赫斯伯斯坦福（Simparameters.ssburst）;％映射占用的子载波和SS突发的传输符号％（以SS定义突发命理）到PDSCH的PRB和符号在％PDSCH BWP /载体命理[mappedPRB, mappedSymbols] = mapNumerology (ssbInfo.OccupiedSubcarriers、ssbInfo.OccupiedSymbols ssbInfo.NRB, simParameters.Carrier.NSizeGrid, ssbInfo.SubcarrierSpacing, simParameters.Carrier.SubcarrierSpacing);％配置PDSCH以保护这些资源，以便PDSCH%传输不重叠的SS爆发预约= nrpdschreservedconfig;预留.symbolset = mappedsymbols;预约.prset = mappedprog;保留.period = simparameters.ssburst.ssbperifyicity *（simparameters.carrier.subcarrierspacing / 15）;%插槽的周期simparameters.pdsch.reservedprb {end + 1} =预留;

％数组存储所有SNR点的最大吞吐量maxthrougpul = zeros（长度（simparameters.snrin），1）;％数组存储所有SNR点的模拟吞吐量simThroughput = 0(长度(simParameters.SNRIn), 1);%根据HARQ配置，设置要使用的冗余版本(RV)序列如果simParameters.PDSCHExtension.EnableHARQ%在RAN WG1 #91会议(R1-1719301)的最终报告中，它是在R1-1717405中观察到的％，如果性能是优先级，[0 2 3 1]应该使用%。如果自解码是优先级，那么它应该是％考虑到码率的上限%每个RV自解码的顺序如下:0>3>2>1rvSeq = [0 2 3 1];别的％HARQ禁用 - 具有RV = 0的单个传输，无重传rvseq = 0;结束％创建DL-SCH编码器系统对象以执行传输通道编码EncodedLsch = nrdlsch;EncodedLsch.multipreharqprocesses = true;EncodedLsch.targetCoderate = SimParameters.PDSchextension.targetCoderate;％创建DL-SCH解码器系统对象以执行传输信道解码％使用分层信仰传播为LDPC解码，有一半的数量％迭代相比的默认置信传播解码DecodedLsch = nrdlschdecoder;DecodedLsch.multipreharqprocesses = true;DecodedLsch.targetCoderate = SimParameters.pdschextension.targetCoderate;DecodedLsch.ldpcdecodingAlgorithm = SimParameters.PDSchextension.LDPCDecodingAlgorithm;解码ledlsch.maximumldpciterationcount = simparameters.pdschextension.maximumLdpciterationCount;为了snridx = 1：numel（simparameters.snrin）%注释掉用于并行计算％parcor snridx = 1：numel（simpparameters.snrin）％解释并行计算％要减少总模拟时间，可以执行此循环% parallel使用并行计算工具箱。注释掉for%语句并取消“parfor”语句的注释。如果并行计算未安装％工具箱，'Parcon'默认为正常'for'语句%将随机数生成器设置为默认值RNG（“默认”);％采取仿真级参数结构的完整副本，以便它们不是使用Parcom时％PCT广播变量simlocal = simparameters;WaveInfolocal = WaveformInfo;％拍摄频道级参数的副本以简单地进行后续参数引用carrier = simlocal.carrier;pdsch = simlocal.pdsch;pdschextra = simlocal.pdschextension;ssburst = simlocal.ssburst;DecodedLschlocal =解码ledlsch;%复制的解码器句柄，以帮助PCT分类变量decodedlschlocal.reset（）;%在每个信噪比点开始时重置译码器pathFilters = [];ssbWaveform = [];%为新的信噪比点准备模拟SNRdB = simLocal.SNRIn (snrIdx);fprintf中（'\n模拟传输方案1 (%dx%d)和SCS=%dkHz， %s通道在%gdB信噪比为%d 10ms帧(s)\n'那......simparameters.ntxants，simparameters.nrxants，carrier.subcarrierspacing，......simlocal.delayprofile，snrdb，simlocal.nframes）;在模拟和分析中使用％初始化变量bittput = [];每次传输成功收到的比特数量txedTrBlkSizes = [];%每次传输的传输信息比特数％指定的顺序周期，我们通过HARQ进程harqsequence = 1：pdschextra.nharqprocesses;％初始化所有HARQ进程的状态HarqProcesses = HNewHarqProcesses（PDSchextra.nharqProcesses，Rvseq，PDSCH.Numcodewords）;harqproccntr = 0;％HARQ过程计数器％重置通道，使每个SNR点会遇到相同的%通道实现复位（信道）;仿真期间中的斜槽总数nslots = simlocal.nframes * carrier.slotsperframe;％索引到当前SS突发样本的开始的开始传送％ssbsampleindex = 1;％获得用于传输的预编码矩阵（WTX）第一个传输块estChannelGrid = getInitialChannelEstimate(载体、simLocal.NTxAnts通道);newWtx = getPrecodingMatrix(载体,pdsch estChannelGrid);％定时偏置，更新在每个时隙为完美的同步和当相关性强直时％是实际同步时的％偏移量= 0;％波形长度的％循环为了nslot = 0: NSlots-1％更新新插槽的运营商插槽号母舰。NSlot = NSlot;在必要时生成新的SS突发如果（ssbsampleindex == 1）nsubframe = carrier.nslot / carrier.slotspersubframe;ssburst.nframe =楼层（nsubframe / 10）;ssburst.nhalfframe = mod（nsubframe / 5,2）;[SSBWAVEFORM，〜，SSBINFO] =奶族（SSBURST）;结束从HARQ索引表中获取当前PDSCH的HARQ进程索引harqprocidx = harqsequence（mod（harqproccntr，长度（harqsequence））+ 1）;％更新当前HARQ进程信息（此更新RV%取决于CRC通过或失败在以前的传输％此HARQ过程）harqProcesses (harqProcIdx) = hUpdateHARQProcess (harqProcesses (harqProcIdx) pdsch.NumCodewords);％计算插槽中的码字的传输块大小[pdschIndices，pdschIndicesInfo] = nrPDSCHIndices（载体，PDSCH）;trBlkSizes = nrTBS（pdsch.Modulation，pdsch.NumLayers，numel（pdsch.PRBSet），pdschIndicesInfo.NREPerPRB，pdschextra.TargetCodeRate，pdschextra.XOverhead）;％HARQ处理％检查每个码字的先前传输的CRC，即需要重传？为了cwIdx = 1: pdsch。NumCodewords newdata = false;如果HarqProcesses（HarqProcidx）.blkerr（cwidx）最后录制解码的％错误如果(harqProcesses (harqProcIdx) .RVIdx (cwIdx) = = 1)％发出RV序列的开始ResetsoftBuffer（解码型雪地，CWIDX-1，HarqProcidx-1）;在这种情况下需要％显式重置newdata = true;结束别的%没有错误newdata = true;结束如果newdata trBlk = randi([0 1]，trBlkSizes(cwIdx)，1);setTransportBlock (encodeDLSCH trBlk、cwIdx-1 harqProcIdx-1);结束结束%编码DL-SCH传输块codedTrBlocks = encodeDLSCH（pdsch.Modulation，pdsch.NumLayers，......pdschIndicesInfo.G harqProcesses (harqProcIdx) .RV harqProcIdx-1);得到在前一个槽中计算的预编码矩阵(wtx)wtx = newwtx;％资源网格阵列pdschGrid = nrResourceGrid(载体、simLocal.NTxAnts);% PDSCH调制和预编码pdschSymbols = nrPDSCH(载体,pdsch codedTrBlocks);[pdschAntSymbols, pdschAntIndices] = hPRGPrecode(大小(pdschGrid)、carrier.NStartGrid pdschSymbols, pdschIndices, mtx);与PDSCH传输周期相关联的电网中的PDSCH映射%pdschGrid (pdschAntIndices) = pdschAntSymbols;％PDSCH DM-RS预编码和映射DMRSSYMBOLS = NRPDSCHDMRS（载波，PDSCH）;dmrsindices = nrpdschdmrsindices（载体，pdsch）;[dmrsantsymbols，dmrsantindices] = hprgprecode（size（pdschgrid），carrier.nstartgrid，dmrssymbols，dmrsindices，wtx）;pdschgrid（dmrsantindices）= dmrsantsymbols;％PDSCH PT-RS预编码和映射ptrssymbols = nrpdschptrs（载体，pdsch）;ptrsindices = nrpdschptrsindices（载体，pdsch）;[ptrsantsymbols，ptrsantindices] = hprgprecode（大小（pdschgrid），carrier.nstartgrid，ptrssymbols，ptrsindices，wtx）;pdschgrid（ptrsantindices）= ptrsantsymbols;关联资源元素的％OFDM调制txwaveform = nrofdmodulate（载体，pdschgrid）;％将SS突发波形的适当部分添加到％传输波形元=大小(txWaveform, 1);txWaveform = txWaveform + ssbWaveform(ssbSampleIndex + (0:Nt-1)，:);ssbSampleIndex = mod(ssbSampleIndex + Nt,size(ssb波形，1));％通过频道模型传递数据。在末尾附加零%发送波形以刷新通道内容。这些零百分比频道引入的任何延迟。这是一个混合多路径延迟和实现延迟的百分比。这个值可能取决于采样率，延迟配置文件和延迟的百分比变更％ 传播txWaveform = [txWaveform;0 (maxChDelay、尺寸(txWaveform 2)));[rxWaveform, pathGains sampleTimes] =通道(txWaveform);%将AWGN添加到接收的时域波形中通过在OFDM调制中使用的IFFT大小％规格化噪声功率，随着OFDM调制器应用此归一化的％%传播波形。也要按照接收的数量进行规范%天线，因为信道模型将这种标准化应用于默认情况下％接收波形信噪比= 10 ^ (SNRdB / 20);％计算线性噪声增益n0 = 1 /（sqrt（2.0 * simlocal.nrxants * double（waveinfolocal.nfft））* snr）;噪声= N0 *复合物（Randn（randn（rxwaveform）），randn（size（rxwaveform）））;rxwaveform = rxwaveform +噪声;如果（Simlocal.PerfectChannEleStimator）％完美同步。使用频道提供的信息％找到最强的多径组件pathfilters = getpathfilters（频道）;得到完美的信道估计路径滤波器[offset，mag] = nrperfecttimingestimate（pathgains，pathfilters）;别的％实际同步。将接收的波形相关联%与PDSCH DM-RS给出时间偏移估计't'和％相关幅度'mag'。函数hskipweaktimingoffset%用于更新接收机定时偏差。如果%相关峰值在'mag'是弱的，当前的时机%估计值't'被忽略，之前的估计值'offset'被忽略使用%[t，mag] = nrtimingestimate（运营商，rxwaveform，dmrsindices，dmrssymbols）;Offset = HskipWeakTimingOffsoffs（偏移，T，MAG）;％显示警告，如果估计的定时偏移量超过％最大信道延迟如果offset> maxchdelay警告（[估计的时间偏移(%d)大于最大信道延迟(%d)。......这将导致解码失败。这可能是低信噪比造成的。”......'或者没有足够的DM-RS符号来成功同步。']，偏移，maxchdelay）;结束结束rxwaveform = rxwaveform（1 + offset：结束，:);％在接收数据上执行OFDM解调以重新创建%资源网格，包括实际事件中的填充%同步导致一个不完整的槽被解调rxGrid = nrOFDMDemodulate（载体，rxWaveform）;[K，L，R] =尺寸（rxGrid）;如果（L 结束如果（Simlocal.PerfectChannEleStimator）％完美的通道估计，使用路径增益的值％由渠道提供。该频道估计没有％包括发射器预编码的效果estChannelGrid = nrPerfectChannelEstimate(载体,pathGains pathFilters,抵消,sampleTimes);％获得完美的噪音估计（来自噪声实现）kotidgrid = nrofdmdemodulate（载波，噪声（1 +偏移：结束，:)）;destemest = var（诊断（:)）;％获取预编码矩阵用于下一时隙newWtx = getPrecodingMatrix(载体,pdsch estChannelGrid);％获取从接收到的电网PDSCH资源要素和%信道估计[pdschRx pdschHest, ~, pdschHestIndices] = nrExtractResources (pdschIndices, rxGrid estChannelGrid);％将预编码应用于频道估计pdschhest = hprgprecode（size（estchannelgrid），carrier.nstartgrid，pdschhest，pdschhestindices，perfute（wtx，[2 1 3]））;别的所接收网格与百分比实际信道估计%每个传输层，使用PDSCH DM-RS为每个％ 层。该频道估计包括效果％变送器预编码[EstchannElgrid，Desement] = NrchanneLestimate（运营商，RxGrid，DMRSindices，DMRSSymbols，'cdmlength', pdsch.DMRS.CDMLengths);％获取从接收到的电网PDSCH资源要素和%信道估计[pdschrx，pdschhest] = nrextractresources（pdschindices，rxgrid，estchannelgrid）;％在预编码之前从EstchannelGrid中删除预编码％矩阵计算estChannelGridPorts = precodeChannelEstimate(载体、estChannelGrid连词(mtx));％获取预编码矩阵用于下一时隙newwtx = getPrecodingMatrix（运营商，PDSCH，EstchannelGridports）;结束％均衡[PDSCHEQ，CSI] = NREQUALIZEMSE（PDSCHRX，PDSCHHEST，噪音）;％常见相位误差（CPE）补偿如果〜isempty（ptrsindices）初始化临时网格以存储相等的符号tempGrid = nrResourceGrid(载体、pdsch.NumLayers);%从接收的网格中提取PT-RS符号并进行估计%通道网格[ptrsrx，ptrshest，〜，〜，ptrshestindices，ptrslayerindices] = nrextractresources（ptrsindices，rxgrid，estchannelgrid，tempgrid）;如果（Simlocal.PerfectChannEleStimator）％将预编码应用于频道估计ptrshest = hprgprecode（size（estchannelgrid），carrier.nstartgrid，ptrshest，ptrshestindices，erfute（wtx，[2 1 3]））;结束%均衡PT-RS符号并将它们映射到tempGridptresq = nrequalizemse（ptrsrx，ptrshest，噪音）;tempgrid（ptrslayerindices）= ptreqe;％估计PT-RS位置处的残差通道％tempgrid.CPE = nrchannelestimate（tempgrid，ptrsindices，ptrssymbols）;在子载波，接收天线和%层。然后，通过取这个角得到CPE%的总和CPE =角度（SUM（CPE，[1 3 4]））;%将均衡的PDSCH符号映射到tempGridtempgrid（pdschindices）= pdscheq;在参考范围内的每个OFDM符号中的％正确CPE％pt-rs OFDM符号symloc = pdschindicsInfo.ptrssymbolset（1）+1：pdschindicsInfo.ptrssymbolset（END）+1;tempgrid（：，symloc，:) = tempgrid（：symloc，:)。* exp（-1i * cpe（symloc））;%提取PDSCH符号pdscheq = tempgrid（pdschindices）;结束％解码PDSCH物理通道[dlschlls，rxsymbols] = nrpdschdecode（载波，pdsch，pdscheq，噪声）;％每层显示EVM，每个插槽和每个RB如果（SIMLOCAL.DISPLYDIAGNOSTICS）PLOTLAYEREVM（NSLOTS，NSLOT，PDSCH，尺寸（PDSCHGRID），PDSCHINDICE，PDSCHSYMBOLS，PDSCHEQ）;结束CSI％SCALE LLRScsi = nrlayerdemap（csi）;％CSI层解映射为了cwidx = 1：pdsch.numcodewords qm = length（dlschllrs {cwidx}）/ length（rxsymbols {cwidx}）;每符号%位CSI {CWIDX} = REPMAT（CSI {CWIDX}。'，QM，1）;%按每个符号的每个位展开dlschLLRs{cwIdx} = dlschLLRs{cwIdx} .* csi{cwIdx}(:);百分率(CSI)结束%解码DL-SCH传输信道%将解码CRC错误写入HARQ进程状态结构中decodeDLSCHLocal.TransportBlockLength = trBlkSizes;[decbits，harqProcesses（harqProcIdx）.blkerr] = decodeDLSCHLocal（dlschLLRs，pdsch.Modulation，pdsch.NumLayers，harqProcesses（harqProcIdx）.RV，harqProcIdx-1）;%存储用于计算吞吐量的值(仅针对活动PDSCH实例)如果（任何（trblksize〜= 0））bittput = [bittput trblksizes。*（1-harqprocesses（harqprociddx）.blkerr）];txedtrblksizes = [txedtrblksizes trblksizes];结束%更新HARQ进程计数器harqproccntr = harqproccntr + 1;显示当前HARQ进程管理的每个码字的传输块CRC错误信息ICR = Trblksize./pdschindicsinfo.g;%瞬时码率csn = mod（nslot，cross.slotsperframe）;帧中的％插槽号fprintf中（'\ n（％3.2f %%）harq proc％d：'，100 *（nslot + 1）/ nslots，harqprocidx）;estrings = {“通过”那'失败的'};RVI = HarqProcesses（HarqProcidx）.rvidx;为了CW = 1：长度（RVI）CWRVI = RVI（CW）;％上创建在RV序列和译码错误的RV状态给定的位置的报告如果Cwrvi == 1 ts = sprintf(最初的传播(NSlot = % d,房车= % d, CR = % f)”、csn rvSeq (cwrvi)、icr (cw));别的ts = sprintf（'重传＃％d（nslot =％d，rv =％d，cr =％f）'，cwrvi-1，csn，rvseq（cwrvi），ICR（CW））;结束fprintf中（'cw％d：％s％s。'，cw-1，ts，estrings {1 + harqprocesses（harqprocidx）.blkerr（cw）}）;结束结束％计算最大和模拟吞吐量maxthrougpul（snridx）= sum（txedtrblksizes）;％最大吞吐量simThroughput (snrIdx) = (bitTput, 2)总和;％模拟吞吐量%在命令窗口中动态显示结果流([['\ n \ n \ nthrougpul（mbps）为'，num2str（simLocal.NFrames）'帧（S）']，......'=％.4f \ n'), 1 e-6 * simThroughput (snrIdx) / (simLocal.NFrames * 10 e - 3));流(['吞吐量（%%）关于'，num2str（simLocal.NFrames）'帧(s) = %.4f\n']，......Simkrouguppul（Snridx）* 100 / maxthrougpul（snridx））;结束

图;情节(simParameters.SNRIn simThroughput * 100. / maxThroughput,'o-。')包含('snr（db）');ylabel（'吞吐量（％）');网格在;标题（Sprintf（'%s (%dx%d) / NRB=%d / SCS=%dkHz'那......simParameters.DelayProfile、simParameters.NTxAnts simParameters.NRxAnts,......simparameters.carrier.nsizegrid，simparameters.carrier.subcarrierspacing））;％捆绑密钥参数并导致录制组合结构simResults。simParameters = simParameters;simResults。simThroughput = simThroughput;

功能validateNumLayers (simParameters)％验证相对于天线几何形状的层数nlayers = simparameters.pdsch.numlayers;ntxants = simparameters.ntxants;nrxants = simparameters.nrxants;Antennadecrounction = Sprintf（'min(NTxAnts,NRxAnts) = min(%d，%d) = %d'，ntxants，nrxants，分钟（ntxants，nrxants））;如果Nlayers > min(ntxants,nrxants) error('层数(%d)必须满足NLayers <= %s'那......nlayers，天平页面）;结束如果通道的最大可能等级等于，则显示警告%层数如果（nlayers> 2）&&（nlayers == min（ntxants，nrxants））警告（[“通道的最大可能等级，由%s给出，等于NLayers (%d)。”......“这可能导致某些频道条件下的解码失败。......'尝试减少层数或增加频道排名'......(使用更多的发射或接收天线)。, antennaDescription nlayers);％＃确定结束结束功能estChannelGrid = getInitialChannelEstimate(载体,nTxAnts propchannel)％获得第一传输之前信道估计。这可以用于％获得第一槽的预编码矩阵。ofdmInfo = nrOFDMInfo（载体）;chInfo = INFO（propchannel）;maxChDelay =小区（MAX（chInfo.PathDelays * propchannel.SampleRate））+ chInfo.ChannelFilterDelay;%临时波形(仅对大小需要)tmpWaveform = 0 ((ofdmInfo.SampleRate / 1000 / carrier.SlotsPerSubframe) + maxChDelay nTxAnts);％过滤通过通道[~, pathGains sampleTimes] = propchannel (tmpWaveform);%完美的时间同步pathFilters = getPathFilters (propchannel);抵消= nrPerfectTimingEstimate (pathGains pathFilters);％完美的频道估计estChannelGrid = nrPerfectChannelEstimate(载体,pathGains pathFilters,抵消,sampleTimes);结束功能WTX = getPrecodingMatrix（载体，PDSCH，hestGrid）％计算预编码所有PRGS矩阵在载体重叠的PDSCH分配的百分比通过载波网格解决的％最大CRBmaxCRB =载体。NStartGrid +载波。NSizeGrid - 1;％prg大小如果（isfield（pdsch，'prgbundlesize'）&&〜的isEmpty（pdsch.PRGBundleSize））Pd_BWP = pdsch.PRGBundleSize;别的pd_bwp = maxcrb + 1;结束载波网格中每个RB的% PRG编号(基于1)NPRG = cell ((maxCRB + 1) / Pd_BWP);prgset = repmat ((1: NPRG)、Pd_BWP 1);prgset = prgset(母舰。NStartGrid + (1: carrier.NSizeGrid)。');(~, ~, R, P) = (hestGrid)大小;wtx = 0 ([pdsch。NumLayers P NPRG]);为了我= 1：nprg当前PRG和PDSCH中的％子载波指数%分配thisprg = find（prgset == i） -  1;thisprg =交叉（thiceprg，pdsch.prbset（:) + cross.nstartgrid，'行');prgSc =（1:12） '+ 12 * thisPRG';prgSc = prgSc（:);如果(~ isempty (prgSc))PRG的％平均信道估计eStallocgrid = hestgrid（prgsc，：：，:);hest =允许（平均值（重塑（exhape（extallocgrid，[]，r，p）），[2 3 1]）;％SVD分解[~ ~ V] =圣言(命令);wtx(:,:我)= V (:, 1: pdsch.NumLayers)。”;结束结束wtx = wtx / sqrt（pdsch.numlayers）;Numlayers的％正常化结束功能estChannelGrid = precodeChannelEstimate(载体、estChannelGrid W)％将预编码矩阵W应用于频道估计的最后一个维度(K, L R P) = (estChannelGrid)大小;estChannelGrid =重塑(estChannelGrid，[K*L R P]);estChannelGrid = hPRGPrecode ([K L R P], carrier.NStartGrid estChannelGrid,重塑(1:元素个数(estChannelGrid), [K * L R P]), W);estChannelGrid =重塑(estChannelGrid K, L, R, []);结束功能[MappedPro，MappedSymbols] = MapNumerology（子载波，符号，NRB，NRBT，FS，FT）％将SSBurst数字映射到PDSCH数字。输出是：% - mappedPRB:载波资源网格的基于0的PRB索引(按列排列)％ -  MapedSymbols：载波资源网格中的插槽中的基于0的OFDM符号索引（连续排列）%输入参数为:％ - 子载波：SSB资源网格的基于1的行下标（排列在列中）％ - 符号：用于SSB资源网格的基于1的列下标（以n×4矩阵排列，为行中的每个传输突发的4个符号，n发送突发）％SSB资源网格使用SSBINFO.nRB，普通CP，跨越5子帧来划分% - nrbs:源(SSB) NRB% - nrbt:目标(载波)NRB％ -  FS：源（SSB）SCS％ - 英尺：靶（载体）SCSMappedPrb =唯一（FIX（（子载波 - （NRB * 6） -  1）* FS /（FT * 12）+ NRBT / 2），'稳定的');符号=符号。';符号=符号（:)。'-  1;如果（FT %如果ft/fs < 1，减少mappedsymbols =唯一（修复（符号* ft / fs），'稳定的');别的% Else, repeat by ft/fsMappedSymbols = REPAPE（（0：（0：（ft / fs-1））'+符号（:)'* ft / fs，1，[]）;结束结束功能Plotlayerevm（NSLOTS，NSLOT，PDSCH，SIZ，PDSCHINDICE，PDSCHSYMBOL，PDSCHEQ）％plot EVM信息执着的slotevm;执着的rbEVM执着的evmPerSlot;如果（nslot == 0）slotEVM = comm.EVM;rbEVM = comm.EVM;evmPerSlot = NaN的（NSLOTS，pdsch.NumLayers）;数字;结束Evmperslot（nslot + 1，:) = slotevm（pdschsymbols，pdscheq）;子图（2,1,1）;图（0：（nslots-1），EVMPERSLOT）;Xlabel（'插槽号');ylabel（'维生素(%)');传奇(“层”+ (1: pdsch.NumLayers),“位置”那'eastoutside');标题('每个插槽每层EVM');次要情节(2,1,2);(k ~ p) = ind2sub (siz pdschIndices);Rbsubs = floor((k-1) / 12);NRB = siz(1) / 12;evmPerRB =南(NRB pdsch.NumLayers);为了NU = 1：pdsch.NumLayers为了RB =唯一（RBSBS）。这个=（rbsubs == rb＆p == nu）;EVMPERRB（RB + 1，NU）= RBEVM（PDSCHSYMBOLS（此），PDSCHEQ（此））;结束结束图（0：（NRB-1），EVMPerrb）;Xlabel（'资源块');ylabel（'维生素(%)');传奇(“层”+ (1: pdsch.NumLayers),“位置”那'eastoutside');标题(['EVM每层每个资源块，slot #'num2str (nslot)]);drawnow;结束