numTrBlks = 4;%模拟传输块数SNRDB = -32：4：0;%信噪比范围，单位为dB[0 5 9];%模拟的重复次数范围

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

安全部队= 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;％NR天线端口数量，应为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;eleesifstrcmpi（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结束

频道=结构;初始化通道配置结构通道。种子= 6;%通道种子通道。NRxAnts = 1;% 1接收天线通道。DelayProfile =“环保署”；%延迟概要channel.dopplerfreq = 5;Hz中％多普勒频率通道。MIMOCorrelation ='低的'；% Multi-antenna相关性通道。NTerms = 16;衰落模型中使用的振荡器通道。ModelType =“GMEDS”；%瑞利衰落模型类型channel.initphase =.“随机”；%随机初始阶段通道。NormalizePathGains =“上”；正规化延迟轮廓功率通道。NormalizeTxAnts =“上”；传送天线规格化

信道估计器行为PerfectChannElestimator = True;

配置信道估计器cec.pilotaverage =.“UserDefined”；导频符号平均值的类型cec。TimeWindow = 1;%时间窗口大小，以REs表示cec。FreqWindow = 25;RES中的％频率窗口大小cec。InterpType =“立方”；二维插值型cec。InterpWindow =“中心”；％插值窗口类型cec。InterpWinSize = 3;插值窗口大小cec。Reference =“关系”；信道估计器参考信号

为repidx = 1：numel（ireps）

npdschinfo = hnpdschinfo;npdschinfo.npdschdatatype = npdschdattype;npdschinfo.isf = ISF;如果Strcmpi（npdschdatype，“SIB1NB”）％npdsch携带sib1-nbnpdschinfo.schedulinginfosib1 = schedulinginfosib1;其他的% NPDSCH不携带SIB1-NBnpdschInfo。IRep= ireps(repIdx);DCI中的％重复编号字段（DCI格式N1或N2）npdschinfo.imcs = imcs;DCI中的％调制和编码方案字段（DCI格式N1或N2）％验证IREP和IMC的输入如果isempty（npdschinfo.tbs）npdschinfo.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;

[~,信息]= 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.displaySubframePattern;结束

