NFrames = 2;帧数%SNRIn = [10.3 12.3 14.3];%信噪比范围，单位为dB

txMode =“TM4”；% tm1, tm2, tm3, tm4, tm6

simulationParameters = [];%清除模拟参数simulationParameters。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 .调制= {16 qam的};发射分集模式(TM2)情况下“TM2”流(\nTM2 -发送分集\n）;simulationParameters.PDSCH.TxScheme =“TxDiversity”；simulationParameters.PDSCH.DCIFormat =“Format1”；simulationParameters。CellRefP = 2;simulationparameters . pdsch .调制= {16 qam的};simulationParameters.PDSCH.NLayers = 2;% CDD模式(TM3)情况下“TM3”流('\nTM3 - CDD\n'）;simulationParameters.PDSCH.TxScheme =CDD的；simulationParameters.PDSCH.DCIFormat =“Format2A”；simulationParameters。CellRefP = 2;simulationparameters . pdsch .调制= {16 qam的，16 qam的};simulationParameters.PDSCH.NLayers = 2;空间多路复用模式(TM4)情况下“TM4”流(\nTM4 -基于码本的空间多路复用\n）;simulationParameters。CellRefP = 2;simulationparameters . pdsch .调制= {16 qam的，16 qam的};simulationParameters.PDSCH.DCIFormat =“Format2”；simulationParameters.PDSCH.TxScheme =“SpatialMux”；simulationParameters.PDSCH.NLayers = 2;%无码本限制simulationParameters.PDSCH.CodebookSubset =”；单层空间多路复用模式(TM6)情况下“TM6”流(.．.\nTM6 -基于码本的单层空间多路复用\n）;simulationParameters。CellRefP = 4;simulationparameters . pdsch .调制= {“正交相移编码”};simulationParameters.PDSCH.DCIFormat =“Format2”；simulationParameters.PDSCH.TxScheme =“SpatialMux”；simulationParameters.PDSCH.NLayers = 1;%无码本限制simulationParameters.PDSCH.CodebookSubset =”；否则错误(“传输模式应该是TM1、TM2、TM3、TM4或TM6中的一种。”）结束%设置适用于所有tm的其他simulationParameters字段simulationParameters。TotSubframes = 1;每次生成一个子帧simulationParameters.PDSCH.CSI =“上”；%软位由CSI加权

enb = lteRMCDL(simulationParameters);

rvSequence = enb.PDSCH.RVSeq;trBlkSizes = enb.PDSCH.TrBlkSizes;

ncw = length(字符串(enb. pdsch .调制));

pmiddelay = 8;

hDisplayENBParameterSummary (enb txMode);

通道。种子= 6;%沟道种子通道。NRxAnts = 2;% 2接收天线通道。DelayProfile =“环保署”；%延迟配置文件通道。DopplerFreq = 5;%多普勒频率通道。MIMOCorrelation =“低”；%多天线相关性通道。NTerms = 16;衰落模型中使用的振荡器通道。ModelType =“GMEDS”；瑞利衰落模型类型通道。InitPhase =“随机”；%随机初始相通道。NormalizePathGains =“上”；正常化延迟配置文件功率通道。NormalizeTxAnts =“上”；用于发射天线归一化

ofdmInfo = lteOFDMInfo(enb);通道。SamplingRate = ofdmInfo.SamplingRate;最大信道延迟(多路径和信道滤波)通道。InitTime = 0;[~，chInfo] = lteFadingChannel(channel,0);获取频道信息maxChDelay = ceil(max(chinfo . pathsampledelayed)) + chInfo.ChannelFilterDelay;

%完美信道估计器标志perfectChanEstimator = false;

配置信道估计器cec。PilotAverage =“UserDefined”；导频符号平均的类型cec。FreqWindow = 41;%频率窗口大小在REscec。TimeWindow = 27;%时间窗口大小在REscec。InterpType =“立方”；% 2D插值类型cec。InterpWindow =“中心”；插补窗口类型cec。InterpWinSize = 1;插值窗口大小%

displaySimulationInformation = true;

%从资源网格中获取发射天线个数P%的维度。'dims'是M乘N乘p其中M是子载波数N%为符号数，P为发射天线数。dims = lteDLResourceGridSize(enb);P = dims(3);初始化模拟和分析中使用的变量数组来存储所有信噪比点的最大吞吐量maxThroughput = 0(长度(SNRIn)，1);数组来存储所有信噪比点的模拟吞吐量simThroughput = 0(长度(SNRIn)，1);获取HARQ处理的HARQ ID序列。这是一个索引列表%用于HARQ进程调度。[~，~，enbOut] = lteRMCDLTool(enb， []);harqProcessSequence = en布特. pdsch . harqProcessSequence;临时变量“enb_init”，“channel_init”和'harqProcessSequence_init'用于优化并行处理(仅%如果使用并行计算工具箱运行示例)Enb_init = enb;Channel_init = channel;harqProcessSequence_init = harqProcessSequence;legendString = [的吞吐量:char (enb.PDSCH.TxScheme)];allRvSeqPtrHistory = cell(1, number (SNRIn));nFFT = ofdmInfo.Nfft;为snrIdx = 1:数字(SNRIn)% parfor snrIdx = 1:数值(SNRIn)允许使用并行计算，以提高注释输出的速度%上面的'for'语句，并取消下面的'parfor'语句的注释。这需要并行计算工具箱。如果没有安装% 'parfor'将默认为普通的'for'语句。如果'parfor'是建议变量'displaySimulationInformation'%以上设置为false，否则模拟信息显示为%每个信噪比点都会重叠。根据循环变量设置随机数生成器种子%，确保独立的随机流rng (snrIdx“combRecursive”）;SNRdB = SNRIn(snrIdx);流(在%g dB信噪比下模拟%d帧(s)\n'、SNRdB NFrames);初始化模拟和分析中使用的变量偏移量= 0;初始化帧偏移值偏移量= 0;初始化无线电帧的帧偏移值blkCRC = [];所有考虑的子帧的块CRCbitTput = [];每个子帧成功接收的比特数txedTrBlkSizes = [];%每子帧传输的比特数Enb = enb_init;初始化RMC配置Channel = channel_init;初始化通道配置harqProcessSequence = harqProcessSequence_init;初始化HARQ进程序列pmiIdx = 0;延迟队列中的PMI索引变量harqPtrTable存储的历史值%指针指向所有HARQ进程的RV序列值。%由于有些子帧没有数据，所以预分配给nanrvSeqPtrHistory = NaN(ncw, NFrames*10);初始化所有HARQ进程的状态harqProcesses = hNewHARQProcess(enb);使用随机pmi的第一' pmiddelay '子帧，直到反馈UE可用百分比;请注意，PMI反馈仅适用于%空间多路复用tm (TM4和TM6)，但是这里的代码是必需的%用于信噪比循环中变量的完全初始化并行计算工具箱。pmidims = ltePMIInfo(enb,enb. pdsch);txPMIs = randi([0 pmidims.]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;data = harqProcesses(harqID).data;%在延迟队列中设置合适的PMI值如果strcmpi (enb.PDSCH.TxScheme“SpatialMux”) pmiIdx = mod(subframeNo, pmiddelay);延迟队列中的PMI索引enb.PDSCH.PMISet = txPMIs(:， pmiIdx+1);%设置PMI结束生成发射波形tx波形= lteRMCDLTool(enb, data);添加maxChDelay示例填充。这是覆盖的范围%的信道建模预期延迟(组合为%实现延迟和信道延迟扩展)tx波形= [tx波形;0 (maxChDelay P));% #好< AGROW >从'enbOut'获取HARQ ID序列用于HARQ处理harqProcessSequence = en布特. pdsch . harqProcessSequence;初始化每个子帧的通道时间通道。InitTime = subframeNo/1000;通过通道模型传递数据rx波形= lteFadingChannel(通道，tx波形);计算噪声增益，包括下行功率补偿%分配信噪比= 10^((SNRdB-enb.PDSCH.Rho)/20);将噪声功率归一化以考虑采样率，即%是OFDM调制中使用的IFFT大小的函数天线数%N0 = 1/(√(2.0*enb.CellRefP*double(nFFT))*信噪比);创建加性高斯白噪声噪声= N0*complex(randn(size(rx波形))，.．.randn(大小(rxWaveform)));将AWGN添加到接收的时域波形中rx波形= rx波形+噪声;每一帧，在子帧0上计算一次新的同步%抵消如果(mod(subframameno,10) == 0) offset = lteDLFrameOffset(enb, rx波形);如果(offset > maxChDelay) offset = offset (end);结束offset = [offset offset];% #好吧结束同步接收到的波形rx波形= rx波形(1+offset:end，:);对接收到的数据进行OFDM解调以重建数据%资源网格rxSubframe = lteOFDMDemodulate(enb, rx波形);%信道估计如果(perfectChanEstimator) estChannelGrid = lteDLPerfectChannelEstimate(enb, channel, offset);% #好吧noiseGrid = lteOFDMDemodulate(enb, noise(1+offset:end，:));noiseEst = var(noiseGrid(:));其他的[estChannelGrid, noiseEst] = lteDLChannelEstimate(.．.enb, enb。PDSCH, cec, rxSubframe);结束%获取PDSCH指数pdschIndices = ltePDSCHIndices(enb, enb。PDSCH enb.PDSCH.PRBSet);从接收的子帧中获取PDSCH资源元素。规模的接收子帧百分比由PDSCH功率因数Rho。PDSCH是%按此量缩放，而单元格引用符号用于%信道估计(用于PDSCH解码阶段)则不是。[pdschRx, pdschHest] = lteExtractResources(pdschIndices, pdschIndices，.．.rxSubframe * (10 ^ (-enb.PDSCH.Rho / 20)), estChannelGrid);解码PDSCHdlschBits = ltePDSCHDecode(.．.enb, enb。PDSCH, pdschRx, pdschest, 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)])结束结束存储用于计算吞吐量的值仅适用于有数据的子帧如果any(trBlk) blkCRC = [blkCRC harqProcesses(harqID).blkerr];% #好< AGROW >bitTput = [bitTput trBlk.*(1- . ..．.harqProcesses (harqID) .blkerr)];% #好< AGROW >txedTrBlkSizes = [txedTrBlkSizes trBlk];% #好< AGROW >结束向eNodeB提供PMI反馈如果strcmpi (enb.PDSCH.TxScheme“SpatialMux”) PMI = ltePMISelect(enb, enb.)PDSCH, estChannelGrid, noiseEst);txpmi (:， pmiIdx+1) = PMI;结束结束计算最大吞吐量和模拟吞吐量maxThroughput(snrIdx) = sum(txedTrBlkSizes);最大可能吞吐量simThroughput(snrIdx) = sum(bitTput,2);%模拟吞吐量在命令窗口中动态显示结果流([[\ n吞吐量(Mbps)num2str (NFrames)“帧”]，.．.“= % .4f \ n”), 1 e-6 * simThroughput (snrIdx) / (NFrames * 10 e - 3));流(['吞吐量(%%)'num2str (NFrames)'帧(s) = %.4f\n']，.．.simThroughput (snrIdx) * 100 / maxThroughput (snrIdx));allRvSeqPtrHistory{snrIdx} = rvSeqPtrHistory;结束绘制所有HARQ进程的RV序列hPlotRVSequence (SNRIn allRvSeqPtrHistory NFrames);

地块吞吐量(SNRIn, simThroughput*100./maxThroughput，* - - - - - -。）;包含(“信噪比(dB)”）;ylabel (的吞吐量(%)）;标题(“吞吐量vs.信噪比”)传说(legendString“位置”，“西北”）;网格在；