numTrBlks = 4;%模拟传输块数SNRdB = - 34:0;%信噪比范围，单位为dBIreps = [0 5 9];%模拟的次数范围

NPDSCHDataType =“NotBCCH”；%允许的值为'SIB1NB'， 'BCCHNotSIB1NB'或'NotBCCH'

Isf = 0;% DCI中的资源分配字段(DCI格式N1或N2)SchedulingInfoSIB1 = 0;MasterInformationBlock-NB中的调度信息字段(MIB-NB)Imcs = 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”；允许的值是“Inband-SamePCI”，“Inband-DifferentPCI”，“Guardband”或“Standalone”如果strcmpi (enb。OperationMode,“Inband-SamePCI”) enb。CellRefP = enb.NBRefP;%允许取值为NBRefP或4enb。NCellID = enb.NNCellID;elseifstrcmpi (enb。OperationMode,“Inband-DifferentPCI”) enb。CellRefP = 4;% Cell 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结束

通道= struct;初始化通道配置结构通道。种子= 6;%沟道种子通道。NRxAnts = 1;% 1接收天线通道。DelayProfile =“环保署”；%延迟配置文件通道。DopplerFreq = 5;%多普勒频率，单位为Hz通道。MIMOCorrelation =“低”；%多天线相关性通道。NTerms = 16;衰落模型中使用的振荡器通道。ModelType =“GMEDS”；瑞利衰落模型类型通道。InitPhase =“随机”；%随机初始相通道。NormalizePathGains =“上”；正常化延迟配置文件功率通道。NormalizeTxAnts =“上”；用于发射天线归一化

%通道估计器行为perfectChannelEstimator = true;

配置信道估计器cec。PilotAverage =“UserDefined”；导频符号平均的类型cec。TimeWindow = 1;%时间窗口大小在REscec。FreqWindow = 25;%频率窗口大小在REscec。InterpType =“立方”；% 2D插值类型cec。InterpWindow =“中心”；插补窗口类型cec。InterpWinSize = 3;插值窗口大小%cec。Reference =“关系”；信道估计器参考信号

为repIdx = 1:数字(ireps)

