loadFromFile = 1;%设置为0，表示在本地生成eNodeB输出

如果loadFromFile负载eNodeBOutput.mat% eNodeB输出的负载I/Q捕获eNodeBOutput = double(eNodeBOutput)/32768;%比例样本Sr = 15.36e6;%加载样本的采样率其他的rmc = lteRMCDL(“R.3”）;% #好< UNRCH >rmc。ncelllid = 17;rmc。TotSubframes = 41;rmc.PDSCH.RNTI = 61;% SIB参数rmc.SIB.Enable =“上”;rmc.SIB.DCIFormat =“Format1A”;rmc.SIB.AllocationType = 0;rmc. sig . vrbstart = 8;rmc. sig . vrlength = 8;rmc.SIB.Gap = 0;% SIB数据字段填充随机位，这不是一个有效的SIB%的消息rmc.SIB.Data = randi([0 1]， 1,176,1);[eNodeBOutput，~，info] = lteRMCDLTool(rmc，[1;0;0;1]);sr = info.SamplingRate;%生成样本的采样率结束

%设置一些内务变量:%分隔符用于命令窗口日志记录分离器=修复“- - -”1、50);%的情节如果(~ (“channelFigure”，“var”) || ~isvalid(channelFigure))“可见”，“关闭”）;结束(简介、synchCorrPlot pdcchConstDiagram] =…hSIB1RecoveryExamplePlots (channelFigure sr);% PDSCH evmpdschEVM = com . evm ();pdschEVM。MaximumEVMOutputPort = true;初始单元格搜索的采样率是使用% lteofdmino为6个资源块配置。enb。设置CyclicPrefix%暂时在调用lteOFDMInfo中抑制默认值%警告(不影响采样率)。Enb = struct;% eNodeB配置结构enb。NDLRB = 6;%资源块数ofdmInfo = lteOFDMInfo(setfield)“CyclicPrefix”，“正常”））;% #好< SFLD >如果(isempty (eNodeBOutput))流(\n接收到的信号不能为空。\n）;返回;结束显示接收到的信号频谱流('绘制接收到的信号频谱…'）;简介(awgn (eNodeBOutput, 100.0));如果(sr ~ = ofdmInfo.SamplingRate)如果(sr < ofdmInfo.SamplingRate)警告(接收信号采样率(%0.3fMs/s)低于小区搜索/ MIB解码所需的采样率(%0.3fMs/s);小区搜索/ MIB解码可能失败。1、sr / 1 e6 ofdmInfo.SamplingRate / e6);结束流(\n从%0.3fMs/s重新采样到%0.3fMs/s用于小区搜索/ MIB解码…\n1、sr / 1 e6 ofdmInfo.SamplingRate / e6);其他的流(不需要重新采样;接收到的信号达到小区搜索/ MIB解码所需的采样率(%0.3fMs/s)。1、sr / e6);结束下采样接收信号nSamples = ceil(ofdmInfo.SamplingRate/round(sr)*size(eNodeBOutput,1));nRxAnts = size(eNodeBOutput, 2);downsampling = 0 (nSamples, nRxAnts);为i=1:nRxAnts downsampling (:，i) = ressample (eNodeBOutput(:，i)， ofdmInfo。SamplingRate圆(sr));结束

