NFrames = 10;SNRdB = 6.0;

enb =结构(“钢筋混凝土”，“R.3”）;%设置RMC r的参数enb = lteRMCDL (enb);enb。CFI = 3;%重新配置控制格式指示器enb。OCNGPDSCHEnable =“上”;%为未分配的PDSCH REs启用OCNGenb。TotSubframes = 1;%重新配置单个子帧enb.PDSCH.RVSeq = 0;%禁用HARQenb.PDSCH.CSIMode =“PUCCH 1 - 0”;%配置CSI上报方式enb.PDSCH.CSI =“上”;软位的% CSI缩放

通道。种子= 10;%通道种子通道。NRxAnts = 2;% 2接收天线通道。DelayProfile =“环保署”;%延迟概要通道。DopplerFreq = 5.0;%多普勒频率通道。MIMOCorrelation =“高”;% Multi-antenna相关性通道。ModelType =“GMEDS”;%瑞利衰落模型类型通道。NormalizeTxAnts =“上”;%对发射天线进行归一化通道。NormalizePathGains =“上”;%归一化延迟配置文件功率通道。InitPhase =“随机”;%随机初始阶段通道。NTerms = 16;衰落模型中使用的振荡器%设置通道模型采样率ofdmInfo = lteOFDMInfo (enb);通道。SamplingRate = ofdmInfo.SamplingRate;

%配置信道估计器cec。PilotAverage =“UserDefined”;%导频符号平均类型cec。FreqWindow = 9;% REs中的频率窗口大小cec。TimeWindow = 9;时间窗口大小(REs)cec。InterpType =“立方”;% 2D插补型cec。InterpWindow =“中心”;%插值窗口类型cec。InterpWinSize = 1;%插值窗口大小%信道估计器行为perfectChanEstimator = true;

cqiDelay = 8;%子帧

初始化用于结果记录的变量CQIReport = [];报告的CQI值SINRReport = [];%对应的SINR值xaxis = [];%对应的子帧编号%对于每个CQI配置(UE报告和中位数):为cqiConfig = 1:2如果(cqiConfig = = 1) cqiConfigStr =“报告问题”;其他的cqiConfigStr =“中值”;结束流(' nsimulation with %s CQI at %gdB SNR for %d Frame(s)\n'，.．.cqiConfigStr、SNRdB NFrames);%初始化CQI值:对于报告的UE，设置为所有1;中位数,%设置为报告运行的UE的CQI值的中位数如果cqiConfig==1) cqbuffer = ones(1, cqiddelay);其他的cqiDelay cqiBuffer = 1 (1) * medianCQI;结束%初始化变量rng (“默认”）;%默认随机数生成器种子totalCRC = [];% CRC值，用于吞吐量计算totalTBS = [];% TBS值，用于吞吐量计算偏移量= 0;%初始化帧偏移值%每个子帧:为subframeNo = 0: (NFrames * 10 - 1)%更新子帧编号enb。NSubframe =国防部(subframeNo 10);%选择CQI，从CQI缓冲区中读取最老的值cqiPtr =国防部(subframeNo cqiDelay);CQI = cqiBuffer (cqiPtr + 1);表a . -1 CSI RMC . MCS% RC.1 FDD (MCS.1)，它定义了CQI之间的关系%指数和MCS指数IMCSTable = [-1 0 0 2 4 6 8 11 13 16 18 21 23 25 27 27];医院药学部imc = IMCSTable (+ 1);%确定TBS和调制顺序(髂胫,调制)= lteMCS (imc);enb.PDSCH.Modulation ={调制};如果(mod(enb.NSubframe,5)==0) TBS =0;其他的TBS =双(lteTBS(大小(enb.PDSCH.PRBSet, 1)、胫));结束enb.PDSCH.TrBlkSizes (enb.NSubframe + 1) = TBS;根据，确定该子帧中是否需要CQI更新%到报告周期性N_pd = 2ms和配置索引% cqi-pmi-ConfigurationIndex = 1 from TS36.101cqiPeriod = 2;%周期性N_pdcqiOffset = 1;从cqi-pmi-ConfigurationIndex派生的%偏移量cqiUpdate = (mod (subframeNo cqiPeriod) = = cqiOffset);%确定PDSCH是否需要执行此子帧%接收，CQI估计或初始定时偏移估计:如果((TBS~=0 && subframeNo>=(cqiddelay +cqiOffset)) ||。.．.(cqiConfig==1 && cqiUpdate) || subframe ==0)%为子帧生成随机位数据= randi([0 1]，TBS,1);%创建包含RMC传输的OFDM资源网格%执行OFDM调制。txWaveform = lteRMCDLTool (enb、数据);%信道建模的初始化时间被每个设置%子帧来模拟连续变化的信道通道。InitTime = subframeNo / 1000;%通过衰落信道模型传递数据。另外的25个样品被添加到样品的末端%波形。这些是为了弥补预期的延误范围%从通道建模(实现的组合%延迟和信道延迟扩展)rxWaveform = lteFadingChannel(频道,.．.[txWaveform;0(25、大小(txWaveform 2))));%计算噪声增益，包括下行链路补偿%功率分配信噪比= 10 ^ ((SNRdB-enb.PDSCH.Rho) / 20);%归一化噪声功率，以考虑采样率%是OFDM调制中IFFT大小的函数天线的数量N0 = 1 /(√(2.0 * enb.CellRefP *双(ofdmInfo.Nfft)) *信噪比);%添加高斯白噪声噪音= N0 *复杂(randn(大小(rxWaveform)),.．.randn(大小(rxWaveform)));%在接收到的时域波形中加入AWGNrx波形= rx波形+噪声;%执行同步的延迟范围内的偏移量%通道建模(实现延迟和% channel delay spread)表示成功如果(mod(subframe,10)==0) offset = lteDLFrameOffset(enb, rx波形);如果offset = offset (end);结束offset = [offset offset];% #好< 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 ~ = 0 & & subframeNo > = (cqiDelay + cqiOffset))%解码PDSCH印第安纳州= ltePDSCHIndices (enb enb.PDSCH enb.PDSCH.PRBSet);pdschRx = lteExtractResources(ind,rxGrid) *.．.(10 ^ (-enb.PDSCH.Rho / 20));pdschChEst = lteExtractResources(印第安纳州,chEstGrid);[rxBits, rxSymbols] = ltePDSCHDecode (enb, enb。PDSCH,.．.pdschRx pdschChEst,噪音);解码DL-SCH(decbits, crc) = lteDLSCHDecode (enb.PDSCH, enb TBS, rxBits);%记录最终吞吐量的CRC和TBS值%计算totalCRC = [totalCRC crc];% #好< AGROW >totalTBS = [totalTBS TBS];% #好< AGROW >结束%医院药学部更新:如果(& & cqiUpdate cqiConfig = = 1)%执行CQI选择[thisCQI, thisSINR] =.．.lteCQISelect (enb enb.PDSCH chEstGrid,噪音);将CQI值返回到UE(在缓冲区中)cqiBuffer (cqiPtr + 1) = thisCQI;%记录用于绘图的值CQIReport = [CQIReport thisCQI];% #好< AGROW >SINRReport = [SINRReport thisSINR];% #好< AGROW >xaxis = [xaxis subframe];% #好< AGROW >结束结束%对于CQI未更新的子帧，重用之前的子帧%值如果cqiConfig==1 && ~cqiUpdate cqiBuffer(cqiPtr+1) = cqiBuffer(mod(cqiPtr-1, cqiddelay)+1);结束结束%显示当前CQI配置结果流('\nResults with %s CQI:\n', cqiConfigStr);tputTotal =总和(totalTBS);如果(cqiConfig = = 1)%计算和显示吞吐量tputUEReported =总和(totalTBS。* (1-totalCRC));流('吞吐量:%d位(%0.2f%%)\n'，.．.tputUEReported tputUEReported / tputTotal * 100);%计算和显示BLERmeasuredBLER =意味着(totalCRC);流('BLER: %0.3f (requirement is >= 0.02)\n', measuredBLER);%计算和显示中间CQImedianCQI =装天花板(中位数(CQIReport));流(“医院药学部中位数:% d \ n”medianCQI ');%计算和显示CQI值的比例中位数的+/- 1measuredAlpha = (sum(CQIReport<(medianCQI-1)) +.．.sum (CQIReport > (medianCQI + 1))) /长度(CQIReport) * 100;流([中位数+/- 1以外的CQI指数百分比:.．.' %0.2f%%(要求是>= 20%%)\n'), measuredAlpha);其他的%计算和显示吞吐量tputMedian =总和(totalTBS。* (1-totalCRC));流('吞吐量:%d位(%0.2f%%)\n'，.．.tputMedian tputMedian / tputTotal * 100);%计算和显示吞吐量比率measuredGamma = tputUEReported / tputMedian;流(['吞吐量比(伽马):%0.3f'.．.'(要求>= 1.05)'), measuredGamma);结束结束

图;次要情节(2,1,1);情节(xaxis SINRReport,“rx - - - - - -”）;包含(“子帧数”）;ylabel (“估计SINR (dB)”）;标题(“估计SINR与子帧数”）;持有在;次要情节(2,1,2);情节(xaxis CQIReport,“bo - - - - - -”）;包含(“子帧数”）;ylabel (“医院药学部选择指数”）;标题('选定的CQI索引与子帧数'）;