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

如果loadFromFile负荷eNodeBOutput.mat％负载I / e节点B输出的Q Capture的eNodeBOutput =双（eNodeBOutput）/ 32768;%样本规模SR = 15.36e6;％采样速率加载的样品其他的RMC = lteRMCDL（'R.3');％＃确定rmc。NCellID = 17;rmc。TotSubframes = 41;rmc.PDSCH。RNTI = 61;％SIB参数rmc.SIB。使='上';rmc.SIB。DCIFormat ='Format1A';rmc.SIB。所有ocationType = 0; rmc.SIB.VRBStart = 8; rmc.SIB.VRBLength = 8; rmc.SIB.Gap = 0;充满随机位％SIB数据字段，这不是有效的SIB％ 信息rmc.SIB.Data =兰迪（[0 1]，176,1）;[eNodeBOutput，〜，信息] = lteRMCDLTool（RMC，[1; 0; 0; 1]）;SR = info.SamplingRate;生成的样品的％采样速率结束

％设置一些看家变量：％分离器，用于命令窗口记录隔板= repmat（' - '，1,50）;地块％如果（〜存在（“channelFigure”，“var”）||〜的isValid（channelFigure））channelFigure =系数（“可见”，“关”);结束[spectrumAnalyzer，synchCorrPlot，pdcchConstDiagram] =...hSIB1RecoveryExamplePlots（channelFigure，SR）;％PDSCH EVMpdschEVM = comm.EVM（）;pdschEVM.MaximumEVMOutputPort = TRUE;％，为初始小区搜索的采样速率，使用既定的％lteOFDMInfo配置为6个资源块。enb.CyclicPrefix是集％暂居调用lteOFDMInfo压制默认值％警告（不影响采样率）。基站=结构;％的eNodeB配置结构enb.NDLRB = 6;资源块的数量％ofdmInfo = lteOFDMInfo (setfield (enb'CyclicPrefix'，'正常'））;％＃确定如果（的isEmpty（eNodeBOutput））fprintf中（“\ nReceived信号不能为空。\ n”);返回;结束％显示接收信号频谱fprintf中（'绘制接收信号频谱…');spectrumAnalyzer（AW​​GN（eNodeBOutput，100.0））;如果（SR〜= ofdmInfo.SamplingRate）如果（SR “接收到的信号采样率（％0.3fMs / S）比为小区搜索所需的采样率下/ MIB解码（％0.3fMs / S）;小区搜索/ MIB解码可能失败“。，SR / 1E6，ofdmInfo.SamplingRate / 1E6）;结束fprintf中（'\ nResampling来自％0.3fMs / s到％0.3fMs / s的小区搜索/ MIB解码... \ N'，SR / 1E6，ofdmInfo.SamplingRate / 1E6）;其他的fprintf中（' \ nResampling不是必需的;接收到的信号以期望的采样率进行小区搜索/ MIB解码(%0.3fMs/s)。，SR / 1E6）;结束％缩减像素采样接收到的信号NSAMPLES =小区（ofdmInfo.SamplingRate /轮（SR）*大小（eNodeBOutput，1））;nRxAnts =尺寸（eNodeBOutput，2）;降频采样=零（NSAMPLES，nRxAnts）;对于I = 1：nRxAnts下采样（：，1）=重采样（eNodeBOutput（：，i）中，ofdmInfo.SamplingRate，圆（SR））;结束

