NFrames = 4;帧数%eqMethod =“MMSE_IRC”；%的患者,MMSE_IRC

SINR = 0.8;% SINR in dBDIP2 = -1.73;细胞2的dB的% DIPDIP3 = -8.66;细胞3的dB的% DIPNoc = -98;% dBm/15kHz平均功率谱密度

设置随机数生成器种子rng (“默认”）;根据R.47设置单元1 eNodeB配置simulationParameters =结构;simulationParameters。NDLRB = 50;simulationParameters。CellRefP = 2;simulationParameters。NCellID = 0;simulationParameters。CFI = 2;simulationParameters。DuplexMode =“FDD”；simulationParameters。TotSubframes = 1;这不是总数%子帧在模拟中使用，只是我们的子帧的数量对波形发生器的每次调用生成%。指定PDSCH配置子结构simulationParameters.PDSCH.TxScheme =“SpatialMux”；simulationParameters.PDSCH.Modulation = {16 qam的};simulationParameters.PDSCH.NLayers = 1;simulationParameters.PDSCH.Rho = 3;simulationParameters.PDSCH.PRBSet = (0:49) ';表A.3.3.2.1-2, TS 36.101simulationParameters.PDSCH.TrBlkSizes = [8760 8760 8760 8760 8760 8760 0 8760 8760 8760 8760 8760];表A.3.3.2.1-2, TS 36.101simulationparameters . pdch . codedtrblksizes = [24768 26400 26400 26400 26400 26400 0 26400 26400 26400];simulationParameters.PDSCH.CSIMode =“PUCCH 1 - 1”；simulationParameters.PDSCH.PMIMode =“宽带”；simulationParameters.PDSCH.CSI =“上”；simulationParameters.PDSCH.W = [];simulationParameters.PDSCH.CodebookSubset =“1111”；指定PDSCH OCNG配置simulationParameters。OCNGPDSCHEnable =“上”；%启用OCNG填充simulationParameters。OCNGPDSCHPower = 3;% OCNG功率与PDSCH Rho相同simulationParameters.OCNGPDSCH.RNTI = 0;虚拟终端simulationParameters.OCNGPDSCH.Modulation =“正交相移编码”；OCNG符号调制simulationParameters.OCNGPDSCH.TxScheme =“TxDiversity”；% OCNG传输模式

enb1 = lteRMCDL (simulationParameters);

%细胞2enb2 = enb1;enb2。NCellID = 1;enb2。OCNGPDSCHEnable =“关闭”；%细胞3enb3 = enb1;enb3。NCellID = 2;enb3。OCNGPDSCHEnable =“关闭”；

eNodeB1到UE传播通道channel1 =结构;通道配置结构。种子= 20;%通道种子。NRxAnts = 2;% 2接收天线。DelayProfile =“爱娃”；%延迟概要。DopplerFreq = 5;%多普勒频率。MIMOCorrelation =“低”；% Multi-antenna相关性。NTerms = 16;在衰落模型中使用的振荡器。ModelType =“GMEDS”；瑞利衰落模型类型。InitPhase =“随机”；%随机初相。没有rmalizePathGains =“上”；规范化延迟配置文件功率。没有rmalizeTxAnts =“上”；%发射天线归一化

ofdmInfo = lteOFDMInfo (enb1);。SamplingRate = ofdmInfo.SamplingRate;eNodeB2(干扰)到UE传播通道channel2 = channel1;channel2。种子= 122;%通道种子eNodeB3(干扰)到UE传播通道channel3 = channel1;channel3。种子= 36;%通道种子

通道估计器行为perfectChanEstimator = true;通道估计器配置结构定义如下。的%的频率和时间平均窗口大小被选择来跨越一个相对%大量的资源元素。选择大窗口尺寸%平均尽可能多的噪音和干扰的资源%的元素。注意，时间和/或频率的平均窗口太大%将导致信息丢失，因为平均出通道变化。这会产生越来越不完美的信道估计，这可能会影响均衡器的性能。。cec =结构;通道估计配置结构cec。PilotAverage =“UserDefined”；导频符号平均的类型cec。FreqWindow = 31;%频率窗口大小cec。TimeWindow = 23;%时间窗口大小cec。InterpType =“立方”；% 2D插值类型cec。InterpWindow =“中心”；插值窗口类型cec。InterpWinSize = 1;插值窗口大小

%通道噪声设置nocLin = 10 ^ (Noc / 10) * (1 e - 3);瓦茨线性百分比%考虑FFT (OFDM)缩放没有=√nocLin /(2 *双(ofdmInfo.Nfft)));%信号和干扰幅值的比例因子计算[K1, K2, K3] = hENBscalingFactors(DIP2, DIP3, Noc, SINR, enb1, enb2, enb3);

