NFrames = 2;%的帧数SNRIn = 1 [4];在dB %信噪比范围

simulationParameters = [];%明显simulationParameterssimulationParameters。NDLRB = 50;simulationParameters.PDSCH。TargetCodeRate = 0.5;simulationParameters.PDSCH。调制= {16 qam的};simulationParameters.PDSCH。PRBSet = (0:9) ';txMode =“TM8”;% TM7, TM8 TM9 TM10DuplexMode =“TDD”;%的FDD, TDD的开关txMode情况下“TM7”%单层(端口5)simulationParameters.PDSCH。DCIFormat =“Format1”;simulationParameters.PDSCH。TxScheme =“Port5”;simulationParameters.PDSCH。NLayers = 1;ntxants = 4;情况下“TM8”%两层(双港7 - 8层,或单一的天线%端口,端口7或8)simulationParameters.PDSCH。DCIFormat =“Format2B”;simulationParameters.PDSCH。TxScheme =“Port7-8”;%’Port7-8”、“Port8”simulationParameters.PDSCH。NLayers = 2;ntxants = 4;simulationParameters.PDSCH。NSCID = 0;情况下“TM9”% 8层(八层端口7 - 14或单身%天线端口,端口7或8)simulationParameters.PDSCH。DCIFormat =“Format2C”;simulationParameters.PDSCH。TxScheme =“Port7-14”;simulationParameters.PDSCH。NLayers = 4;ntxants = 8;simulationParameters。CSIRefP = ntxants;情况下“TM10”% 8层(端口7 - 14或8层%单天线端口,端口7或8)simulationParameters.PDSCH。DCIFormat =“Format2D”;simulationParameters.PDSCH。TxScheme =“Port7-14”;simulationParameters.PDSCH。NLayers = 4;ntxants = 8;simulationParameters。CSIRefP = ntxants;否则错误(的传输方式应该是TM7, TM8, TM9或TM10”。)结束% 2码字(超过4层),调制领域需要两个值,每一个码字如果(simulationParameters.PDSCH。NLayers> 4) && (length(simulationParameters.PDSCH.Modulation)<2) simulationParameters.PDSCH.Modulation = repmat(simulationParameters.PDSCH.Modulation,1,2);结束

%初始化W为零simulationParameters.PDSCH。W= zeros(simulationParameters.PDSCH.NLayers, ntxants);

simulationParameters。DuplexMode =DuplexMode; simulationParameters.TotSubframes = 1;

enb = lteRMCDL (simulationParameters);

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

作战飞机=长度(string (enb.PDSCH.Modulation));

hDisplayENBParameterSummary (enb txMode);

频道=结构;%通道配置结构通道。种子= 6;%通道种子通道。NRxAnts = 8;%的接收天线通道。DelayProfile =“环保署”;%延迟概要通道。DopplerFreq = 5;%多普勒频率通道。MIMOCorrelation =“低”;% Multi-antenna相关性通道。NTerms = 16;%振荡器用于衰落模型通道。ModelType =“GMEDS”;%瑞利衰落模型类型通道。InitPhase =“随机”;%随机初始阶段通道。NormalizePathGains =“上”;%正常化延迟概要文件的权力通道。NormalizeTxAnts =“上”;%正常化传输天线%通道采样率取决于用于OFDM FFT大小%调制器。这可以使用函数lteOFDMInfo获得。ofdmInfo = lteOFDMInfo (enb);通道。SamplingRate = ofdmInfo.SamplingRate;%最大信道延迟(多路径和通道过滤)通道。InitTime = 0;[~,chInfo] = lteFadingChannel(频道,0);%获得信道信息maxChDelay =装天花板(max (chInfo.PathSampleDelays)) + chInfo.ChannelFilterDelay;

cec =结构;%信道估计配置结构cec。PilotAverage =“UserDefined”;%的试点,平均象征cec。FreqWindow = 41;%频率窗口大小cec。TimeWindow = 27个;%的时间窗口大小cec。InterpType =“立方”;% 2 d插值类型cec。InterpWindow =“中心”;%插值窗口类型cec。InterpWinSize = 1;%插值窗口大小cec。参考=dmr的;%使用dmr作为参考信号

displaySimulationInformation = true;