fprintf中（'\ nPerforming小区搜索... \ N');%设置双工模式和循环前缀长度组合进行搜索;如果这些参数中的任何一个在|enb|中配置，则其值为％假设是正确的如果(~ isfield (enb'DuplexMode'））duplexModes = {“TDD”'FDD'};其他的duplexModes = {enb.DuplexMode};结束如果(~ isfield (enb'CyclicPrefix'））cyclicPrefixes = {'正常'“扩展”};其他的cyclicPrefixes = {enb.CyclicPrefix};结束％进行跨越双工模式和循环前缀长度的小区搜索％的组合和记录与所述最大相关的组合;如果%多单元搜索已配置，本例将解码第一个单元%(最强)检测到的细胞searchalg。MaxCellCount = 1;searchalg。SSSDetection =“PostFFT”;peakMax = -Inf;对于duplexMode = duplexModes对于cyclicPrefix = cyclicPrefixesDuplexMode = DuplexMode {1};enb。CyclicPrefix = CyclicPrefix {1};[enb。NCellID, offset, peak] = lteCellSearch(enb, down, searchalg);enb。NCellID = enb.NCellID (1);抵消=抵消(1);峰=峰(1);如果（峰值> peakMax）enbMax =基站;offsetMax =偏移;peakMax =峰;结束结束结束％使用小区标识，循环前缀长度，双工模式和定时％偏置这给了小区搜索时的最大相关基站= enbMax;偏移量= offsetMax;％计算每个所述三种可能的主小区的相关性％的身份;对于小区身份的相关性建立的峰值％以上与相关的其它两个峰相比％，以主小区标识建立的质量％的相关性。科尔=细胞（1,3）;idGroup =地板（enbMax.NCellID / 3）;对于I = 0：2 enb.NCellID = idGroup * 3 + MOD（enbMax.NCellID + 1，3）;[〜，科尔{I + 1}] = lteDLFrameOffset（ENB，下采样）;科尔{I + 1} =总和（相应{I + 1}，2）;结束阈值= 1.3 * max([corr{2};corr {3}]);根据经验得出的乘数为1.3如果（最大（相应{1}）<阈值）警告（同步信号相关性弱;检测到的细胞身份可能不正确。”);结束%返回到原来检测到的单元格标识enb.NCellID = enbMax.NCellID;%绘制PSS/SSS相关性和阈值synchCorrPlot。max([corr{1};阈值)* 1.1);{1}synchCorrPlot ([corr阈值*的(大小{1})(corr)]);％执行定时同步fprintf中（帧开始的定时偏移:%d采样\n'，偏移量）;下采样下采样=（1 +偏移：端，:);enb.NSubframe = 0;％显示单元格范围的设置fprintf中（“小区搜索后的细胞级设置：\ N”);DISP（ENB）;

fprintf中（'\ nPerforming频偏估计... \ N');％对于TDD，TDDConfig和SSC被默认为0。这些参数不％建立的系统，直到SIB1解码，所以在这个阶段，0％值做出最保守的假设（最少下行％子帧和最短的特殊子帧）。如果（strcmpi（enb.DuplexMode，“TDD”enb))。TDDConfig = 0;enb。SSC = 0;结束delta_f = lteFrequencyOffset（ENB，下采样）;fprintf中（“频率偏移:% 0.3 fhz \ n”,delta_f);下采样= ltefrequency - correct (enb，下采样，delta_f);

％信道估计的配置cec.PilotAverage =“UserDefined”;导频平均化的％类型cec。FreqWindow = 13;频率窗口大小cec.TimeWindow = 9;％时间窗口大小cec。InterpType =“立方”;％2D插值类型cec.InterpWindow =“中心”;％插值窗型cec.InterpWinSize = 1;插值窗口大小％假设为初始解码尝试4小区特定参考信号;％，确保信道估计可用于所有小区特定的参考％的信号enb.CellRefP = 4;fprintf中（“执行OFDM解调。\ n \ n”);griddims = lteResourceGridSize（ENB）;％资源网格尺寸L = griddims（2）;在子帧中的OFDM符号的数量％％OFDM解调信号rxgrid = lteOFDMDemodulate（ENB，下采样）;如果（的isEmpty（rxgrid））fprintf中（在定时同步后，信号短于一个子帧，因此不会进行进一步的解调。);返回;结束％进行信道估计[HEST，巢] = lteDLChannelEstimate（ENB，CEC，rxgrid（：，1：L，:)）;

