numFrames = 1;% 10 ms的框架信噪比= 20;在dB %信噪比

%创建问题/载波配置问题= nrCarrierConfig;问题。NSizeGrid = 52个;%的带宽数量的资源块(52苏格兰皇家银行在15 khz SCS 10 mhz BW)问题。SubcarrierSpacing = 15;% 15、30、60、120、240(赫兹)问题。CyclicPrefix =“正常”;%“正常”或“扩展”nTxAnts = 2;%的发送天线数量(1、2、4)nRxAnts = 2;%的接收天线nLayers = min (nTxAnts nRxAnts);%配置一个周期多口SRS和支持跳频srs = nrSRSConfig;srs。NumSRSSymbols = 4;%的OFDM符号分配每槽(1、2、4)srs。SymbolStart = 8;%开始OFDM符号在一个插槽srs。NumSRSPorts = nTxAnts;%的SRS天线端口(1、2、4)。srs。FrequencyStart = 0;%的频率位置SRS BWP苏格兰皇家银行srs。NRRC = 0;%额外抵消FreqStart 4块规定的伪随机位序列(0。67)srs。csr= 14;%的带宽配置C_SRS (0。63)。它控制分配带宽SRSsrs。建筑= 0;%的带宽配置B_SRS (0。3)。它控制分配带宽SRSsrs。BHop= 0;%跳频配置(0。3)。设置BHop <建筑实现跳频srs。KTC = 2;%梳数(2、4)。副载波频率密度srs。重复= 2;%(1、2、4)重复。它禁用跳频块| |符号重复srs。SRSPeriod = [2 0];在槽%周期性和偏移量。必须< SRSPeriod SRSPeriod (2) (1)srs。ResourceType =“周期”;%的资源类型(“周期”、“semi-persistent”,“非周期”)。使用“非周期”禁用inter-slot跳频

practicalSynchronization = true;

csiSubbandSize = 4;%每个部分波段的苏格兰皇家银行的数量

频道= nrTDLChannel;通道。DelayProfile =“TDL-C”;通道。DelaySpread = 40 e-9;通道。MaximumDopplerShift = 30;通道。NumTransmitAntennas = nTxAnts;通道。NumReceiveAntennas = nRxAnts;通道。Seed = 5;%设置通道采样率ofdmInfo = nrOFDMInfo(问题);通道。SampleRate = ofdmInfo.SampleRate;

chInfo =信息(渠道);maxChDelay =装天花板(max (chInfo.PathDelays * channel.SampleRate));maxChDelay = maxChDelay + chInfo.ChannelFilterDelay;%重置随机发生器再现性rng (“默认”);

