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

如果loadFromFile负载eNodeBOutput.mat% I / Q capture eNodeB输出负载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。使=“上”;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 =图(“可见”,“关闭”);结束(简介、synchCorrPlot pdcchConstDiagram] =…hSIB1RecoveryExamplePlots (channelFigure sr);% PDSCH维生素pdschEVM = comm.EVM ();pdschEVM。MaximumEVMOutputPort = true;%的采样率初始细胞建立搜索使用% lteOFDMInfo资源配置为6块。enb。CyclicPrefix设置%的临时调用lteOFDMInfo压制一个默认值%的警告(它不会影响采样率)。enb =结构;% eNodeB配置结构enb.NDLRB= 6;%资源块的数量ofdmInfo = lteOFDMInfo (setfield (enb“CyclicPrefix”,“正常”));% #好< SFLD >如果(isempty (eNodeBOutput))流(”不能为空。\ n \ nReceived信号”);返回;结束%显示接收信号频谱流(“\ nPlotting收到信号频谱…\ n”);简介(awgn (eNodeBOutput, 100.0));如果(sr ~ = ofdmInfo.SamplingRate)如果(sr < ofdmInfo.SamplingRate)警告(接收信号的采样率(% 0.3 fms / s)低于细胞搜索/ MIB译码所需的采样率(% 0.3 fms / s);手机搜索/ MIB解码可能会失败。”1、sr / 1 e6 ofdmInfo.SamplingRate / e6);结束流(的\ nResampling % 0.3 fms / s % 0.3 fms / s细胞搜索/ MIB解码…\ n '1、sr / 1 e6 ofdmInfo.SamplingRate / e6);其他的流(' \ nResampling不是必需的;接收信号在细胞搜索/ MIB译码所需的采样率(% 0.3 fms / s)。\ n”1、sr / e6);结束% Downsample接收信号nSamples =装天花板(ofdmInfo.SamplingRate /轮(sr) *大小(eNodeBOutput, 1));nRxAnts =大小(eNodeBOutput, 2);downsampled = 0 (nSamples nRxAnts);为i = 1: nRxAnts downsampled (:, i) =重新取样(eNodeBOutput(:,我),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 searchalg);enb。NCellID = enb.NCellID (1);抵消=抵消(1);峰=峰(1);如果(峰值> peakMax) enbMax = enb;offsetMax =抵消;peakMax =峰值;结束结束结束%使用细胞的身份,循环前缀长度、双工模式和时机%抵消使最大相关细胞搜索enb = enbMax;抵消= offsetMax;%计算相关的三种可能的主要细胞%的身份;细胞的相关身份的峰值%以上的峰值相比其他两个相关性%原代细胞的身份为了建立的质量%的相关性。相关系数=细胞(1、3);idGroup =地板(enbMax.NCellID / 3);为i = 0:2 enb。NCellID = idGroup * 3 +国防部(enbMax。NCellID +我,3);{i + 1}] [~, corr = lteDLFrameOffset (enb downsampled);corr {i + 1} =总和(corr {i + 1}, 2);结束阈值= 1.3 *马克斯([corr {2};corr {3}]);%的乘数为1.3经验如果{1}(max (corr) <阈值)警告(的同步信号相关性弱;检测细胞的身份可能是不正确的。”);结束%回到最初发现细胞的身份enb。NCellID = enbMax.NCellID;%情节PSS /瑞士相关性和阈值synchCorrPlot。YLimits =[0马克斯([corr {1};阈值)* 1.1);{1}synchCorrPlot ([corr阈值*的(大小{1})(corr)]);%完成定时同步流(帧的定时偏移量开始:\ n % d样本”,抵消);downsampled = downsampled(1 +抵消:最终,);enb。NSubframe = 0;%显示单元范围的设置流(“它是设置后细胞搜索:\ n”);disp (enb);

流(“……\ n \ nPerforming频差估计”);%对TDD, TDDConfig和SSC违约为0。这些参数不%建立于系统直到SIB1解码,所以在这个阶段%值0做最保守的假设(下行最少%子帧和最短的特殊的子帧)。如果(strcmpi (enb.DuplexMode“TDD”enb))。TDDConfig = 0;enb。SSC = 0;结束delta_f = lteFrequencyOffset (enb downsampled);流(“频率偏移:% 0.3 fhz \ n”,delta_f);downsampled = lteFrequencyCorrect (enb downsampled delta_f);