％在解码MIB从第一对应于PBCH％提取资源元素（RE）%子帧通过所有接收天线和信道估计fprintf中（'表演MIB解码... \ N');pbchIndices = ltePBCHIndices（ENB）;[pbchRx，pbchHest] = lteExtractResources（...pbchIndices，rxgrid（：，1：L，:)，HEST（：，1：L，：，:)）;％解码PBCH[bchBits，pbchSymbols，nfmod4，MIB，enb.CellRefP] = ltePBCHDecode（...基站，pbchRx，pbchHest，巢）;％解析MIB位enb = lteMIB(mib, enb);将来自ltePBCHDecode函数的nfmod4值输出合并为％从MIB建立NFrame值是系统帧号％（SFN）模4（它是存储在MIB地板（SFN / 4））enb.NFrame = enb.NFrame + nfmod4;MIB解码之后％显示细胞范围的设置fprintf中（' MIB解码后的小区范围设置:\n');DISP（ENB）;如果(enb.CellRefP = = 0)流('MIB解码失败（enb.CellRefP = 0）。\ n \ n');返回;结束如果（enb.NDLRB == 0）fprintf中（'MIB解码失败（enb.NDLRB = 0）。\ n \ n');返回;结束

fprintf中（'现在重新启动接收带宽（NDLRB =％d）是已知的... \ N'，enb.NDLRB）;％重新取样现在我们知道真正带宽ofdmInfo = lteOFDMInfo（ENB）;如果（SR〜= ofdmInfo.SamplingRate）如果（SR “接收到的信号采样率（％0.3fMs / S）是比NDLRB =％d所需的采样率下（％0.3fMs / S）;PDCCH搜索/ SIB1解码可能失败“。，SR / 1E6，enb.NDLRB，ofdmInfo.SamplingRate / 1E6）;结束fprintf中（'\ nResampling从％0.3fMs / s到％0.3fMs /秒... \ N'，SR / 1E6，ofdmInfo.SamplingRate / 1E6）;其他的fprintf中（' \ nResampling不是必需的;接收的信号是在用于NDLRB =％d所需的采样率（％0.3fMs / S）。\ N”，enb.NDLRB，SR / 1E6）;结束NSAMPLES =小区（ofdmInfo.SamplingRate /轮（SR）*大小（eNodeBOutput，1））;重采样=零（NSAMPLES，nRxAnts）;对于i = 1:nRxAnts resam(:，i) = resample(eNodeBOutput(:，i)， ofdmInfo。SamplingRate圆(sr));结束％进行频率偏移估计和校正fprintf中（'\ nPerforming频偏估计... \ N');delta_f = lteFrequencyOffset（ENB，重采样）;fprintf中（“频率偏移:% 0.3 fhz \ n”,delta_f);重采样= ltefrequency - correct (enb，重采样，delta_f);％查找帧的开始fprintf中（执行定时偏移估计…);偏移量= lteDLFrameOffset（ENB，重采样）;fprintf中（帧开始的定时偏移:%d采样\n'，偏移量）;调整信号与帧的开始重新采样的再采样=（1 +偏移：端，:);% OFDM解调fprintf中（“\ nPerforming OFDM解调。\ n \ n”);rxgrid = lteOFDMDemodulate（ENB，重采样）;

