loadFromFile = 1;％设置为0以在本地生成eNodeB输出

如果loadFromFile负载eNodeBOutput.mat%负载eNodeB输出的I/Q捕获eNodeBOutput = (eNodeBOutput) / 32768年翻一番;％缩放样本sr = 15.36 e6;载样取样率%别的rmc = lteRMCDL (“R.3”）;% #好< UNRCH >rmc.ncellid = 17;rmc.totsubframes = 41;rmc.pdsch.rnti = 61;% SIB参数rmc.sib.enable =.“上”；rmc.sib.dcormat =“Format1A”；rmc.sib.allocationtype = 0;rmc.sib.vrbstart = 8;rmc.sib.vrblength = 8;rmc.sib.gap = 0;% SIB数据字段填充随机位，这不是一个有效的SIB%的消息rmc.SIB.Data = randi([0 1]，176,1);[eNodeBOutput ~,信息]= lteRMCDLTool (rmc中,[1,0,0,1]);sr = info.SamplingRate;%生成样本的采样率结尾

%设置一些家政变量:命令窗口日志记录的%分隔符分隔符= repmat (“- - -”1、50);%的情节如果(~ ('Channelfigure'，'var') || ~isvalid(channelFigure))'可见的'，“关闭”）;结尾(简介、synchCorrPlot pdcchConstDiagram] =.．.hSIB1RecoveryExamplePlots (channelFigure sr);% PDSCH维生素pdschEVM = comm.EVM ();pdschEVM。MaximumEVMOutputPort = true;初始单元搜索的抽样率是用为6个资源块配置% lteOFDMInfo。enb。CyclicPrefix设置%在调用lteOFDMInfo时暂时抑制默认值%警告(不影响采样率)。enb =结构;% eNodeB配置结构enb。NDLRB = 6;%资源块数OFDMINFO = LTEOFDMINFO（SETFIELD（eNB，“CyclicPrefix”，“正常”))；% #好< SFLD >如果(isempty (eNodeBOutput))流(接收到的信号不能为空。）;返回；结尾%显示接收信号频谱流('\ npplotting接收信号谱... \ n'）;简介(awgn (eNodeBOutput, 100.0));如果(sr ~ = ofdmInfo.SamplingRate)如果(sr < ofdmInfo.SamplingRate)警告('接收信号采样率(%0.3fMs/s)低于cell搜索/ MIB解码所需的采样率(%0.3fMs/s);cell搜索/ MIB解码可能失败。1、sr / 1 e6 ofdmInfo.SamplingRate / e6);结尾流('\nResampling from %0.3fMs/s to %0.3fMs/s for cell search / MIB decoding…\n'1、sr / 1 e6 ofdmInfo.SamplingRate / e6);别的流(' \ nResampling不是必需的;收到的信号是用于小区搜索/ MIB解码的所需采样率（％0.3fms / s）。\ n'1、sr / e6);结尾% Downsample收到信号nSamples =装天花板(ofdmInfo.SamplingRate /轮(sr) *大小(eNodeBOutput, 1));nrxants = size（eNodeboutput，2）;DownSampled = Zeros（nsamples，nrxants）;为了i=1:nRxAnts downsampling (:，i) = resample(eNodeBOutput(:，i)， ofdmInfo. .SamplingRate圆(sr));结尾