流(\n正在执行单元格搜索…\n）;%设置双工模式和循环前缀长度组合进行搜索;如果%在|enb|中配置了这些参数中的任何一个，则值为%假定是正确的如果(~ isfield (enb“DuplexMode”)) duplexModes = {“TDD”“FDD”}；其他的= {enb.DuplexMode};结束如果(~ isfield (enb“CyclicPrefix”)) cyclicPrefixes = {“正常”“扩展”}；其他的cyclicprefix = {enb.CyclicPrefix};结束%跨双工模式和循环前缀长度执行小区搜索%组合，记录相关性最大的组合;如果如果配置了多个单元搜索，本例将解码第一个%(最强)检测细胞searchalg。MaxCellCount = 1;searchalg。SSSDetection =“PostFFT”;peakMax = -Inf;为duplexMode = duplexModes为cyclicPrefix = cyclicPrefixes enb。DuplexMode = DuplexMode {1};enb。CyclicPrefix = CyclicPrefix {1};[enb。ncelllid, offset, peak] = lteCellSearch(enb, downsampling, searchalg);enb。NCellID = enb.NCellID(1);Offset = Offset (1);Peak =峰值(1);如果(峰值>峰值max) enbMax = enb;offsetMax = offset;peakMax =峰值;结束结束结束使用小区标识、循环前缀长度、双工模式和定时在单元格搜索期间给出最大相关性的%偏移量enb = enbMax;offset = offsetMax;计算三个可能的主单元格的相关性%的身份;建立了细胞身份的相关峰%与其他两个的相关性峰值进行比较%原代细胞鉴定，以建立质量%的相关性。Corr = cell(1,3);idGroup = floor(enbMax.NCellID/3);为I = 0:2。NCellID = idGroup*3 + mod(enbMax. 3)ncelllid + i,3);[~，corr{i+1}] = lteDLFrameOffset(enb, downsampling);Corr {i+1} = sum(Corr {i+1}，2);结束阈值= 1.3 * max([corr{2};corr {3}]);经验得出的%乘数为1.3如果{1} (max (corr) <阈值)警告(同步信号相关性较弱;检测到的细胞身份可能不正确。）;结束返回到最初检测到的单元格标识enb。NCellID = enbMax.NCellID;%绘制PSS/SSS相关性和阈值synchCorrPlot。YLimits = [0 max([corr{1};阈值)* 1.1);{1} synchCorrPlot ([corr阈值*的(大小{1})(corr)]);%执行定时同步流('帧开始的定时偏移量:%d采样\n',抵消);downsampling = downsampling (1+offset:end，:);enb。NSubframe = 0;%显示单元格范围的设置流(“搜索单元格后的单元格范围设置:\n”）;disp (enb);

流(“\n执行频率偏移估计…\n”）;%对于TDD，“TDDConfig”和“SSC”默认为“0”。这些参数不是%在系统中建立，直到SIB1被解码，所以在这个阶段%值为0做出最保守的假设(最少的下行链路%子帧和最短特殊子帧)。如果(strcmpi (enb。DuplexMode,“TDD”enb))。TDDConfig = 0;enb。SSC = 0;结束delta_f = lteFrequencyOffset(enb, downsampling);流(“频率偏移:%0.3fHz\n”, delta_f);downsampling = lteFrequencyCorrect(enb, downsampling, delta_f);

%信道估计器配置cec。PilotAverage =“UserDefined”;%先导平均类型cec。FreqWindow = 13;%频率窗口大小cec。TimeWindow = 9;%时间窗口大小cec。InterpType =“立方”;% 2D插值类型cec。InterpWindow =“中心”;%插值窗口类型cec。InterpWinSize = 1;%插值窗口大小%假设有4个cell特有的参考信号用于初始解码尝试;%确保信道估计可用于所有特定于小区的参考%的信号enb。CellRefP = 4;流(“正在进行OFDM解调…\n\n”）;griddims = lteResourceGridSize(enb);%资源网格维度L = griddims(2);%子帧中OFDM符号的数目% OFDM解调信号rxgrid = lteofdmdemomodulation (enb, downsampling);如果(isempty (rxgrid))流(在时间同步后，信号短于一个子帧，因此不会进行进一步的解调。）;返回;结束%执行信道估计[est, nest] = lteDLChannelEstimate(enb, cec, rxgrid(:，1:L，:));

%解码MIB%从第一个PBCH中提取对应的资源元素(REs)%子帧跨所有接收天线和信道估计流('正在执行MIB解码…\n'）;pbchIndices = ltePBCHIndices(enb);[pbchRx, pbchHest] = lteExtractResources(…pbchIndices, rxgrid(:，1:L，:)， hest(:，1:L，:，:));%解码PBCH[bchBits, pbchSymbols, nfmod4, mib, enb.]CellRefP] = ltePBCHDecode(…enb, pbchRx, pbchhst, nest);%解析MIB位enb = lteMIB(mib, enb);%合并函数ltePBCHDecode的nfmod4值输出，如%从MIB建立的NFrame值是系统帧号% (SFN)模4(在MIB中存储为floor(SFN/4))enb。NFrame = enb.NFrame+nfmod4;显示MIB解码后的单元格宽度设置流(' MIB解码后的小区范围设置:\n'）;disp (enb);如果(enb.CellRefP = = 0)流('MIB解码失败(enb.CellRefP=0).\n\n'）;返回;结束如果(enb.NDLRB = = 0)流('MIB解码失败(enb.NDLRB=0).\n\n'）;返回;结束

流('重新开始接收，现在带宽(NDLRB=%d)已知…\n', enb.NDLRB);现在我们知道了真正的带宽ofdmInfo = lteOFDMInfo(enb);如果(sr ~ = ofdmInfo.SamplingRate)如果(sr < ofdmInfo.SamplingRate)警告(接收信号采样率(%0.3fMs/s)低于NDLRB=%d (%0.3fMs/s)的期望采样率;PDCCH搜索/ SIB1解码可能失败。sr / 1 e6, enb.NDLRB ofdmInfo.SamplingRate / 1 e6);结束流(\n从%0.3fMs/s重新采样到%0.3fMs/s…\n1、sr / 1 e6 ofdmInfo.SamplingRate / e6);其他的流(不需要重新采样;接收到的信号在NDLRB=%d (%0.3fMs/s)的期望采样率下。1、enb.NDLRB sr / e6);结束nSamples = ceil(ofdmInfo.SamplingRate/round(sr)*size(eNodeBOutput,1));resampled = 0 (nSamples, nRxAnts);为i = 1:nRxAnts resampled(:，i) = ressample (eNodeBOutput(:，i)， ofdmInfo。SamplingRate圆(sr));结束执行频率偏移估计和校正流(“\n执行频率偏移估计…\n”）;delta_f = lteFrequencyOffset(enb, ressampling);流(“频率偏移:%0.3fHz\n”, delta_f);resampled = lteFrequencyCorrect(enb, resampled, delta_f);找到帧的开头流(“正在执行时序偏移估计…\n”）;offset = lteDLFrameOffset(enb, resampled);流('帧开始的定时偏移量:%d采样\n',抵消);将信号与帧的起始位置对齐Resampled = Resampled (1+offset:end，:);% OFDM解调流(“正在进行OFDM解调…\n\n”）;rxgrid = lteofdm解调(enb，重采样);

检查此帧是否包含SIB1，如果没有提前1帧%有足够的数据，否则终止。如果(mod (enb.NFrame, 2) ~ = 0)如果(大小(rxgrid, 2) > = (L * 10)) rxgrid (:, 1: (L * 10 ),:) = [];流('跳过帧%d(奇数帧号不包含SIB1)。\n\n', enb.NFrame);其他的Rxgrid = [];结束enb。NFrame = enb。NFrame + 1;结束前进到子帧5，如果子帧数少于5则终止如果(size(rxgrid,2)>=(L*5)) rxgrid(:，1):(L*5)，:) = [];%删除子帧0到4其他的Rxgrid = [];结束enb。NSubframe = 5;如果(isempty (rxgrid))流('接收到的信号不包含携带SIB1的子帧。\n\n'）;结束%重置HARQ缓冲区decState = [];当我们有更多的数据剩余时，尝试解码SIB1而(size(rxgrid,2) >) fprintf(“% s \ n”、分离器);流('帧%d\n的SIB1解码',国防部(enb.NFrame, 1024));流(“% s \ n \ n”、分离器);重置HARQ缓冲区，每8帧为SIB1% info可能不同如果(mod (enb.NFrame 8) = = 0)流(“重置HARQ缓冲区。\n\n”）;decState = [];结束提取当前子帧rxsubframe = rxgrid(:，1:L，:);%执行信道估计[est,nest] = lteDLChannelEstimate(enb, cec, rxsubframe);% pcfi解调，CFI解码。CFI现在被解调并且%解码使用类似的资源提取和解码功能来%已显示用于BCH接收的内容。lteExtractResources用于从接收到的子帧中提取与PCFICH相对应的REs% rx子帧和信道估计测试。流(“解码CFI。\ n \ n”）;pcfichinices = ltepcfichinices (enb);%获取PCFICH索引[pcfichRx, pcfichest] = lteExtractResources(pcfichIndices, rxsubframe, hest);解码PCFICHcfiBits = ltePCFICHDecode(enb, pcfichRx, pcfichest, nest);cfi = lteCFIDecode(cfiBits);%获得CFI如果(isfield (enb“CFI”) && cfi~=enb.CFI) release(pdcchConstDiagram);结束enb。CFI = CFI;流('已解码的CFI值:%d\n\n', enb.CFI);%对于TDD，必须对PDCCH的可能值进行盲解码在TS36.211表6.9-1中，PHICH配置因子m_i (0,1,2) %% m_i = 0,1和2可以通过配置TDD实现%上行-下行配置分别为1、6和0。如果(strcmpi (enb。DuplexMode,“TDD”tddConfigs = [1 60 0];其他的tddConfigs = 0;%不用于FDD，只用于控制while循环结束Alldci = {};而(isempty(alldci) && ~isempty(tddConfigs))%配置TDD上行下行配置如果(strcmpi (enb。DuplexMode,“TDD”enb))。TDDConfig = tddConfigs(1);结束tddConfigs(1) = [];% PDCCH解调。现在对PDCCH进行解调和解码%使用类似的资源提取和解码功能BCH和CFI接收的百分比已显示pdcchIndices = ltePDCCHIndices(enb);%获取PDCCH索引[pdcchRx, pdcchHest] = lteExtractResources(pdcchIndices, rxsubframe, hest);解码PDCCH并绘制星座[dciBits, pdcchSymbols] = ltePDCCHDecode(enb, pdcchRx, pdcchHest, nest);pdcchConstDiagram (pdcchSymbols);% PDCCH盲搜索系统信息(SI)和DCI解码。LTE工具箱提供了全盲查找的PDCCH来查找%具有指定RNTI的任何DCI消息，在本例中为SI-RNTI。流(“PDCCH搜索SI-RNTI…\n\n”）;PDCCH = struct(“RNTI”, 65535);pdcch。ControlChannelType =“PDCCH”;pdcch。EnableCarrierIndication =“关闭”;pdcch。SearchSpace =“普通”;pdcch。EnableMultipleCSIRequest =“关闭”;pdcch。EnableSRSRequest =“关闭”;pdcch。NTxAnts = 1;alldci = ltePDCCHSearch(enb, pdcch, dciBits);%在PDCCH中搜索DCI结束%如果DCI已解码，则继续解码PDSCH / DL-SCH为I = 1:num (alldci) dci = alldci{I};流('带有SI-RNTI的DCI消息:\n'）;disp (dci);%从DCI获取PDSCH配置[pdsch, trblklen] = hPDSCHConfiguration(enb, dci, pdsch . rnti);%如果创建了PDSCH配置，继续解码PDSCH% / dl-sch如果~ isempty (pdsch) pdsch。NTurboDecIts = 5;流(' DCI解码后的PDSCH设置:\n'）;disp (pdsch);% PDSCH解调和DL-SCH解码恢复SIB位。DCI消息现在被解析为给出的配置%对应的携带SIB1的PDSCH，则PDSCH为%解调，最后对接收到的比特进行DL-SCH解码%产生SIB1位。流(“解码SIB1。\ n \ n”）;%获取PDSCH索引[pdschIndices, pdschindesinfo] = ltePDSCHIndices(enb, pdsch, pdsch. prbset);[pdschRx, pdschHest] = lteExtractResources(pdschIndices, rxsubframe, hest);%解码PDSCH[dlschBits,pdschSymbols] = ltePDSCHDecode(enb, pdschch, pdschRx, pdschHest, nest);解码DL-SCH与软缓冲输入/输出HARQ组合如果~ isempty (decState)流(“与之前的传输重新组合。\n\n”）;结束[sib1, crc, decState] = lteDLSCHDecode(enb, pdsch, trblklen, dlschBits, decState);%计算PDSCH EVMrecoded = lteDLSCH(enb, pdsch, pdschIndicesInfo.)G, sib1);remod = ltePDSCH(enb, pdsch, recoded);[~，refSymbols] = ltePDSCHDecode(enb, pdsch, remod);[rmsevm,peakevm] = pdschEVM(refSymbols{1}， pdschSymbols{1});流('PDSCH RMS EVM: %0.3f%%\n', rmsevm);流('PDSCH峰值EVM: %0.3f%%\n\n', peakevm);流('SIB1 CRC: %d\n', crc);如果CRC == 0“成功恢复SIB1。\n\n”）;其他的流('SIB1解码失败。\n\n'）;结束其他的%表示从DCI创建PDSCH配置%消息失败流('从DCI消息创建PDSCH配置失败。\n\n'）;结束结束如果(元素个数(alldci) = = 0)% DCI解码失败流('DCI解码失败。\n\n'）;结束%更新信道估计图图(channelFigure);冲浪(abs(命令(:,:,1,1)));hSIB1RecoveryExamplePlots (channelFigure);channelFigure.CurrentAxes.XLim = [0 size(hest,2)+1];channelFigure.CurrentAxes.YLim = [0 size(hst,1)+1];跳过2帧并再次尝试SIB1解码，如果失败则终止剩下的帧数少于2帧。如果(大小(rxgrid, 2) > = (L * 20)) rxgrid (:, 1: (L * 20 ),:) = [];再移除2帧其他的Rxgrid = [];%剩余少于2帧结束enb。NFrame = mod(enb);NFrame + 2,1024);结束