%检查此帧是否包含SIB1，如果我们没有提前1帧％有足够的数据，否则终止。如果（MOD（enb.NFrame，2）〜= 0）如果（大小（rxgrid，2）> =（L * 10））rxgrid（：，1：（L * 10），:) = [];fprintf中（'跳过帧%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））fprintf中（“接收到的信号中不包含的子帧携带SIB1。\ n \ n”);结束％重置HARQ缓冲器decState = [];％，而我们有更多的数据离开了，尝试解码SIB1而（大小（rxgrid，2）> 0）fprintf中（“% s \ n”、分离器);fprintf中（'SIB1解码为帧％d \ N'，MOD（enb.NFrame，1024））;fprintf中（'％s的\ n \ n'、分离器);用每组新的8帧作为SIB1重置HARQ缓冲区％的信息可能会有所不同如果（MOD（enb.NFrame，8）== 0）fprintf中（“重置HARQ缓冲区。\ n \ n”);decState = [];结束提取当前子帧rxsubframe = rxgrid（：，1：L，:);％进行信道估计[HEST，巢] = lteDLChannelEstimate（ENB，CEC，rxsubframe）;％PCFICH解调，CFI进行解码。原讼法庭现在解调%解码使用类似的资源提取和解码功能％那些BCH接收已经示出。lteExtractResources用于从接收的子帧中提取与PCFICH对应的REs％rxsubframe和信道估计HEST。fprintf中（'解码CFI ... \ n \ n');pcfichIndices = ltePCFICHIndices（ENB）;%得到PCFICH指数[pcfichRx，pcfichHest] = lteExtractResources（pcfichIndices，rxsubframe，HEST）;％解码PCFICHcfiBits = ltePCFICHDecode(enb, pcfichRx, pcfichHest, nest);cfi = lteCFIDecode (cfiBits);%得到CFI如果(isfield (enb'CFI'）&& CFI〜= enb.CFI）释放（pdcchConstDiagram）;结束enb.CFI = CFI;fprintf中（'解码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 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，pdcchHest] = lteExtractResources（pdcchIndices，rxsubframe，HEST）;％PDCCH进行解码和情节的星座[dciBits，pdcchSymbols] = ltePDCCHDecode（ENB，pdcchRx，pdcchHest，巢）;pdcchConstDiagram（pdcchSymbols）;％PDCCH盲目搜索系统信息（SI）和DCI解码。％的LTE工具箱提供全盲搜索PDCCH的发现带有指定RNTI(在本例中为SI-RNTI)的任何DCI消息。fprintf中（'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）;％搜索PDCCH的DCI结束％如果DCI被解码，则继续PDSCH / DL-SCH的解码对于I = 1：numel（alldci）DCI = alldci {I};fprintf中（带SI-RNTI的DCI信息:\n');DISP（DCI）;从DCI获取PDSCH配置[pdsch, trblklen] = hPDSCHConfiguration(enb, dci, pdcch.RNTI);％如果PDSCH创建配置的，则继续进行解码PDSCH% / DL-SCH如果〜的isEmpty（PDSCH）pdsch.NTurboDecIts = 5;fprintf中（DCI解码后的PDSCH设置:\n');DISP（PDSCH）;％PDSCH解调和DL-SCH的解码，以恢复SIB比特。％的DCI消息现在解析，得到的结构%对应的PDSCH携带SIB1，该PDSCH为％解调，最后所接收的位是DL-SCH的解码%生成SIB1位。fprintf中（'解码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，巢）;％译码DL-SCH与HARQ组合的软缓冲器的输入/输出如果〜的isEmpty（decState）fprintf中（'与先前的传输重组。\ n \ n');结束[SIB1，CRC，decState] = lteDLSCHDecode（ENB，PDSCH，trblklen，dlschBits，decState）;％计算PDSCH EVM重新编码的= lteDLSCH（ENB，PDSCH，pdschIndicesInfo.G，SIB1）;remod = ltePDSCH（ENB，PDSCH，重新编码）;[〜，refSymbols] = ltePDSCHDecode（ENB，PDSCH，remod）;[rmsevm，peakevm] = pdschEVM（refSymbols {1}，pdschSymbols {1}）;fprintf中（'PDSCH RMS EVM: %0.3f%%\n'，rmsevm）;fprintf中（'PDSCH峰值EVM：％0.3f %% \ n \ n'，peakevm）;fprintf中（'SIB1 CRC：％d \ n'CRC）;如果CRC == 0 fprintf中（'成功SIB1恢复。\ n \ n');其他的fprintf中（'SIB1解码失败。\ n \ n');结束其他的表示从DCI创建一个PDSCH配置％消息失败fprintf中（'从DCI消息创建PDSCH的配置失败。\ n \ n');结束结束如果(元素个数(alldci) = = 0)％表示DCI解码失败fprintf中（“DCI解码失败。\ n \ n”);结束％更新信道估计情节图（channelFigure）;冲浪（ABS（HEST（：，：，1,1）））;hSIB1RecoveryExamplePlots（channelFigure）;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 = MOD（enb.NFrame + 2,1024）;结束