NFrames = 10;SNRdB = 6.0;

enb =结构(“钢筋混凝土”那'r.3'）;％设置RMC R.3的参数eNB = ltermcdl（eNB）;enb.cfi = 3;％重新配置控制格式指示符enb.ocngpdschenable ='在';％为未分配的PDSCH RES启用OCNGenb。TotSubframes = 1;单个子帧的％重新配置enb.pdsch.rvseq = 0;％禁用HARQ.eNB.PDSCH.CSIMODE ='PUCCH 1-0';％配置CSI报告模式eNB.PDSCH.CSI ='在';软比特的％CSI缩放

channel.seed = 10;％沟道种子channel.nrxants = 2;％2接收天线channel.delayprofile ='EPA';％延迟配置文件channel.dopplerfreq = 5.0;％多普勒频率通道。MIMOCorrelation ='高的';％多天线相关性channel.modeltype =.“GMEDS”;％瑞利衰落模型类型channel.normalizetxants =.'在';传输天线的％标准化channel.normalizepathgains ='在';％标准化延迟配置文件channel.initphase =.'随机的';％随机初始阶段通道。NTerms = 16;褪色模型中使用的％振荡器％设置频道模型采样率OFDMINFO = LTEOFDMINFO（eNB）;Channel.SamplingRate = OFDMINFO.SAMPLINGRATE;

％配置信道估计器cec.pilotaverage =.'用户自定义';％类型的导频符号平均cec.freqwindow = 9;RES中的％频率窗口大小cec.timewindow = 9;％时间窗口大小cec.interpype =.'立方体';％2D插值类型cec。InterpWindow ='中心';％插值窗口类型cec.Interpwinsize = 1;％插值窗口大小％信道估计行为PerfectChanestimator = True;

cqiDelay = 8;%子帧

％初始化用于结果录制的变量cqireport = [];％报告​​的CQI值sinrreport = [];％对应的SINR值Xaxis = [];%对应的子帧号每个CQI配置的％（UE报告和中位数）：为了cqiconfig = 1：2如果(cqiConfig = = 1) cqiConfigStr ='UE报道';别的cqiConfigStr ='中位';结尾fprintf（'\以％GB SNR为％GB SNR为％D帧（s）\ n'那......cqiConfigStr、SNRdB NFrames);％初始化CQI值：对于UE报告，设置为所有人;中位数，％设置为UE报告运行的CQI值的中位数如果（CQICONFIG == 1）CQIBUFFER = ONE（1，CQIDELAY）;别的cqibuffer = =（1，cqidelay）* mediancqi;结尾%初始化变量RNG（'默认'）;%默认随机数生成器种子totalcrc = [];％CRC值，用于吞吐量计算totalTBS = [];％TBS值，用于吞吐量计算偏移= 0;％初始化帧偏移值%每个子帧:为了subframeNo = 0: (NFrames * 10 - 1)％更新子帧号enb.nsubframe = mod（subframeno，10）;％选择CQI，从CQI缓冲区读取最旧的值CQIPTR = MOD（Subframeno，Cqidelay）;cqi = cqibuffer（CQIPTR + 1）;根据CQI使用TS36.101表A.4-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 27 27];imcs = imcstable（CQI + 1）;％确定TB和调制顺序[itbs，调制] = LTEMCS（IMC）;enb.pdsch.modulation = {调制};如果(mod(enb.NSubframe,5)==0) TBS =0;别的TBS =双(lteTBS(大小(enb.PDSCH.PRBSet, 1)、胫));结尾eNB.PDSCH.trblsze（eNB.NSUBFRAME + 1）= TBS;%确定这个子帧是否需要CQI更新，根据报告周期性n_pd = 2ms和配置索引％CQI-PMI-ConfigurationIndex = 1来自TS36.101表9.3.2.1.1-1cqiperiod = 2;％周期性n_pd.cqioffset = 1;由cqi-pmi-ConfigurationIndex导出的偏移量cqiUpdate = (mod (subframeNo cqiPeriod) = = cqiOffset);％建立此子帧是否实际需要为PDSCH执行％接收，CQI估计或初始定时偏移估计：如果（（tbs〜= 0 &&子帧> =（cqidelay + cqioffset））||......（cqiconfig == 1 && cqiupdate）||subframeno == 0）％为子帧生成随机位DATA = RANDI（[0 1]，TBS，1）;％创建包含RMC传输的DM资源网格和％执行OFDM调制。txwaveform = ltermcdltool（eNB，数据）;％频道建模的初始化时间设置为每个％子帧模拟连续变化的通道channel.inittime = subframeno / 1000;将数据通过衰落信道模型。％额外的25个样本被添加到末端%波形。这是为了弥补预期的延误频道建模的％（实现的组合％延迟和渠道延迟传播）rxWaveform = lteFadingChannel(频道,......[txwaveform;零（25，尺寸（txwaveform，2））]）;％计算噪声增益，包括补偿下行链路%功率分配SNR = 10 ^（（snrdb-eNb.pdsch.rho）/ 20）;％归一化噪声功率，以考虑采样率，这％是OFDM调制中使用的IFFT大小的函数，以及％天线的数量n0 = 1 /（SQRT（2.0 * enb.cellrefp * double（Ofdminfo.nfft））* SNR）;％创建添加剂白色高斯噪声噪声= n0 *复杂（Randn（randn（rxwaveform）），......randn(大小(rxWaveform)));%将AWGN添加到接收的时域波形中rx波形= rx波形+噪声;％执行同步％在预期的延迟范围内的偏移量％信道建模（实现延迟的组合和％信道延迟扩展）表明成功如果（mod（subframeno，10）== 0）offset = ltedlframeoffset（eNB，rxwaveform）;如果(offset > 25) offset = offset (end);结尾offset = [offset];％＃OK 结尾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 &&子阵容> =（CQIDELAY + CQIOFFSET））％解码PDSCHIND = 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，PDSchchchest，噪音）;％解码DL-SCH[decbits，crc] = LTEDLSCHDECODE（eNB，eNB.PDSCH，TBS，RXBITS）;％记录CRC和TBS值的最终吞吐量％ 计算TotalCRC = [TotalCRC CRC];％＃OK totaltbs = [totaltbs tbs];％＃OK 结尾％更新CQI：如果（cqiconfig == 1 && cqiupdate）％执行CQI选择[thiscqi，thissinr] =......LTECQISELECT（eNB，eNB.PDSCH，CHESTGRID，噪音）;％将CQI值送回UE（在缓冲区中）cqibuffer（CQIPTR + 1）= Thiscqi;％记录绘图的值cqireport = [cqireport topcqi];％＃OK SINRReport = [SINRReport thisSINR];％＃OK xaxis = [xaxis subframeno];％＃OK 结尾结尾CQI未更新的子帧的百分比，重新使用上一个缓冲区中的％值如果（cqiconfig == 1 &&〜cqiupdate）cqibuffer（cqiptr + 1）= cqibuffer（mod（cqiptr-1，cqidelay）+1）;结尾结尾％显示当前CQI配置的结果fprintf（'\ nresults pyl％s cqi：\ n'，cqiconfigstr）;tputtotal = sum（totaltbs）;如果（Cqiconfig == 1）％compute和显示吞吐量tputuereported = sum（totaltbs。*（1-totalcrc））;fprintf（'吞吐量:%d位(%0.2f%%)\n'那......tputueEreported，tputuereported / tputtotal * 100）;％Compute和Display Bler测量表=均值（TotalCRC）;fprintf（'Bler：％0.3F（要求> = 0.02）\ n'，测量夹）;％compute和显示中位数cqiMediancqi = CEIL（中位数（CQIREPORT））;fprintf（“医院药学部中位数:% d \ n”，中美奇'）;CQI值的％计算和显示比例％以外的+/- 1中位数MotoottalAlpha =（Sum（CQireport <（Mediancqi-1））+......总和（cqireport>（mediancqi + 1）））/长度（cqireport）* 100;fprintf（[' CQI指数超出中位数+/- 1的百分比'......'％0.2f %%（要求> = 20 %%）\ n']，测量的alpha）;别的％compute和显示吞吐量tputMedian =总和(totalTBS。* (1-totalCRC));fprintf（'吞吐量:%d位(%0.2f%%)\n'那......tputMedian tputMedian / tputTotal * 100);％计算和显示吞吐量比测量扬声器= tputueEnterported / tprutmedian;fprintf（['吞吐率比（伽玛）：％0.3f'......'（要求> = 1.05）']，测量伽玛）;结尾结尾

数字;子图（2,1,1）;情节（Xaxis，SinRReport，'rx-'）;Xlabel（'子帧号'）;ylabel（“估计SINR (dB)”）;标题（'估计SINR与子帧号码'）;抓住在;子图（2,1,2）;plot（xaxis，cqireport，'bo-'）;Xlabel（'子帧号'）;ylabel（'选定的CQI索引'）;标题（'选定的CQI索引与子帧编号'）;