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.Modulation = {16 qam的};发射分集模式(TM2)情况下“TM2”流('\nTM2 -传输分集\n'）;simulationParameters.PDSCH.TxScheme =“TxDiversity”；simulationParameters.PDSCH.DCIFormat =“Format1”；simulationParameters。CellRefP = 2;simulationParameters.PDSCH.Modulation = {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.Modulation = {16 qam的，16 qam的};simulationParameters.PDSCH.NLayers = 2;%空间复用模式(TM4)情况下“TM4”流(基于码本的空间多路复用\n′）;simulationParameters。CellRefP = 2;simulationParameters.PDSCH.Modulation = {16 qam的，16 qam的};simulationParameters.PDSCH.DCIFormat =“Format2”；simulationParameters.PDSCH.TxScheme =“SpatialMux”；simulationParameters.PDSCH.NLayers = 2;%没有码本限制simulationParameters.PDSCH.CodebookSubset =＇＇；%单层空间复用模式(TM6)情况下“TM6”流(．..基于码本的单层空间复用\n′）;simulationParameters。CellRefP = 4;simulationParameters.PDSCH.Modulation = {“正交相移编码”};simulationParameters.PDSCH.DCIFormat =“Format2”；simulationParameters.PDSCH.TxScheme =“SpatialMux”；simulationParameters.PDSCH.NLayers = 1;%没有码本限制simulationParameters.PDSCH.CodebookSubset =＇＇；否则错误(“传输模式应该是TM1, TM2, TM3, TM4或TM6之一。”）结束%设置适用于所有TMs的其他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 =“上”；传送天线规格化

ofdmInfo = lteOFDMInfo (enb);通道。SamplingRate = ofdmInfo.SamplingRate;

完美信道估计标志perfectChanEstimator = false;

配置信道估计器cec。PilotAverage =“UserDefined”；导频符号平均值的类型cec。FreqWindow = 41;%频率窗口大小，以REs表示cec。TimeWindow = 27个;%时间窗口大小，以REs表示cec。InterpType =“立方”；二维插值型cec。InterpWindow =“中心”；插值窗口类型cec。InterpWinSize = 1;插值窗口大小

displaySimulationInformation = true;

%发射天线数量P从资源网格中获取%的维度。'dim '是M乘N乘p，其中M是子载波数N%为符号数，P为发射天线数。= lteDLResourceGridSize dim (enb);P = dim (3);%初始化在模拟和分析中使用的变量%阵列存储所有信噪比点的最大吞吐量maxThroughput = 0(长度(SNRIn), 1);%阵列存储所有信噪比点的模拟吞吐量simThroughput = 0(长度(SNRIn), 1);%获取HARQ ID序列，进行HARQ处理。这是一个索引列表%用于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:numel(SNRIn)%要启用并行计算以提高速度，请注释掉%上面的'for'语句并取消下面的'parfor'语句的注释。%这需要并行计算工具箱。如果没有安装% 'parfor'将默认为普通的'for'语句。如果“parfor”建议使用变量“displaySimulationInformation”上面的%设置为false，否则将显示模拟信息%每个信噪比点将重叠。%根据循环变量设置随机数生成器种子%以确保独立的随机流rng (snrIdx“combRecursive”）;SNRdB = SNRIn (snrIdx);流('\ n模拟%g dB信噪比为%d帧\n'、SNRdB NFrames);%初始化在模拟和分析中使用的变量偏移量= 0;%初始化帧偏移值抵消= 0;初始化无线电帧的帧偏移值blkCRC = [];%阻塞所有考虑的子帧的CRCbitTput = [];%每子帧成功接收比特数txedTrBlkSizes = [];%每子帧传输的比特数enb = enb_init;%初始化RMC配置频道= channel_init;初始化通道配置harqProcessSequence = harqProcessSequence_init;%初始化HARQ进程序列pmiIdx = 0;延迟队列中的% PMI索引变量harqPtrTable存储的是%指针指向所有HARQ进程的RV序列值。%预分配与nan，因为一些子帧没有数据rvSeqPtrHistory = NaN(ncw, NFrames*10);％HARQ进程的％初始化状态harqProcesses = hNewHARQProcess (enb);%对第一个'pmiDelay'子帧使用随机pmi，直到反馈%可从UE获得;注意PMI反馈只适用于%空间复用TMs (TM4和TM6)，但是这里的代码是必需的%表示在使用时，信噪比循环中变量的完整初始化%并行计算工具箱。pmidims = ltePMIInfo (enb enb.PDSCH);txPMIs = randi([0 pmidimms . txt])MaxPMI], pmidims。NSubbands pmiDelay);为subframeNo = 0: (NFrames * 10 - 1)%更新子帧号enb。NSubframe = subframeNo;%从HARQ进程序列中获取子帧的HARQ进程IDharqID = harqProcessSequence(mod(subframeNo, length(harqProcessSequence))+1);%当前子帧中调度了一个传输块%(由非零的'harqID'表示)，执行传输和%接待。否则继续下一帧如果harqID = = 0继续；结束%更新当前HARQ进程harqProcesses (harqID) = hHARQScheduling (．..harqProcesses (harqID)、subframeNo rvSequence);提取当前子帧传输块大小trBlk = trBlkSizes(:， mod(subframeNo, 10)+1).';%显示运行时信息如果displaySimulationInformation disp (' '）;disp ([的子帧:num2str (subframeNo)．..”。HARQ进程ID: 'num2str (harqID)]);结束%更新RV序列指针表rvSeqPtrHistory (:, subframeNo + 1) =．..harqProcesses harqID .txConfig.RVIdx。';%用HARQ进程状态更新PDSCH传输配置enb。PDSCH = harqProcesses (harqID) .txConfig;data = harqProcesses (harqID) . data;%在延迟队列中设置相应的PMI值如果strcmpi (enb.PDSCH.TxScheme“SpatialMux”) pmiIdx = mod(subframeNo, pmiDelay);延迟队列中的% PMI索引enb. pdsh . pmiset = txPMIs(:， pmiIdx+1);%设置采购经理人指数结束%产生发射波形txWaveform = lteRMCDLTool(enb, data);%添加25个样品填充。这是为了覆盖延迟的范围%期望从渠道建模(组合实现延迟和信道延迟扩展%)txWaveform = [txWaveform;0(25日P)];% #好< AGROW >%从'enbOut'获取HARQ ID序列进行HARQ处理harqProcessSequence = enbOut.PDSCH.HARQProcessSequence;%初始化每个子帧的信道时间通道。InitTime = subframeNo / 1000;%通过通道模型传递数据rx波形= lteFadingChannel(channel, tx波形);计算噪声增益，包括下行功率的补偿%分配信噪比= 10 ^ ((SNRdB-enb.PDSCH.Rho) / 20);将噪声功率归一化以考虑采样率，即%是OFDM调制中使用的IFFT大小的函数天线数%N0 = 1 /(√(2.0 * enb.CellRefP *双(nFFT)) *信噪比);创建加性高斯白噪声噪音= N0 *复杂(randn(大小(rxWaveform)),．..randn(大小(rxWaveform)));%将AWGN添加到接收的时域波形中rx波形= rx波形+噪声;在子帧0上，每帧计算一次新的同步%抵消如果(mod(subframeNo,10) == 0) offset = lteDLFrameOffset(enb, rxWaveform); / /帧偏移如果(offset > 25) offset = offset (end);结束offset = [offset];% #好吧结束%同步接收的波形rx波形= rx波形(1+偏移:结束，:);%对接收到的数据进行OFDM解调，重新创建%资源网格rxSubframe = lteofdm解调器(enb, rx波形);%信道估计如果(perfectChanEstimator) estChannelGrid = lteDLPerfectChannelEstimate(enb, channel, offset);% #好吧noiseGrid = lteOFDMDemodulate(enb, noise(1+offset:end，:));噪音= var (noiseGrid (:));其他的[estChannelGrid, noiseEst] = lteDLChannelEstimate(．..enb, enb。PDSCH、cec rxSubframe);结束%获取PDSCH索引pdschIndices = ltePDSCHIndices(enb, enb.)PDSCH enb.PDSCH.PRBSet);从接收的子帧中获取PDSCH资源元素。规模的由PDSCH功率因数Rho接收的子帧。PDSCH是的单元格参考符号%信道估计(用于PDSCH解码阶段)不是。[pdschRx, pdschHest] = lteExtractResources(pdschIndices, pdschRx, pdschHest)．..rxSubframe * (10 ^ (-enb.PDSCH.Rho / 20)), estChannelGrid);%解码PDSCHdlschBits = ltePDSCHDecode (．..enb, enb。PDSCH, pdschRx, pdschhester, noiseEst);解码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提供巡检反馈信息如果strcmpi (enb.PDSCH.TxScheme“SpatialMux”) PMI = ltePMISelect(enb, enb.)PDSCH estChannelGrid,噪音);txPMIs(:， pmiIdx+1) = PMI;结束结束计算最大和模拟吞吐量maxThroughput (snrIdx) =总和(txedTrBlkSizes);最大可能吞吐量%simThroughput (snrIdx) = (bitTput, 2)总和;%模拟吞吐量%在命令窗口中动态显示结果流([[' \ nThroughput (Mbps)的num2str (NFrames)“帧(s)”]，．..“= % .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“位置”，“西北”）;网格在；