% Configuration TS 36.101 25 REs (5 MHz)， 64-QAM, full allocationrmc = lteRMCDL (“R.6”）;rmc。OCNGPDSCHEnable =“上”；%使用固定PDSCH数据创建eNodeB传输rng (2);固定随机种子(任意)data = randi([0 1]， sum(rmc.PDSCH.TrBlkSizes)，1); / /将rmc.PDSCH.TrBlkSizes赋值%生成1帧，重复模拟总共N帧[tx， ~， info] = lteRMCDLTool(rmc, data);% 1帧%计算采样周期和帧的长度。SamplePeriod=1/info.SamplingRate；帧长度=长度（tx）；

带限插值滤波器FiltOrd = 32;h = firpm(FiltOrd，[0 2.25e6*2*SamplePeriod 2.7e6*2*SamplePeriod 1]，[1 1 0 0]);FilterDelaySamples = FiltOrd / 2;滤波器组延迟%%传播模型SNRdB = 57;%Es/Noc（单位：dB）NocdBm = -98;% Noc在dBm/15kHzNocdBW = NocdBm - 30;% Noc在dBW/15kHz57信噪比=10^（SNRdB/10）；%线性Es / NocEs =信噪比* (10 ^ (NocdBW / 10));每RE线性Es %FFTLength = info.Nfft;SymbolPower = Es /双(FFTLength);%模拟帧数N>=1N = 3;%为N-1帧的结果预分配向量% EVM在第一帧不测量，以避免瞬变效应evmpeak = 0 (n - 1, 1);%预分配结果evmrms = 0 (n - 1, 1);%预分配结果

%指定并打开Simulink模型金宝app模型=“RFLTEReceiverModel”；disp(“开始仿真软件”金宝app); 开放式系统（模型）；

%生成射频接收机测试数据时间= (0:FrameLength + FilterDelaySamples) * SamplePeriod;在帧的末尾添加足够的采样以补偿延迟过滤器百分比txWaveform = timeseries ([tx;tx (1: FilterDelaySamples + 1),时间);%模拟RF RX的RF块集模型set_param(模型,“LoadInitialState”，“关闭”）;disp (“模拟LTE第一帧……”）;sim(模型、时间(结束));%将模型的最终状态保存在xInitial中以供下一帧处理xInitial = xFinal;与接收的波形同步抵消= lteDLFrameOffset (rmc中,挤压(rx),“TestEVM”）;%在这种情况下，Offset = FilterDelaySamples，因此如下所示%帧不需要同步

前一帧执行后的加载状态。因为我们在重复%同一帧每次执行帧后，模型状态将相同。set_param(模型,“LoadInitialState”，“上”，“InitialState”，“xInitial”）;%修改输入向量以考虑带宽限制的延迟%过滤器RepeatFrame = (tx (FilterDelaySamples + 1:结束);tx (1: FilterDelaySamples + 1)];EVMalg。EnablePlotting =“关闭”；cec。PilotAverage =“TestEVM”；为n = 2: n%为所有剩余帧%生成数据时间=（（n-1）*帧长+（0:FrameLength）+过滤器延迟采样）*采样周期；txWaveform=时间序列（RepeatFrame，time）；%执行Simulin金宝appk RF Blockset测试disp ([“模拟LTE帧”num2str (n),' ...']）;sim(模型、时间(结束));xInitial = xFinal;保存模型状态计算和显示EVM测量值evmmeas = hPDSCHEVM (rmc cec,挤压(rx) EVMalg);evmpeak (n - 1) = evmmeas.Peak;evmrms (n - 1) = evmmeas.RMS;结束

图;情节(100 * evmpeak (2: N),“啊——”)标题(“维生素与峰值%”）;包含(“帧数”）;图;情节(100 * evmrms (2: N),“啊——”); 头衔(“EVM RMS%”）;包含(“帧数”）;

bdclose（模型）；clear（模型，“_acc”]）;