％配置NDLRB值集以描述要的载波％汇总ndlrb = [50 75 100];％建立了组件载体的数量numcc = numel（ndlrb）;如果numcc <2错误（'请指定多个CC带宽值。'）;结尾

％配置要生成的子帧的数量numsubframes = 10;％CC配置eNB = CELL（1，NUMCC）;为了i = 1：numcc转变ndlrb（i）案件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 {i} .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 = struct;％信道估计配置结构cec.pilotaverage ='用户自定义';％类型的导频符号平均cec.freqwindow = 15;％频率窗口大小cec.timewindow = 15;％时间窗口大小cec.interpype =.'立方体';％2D插值类型cec.interpwindow ='中心';％插值窗口类型cec.Interpwinsize = 1;％插值窗口大小

％定义MHz中的Delta_F1参数，用于标称保护频段和频率％偏移量计算。在下行链路delta_f1中是子载波间距％（TS 36.101，表5.6A-1）delta_f1 = 0.015;MHz中的％在MHz中，LTE子载波间距在MHz中maxbw = hnrbtobandwidth（max（ndlrb））;f_c = zeros（1，numcc）;ccspacing =零（1，numcc-1）;％CC间距％计算标称卫频带，TS 36.101 5.6a-1名词威达带= 0.05 * maxbw-0.5 * delta_f1;％最初假设较低的载波频率是基带f_c（1）= 0;％计算CC间距和载波值为了K = 2：NUMCC CCSPACING（K-1）= HCARRAGGREGATIONCHANNELSPACING（......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（结束）+ f_offset_high;％计算聚合信道带宽，TS 36.101 5.6abw_channel_ca = f_edge_high  -  f_edge_low;fprintf（'bw_channel_ca：％0.4f 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;％显示频带边缘fprintf（'f_edge_low：％0.4f mhz \ n'，f_edge_low）;fprintf（'f_edge_high：％0.4f mhz \ n'，f_edge_high）;fprintf（'f_offset_low：％0.4f mhz \ n'，f_offset_low）;fprintf（'f_offset_high：％0.4f mhz \ n'，f_offset_high）;％显示载波频率fprintf（'\ n'）;为了i = 1：numcc fprintf（'组分载体％d：\ n'，一世）;fprintf（'fc：％0.4f mhz \ n'，f_c（i））;结尾

85％的百分比利用率Bwfraction = 0.85;％计算分量载波的采样率CCSR =零（1，NUMCC）;为了i = 1：numcc info = lteofdminfo（eNB {i}）;ccsr（i）= info.samplingrate;结尾％计算聚合信号的整体采样率％计算过采样比率（最大的BW CC）以确保％信号占总带宽的85％（BWFraction）OSR =（BW_CHANNEL_CA / BWFRACTECT）/（MAX（CCSR）/ 1E6）;％为了简化重采样操作，选择过采样率％是2的力量：计算OSR上方两个的下一个功率osr = 2 ^ ceil（log2（osr））;sr = osr * max（ccsr）;fprintf（'\ noutput采样率：％0.4f ms / s \ n \ n'，SR / 1E6）;％计算组件载体的单个过采样因子osrs = sr./ccsr;

％生成分量载波tx = cell（1，numcc）;为了i = 1：numcc tx {i} = ltermcdltool（eNB {i}，randi（[0 1]，1000,1））;tx {i} =重组（tx {i}，osrs（i），1）/ osrs（i）;tx {i} = hcArrieraggregationmodulate（tx {i}，sr，f_c（i）* 1e6）;结尾％叠加分量载体波形= Tx {1};为了i = 2：numcc waveform = waveform + tx {i};结尾

Specplot = HcArrieraggregationPlotsPectrum（波形，SR，......'载波聚合的功率范围'，{'信号谱'}）;

％选择CC解调ccofinterest = 1;如果ccofinterest> numcc ||ccofinterest <= 0错误（'无法解调CC号％d，有％d ccs \ n'那......CCOFINTEREST，NUMCC）;结尾％定义下采样过滤器顺序FilterOrder = 201;所有CC的％预先烫伤通带和阻带值FiRPassBandVec =（0.18 * NDLRB + DELTA_F1）/ 2 /（SR / 1E6 / 2）;firstopbandvec = hcarrieraggregationstopband（ccspacing，ndlrb，sr）;％定义了感兴趣的载体的通带和停车FirPassBand = FirPassBandVec（CCOFInterest）;firstopband = firstopbandvec（ccofinterest）;％提取物和解码兴趣的CCFPRINTF（1，'提取CC号％d：\ n'，ccofinterest）;％PAD信号用零点考虑过滤器瞬态响应％ 长度波形= [波形;零（FilterOrder + 1，尺寸（波形，2））];％解调感兴趣的载体DemodulatedCC =......hcArrieraggregationmodulate（波形，SR，-F_C（CCOFInterest）* 1E6）;如果过滤器太窄，则％下采样在两个阶段完成。这％简化过滤器设计要求。如果这是一个初始的情况施加％下采样系数为4。你可能想考虑一个不同的如果产生的质量，在这个初始阶段中的％过滤器设计％信号是不可接受的。如果（FirstopBand <0.1）百分比下样品4SRC = 4;DemodulatedCC =重组（DemodulatedCC，1，SRC）;％更新传递和停机带到首先进行下式采样FirePassBand = 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，'分析'那'弗里克'）;

％过滤信号以提取感兴趣的组件rxwaveform = firfilter（demodulatedcc）;％绘制解调和过滤的信号谱filteredspecplot =.......hcarrieraggregationplotspectrum（[demodulatedcc，rxwaveform]，sr，......'解调和过滤波形功率谱'那......{'载波聚合信号''过滤信号'}）;

rxwaveform = downsample（rxwaveform，osrs（ccofinterest）/ src）;

感兴趣的CC的％参数eNodeB = eNB {CCOFINTEREST};PDSCH = eNodeB.PDSCH;％同步接收波形offset = ltedlframeoffset（eNodeb，rxwaveform）;rxwaveform = rxwaveform（1 + offset：结束，:);

EVM测量的％信道估计结构cecevm = cec;cecevm.pilotaverage =.'testevm';[EVMMEAS，PLOTS] = HPDSCHEVM（eNB {CCOFINTEREST}，CECEVM，RXWAVEFORM）;％OFDM解调rxgrid = lteofdmdemodulate（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，噪声] = LTedLchanneLeStimate（......eNodeB，CEC，RXSUBFRAME）;％执行贬低，层解调，解调和解调使用频道估计的收到数据的％[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通过了']）;结尾结尾