%的帧数来模拟NFrames = 4;%设置一系列涵盖高和低信噪比分提单操作%的条件SNRIn = (5 0 5);

sfGap = 8;%时间之间的子帧的初始和重传创科实业= 4;%的时间在不同的HARQ过程之间的子帧如果创科实业> sfGap错误(“不能大于sfGap创科实业”);结束%的数量分配资源块(NPRB)和传输块大小(TBS)%如表A.8.3-1 TS 36.101给出的。注意,TBS必须选择%根据TS 36.213节14.1.1和NPRB TS 36.213% 14.1.1.4C节。如果指定一个无效的TBS / NPRB,传播% ModCoding SCI的价值信息是不正确的,这可能%导致SL-SCH解码失败。一个无效的NPRB也可能的结果%在SCI消息和相应PSSCH传输%不同于这里的NPRB指定分配资源。NPRB = 48;%的数量分配资源块(NPRB)TBS = 3496;%的传输块大小%定义的起始RB资源分配。最初的和%的后续传播可以有不同的配置rbStart = [0 0];%定义问题的配置参数ueConfig =结构(“SidelinkMode”,“V2X”);%释放14 V2X模式ueConfig。NSLRB = 50;% 10 mhz带宽ueConfig。DuplexMode =“FDD”;%双工模式ueConfig。CyclicPrefixSL =“正常”;%循环前缀长度ueConfig。调制=“正交相移编码”;%符号调制(“正交相移编码”、“16 qam”)

频道=结构;%通道配置结构通道。种子= 6;%通道种子通道。NRxAnts = 2;%的接收天线通道。DelayProfile =“爱娃”;%延迟概要通道。DopplerFreq = 500;在赫兹%多普勒频率通道。MIMOCorrelation =“低”;% Multi-antenna相关性通道。NTerms = 16;%振荡器用于衰落模型通道。ModelType =“GMEDS”;%瑞利衰落模型类型通道。InitPhase =“随机”;%随机初始阶段通道。NormalizePathGains =“上”;%正常化延迟概要文件的权力通道。NormalizeTxAnts =“上”;%正常化传输天线%通道采样率取决于用于OFDM FFT大小%调制器。这可以使用函数lteSLSCFDMAInfo获得。ofdmInfo = lteSLSCFDMAInfo (ueConfig);通道。SamplingRate = ofdmInfo.SamplingRate;

%完美的信道估计值的旗帜perfectChanEstimator = false;%定义实际信道估计量的配置结构。请注意,%根据通道延迟资料和多普勒频率,%的操作信噪比,PSSCH调制秩序和分配的数量%的资源块,不同的信道估计参数可能会更好%的性能。cec =结构;%信道估计配置结构cec。PilotAverage =“UserDefined”;%的试点,平均象征cec。FreqWindow = 27个;%频率窗口大小cec。TimeWindow = 1;%的时间窗口大小cec。InterpType =“立方”;% 2 d插值类型

measureBLERForSCI = true;measureTputForSLSCH = true;

sciAssumed = false;

displaySimulationInformation = true;

resourcePool = V2XSidelinkResourcePool;resourcePool.Config。NSLRB = ueConfig.NSLRB;resourcePool.Config。DuplexMode =ueConfig。DuplexMode;如果displaySimulationInformation resourcePool.displayPeriod;drawnow;结束

