simParameters = [];％清除SimParameters变量simParameters。NFrames = 2;% 10ms帧数simParameters。SNRIn = [-5 0 5];%信噪比范围(dB)

displaysimulationInformation = true;

PerfectChannElestimator = True;

%带宽、命理(SCS和CP类型)等通用参数simparameters.nrb = 52;资源块数量的％带宽（10MHz BW为15kHz SCS的52RB）simParameters.SubcarrierSpacing=15；％15,30,60,120,240（kHz）simParameters。CyclicPrefix =“正常”；% '普通'或'扩展'simParameters。NCellID = 0;%细胞身份simparameters.ntxants = 1;%发送天线数simParameters。NRxAnts = 2;%接收天线数％UL-SCH / PUSCH参数simParameters.PUSCH.TargetCodeRate = 193 / 1024;%用于计算传输块大小的码率simParameters.PUSCH.PRBSet = (0: simParameters.NRB-1);％PUSCH PRB分配simParameters.PUSCH.SymbolSet = 0:13;每个槽中的PUSCH符号分配simparameters.pusch.nohprb = 0;%每个PRB额外的RE开销simParameters.PUSCH.EnableHARQ = true;％启用/禁用HARQ，如果禁用，则具有RV = 0的单个传输，即无重传simParameters.PUSCH.Modulation ='QPSK'；% 'pi/2-BPSK'， 'QPSK'， '16QAM'， '64QAM'， '256QAM'simParameters.PUSCH.NLayers = 1;% push层数simParameters.PUSCH.RNTI = 1;%无线网络临时标识符simParameters.PUSCH.TransformPrecoding = false;％启用/禁用变换预编码simParameters.PUSCH.TxScheme =“nonCodebook”；%传输方案('nonCodebook'，'codebook')simParameters.PUSCH.NAntennaPorts = 1;%用于码本预编码的天线端口数simParameters.PUSCH.TPMI = 0;基于码本的预编码预测矩阵指示% push DM-RS配置simParameters.PUSCH.PUSCHMappingType =“一个”；% PUSCH映射类型('A'(slot-wise)，'B'(non - slot-wise))simparameters.pusch.dmrstypeaposition = 2;%仅映射类型A。第一个DM-RS符号位置(2,3)simparameters.pusch.dmrslength = 1;%前置DM-RS符号数(1个(单符号)，2个(双符号))simparameters.pusch.dmrsadditionalposition = 1;%附加DM-RS符号位置(最大范围0…3)simparameters.pusch.dmrsconfigurytype = 1;％DM-RS配置类型（1,2）simParameters.PUSCH.NumCDMGroupsWithoutData = 2;%无数据的CDM组simParameters.PUSCH.NIDNSCID = 0;%置乱标识(0…65535)simParameters.PUSCH.NSCID = 0;%置乱初始化(0,1)simParameters.PUSCH.NRSID = 0;低papr序列的打乱ID(0…1007)simparameters.pusch.grouphopping =.'禁用'；％跳跃类型（'启用'，'禁用'）％定义传播信道类型simParameters。ChannelType ='tdl'；％'cdl'或'tdl'

问题= simParameters;pusch = simParameters.PUSCH;

snrIn = simParameters.SNRIn;nTxAnts = simParameters.NTxAnts;nRxAnts = simParameters.NRxAnts;channelType = simParameters.ChannelType;

