NFrames = 2;%的帧数SNRIn = (10.3 12.3 14.3);在dB %信噪比范围

txMode =“TM4”;% TM1, TM2、TM3 TM4, TM6

simulationParameters = [];%明显simulationParameterssimulationParameters。NDLRB = 50;simulationParameters.PDSCH。TargetCodeRate = 0.5;simulationParameters.PDSCH。PRBSet = (0:49) ';开关txMode%单天线(Port0)模式(TM1)情况下“TM1”流(“\ nTM1——单天线(端口0)\ n ');simulationParameters.PDSCH。TxScheme =“Port0”;simulationParameters.PDSCH。DCIFormat =“Format1”;simulationParameters。CellRefP = 1;simulationParameters.PDSCH。调制= {16 qam的};%发射分集模式(TM2)情况下“TM2”流(“\ nTM2——发射分集\ n”);simulationParameters.PDSCH。TxScheme =“TxDiversity”;simulationParameters.PDSCH。DCIFormat =“Format1”;simulationParameters。CellRefP = 2;simulationParameters.PDSCH。调制= {16 qam的};simulationParameters.PDSCH。NLayers = 2;% CDD模式(TM3)情况下“TM3”流(“\ nTM3 CDD \ n”);simulationParameters.PDSCH。TxScheme =CDD的;simulationParameters.PDSCH。DCIFormat =“Format2A”;simulationParameters。CellRefP = 2;simulationParameters.PDSCH。调制= {16 qam的,16 qam的};simulationParameters.PDSCH。NLayers = 2;%的空间复用模式(TM4)情况下“TM4”流(“\ nTM4——基于码书的空间复用\ n”);simulationParameters。CellRefP = 2;simulationParameters.PDSCH。调制= {16 qam的,16 qam的};simulationParameters.PDSCH。DCIFormat =“Format2”;simulationParameters.PDSCH。TxScheme =“SpatialMux”;simulationParameters.PDSCH。NLayers = 2;%没有速率限制simulationParameters.PDSCH。CodebookSubset =”;%单层空间复用模式(TM6)情况下“TM6”流(…“\ nTM6——基于码书的空间复用与单层\ n”);simulationParameters。CellRefP = 4;simulationParameters.PDSCH。调制= {“正交相移编码”};simulationParameters.PDSCH。DCIFormat =“Format2”;simulationParameters.PDSCH。TxScheme =“SpatialMux”;simulationParameters.PDSCH。NLayers = 1;%没有速率限制simulationParameters.PDSCH。CodebookSubset =”;否则错误(的传输方式应该是TM1, TM2、TM3, TM4或TM6”。)结束%设置其他simulationParameters领域应用经颅磁刺激simulationParameters。TotSubframes = 1;%生成一副框架simulationParameters.PDSCH。CSI =“上”;%软比特被CSI加权

enb = lteRMCDL (simulationParameters);

rvSequence = enb.PDSCH.RVSeq;trBlkSizes = enb.PDSCH.TrBlkSizes;

作战飞机=长度(string (enb.PDSCH.Modulation));

pmiDelay = 8;

hDisplayENBParameterSummary (enb txMode);

通道。种子= 6;%通道种子通道。NRxAnts = 2;% 2接收天线通道。DelayProfile =“环保署”;%延迟概要通道。DopplerFreq = 5;%多普勒频率通道。MIMOCorrelation =“低”;% Multi-antenna相关性通道。NTerms = 16;%振荡器用于衰落模型通道。ModelType =“GMEDS”;%瑞利衰落模型类型通道。InitPhase =“随机”;%随机初始阶段通道。NormalizePathGains =“上”;%正常化延迟概要文件的权力通道。NormalizeTxAnts =“上”;%正常化传输天线%通道模型的采样率设置使用返回的值从< docid: lte_ref # % bt1bgp9 lteOFDMInfo >。ofdmInfo = lteOFDMInfo (enb);通道。SamplingRate = ofdmInfo.SamplingRate;

%完美的信道估计值的旗帜perfectChanEstimator = false;

%配置信道估计cec。PilotAverage =“UserDefined”;%的试点,平均象征cec。FreqWindow = 41;在REs %频率窗口大小cec。TimeWindow = 27个;在REs %时间窗口大小cec。InterpType =“立方”;% 2 d插值类型cec。InterpWindow =“中心”;%插值窗口类型cec。InterpWinSize = 1;%插值窗口大小

displaySimulationInformation = true;