%槽数来模拟numSlots = numFrames * ue.SlotsPerFrame;%的总数每槽副载波和符号K =问题。NSizeGrid * 12;L = ue.SymbolsPerSlot;%初始化数组存储信道估计allTxGrid = 0 ([K L * numSlots nTxAnts]);slotGridSize = [K L nRxAnts nTxAnts];命令= 0 (slotGridSize);hestInterp = 0 (slotGridSize);hEstUpdate = 0 (slotGridSize);totalGridSize = [K L * numSlots nRxAnts nTxAnts);allHest = 0 (totalGridSize);allHestPerfect = 0 (totalGridSize);allHestInterp = 0 (totalGridSize);%初始化噪声功率估计据nvar = 0;%计算承运人CSI部分波段的数目numCSISubbands =装天花板(ue.NSizeGrid / csiSubbandSize);%初始化每部分波段SINR,槽,PMImaxPMI = hMaxPUSCHPrecodingMatrixIndicator (nLayers nTxAnts);sinrSubband = 0 ([numCSISubbands numSlots maxPMI + 1]);%初始化采购经理人指数矩阵和SINR的损失采购经理人指数=南(numCSISubbands numSlots);pmiPerfect =采购经理人指数;损失= 0(大小(pmi));%初始化时间估计抵消抵消= chInfo.ChannelFilterDelay;%计算SRS CDM的长度cdmLengths = hSRSCDMLengths (srs);% OFDM符号用于CSI收购csiSelectSymbols = srs。SymbolStart + (1: srs.NumSRSSymbols);为nSlot = 0: numSlots-1%更新位置计数器问题。NSlot = NSlot;%生成SRS和映射到网格[srsIndices, srsIndInfo] = nrSRSIndices(问题,srs);srsSymbols = nrSRS(问题,srs);%创建一个SRS slot-wise资源网格空网格和地图符号nTxAnts txGrid = nrResourceGrid(问题);txGrid (srsIndices) = srsSymbols;%确定槽包含SRSisSRSSlot = ~ isempty (srsSymbols);% OFDM调制[txWaveform, waveformInfo] = nrOFDMModulate(问题,txGrid);txWaveform = [txWaveform;0 (maxChDelay、尺寸(txWaveform 2)));%要求后冲洗通道滤波器得到接收信号%通过信道传输[rxWaveform, pathGains] =通道(txWaveform);%添加情况下接收到的时域波形%正常化噪声功率的采样率,这是%传输线大小的函数用于OFDM调制。的信噪比% /重新定义为每个接收天线(TS 38.101 - 4)。信噪比= 10 ^(信噪比/ 10);N0 = 1 /√(2.0 * nRxAnts *双(waveformInfo.Nfft) *信噪比);噪音= N0 *复杂(randn(大小(rxWaveform)), randn(大小(rxWaveform)));rxWaveform = rxWaveform +噪声;%执行定时偏移估计pathFilters = getPathFilters(渠道);%的时间估计是只执行SRS的插槽中%传播保持完美的和实际的信道估计%同步。如果isSRSSlot如果practicalSynchronization%实际同步。将接收到的波形%的SRS给定时偏移估计抵消= nrTimingEstimate(问题、rxWaveform srsIndices, srsSymbols);其他的抵消= nrPerfectTimingEstimate (pathGains pathFilters);结束结束%完成OFDM解调rxGrid = nrOFDMDemodulate(问题,rxWaveform(1 +抵消:最终,));%执行实际的信道估计%更新信道估计仅在符号包含SRS和苏格兰皇家银行%在这个位置,前面的位置的估计%的位置。据nvar不是当没有更新SRS传播。使用%的时间平均窗口覆盖所有的SRS符号传播。如果isSRSSlot%这个位置包含一个SRS传播[命令,据nvar] = nrChannelEstimate(问题、rxGrid srsIndices, srsSymbols,“AveragingWindow”[0 7],“CDMLengths”,cdmLengths);%使用从先前的槽的OFDM信道估计符号之前第一个SRS的象征hestInterp = repmat (hestInterp (:,:,:), 1, ue.SymbolsPerSlot);%更新信道估计在OFDM符号和RB SRS%现在和保存所有信道估计在年底前槽firstSymbol = srs.SymbolStart + 1;lastSymbol = srs。SymbolStart + srs.NumSRSSymbols;firstSymbol hEstUpdate (:: lastSymbol、::) =命令(:,firstSymbol: lastSymbol,:,:);lastSymbol hEstUpdate (::,:,:) = repmat(命令(lastSymbol,::,:), 1, ue.SymbolsPerSlot-lastSymbol + 1);idxHEstUpdate = hEstUpdate ~ = 0;%更新信道估计的指标hestInterp (idxHEstUpdate) = hEstUpdate (idxHEstUpdate);其他的%保存以前的信道估计,如果这个位置不包含SRShestInterp = repmat (hestInterp (:,:,:), 1, ue.SymbolsPerSlot);结束% PMI选择%选择预编码器矩阵指标层使用%的插值信道估计OFDM符号规定% csiSelectSymbols。采购经理人指数的估计每CSI部分波段pmi (:, nSlot + 1), sinrSubband (nSlot + 1,::), subbandIndices] = hPMISelect (nLayers, hestInterp (csiSelectSymbols,::,:),据nvar, csiSubbandSize);% PMI选择SINR的损失%计算性能损失的比SINR预编码%与采购经理人指数选择使用一个实际的信道估计和SINR%与采购经理人指数预编码后选择使用一个完美的信道估计。%计算完美的PMI选择完美的信道估计hEstPerfect = nrPerfectChannelEstimate (pathFilters pathGains退居其次,抵消);%完美的噪声估计噪声的实现noiseGrid = nrOFDMDemodulate(问题、噪声(1 +抵消:最终,));nvarPerfect = var (noiseGrid (:));[损失:nSlot + 1), pmiPerfect (:, nSlot + 1)] = hPMISelectionSINRLoss (pmi (:, nSlot + 1), nLayers, hEstPerfect csiSelectSymbols,::,:, nvarPerfect);%保存一个副本的OFDM传输网格和信道估计%显示目的thisSlot = nSlot * L + (1: L);%符号的当前位置thisSlot, allTxGrid (::) = txGrid;allHest (thisSlot,::,:) =命令;allHestInterp (thisSlot,::,:) = hestInterp;allHestPerfect (thisSlot,::,:) = hEstPerfect;结束

