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 = {'16QAM'};发射分集模式(TM2)案件“TM2”流('\nTM2 -传输分集\n'）;simulationParameters.PDSCH.TxScheme =“TxDiversity”；simulationParameters.PDSCH.DCIFormat ='format1'；SimulationParameters.cellrefp = 2;simulationParameters.PDSCH.Modulation = {'16QAM'};simulationParameters.PDSCH.NLayers = 2;% CDD模式(TM3)案件“TM3”流(“\ nTM3 CDD \ n”）;simulationParameters.PDSCH.TxScheme ='CDD'；simulationParameters.PDSCH.DCIFormat =“Format2A”；SimulationParameters.cellrefp = 2;simulationParameters.PDSCH.Modulation = {'16QAM'，'16QAM'};simulationParameters.PDSCH.NLayers = 2;％空间多路复用模式（TM4）案件“TM4”流('\ ntm4  - 基于码本的空间多路复用\ n'）;SimulationParameters.cellrefp = 2;simulationParameters.PDSCH.Modulation = {'16QAM'，'16QAM'};simulationParameters.PDSCH.DCIFormat ='format2'；simulationParameters.PDSCH.TxScheme ='spatialmux'；simulationParameters.PDSCH.NLayers = 2;%没有码本限制simulationParameters.PDSCH.CodebookSubset =＇＇；%单层空间复用模式(TM6)案件“TM6”流(......'\ ntm6  - 基于码本的空间复用与单层\ 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中的一种。）结束％设置应用于所有TM的其他模拟参数字段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;%通道种子channel.nrxants = 2;% 2接收天线通道。DelayProfile =“环保署”；%延迟概要channel.dopplerfreq = 5;%多普勒频率channel.mimocorrelation =.'低的'；％多天线相关性channel.nterms = 16;衰落模型中使用的振荡器通道。ModelType ='gmeds'；%瑞利衰落模型类型通道。InitPhase =“随机”；%随机初始阶段通道。NormalizePathGains =“上”；正规化延迟轮廓功率通道。NormalizeTxAnts =“上”；传送天线规格化％使用返回的值设置通道模型的采样率% from 。OFDMINFO = LTEOFDMINFO（eNB）;Channel.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从资源网格中获取％ 方面。'dims'是m-by-n by-p，其中m是子载波的数量，n%为符号数，P为发射天线数。dims = ltedlresourcegridsize（eNB）;p =昏暗（3）;%初始化在模拟和分析中使用的变量%阵列存储所有信噪比点的最大吞吐量maxThroughput = 0(长度(SNRIn), 1);％数组存储所有SNR点的模拟吞吐量simthrougupput = zeros（长度（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;letcendstring = [的吞吐量:char (enb.PDSCH.TxScheme)];allRvSeqPtrHistory =细胞(1,元素个数(SNRIn));nFFT = ofdmInfo.Nfft;为snrIdx = 1:元素个数(SNRIn)% parfor snrIdx = 1:numel(SNRIn)％要使用并行计算的使用，以增加速度评论%上面的'for'语句并取消下面的'parfor'语句的注释。%这需要并行计算工具箱。如果没有安装％'parcon'默认为正常的'for'语句。如果是'parfor'百分比使用它建议变量'displaysimulationInformation'上面的%设置为false，否则将显示模拟信息％每个SNR点都将重叠。%根据循环变量设置随机数生成器种子%以确保独立的随机流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），但此处的代码是必需的使用时为SNR循环中的变量完全初始化的％%并行计算工具箱。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;数据= HarqProcesses（Harqid）.data;％将PMI设置为延迟队列中的适当值如果Strcmpi（eNB.PDSCH.TXScheme，'spatialmux'）PMIIDX​​ = MOD（子阵容，PMIDELAY）;延迟队列中％PMI索引enb. pdsh . pmiset = txPMIs(:， pmiIdx+1);％设置PMI结束％生成传输波形txWaveform = lteRMCDLTool(enb, data);％加入25个样品填料。这是为了涵盖延迟范围%期望从渠道建模(组合实现延迟和信道延迟扩展%)txwaveform = [txwaveform;零（25，p）];% #好< AGROW >%从'enbOut'获取HARQ ID序列进行HARQ处理harqProcessSequence = enbOut.PDSCH.HARQProcessSequence;％初始化每个子帧的通道时间通道。InitTime = subframeNo / 1000;％通过频道模型传递数据rxwaveform = ltefadingchannel（通道，txwaveform）;计算噪声增益，包括下行功率的补偿%分配信噪比= 10 ^ ((SNRdB-enb.PDSCH.Rho) / 20);将噪声功率归一化以考虑采样率，即%是OFDM调制中使用的IFFT大小的函数天线数%n0 = 1 /（SQRT（2.0 * enb.cellrefp * double（nfft））* snr）;创建加性高斯白噪声噪声= n0 *复杂（Randn（randn（rxwaveform）），......Randn（尺寸（rxwaveform））））;％将AWGN添加到接收的时域波形rxwaveform = rxwaveform +噪声;在子帧0上，每帧计算一次新的同步％ 抵消如果（mod（subframeno，10）== 0）offset = ltedlframeoffset（eNB，rxwaveform）;如果（offset> 25）offset =偏移（结束）;结束偏移= [偏移偏移];％＃行结束%同步接收的波形rxwaveform = rxwaveform（1 + offset：结束，:);%对接收到的数据进行OFDM解调，重新创建%资源网格rxSubframe = lteofdm解调器(enb, rx波形);％信道估计数如果(perfectChanEstimator) estChannelGrid = lteDLPerfectChannelEstimate(enb, channel, offset);％＃行noiseGrid = lteOFDMDemodulate(enb, noise(1+offset:end，:));噪音= var (noiseGrid (:));其他的[EstchannelGrid，噪声] = LTedLchanneLeStimate（......eNB，eNB.PDSCH，CEC，RXSUBFRAME）;结束％获得PDSCH指数pdschindices = ltepdschindices（eNB，enb.pdsch，enb.pdsch.prbset）;％从接收的子帧获取PDSCH资源元素。缩放由PDSCH功率因数Rho接收的子帧。PDSCH是的单元格参考符号%信道估计(用于PDSCH解码阶段)不是。[pdschrx，pdschhest] = lteextractresources（pdschindices，......rxsubframe *（10 ^（ -  enb.pdsch.rho / 20）），Estchannelgrid）;％解码PDSCH.dlschbits = 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;结束结束％计算最大和模拟吞吐量maxthrougpul（snridx）= sum（txedtrblksizes）;％最大吞吐量simThroughput (snrIdx) = (bitTput, 2)总和;％模拟吞吐量%在命令窗口中动态显示结果fprintf（[[' \ nThroughput (Mbps)的num2str (NFrames)“帧(s)”]，......“= % .4f \ n”), 1 e-6 * simThroughput (snrIdx) / (NFrames * 10 e - 3));流([的吞吐量(% %)num2str (NFrames)'帧(s) = %.4f\n']，......Simkrouguppul（Snridx）* 100 / maxthrougpul（snridx））;allrvseqptrhistory {snridx} = rvseqptrhistory;结束%绘制所有HARQ进程的RV序列hplotrvsequence（snrin，allrvseqptrhistory，nframes）;

%绘制吞吐量图绘制(SNRIn simThroughput * 100. / maxThroughput,* - - - - - -。）;包含('snr（db）'）;ylabel ('吞吐量（％）'）;标题（'吞吐量与snr')传说(legendString“位置”，“西北”）;网格在；