如果strcmpi (channelTypeCDL的) channel = nrCDLChannel;通道。DelayProfile =“CDL-A”；[txsize, rxsize] = hArrayGeometry (nTxAnts nRxAnts,'上行链路'）;channel.TransmitAntennaArray.Size = txsize;channel.ReceiveAntennaArray.Size = rxsize;其他的频道= nrTDLChannel;通道。DelayProfile =“TDL-A”；channel.numtransmitantennas = ntxants;channel.numreceiveantennas = nrxants;结束channel.delayspread = 30e-9;%在几秒钟内通道。MaximumDopplerShift = 10;%在赫兹

waveformInfo = hOFDMInfo(问题);通道。SampleRate = waveformInfo.SamplingRate;

chInfo =信息(渠道);maxChDelay =装天花板(max (chInfo.PathDelays * channel.SampleRate));maxChDelay = maxChDelay + chInfo.ChannelFilterDelay;

％数组存储所有SNR点的最大吞吐量maxthrougpul = zeros（长度（snrin），1）;%阵列存储所有信噪比点的模拟吞吐量simThroughput = 0(长度(snrIn), 1);%设置冗余版本(RV)序列，HARQ进程数和%表示使用HARQ进程的顺序如果pusch。EnableHARQ%在RAN WG4中PUSCH解调要求满足#88bis% (r4 - 1814062)rvSeq = [0 2 3 1];其他的% HARQ禁用-单次传输，RV=0，没有重传rvseq = 0;结束%创建UL-SCH编码器系统对象encodeULSCH = nrULSCH;encodeULSCH。MultipleHARQProcesses = true;encodeULSCH。TargetCodeRate = pusch.TargetCodeRate;创建UL-SCH解码器系统对象% LDPC译码采用分层信念传播，次数减半与信念传播解码的默认值相比，%迭代decodeULSCH = nrULSCHDecoder;decodeULSCH。MultipleHARQProcesses = true;decodeULSCH。TargetCodeRate = pusch.TargetCodeRate;decodeULSCH。LDPCDecodingAlgorithm =“分层信念传播”；decodeULSCH。MaximumLDPCIterationCount = 6;％临时变量'ue_init'，'pusch_init'和'decodeulsch_init'%用于创建临时变量ue、push和% 'decodeULSCH'在信噪比循环中创建独立实例并行模拟%ue_init =问题;pusch_init = pusch;decodeULSCH_init =克隆(decodeULSCH);为snridx = 1：numel（snrin）%注释掉用于并行计算% parfor snrIdx = 1:numel(snrIn) % uncomment用于并行计算为了减少总的模拟时间，你可以执行这个循环% parallel使用并行计算工具箱。注释掉for%语句并取消“parfor”语句的注释。如果并行计算% Toolbox未安装，'parfor'默认为普通的'for'语句。因为parfor-loop迭代是并行执行的%不确定顺序，模拟信息显示为每个信噪比百分点可以交织在一起。要关闭模拟信息显示，％将上面的“DisplaySimulationInmation”变量设置为FALSE重置随机数生成器和信道，使每个信噪比点%将经历相同的噪声和通道实现rng (“默认”）;重置(渠道);%初始化这个SNR点的变量，需要初始化使用并行计算工具箱时变量的百分比问题= ue_init;pusch = pusch_init;pathFilters = [];%指定HARQ进程循环的顺序NHARQProcesses = 16;harqSequence = 1: NHARQProcesses;％初始化所有HARQ进程的状态并重置UL-SCH%译码器harqProcesses = hNewHARQProcesses (NHARQProcesses rvSeq 1);harqProcCntr = 0;％HARQ过程计数器decodeULSCH =克隆(decodeULSCH_init);SNRdB = snrIn (snrIdx);流('\n模拟基于%s的传输方案(%dx%d)和SCS=%dkHz， %s通道在%gdB SNR for %d 10ms frame(s)\n', nTxAnts pusch.TxScheme nRxAnts、ue.SubcarrierSpacing channelType, SNRdB, ue.NFrames);%仿真周期内OFDM符号总数waveformInfo = hOFDMInfo(问题);NSymbols =问题。NFrames * 10 * waveformminfo . symbolspersubframe;与每次PUSCH传输开始相关的% OFDM符号号UE.NSYMBOL = 0;PUSCH传输实例数量的％运行计数器%模拟将使用此计数器作为每个槽位号% PUSCHpusch。NSlot = 0;%定时偏移，更新在每个插槽完美的同步和%，用于实际的同步抵消= 0;虽然问题。NSymbol < NSymbols％计算此插槽的传输块大小[PUSCHINDICES，DMRSINDICES，DMRSSYMBOLS，PUSCHINDISSINFO] = HPUSCHRESOUSE（UE，PUSCH）;TBS = HPUSCHTBS（PUSCH，PUSCHINDICESFO.NREPERPRB  -  PUSCH.NOHPRB）;从HARQ索引获取当前PUSCH的HARQ进程索引％ 桌子harqProcIdx = harqSequence (mod (harqProcCntr长度(harqSequence)) + 1);%更新当前HARQ进程信息(这将更新RV%取决于CRC通过或失败在以前的传输%此HARQ进程)HarqProcesses（HarqProcidx）= hupdateharqprocess（HarqProcesses（HarqProciddx），1）;% HARQ:检查以前传输的CRC，即是a%需要重传?NDI = false;如果harqProcesses .blkerr (harqProcIdx)％错误如果(harqProcesses (harqProcIdx) .RVIdx = = 1)rvSeq的% endresetSoftBuffer (decodeULSCH harqProcIdx-1);NDI = true;结束其他的%没有错误NDI = true;结束如果NDI trBlk = randi([0 1]，TBS,1); / /输出setTransportBlock (encodeULSCH trBlk harqProcIdx-1);结束% UL-SCH编码codedTrBlock = encodeULSCH (pusch.Modulation pusch.NLayers、puschIndicesInfo.G harqProcesses (harqProcIdx) .RV harqProcIdx-1);% PUSCH调制，包括基于码本的MIMO预编码if% TxScheme = 'codebook'MRB =元素个数(pusch.PRBSet);puschSymbols = nrPUSCH (codedTrBlock pusch.Modulation、pusch.NLayers ue.NCellID, pusch.RNTI, pusch.TransformPrecoding, MRB, pusch.TxScheme, pusch.NAntennaPorts, pusch.TPMI);创建与PUSCH传输周期相关的资源网格puschGrid = 0 (waveformInfo.NSubcarriers waveformInfo.SymbolsPerSlot nTxAnts);％实现特定于PUSCH MIMO预编码和映射。这个％MIMO预编码步骤是基于任何码本的% MIMO预编码做PUSCH调制以上如果（strcmpi（pusch.txscheme，“码”）)%码本的MIMO预编码，F预编码之间的PUSCH%发射天线端口和发射天线F =眼睛(pusch.NAntennaPorts nTxAnts);其他的基于非码本的MIMO预编码，PUSCH之间的F比%层和发射天线F =眼睛(pusch.NLayers nTxAnts);结束[~, puschAntIndices] = nrExtractResources (puschIndices puschGrid);puschGrid(puschAntIndices) = puschSymbols * F;%实现特定的PUSCH DM-RS MIMO预编码和映射。%在hPUSCHResources上面的DM-RS创建包括代码本%基于MIMO预编码(如适用)为p = 1:size(dmrsSymbols,2) [~，dmrsAntIndices] = nrExtractResources(dmrsIndices(:，p)，puschGrid);puschGrid(dmrsAntIndices) = puschGrid(dmrsAntIndices) + dmrsSymbols(:，p) * F(p，:);结束% OFDM调制puschGrid txWaveform = hOFDMModulate(问题);%通过通道模型传递数据。的末尾加0%发送波形以刷新通道内容。这些零%考虑信道中引入的任何延迟。这是一种混合%多径时延与实现时延的关系。该值可能会%的变化取决于采样率，延迟轮廓和延迟%的传播txWaveform = [txWaveform;0 (maxChDelay、尺寸(txWaveform 2)));% #好< AGROW >[rxWaveform, pathGains sampleTimes] =通道(txWaveform);%将AWGN添加到接收的时域波形中根据OFDM调制中使用的IFFT尺寸对噪声功率进行归一化，%时，OFDM调制器将这种归一化应用到%传播波形。也要按照接收的数量进行规范%天线，因为信道模型的默认行为是%将此归一化应用于接收的波形信噪比= 10 ^ (SNRdB / 20);N0 = 1 /(√(2.0 * nRxAnts *双(waveformInfo.Nfft)) *信噪比);噪音= N0 *复杂(randn(大小(rxWaveform)), randn(大小(rxWaveform)));rx波形= rx波形+噪声;如果（PerfectChanneLestimator）％完美同步。使用提供的信息%通道以找到最强的多路径组件pathFilters = getPathFilters(渠道);(抵消,mag) = nrPerfectTimingEstimate (pathGains pathFilters);其他的%实际同步。关联接收波形%与PUSCH DM-RS给出时间偏移估计't'和%相关幅度'mag'。这个函数% hskipweaktiminggoffset用于更新接收机定时％ 抵消。如果“mag”中的相关峰值弱，则当前%计时估计't'被忽略，之前的估计使用％'offset'(t, mag) = nrTimingEstimate (rxWaveform、ue.NRB ue.SubcarrierSpacing, pusch.NSlot, dmrsIndices, dmrsSymbols,“CyclicPrefix”, ue.CyclicPrefix);% #好< UNRCH >抵消= hSkipWeakTimingOffset(抵消t杂志);结束rxWaveform = rxWaveform(1 +抵消:最终,);%对接收到的数据执行OFDM解调以重新创建％资源网格，包括实用的事件中的填充%同步导致一个不完整的槽被解调rxgrid = hofdmdemodulate（ue，rxwaveform）;[k，l，r] =尺寸（rxgrid）;如果（l 结束如果（PerfectChanneLestimator）%完美信道估计，使用路径增益的值%由通道提供estChannelGrid = nrPerfectChannelEstimate (pathGains pathFilters、ue.NRB ue.SubcarrierSpacing, pusch.NSlot,抵消,sampleTimes, ue.CyclicPrefix);％获得完美的噪音估计（来自噪声实现）noiseGrid = hOFDMDemodulate(问题、噪声(1 +抵消:最终,));噪音= var (noiseGrid (:));%对estChannelGrid应用MIMO预编码给出估计每传输层%K =大小(estChannelGrid, 1);estChannelGrid =重塑(estChannelGrid K * waveformInfo.SymbolsPerSlot * nRxAnts nTxAnts);estChannelGrid = estChannelGrid * F.';如果（strcmpi（pusch.txscheme，“码”)) W = nrPUSCHCodebook(pusch.NLayers,pusch. ntennapports,pusch.TPMI,pusch.TransformPrecoding);estChannelGrid = estChannelGrid * W.';结束estChannelGrid =重塑(estChannelGrid K waveformInfo.SymbolsPerSlot nRxAnts, []);其他的%实际信道估计之间的接收网格和%每一传输层，使用PUSCH DM-RS为每一层[〜，dmrslayerindices，dmrslayersymbols] = hpuschresources（Ue，setfield（pusch，“TxScheme”，“nonCodebook”))；% #好< UNRCH >[estChannelGrid,噪音]= nrChannelEstimate (rxGrid、dmrsLayerIndices dmrsLayerSymbols,“CyclicPrefix”,问题。CyclicPrefix,'cdmlength', puschIndicesInfo.CDMLengths);结束％从收到的网格获取PUSCH资源元素[puschRx, puschHest] = nrExtractResources (puschIndices、rxGrid estChannelGrid);%均衡[puschEq，csi]=nrEqualizeMMSE（puschRx，puschHest，noiseEst）；％解码PUSCH物理通道[ulschLLRs, rxSymbols] = nrPUSCHDecode (puschEq、pusch.Modulation ue.NCellID, pusch.RNTI,噪音,pusch.TransformPrecoding, MRB);%应用均衡器产生的信道状态信息(CSI)，%如果启用，包括转换预编码的效果如果(pusch.TransformPrecoding) MSC = MRB * 12;csi = nrTransformDeprecode(csi,MRB) / sqrt(MSC);csi = repmat (csi ((1: MSC:结束)。”),1,MSC)。“;csi =重塑(csi,大小(rxSymbols));结束csi = nrLayerDemap (csi);Qm = length(ulschLLRs) / length(rxSymbols);csi =重塑(repmat (csi{1}”。Qm, 1), [], 1);ulschLLRs = ulschLLRs .* csi;%译码UL-SCH传输信道decodeULSCH。TransportBlockLength = TBS;[decbits harqProcesses (harqProcIdx)。blkerr] = decodeULSCH (ulschLLRs、pusch.Modulation pusch.NLayers, harqProcesses (harqProcIdx) .RV harqProcIdx-1);%存储计算吞吐量的值simThroughput(snrIdx) = simThroughput(snrIdx) + (~harqProcesses(harqProcIdx))。blkerr * TBS);maxThroughput(snrIdx) = maxThroughput(snrIdx) + TBS;显示传输块错误信息如果(displaySimulationInformation)流('\ n（％3.2f %%）harq proc％d：'，100 *（UE.NSYMBOL +大小（PUSCHGRID，2））/ nsymbols，HarqProcidx）;estrings = {“通过”，“失败”};rvi = harqProcesses (harqProcIdx) .RVIdx;如果Rvi == 1 ts = sprintf('初始传输（RV=%d）'，rvseq（RVI））;其他的ts = sprintf (“重传# % d (RV = % d) '、rvi-1 rvSeq (rvi));结束流('％s％s。'ts, estrings {1 + harqProcesses (harqProcIdx) .blkerr});结束%更新下一个PUSCH传输的开始符号号UE.NSYMBOL = UE.NSYMBOL +大小（PUSCHGRID，2）;% PUSCH传输总数的更新计数pusch。NSlot = pusch。NSlot + 1;%更新HARQ进程计数器harqProcCntr = harqProcCntr + 1;结束%在命令窗口中动态显示结果如果(displaySimulationInformation)流(' \ n '）;结束流([[“\n的总吞吐量（Mbps）”num2str（UE.nframes）“帧(s)”]，“= % .4f \ n”]，1E-6 * SimThrougupul（SNRIDX）/（UE.NFRAMES * 10E-3））;fprintf（[的吞吐量(% %)num2str（UE.nframes）'帧(s) = %.4f\n'), simThroughput (snrIdx) * 100 / maxThroughput (snrIdx));结束

图;情节(snrIn simThroughput * 100. / maxThroughput,'o-。'）Xlabel（“信噪比(dB)”）;ylabel（的吞吐量(%)）;网格在；如果(pusch_init.TransformPrecoding) ofdmType =“DFT-s-OFDM”；其他的ofdmType =“CP-OFDM”；结束标题(sprintf ('%s / NRB=%d / sc =%dkHz / %s %d/1024 / %dx%d'，．..ofdmType、ue_init.NRB ue_init.SubcarrierSpacing,．..pusch_init。调制,．..轮(pusch_init.TargetCodeRate * 1024)、nTxAnts nRxAnts));simResults。simParameters = simParameters;simResults。simThroughput = simThroughput;simResults。maxThroughput = maxThroughput;