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

如果loadFromFile负载eNodeBOutput.mat%负载eNodeB输出的I/Q捕获eNodeBOutput = (eNodeBOutput) / 32768年翻一番;%鳞片样品sr = 15.36 e6;载样取样率%其他的rmc = lteRMCDL (“R.3”）;%#好的rmc.NCellID=17；rmc.TotSubframes=41；rmc.PDSCH.RNTI=61；% SIB参数rmc.SIB.Enable=“上”; rmc.SIB.DCIFormat=“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））channelFigure=figure(“可见”,“关”）;终止(简介、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））fprintf(“\n收到的信号不能为空。\n”）;返回；终止%显示接收信号频谱流(“\n正在设置接收信号频谱…\n”）;简介(awgn (eNodeBOutput, 100.0));如果(sr ~ = ofdmInfo.SamplingRate)如果(sr < ofdmInfo.SamplingRate)警告('接收信号采样率(%0.3fMs/s)低于cell搜索/ MIB解码所需的采样率(%0.3fMs/s);cell搜索/ MIB解码可能失败。，sr/1e6，ofdmInfo.SamplingRate/1e6）；终止流(“\n正在从%0.3fMs/s重新采样到%0.3fMs/s以进行单元搜索/MIB解码…\n”，sr/1e6，ofdmInfo.SamplingRate/1e6）；其他的流(“\n不需要重新采样；接收到的信号处于小区搜索/MIB解码所需的采样率（%0.3fMs/s）。\n1、sr / e6);终止% Downsample收到信号nSamples =装天花板(ofdmInfo.SamplingRate /轮(sr) *大小(eNodeBOutput, 1));nRxAnts = size(eNodeBOutput, 2);downsampling = 0 (nSamples, nRxAnts);对于i=1:nRxAnts downsampling (:，i) = resample(eNodeBOutput(:，i)， ofdmInfo. .SamplingRate圆(sr));终止

流(“\n正在执行单元格搜索…\n”）;%设置搜索的双工模式和循环前缀长度组合;如果%在| enb |中配置这些参数中的任何一个，则该值为%被认为是正确的如果(~ isfield (enb“双工模式”)) 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，offset，peak]=lteCellSearch（enb，下采样，searchalg）；enb.NCellID=enb.NCellID（1）；偏移量=偏移量（1）；峰值=峰值（1）；如果（峰值>峰值）enbMax=enb；offsetMax=偏移量；peakMax=峰值；终止终止终止%使用小区标识、循环前缀长度、双工模式和定时在单元格搜索期间给出最大相关性的偏移量enb = enbMax;抵消= offsetMax;%计算三个可能的主单元的相关性%的身份;建立了细胞同一性的相关性峰值%以上与其他两个相关的峰值进行比较%初级细胞的身份，以建立质量%的相关性。相关系数=细胞(1、3);idGroup =地板(enbMax.NCellID / 3);对于i=0:2 enb.NCellID=idGroup*3+mod（enbMax.NCellID+i，3）；[~，corr{i+1}]=lteDLFrameOffset（enb，下采样）；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',抵消);downsampled = downsampled(1 +抵消:最终,);enb。NSubframe = 0;显示单元格范围设置流(' cell -wide settings after cell search:\n')(九); disp(enb);；

流('\n '执行频率偏移估计…\n'）;%对于TDD，TDDConfig和SSC默认为0。这些参数不适用%在系统中建立，直到SIB1被解码，因此在这个阶段%值为0表示最保守的假设(最少的下行链路)%子帧和最短的特殊子帧)。如果(strcmpi (enb。DuplexMode,“TDD”))enb.TDDConfig=0；enb.SSC=0；终止delta_f=LTE频率偏移（enb，下采样）；fprintf(“频率偏移:% 0.3 fhz \ n”, delta_f);downsampled = lteFrequencyCorrect(enb, downsampled, delta_f); / /采样

通道估计器配置cec。PilotAverage =“用户定义”；%导频平均类型cec。FreqWindow = 13;频率窗口大小cec。TimeWindow = 9;%时间窗口大小cec.InterpType=“立方”；二维插值型电子窗口=“居中”；插值窗口类型cec。InterpWinSize = 1;插值窗口大小%假设初始解码尝试使用4个细胞特定的参考信号;%确保信道估计值可用于所有特定于小区的参考%信号enb.CellRefP=4；fprintf(“执行OFDM解调。\ n \ n”）;griddims = lteResourceGridSize (enb);%资源网格维度L=griddims（2）；%子帧中的OFDM符号数%OFDM解调信号rxgrid = lteofdm解调器(enb，下采样);如果(isempty (rxgrid))流('在定时同步后，信号小于一个子帧，因此不会执行进一步的解调'）;返回；终止执行信道估计[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解码后的Cell-wide设置:\n')(九); disp(enb);；如果(enb.CellRefP = = 0)流('MIB解码失败(enb.CellRefP=0).\n\n'）;返回；终止如果（enb.NDLRB==0）fprintf('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'，sr/1e6，ofdmInfo.SamplingRate/1e6）；其他的流(“\n不需要重新采样；接收到的信号为NDLRB=%d（%0.3fMs/s）的所需采样率。\n”，enb.NDLRB，sr/1e6）；终止nSamples =装天花板(ofdmInfo.SamplingRate /轮(sr) *大小(eNodeBOutput, 1));resampled = 0 (nSamples, nRxAnts);对于i = 1:nRxAnts resampled(:，i) = resample(eNodeBOutput(:，i)， ofdmInfo. .SamplingRate圆(sr));终止%进行频率偏差估计和校正流('\n '执行频率偏移估计…\n'）;delta_f = lteFrequencyOffset(enb, resampled); / /重新采样流(“频率偏移:% 0.3 fhz \ n”，delta_f）；重采样=LTE频率正确（enb，重采样，delta_f）；找到帧的开始流(“\n正在执行定时偏移估计…\n”）;offset = lteDLFrameOffset(enb, resampled);流('帧开始的时间偏移：%d个样本\n',抵消);%将信号与帧的开始对齐重新取样=重新取样(1 +抵消:最终,);%正交频分复用解调流(“\n正在执行OFDM解调…\n\n”）;rxgrid = lteofdm解调器(enb，重采样);

% Check this frame contains SIB1, if not advance by 1 frame provided we . %检查这个帧包含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”，分隔符）；fprintf('SIB1解码帧%d\n',国防部(enb.NFrame, 1024));流(“% s \ n \ n”、分离器);%使用每一组新的8帧作为SIB1重置HARQ缓冲区% 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的REs% rxsubframe和channel estimate hest。流(“解码CFI。\ n \ n”）;pcfichIndices = ltePCFICHIndices (enb);%获取PCFICH指数[pcfichRx，pcfichHest]=lteExtractResources（pcfichIndices，rxsubframe，hest）；%解码PCFICHcfiBits=ltePCFICHDecode（enb、pcfichRx、pcfichHest、nest）；cfi=lteCFIDecode（cfiBits）；%得到CFI如果（伊斯菲尔德（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。流('PDCCH搜索SI-RNTI…\n\n'）;pdcch =结构(“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};流('带有SI-RNTI的DCI消息：\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，pdschIndicesInfo]=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);如果crc==0 fprintf(“成功SIB1复苏。\ n \ n”）;其他的流(“SIB1解码失败。\ n \ n”）;终止其他的%表示从DCI创建PDSCH配置%消息失败流(“从DCI消息创建PDSCH配置失败。\n\n”）;终止终止如果（numel（所有DCI）==0）%指示DCI解码失败流(“DCI解码失败。\ n \ n”）;终止%更新信道估计图图(channelFigure);冲浪(abs(命令(:,:,1,1)));hSIB1RecoveryExamplePlots (channelFigure);channelfigel . currentaxis . xlim = [0 size(hst,2)+1];channelfigel . currentaxis . ylim = [0 size(hst,1)+1]; / /当前位置%跳过2帧并再次尝试SIB1解码，或者终止%剩下的帧少于2帧。如果（大小（rxgrid，2）>=（L*20））rxgrid（：，1:（L*20），：）=[]；%再移除2个框架其他的rxgrid = [];%少于2帧剩余终止enb。NFrame =国防部(enb。NFrame + 2, 1024);终止