numTrBlks = 4;%的模拟传输块SNRdB = 32:4:0;在dB %信噪比范围ireps = [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;%的数量关系,天线端口,应该是1或2enb.OperationMode=“Inband-DifferentPCI”;%允许的值是‘Inband-SamePCI’,‘Inband-DifferentPCI’,‘Guardband’或‘独立’如果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结束

频道=结构;%初始化通道配置结构通道。种子= 6;%通道种子通道。NRxAnts = 1;% 1接收天线通道。DelayProfile =“环保署”;%延迟概要通道。DopplerFreq = 5;在赫兹%多普勒频率通道。MIMOCorrelation =“低”;% Multi-antenna相关性通道。NTerms = 16;%振荡器用于衰落模型通道。ModelType =“GMEDS”;%瑞利衰落模型类型通道。InitPhase =“随机”;%随机初始阶段通道。NormalizePathGains =“上”;%正常化延迟概要文件的权力通道。NormalizeTxAnts =“上”;%正常化传输天线

%信道估计的行为perfectChannelEstimator = true;

%配置信道估计cec。PilotAverage =“UserDefined”;%的试点,平均象征cec。TimeWindow = 1;在REs %时间窗口大小cec。FreqWindow = 25;在REs %频率窗口大小cec。InterpType =“立方”;% 2 d插值类型cec。InterpWindow =“中心”;%插值窗口类型cec。InterpWinSize = 3;%插值窗口大小cec。Reference =“关系”;%信道估计的参考信号

为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。imc= IMCS;%调制和编码方案在DCI (DCI格式N1和N2)%验证IRep和imc的输入如果npdschInfo isempty (npdschInfo.TBS)。TBSTable错误([“无效(髂胫,安全部队)(髂胫= imc = 'num2str (imc)…”,安全部队= 'num2str(安全部队)”),返回空TBS,检查有效的对在上面的表或3 gpp 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;%初始化提单和吞吐量的结果maxThroughput = 0(长度(SNRdB), 1);simThroughput = 0(长度(SNRdB), 1);提单= 0(1,元素个数(SNRdB));%的临时变量“enb_init”和“channel_init”是用于创建%的临时变量' enb '和信噪比回路中的“通道”来创建%独立模拟循环parfor的循环enb_init = enb;channel_init =通道;为snrIdx = 1:元素个数(SNRdB)% parfor snrIdx = 1:元素个数(SNRdB)%,使并行计算的使用增加速度注释掉%以上”的声明和取消下面的“parfor”声明。%这需要并行计算工具箱。如果这是没有安装%’parfor将默认为正常”的声明。%设置随机数生成器的种子不同循环变量%,以确保独立随机流rng (snrIdx“combRecursive”);流(在% \ nSimulating % d运输块gdB信噪比\ 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= mod(subframeIdx,10); enb.NFrame = floor((subframeIdx)/10);%生成nps符号和指数npssSymbols = lteNPSS (enb);npssIndices = lteNPSSIndices (enb);%符号映射到子帧网格subframeGrid (npssIndices) = npssSymbols;%生成、符号和指数nsssSymbols = lteNSSS (enb);nsssIndices = lteNSSSIndices (enb);%符号映射到子帧网格subframeGrid (nsssIndices) = nsssSymbols;%建立如果nps或者、传播,如果是这样,%不传输NPDSCH子帧isDataSubframe = isempty (npssSymbols) & & isempty (nsssSymbols);%时创建一个新的交通阻塞和编码%传输码字是空的。接收方设置码字%空信号,所有的子帧一个包%接收第一传输之前(这也是空)如果isempty (txcw) txTrBlk =兰迪([0 1]trblklen 1);txcw = lteNDLSCH (rmoutlen txTrBlk);结束如果(isDataSubframe)%为子帧生成NPDSCH符号和指数[txNpdschSymbols,房地产]= lteNPDSCH (txcw, enb, npdsch房地产);npdschIndices = lteNPDSCHIndices (enb npdsch);%符号映射到子帧网格subframeGrid (npdschIndices) = txNpdschSymbols;%生成关系,符号和指数nrsSymbols = lteNRS (enb);nrsIndices = lteNRSIndices (enb);%符号映射到子帧网格subframeGrid (nrsIndices) = nrsSymbols;结束%完成OFDM调制生成时域波形[txWaveform, ofdmInfo] = nbOFDMModulate (enb subframeGrid);% 25样品添加填充。这是延迟的范围%预期从信道建模(的组合%实现延迟和延迟扩展频道)txWaveform = [txWaveform;enb.NBRefP 0(25日)];% #好< AGROW >%初始化通道为每个子帧时间通道。InitTime = subframeIdx / 1000;%通过通道模型传递数据通道。SamplingRate = ofdmInfo.SamplingRate;[rxWaveform, fadingInfo] = lteFadingChannel(频道,txWaveform);%计算噪声增益包括下行功率补偿%分配信噪比= 10 ^ (SNRdB (snrIdx) / 20);%正常化噪声功率的采样率,%是传输线大小的函数用于OFDM调制、和%的天线数量N0 = 1 / (√(2.0 * enb.NBRefP *双(ofdmInfo.Nfft)) *信噪比);%建立加性高斯白噪声噪音= N0 *复杂(randn(大小(rxWaveform)),…randn(大小(rxWaveform)));%添加情况下接收到的时域波形rxWaveform = rxWaveform +噪声;% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -%接收% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -%完成定时同步,提取适当的%的子帧接收到的波形,并执行OFDM%解调如果(perfectChannelEstimator)抵消= hPerfectTimingEstimate (fadingInfo);其他的%在这个例子中,子帧偏移量的计算% nps在子帧5,所以我们需要垫%子帧之前返回的帧偏移% lteNBDLFrameOffset 5是子帧的偏移sfTsamples = ofdmInfo.SamplingRate * 1 e - 3;如果(enb.NSubframe = = 5)填充= 0 ([sfTsamples * 5、大小(rxWaveform, 2)));抵消= lteNBDLFrameOffset (enb[填充;rxWaveform]);如果(抵消> 25)| |(抵消< 0)抵消= lastOffset;结束lastOffset =抵消;结束结束%同步接收到的波形rxWaveform = rxWaveform(1 +抵消:最终,);%完成OFDM解调接收的数据重新创建%资源网格rxSubframe = nbOFDMDemodulate (enb rxWaveform);%信道估计如果(perfectChannelEstimator)%完美的信道估计estChannelGrid = nbDLPerfectChannelEstimate (enb,通道,抵消);noiseGrid = nbOFDMDemodulate (enb、噪音(1 +抵消:最终,));噪音= var (noiseGrid (:));其他的[estChannelGrid,噪音]= lteDLChannelEstimate (…cec, enb rxSubframe);结束如果(isDataSubframe)%得到NPDSCH指数npdschIndices = lteNPDSCHIndices (enb npdsch);%得到PDSCH资源元素从接收到的子帧。规模的%接收子帧PDSCH功率因数ρ。PDSCH是%比例的金额,而细胞符号用于参考%信道估计(用于PDSCH解码阶段)。[rxNpdschSymbols, npdschHest] = lteExtractResources (npdschIndices,…rxSubframe estChannelGrid);%解码NPDSCH[rxcw dstate,符号]= lteNPDSCHDecode (…npdsch, enb rxNpdschSymbols npdschHest,噪音,dstate);%解码传输块当所有的子帧包%已收到如果dstate。EndOfTx [trblkout blkerr] = lteNDLSCHDecode (trblklen rxcw);numBlkErrors = numBlkErrors + blkerr;numRxTrBlks = numRxTrBlks + 1;%初始化,使一个新的传输块的传播txcw = [];结束结束subframeIdx = subframeIdx + 1;结束%计算块错误率提单(snrIdx) = numBlkErrors / numTrBlks;流(“NPDSCH提单= %。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 =图;网格在;持有在;包含(“信噪比(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端口(s)”]);传奇(legendstr);

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