问题=结构();％初始化结构问题。NULRB = 50;%的带宽问题。DuplexMode =“FDD”；％双工模式问题。TDDConfig = 1;% UL/DL配置，如果TDD双工模式问题。CyclicPrefixUL =“正常”；％循环前缀长度问题。NCellID = 1;％细胞标识UE.rnti = 1;% RNTI价值问题。NFrame = 0;%帧数问题。NSubframe = 0;%子帧数问题。NTxAnts = 1;%发送天线数问题。缩短= 1;%最后一个符号可用性(SRS分配)%设置跳跃的具体参数问题。HoppingOffset = 1;一个窄带和下一个窄带之间％窄带偏移%一个PUSCH跳到，表示为上行窄带数UE.Nchulnb = 2;%所覆盖的连续绝对子帧数％PUCCH或PUSCH保持在同​​一个pusch =结构();pusch。CEMode =“一个”；% CE模式A或CE模式Bpusch.hopping.= true;％启用/禁用跳频pusch。NRepPUSCH = 8;% PUSCH总重复次数pusch.modulation =.'QPSK'；%符号调制pusch。房车= 0;% UL-SCH处理的冗余版本pusch.nlayers = 1;％的层数pusch.trblksizes = 100;％运输块大小

%为所有情况指定窄带内基于1的RBs相对索引%除了5MHz Cat-M2 CE模式A，如果5MHz Cat-M2 CE模式A，这些是绝对PRBs用于传播的百分比pusch。InitPRBSet = (2:3) ';%用于传输的窄带(非跳频，非5mhz)pusch。InitNarrowbandIndex = 1;%指定PUSCH的功率缩放，以dB为单位，PUSCH DM-RSpusch.puschpower = 30;pusch.puschdmrspower = 100;%如果分配跨越多个窄带，关闭跳变如果numel(pusch. initprbset) > 6 pusch. initprbset跳= false;结束

