tbs = 144;％传输块大小totnumblks = 1;%模拟传输块数问题=结构();%初始化UE结构UE.nbulsubcarierspacing ='15khz'；％3.75khz，15khzUE.nncellid = 0;％窄带细胞标识chs = struct（）;％NPUSCH带有数据或控制信息chs.npuschformat ='数据'；%负载类型(数据或控制)用于NPUSCH ' nscu '的子载波数量取决于NPUSCH％格式和子载波间隔'nbulsubcarierspacing'如ts所示％36.211表10.1.2.3-1。有1,3,6或12个连续的子载波％npusch.chs.nbulsubcarierset = 0;％范围为0-11（15kHz）;0-47（3.75kHz）chs.nrusc =长度（chs.nbulsubcarrierset）;chs.cyclicshift = 0;当NRUsc = 3或6时需要的循环移位%chs.rnti = 0;% RNTI价值chs.nlayers = 1;％的层数chs。NRU= 2;%资源单位数chs。NRep = 4;％NPUSCH的重复次数Chs.slotidx = 0;％启动捆绑中的插槽索引％符号调制取决于NPUSCH格式和NSCRU由TS 36.211给出的％表10.1.3.2-1chs。调制='QPSK'；rvDCI = 0;通过DCI信号的% RV偏移量(见36.213 16.5.1.2)％在DB中指定NPUSCH和DM-RS功率缩放，用于绘图可视化chs.npuschpower = 30;chs.npuschdrspower = 34;

Chs.seqgrouphopping ='在'；％启用/禁用序列组跳跃的UEchs.seqgroup = 0;％delta_ss。高层参数Groupassignmentnpusch％获取资源单元纽尔斯中的时间插槽数量％TS 36.211表10.1.2.3-1如果Strcmpi（Chs.npuschformat，'数据'）如果chs.nrusc == 1 nulslots = 16;eleesif任何(chs。NRUsc == [3 6 12]) NULSlots = 24/chs.NRUsc;别的错误（'子载波的数量无效。NRUSC必须是1,3,6,12'中的一个）;结尾eleesifStrcmpi（Chs.npuschformat，'控制'）Nulslots = 4;别的错误（'无效的NPUSCH格式(%s)。NPUSCHFormat必须是" Data "或" Control "，chs.npuschformat）;结尾chs。NULSlots = NULSlots;NSlotsPerBundle = chs.NRU * chs.NULSlots * chs.NRep;码字捆绑中的次数totnslots = totnumblks * nslotsperbundle;模拟插槽的％总数

％将随机生成器初始化为默认状态RNG（'默认'）;％获取每个框架的插槽网格和插槽数emptySlotGrid = lteNBResourceGrid(问题);slotGridSize =大小(emptySlotGrid);NSlotsPerFrame = 20 / (slotGridSize (1) / 12);状态= [];％NPUSCH编码器和DM-RS状态，在功能中自动重新初始化trblk = [];％初始化传输块txgrid = [];全网格初始化％显示正在生成的时隙数fprintf（'\n生成%d个槽，对应%d个传输块(s)\n'，totnslots，totnumblks）;为了slotidx = 0+（0：totnslots-1）％计算帧内的帧号和插槽号UE.NFRAME =修复（SLOTIDX / NSLOTSPERFREME）;ue.nslot = mod（slotidx，nslotsperframe）;如果isempty（trblk）如果Strcmpi（Chs.npuschformat，'数据'）对于用于的两个RV值来说，完成％UL-SCH编码％发送码字。使用的RV序列确定从DCI中发出的RVDCI值的％并交替在TS 36.213第16.5.1.2节中给出的0和2之间的％％定义将被编码以创建的传输块不同RV的％码字trblk = randi（[0 1]，TBS，1）;％确定编码传输块大小[〜，INFO] = LTEPUSCHINDICES（UE，CHS）;outblklen = info.g;％创建与使用的两个RV值对应的码字％在第一和第二块中，这将重复直到所有传输％块chs.rv = 2 * mod（rvdci + 0,2）;第一个块的％RV连续波= lteNULSCH (chs、outblklen trblk);% CRC和Turbo编码chs.rv = 2 * mod（rvdci + 1,2）;第二块的％RVCW = [CW LTTULSCH（CHS，OUTBLKLEN，TRBLK）];％＃OK ％CRC和Turbo编码重复别的trblk = randi（[0 1]，1）;％1位ACK%对于ACK，每个块传输相同的码字TS 36.212第6.3.3节中定义的％cw = ltenulsch（trblk）;结尾blockIdx = 0;%要传输的第一个块结尾％初始化网格slotgrid = emptyslotgrid;％NPUSCH编码和映射到插槽网格Txsym = Ltenusch（UE，CHS，CW（：，mod（blockidx，尺寸（cw，2））+ 1），状态）;％映射插槽网格中的NPUSCH符号indicsnpusch = ltenuschindices（Ue，chs）;SlotGrid（indicsnpusch）= txsym * db2mag（chs.npuschpower）;%创建DM-RS序列并映射到槽网格(dmr、州)= lteNPUSCHDRS(问题、chs、状态);chs indicesDMRS = lteNPUSCHDRSIndices(问题);slotGrid (indicesDMRS) = dmr * db2mag (chs.NPUSCHDRSPower);％将此插槽连接到插槽网格txgrid = [txgrid slotgrid];％＃OK ％如果传输完整块，则递增时钟计数器％可以选择正确的码字如果state.endofblk blockidx = blockidx + 1;结尾% TRBLK错误计数和重新初始化如果状态。EndOfTx％重新初始化以启用新传输块的传输trblk = [];结尾结尾％执行SC-FDMA调制以创建时域波形UE.CyclicPrefixul ='普通的'；NB-IOT的％常规循环前缀长度[波形，scfdmainfo] = LTESCFDMAMODULES（UE，CHS，TxGrid）;

％创建了整个资源网格的图像图im = image(abs(txgrid));提出= parula (64);colormap (im.Parent提出);轴xy；标题(sprintf ('NB-IOT上行链路RE网格（NREP =％D，NRUSC =％D，NRU =％d）', chs.NRUsc chs.NRep chs.NRU)包含('OFDM符号'）ylabel（'子载波'）%创建图例框以指示与REs关联的通道/信号类型重命名= {'NPUSCH'；“DM-RS”};CLEVELS = ROUND（DB2MAG（[CHS.NPUSCHPOWER CHS.NPUSCHDRSPOWER]））;n = numel（重命名）;l =行（那些（n），一个（n），“线宽”，8）;％生成线条％根据CMAP设置颜色集(L, {'颜色'}，mat2cell（cmap（min（min（1 + clevels，长度（cmap）），:)，（1，n），3））;传奇（重命名{：}）;