%初始化变量用于模拟和分析maxThroughputSLSCH = 0(长度(SNRIn), 1);simThroughputSLSCH = 0(长度(SNRIn), 1);simBLERSCI = 0(长度(SNRIn), 1);%得到所需的HARQ进程数量nHARQProcesses =地板(sfGap /创科实业);%计算HARQ ID为每个传输子帧序列harqProcessSequence = 0 (1, NFrames * 10);为h = 1: nHARQProcesses harqProcessSequence(1 +创科实业* (h): sfGap: NFrames * 10) = h;结束%创建部分SCI消息的参数有固定值%的模拟。其他参数,变量值%制定进一步的仿真循环sciMessage =结构(“SCIFormat”,“Format1”);sciMessage。TimeGap= sfGap;%设置根据TBS PSSCH MCS和调制方案[胫,modn] = lteMCS (0:28,“PUSCH”);%根据TS 36.213节14.1.1,imc 28(0),调制%顺序设置为Q ' = min (4, Qm), Qm可以2、4或6。所以改变%的情况Qm ' = 6 (64 qam) Qm”= 4 (16 qam的)modn (strcmpi (modn,64 qam))= {16 qam的};%滤波器有效的传输块大小根据使用的调制方案possibleTBS = lteTBS (NPRB、髂胫);possibleTBS (~ strcmpi (modn ueConfig.Modulation)) = 0;%现在设置根据TBS MCS指数sciMessage。ModCoding= find(possibleTBS==TBS,1) - 1;%图显示创建新的资源分配,如果需要的话如果displaySimulationInformation图;结束%如果SCI提单没有被测量,必须设置,以便sciAssumed%接收机可以假设的科学知识如果(~ measureBLERForSCI) sciAssumed = true;结束rng (“默认”);为snrIdx = 1:元素个数(SNRIn)%设置随机数生成器的种子取决于循环变量%,以确保独立随机流rng (snrIdx);%初始化所有HARQ缓冲区的状态harqProcess =结构(“RVIdx”,1“数据”[],“decState”[]);harqProcesses (1: nHARQProcesses) = harqProcess;RVSeq = [0 2];%房车首次传输和重传%设置变量的主循环lastOffset = 0;%初始化总体框架定时偏移量frameOffset = 0;%初始化帧定时偏移量sharedbitTput = [];%成功收到位/子帧的数量txedTrBlkSizes = [];%每个子帧的传输比特数controlErrors = [];%的SCI块错误为subframeNo = 0: (NFrames * 1) harqID = harqProcessSequence (subframeNo + 1);如果harqID = = 0%如果没有HARQ过程传输电流%子帧,继续下一个继续结束%更新当前HARQ过程与新的传输数据和重置%接收缓冲区,如果是最初的传播。PSSCH是%传播两次的差距sfGap子帧不分%的成功第一个传播如果harqProcesses (harqID)。RVIdx = = 1 harqProcesses (harqID)。data =兰迪([0,1],TBS, 1);harqProcesses (harqID) .decState = [];结束%设置房车ueConfig。房车= RVSeq (harqProcesses (harqID) .RVIdx);%设置子帧数ueConfig。NSubframePSSCH = subframeNo;%显示运行时信息如果displaySimulationInformation流(“子帧:% d。HARQ过程指数:% d。冗余版本:% d \ n ',harqID subframeNo ueConfig.RV)结束%通道当前子帧的时间通道。InitTime = subframeNo / 1000;%创建一个空的网格资源slgrid = lteSLResourceGrid (ueConfig);%更新传输配置与RB开始startingRB = rbStart (harqProcesses (harqID) .RVIdx);resourcePool.Config。startRB_Subchannel_r14 = startingRB;%计算RIV分配子信道数的和%的起始子通道ueConfig。FirstSubchannelIdx = 0;β= 0;如果strcmpi (resourcePool.Config.adjacencyPSCCH_PSSCH_r14“上”)β= 2;结束lsubCH =装天花板((NPRB +β)/ resourcePool.Config.sizeSubchannel_r14);sciMessage。RIV= resourcePool.encodeRIV(lsubCH,ueConfig.FirstSubchannelIdx);%设置传输指数sciMessage。RetransmissionIdx = harqProcesses (harqID)。RVIdx - 1;%创建SCI和消息sciConfig = ueConfig;sciConfig。PSSCHNSubchannels = resourcePool.Config.numSubchannel_r14;[sciMessage, sciBits] = lteSCI (sciConfig sciMessage);%执行SCI编码(连续波,crc) = lteSCIEncode (sciConfig sciBits);% NXID是相关的16位CRC PSCCH SCI格兰特和%作为PSSCH匆忙的身份ueConfig。NXID= crc;% NXID值存储在“expectedNXID”,这样接收器不会%访问“sciMessage”expectedNXID = ueConfig.NXID;%显示传输的子帧和资源块。请注意%,所有HARQ进程的配置都是一样的%情节没有显示的子帧抵消HARQ进程以外的% 1如果(displaySimulationInformation) & & (harqID = = 1)%指定其他参数所需的资源%池显示sciMessage。FirstSubchannelIdx = ueConfig.FirstSubchannelIdx;sciMessage。SLIndex = 0;resourcePool.displayPeriod (sciMessage 10);drawnow;结束%创建PSCCH码字pscchSymbols = ltePSCCH (cw);%得到PSCCH的资源分配[pscchsf, pscchrb] = resourcePool.getPSCCHResources (setfield (sciMessage,“FirstSubchannelIdx”ueConfig.FirstSubchannelIdx));% #好< SFLD >%选择资源传输sciConfig。PRBSet= pscchrb(:,harqProcesses(harqID).RVIdx);%得到PSCCH象征指数pscchIndices = ltePSCCHIndices (sciConfig);%的PSCCH符号映射到网格slgrid (pscchIndices) = pscchSymbols;%定义一组随机循环移位从{0、3、6、9}% PSCCH传输sciConfig。CyclicShift =兰迪([0 3])* 3;%创建DM-RS符号和指数pscchdrssymbols = ltePSCCHDRS (sciConfig);pscchdrsindices = ltePSCCHDRSIndices (sciConfig);%将PSCCH DRS符号映射到网格slgrid (pscchdrsindices) = pscchdrssymbols;%得到两个传输的资源分配[psschsf, psschrb] = resourcePool.getPSSCHResources (setfield (sciMessage,“FirstSubchannelIdx”ueConfig.FirstSubchannelIdx));% #好< SFLD >%现在选择的资源传播ueConfig。PRBSet= psschrb(:,harqProcesses(harqID).RVIdx);%计算PSSCH资源指数[psschIndices, psschinfo] = ltePSSCHIndices (ueConfig);% SL-SCH和PSSCH连续波= lteSLSCH (ueConfig、psschinfo.G harqProcesses (harqID) . data);psschSymbols = ltePSSCH (ueConfig, cw);slgrid (psschIndices) = psschSymbols;% PSSCH DRS和DRS指数psschdrssymbols = ltePSSCHDRS (ueConfig);psschdrsindices = ltePSSCHDRSIndices (ueConfig);slgrid (psschdrsindices) = psschdrssymbols;% SC-FDMA调制不会使用窗口由于%分段波形结构的性质窗口= 0;%执行sidelink SC-FDMA调制[txWaveform, scfdmaInfo] = lteSLSCFDMAModulate (ueConfig slgrid,窗口);% 25样品添加填充。这是延迟的范围%预期从信道建模(的组合%实现延迟和延迟扩展频道)ntxants =大小(txWaveform, 2);txWaveform = [txWaveform;ntxants 0(25日)];% #好< AGROW >%通过通道模型传递数据rxNoiselessWaveform = lteFadingChannel(通道、txWaveform);%计算噪声增益信噪比= 10 ^ (SNRIn (snrIdx)) / 20);%正常化噪声功率的采样率,%是传输线大小的函数用于OFDM%调制,天线的数量N0 = 1 / (√(2.0 * ntxants *双(ofdmInfo.Nfft)) *信噪比);%建立情况噪音= N0 *复杂(randn(大小(rxNoiselessWaveform)),…randn(大小(rxNoiselessWaveform)));%添加情况下接收到的时域波形rxWaveform = rxNoiselessWaveform +噪声;% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -%接收% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -%定义接收机顶级配置rxConfig = [];rxConfig。SidelinkMode = ueConfig.SidelinkMode;rxConfig。CyclicPrefixSL = ueConfig.CyclicPrefixSL;rxConfig。FirstSubchannelIdx = ueConfig.FirstSubchannelIdx;rxConfig。NSLRB = resourcePool.Config.NSLRB;rxConfig。DuplexMode = resourcePool.Config.DuplexMode; rxConfig.TDDConfig = resourcePool.Config.TDDConfig; rxConfig.PSSCHNSubchannels = resourcePool.Config.numSubchannel_r14; rxConfig.PSSCHSubchannelSize = resourcePool.Config.sizeSubchannel_r14; rxConfig.PSSCHAdjacency = resourcePool.Config.adjacencyPSCCH_PSSCH_r14; rxConfig.PSSCHSubchannelPRBStart = resourcePool.Config.startRB_Subchannel_r14; rxConfig.NSubframePSSCH = subframeNo;%设置子帧数%如果SCI提单测量配置:sciDecoded = false;如果measureBLERForSCI%控制信道的盲解码%我们知道当前的子帧包含控制和数据,%,但不知道位置。检查中的SCI消息%的资源池循环直到发现或转变%搜索整个资源池pscchPool = resourcePool.PSCCHResourceBlockPool;%为可能的复审委员会对重建资源pscchPool =重塑(pscchPool 2元素个数(pscchPool) / 2);p = 1;而p < = (pscchPool, 2) & & ~ sciDecoded大小% PSCCH DRS V2X传播在两个连续的苏格兰皇家银行rxConfig。PRBSet= pscchPool(:,p);为c = [0 3 6 9] rxConfig。CyclicShift = c;(frameOffset, corr) = lteSLFrameOffsetPSCCH (rxConfig rxWaveform);%执行子帧同步如果(frameOffset > 25) | | (frameOffset < 0) frameOffset = lastOffset;结束lastOffset = frameOffset;%对接收的数据获得执行SC-FDMA解调%的资源网格rxSubframe = lteSLSCFDMADemodulate (rxConfig rxWaveform (1 + frameOffset:最终,));%进行信道估计如果(perfectChanEstimator)[命令,巢]= perfectChannelEstimate (rxConfig、通道、噪音frameOffset);其他的[命令,巢]= lteSLChannelEstimatePSCCH (cec, rxConfig rxSubframe);结束%得到当前PRBSet PSCCH候选人,%提取相应的符号和接收%信道估计和均衡执行pscchCandidate = ltePSCCHIndices (rxConfig);[pscchRx, pscchHest] = lteExtractResources (pscchCandidate rxSubframe,命令);pscchSymbols = lteEqualizeMMSE (pscchRx pscchHest,巢);%解调PSCCHcodedsciBits = ltePSCCHDecode (pscchSymbols);%得到释放的有效载荷大小14 V2X SCI消息scilengthinfo = lteSCIInfo (rxConfig);%解码SCI和恢复消息[rxinfo, rxSCIErr rxcrc] = lteSCIDecode (scilengthinfo.Format1 codedsciBits);如果~ rxSCIErr & & (rxcrc = = expectedNXID) sciMessageRx = lteSCI (rxConfig rxinfo);rxConfig。NXID= rxcrc; sciDecoded = true;如果(displaySimulationInformation)流(“SCI解码\ n”);结束打破;结束结束%增加池索引p = p + 1;结束如果~ sciDecoded如果(displaySimulationInformation)流(“SCI解码失败\ n”);结束结束%计算提单所需存储值controlErrors = [controlErrors ~ sciDecoded];% #好< AGROW >结束%如果SCI解码失败和SL-SCH SCI假定%解码,解码SCI消息等于传播% SCI消息。还设置NXID期望值如果(~ sciDecoded) & & (sciAssumed) sciMessageRx = sciMessage;rxConfig。NXID= expectedNXID;结束%如果SL-SCH吞吐量测量配置和科学%成功解码或假设如果SCI解码成功,%执行PSSCH接待和SL-SCH解码如果measureTputForSLSCH如果(sciDecoded | | sciAssumed)%得到PSSCH解码SCI的资源分配%的消息和池配置rxConfig。PRBSet= lteSCIResourceAllocation(rxConfig,sciMessageRx);%执行子帧同步frameOffset = lteSLFrameOffsetPSSCH (rxConfig rxWaveform);如果(frameOffset > 25) | | (frameOffset < 0) frameOffset = lastOffset;结束lastOffset = frameOffset;%对接收的数据获得执行SC-FDMA解调%的资源网格rxSubframe = lteSLSCFDMADemodulate (rxConfig rxWaveform (1 + frameOffset:最终,));%进行信道估计,提取PSSCH收到%符号和相应的信道估计%进行均衡如果(perfectChanEstimator)[命令,巢]= perfectChannelEstimate (rxConfig、通道、噪音frameOffset);其他的[命令,巢]= lteSLChannelEstimatePSSCH (cec, rxConfig rxSubframe);结束[psschRx, psschHest] = lteExtractResources (psschIndices rxSubframe,命令);psschSymbols = lteEqualizeMMSE (psschRx psschHest,巢);%解调PSSCH[~,rxConfig。调制]= lteMCS (sciMessageRx.ModCoding“PUSCH”);如果(strcmpi (rxConfig.Modulation64 qamrxConfig))。调制=16 qam的;结束codedSlschBits = ltePSSCHDecode (rxConfig psschSymbols);%解码SL-SCH,包括软结合到%接收缓冲区并检查CRCrxConfig。NTurboDecIts = 5;%涡轮迭代译码器rxConfig。房车= RVSeq (harqProcesses (harqID) .RVIdx);[~,slschCRC harqProcesses (harqID) .decState] =…codedSlschBits lteSLSCHDecode (rxConfig TBS,…harqProcesses (harqID) .decState);slschDecoded = ~ slschCRC;其他的%如果SCI解码失败,SL-SCH解码也失败了slschDecoded = false;结束%存储值计算所需的吞吐量sharedbitTput = [sharedbitTput TBS * slschDecoded];% #好< AGROW >txedTrBlkSizes = [txedTrBlkSizes TBS];% #好< AGROW >如果(displaySimulationInformation)如果slschDecoded流(“SL-SCH解码\ n”);其他的流(“SL-SCH解码失败\ n”);结束结束结束%更新RV序列索引下传播harqProcesses (harqID)。RVIdx =国防部(harqProcesses harqID .RVIdx,大小(RVSeq, 2)) + 1;结束如果measureBLERForSCI%计算SCI提单simBLERSCI (snrIdx) = 100 *意味着(controlErrors, 2);%动态在命令窗口中显示结果流(' \ nSNR = %。2 f dB。SCI提单(% %)% d帧(s) = %。4 f % % \ n”,…SNRIn (snrIdx)、NFrames simBLERSCI (snrIdx));结束如果measureTputForSLSCH%计算和模拟的最大吞吐量maxThroughputSLSCH (snrIdx) =总和(txedTrBlkSizes);%最大可能的吞吐量simThroughputSLSCH (snrIdx) =总和(sharedbitTput);%模拟吞吐量%动态在命令窗口中显示结果流(' \ nSNR = %。2 f dB。SL-SCH吞吐量(s) = % % d帧。4 f Mbps \ n”,…SNRIn (snrIdx) NFrames 1 e-6 * simThroughputSLSCH (snrIdx) / (NFrames * 10 e - 3));流('信噪比= %。2 f dB。SL-SCH吞吐量(% %)% d帧(s) = %。4 f % % \ n”,…SNRIn (snrIdx)、NFrames simThroughputSLSCH (snrIdx) * 100 / maxThroughputSLSCH (snrIdx));结束结束

