％初始化nbiotdownlinkwaveformgenerator对象％％使用'r.nb.5'和'r.nb.6'的'rc'值进行MIB解码。注意％rmcs'r.nb.5-1'，'r.nb.6-1'和'r.nb.7'对应于非锚型％不包含所需的NPS，NSS和NPBCH的载体。rc ='r.nb.5'；ngen = nbiotdownlinkwaveformgenerator（RC）;ngen.config.nncellid = 120;ngen.config.nbrefp = 1;ngen.config.nframe = 640;ngen.config.totsubframes = 50;ngen.config.npbch.datasource ='mib'；%生成时域窄带物联网波形[eNodeBoutput，网格，OFDMINFO] =生成WaveForm（ngen）;sr = ofdminfo.samplingrate;％引入延迟以模拟相关的未知定时对齐百分比是任意信号捕获。这在接收器处得到补偿％同时同步接收信号。延迟= 50;波形= circshift(eNodeBOutput, delay);％初始化图如果（〜存在（'ChannelFigure'，'var'）||〜IsValid（ChannelFigure））ChannEvigure = Figure（'可见的'，'离开'）;结尾[Spectrumanalyzer，Synchcorlplot，PdcchConstdiagram] =.．.hSIB1RecoveryExamplePlots (channelFigure sr);％显示接收信号频谱流('\ npplotting接收信号谱... \ n'）;简介(awgn(波形,100.0));

%细胞搜索流('\ n informing小区搜索... \ n'）;alg。SSSDetection ='postfft'；alg.maxcellcount = 1;[nncellid，偏移，峰] = hnbcellsearch（波形，all）;nbenb = struct;nbenb.nncellid = nncellid;[〜，corr] = ltenbdlframeoffset（setfield（nbenb，'操作模式'，“独立”），波形）;% #好< SFLD >％作为NB-IOT具有与LTE波形相似的波形结构，使用功能与LTE有关的％与NB-IOT一起使用。为此，创建一个设置结构enb。enb = nbenb;％plot npss / nsss相关和阈值Synchcorlplot.ylimits = [0 max（corr（：））* 1.1];同步载推（Corr）;％执行时序同步流(“时间偏移到帧开始：％d样本\ n',抵消);timesynced =波形(1 +抵消:最终,);%显示单元格范围的设置流(单元格搜索后的单元格范围设置:\n）;DISP（NBENB）;

频率偏移估计流('\n执行频率偏移估计…\n'）;enb。NULRB = 6;频率偏移计算的％enb。DuplexMode =“FDD”；foffset = ltefrequencyOffset（eNB，倍增，0）;流('频率偏移：％0.3FHz \ n', foffset);%补偿频率偏差freqsynced = ltefrequencycorect（eNB，乘时，foffset）;

％信道估计器配置cec。PilotAverage ='用户自定义'；％类型的飞行员平均cec.freqwindow = 13;％频率窗口大小cec。TimeWindow = 9;％时间窗口大小cec.interpype ='立方体'；％2D插值类型cec.interpwindow ='中心'；%插值窗口类型cec.Interpwinsize = 1;％插值窗口大小cec.reference =.'nrs'；％NB-IOT下行链路信道估计％假设4个小区特定和2个窄带参考信号，用于初始％解码尝试。这确保了所有可用的通道估计%的参考信号enb.cellrefp = 4;enb.nbrefp = 2;eNB.ntxants = 2;nbenb.nbrefp = 2;流('执行OFDM解调...... \ n \ n'）;griddims = lteresourcegridsize（eNB）;％资源网格尺寸l = griddims（2）;子帧中的％OFDM符号数enb.nbulsubcarierspacing =“15 khz”；NB-IOT下行链路的％需要DDM解调rxgrid = LTESCFDMADEMODULES（eNB，Freqsynced）;NB-IOT中OFDM解调的百分比如果（isempty（rxgrid））fprintf（“定时同步后，信号短于一个子帧，因此不会执行进一步的解调。\ n'）;返回；结尾%进行信道估计enb。NSubframe = 0;％初始化nsubframe以进行信道估计[hest, nest] = lteDLChannelEstimate(enb, cec, rxgrid(:，1:L，:));

％解码MIB%从第一个NPBCH中提取对应的资源元素(REs)所有接收天线的％子帧和信道估计分隔符= repmat (“- - -”，1,50）;流('％s \ n'，分隔器）;流(“执行MIB解码…\ n”）;流('％s \ n \ n'，分隔器）;npbchindices = ltenbchindices（nbenb）;[npbchrx，npbchht] = lteextractresouse（npbchindices，.．.rxgrid（1：12,1：l，:)，hest（：，1：l：，:)）;%解码NPBCHdstate = [];[bchBits, dstate npbchSymbols nfmod64, mib, nbenb。NBRefP] = lteNPBCHDecode (.．.NBENB，NPBCHRX，NPBCHHEST，NEST，DSTATE）;％parse mib bitsnbenb nbenb = hNBMIB (mib);如果（nbenb.nbrefp == 0）fprintf（'MIB解码失败（nbenb.nbrefp = 0）。\ n \ n'）;返回；结尾%将lteNPBCHDecode函数中的nfmod64值合并为％从MIB建立的NFRAME值是SFN Modulo 64（它是%作为floor(SFN/64)存储在MIB中。如果nbenb.nbrefp.如果传递了NPBCH解码，则会设置nframe值nbenb.nframe = nbenb.nframe + nfmod64;结尾如果Strcmpi（nbenb.operationMode，“Inband-SamePCI”) nbenb。NCellID = NNCellID;结尾% MIB解码后显示整个单元的设置流('MIB解码后的单元格设置：\ n'）;DISP（NBENB）;