%配置一组NDLRB值，以描述要使用的载波%聚合NDLRB = [50 75 100];％建立了组件载体的数量numCC=numel（NDLRB）；如果numCC < 2错误('请指定多个CC带宽值。')；结尾

%配置要生成的子帧数numSubframes = 10;% CC配置enb =细胞(1、numCC);为了i=1:numCC开关NDLRB(我)案例6 enb{i}=lteRMCDL('r.4')；案例15 enb{i} = lteRMCDL('r.5')；案例25 enb{i} = lteRMCDL('r.6')；案例50 enb{i} = lteRMCDL('r.7')；案例75 enb{i}=lteRMCDL(“R.8”)；案例100 enb{i} = lteRMCDL('r.9')；否则fprintf（'没有有效的资源块数量:%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 =“立方”；二维插值型电子窗口=“中心”；插值窗口类型cec。InterpWinSize = 1;插值窗口大小

%为额定保护频带和频率定义MHz中的Delta_f1参数％偏移量计算。在下行链路delta_f1中是子载波间距％（TS 36.101，表5.6A-1）delta_f1 = 0.015;%以MHz为单位，LTE子载波间隔以MHz为单位maxBW = hNRBToBandwidth (max (NDLRB));numCC F_c = 0 (1);numCC-1 ccSpacing = 0 (1);% CC间距计算公称保护带，TS 36.101 5.6A-1名词威达带= 0.05 * maxbw-0.5 * delta_f1;%最初假设基带的载波频率较低F_c (1) = 0;％计算CC间距和载波值为了numCC ccSpacing(k-1) = hCarrierAggregationChannelSpacing(...eNB {K-1} .BandWidth，eNB {k} .BandWidth）;f_c（k）= f_c（k-1）+ ccspacing（k-1）;结尾%计算较低和较高频率偏移，TS 36.101 5.6Af_offset_low =（0.18 * ndlrb（1）+ delta_f1）/ 2 +名义根带;f_offset_high =（0.18 * ndlrb（end）+ delta_f1）/ 2 +名称是;%计算较低和较高频率边缘，TS 36.101 5.6AF_edge_low = F_c(1) - F_offset_low;F_edge_high = F_c(end) + F_offset_high;％计算聚合信道带宽，TS 36.101 5.6aBW\u信道\u CA=F\u边缘\u高-F\u边缘\u低；fprintf(“BW_channel_CA: f % 0.4 MHz \ n”，BW_channel_CA）；计算移动到中心基带shiftTocenter = -1 *（bw_channel_ca / 2 + f_edge_low）;％聚合带宽以基带为中心F_c=F_c+移位中心；F_边缘低=F_c（1）-F_偏移低；F_边缘高=F_c（末端）+F_偏移高；％显示频带边缘fprintf（'频率边缘低：%0.4f MHz\n', F_edge_low);fprintf（“F_edge_high: f % 0.4 MHz \ n”, F_edge_high);fprintf（“F_offset_low: f % 0.4 MHz \ n”, F_offset_low);fprintf（“F_offset_high: f % 0.4 MHz \ n”，f_offset_high）;显示载波频率fprintf（' \ n ')；为了i = 1:numCC fprintf(“分量载波% d: \ n”，i）；fprintf('Fc:%0.4f MHz\n'，F_c（i））；结尾

%带宽利用率为85%bw分数=0.85；％计算分量载波的采样率numCC CCSR = 0 (1);为了i = 1：numcc info = lteofdminfo（eNB {i}）;ccsr（i）= info.samplingrate;结尾%计算聚合信号的总体采样率计算过采样比率(最大BW CC)，以确保%信号占总带宽的85%（BW分数）OSR =（BW_CHANNEL_CA / BWFRACTECT）/（MAX（CCSR）/ 1E6）;为了简化重采样操作，选择一个过采样比%是2的幂：计算OSR上方2的下一个幂osr = 2 ^ ceil（log2（osr））;sr = osr * max（ccsr）;fprintf（'\nOutput采样率:%0.4f Ms/s\n\n'，SR/1e6）；％计算组件载体的单个过采样因子老osr = / CCSR;

生成组件载波tx = cell（1，numcc）;为了i = 1:numCC tx{i} = lteRMCDLTool(enb{i}，randi([0 1]，1000,1));tx{我}=重新取样(tx {}, osr(我),1)/ osr(我);tx{我}= hCarrierAggregationModulate (tx{我}、SR、F_c(我)* 1 e6);结尾%叠加组件载体波形= Tx {1};为了i = 2:numCC波形=波形+ tx{i};结尾

specPlot = hCarrierAggregationPlotSpectrum(波形、SR、...“载波聚合的功率谱”, {信号频谱的});

%选择CC进行解调CCofInterest=1；如果CCofInterest>numCC | | CCofInterest<=0错误(“无法解调CC号%d，有%d CC \n”那...CCOFINTEREST，NUMCC）;结尾％定义下采样过滤器顺序过滤器顺序=201；预计算所有CC的通频带和阻频带值firPassbandVec=（0.18*NDLRB+Delta_f1）/2/（SR/1e6/2）；firStopbandVec=hCarrierAggregationTopBand（ccSpacing，NDLRB，SR）；定义感兴趣载波的通频带和阻频带firPassband = firPassbandVec (CCofInterest);firStopband = firStopbandVec (CCofInterest);％提取物和解码兴趣的CCfprintf（1，'正在提取抄送号码%d:\n'，cFinterest）；%填充零信号以考虑滤波器瞬态响应%的长度波形=[波形；零（滤波器顺序+1，大小（波形，2））；%解调感兴趣的载波DemodulatedCC =...hcArrieraggregationmodulate（波形，SR，-F_C（CCOFInterest）* 1E6）;如果滤镜太窄，下采样分两阶段进行。这%简化过滤器设计要求。如果是这样的话，首字母缩写%的下采样系数为4。滤光片的设计在这一初始阶段若质量不佳%信号不可接受。如果(firStopband < 0.1)百分比下样品4src = 4;DemodulatedCC =重组（DemodulatedCC，1，SRC）;％更新传递和停机带到首先进行下式采样firPassband=firPassband*SRC；firStopband=firStopband*SRC；其他的%不预滤器SRC = 1;结尾%设计低通滤波器以滤除感兴趣的CCFredges = [0 FirepassBand FirstopBand 1];firfilter = dsp.firfilter;firfilter.numerator = FILPM（FilterOrder，Fredges，[1 1 0 0]）;

fvtool (firFilter“分析”那'弗里克')；

对信号进行滤波，提取感兴趣的分量rx波形=滤波器（解调的CC）；%绘制解调滤波信号频谱过滤样图=...hCarrierAggregationPlotSpectrum ([demodulatedCC rxWaveform],老...“解调和滤波波形功率谱”那...{'载波聚合信号'“过滤信号”});

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

%利息CC参数eNodeB=enb{CCofInterest}；PDSCH=eNodeB.PDSCH；%同步接收波形offset = ltedlframeoffset（eNodeb，rxwaveform）;rxwaveform = rxwaveform（1 + offset：结束，:);

%用于EVM测量的信道估计结构cecEVM = cec;cecEVM。PilotAverage =“TestEVM”；[evmmeas, plot] = hPDSCHEVM(enb{CCofInterest}，cecEVM, rx波形);% OFDM解调rxGrid=LTEOFDMDDemodulate（eNodeB，rxWaveform）；%获取每个子帧接收的子帧和OFDM符号的数量dims = lteofdminfo（eNodeB）;samplespserubframe = dims.samplingrate / 1000;nrxsubframes = bloor（size（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，noiseEst]=lteDLChannelEstimate(...eNodeB、cec、rxSubframe）；％执行贬低，层解调，解调和解调%在接收数据上使用信道估计[rxEncodedBits，pdschsymbols]=ltePDSCHDecode（eNodeB，PDSCH，...rxSubframe*（10^（-PDSCH.Rho/20）），estChannelGrid，noiseEst）；%解码下行链路共享信道（DL-SCH）[十进制位，crc]=...lteDLSCHDecode (eNodeB PDSCH、trBlkSize rxEncodedBits {1});如果crc显示([“子帧”num2str (n)'：CRC失败']）;其他的disp([“子帧”num2str (n): CRC传递的]）;结尾结尾