如果measureBLERForSCI%情节PSCCH / SCI提单plotResults (“科学”,“提单”、通道、SNRIn simBLERSCI);结束如果measureTputForSLSCH%绘制PSSCH / SL-SCH吞吐量plotResults (“SL-SCH”,的吞吐量、通道、SNRIn simThroughputSLSCH * 100. / maxThroughputSLSCH);结束

%计算完美的信道估计和噪声估计函数[命令,巢]= perfectChannelEstimate (rxConfig、通道、噪音、frameOffset)%使用完美上行信道估计量来计算%都使用SC-FDMA sidelink估计rxConfig。NTxAnts = 1;rxConfig。NULRB = rxConfig.NSLRB;命令= lteULPerfectChannelEstimate (rxConfig、通道frameOffset);noiseGrid = lteSLSCFDMADemodulate (rxConfig、噪音(1 + frameOffset:最终,));巢= var (noiseGrid (:));结束%绘制SCI提单或SL-SCH吞吐量函数plotResults (chName simOpt、通道、xvalues yvalues) fadingDescription = sprintf (' % s % s '、channel.DelayProfile num2str (channel.DopplerFreq));图;情节(xvalues yvalues,* - - - - - -。);v =轴;轴([v(1)(2) -10 110])标题(sprintf (“V2X Sidelink % s % s (% %) % s衰落和AWGN”,simOpt chName fadingDescription));包含(“信噪比(dB)”);ylabel ([simOpt“(%)”]);网格在;结束