numTrBlks = 4;%模拟传输块的数量SNRdB = -32:4:0;%信噪比范围，单位为dBirep = [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物理cell 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。nncellid = 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结束

Channel = 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;REs中的时间窗口大小cec。FreqWindow = 25;% REs中的频率窗口大小cec。InterpType =“立方”；% 2D插值类型cec。InterpWindow =“中心”；插值窗口类型cec。InterpWinSize = 3;插值窗口大小cec。Reference =“关系”；信道估计参考信号

为repIdx = 1:numel(ireps)

npdschInfo = hNPDSCHInfo;npdschInfo。NPDSCHDataType = NPDSCHDataType;npdschInfo。Isf = Isf;如果strcmpi (NPDSCHDataType“SIB1NB”）% NPDSCH携带SIB1-NBnpdschInfo。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 error(['DL-SCH编码率('num2str(右)')大于或等于1的配置参数。']);结束

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

%模拟开始时的绝对子帧数NSubframe = enb.NFrame*10+enb.NSubframe;初始化BLER和吞吐量结果maxThroughput = 0 (length(SNRdB)，1);simThroughput = 0 (length(SNRdB)，1);bler = 0 (1,numel(SNRdB));临时变量'enb_init'和'channel_init'用于创建%临时变量'enb'和'channel'在信噪比循环内创建%独立的模拟循环为'parfor'循环Enb_init = enb;Channel_init = channel;为snrIdx = 1:数值(SNRdB)SNRdB = 1:数值(SNRdB)启用并行计算以提高注释输出的速度%删除上面的'for'语句并取消注释下面的'parfor'语句。这需要并行计算工具箱。。如果没有安装% 'parfor'将默认为正常的'for'语句。根据循环变量设置随机数生成器种子%以确保独立的随机流rng (snrIdx“combRecursive”）;流('\n模拟%d传输块在%gdB信噪比\n'、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);%时创建一个新的传输块并对其进行编码%传输码字为空。接收方设置码字%为空，表示bundle中的所有子帧都已被替换%已接收(在第一次传输之前也是空的)如果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);将符号映射到子帧网格subframgrid (npdschIndices) = txNpdschSymbols;生成NRS符号和索引nrsSymbols = lteNRS(enb);nrsIndices = lteNRSIndices(enb);将符号映射到子帧网格subframeGrid(nrsIndices) = nrsSymbols;结束执行OFDM调制生成时域波形[tx波形，ofdmInfo] = nbofdmmodululate (enb, subframgrid);添加25个样本填充。这是为了涵盖延迟的范围%期望从通道建模(的组合%实现延迟和通道延迟扩散)tx波形= [tx波形;enb.NBRefP 0(25日)];% #好< AGROW >初始化每个子帧的通道时间通道。InitTime = subframeIdx/1000;通过通道模型传递数据通道。SamplingRate = ofdmInfo.SamplingRate;[rx波形，fadingInfo] = lteFadingChannel(channel, tx波形);计算包括下行功率补偿在内的噪声增益%分配SNR = 10^(SNRdB(snrIdx)/10);%归一化噪声功率以考虑抽样率，其中%是OFDM调制中使用的IFFT大小的函数，和%天线数N0 = 1/√(2.0*enb.NBRefP*double(ofdmInfo.Nfft)*SNR);创建加性高斯白噪声噪声= N0*complex(randn(size(rx波形))，.．.randn(大小(rxWaveform)));在接收的时域波形中添加AWGNrx波形= rx波形+噪声;%------------------------------------------------------------------%接收%------------------------------------------------------------------%执行定时同步，提取适当的%子帧的接收波形，并执行OFDM%解调如果(perfectChannelEstimator) offset = hPerfectTimingEstimate(fadingInfo);其他的在本例中，子帧偏移量计算依赖于NPSS的%出现在子帧5中，所以我们需要填充%子帧之前，使帧偏移量返回% lteNBDLFrameOffset是子帧5的偏移量sfsamples = ofdmInfo.SamplingRate*1e-3;如果(enb.NSubframe==5) padding = zero ([sfTsamples*5,size(rx波形，2)]);offset = lteNBDLFrameOffset(enb， [padding;]rxWaveform]);如果(offset > 25) || (offset < 0) offset = lastOffset;结束lastOffset =偏移量;结束结束同步接收到的波形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功率因数。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)]};其他的Legendstr = [Legendstr]'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);

物联网DL OFDM调制器函数[波形，信息]= nbofdmmodululate (enb,grid)根据TS 36.104表E.5.1-1a应用默认窗口大小如果(~ isfield (enb“窗口”enb))。Windowing = 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);结束