流(“搜索…\ n \ nPerforming细胞”）;％设置双工模式和循环前缀长度组合进行搜索;如果％在| eNB |配置了这些参数中的任何一个那么值是%被认为是正确的如果（〜isfield（eNB，“DuplexMode”)) duplexModes = {'TDD'“FDD”}；别的duplexmodes = {enb.duplexmode};结尾如果（〜isfield（eNB，“CyclicPrefix”)) cyclicPrefixes = {“正常”'延伸'}；别的cyclicPrefixes = {enb.CyclicPrefix};结尾％跨双面模式和循环前缀长度执行小区搜索%组合，并记录相关性最大的组合;如果％多个单元搜索配置，此示例将解码第一个％（最强）检测到的细胞searchalg.maxcellcount = 1;searchalg.sssdetection ='postfft'；peakMax =无穷;为了duplexMode = duplexModes为了cyclicprefix = cyclicprefixes enb.duplexmode = duplexmode {1};enb.cyclicprefix = cyclicprefix {1};[enb.ncellid，偏移，峰值] = LTECellSearch（eNB，DownSampled，SearchaL）;eNB.ncellid = eNB.ncellid（1）;offset =偏移（1）;峰=峰（1）;如果(peak>peakMax) enbMax = enb;offsetMax =抵消;peakMax =峰值;结尾结尾结尾%使用小区标识、循环前缀长度、双工模式和定时在单元格搜索期间给出最大相关性的偏移量enb = enbMax;抵消= offsetMax;%计算三个可能的主单元的相关性%的身份;建立了细胞同一性的相关性峰值%以上与其他两个相关的峰值进行比较%初级细胞的身份，以建立质量%的相关性。相关系数=细胞(1、3);idGroup =地板(enbMax.NCellID / 3);为了I = 0:2 enb。NCellID = idGroup*3 + mod(enbMax。NCellID +我,3);[~，corr{i+1}] = lteDLFrameOffset(enb, downsampling);corr {i + 1} =总和(corr {i + 1}, 2);结尾阈值= 1.3 * max（[corr {2}; corr {3}]）;经验获得的1.3％乘数如果{1} (max (corr) <阈值)警告('同步信号相关性弱;检测到的细胞标识可能是不正确的。）;结尾％返回最初检测到的细胞标识enb。NCellID = enbMax.NCellID;％plot pss / sss相关性和阈值synchCorrPlot。YLimits = [0 max([corr{1};阈值)* 1.1);{1} synchCorrPlot ([corr阈值*的(大小{1})(corr)]);%执行定时同步流('定时偏移到帧开始：％d样本\ n',抵消);downsampled = downsampled(1 +抵消:最终,);enb。NSubframe = 0;显示单元格范围设置流(' cell -wide settings after cell search:\n'）;disp (enb);

流('\ n qualforming频率偏移估计... \ 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）;DownSampled = LtefRequencyCorrect（eNB，DownPled，Delta_F）;

通道估计器配置cec。PilotAverage ='用户自定义'；%引航平均的类型cec。FreqWindow = 13;％频率窗口大小cec。TimeWindow = 9;%时间窗口大小cec.interptype =.'立方体'；二维插值型cec。InterpWindow ='中心'；插值窗口类型cec。InterpWinSize = 1;％插值窗口大小%假设初始解码尝试使用4个细胞特定的参考信号;%确保通道估计可用于所有特定单元格的参考%的信号enb.cellrefp.= 4; fprintf('执行OFDM解调... \ n \ n'）;griddims = lteResourceGridSize (enb);%资源网格维度L = griddims (2);%一帧OFDM符号数% OFDM解调信号rxgrid = lteofdm解调器(enb，下采样);如果(isempty (rxgrid))流('在定时同步后，信号短于一个子帧，因此不会执行进一步的解调。\ n'）;返回；结尾执行信道估计[hest, nest] = lteDLChannelEstimate(enb, cec, rxgrid(:，1:L，:));

%解码MIB从第一个中提取PBCH对应的资源元素(REs)%子帧横跨所有接收天线和信道估计流(“执行MIB解码…\ n”）;pbchIndices = ltePBCHIndices (enb);[pbchRx, pbchHest] = lteExtractResources(.．.pbchIndices rxgrid (::, 1: L),命令(:,:,1:L:));%解码PBCH[bchBits, pbchSymbols, nfmod4, mib, enb.]CellRefP] = ltePBCHDecode (.．.enb, pbchRx, pbchest, nest);%解析MIB位eNB = LTEMIB（MIB，eNB）;％从函数LTEPBCHDECODE中加入NFMOD4值输出，如%从MIB中建立的NFrame值就是系统帧号% (SFN) modulo 4(在MIB中存储为floor(SFN/4))enb。NFrame = enb.NFrame + nfmod4;% MIB解码后显示单元宽设置流('MIB解码后的单元格设置：\ n'）;disp (enb);如果（enb.cellrefp == 0）fprintf（'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);结尾流('\nResampling from %0.3fMs/s to %0.3fMs/s…\n'1、sr / 1 e6 ofdmInfo.SamplingRate / e6);别的流(' \ nResampling不是必需的;接收信号在NDLRB=%d (%0.3fMs/s)的期望采样率下。\n'1、enb.NDLRB sr / e6);结尾nSamples =装天花板(ofdmInfo.SamplingRate /轮(sr) *大小(eNodeBOutput, 1));resampled = 0 (nSamples, nRxAnts);为了i = 1：nrxants重新采样（：，i）=重组（eNodeboutput（：，i），Ofdminfo.samplingRate，Round（SR））;结尾%进行频率偏差估计和校正流('\ n qualforming频率偏移估计... \ n'）;delta_f = lteFrequencyOffset(enb, resampled); / /重新采样流('频率偏移：％0.3FHz \ n'，delta_f）;重新采样= ltefrequencycorect（eNB，重新采样，delta_f）;找到帧的开始流('\ nuperforming定时偏移估计... \ n'）;offset = lteDLFrameOffset(enb, resampled);流('定时偏移到帧开始：％d样本\ n',抵消);用帧的开始％对准信号重新取样=重新取样(1 +抵消:最终,);％OFDM解调流('\ n ofdming 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个子帧终止如果(大小(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) > 0) fprintf('％s \ n'、分离器);流('SIB1解码帧%d\n',国防部(enb.NFrame, 1024));流(“% s \ n \ n”、分离器);％将HARQ缓冲区重置为SIB1的每个新的8帧% info可能不同如果(mod (enb.NFrame 8) = = 0)流('重置HARQ缓冲区。\ n \ n'）;decState = [];结尾％提取电流子帧rxsubframe = rxgrid (::, 1: L);执行信道估计[hest,nest] = lteDLChannelEstimate(enb, cec, rxsubframe);% PCFICH解调，CFI解码。CFI现在被解调使用类似的资源提取和解码功能进行解码对于BCH接收而显示的％。lteextractresources用于％从接收的子帧的PCFICH对应的％提取物% rxsubframe和channel estimate hest。流(“解码CFI。\ n \ n”）;pcfichIndices = ltePCFICHIndices (enb);％获得PCFICH指数[pcfichRx, pcfichest] = lteExtractResources(pcfichIndices, rxsubframe, hest);%解码PCFICHcfibits = ltepcfichdecode（eNB，pcfichrx，pcfichhest，巢）;CFI = LTECFIDECODE（CFIBITS）;％获得CFI.如果（ISFIELD（eNB，“CFI”) && cfi~=enb.CFI)释放(pdcchConstDiagram);结尾enb。CFI = CFI;流('Decoded CFI值:%d\n\n', enb.CFI);%对于TDD，必须盲译码PDCCH的可能值% the PHICH configuration factor m_i (0,1,2) in TS36.211% m_i = 0,1和2的值可以通过配置TDD来实现%上行-下行配置分别为1、6和0。如果(strcmpi (enb。DuplexMode,'TDD'tddConfigs = [1 6 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, pdcchhst] = lteExtractResources(pdcchIndices, rxsubframe, hst); / /输出译码PDCCH并绘制星座图[dciBits, pdcchSymbols] = ltePDCCHDecode(enb, pdcchRx, pdcchest, nest);pdcchConstDiagram (pdcchSymbols);% PDCCH盲搜索系统信息(SI)和DCI解码。LTE工具箱提供了PDCCH的完全盲搜索来查找具有指定RNTI的任何DCI消息，在这种情况下，SI-RNTI。流(' pdch search for 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); / /查询所有数据%在pcch中搜索DCI结尾％如果解码DCI，则继续解码PDSCH / DL-SCH为了I = 1:numel(alldci) dci = alldci{I};流('DCI消息与SI-RNTI：\ n'）;disp (dci);％从DCI获取PDSCH配置[PDSCH，TRBLKLEN] = HPDSCHCONFIGURATION（eNB，DCI，PDCCH.rnti）;%如果已创建PDSCH配置，则继续解码PDSCH％/ DL-SCH如果~ isempty (pdsch) pdsch。NTurboDecIts = 5;流('DCI解码后的PDSCH设置：\ n'）;disp (pdsch);% PDSCH解调和DL-SCH译码恢复SIB位。%现在解析DCI消息以给出配置％PDSCH携带SIB1，PDSCH是%解调，最后对接收的比特进行DL-SCH解码%来生成SIB1位。流('解码sib1 ... \ n \ n'）;%获取PDSCH索引[pdschindices，pdschindicsinfo] = ltepdschindices（eNB，pdsch，pdsch.prbset）;[pdschrx，pdschhest] = lteextractresources（pdschindices，rxsubframe，hest）;％解码PDSCH.[dlschBits,pdschSymbols] = ltePDSCHDecode(enb, pdsch, pdschRx, pdschHest, nest); / /将pdsch解码解码DL-SCH与软缓冲输入/输出HARQ结合如果~ isempty (decState)流(“与之前的传输重组。”）;结尾[sib1, crc, decState] = lteDLSCHDecode(enb, pdsch, trblklen, dlschBits, decState);%计算PDSCH EVMencoded = lteDLSCH(enb, pdsch, pdschIndicesInfo. txt)。G, sib1);redd = ltePDSCH(enb, pdsch, recoded); / /重新编码[~，refSymbols] = ltePDSCHDecode(enb, pdsch, remod);[rmsevm,peakevm] = pdschEVM(refSymbols{1}， pdschSymbols{1});流('PDSCH RMS EVM：％0.3F %% \ n', rmsevm);流('PDSCH Peak EVM: %0.3f%%\n\n', peakevm);流(“SIB1 CRC: % d \ n”, crc);如果fprintf == 0 (“成功SIB1复苏。\ n \ n”）;别的流(“SIB1解码失败。\ n \ n”）;结尾别的％表示从DCI创建PDSCH配置%消息失败流(“从DCI消息创建PDSCH配置失败。\n\n”）;结尾结尾如果（NUMER（ALLDCI）== 0）% DCI解码失败流(“DCI解码失败。\ n \ n”）;结尾%更新信道估计图图（ChannFeigure）;冲浪（ABS（hest（：：，1,1））））;hsib1recoveryexampleplots（ChannFogigure）;ChannelFigure.currentaxes.xlim = [0大小（hest，2）+1];channelfigure.currentaxes.ylim = [0大小（hest，1）+1];％跳过2帧并再次尝试SIB1解码，或者如果我们终止％剩余少于2帧。如果（大小（rxgrid，2）> =（l * 20））rxgrid（：，1：（l * 20），:) = [];％删除了2个框架别的rxgrid = [];％少于2帧结尾enb。NFrame =国防部(enb。NFrame + 2, 1024);结尾