NFrames = 10;SNRdB = 6.0;

enb =结构(“钢筋混凝土”，“R.3”）;%设置RMC R.3参数enb=lteRMCDL（enb）；enb.CFI=3；%重新配置控制格式指示灯enb.OCNGPDSCHEnable=“上”；%为未分配的PDSCH资源启用OCNGenb。TotSubframes = 1;%为单个子帧重新配置enb.PDSCH.RVSeq=0；％禁用HARQ.enb.PDSCH.CSIMode='PUCCH 1-0'；%配置CSI报告模式enb.PDSCH.CSI =“上”；%软比特的CSI缩放

通道。种子= 10;%通道种子channel.nrxants = 2;% 2接收天线频道延迟配置文件=“环保署”；%延迟曲线channel.DopplerFreq=5.0；%多普勒频率通道。MIMOCorrelation =“高”；％多天线相关性channel.ModelType='gmeds'；%瑞利衰落模型类型通道。NormalizeTxAnts =“上”；%发射天线的标准化通道。NormalizePathGains =“上”；正规化延迟轮廓功率通道初始化阶段=“随机”；%随机初始相位通道。NTerms = 16;褪色模型中使用的％振荡器%设置信道模型采样率OFDMINFO = LTEOFDMINFO（eNB）;Channel.SamplingRate = OFDMINFO.SAMPLINGRATE;

配置信道估计器飞行员平均数=“UserDefined”；导频符号平均值的类型cec.FreqWindow=9；%频率窗口大小（以REs为单位）cec.TimeWindow=9；%时间窗口大小，以REs表示cec。InterpType =“立方”；%二维插值类型cec。InterpWindow =“中心”；%插值窗口类型cec.InterpWinSize=1；％插值窗口大小%信道估计器行为估计值=真；

cqiDelay = 8;%子帧

%初始化用于结果记录的变量CQIReport = [];%报告的CQI值SINRReport = [];%对应SINR值xaxis = [];%对应的子帧号对于每个CQI配置(UE报告和中位数):为cqiConfig=1:2如果(cqiConfig = = 1) cqiConfigStr =“报告问题”；其他的cqiConfigStr =“中位数”；结束流(“\n正在用%s CQI在%gdB SNR下模拟%d帧\n”，...cqiConfigStr、SNRdB NFrames);%初始化CQI值：对于报告的UE，设置为所有值；对于中值，％设置为UE报告运行的CQI值的中位数如果（CQICONFIG == 1）CQIBUFFER = ONE（1，CQIDELAY）;其他的cqiDelay cqiBuffer = 1 (1) * medianCQI;结束%初始化变量RNG（'默认'）;%默认随机数生成器种子总CRC=[]；% CRC值，用于吞吐量计算totalTBS = [];% TBS值，用于吞吐量计算偏移量=0；％初始化帧偏移值%每个子帧:为subframeNo = 0: (NFrames * 10 - 1)%更新子帧号enb。NSubframe =国防部(subframeNo 10);%选择CQI，从CQI缓冲区读取最老的值cqiPtr＝mod（子帧，cqiDelay）；CQI＝cqiBuffer（cqiPtr+1）；%使用TS36.101根据CQI选择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);％确定TB和调制顺序[ITBS，调制]=lteMCS（IMCS）；enb.PDSCH.Modulation={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_pd.cqioffset = 1;由cqi-pmi-ConfigurationIndex导出的偏移量cqiUpdate = (mod (subframeNo cqiPeriod) = = cqiOffset);%确定这个子帧是否需要执行PDSCH接收百分比，CQI估计或初始时间偏移估计:如果（（tbs〜= 0 &&子帧> =（cqidelay + cqioffset））||...(cqiConfig==1 && cqupdate) || subframeNo==0％为子帧生成随机位数据= randi([0 1]，TBS,1);％创建包含RMC传输的DM资源网格和%进行OFDM调制。txWaveform=LTERMCDLTOL（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（大小（Rx波形）），...randn(大小(rxWaveform)));％将AWGN添加到接收的时域波形rx波形= rx波形+噪声;%执行同步%在系统预期延迟范围内的偏移量%通道建模（实现延迟和%通道延迟扩展)表示成功如果（mod（子帧，10）=0）offset=lteDLFrameOffset（enb，rxWaveform）；如果(offset > 25) offset = offset (end);结束offset = [offset];%#好的结束rxWaveform = rxWaveform(1 +抵消:最终,);%对接收到的数据执行OFDM解调以创建%接收的资源网格rxGrid=LTEOFDMDdemodulate（enb，rxWaveform）；%信道估计如果(perfectChanEstimator) chEstGrid =...lteDLPerfectChannelEstimate (enb,通道,抵消);n = lteOFDMDemodulate (enb噪声(1 +抵消:最终,));噪音= var(重塑(n元素个数(n), 1));其他的[Chestgrid，噪音] =...LTEDL信道估计（enb，enb.PDSCH，...cec、rxGrid);% #好< UNRCH >结束％如果此子帧需要PDSCH接收：如果（TBS~=0&&子帧>=（cqiDelay+cqiOffset））%解码PDSCH印第安纳州= ltePDSCHIndices (enb enb.PDSCH enb.PDSCH.PRBSet);pdschRx = lteExtractResources(ind,rxGrid) * . pdschRx = lteExtractResources(ind,rxGrid) *...(10 ^ (-enb.PDSCH.Rho / 20));pdschChEst = lteExtractResources(印第安纳州,chEstGrid);[rxBits, rxSymbols] = ltePDSCHDecode (enb, enb。PDSCH,...pdschRx、PDSCHCURT、noiseEst）；解码DL-SCH[decbits，crc] = LTEDLSCHDECODE（eNB，eNB.PDSCH，TBS，RXBITS）;％记录CRC和TBS值的最终吞吐量%计算totalCRC=[totalCRC]；%#好的totalTBS = [totalTBS TBS];%#好的结束％更新CQI：如果（cqiconfig == 1 && cqiupdate）%执行CQI选择[thisCQI, thisSINR] =...lteCQISelect (enb enb.PDSCH chEstGrid,噪音);%将CQI值反馈回UE（在缓冲区中）Cqibuffer（CQIPTR + 1）= Thiscqi;%记录绘图值CQIReport=[CQIReport thisCQI]；%#好的SINRReport = [SINRReport thisSINR];%#好的xaxis=[xaxis子帧]；%#好的结束结束CQI未更新的子帧的百分比，重新使用上一个%缓冲区中的值如果（cqiConfig==1&&~cqiUpdate）cqiBuffer（cqiPtr+1）=cqiBuffer（mod（cqiPtr-1，cqiDelay）+1）；结束结束%显示当前CQI配置的结果流('\nResults with %s CQI:\n'，cqiConfigStr）；tputTotal=总和（总TBS）；如果（cqiConfig==1）%计算和显示吞吐量tputUEReported=sum（totalTBS.*（1-totalCRC））；fprintf('吞吐量:%d位(%0.2f%%)\n'，...tputUEReported，tputUEReported/tputTotal*100）；%计算和显示BLER测量布尔=平均值（总CRC）；fprintf('Bler：％0.3F（要求> = 0.02）\ n'，被测量者）；%计算并显示中值CQImedianCQI =装天花板(中位数(CQIReport));流(“医院药学部中位数:% d \ n”中成药'；%计算并显示CQI值的比例%中间带+/-1之外MotoottalAlpha =（Sum（CQireport <（Mediancqi-1））+...总和（CQIReport>（medianCQI+1））/长度（CQIReport）*100；fprintf([' CQI指数超出中位数+/- 1的百分比'...%0.2f%%(要求是>= 20%%)\n']，被测α）；其他的%计算和显示吞吐量tputMedian =总和(totalTBS。* (1-totalCRC));流('吞吐量:%d位(%0.2f%%)\n'，...tputMedian tputMedian / tputTotal * 100);%计算和显示吞吐量比率测量扬声器= tputueEnterported / tprutmedian;fprintf（['吞吐量比率（伽马）：%0.3f'...'（要求> = 1.05）']，测量伽玛）;结束结束

图；子地块（2,1,1）；地块（xaxis、SINRReport、，“rx-”）;包含(“子帧数”)；伊莱贝尔(“估计SINR (dB)”）;标题（“估计SINR与子帧数”); 持有在; 子批次（2,1,2）；绘图（xaxis、CQIReport、，“bo - - - - - -”）;包含(“子帧数”)；伊莱贝尔(“选定CQI索引”）;标题（'选定的CQI索引与子帧编号'）;