所有HARQ进程的初始化状态harqProcesses = hNewHARQProcess (enb1);初始化HARQ进程id为1作为第一个非零传输%块将始终使用第一个HARQ进程传输。这%将在之后用lteRMCDLTool输出的完整序列更新%对函数的第一次调用harqProcessSequence = 1;为主循环设置变量lastOffset = 0;初始化上一帧的帧时间偏移量frameOffset = 0;初始化帧定时偏移量nPMI = 0;初始化预编码器矩阵的数量指示% (PMI)已计算集blkCRC = [];为所有考虑的子帧阻塞CRCbitTput = [];每子帧成功接收的比特数txedTrBlkSizes = [];每子帧传输的比特数%传输的总比特数的向量，为每个子帧计算。runningMaxThPut = [];用于计算的，存储成功接收比特数的向量每个子帧%。runningSimThPut = [];%获取发射天线数。= lteDLResourceGridSize dim (enb1);P = dim (3);为每个码字和传输指定冗余版本序列。每个子帧的%块大小rvSequence = enb1.PDSCH.RVSeq;trBlkSizes = enb1.PDSCH.TrBlkSizes;设置闭环空间复用的PMI延迟时间pmiDelay = 8;如表8.2.1.4.1B-1中TS 36.101所示初始化第一个'pmiDelay'子帧的pmipmiDims = ltePMIInfo (enb1 enb1.PDSCH);(pmiDims txPMIs = 0。NSubbands pmiDelay);%标志，指示是否有有效的PMI反馈从UE可用pmiReady = false;

