TBS = 144;%传输块大小totnumblks = 1;百分比的模拟传输块UE = struct（）;％初始化UE结构UE.nbulsubcarierspacing =“15 khz”;% 3.75 khz, 15 khzUE.nncellid = 0;窄带单元识别chs = struct（）;％NPUSCH带有数据或控制信息chs。NPUSCHFormat ='数据';％有效载荷类型（数据或控制）用于NPUSCH ' nscu '的子载波数量取决于NPUSCH％格式和子载波间隔'nbulsubcarierspacing'如ts所示表10.1.2.3-1。有1、3、6或12个连续的子载波％npusch.chs.nbulsubcarierset.= 0;%量程0-11 (15kHz);0-47 (3.75 khz)chs.nrusc.= length(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表10.1.3.2-1chs.modulation =.“正交相移编码”;rvdci = 0;％RV偏移通过DCI信号信号（参见36.213 16.5.1.2）%指定NPUSCH和DM-RS功率缩放在dB中用于绘图可视化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;elseif任何（chs.nrusc == [3 6 12]）Nulslots = 24 / chs.nrusc;别的错误('子载波数量无效。NRUsc必须是1 3 6 12'中的一个）;结尾elseifStrcmpi（Chs.npuschformat，'控制'NULSlots = 4;别的错误('无效的NPUSCH格式（％s）。npuschformat必须是''数据''或''''''''，chs.npuschformat）;结尾chs.nulslots =鼻子;nslotsperbundle = chs.nru * chs.nulslots * chs.nrep;%码字包中的槽数totnslots = totnumblks * nslotsperbundle;模拟插槽的％总数

％将随机生成器初始化为默认状态RNG（'默认'）;％获取每个框架的插槽网格和插槽数emptyslotgrid = ltenbresourcegrid（UE）;slotgridsize = size（emptyslotgrid）;nslotsperframe = 20 /（SlotGridSize（1）/ 12）;国家= [];％NPUSCH编码器和DM-RS状态，在功能中自动重新初始化trblk = [];％初始化传输块txgrid = [];％全网格初始化%显示正在生成的槽位数fprintf（'\ ngeneration％d槽对应于％d传输块（s）\ n'、TotNSlots totNumBlks);为了slotidx = 0+（0：totnslots-1）％计算帧内的帧号和插槽号问题。NFrame =修复(slotIdx / NSlotsPerFrame);问题。NSlot =国防部(slotIdx NSlotsPerFrame);如果isempty（trblk）如果Strcmpi（Chs.npuschformat，'数据'）对于用于的两个RV值来说，完成％UL-SCH编码％发送码字。使用的RV序列确定%从rvDCI值在DCI和交替信号在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。房车= 2 *国防部(rvDCI + 0, 2);第一个块的％RVCw = Ltenulsch（Chs，Outblklen，Trblk）;％CRC和Turbo编码chs。房车= 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符号chs indicesNPUSCH = lteNPUSCHIndices(问题);slotGrid (indicesNPUSCH) = txsym * db2mag (chs.NPUSCHPower);％创建DM-RS序列并映射到插槽网格[DMRS，状态] = LTEPUSCHDRS（UE，CHS，状态）;indicsdmrs = ltenuschdrsindices（UE，chs）;SlotGrid（indicsdmrs）= dmrs * db2mag（chs.npuschdrspower）;％将此插槽连接到插槽网格txgrid = [txgrid slotGrid];％＃OK ％如果传输完整块，则递增时钟计数器％可以选择正确的码字如果状态。endofblock blockIdx = blockIdx + 1;结尾％trblk err count和重新初始化如果State.endoftX.%重新初始化以启用新传输块的传输trblk = [];结尾结尾％执行SC-FDMA调制以创建时域波形UE.CyclicPrefixul ='普通的';NB-IOT的％常规循环前缀长度[波形，scfdmainfo] = LTESCFDMAMODULES（UE，CHS，TxGrid）;

％创建了整个资源网格的图像图IM =图像（ABS（TxGrid））;CMAP = Parula（64）;Colormap（IM.Parent，CMAP）;轴XY.;标题（Sprintf（'NB-IOT上行链路RE网格（NREP =％D，NRUSC =％D，NRU =％d）'，chs.nrep，chs.nrusc，chs.nru））xlabel（'OFDM符号') ylabel (副载波的）％创建图例框以指示与res关联的频道/信号类型重命名= {'NPUSCH';'dm-rs'};CLEVELS = ROUND（DB2MAG（[CHS.NPUSCHPOWER CHS.NPUSCHDRSPOWER]））;n = numel（重命名）;l =行（那些（n），一个（n），'行宽'，8）;%生成线％根据CMAP设置颜色set（l，{“颜色”}，mat2cell（cmap（min（min（1 + clevels，长度（cmap）），:)，（1，n），3））;传奇（重命名{：}）;