tbs = 144;传输块大小totNumBlks = 1;模拟传输块的数目问题=结构();初始化UE结构问题。NBULSubcarrierSpacing =“15 khz”;% 3.75 khz, 15 khz问题。NNCellID = 0;窄带单元标识chs =结构();% NPUSCH携带数据或控制信息chs。NPUSCHFormat =“数据”;负载类型(数据或控制)NPUSCH 'NscRU'使用的子载波数量取决于NPUSCH%格式和子载波间距' nbulsubcarrier行距'，如TS所示% 36.211表10.1.2.3-1有1个，3个，6个或12个相邻的子载波% NPUSCHchs.NBULSubcarrierSet= 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;重复次数%的NPUSCHchs。SlotIdx = 0;在一个包中启动插槽索引符号调制依赖于NPUSCH格式和NscRU as%由TS 36.211给出，表10.1.3.2-1chs。调制=“正交相移编码”;rvDCI = 0;% RV偏移通过DCI信号(见36.213 16.5.1.2)%在数据库中指定NPUSCH和DM-RS功率比例，以用于绘图可视化chs。NPUSCHPower = 30;chs。NPUSCHDRSPower = 34;

chs。SeqGroupHopping =“上”;为UE启用/禁用序列组跳转chs。SeqGroup = 0;% Delta_SS。更高层次参数groupAssignmentNPUSCH%获取资源单元空槽中的时间槽数表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“控制”) nulslot = 4;其他的错误('无效的NPUSCH格式(%s)。NPUSCHFormat必须是“数据”或“控制”,chs.NPUSCHFormat);结束chs。NULSlots = NULSlots;NSlotsPerBundle = chs.NRU * chs.NULSlots * chs.NRep;一个码字包中的槽数TotNSlots = totNumBlks * NSlotsPerBundle;%模拟槽的总数

初始化随机生成器为默认状态rng (“默认”);%得到插槽网格和每个帧的插槽数如果strcmpi (ue.NBULSubcarrierSpacing“15 khz”) slotGrid0 = 0 (12,7);NSlotsPerFrame = 20;%槽0…19其他的slotGrid0 = 0(48岁,7);NSlotsPerFrame = 5;%插槽0…4结束状态= [];% NPUSCH编码器和DM-RS状态，自动重新初始化的功能trblk = [];初始化传输块txgrid = [];%全网格初始化显示正在生成的槽的数量流('\ ngenerate %d slot to %d transport block(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] = lteNPUSCHIndices(ue,chs);outblklen = info.G;创建与所使用的两个RV值对应的码字在第一个和第二个块中，将重复此操作，直到所有传送%块chs。房车= 2 *国防部(rvDCI + 0, 2);% RV为第一个块连续波= lteNULSCH (chs、outblklen trblk);% CRC和Turbo编码chs。房车= 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节中定义连续波= lteNULSCH (trblk);结束blockIdx = 0;要传输的第一个块结束%初始化网格slotGrid = slotGrid0;% NPUSCH编码并映射到槽网格txsym = lteNPUSCH (chs问题,连续波(:,国防部(blockIdx、大小(cw, 2)) + 1),状态);将NPUSCH符号映射到槽的网格中chs indicesNPUSCH = lteNPUSCHIndices(问题);slotGrid (indicesNPUSCH) = txsym * db2mag (chs.NPUSCHPower);创建DM-RS序列并映射到槽网格(dmr、州)= lteNPUSCHDRS(问题、chs、状态);chs indicesDMRS = lteNPUSCHDRSIndices(问题);slotGrid (indicesDMRS) = dmr * db2mag (chs.NPUSCHDRSPower);将此插槽连接到插槽网格txgrid = [txgrid slotGrid];% #好< AGROW >如果传输了一个完整的块，则增加时钟计数器，以便可以选择正确的码字如果状态。EndOfBlk blockIdx = blockIdx + 1;结束错误计数和重新初始化如果state.EndOfTx重新初始化以启用新传输块的传输trblk = [];结束结束执行SC-FDMA调制以创建时域波形问题。CyclicPrefixUL =“正常”;NB-IoT的正常循环前缀长度%[波形,scfdmaInfo] = lteSCFDMAModulate (chs退居其次,txgrid);

创建一个整体资源网格的映像figure 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 =圆(db2mag ([chs。NPUSCHPower chs.NPUSCHDRSPower])); N = numel(reNames); L = line(ones(N),ones(N),“线宽”8);%生成线%根据cmap设置颜色集(L, {“颜色”},mat2cell(提出(min (1 + clevels长度城市规划机构(cmap)),:), (1, N), 3));传奇(重命名{:});