numTrBlks = 4;模拟传输块的数目SNRdB = 32:4:0;％SNR范围以dB为单位ireps = [0 5 9];代表的范围％模拟

NPDSCHDataType ='NotBCCH';％的允许值是“SIB1NB”，“BCCHNotSIB1NB”或“NotBCCH”

ISF = 0;% DCI资源分配字段(DCI格式N1或N2)SchedulingInfoSIB1 = 0;MasterInformationBlock-NB (MIB-NB)中的调度信息字段%imc = 4;DCI中%调制和编码方案域(DCI格式N1或N2)

enb.NFrame = 0;％仿真起始帧编号enb.NSubframe = 0;%模拟启动子帧数enb.NNCellID = 0;％NB-的IoT的物理小区IDenb.NBRefP = 2;NRS天线端口的数目，应该是1或2enb.OperationMode=“Inband-DifferentPCI”;％的允许值是“带内-SamePCI”，“带内-DifferentPCI”，“防护频带”或“独立”如果strcmpi（enb.OperationMode，“Inband-SamePCI”）enb.CellRefP = enb.NBRefP;允许的值是NBRefP或4enb。NCellID = enb.NNCellID;elseifstrcmpi（enb.OperationMode，“Inband-DifferentPCI”)enb。CellRefP = 4;细胞的％数Rs天线端口（必须等于NBRefP或4）enb。NCellID = 1;结束如果(strcmpi (NPDSCHDataType'BCCHNotSIB1NB')| | strcmpi (NPDSCHDataType'NotBCCH'））&&...(strcmpi (enb.OperationMode“Inband-SamePCI”）||strcmpi（enb.OperationMode，“Inband-DifferentPCI”））enb.ControlRegionSize = 3;％的允许值为0 ... 13结束

频道=结构;％初始化信道的配置结构channel.Seed = 6;％频道种子通道。NRxAnts = 1;％1接收天线channel.DelayProfile ='EPA';％延迟分布通道。DopplerFreq = 5;%多普勒频率(Hz)通道。MIMOCorrelation =“低”;% Multi-antenna相关性通道。NTerms = 16;在衰落模型用于％振荡器channel.ModelType =“GMEDS”;％瑞利衰落模型类型通道。InitPhase ='随机';％随机初始阶段channel.NormalizePathGains ='上';％正常化延迟分布电源channel.NormalizeTxAnts ='上';用于发射天线％规格化

信道估计器行为perfectChannelEstimator = true;

％配置信道估计器cec。PilotAverage ='用户自定义';导频码元的平均％类型cec。TimeWindow = 1;在res％的时间窗口大小cec.FreqWindow = 25;频率窗口大小(以REs为单位)cec.InterpType ='立方体';％2D插值类型cec。InterpWindow =“中心”;插值窗口类型cec.InterpWinSize = 3;％插值窗口大小cec。[Reference ='NRS';％信道估计的参考信号

对于repIdx = 1:元素个数(ireps)

npdschInfo = hNPDSCHInfo;npdschInfo。NPDSCHDataType = NPDSCHDataType;npdschInfo。安全部队=安全部队;如果strcmpi (NPDSCHDataType'SIB1NB'）% NPDSCH携带SIB1-NBnpdschInfo。SchedulingInfoSIB1 = SchedulingInfoSIB1;其他% NPDSCH不携带SIB1-NBnpdschInfo.IRep = ireps（repIdx）;DCI中%重复数字段(DCI格式N1或N2)npdschInfo。IMCS= IMCS;DCI中%调制和编码方案域(DCI格式N1或N2)%验证IRep和IMCS的输入如果npdschInfo isempty (npdschInfo.TBS)。TBSTable错误(['无效[ITBS，ISF]（其中ITBS = IMCS ='num2str (imc)...'ISF ='num2str(安全部队)“）对，则返回空TBS，检查有效对在上表中或3GPP TS 36.213表16.4.1.5.1-1”]）;结束结束

npdsch.NSF = npdschInfo.NSF;npdsch.NRep = npdschInfo.NRep;npdsch.NPDSCHDataType = NPDSCHDataType;npdsch.RNTI = 1;

[〜，信息] = lteNPDSCHIndices（ENB，npdsch）;rmoutlen = info.G;码率匹配后的比特长度，即码字长度trblklen = npdschInfo.TBS;传输块大小%R = (trblklen + 24) / rmoutlen;% DL-SCH信道编码率，24为CRC比特数如果R> = 1个错误（['DL-SCH的编码率（'num2str（R）“）大于或等于1的配置的参数大”。]）;结束

当前配置的NPDSCH重复模式为%显示下面displayPattern = FALSE;％显示NPDSCH重复图形如果displayPattern == true npdschInfo.displaySubframePattern;结束

在仿真起始点的绝对子帧数NSubframe = enb.NFrame * 10 + enb.NSubframe;初始化BLER和吞吐量结果maxThroughput =零（长度（SNRdB），1）;simThroughput =零（长度（SNRdB），1）;BLER =零（1，numel（SNRdB））;％的临时变量“enb_init”和“channel_init”用于创建在SNR循环中创建临时变量'enb'和'channel'%独立的模拟循环的'parfor'循环enb_init = enb;channel_init =通道;对于snrIdx = 1:元素个数(SNRdB)％PARFOR snrIdx = 1：numel（SNRdB）允许使用并行计算来提高注释输出速度％的“对”语句的上方和下方取消注释“PARFOR”语句。％这需要并行计算工具箱。如果不安装％“PARFOR”将默认为正常的“For”语句。％，取决于给循环变量设置随机数生成器种子％，以确保独立的随机流RNG（snrIdx，'combRecursive');fprintf中（'\ nSimulating％d原子％GDB SNR的传输块\ N'、numTrBlks SNRdB (snrIdx));enb = enb_init;％初始化的eNodeB配置信道= channel_init;％初始化衰落信道配置txcw = [];初始化传输的码字numBlkErrors = 0;有错误的传输块的数目地产= [];％初始化NPDSCH编码器状态dstate = [];％初始化NPDSCH解码器状态lastOffset = 0;％初始化整体帧时序偏移偏移量= 0;初始化帧定时偏移量enb.NBRefP subframeGrid = 0 (12、14);％初始化时隙格subframeIdx = NSubframe;numRxTrBlks = 0;而（numRxTrBlks ％设定电流和子帧号enb.NSubframe = MOD（subframeIdx，10）;enb.NFrame =地板（（subframeIdx）/ 10）;％生成NPSS符号和标记npssSymbols = lteNPSS（ENB）;npssIndices = lteNPSSIndices（ENB）;％映像符号副车架网格subframeGrid（npssIndices）= npssSymbols;%生成NSSS符号和索引nsssSymbols = lteNSSS (enb);nsssIndices = lteNSSSIndices (enb);％映像符号副车架网格subframeGrid（nsssIndices）= nsssSymbols;%确定是否传播了NPSS或NSSS，如果是，在此子帧中不传输NPDSCHisDataSubframe =的isEmpty（npssSymbols）&&的isEmpty（nsssSymbols）;％创建一个新的传输块，并且当它编码％发送的码字是空的。接收机设定码%为空，表示包中的所有子帧都已完成%已接收(在第一次传输之前也是空的)如果的isEmpty（txcw）txTrBlk =兰迪（[0 1]，trblklen，1）;txcw = lteNDLSCH（rmoutlen，txTrBlk）;结束如果(isDataSubframe)为一个子帧生成NPDSCH符号和索引[txNpdschSymbols，房地产] = lteNPDSCH（ENB，npdsch，txcw，房地产）;npdschIndices = lteNPDSCHIndices（ENB，npdsch）;％映像符号副车架网格subframeGrid (npdschIndices) = txNpdschSymbols;％生成NRS符号和标记nrsSymbols = lteNRS (enb);nrsIndices = lteNRSIndices (enb);％映像符号副车架网格subframeGrid (nrsIndices) = nrsSymbols;结束％进行OFDM调制，以生成时域波形[txWaveform，ofdmInfo] = nbOFDMModulate（ENB，su​​bframeGrid）;增加25个样品填充。这是为了覆盖延迟的范围从信道建模的组合（％预期％执行延迟和信道延迟扩展）txWaveform = [txWaveform;零（25，enb.NBRefP）];％＃确定为每个子帧初始化通道时间通道。InitTime = subframeIdx / 1000;通过通道模型传递数据channel.SamplingRate = ofdmInfo.SamplingRate;[rxWaveform，fadingInfo] = lteFadingChannel（信道，txWaveform）;％计算噪声增益包括用于下行链路功率补偿％分配信噪比= 10 ^ (SNRdB (snrIdx) / 20);％规格化噪声功率，以考虑采样率，其中的％是在OFDM调制中使用的IFFT大小的函数，并％的天线的数量N0 = 1 /（SQRT（2.0 * * enb.NBRefP双（ofdmInfo.Nfft））* SNR）;％创建加性高斯白噪声噪音= N0 *复杂(randn(大小(rxWaveform)),...randn(大小(rxWaveform)));在接收到的时域波形中加入AWGNrx波形= rx波形+噪声;％------------------------------------------------------------------％接收机％------------------------------------------------------------------％执行定时同步，提取相应的%子帧接收的波形，并执行OFDM％解调如果(perfectChannelEstimator) offset = hPerfectTimingEstimate(fadingInfo);其他％。在该示例中，子帧偏移计算依赖％上NPSS存在于子帧5，所以我们需要垫％的子帧之前，以使帧偏移通过返回% lteNBDLFrameOffset是子帧5的偏移量sfTsamples = ofdmInfo.SamplingRate * 1 e - 3;如果（enb.NSubframe == 5）填充=零（[sfTsamples * 5，尺寸（rxWaveform，2）]）;偏移量= lteNBDLFrameOffset（ENB，[填充; rxWaveform]）;如果(offset > 25) || (offset < 0) offset = lastOffset;结束lastOffset =抵消;结束结束％同步接收波形rx波形= rx波形(1+偏移量:end，:);％对接收的数据执行OFDM解调以重新创建％资源网格rxSubframe = nbOFDMDemodulate (enb rxWaveform);%信道估计如果(perfectChannelEstimator)完美信道估计estChannelGrid = nbDLPerfectChannelEstimate(enb, channel, offset);噪声组= nbOFDMDemodulate(enb, noise(1+偏移:end，:));噪音= var (noiseGrid (:));其他[estChannelGrid, noise eest] = lteDLChannelEstimate(...基站，CEC，rxSubframe）;结束如果(isDataSubframe)%得到NPDSCH指数npdschIndices = lteNPDSCHIndices(enb, npdsch);从接收的子帧中获取PDSCH资源元素。规模的％的由PDSCH功率因数的Rho接收子帧。PDSCH是％缩放通过该量，而用于小区的参考符号％的信道估计（在PDSCH解码阶段中使用）则不是。[rxNpdschSymbols, npdschHest] = lteExtractResources(npdschIndices，...rxSubframe，estChannelGrid）;%解码NPDSCH[rxcw dstate,符号]= lteNPDSCHDecode (...enb, npdsch, rxNpdschSymbols, npdschHest, noise eest,dstate);％解码的传输块时，在一个捆绑所有子帧%已收到如果dstate.EndOfTx [trblkout，blkerr] = lteNDLSCHDecode（trblklen，rxcw）;numBlkErrors = numBlkErrors + blkerr;numRxTrBlks = numRxTrBlks + 1;重新初始化以启用新传输块的传输txcw = [];结束结束subframeIdx = subframeIdx + 1;结束%计算块错误率提单(snrIdx) = numBlkErrors / numTrBlks;fprintf中（'NPDSCH BLER =％.4f \ N'，BLER（snrIdx））;%计算最大吞吐量和模拟吞吐量maxThroughput (snrIdx) = trblklen * numTrBlks;最大可能吞吐量simThroughput (snrIdx) = trblklen * (numTrBlks-numBlkErrors);%模拟吞吐量fprintf中（'NPDSCH吞吐量（%%）=％.4f %% \ N'，simThroughput（snrIdx）* 100 / maxThroughput（snrIdx））;结束

如果repIdx == 1 FH =图。格上;持有上;包含(“信噪比(dB)”);ylabel (“提单”);legendstr = {[“NRep = 'num2str (npdsch.NRep)]};其他勒让德str =[勒让德str]“NRep = 'num2str (npdsch.NRep)];％＃确定结束图（FH）;图（SNRdB，BLER，'-o');

结束％设定人物称号如果strcmpi (NPDSCHDataType'SIB1NB'）npdsch.NSF = 8;结束标题([' 'CHAR（npdsch.NPDSCHDataType）”:TBS = 'num2str（trblklen）...”;NSF = 'num2str（npdsch.NSF）”;“num2str (enb_init.NBRefP)“NRS端口(s)”]）;图例（legendstr）;

％NB-的IoT DL OFDM调制器功能(波形信息]= nbOFDMModulate (enb、网格)％根据TS 36.104表E.5.1-1a应用默认窗口大小如果（〜isfield（ENB，“窗口”enb))。窗口= 6;结束％用途NB-的IoT SC-FDMA来获得关于OFDM调制的1/2子载波移位enb。NBULSubcarrierSpacing =“15 khz”;[波形，信息] = lteSCFDMAModulate（ENB，网格）;结束％NB-的IoT DL OFDM解调器功能格= nbOFDMDemodulate（ENB，rxWaveform）％用途NB-的IoT SC-FDMA来获得关于OFDM调制的1/2子载波移位enb。NBULSubcarrierSpacing =“15 khz”;网格= lteSCFDMADemodulate (enb, rxWaveform, 0.55);％0.55 CP分数结束完美的信道估计器功能H = nbDLPerfectChannelEstimate (enb、通道timefreqoffset)重新配置NB-IoT UL完美通道估计器来执行DL完美％的信道估计enb。NBULSubcarrierSpacing =“15 khz”;enb。NTxAnts = enb.NBRefP;enb。TotSlots = 2;H = lteULPerfectChannelEstimate(enb, channel,timefreqoffset);结束