TBS = 144;%传输块大小totNumBlks = 1;%模拟传输块数Ue = struct();初始化UE结构问题。NBULSubcarrierSpacing =“15 khz”；% 3.75kHz, 15kHz问题。NNCellID = 0;%窄带单元识别CHS = struct();% NPUSCH携带数据或控制信息chs。NPUSCHFormat =“数据”；有效载荷类型(数据或控制)用于NPUSCH 'NscRU'的子载波数量取决于NPUSCH%格式和子载波间距'NBULSubcarrierSpacing'，如TS所示表10.1.2.3-1. % 36.211。有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;% NPUSCH的重复次数chs。SlotIdx = 0;包中的起始槽索引符号调制取决于NPUSCH格式和NscRU格式表10.1.3.2-1chs。调制=“正交相移编码”；rvDCI = 0;通过DCI信号的RV偏移百分比(见36.213 16.5.1.2)在dB中指定NPUSCH和DM-RS功率缩放，用于图形可视化chs。NPUSCHPower = 30;chs。npuschdrpower = 34;

chs。SeqGroupHopping =“上”；启用/禁用终端序列组跳转功能chs。SeqGroup = 0;% Delta_SS。上层参数groupAssignmentNPUSCH获取资源单元NULSlots中的时隙个数表10.1.2.3-1 . % TS 36.211如果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必须为" Data "或" Control ", chs.NPUSCHFormat);结束chs。NULSlots = NULSlots;nsslotsperbundle = chs.NRU*chs.NULSlots*chs.NRep;码字包中的槽数TotNSlots = totNumBlks*NSlotsPerBundle;%模拟插槽总数

将随机生成器初始化为默认状态rng (“默认”）;获取每帧的槽网格和槽数emptySlotGrid = lteNBResourceGrid(ue);slotGridSize = size(emptySlotGrid);NSlotsPerFrame = 20/(slotGridSize(1)/12);State = [];% NPUSCH编码器和DM-RS状态，自动重新初始化的功能TRBLK = [];初始化传输块Txgrid = [];%完全网格初始化显示正在生成的槽位数流(生成%d个插槽，对应%d个传输块\n'、TotNSlots totNumBlks);为slotIdx = 0+(0:TotNSlots-1)%计算帧号和帧内槽号问题。NFrame = fix(slotIdx/NSlotsPerFrame);问题。NSlot = mod(slotIdx,NSlotsPerFrame);如果isempty (trblk)如果strcmpi (chs。NPUSCHFormat,“数据”）% UL-SCH编码用于的两个RV值%传输码字。所使用的RV序列已确定%从rvDCI值在DCI和交替信号%在0和2之间，如TS 36.213节16.5.1.2所述定义传输块，该传输块将被编码以创建%码字用于不同RVTRBLK = randi([0 1]，tbs,1);确定编码的传输块大小[~， info] = lteNPUSCHIndices(ue,chs);outblklen = info.G;创建对应于所使用的两个RV值的码字%在第一个和第二个块，这将重复直到所有传输块的百分比chs。RV = 2*mod(rvDCI+0,2);% RV为第一个块cw = 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.3cw = lteNULSCH(trblk);结束blockIdx = 0;要传输的第一个数据块结束初始化网格slotGrid = emptySlotGrid;% NPUSCH编码并映射到插槽网格txsym = lteNPUSCH(ue,chs,cw(:，mod(blockIdx,size(cw,2))+1)，state);将NPUSCH符号映射到插槽的网格中indicesNPUSCH = lteNPUSCHIndices(ue,chs);slotGrid(indicesNPUSCH) = txsym*db2mag(chs.NPUSCHPower);创建DM-RS序列并映射到插槽网格[dmrs,state] = lteNPUSCHDRS(ue,chs,state);indicesDMRS = lteNPUSCHDRSIndices(ue,chs);slotGrid(indicesDMRS) = dmrs*db2mag(chs. npuschdrpower);将此插槽连接到插槽网格txgrid = [txgrid slotGrid];% #好< AGROW >%如果传输了一个完整的块，增加时钟计数器，使%可以选择正确的码字如果状态。EndOfBlk blockIdx = blockIdx + 1;结束% TRBLK错误计数和重新初始化如果状态。EndOfTx重新初始化以启用新的传输块的传输TRBLK = [];结束结束%执行SC-FDMA调制以创建时域波形问题。CyclicPrefixUL =“正常”；NB-IoT正常循环前缀长度[波形，scfdmaInfo] = lteSCFDMAModulate(ue,chs,txgrid);

%创建一个整体资源网格的图像图im = image(abs(txgrid));Cmap = 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(reNames); L = line(ones(N),ones(N),“线宽”8);%生成行根据cmap设置颜色集(L, {“颜色”}, mat2cell(提出(min (1 + clevels长度城市规划机构(cmap)),:), (1, N), 3));传奇(重命名{:});