%建立轴和轴向量符号= 0:(ue.NSlot + 1) * ue.SymbolsPerSlot-1;槽= 0:ue.NSlot;副载波= 1:ue.NSizeGrid * 12;resourceBlocks = 1: ue.NSizeGrid;

图显示亮度图像(符号、副载波、abs (allTxGrid (:,:, 1, 1)));包含(OFDM符号的);ylabel (副载波的);轴xy;标题(“传播SRS(端口1)”);

%删除第一个OFDM符号不包含SRS改进的可视化%信道估计误差hEstInterp = allHestInterp;idx = 1:(srs.SRSPeriod(2) *问题。SymbolsPerSlot + srs.SymbolStart);hEstInterp (idx,::,:) =南;hEstInterp ((srs.NRB * 12 + 1):,:,:,:) =南;hChannelEstimationPlot(符号、副载波allHestPerfect hEstInterp);

%首次扩大亏损部分波段为苏格兰皇家银行用于显示目的pmiRB = hExpandSubbandToRB (pmi csiSubbandSize ue.NSizeGrid);pmiPerfectRB = hExpandSubbandToRB (pmiPerfect csiSubbandSize ue.NSizeGrid);lossRB = hExpandSubbandToRB(损失,csiSubbandSize ue.NSizeGrid);hPMIPlot(槽resourceBlocks pmiRB、pmiPerfectRB lossRB);