npdschInfo = hNPDSCHInfo;npdschInfo。NPDSCHDataType = NPDSCHDataType;npdschInfo。Isf = Isf;如果strcmpi (NPDSCHDataType“SIB1NB”）携带SIB1-NB的% NPDSCHnpdschInfo。SchedulingInfoSIB1 = SchedulingInfoSIB1;其他的% NPDSCH不携带SIB1-NBnpdschInfo。IRep= ireps(repIdx);DCI中重复数字段的百分比(DCI格式N1或N2)npdschInfo。imc= IMCS;DCI中调制和编码方案字段(DCI格式N1或N2)验证IRep和IMCS的输入如果npdschInfo isempty (npdschInfo.TBS)。TBSTable错误([无效[ITBS,ISF](其中ITBS=IMCS='num2str (imc).．.”,安全部队= '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;

[~，info] = lteNPDSCHIndices(enb,npdsch);rmoutlen = info.G;速率匹配后的比特长度，即码字长度trblklen = npdschInfo.TBS;%传输块大小R = (trblklen+24)/rmoutlen;% DL-SCH信道编码率，24表示CRC比特数如果R >= 1错误(['DL-SCH编码率('num2str(右)')大于或等于1的配置参数。']);结束

%当前配置的NPDSCH重复模式为%显示如下displayPattern = false;显示NPDSCH重复模式如果displayPattern == true npdschinfo . displaysubframpattern;结束

模拟起始点的绝对子帧数NSubframe = enb.NFrame*10+enb.NSubframe;初始化BLER和吞吐量结果maxThroughput = 0 (length(SNRdB)，1);simThroughput = 0(长度(SNRdB)，1);bler = 0(1，数字(SNRdB));临时变量'enb_init'和'channel_init'用于创建在信噪比循环中创建临时变量'enb'和'channel''parfor'循环的%独立模拟循环Enb_init = enb;Channel_init = channel;为snrIdx = 1:数字(SNRdB)% parfor snrIdx = 1:数值(SNRdB)允许使用并行计算，以提高注释输出的速度%上面的'for'语句，并取消下面的'parfor'语句的注释。这需要并行计算工具箱。如果没有安装% 'parfor'将默认为普通的'for'语句。根据循环变量设置随机数生成器种子%，确保独立的随机流rng (snrIdx“combRecursive”）;流(在%gdB信噪比\n'处模拟%d个传输块、numTrBlks SNRdB (snrIdx));Enb = enb_init;初始化eNodeB配置Channel = channel_init;初始化衰落信道配置TXCW = [];初始化传输的码字numBlkErrors = 0;出错传输块的百分比地产= [];初始化NPDSCH编码器状态Dstate = [];初始化NPDSCH解码器状态lastOffset = 0;初始化整体帧定时偏移量偏移量= 0;初始化帧定时偏移量subframeGrid = lteNBResourceGrid(enb);初始化子帧网格subframeIdx = NSubframe;numRxTrBlks = 0;而(numRxTrBlks < numTrBlks)设置当前子帧和帧号enb。NSubframe = mod(subframeIdx,10);enb。NFrame = floor((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 = randi([0 1]，trblklen,1);txcw = lteNDLSCH(rmoutlen,txTrBlk);结束如果(isDataSubframe)为子帧生成NPDSCH符号和索引[txNpdschSymbols,estate] = lteNPDSCH(enb,npdsch,txcw,estate);npdschIndices = lteNPDSCHIndices(enb,npdsch);将符号映射到子帧网格subframeGrid(npdschIndices) = txNpdschSymbols;生成NRS符号和索引nrsSymbols = lteNRS(enb);nrsIndices = lteNRSIndices(enb);将符号映射到子帧网格subframeGrid(nrsIndices) = nrsSymbols;结束%执行OFDM调制以生成时域波形[tx波形，ofdmInfo] = nbofdmmodulation (enb,subframeGrid);添加25个样本填充。这是为了覆盖延迟的范围%的预期从通道建模(组合%实现延迟和信道延迟扩展)tx波形= [tx波形;enb.NBRefP 0(25日)];% #好< AGROW >初始化每个子帧的通道时间通道。InitTime = subframeIdx/1000;通过通道模型传递数据通道。SamplingRate = ofdmInfo.SamplingRate;[rx波形，fadingInfo] = lteFadingChannel(通道，tx波形);计算噪声增益，包括下行功率补偿%分配信噪比= 10^(SNRdB(snrIdx)/20);%归一化噪声功率，以考虑采样率，其中%是用于OFDM调制的IFFT大小的函数，并且%天线数N0 = 1/(√(2.0*enb.NBRefP*double(ofdmInfo.Nfft))*信噪比);创建加性高斯白噪声噪声= N0*complex(randn(size(rx波形))，.．.randn(大小(rxWaveform)));将AWGN添加到接收的时域波形中rx波形= rx波形+噪声;%------------------------------------------------------------------%接收%------------------------------------------------------------------%执行定时同步，提取相应的%子帧的接收波形，并执行OFDM%解调如果(perfectChannelEstimator) offset = hPerfectTimingEstimate(fadingInfo);其他的在本例中，子帧偏移量计算依赖于%的NPSS出现在子帧5中，所以我们需要填充%子帧在它前面，使帧偏移量返回由% lteNBDLFrameOffset是子帧5的偏移量sfTsamples = ofdmInfo.SamplingRate*1e-3;如果(enb.NSubframe==5) padding = 0 ([sfTsamples*5,size(rx波形，2)]);offset = lteNBDLFrameOffset(enb，[填充;rxWaveform]);如果(offset > 25) || (offset < 0) offset = lastOffset;结束lastOffset = offset;结束结束同步接收到的波形rx波形= rx波形(1+offset:end，:);对接收到的数据进行OFDM解调以重建数据%资源网格rxSubframe = nbOFDMDemodulate(enb, rx波形);%信道估计如果(perfectChannelEstimator)%完美信道估计estChannelGrid = nbDLPerfectChannelEstimate(enb, channel, offset);noiseGrid = nbOFDMDemodulate(enb, noise(1+offset:end，:));noiseEst = var(noiseGrid(:));其他的[estChannelGrid, noiseEst] = lteDLChannelEstimate(.．.enb, cec, rxSubframe);结束如果(isDataSubframe)获取NPDSCH指数npdschIndices = lteNPDSCHIndices(enb, npdsch);从接收的子帧中获取PDSCH资源元素。规模的接收子帧百分比由PDSCH功率因数Rho。PDSCH是%按此量缩放，而单元格引用符号用于%信道估计(用于PDSCH解码阶段)则不是。[rxNpdschSymbols, npdschest] = lteExtractResources(npdschIndices，.．.rxSubframe estChannelGrid);解码NPDSCH[rxcw,dstate,symbols] = lteNPDSCHDecode(.．.enb, npdsch, rxNpdschSymbols, npdschest, noiseEst,dstate);解码传输块时，所有的子帧在一个捆绑%已收到如果dstate。EndOfTx [trblkout,blkerr] = lteNDLSCHDecode(trblklen,rxcw);numBlkErrors = numBlkErrors + blkerr;numRxTrBlks = numRxTrBlks + 1;重新初始化以启用新的传输块的传输TXCW = [];结束结束subframeIdx = subframeIdx + 1;结束计算块错误率bler(snrIdx) = numBlkErrors/numTrBlks;流(' npdsch bler = %。4 f \ n”提单(snrIdx));计算最大吞吐量和模拟吞吐量maxThroughput(snrIdx) = trblklen*numTrBlks;最大可能吞吐量simThroughput(snrIdx) = trblklen*(numTrBlks-numBlkErrors);%模拟吞吐量流('NPDSCH吞吐量(%%)= %。4 f % % \ n”simThroughput (snrIdx) * 100 / maxThroughput (snrIdx));结束

如果repIdx == 1 fh = figure;网格在；持有在；包含(“信噪比(dB)”）;ylabel (“提单”）;Legendstr = {['NRep = 'num2str (npdsch.NRep)]};其他的传奇str =[传奇str]'NRep = 'num2str (npdsch.NRep)];% #好< AGROW >结束图(跳频);情节(SNRdB,提单,“o”）;

结束%设置图形标题如果strcmpi (NPDSCHDataType“SIB1NB”) npdsch。NSF = 8;结束标题([' 'char (npdsch.NPDSCHDataType)”:TBS = 'num2str (trblklen).．.”;NSF = 'num2str (npdsch.NSF)“;”num2str (enb_init.NBRefP)“NRS端口”]);传奇(legendstr);

NB-IoT DL OFDM调制器函数[波形，信息]= nbofdmmodulation (enb,grid)根据TS 36.104表E.5.1-1a应用默认窗口大小如果(~ isfield (enb“窗口”enb))。窗口= 6;结束使用NB-IoT SC-FDMA在OFDM调制上获得1/2子载波移位enb。NBULSubcarrierSpacing =“15 khz”；[波形，信息]= lteSCFDMAModulate(enb,grid);结束NB-IoT DL OFDM解调器函数grid = nbOFDMDemodulate(enb, rx波形)使用NB-IoT SC-FDMA在OFDM调制上获得1/2子载波移位enb。NBULSubcarrierSpacing =“15 khz”；grid = lteSCFDMADemodulate(enb, rx波形，0.55);% CP分数为0.55结束% NB-IoT DL完美信道估计器函数H = nbDLPerfectChannelEstimate(enb,channel,timefreqoffset)重新配置NB-IoT UL完美的信道估计器，以执行完美的DL%信道估计enb。NBULSubcarrierSpacing =“15 khz”；enb。NTxAnts = enb.NBRefP;enb。TotSlots = 2;H = lteULPerfectChannelEstimate(enb, channel,timefreqoffset);结束