％识别帧中的所有上行链路子帧信息= arrayfun (@ (x) lteDuplexingInfo (setfield(问题,“NSubframe”, x)), 0:9);ulsfs = arrayfun (@ (x) strcmpi (x。SubframeType,“上行”），信息）;％在这个例子中，我们假设第一个绝对子帧% PUSCH传输的是第一个可用的上行子帧pusch。InitNSubframe =找到(ulsfs, 1) 1;计算分配pusch。pusch PRBSet = getPUSCHAllocation(问题);ueTemp =问题;％为所有符号创建编码传输块如果strcmpi (pusch。CEMode,“B”) & &问题。缩短ueTemp。缩短= 0;结束[~,信息]= ltePUSCHIndices (ueTemp pusch);％定义UL-SCH消息比特TRDATA = ON（PUSCH.TRBLKSIZE（1），1）;％创建编码的UL-SCH位pusch。BetaCQI = 2.0;pusch。BetaRI = 2.0;pusch。BetaACK = 2.0;codedTrBlock = lteULSCH (ueTemp pusch、info.G trData);

%置乱块中的子帧数Nacc = 1;如果strcmpi（ue.duplexmode，“FDD”）&& strcmpi（pusch.code，“B”) Nacc = 4;elseifstrcmpi（ue.duplexmode，'TDD'）&& strcmpi（pusch.code，“B”) Nacc = 5;结束要模拟的％BL / CE子帧（所有上行链路子帧都是BL / CE％子帧）和PUSCH在没有任何子帧间隙的情况下发送totmtcsubframes = pusch.nreppusch;要模拟的绝对子帧总数startSubframe = ue.NFrame * 10 + ue.NSubframe;%初始绝对子帧数Lastabssf = getlastabssf（ulsfs，pusch.initnsubframe，totmtcsubframes）;totsubframes = lastabssf-startsubframe + 1;％为整个传输创建资源网格。脓液和% DM-RS符号将被映射到这个数组中子帧= lteulresourcegridsize（UE）;sfgrid = zeros（[子帧大小（1）子帧（2）* totsubframes子帧（3：结束）]）;mpuschsym = [];%初始化PUSCH符号为sf = startSubframe + (0: tosubframes -1)%设置当前的绝对子帧和帧数ue.nsubframe = mod（sf，10）;UE.NFRAME =楼层（（SF）/ 10）;%如果这不是上行子帧，跳过处理duplexInfo = lteDuplexingInfo(问题);如果~ strcmpi (duplexInfo。SubframeType,“上行”）继续结束％计算当前子帧中使用的PRBSETprbset = getPuschallocation（UE，PUSCH）;%计算当前子帧的PDSCH索引。提单/ CE问题%在CE模式B中，a中最后一个SC-FDMA符号中的资源元素配置有单元特定SRS的%子帧将被计算在% PUSCH映射，但不用于PUSCH的传输pusch。PRBSet = PRBSet;pusch mpuschIndices = ltePUSCHIndices(问题);%创建一个空的子帧网格子帧= lteULResourceGrid(问题);％从码字中编码PUSCH符号在重复的情况下，相同的符号在每一个中重复%一个NRepPUSCH子帧块。根据需要进行跳频如果~mod(sf,Nacc) || isempty(mpuschSym) ueTemp = ue;如果strcmpi (pusch。CEMode,“B”) & &问题。缩短ueTemp。缩短= 0;%为完整的子帧创建符号结束mpuschSym = ltePUSCH (ueTemp pusch codedTrBlock) * db2mag (pusch.PUSCHPower);结束%将SRS戳穿的PUSCH符号映射到子帧网格子帧(mpuschIndices) = mpuschSym(1:元素个数(mpuschIndices));%创建并映射DMRS符号。问题。跳=“关闭”；% DRS跳跃问题。SeqGroup = 0;% PUSCH序列组问题。CyclicShift = 0;%用于n1DMRSLTE-M UE的％，应当假设循环移位字段何时%确定n2DMRS从表5.5.2.1.1 TS 36.211pusch。DynCyclicShift = 0;% n2DMRS循环移位'000'pusch.torthcover =“关闭”；%无正交覆盖序列pusch mpuschDrs = ltePUSCHDRS(问题)* db2mag (pusch.PUSCHDMRSPower);pusch mpuschDrsIndices = ltePUSCHDRSIndices(问题);子帧(mpuschDrsIndices) = mpuschDrs;%现在将当前的子帧分配到整个网格中sfgrid (:, (1: subframeSize(2)) +科幻* subframeSize(2):) =子帧;结束

WaveForm = LTESCFDMAMODULUTS（UE，SFGRID）;

创建一个整体资源网格的映像。自从PUSCH开始％变换预编码，我们需要为全部分配单个功率电平％符号以在绘图中可视化plotgrid = abs（sfgrid（：，：，1））;%获取DM-RS位置drspos = (plotgrid = = db2mag (pusch.PUSCHDMRSPower));plotgrid (drspos) = 0;％现在将所有PUSCH符号设置为一个功率电平进行绘图plotgrid (plotgrid ~ = 0) = db2mag (pusch.PUSCHPower);%现在写回DRS和绘图plotgrid (drspos) = db2mag (pusch.PUSCHDMRSPower);图im = image(plotgrid);提出= parula (64);colormap (im.Parent提出);轴xy；标题(sprintf ('LTE-M CEMode%s Uplink RE Grid (NRepPUSCH = %d)'，pusch.ceode，pusch.nreppusch））xlabel（OFDM符号的）ylabel（'子载波'）%创建图例框以指示与REs关联的通道/信号类型重命名= {'pusch'；'pusch drs'};clevels =圆(db2mag ([pusch。PUSCHPower pusch.PUSCHDMRSPower])); N = numel(reNames); L = line(ones(N),ones(N),“线宽”8);%生成线%根据cmap设置颜色集(L, {'颜色'}，mat2cell（cmap（min（min（1 + clevels，length（cmap）），:)，一个（1，n），3））;%根据cmap设置颜色传奇(重命名{:});

％计算LTE载波带宽的宽带，窄带和PRBSETS功能[prbsets，nnb，nwb] = calcnarrowbandprbsets（nulrb）窄带和宽带(见36.211章节5.2.4)NULNB =地板(NULRB / 6);nNB = 0: (NULNB-1);%窄带如果nulnb> = 4 nulwb =地板（nulnb / 4）;其他的NULWB = 1;结束nwb = 0 :( nulwb-1）;%宽带% prb窄带2 = 0:5;i0 = floor(NULRB/2) - 6*(NULNB/2);prbsets = 0(6元素个数(nNB));为nb = 1：numel（nnb）如果国防部(NULRB, 2) & & nNB (nb) > = (NULNB / 2) prbsets (:, nb) = 6 * (nNB (nb)) + ii0 + 2 + 1;其他的prbsets (:, nb) = 6 * (nNB (nb)) + ii0 + 2;结束结束结束%计算子帧中为PUSCH分配的资源块功能chs prbset = getPUSCHAllocation(问题)%如果5MHz模式(最多可以使用24个prb)，分配是相同的％是initprbset.如果numel(chs.InitPRBSet) > 6 prbset = chs.InitPRBSet;返回；结束%获取窄带和相应的资源[prbsets, nNB] = calcNarrowbandPRBSets (ue.NULRB);如果max(chs. init窄bandindex) > max(nNB) error('指定的窄带无效。从0 ...％d'的带宽中只有％d窄带，NNB + 1，NNB）;结束%如果跳频被禁用，分配是相同的%每个子帧如果~ chs。跳跃prbset = prbset (chs.InitPRBSet,chs.InitNarrowbandIndex+1);返回结束％跳窄带计算根据TS 36.211第5.3节j₀地板= ((chs.InitNSubframe) / ue.NChULNB);%计算该子帧的窄带问题。NSubframe = ue.NFrame * 10 + ue.NSubframe;%获取绝对子帧数如果mod(floor(ue.NSubframe/ue.NChULNB-j0)，2) == 0 nnBi = chs. init窄带索引;其他的nnBi =国防部(chs.InitNarrowbandIndex + ue.HoppingOffset元素个数(nNB));结束％计算此子帧的PRBSET，它们位于同一RB上%在窄带内[rbstartIndex,nbstartIndex] = find(prbsets == chs.InitPRBSet(1));[rbendIndex,nbendIndex] = find(prbsets == chs.InitPRBSet(end));如果（isempty（rbstartindex）|| isempty（rbendindex））||（nbstartIndex〜= nbendIndex）错误（'指定的PRBSet无效，必须是单一窄带内的资源')；结束prbset = prbsets (rbstartIndex: rbstartIndex +元素个数(chs.InitPRBSet) 1, nnBi + 1);结束%获取最后一次传输使用的绝对子帧数信道百分比功能lastabssf = getlastabsSF(ulsfs,InitNSubframe,totmtcSubframes) numulsfsinFrame = sum(ulsfs);%每帧中活动的sfs数量ulsfsinFrame =找到(ulsfs);帧中的% UL子帧(基于1)找到第一个绝对子帧和帧initabssf = mod（ini​​tnsubframe，10）;initabsf =楼层（initnsubframe / 10）;startidxwithinframe = initabssf + 1;基于% 1的索引到UL sf如果~ ulsfs (startIdxwithinFrame)错误('第一个上行子帧的绝对子帧数无效'，．..'用于指定的PUSCH (%d)。这不是上行子帧'】,InitNSubframe)结束sfslastframe = mod（（find（ulsfsinframe == startidxwithInframe）-1）+ totmtcsubframes，numulesfsinframe）;%子帧到最后一帧的tx如果sfslastFrame％找到与上次子帧相对应的子帧号sfsnumlastFrame =找到(ulsfs sfslastFrame) 1;sfsnumlastFrame = sfsnumlastFrame(结束);其他的%不需要部分帧sfsnumlastFrame = 0;结束lastabssf = (initabsf + floor(((find(ulsfsinFrame==startIdxwithinFrame)-1)+totmtcSubframes)/numulsfsinFrame)) * 10 + sfsnumlastFrame;结束