NFrames = 10;SNRdB = 0.0;FixedRI = 2;

通道。种子= 1;%通道种子通道。NRxAnts = 2;% 2接收天线通道。DelayProfile =“环保署”;%延迟概要通道。DopplerFreq = 5.0;%多普勒频率通道。MIMOCorrelation =“低”;% Multi-antenna相关性通道。ModelType =“GMEDS”;%瑞利衰落模型类型通道。NormalizeTxAnts =“上”;%正常化传输天线通道。NormalizePathGains =“上”;%正常化延迟概要文件的权力通道。InitPhase =“随机”;%随机初始阶段通道。NTerms = 16;%振荡器用于衰落模型

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

riDelay = 9;%子帧cqipmiDelay = 8;%子帧

codebookSubsetRI1 =“000011”;codebookSubsetRI2 =“010000”;codebookSubsetUEReported =“010011”;

%检查固定RI有效1或2如果(~ (FixedRI = =(1 - 2)))错误(指定固定的RI值必须是1或2。);结束%为每个国际扶轮配置(固定RI和问题报告):为riConfig = 1:2如果(riConfig = = 1) riConfigStr = sprintf (“固定RI = % d ',FixedRI);其他的riConfigStr =“问题报告RI”;结束流(与在% % s \ nSimulating gdB信噪比% d (s) \ n帧的,…riConfigStr、SNRdB NFrames);%建立eNodeB设置2码字enb =结构(“钢筋混凝土”,“R.3”);enb。CellRefP = 2;enb。CFI = 3;enb。OCNGPDSCHEnable =“上”;enb。TotSubframes = 1;enb.PDSCH。RVSeq = 0;enb.PDSCH。CSIMode =“PUCCH 1 - 1”;enb.PDSCH。TxScheme =“SpatialMux”;enb.PDSCH。NLayers = 2;enb.PDSCH。ρ= -3.0;enb.PDSCH。CSI =“上”;作战飞机= 2;enb = lteRMCDL (enb ncw);%设置通道模型采样率ofdmInfo = lteOFDMInfo (enb);通道。SamplingRate = ofdmInfo.SamplingRate;%配置适当的约束和目标码子集%的吞吐量比如果(riConfig = = 1)如果enb.PDSCH (FixedRI = = 1)。CodebookSubset = codebookSubsetRI1;targetGamma = 1.05;其他的enb.PDSCH。CodebookSubset = codebookSubsetRI2;targetGamma = 1.0;结束其他的enb.PDSCH。CodebookSubset = codebookSubsetUEReported;结束%初始化变量totalCRC = [];% CRC值,用于吞吐量计算totalTBS = [];% TBS值,用于吞吐量计算RIReport = [];%报道RI值riXaxis = [];%对应的子帧数PMIReport = [];% PMI报告值CQIReport = [];医院药学部%报道值SINRReport = [];%的SINR值cqipmiXaxis = [];%对应的子帧数偏移量= 0;%初始化帧偏移值rng (“默认”);%默认随机数生成器的种子%初始化医院药学部RI / PMI /缓冲区;注意,PDSCH吞吐量%结果从这些初始值将被忽略,吞吐量%记录等待,直到缓冲区已满是有意义的%的报告riDelay riBuffer = 1 (1);cqipmiDelay pmiBuffer = 1 (1);cqipmiDelay cqiBuffer = 1 (ncw);%为每个子帧:为subframeNo = 0: (NFrames * 10 - 1)%更新子帧数enb。NSubframe =国防部(subframeNo 10);%更新PMI,阅读最古老的价值从PMI缓冲区cqipmiPtr =国防部(subframeNo cqipmiDelay);enb.PDSCH。PMISet = pmiBuffer (cqipmiPtr + 1);%更新数量的传输层基于国际扶轮,阅读从国际扶轮缓冲%最大的价值riPtr =国防部(subframeNo riDelay);enb.PDSCH。NLayers = riBuffer (riPtr + 1);%更新医院药学部MCS根据阅读最古老的价值医院药学部%缓冲CQI = cqiBuffer (:, cqipmiPtr + 1);%选择医院药学部MCS根据使用RMC A.4-1 CSI TS36.101表% RC.2 FDD (MCS.2),它定义了医院药学部之间的关系%指标和MCS指标IMCSTable =[1 0 0 2 4 6 8 11 13 15 18 20 22日24日26日27日);医院药学部imc = IMCSTable (+ 1);%确定TBS和调制(髂胫,调制)= lteMCS (imc);如果(mod (enb.NSubframe 5) = = 0) TBS = [0;0);其他的TBS =双(lteTBS(大小(enb.PDSCH.PRBSet, 1)、胫));结束enb.PDSCH.TrBlkSizes (:, enb.NSubframe + 1) = TBS;%确定一个医院药学部国际扶轮或PMI /更新是必需的%子帧,根据报告N_pd = 5 ms和周期性%的配置指标cqi-pmi-ConfigurationIndex = 6从表9.5.1.1-1 TS36.101 % ri-ConfigurationInd = 1riPeriod = 5;%周期性N_pdriOffset = 3;%抵消源于配置指标riUpdate = (mod (subframeNo riPeriod) = = riOffset);cqipmiPeriod = 5;%周期性N_pdcqipmiOffset = 4;%从cqi-pmi-ConfigurationIndex抵消推导cqipmiUpdate = (mod (subframeNo cqipmiPeriod) = = cqipmiOffset);%建立PDSCH如果这实际上子帧需要执行%接待,RI / PMI医院药学部/估计或初始时间偏移量%的评估:firstRxSubframe = max (cqipmiDelay + cqipmiOffset riDelay + riOffset);如果(((TBS) & & subframeNo > = firstRxSubframe) | |…(riUpdate | | cqipmiUpdate) | | subframeNo = = 0)%为传输配置数据%(1码字在1层,2层2日密语)data ={兰迪([0,1],TBS (1), 1)};如果(enb.PDSCH.NLayers = = 2) data ={1}{数据兰迪([0,1],TBS (2), 1)};结束enb.PDSCH。调制=调制(1:元素个数(数据);%创建OFDM资源包含RMC传输和网格%完成OFDM调制txWaveform = lteRMCDLTool (enb、数据);%初始化时间建模是设置每个频道%子帧来模拟连续不同的频道通道。InitTime = subframeNo / 1000;%通过衰落信道模型传递数据。%额外的25样品被添加到年底%波形。这些都是延迟预期的范围从信道建模(%的组合实现%延迟和延迟扩展频道)rxWaveform = lteFadingChannel(频道,…[txWaveform;0(25、大小(txWaveform 2))));%计算噪声增益包括下行补偿%功率分配信噪比= 10 ^ ((SNRdB-enb.PDSCH.Rho) / 20);%正常化噪声功率的采样率,%是传输线大小的函数用于OFDM调制、和%的天线数量N0 = 1 / (√(2.0 * enb.CellRefP *双(ofdmInfo.Nfft)) *信噪比);%建立加性高斯白噪声噪音= N0 *复杂(randn(大小(rxWaveform)),…randn(大小(rxWaveform)));%添加情况下接收到的时域波形rxWaveform = rxWaveform +噪声;%执行同步%的抵消延迟预期的范围内%信道建模(实现延迟和的组合%通道延迟扩展)表示成功如果(国防部(subframeNo 10) = = 0)抵消= lteDLFrameOffset (enb rxWaveform);如果(抵消> 25)抵消=补偿(结束);结束偏移量=(补偿抵消);% #好< AGROW >结束rxWaveform = rxWaveform(1 +抵消:最终,);%完成OFDM解调接收的数据来创建%接收资源网格rxGrid = lteOFDMDemodulate (enb rxWaveform);%信道估计如果(perfectChanEstimator) chEstGrid =…lteDLPerfectChannelEstimate (enb,通道,抵消);n = lteOFDMDemodulate (enb噪声(1 +抵消:最终,));噪音= var(重塑(n元素个数(n), 1));其他的[chEstGrid,噪音]=…lteDLChannelEstimate (enb enb.PDSCH,…cec、rxGrid);% #好< UNRCH >结束%如果这子帧需要PDSCH接待:如果(任何(TBS) & & subframeNo > = firstRxSubframe)%解码PDSCH印第安纳州= ltePDSCHIndices (enb enb.PDSCH enb.PDSCH.PRBSet);rxGrid pdschRx = lteExtractResources(印第安纳州)*…(10 ^ (-enb.PDSCH.Rho / 20));pdschChEst = lteExtractResources(印第安纳州,chEstGrid);[rxBits, rxSymbols] = ltePDSCHDecode (enb enb.PDSCH,…pdschRx pdschChEst,噪音);%解码DL-SCH(decbits, crc) = lteDLSCHDecode (enb enb.PDSCH,…enb.PDSCH.TrBlkSizes(1:元素个数(rxBits) enb.NSubframe + 1), rxBits);%记录CRC和TBS值最后的吞吐量%计算totalCRC = (totalCRC crc);% #好< AGROW >totalTBS = [totalTBS TBS(1:元素个数(crc)。');% #好< AGROW >结束%更新RI:如果(riUpdate)%更新国际扶轮thisRI = lteRISelect (enb enb.PDSCH chEstGrid,噪音);%饲料价值回归问题(缓冲)riBuffer (riPtr + 1) = thisRI;%绘制的记录值如果(riConfig = = 2) RIReport = [RIReport thisRI);% #好< AGROW >riXaxis = [riXaxis subframeNo];% #好< AGROW >结束结束医院药学部%更新和PMI:如果(cqipmiUpdate)%更新PMI最近报道国际扶轮的条件,%或固定国际扶轮扶轮报告不是配置如果enb.PDSCH (isempty (RIReport))。NLayers = FixedRI;其他的enb.PDSCH。NLayers = RIReport(结束);结束thisPMI = ltePMISelect (enb enb.PDSCH chEstGrid,噪音);医院药学部%更新条件在最近发布的采购经理人指数%,国际扶轮;这包括配置码字的数量基于国际扶轮%医院药学部的选择enb.PDSCH。PMISet = thisPMI;enb.PDSCH。NCodewords = enb.PDSCH.NLayers;[thisCQI, thisSINR] =…lteCQISelect (enb enb.PDSCH chEstGrid,噪音);%饲料回到问题的值(缓冲)pmiBuffer (:, cqipmiPtr + 1) = thisPMI;% CSIMode =“PUCCH 1 - 1”,医院药学部的第二码字%报告作为微分从第一个码字;%这里我们医院药学部转换为绝对指数%密语。如果(元素个数(thisCQI) = = 2) thisCQI (2) = thisCQI (1) - thisCQI (2);结束医院药学部%传输与医院药学部最低= 1时选择报告医院药学部% = 0(的范围)thisCQI (thisCQI = = 0) = 1;cqiBuffer (:, cqipmiPtr + 1) = thisCQI;%绘制的记录值如果(riConfig = = 2) PMIReport = [PMIReport thisPMI);% #好< AGROW >CQIReport = [CQIReport thisCQI (1)];% #好< AGROW >SINRReport = [SINRReport thisSINR (1)];% #好< AGROW >cqipmiXaxis = [cqipmiXaxis subframeNo];% #好< AGROW >结束结束结束%为子帧RI没有更新,不重用的前一个值%在缓冲如果(~ riUpdate) riBuffer (riPtr + 1) = riBuffer (mod (riPtr-1 riDelay) + 1);结束%的子帧PMI医院药学部/没有更新,重用以前的%值的缓冲区如果(~ cqipmiUpdate) pmiBuffer (:, cqipmiPtr + 1) =…pmiBuffer(:,国防部(cqipmiPtr-1 cqipmiDelay) + 1);cqiBuffer (:, cqipmiPtr + 1) =…cqiBuffer(:,国防部(cqipmiPtr-1 cqipmiDelay) + 1);结束结束%显示结果为当前国际扶轮的配置流(“\ % s nResults: \ n”,riConfigStr);tputTotal =总和(totalTBS);如果(riConfig = = 1)%计算和显示的吞吐量tputFixedRI =总和(totalTBS。* (1-totalCRC));流(吞吐量:% d位(% 0.2 f % %) \ n”,…tputFixedRI tputFixedRI / tputTotal * 100);其他的%计算和显示的吞吐量tputUEReported =总和(totalTBS。* (1-totalCRC));流(吞吐量:% d位(% 0.2 f % %) \ n”,…tputUEReported tputUEReported / tputTotal * 100);%计算和显示吞吐量比measuredGamma = tputUEReported / tputFixedRI;流([的吞吐量比(γ):% 0.3 f '…(要求是> = % 0.2 f)”,measuredGamma targetGamma);结束结束

图;次要情节(2,1,1);情节(riXaxis RIReport,“bo - - - - - -”);包含(“子帧数”);ylabel (“选定的国际扶轮”);标题(“选定的国际扶轮和子帧数”);持有在;次要情节(2,1,2);情节(cqipmiXaxis PMIReport,“ks -”);包含(“子帧数”);ylabel (“选择采购经理人指数”);标题(“选定的PMI和子帧数”);图;次要情节(2,1,1);情节(cqipmiXaxis SINRReport,“rx - - - - - -”);包含(“子帧数”);ylabel (“估计SINR (dB)”);标题(估计SINR(1码字)和子帧数的);持有在;次要情节(2,1,2);情节(cqipmiXaxis CQIReport,“bo - - - - - -”);包含(“子帧数”);ylabel (“医院药学部选择指数”);标题(“医院药学部选择指数(1码字)与子帧数的);