流(“\ nSimulating % d帧(s) \ n”, NFrames);% Main for循环:用于所有子帧为subframeNo = 0: (NFrames * 10 - 1)为每个子帧重新初始化通道种子以增加可变性。种子= 1 + subframeNo;channel2。种子= 1 + subframeNo + (NFrames * 10);channel3。种子= 1 + subframeNo + 2 * (NFrames * 10);更新子帧号enb1。NSubframe = subframeNo;enb2。NSubframe = subframeNo;enb3。NSubframe = subframeNo;duplexInfo = lteDuplexingInfo (enb1);如果duplexInfo。NSymbolsDL ~ = 0%目标仅下行子帧从HARQ进程序列中获取子帧的HARQ进程IDharqID = harqProcessSequence(mod(subframe, length(harqProcessSequence))+1);%如果在当前子帧中有一个调度的传输块%(由非零'harqID'表示)，执行传输和%接待。否则，继续到下一个子帧。如果harqID = = 0继续；结束更新当前HARQ进程harqProcesses (harqID) = hHARQScheduling (.．.harqProcesses (harqID)、subframeNo rvSequence);提取当前子帧传输块大小(s)trBlk = trBlkSizes(:， mod(subframino, 10)+1).';使用HARQ进程状态更新PDSCH传输配置enb1.PDSCH。RVSeq = harqProcesses (harqID) .txConfig.RVSeq;enb1.PDSCH。房车= harqProcesses (harqID) .txConfig.RV;dlschTransportBlk = harqProcesses (harqID) . data;%将PMI设置为延迟队列中的适当值如果strcmpi (enb1.PDSCH.TxScheme“SpatialMux”) pmiIdx = mod(subframeNo, pmiDelay);延迟队列中的PMI索引enb1.PDSCH。PMISet = txPMIs(:， pmiIdx+1);%设置采购经理人指数结束%创建发送波形[tx，~，enbOut] = lteRMCDLTool(enb1, dlschTransportBlk);Pad 25个样品，以覆盖来自通道的预期延迟范围。%建模(实现延迟和通道延迟的组合%传播)txWaveform1 = [tx;0(25日P)];从' enabout '获取HARQ ID序列用于HARQ处理harqProcessSequence = enbOut.PDSCH.HARQProcessSequence;根据TS 36.101, B.5.3生成干扰模型。的htm4intermodel生成干扰发射信号。txWaveform2 = [hTM4InterfModel (enb2);0(25日P)];txWaveform3 = [hTM4InterfModel (enb3);0(25日P)];指定当前子帧的通道时间。InitTime = subframeNo / 1000;channel2。InitTime = channel1.InitTime;channel3。InitTime = channel1.InitTime;通过通道传递数据rxWaveform1 = lteFadingChannel (channel1 txWaveform1);rxWaveform2 = lteFadingChannel (channel2 txWaveform2);rxWaveform3 = lteFadingChannel (channel3 txWaveform3);%产生噪音噪音=没有*复杂(randn(大小(rxWaveform1)),.．.randn(大小(rxWaveform1)));在接收的时域波形中添加AWGNrx波形= K1*rxWaveform1 + K2*rxWaveform2 + K3*rxWaveform3 + noise;%接收每一帧，在子帧0上计算一个新的同步%抵消如果(mod(subframeamo,10) == 0) frameOffset = lteDLFrameOffset(enb1, rx波形);如果(frameOffset = lastOffset;结束lastOffset = frameOffset;结束同步接收到的波形rx波形= rx波形(1+frameOffset:end，:);%缩放rx波形1/K1，以避免数值问题%信道解码级rxWaveform = rxWaveform / K1;对接收到的数据进行OFDM解调%资源网格rxSubframe = lteOFDMDemodulate(enb1, rx波形);%执行通道估计如果(perfectChanEstimator) estChannelGrid = lteDLPerfectChannelEstimate(enb1, channel1, frameOffset);noiseInterf = K2*rxWaveform2 + K3*rxWaveform3 + noise;noiseInterf = noiseInterf / K1;noiseGrid = lteOFDMDemodulate(enb1, noiseInterf(1+frameOffset:end，:));噪音= var (noiseGrid (:));其他的[estChannelGrid, noiseEst] = lteDLChannelEstimate(.．.enb1 enb1。PDSCH,cec, rxSubframe);结束%获取PDSCH索引pdschIndices = ltePDSCHIndices (enb1 enb1.PDSCH enb1.PDSCH.PRBSet);获取PDSCH资源元素。缩放接收到的子帧PDSCH功率因数。[pdschRx, pdschest] = lteExtractResources(pdschIndices，.．.rxSubframe * (10 ^ (-enb1.PDSCH.Rho / 20)), estChannelGrid);执行均衡和解预编码如果比较字符串(eqMethod“患者”）% MIMO均衡和预编码(基于MMSE)[rxDeprecoded, csi] = lteEqualizeMIMO (enb1 enb1。PDSCH,.．.pdschRx pdschHest,噪音);其他的% MIMO均衡和预编码(基于MMSE-IRC)[rxDeprecoded, csi] = hEqualizeMMSEIRC (enb1 enb1。PDSCH,.．.rxSubframe estChannelGrid,噪音);结束执行层解映射，解调和解乱水煤浆= ltePDSCHDecode (enb1 setfield (enb1。PDSCH,“TxScheme”，.．.“Port7-14”), rxDeprecoded);PDSCH传输方案为%修改为端口7-14，以便跳过解压缩操作%根据CSI缩放llr水煤浆= hCSIscaling (enb1.PDSCH,水煤浆,csi);%解码DL-SCH[decbits harqProcesses (harqID) .blkerr harqProcesses (harqID) .decState] =.．.lteDLSCHDecode (enb1 enb1。PDSCH trBlk,水煤浆,.．.harqProcesses (harqID) .decState);用于计算吞吐量的存储值%仅用于有数据和有效PMI反馈的子帧如果any(trBlk) && pmiReady blkCRC = [blkCRC harqProcesses(harqID).blkerr];txedTrBlkSizes = [txedTrBlkSizes trBlk];bitTput = [bitTput trBlk.*(1-harqProcesses(harqID).blkerr)];结束runningSimThPut = [runningSimThPut sum(bitTput,2)];runningMaxThPut = [runningMaxThPut sum(txedTrBlkSizes,2)];向eNodeB提供PMI反馈如果strcmpi (enb1.PDSCH.TxScheme“SpatialMux”) PMI = ltePMISelect(enb1, enb1.)PDSCH,estChannelGrid, noiseEst); txPMIs(:, pmiIdx+1) = PMI; nPMI = nPMI+1;如果nPMI>=pmiDelay pmiReady = true;结束结束结束结束

maxThroughput =总和(txedTrBlkSizes);%最大可能吞吐量simThroughput = (bitTput, 2)之和;%模拟吞吐量显示实现吞吐量百分比disp ([“实现吞吐量”num2str (simThroughput * 100 / maxThroughput)“%”])% Plot运行吞吐量人物,情节(runningSimThPut * 100. / runningMaxThPut) ylabel (的吞吐量(%)）;包含(模拟的子帧的）;标题(的吞吐量）;