numLayers = min(大小(allHestInterp [3 - 4]));如果numLayers ~ = 1 pmiErr =总和(abs pmiPerfect (pmi) > 0, [1 - 2])。/笔(~ isnan (pmi), [1 - 2]);TotPMIEst =总和(~ isnan (pmi), [1 - 2]);流(的数量估计PMI: % d \ n ',TotPMIEst);流(数量错了PMI: % d \ n '装天花板(pmiErr * TotPMIEst));流(相对误差:%。1 f (% %) \ n ',pmiErr * 100);其他的流(的一层,PMI总是0。\ n”);结束

hBestSINRPlot(槽resourceBlocks sinrSubband pmi, csiSubbandSize);

%显示完美的和实际的信道估计和信道%估计误差第一传输和接收端口。英吉利海峡%估计误差被定义为不同的绝对值%之间完美的和实际的信道估计。函数hChannelEstimationPlot(符号、副载波allHestPerfect allHestInterp)图次要情节(311)显示亮度图像(符号、副载波、abs (allHestPerfect (:,:, 1, 1)));轴xy;包含(OFDM符号的);ylabel (副载波的);colorbar;标题(“完美信道估计(TxAnt = 1, RxAnt = 1)的);次要情节(312)显示亮度图像(符号、副载波、abs (allHestInterp (:,:, 1, 1)),…“AlphaData”、~ isnan (allHestInterp(:,:, 1, 1)))轴xy;包含(OFDM符号的);ylabel (副载波的);colorbar;标题(“SRS-based实际信道估计(TxAnt = 1, RxAnt = 1)”);%显示信道估计错误,定义为之间的差异% SRS-based和完美的信道估计次要情节(313)hestErr = abs (allHestInterp - allHestPerfect);显示亮度图像(符号、副载波hestErr (:,:, 1, 1),…“AlphaData”、~ isnan (hestErr (:,:, 1, 1)));轴xy;包含(OFDM符号的);ylabel (副载波的);colorbar;标题(的信道估计误差(TxAnt = 1, RxAnt = 1)的);结束%显示PMI进化和PMI估计随着时间的推移和SINR损失%的频率。的SINR损失被定义为一个比SINR%预编码估计和完美的采购经理人指数。预计采购经理人指数%获得使用一个实际的信道估计和完美的采购经理人指数%选择使用一个完美的信道估计。函数hPMIPlot(槽resourceBlocks pmiRB、pmiPerfectRB lossRB)图次要情节(311)显示亮度图像(插槽,resourceBlocks pmiPerfectRB,“AlphaData”~ isnan (pmiPerfectRB));轴xy;c = caxis;厘米= colormap;colormap(厘米(1:地板(尺寸(厘米,1)/ (c (2) - c (1)) 1):,:));%调整colormap PMI离散值colorbar包含(“槽”);ylabel (“资源块”)、标题(PMI使用完美的信道估计来选择的)次要情节(312)显示亮度图像(插槽,resourceBlocks pmiRB,“AlphaData”~ isnan (pmiRB));轴xy;colorbar,包含(“槽”);ylabel (“资源块”)、标题(PMI使用SRS来选择的)次要情节(313)显示亮度图像(插槽,resourceBlocks lossRB,“AlphaData”~ isnan (lossRB));甘氨胆酸colormap(厘米)包含(“槽”);ylabel (“资源块”);轴xy;colorbar;标题(“PMI估计SINR损失(dB)”)结束%显示每个资源块的SINR的预编码后得到% PMI最大化每个部分波段SINR。函数hBestSINRPlot(槽resourceBlocks sinrSubband pmi, csiBandSize)%显示SINR预编码后最好的采购经理人指数bestSINRPerSubband =南(大小(sinrSubband [1 - 2]));%得到SINR每部分波段和槽使用最好的采购经理人指数(某人,nslot) =找到(~ isnan (pmi));为i = 1:长度(某人)bestSINRPerSubband(某人(我),nslot(我))= sinrSubband(某人(我),nslot(我),pmi(某人(我),nslot(我))+ 1);结束%的第一个扩展SINR从次能带到苏格兰皇家银行用于显示目的bestSINRPerRB = hExpandSubbandToRB (bestSINRPerSubband csiBandSize、长度(resourceBlocks));图sinrdb = 10 * log10 (abs (bestSINRPerRB));显示亮度图像(插槽,resourceBlocks sinrdb,“AlphaData”~ isnan (sinrdb));轴xy;colorbar;包含(“槽”);ylabel (“资源块”)标题(每部分波段的平均SINR,槽与最好的PMI (dB)预编码后的)结束%每部分波段扩展一个二维矩阵的值在第一维度成%每个资源块矩阵的值。函数rbValues = hExpandSubbandToRB (subbandValues bandSize NRB) lastBandSize =国防部(NRB bandSize);lastBandSize = lastBandSize + bandSize * (lastBandSize = = 0);rbValues = [kron (subbandValues (1: end-1,:) (bandSize, 1));…subbandValues(最终:)。*的(lastBandSize, 1)];结束