%模拟起始点的绝对子帧数NSubframe = enb.NFrame * 10 + enb.NSubframe;％初始化BLER和吞吐量结果maxThroughput = 0(长度(SNRdB), 1);simThroughput = 0(长度(SNRdB), 1);提单= 0(1,元素个数(SNRdB));%临时变量“enb_init”和“channel_init”用于创建在SNR循环中创建临时变量“enb”和“channel”'parfor'循环的独立仿真循环enb_init = eNB;channel_init =频道;为snridx = 1：numel（snrdb）% parfor snrIdx = 1:numel(SNRdB)％要使用并行计算的使用，以增加速度评论%上面的'for'语句并取消下面的'parfor'语句的注释。%这需要并行计算工具箱。如果没有安装% 'parfor'将默认为普通的'for'语句。%根据循环变量设置随机数生成器种子%以确保独立的随机流rng (snrIdx“combRecursive”）;流(在%gdB信噪比下模拟%d传输块\n'，numtrblks，snrdb（snridx））;eNB = enb_init;%初始化eNodeB配置频道= channel_init;%初始化衰落信道配置txcw = [];％初始化传输的码字numBlkErrors = 0;百分比具有错误的传输块数房地产= [];%初始化NPDSCH编码器状态dstate = [];%初始化NPDSCH解码器状态lastOffset = 0;初始化整体帧定时偏移抵消= 0;初始化帧定时偏移subframeGrid = lteNBResourceGrid (enb);％初始化子帧网格subframeIdx = NSubframe;numRxTrBlks = 0;而(numRxTrBlks < numTrBlks)%设置当前子帧和帧号enb。NSubframe =国防部(subframeIdx 10);enb。NFrame =地板((subframeIdx) / 10);%生成NPSS符号和索引npssSymbols = lteNPSS (enb);npssIndices = lteNPSSIndices (enb);将符号映射到子帧网格subframeGrid (npssIndices) = npssSymbols;％生成NSSS符号和索引nssssymbols = ltenss（eNB）;nsssindices = ltenssindices（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,房地产]= lteNPDSCH (txcw, enb, npdsch房地产);npdschIndices = lteNPDSCHIndices (enb npdsch);将符号映射到子帧网格subframeGrid (npdschIndices) = txNpdschSymbols;%生成NRS符号和索引nrsSymbols = lteNRS (enb);nrsIndices = lteNRSIndices (enb);将符号映射到子帧网格subframegrid（nrsindices）= nrssymbols;结束%进行OFDM调制以产生时域波形[txWaveform, ofdmInfo] = nbOFDMModulate (enb subframeGrid);%添加25个样品填充。这是为了覆盖延迟的范围%期望从渠道建模(组合实现延迟和信道延迟扩展%)txWaveform = [txWaveform;enb.NBRefP 0(25日)];% #好< AGROW >%初始化每个子帧的信道时间Channel.Inittime = SubframeDX / 1000;％通过频道模型传递数据通道。SamplingRate = ofdmInfo.SamplingRate;[rx波形，fadingInfo] = lteFadingChannel(channel, tx波形);计算噪声增益，包括下行功率的补偿%分配SNR = 10 ^（SNRDB（SNRIDX）/ 20）;将噪声功率归一化，以考虑采样率%是OFDM调制中使用的IFFT大小的函数%天线数量N0 = 1 /(√(2.0 * enb.NBRefP *双(ofdmInfo.Nfft)) *信噪比);创建加性高斯白噪声噪声= n0 *复杂（Randn（randn（rxwaveform）），......Randn（尺寸（rxwaveform））））;%将AWGN添加到接收的时域波形中rxwaveform = rxwaveform +噪声;%------------------------------------------------------------------%接收%------------------------------------------------------------------%执行定时同步，提取适当的收到波形的％子帧，并执行OFDM%解调如果(perfectChannelEstimator) offset = hPerfectTimingEstimate(fadingInfo);其他的在这个例子中，子帧偏移计算依赖于% on NPSS present in subframe 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 =抵消;结束结束%同步接收的波形rx波形= rx波形(1+偏移:结束，:);%对接收到的数据进行OFDM解调，重新创建%资源网格rxSubframe = nbOFDMDemodulate (enb rxWaveform);％信道估计数如果（PerfectChanneLestimator）完美信道估计百分比EstchannelGrid = nbdlperfectChanneLestimate（eNB，频道，偏移）;kotidgrid = nbofdmdemodulate（eNB，噪声（1 +偏移：结束，:)）;destemest = var（诊断（:)）;其他的[estChannelGrid, noiseEst] = lteDLChannelEstimate(......cec, enb 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, 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));结束

如果reidx == 1 fh = figure;网格在；持有在；Xlabel（“信噪比(dB)”）;ylabel（'bler'）;legendstr = {('nrep ='num2str（npdsch.nrep）]};其他的Legendstr = [letcendstr ['nrep ='num2str（npdsch.nrep）]];% #好< AGROW >结束图(跳频);情节(SNRdB,提单,“o”）;

结束设置数字标题如果Strcmpi（npdschdatype，“SIB1NB”) npdsch。NSF = 8;结束标题([''char (npdsch.NPDSCHDataType)'：TBS ='num2str (trblklen)......”;NSF = 'num2str (npdsch.NSF)“;”num2str（enb_init.nbrefp）'NRS端口']);传奇(legendstr);

% NB-IoT - DL OFDM调制器函数[波形，信息] = nbofdmodulate（eNB，网格）%应用默认窗口大小根据TS 36.104表E.5.1-1a如果(~ isfield (enb“窗口”）））eNB.Windowing = 6;结束%使用NB-IoT SC-FDMA在OFDM调制上获得1/2副载波移位enb。NBULSubcarrierSpacing ='15khz'；(波形信息]= lteSCFDMAModulate (enb、网格);结束% NB-IoT - DL OFDM解调器函数网格= nbOFDMDemodulate (enb rxWaveform)%使用NB-IoT SC-FDMA在OFDM调制上获得1/2副载波移位enb。NBULSubcarrierSpacing ='15khz'；Grid = LTESCFDMADEMODULES（eNB，RxWaveForm，0.55）;% CP含量为0.55结束% NB-IoT DL完美信道估计器函数h = nbdlperfectChanneLestimate（eNB，频道，timefreqoffset）%重新配置NB-IoT UL完美信道估计器，以执行DL完美%信道估计enb。NBULSubcarrierSpacing ='15khz'；enb。NTxAnts = enb.NBRefP;enb。TotSlots = 2;H = lteULPerfectChannelEstimate(enb, channel,timefreqoffset);结束