tbs = 144;％传输块大小totNumBlks = 1;%模拟传输块数问题=结构();%初始化UE结构问题。NBULSubcarrierSpacing ='15khz'；％3.75khz，15khz问题。NNCellID = 0;％窄带细胞标识chs =结构();% NPUSCH携带数据或控制信息chs.npuschformat =“数据”；%负载类型(数据或控制)用于NPUSCH ' nscu '的子载波数量取决于NPUSCH%格式和子载波间距'NBULSubcarrierSpacing'，如TS所示％36.211表10.1.2.3-1。有1,3,6或12个连续的子载波% NPUSCHchs.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;%在一个bundle中启动插槽索引符号调制取决于NPUSCH格式和nscu由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 =“上”；%启用/禁用终端顺序组跳变功能chs。SeqGroup = 0;% Delta_SS。更高层次参数groupAssignmentNPUSCH%获取资源单元的时间槽数表10.1.2.3-1如果strcmpi (chs。NPUSCH.Format,“数据”）如果chs.nrusc == 1 nulslots = 16;eleesif.任何(chs。NRUsc== [3 6 12]) NULSlots = 24/chs.NRUsc;其他的错误（'子载波的数量无效。NRUSC必须是1,3,6,12'中的一个）;结束eleesif.strcmpi (chs。NPUSCH.Format,“控制”）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 = [];全网格初始化％显示正在生成的时隙数流('\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。NPUSCH.Format,“数据”）% UL-SCH对所使用的两个RV值进行编码%传送密码。确定使用的RV序列从DCI中发出的RVDCI值的％并交替%在TS 36.213第16.5.1.2节中给出的0至2之间定义传输块，它将被编码以创建不同RV的%码字TRBLK = randi([0 1]，tbs,1);确定编码传输块大小[~， info] = lteNPUSCHIndices(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）;% RV为第二个区块cw = [cw lteNULSCH(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 = lteNPUSCH (chs问题,连续波(:,国防部(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];% #好< AGROW >%如果传输的是一个完整的块，则增加时钟计数器%可以选择正确的码字如果state.endofblk blockidx = blockidx + 1;结束% TRBLK错误计数和重新初始化如果状态。EndOfTx％重新初始化以启用新传输块的传输trblk = [];结束结束%执行SC-FDMA调制来创建时域波形问题。CyclicPrefixUL =“正常”；%窄带物联网正常循环前缀长度[波形,scfdmaInfo] = lteSCFDMAModulate (chs退居其次,txgrid);

创建一个整体资源网格的映像图im = image(abs(txgrid));提出= parula (64);colormap (im.Parent提出);轴xy；标题(sprintf ('NB-IoT Uplink RE Grid (NRep = %d, NRUsc = %d, NRU = %d)', chs.NRUsc chs.NRep chs.NRU)包含(OFDM符号的）ylabel（'子载波'）%创建图例框以指示与REs关联的通道/信号类型重命名= {“NPUSCH”；“DM-RS”};clevels =圆(db2mag ([chs。NPUSCH.Power chs.NPUSCHDRSPower])); N = numel(reNames); L = line(ones(N),ones(N),“线宽”8);％生成线条%根据cmap设置颜色集(L, {'颜色'}, mat2cell(提出(min (1 + clevels长度城市规划机构(cmap)),:), (1, N), 3));传奇(重命名{:});