%配置组NDLRB值来描述运营商%聚合NDLRB = (50 75 100);%建立组件的航空公司的数量numCC =元素个数(NDLRB);如果numCC < 2错误(“请指定超过一个CC带宽值。);结束

%配置子帧生成的数量numSubframes = 10;% CC配置enb =细胞(1、numCC);为我= 1:numCC开关NDLRB(我)情况下6 enb{我}= lteRMCDL (“R.4”);情况下15 enb{我}= lteRMCDL (“R.5”);情况下25 enb{我}= lteRMCDL (“R.6”);情况下50 enb{我}= lteRMCDL (“R.7”);情况下75年enb{我}= lteRMCDL (“R.8”);情况下100年enb{我}= lteRMCDL (“R.9”);否则流(”不是一个有效的资源块:% d \ n ',…NDLRB (i));返回结束enb {}。NDLRB= NDLRB(i); enb{i}.Bandwidth = hNRBToBandwidth(NDLRB(i)); enb{i}.TotSubframes = numSubframes; enb{i}.PDSCH.PRBSet = (0:enb{i}.NDLRB-1).'; enb{i}.PDSCH.RVSeq = 0; enb{i}.NCellID = 10;结束

cec =结构;%信道估计配置结构cec。PilotAverage =“UserDefined”;%的试点,平均象征cec。FreqWindow = 15;%频率窗口大小cec。TimeWindow = 15;%的时间窗口大小cec。InterpType =“立方”;% 2 d插值类型cec。InterpWindow =“中心”;%插值窗口类型cec。InterpWinSize = 1;%插值窗口大小

%为名义保护带和定义Delta_f1参数MHz的频率%偏移量的计算。在下行Delta_f1副载波间距% (TS 36.101,表5.6 a - 1)Delta_f1 = 0.015;%在兆赫,LTE副载波间距MHzmaxBW = hNRBToBandwidth (max (NDLRB));numCC F_c = 0 (1);numCC-1 ccSpacing = 0 (1);% CC间距%计算名义保护带,TS 36.101 - 5.6 a - 1nominalGuardBand = 0.05 * maxbw - 0.5 * Delta_f1;%最初假设降低载波频率在基带F_c (1) = 0;%计算CC间距和载体值为k = 2: numCC ccSpacing (k - 1) = hCarrierAggregationChannelSpacing (…enb {k - 1}。带宽,enb {k} .Bandwidth);F_c (k) = F_c (k - 1) + ccSpacing (k - 1);结束%计算低和更高的频率偏移,TS 36.101 - 5.6F_offset_low = (0.18 * NDLRB (1) + Delta_f1) / 2 + nominalGuardBand;F_offset_high = (0.18 * NDLRB(结束)+ Delta_f1) / 2 + nominalGuardBand;%计算低和高频率边缘,TS 36.101 - 5.6F_edge_low = F_c (1) - F_offset_low;F_edge_high = F_c(结束)+ F_offset_high;%计算聚合通道带宽,TS 36.101 - 5.6BW_channel_CA = F_edge_high - F_edge_low;流(“BW_channel_CA: f % 0.4 MHz \ n”,BW_channel_CA);%计算转向中心基带shiftToCenter = 1 * (BW_channel_CA / 2 + F_edge_low);集中在基带%聚合带宽F_c = F_c + shiftToCenter;F_edge_low = F_c (1) - F_offset_low;F_edge_high = F_c(结束)+ F_offset_high;%显示频带边缘流(“F_edge_low: f % 0.4 MHz \ n”,F_edge_low);流(“F_edge_high: f % 0.4 MHz \ n”,F_edge_high);流(“F_offset_low: f % 0.4 MHz \ n”,F_offset_low);流(“F_offset_high: f % 0.4 MHz \ n”,F_offset_high);%显示载波频率流(' \ n ');为i = 1: numCC流(“分量载波% d: \ n”,我);流(舰队指挥官:f % 0.4 MHz \ n的F_c(我));结束

%的带宽利用率为85%bwfraction = 0.85;%计算组件运营商的采样率numCC CCSR = 0 (1);为i = 1: numCC信息= lteOFDMInfo (enb{我});CCSR (i) = info.SamplingRate;结束%计算聚合整体采样率信号%计算过采样率(最大的BW CC)确保占总数的%信号占85% (bwfraction)带宽OSR = (BW_channel_CA / bwfraction) / (max (CCSR) / 1 e6);%,以简化重采样操作选择的过采样率%是2的乘方:计算下一两个以上OSR的力量OSR = 2 ^装天花板(log2 (OSR));SR = OSR *马克斯(CCSR);流(女士' \ nOutput采样率:% 0.4 f / s \ n \ n”1、SR / e6);%计算单个组件的过采样因素运营商老osr = / CCSR;

%生成组件运营商tx =细胞(1、numCC);为i = 1: numCC tx{我}= lteRMCDLTool (enb{},兰迪([0,1],1000,1));tx{我}=重新取样(tx {}, osr(我),1)/ osr(我);结束%聚合所有CC. comm.MultibandCombiner适用频率偏移量%并添加一起产生的信号。sigAggregator = comm.MultibandCombiner (InputSampleRate =老,…e6 FrequencyOffsets = F_c * 1,…OutputSampleRateSource =“属性”,…OutputSampleRate = SR);波形= sigAggregator ((tx {:}));

specPlot = hPlotSpectrum(波形、SR、…载波聚合的功率谱,{信号频谱的});

%选择CC解调CCofInterest = 1;如果CCofInterest > numCC | | CCofInterest < = 0错误(”无法解调CC % d,有% d CCs \ n ',…CCofInterest numCC);结束%定义将采样滤波器秩序filterOrder = 201;% Precalculate通带和阻带所有CC值firPassbandVec = (0.18 * NDLRB + Delta_f1) / 2 / (SR / 1 e6/2);firStopbandVec = hCarrierAggregationStopband (ccSpacing NDLRB, SR);%定义通带和阻带载体的兴趣firPassband = firPassbandVec (CCofInterest);firStopband = firStopbandVec (CCofInterest);%提取和解码CC的兴趣流(1,“提取CC % d: \ n”,CCofInterest);%垫与零信号,考虑滤波器的瞬态响应%的长度波形=[波形;0 (filterOrder + 1、大小(波形,2)));%解调载波的兴趣demodulatedCC =…hCarrierAggregationModulate(波形、SR、-F_c (CCofInterest) * 1 e6);%将采样是在两个阶段完成,如果过滤器太窄。这%简化了滤波器设计要求。如果这是一个初始的情况%将采样因子4。你可能想要考虑不同%滤波器设计在这个初期如果结果的质量%信号是不能接受的。如果(firStopband < 0.1)% Down-sample除以4SRC = 4;demodulatedCC =重新取样(demodulatedCC 1 SRC);%更新传递和阻带考虑首先将采样firPassband = firPassband * SRC;firStopband = firStopband * SRC;其他的%不预滤器SRC = 1;结束%设计低通滤波器过滤掉CC的兴趣frEdges = [0 firPassband firStopband 1];firFilter = dsp.FIRFilter;firFilter。分子= firpm (filterOrder、frEdges [1 1 0 0]);

fvtool (firFilter“分析”,“频率”);

%过滤信号提取感兴趣的组件rxWaveform = firFilter (demodulatedCC);%画出解调和过滤后的光谱信号filteredSpecPlot =…hPlotSpectrum ([demodulatedCC rxWaveform],老…“解调和过滤后的波形功率谱”,…{“载波聚合信号”“过滤信号”});

rxWaveform = downsample (rxWaveform osr (CCofInterest) / SRC);

% CC感兴趣的参数eNodeB = enb {CCofInterest};PDSCH = eNodeB.PDSCH;%同步接收波形抵消= lteDLFrameOffset (eNodeB rxWaveform);rxWaveform = rxWaveform(1 +抵消:最终,);

%为维生素与信道估计结构测量cecEVM = cec;cecEVM。PilotAverage =“TestEVM”;[evmmeas,情节]= hPDSCHEVM (enb {CCofInterest}, cecEVM, rxWaveform);% OFDM解调rxGrid = lteOFDMDemodulate (eNodeB rxWaveform);%的数量获得每个子帧子帧和OFDM符号= lteOFDMInfo dim (eNodeB);samplesPerSubframe = dims.SamplingRate / 1000;nRxSubframes =地板(大小(rxWaveform 1) / samplesPerSubframe);eNodeB。TotSubframes = 1;resGridSize = lteResourceGridSize (eNodeB);L = resGridSize (2);disp (传播子帧,并检查CRC解码。);为n = 0: nRxSubframes-1%提取子帧rxSubframe = rxGrid (:, (1: L) + (n * L):);%当前子帧传输块大小eNodeB。NSubframe = n;trBlkSize = PDSCH.TrBlkSizes (n + 1);%信道估计[estChannelGrid,噪音]= lteDLChannelEstimate (…eNodeB、cec rxSubframe);%执行deprecoding层demapping解调和descrambling%接收的数据使用信道估计[rxEncodedBits, pdschsymbols] = ltePDSCHDecode (eNodeB PDSCH,…rxSubframe * (10 ^ (-PDSCH.Rho / 20)), estChannelGrid,噪音);%下行共享信道解码(DL-SCH)(decbits, crc) =…lteDLSCHDecode (eNodeB PDSCH、trBlkSize rxEncodedBits {1});如果crc disp ([“子帧”num2str (n)“:CRC失败”]);其他的disp ([“子帧”num2str (n):CRC传递的]);结束结束