%初始化变量用于模拟和分析%数组来存储所有信噪比点的最大吞吐量maxThroughput = 0(长度(SNRIn), 1);%数组来存储所有信噪比点的模拟吞吐量simThroughput = 0(长度(SNRIn), 1);%数组来存储房车序列指数的历史allRvSeqIdxHistory =细胞(1,元素个数(SNRIn));%的渠道结构优化并行处理(只有%运行并行计算工具箱的例子)channelInit =通道;%在enb模拟一些字段将被更新,复制%的时候重新启动模拟每一个信噪比点enbInit = enb;% TDD precalculate向量的子帧类型:D, S和U下行,%分别特殊和上行如果strcmpi (DuplexMode“TDD”1)subframeType = char(10日);initialSubframeNo = enb.NSubframe;为sNo = 0:9%为所有子帧的帧enb。NSubframe = sNo;duplexInfo = lteDuplexingInfo (enb);subframeType (sNo + 1) = duplexInfo.SubframeType (1);%的第一个字符:D S或U结束enb。NSubframe = initialSubframeNo;结束% CFI可以是一个标量或矢量,创建一个本地副本CFI的向量%(每个子帧一个值)如果元素个数(enb.CFI) = = 1 CFI = repmat (enb.CFI 1 10);其他的CFI = enb.CFI;结束%得到HARQ的HARQ ID序列处理。这是一个索引的列表% HARQ进程调度。[~,~,enbOut] = lteRMCDLTool (enb, []);harqProcessSequenceInit = enbOut.PDSCH.HARQProcessSequence;为snrIdx = 1:元素个数(SNRIn)%的评论进行并行计算% parfor snrIdx = 1:元素个数为并行计算(SNRIn) %取消%设置随机数生成器的种子取决于循环变量%,以确保独立随机流rng (snrIdx“combRecursive”);%初始化变量(他们在每一个信噪比点被修改%模拟)enb = enbInit;频道= channelInit;harqProcessSequence = harqProcessSequenceInit;%初始化所有HARQ进程的状态harqProcesses = hNewHARQProcess (enb);%设置变量的主循环lastOffset = 7;%初始化总体框架定时偏移量%(设置等于信道实现延迟)frameOffset = 7;%初始化帧定时偏移量%(设置等于信道实现延迟)blkCRC = [];%块CRC认为子帧bitTput = [];%成功收到位/子帧的数量txedTrBlkSizes = [];%每个子帧的传输比特数precodingMatrix = [];%预编码矩阵rxSymbols = [];星座图% DL-SCH符号(码字1)rxSymbols2 = [];星座图% DL-SCH符号(码字2)%变量rvSeqIdxHistory商店的历史价值%房车指数序列的HARQ进程rvSeqIdxHistory =南(ncw, NFrames * 10);%国旗表明如果一个预编码矩阵W是可用的。这是非常有用的%如果没有W计算仿真的开始。isWready = false;%标志表明是子帧处理。如果设置为真%有数据处理子帧,即DL子帧或%零传输块大小。processSubframe = false;%主要for循环:子帧为subframeNo = 0: (NFrames * 10 - 1)%更新子帧数enb。NSubframe = subframeNo;%负荷CFI当前子帧enb。CFI = CFI (mod (subframeNo长度(CFI)) + 1);%通道当前子帧的时间通道。InitTime = subframeNo / 1000;%从HARQ过程得到HARQ进程ID为子帧序列harqID = harqProcessSequence (mod (subframeNo长度(harqProcessSequence)) + 1);%提取当前子帧传输块大小(s)trBlk = trBlkSizes(:,国防部(subframeNo, 10) + 1)”;%在之前的子帧,如果加载预编码矩阵计算%的存在,如果有数据传输。设置标志来触发%子帧处理。如果isWready & &任何(trBlk) harqProcesses .txConfig (harqID)。W = precodingMatrix;processSubframe = true;其他的processSubframe = false;结束%预编码矩阵计算如果strcmpi (DuplexMode“TDD”)%估计通道在UL子帧如果比较字符串(subframeType (mod (subframeNo, 10) + 1),“U”)processSubframe = false;% UL子帧,没有DL数据%只有执行信道估计下一子帧是DL如果比较字符串(subframeType (mod ((subframeNo + 1), 10) + 1),' D ')precodingMatrix = hCalculatePrecodingMatrix (enb、通道);isWready = true;结束结束其他的% FDD%获得下一子帧传输块trBlkNext = trBlkSizes(:,国防部(subframeNo + 1, 10) + 1)”;%计算下的预编码矩阵子帧只有它%进行数据(即非零trBlkNext)如果任何(trBlkNext) precodingMatrix = hCalculatePrecodingMatrix (enb、通道);isWready = true;其他的isWready = false;结束结束%子帧处理如果processSubframe%更新当前HARQ进程harqProcesses (harqID) = hHARQScheduling (…harqProcesses (harqID)、subframeNo rvSequence);%显示运行时信息如果任何(trBlk) disp displaySimulationInformation & & (' ');disp ([的子帧:num2str (subframeNo)”。HARQ过程指数:num2str (harqID)]);结束%更新房车序列索引表rvSeqIdxHistory (:, subframeNo + 1) =…harqProcesses harqID .txConfig.RVIdx。';%更新PDSCH传输配置HARQ%流程状态enb。PDSCH = harqProcesses (harqID) .txConfig;%细胞有效载荷dlschTransportBlk = harqProcesses (harqID) . data;%建立传输波形txWaveform = lteRMCDLTool (enb dlschTransportBlk);%添加maxChDelay样本填充。这是延迟的范围%预期从信道建模(的组合%实现延迟和延迟扩展频道)txWaveform = [txWaveform;0 (maxChDelay ntxants)];% #好< AGROW >%通过通道模型传递数据rxNoiselessWaveform = lteFadingChannel(通道、txWaveform);%计算噪声增益包括下行补偿%功率分配信噪比= 10 ^ (SNRIn (snrIdx) -enb.PDSCH.Rho) / 20);%正常化噪声功率的采样率,%是传输线大小的函数用于OFDM%调制,天线的数量N0 = 1 / (√(2.0 * ntxants *双(ofdmInfo.Nfft)) *信噪比);%建立加性高斯白噪声噪音= N0 *复杂(randn(大小(rxNoiselessWaveform)),…randn(大小(rxNoiselessWaveform)));%增加情况下接收到的时域波形和规模%的需用功率rxWaveform = rxNoiselessWaveform +噪声;%接收%每帧,在子帧0,计算一个新的%同步偏移量。一个偏移量的范围之内%延误预计从信道建模(组合%的实现延迟和延迟扩展频道)%表示成功如果(国防部(subframeNo 10) = = 0) frameOffset = lteDLFrameOffset (enb rxWaveform);如果(frameOffset > maxChDelay) frameOffset = lastOffset;结束lastOffset = frameOffset;结束rxWaveform = rxWaveform (1 + frameOffset:最终,);%完成OFDM解调接收的数据获取%的资源网格rxSubframe = lteOFDMDemodulate (enb rxWaveform);%信道估计[estChannelGrid,噪音]=…lteDLChannelEstimate (enb.PDSCH, enb cec rxSubframe);%进行均衡,deprecoding demapping层,%解调和descrambling接收的数据使用%信道估计。%得到PDSCH指数pdschIndices = ltePDSCHIndices (enb enb.PDSCH enb.PDSCH.PRBSet);%得到PDSCH资源元素。规模接收到的子帧% PDSCH功率因数ρ。PDSCH是按比例缩小的%ρ,而细胞参考符号用于通道%估计(用于PDSCH解码阶段)。[pdschRx, pdschHest] = lteExtractResources (pdschIndices,…rxSubframe * (10 ^ (-enb.PDSCH.Rho / 20)), estChannelGrid);%解码PDSCH[dlschBits, dlschSymbols] = ltePDSCHDecode (enb enb.PDSCH,…pdschRx pdschHest,噪音);%存储解码DLSCH星座符号%绘制rxSymbols = [rxSymbols;dlschSymbols {1} (:));% #好< AGROW >如果作战飞机> 1 rxSymbols2 = [rxSymbols2;dlschSymbols {2} (:));% #好< AGROW >结束%解码DL-SCH[decbits harqProcesses (harqID) .blkerr harqProcesses (harqID) .decState] =…lteDLSCHDecode (enb.PDSCH, enb 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 >结束结束结束%的阴谋得到星座象征如果displaySimulationInformation图;情节(rxSymbols,“r”);标题([的星座。信噪比的…num2str (SNRIn (snrIdx))dB。码字1 '])包含(“同相振幅”);ylabel (“正交振幅”);如果作战飞机> 1图;情节(rxSymbols2“r”);标题([的星座。信噪比的…num2str (SNRIn (snrIdx))dB。码字2》])包含(“同相振幅”);ylabel (“正交振幅”);结束结束%计算和模拟的最大吞吐量maxThroughput (snrIdx) =总和(txedTrBlkSizes);%最大可能的吞吐量simThroughput (snrIdx) = (bitTput, 2)总和;%模拟吞吐量%动态在命令窗口中显示结果流(' \ nSNR = %。2 f dB。吞吐量(s) = % % d帧。4 f Mbps \ n”,…SNRIn (snrIdx) NFrames 1 e-6 * simThroughput (snrIdx) / (NFrames * 10 e - 3));流('信噪比= %。2 f dB。吞吐量(% %)% d帧(s) = %。4 f % % \ n”,…SNRIn (snrIdx)、NFrames simThroughput (snrIdx) * 100 / maxThroughput (snrIdx));allRvSeqIdxHistory {snrIdx} = rvSeqIdxHistory;结束

hPlotRVSequence (SNRIn allRvSeqIdxHistory NFrames);

%绘制吞吐量图legendString = (char (txMode)“:”num2str (simulationParameters.PDSCH.NLayers)…层(s),num2str (ntxants)“蚁(s)”];情节(SNRIn simThroughput * 100. / maxThroughput,* - - - - - -。);包含(“信噪比(dB)”);ylabel (的吞吐量(%));标题(吞吐量和信噪比的);传奇(legendString“位置”,“西北”);网格在;