%的数量传输天线P是获得资源网格%的维度。“暗”M-by-N-by-P M是副载波的数量,N%是符号的数量和P是发射天线的数量。= lteDLResourceGridSize dim (enb);P = dim (3);%初始化变量用于模拟和分析%数组来存储所有信噪比点的最大吞吐量maxThroughput = 0(长度(SNRIn), 1);%数组来存储所有信噪比点的模拟吞吐量simThroughput = 0(长度(SNRIn), 1);%得到HARQ的HARQ ID序列处理。这是一个索引的列表% HARQ进程调度。[~,~,enbOut] = lteRMCDLTool (enb, []);harqProcessSequence = enbOut.PDSCH.HARQProcessSequence;%临时变量“enb_init”,和“channel_init”%的harqProcessSequence_init用于优化并行处理(只有%如果运行示例与并行计算工具箱)enb_init = enb;channel_init =通道;harqProcessSequence_init = harqProcessSequence;legendString = [的吞吐量:char (enb.PDSCH.TxScheme)];allRvSeqPtrHistory =细胞(1,元素个数(SNRIn));nFFT = ofdmInfo.Nfft;为snrIdx = 1:元素个数(SNRIn)% parfor snrIdx = 1:元素个数(SNRIn)%,使并行计算的使用增加速度注释掉%以上”的声明和取消下面的“parfor”声明。%这需要并行计算工具箱。如果这是没有安装%’parfor将默认为正常”的声明。如果“parfor”建议使用%变量“displaySimulationInformation”%以上是设置为false,否则仿真信息显示%每个信噪比点会重叠。%设置随机数生成器的种子不同循环变量%,以确保独立随机流rng (snrIdx“combRecursive”);SNRdB = SNRIn (snrIdx);流(在% g \ nSimulating dB信噪比% d (s) \ n帧的、SNRdB NFrames);%初始化变量用于模拟和分析偏移量= 0;%初始化帧偏移值抵消= 0;%初始化帧偏移值广播帧blkCRC = [];%块CRC认为子帧bitTput = [];%成功收到位/子帧的数量txedTrBlkSizes = [];%每个子帧的传输比特数enb = enb_init;%初始化RMC配置频道= channel_init;%初始化通道配置harqProcessSequence = harqProcessSequence_init;%初始化HARQ过程序列pmiIdx = 0;%的PMI指数延迟队列%变量harqPtrTable商店的历史价值%指针RV序列值的HARQ进程。%与nan预先分配一些子帧没有数据rvSeqPtrHistory =南(ncw, NFrames * 10);%初始化所有HARQ进程的状态harqProcesses = hNewHARQProcess (enb);%使用随机采购经理人指数第一pmiDelay子帧,直到反馈%可以从问题;注意,只适用于PMI反馈%的空间复用经颅磁刺激(TM4和TM6),但这里的代码是必需的%完成初始化的变量的信噪比循环使用%的并行计算工具。pmidims = ltePMIInfo (enb enb.PDSCH);txPMIs =兰迪([0 pmidims。MaxPMI], pmidims。NSubbands pmiDelay);为subframeNo = 0: (NFrames * 10 - 1)%更新子帧数enb。NSubframe = subframeNo;%从HARQ过程得到HARQ进程ID为子帧序列harqID = harqProcessSequence (mod (subframeNo长度(harqProcessSequence)) + 1);%如果有一块运输计划在当前的子帧%(零harqID),执行和传播%接待。否则继续接下来的子帧如果harqID = = 0继续;结束%更新当前HARQ进程harqProcesses (harqID) = hHARQScheduling (…harqProcesses (harqID)、subframeNo rvSequence);%提取当前子帧传输块大小(s)trBlk = trBlkSizes(:,国防部(subframeNo, 10) + 1)”;%显示运行时信息如果displaySimulationInformation disp (' ');disp ([的子帧:num2str (subframeNo)…”。HARQ进程ID:“num2str (harqID)]);结束%更新房车序列指针表rvSeqPtrHistory (:, subframeNo + 1) =…harqProcesses harqID .txConfig.RVIdx。';%更新PDSCH传输配置与HARQ进程状态enb。PDSCH = harqProcesses (harqID) .txConfig;data = harqProcesses (harqID) . data;%将PMI设置为适当的值在延迟队列如果strcmpi (enb.PDSCH.TxScheme“SpatialMux”)pmiIdx =国防部(subframeNo pmiDelay);%的PMI指数延迟队列enb.PDSCH。PMISet = txPMIs (:, pmiIdx + 1);%设置采购经理人指数结束%产生传输波形txWaveform = lteRMCDLTool (enb、数据);% 25样品添加填充。这是延迟的范围%预期从信道建模(的组合%实现延迟和延迟扩展频道)txWaveform = [txWaveform;0(25日P)];% #好< AGROW >%得到HARQ ID序列的enbOut HARQ处理harqProcessSequence = enbOut.PDSCH.HARQProcessSequence;%初始化通道为每个子帧时间通道。InitTime = subframeNo / 1000;%通过通道模型传递数据rxWaveform = lteFadingChannel(通道、txWaveform);%计算噪声增益包括下行功率补偿%分配信噪比= 10 ^ ((SNRdB-enb.PDSCH.Rho) / 20);%正常化噪声功率的采样率,这是%传输线大小的函数用于OFDM调制,以及%的天线N0 = 1 / (√(2.0 * enb.CellRefP *双(nFFT)) *信噪比);%建立加性高斯白噪声噪音= N0 *复杂(randn(大小(rxWaveform)),…randn(大小(rxWaveform)));%添加情况下接收到的时域波形rxWaveform = rxWaveform +噪声;%每帧,在子帧0,计算一个新的同步%抵消如果(国防部(subframeNo 10) = = 0)抵消= lteDLFrameOffset (enb rxWaveform);如果(抵消> 25)抵消=补偿(结束);结束偏移量=(补偿抵消);% #好吧结束%同步接收到的波形rxWaveform = rxWaveform(1 +抵消:最终,);%完成OFDM解调接收的数据重新创建%资源网格rxSubframe = lteOFDMDemodulate (enb rxWaveform);%信道估计如果(perfectChanEstimator) estChannelGrid = lteDLPerfectChannelEstimate (enb,通道,抵消);% #好吧noiseGrid = lteOFDMDemodulate (enb、噪音(1 +抵消:最终,));噪音= var (noiseGrid (:));其他的[estChannelGrid,噪音]= lteDLChannelEstimate (…enb, enb。PDSCH、cec rxSubframe);结束%得到PDSCH指数pdschIndices = ltePDSCHIndices (enb, enb。PDSCH enb.PDSCH.PRBSet);%得到PDSCH资源元素从接收到的子帧。规模的%接收子帧PDSCH功率因数ρ。PDSCH是%比例的金额,而细胞符号用于参考%信道估计(用于PDSCH解码阶段)。[pdschRx, pdschHest] = lteExtractResources (pdschIndices,…rxSubframe * (10 ^ (-enb.PDSCH.Rho / 20)), estChannelGrid);%解码PDSCHdlschBits = ltePDSCHDecode (…enb, enb。PDSCH、pdschRx pdschHest、噪音);%解码DL-SCH[decbits harqProcesses (harqID) .blkerr harqProcesses (harqID) .decState] =…lteDLSCHDecode (enb enb。PDSCH、trBlk dlschBits,…harqProcesses (harqID) .decState);%显示块的错误如果displaySimulationInformation如果任何(harqProcesses (harqID) .blkerr) disp ([“块错误。房车指数:num2str (harqProcesses (harqID) .txConfig.RVIdx)…“CRC:”num2str (harqProcesses (harqID) .blkerr)])其他的disp ([“没有错误。房车指数:num2str (harqProcesses (harqID) .txConfig.RVIdx)…“CRC:”num2str (harqProcesses (harqID) .blkerr)])结束结束%存储值计算的吞吐量%的子帧数据如果任何(trBlk) blkCRC = [blkCRC harqProcesses (harqID) .blkerr];% #好< AGROW >bitTput = [bitTput trBlk。* (1 -…harqProcesses (harqID) .blkerr)];% #好< AGROW >txedTrBlkSizes = [txedTrBlkSizes trBlk];% #好< AGROW >结束% eNodeB提供PMI的反馈如果strcmpi (enb.PDSCH.TxScheme“SpatialMux”)采购经理人指数= ltePMISelect (enb, enb。PDSCH estChannelGrid,噪音);txPMIs (:, pmiIdx + 1) =采购经理人指数;结束结束%计算和模拟吞吐量最大maxThroughput (snrIdx) =总和(txedTrBlkSizes);%最大可能的吞吐量simThroughput (snrIdx) = (bitTput, 2)总和;%模拟吞吐量%动态在命令窗口中显示结果流([[' \ nThroughput (Mbps)的num2str (NFrames)框架(年代)”),…“= % .4f \ n”),1 e-6 * simThroughput (snrIdx) / (NFrames * 10 e - 3));流([的吞吐量(% %)num2str (NFrames)的帧(s) = % .4f \ n '),…simThroughput (snrIdx) * 100 / maxThroughput (snrIdx));allRvSeqPtrHistory {snrIdx} = rvSeqPtrHistory;结束为所有HARQ进程%情节RV序列hPlotRVSequence (SNRIn allRvSeqPtrHistory NFrames);

%绘制吞吐量图绘制(SNRIn simThroughput * 100. / maxThroughput,* - - - - - -。);包含(“信噪比(dB)”);ylabel (的吞吐量(%));标题(吞吐量和信噪比的)传说(legendString“位置”,“西北”);网格在;