%信道估计配置cec。PilotAverage =“UserDefined”;%的试点,平均cec。FreqWindow = 13;%频率窗口大小cec。TimeWindow = 9;%的时间窗口大小cec。InterpType =“立方”;% 2 d插值类型cec。InterpWindow =“中心”;%插值窗口类型cec。InterpWinSize = 1;%插值窗口大小%假设4特异性首次解码尝试参考信号;%保证信道估计对所有特异性可用参考%的信号enb.CellRefP= 4; fprintf(“执行OFDM解调。\ n \ n”);griddims = lteResourceGridSize (enb);%资源网格尺寸L = griddims (2);%的OFDM符号数子帧% OFDM解调信号rxgrid = lteOFDMDemodulate (enb downsampled);如果(isempty (rxgrid))流(的定时同步后,信号比一副框架短所以没有进一步解调将被执行。\ n”);返回;结束%进行信道估计[命令,巢]= lteDLChannelEstimate (cec, enb rxgrid (::, 1: L));

%解码MIB%提取资源元素(REs)对应PBCH从第一%子帧所有接收天线和信道估计流(“执行MIB解码…\ n”);pbchIndices = ltePBCHIndices (enb);[pbchRx, pbchHest] = lteExtractResources (…pbchIndices rxgrid (::, 1: L),命令(:,:,1:L:));%解码PBCH[bchBits, pbchSymbols nfmod4 mib, enb。CellRefP] = ltePBCHDecode (…pbchRx, enb pbchHest、巢);%解析MIB比特enb = lteMIB (mib, enb);%将nfmod4 ltePBCHDecode输出函数值,% NFrame值建立MIB的系统帧数% (SFN)模4(存储在MIB地板(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.3 fms / s)低于所需的采样率NDLRB = % d (% 0.3 fms / s);PDCCH搜索/ SIB1解码可能会失败。”sr / 1 e6, enb.NDLRB ofdmInfo.SamplingRate / 1 e6);结束流(“\ nResampling从% 0.3 fms / s % 0.3 fms / s…\ n '1、sr / 1 e6 ofdmInfo.SamplingRate / e6);其他的流(' \ nResampling不是必需的;接收信号在NDLRB = % d所需的采样率(% 0.3 fms / s)。\ n”1、enb.NDLRB sr / e6);结束nSamples =装天花板(ofdmInfo.SamplingRate /轮(sr) *大小(eNodeBOutput, 1));重新取样= 0 (nSamples nRxAnts);为i = 1: nRxAnts重新取样(:,i) =重新取样(eNodeBOutput(:,我),ofdmInfo。SamplingRate圆(sr));结束%进行频率偏移估计和修正流(“……\ n \ nPerforming频差估计”);delta_f = lteFrequencyOffset (enb、重新取样);流(“频率偏移:% 0.3 fhz \ n”,delta_f);重新取样= lteFrequencyCorrect(重新取样,enb delta_f);%找到帧的开始流(“\ nPerforming定时偏移估计…\ n”);抵消= lteDLFrameOffset (enb、重新取样);流(帧的定时偏移量开始:\ n % d样本”,抵消);%将信号与帧的开始重新取样=重新取样(1 +抵消:最终,);% OFDM解调流(“\ nPerforming OFDM解调。\ n \ n”);rxgrid = lteOFDMDemodulate (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而(大小(rxgrid, 2) > 0)流(“% s \ n”、分离器);流(“SIB1解码帧% d \ n”,国防部(enb.NFrame, 1024));流(“% s \ n \ n”、分离器);%重置HARQ缓冲与每个新的SIB1 8帧%的信息可能会有所不同如果(mod (enb.NFrame 8) = = 0)流(“重置HARQ缓冲区。\ n \ n”);decState = [];结束%提取当前子帧rxsubframe = rxgrid (::, 1: L);%进行信道估计[命令,巢]= lteDLChannelEstimate (cec, enb rxsubframe);% PCFICH解调,CFI解码。CFI现在解调%解码使用类似的资源开采和解码功能%显示已经对BCH接待。lteExtractResources用于%从接收到的子帧中提取REs PCFICH对应% rxsubframe命令和信道估计。流(“解码CFI。\ n \ n”);pcfichIndices = ltePCFICHIndices (enb);%得到PCFICH指数[pcfichRx, pcfichHest] = lteExtractResources (pcfichIndices rxsubframe,命令);%解码PCFICHcfiBits = ltePCFICHDecode (pcfichHest, enb, pcfichRx巢);cfi = lteCFIDecode (cfiBits);%得到CFI如果(isfield (enb“CFI”)& & cfi ~ = enb.CFI)释放(pdcchConstDiagram);结束enb。CFI = CFI;流(“解码CFI值:% d \ n \ n ',enb.CFI);%对TDD, PDCCH必须解码盲目地在可能的值% PHICH配置因素m_i (0, 1, 2) TS36.211表6.9 - 1。% m_i值= 0,1和2可以通过配置TDD% uplink-downlink配置1,6 - 0。如果(strcmpi (enb.DuplexMode“TDD”)tddConfigs = (1 6 0);其他的tddConfigs = 0;%不用于FDD,只用于控制while循环结束alldci = {};而(isempty (alldci) & & ~ isempty (tddConfigs))%配置TDD uplink-downlink配置如果(strcmpi (enb.DuplexMode“TDD”enb))。TDDConfig = tddConfigs (1);结束tddConfigs (1) = [];% PDCCH解调。现在PDCCH解调和解码%使用类似的资源开采和解码功能%显示已经对BCH和CFI接待pdcchIndices = ltePDCCHIndices (enb);%得到PDCCH指数[pdcchRx, pdcchHest] = lteExtractResources (pdcchIndices rxsubframe,命令);%解码PDCCH和星座[dciBits, pdcchSymbols] = ltePDCCHDecode (pdcchHest, enb, pdcchRx巢);pdcchConstDiagram (pdcchSymbols);% PDCCH盲目搜索系统(SI)和DCI解码的信息。% PDCCH的LTE工具箱提供了完整的盲目搜索% DCI RNTI指定的消息,在这种情况下,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);% PDCCH搜索DCI结束%如果DCI解码,进行解码PDSCH / DL-SCH为i = 1:元素个数(alldci) dci = alldci {};流(“与SI-RNTI DCI消息:\ n”);disp (dci);%得到DCI PDSCH配置[pdsch, trblklen] = hPDSCHConfiguration (dci, enb 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 (pdsch, enb pdsch.PRBSet);[pdschRx, pdschHest] = lteExtractResources (pdschIndices rxsubframe,命令);%解码PDSCH[dlschBits, pdschSymbols] = ltePDSCHDecode (pdsch, enb pdschRx, pdschHest,巢);%解码DL-SCH HARQ结合软缓冲输入/输出如果~ isempty (decState)流(的重组与先前的传播。\ n \ n”);结束[sib1, crc, decState] = lteDLSCHDecode (pdsch, enb trblklen, dlschBits, decState);%计算PDSCH维生素记录= lteDLSCH (pdsch, enb pdschIndicesInfo。G, sib1);remod = ltePDSCH (enb pdsch记录);[~,refSymbols] = ltePDSCHDecode (pdsch, enb remod);[rmsevm, peakevm] = pdschEVM (refSymbols {1}, pdschSymbols {1});流(“PDSCH RMS维生素:% 0.3 f % % \ n ',rmsevm);流(' PDSCH高峰维生素:% 0.3 f % % \ n \ n”,peakevm);流(“SIB1 CRC: % d \ n”,crc);如果crc = = 0流(“成功SIB1复苏。\ n \ n”);其他的流(“SIB1解码失败。\ n \ n”);结束其他的%显示创建一个PDSCH DCI的配置%消息失败流(“创建从DCI PDSCH配置消息失败。\ n \ n”);结束结束如果(元素个数(alldci) = = 0)%显示DCI解码失败了流(“DCI解码失败。\ n \ n”);结束%更新信道估计阴谋图(channelFigure);冲浪(abs(命令(:,:,1,1)));hSIB1RecoveryExamplePlots (channelFigure);channelFigure.CurrentAxes。XLim =[0大小(命令,2)+ 1);channelFigure.CurrentAxes。YLim =[0大小(命令,1)+ 1];%跳过2帧,并尝试再次SIB1解码,如果我们或终止%剩下不到2帧。如果(大小(rxgrid, 2) > = (L * 20)) rxgrid (:, 1: (L * 20):) = [];%去除2更多的帧其他的rxgrid = [];%小于2帧了结束enb。NFrame =国防部(enb。NFrame + 2, 1024);结束