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

%限带插补滤波器FiltOrd = 32;h = firpm(FiltOrd，[0 2.25e6*2*SamplePeriod 1]，[1 1 0 0]);FilterDelaySamples = FiltOrd / 2;%滤波器组延迟%传播模型SNRdB = 57;dB中的% Es/NocNocdBm = -98;% Noc在dBm/15kHznodbw = nodbm - 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);open_system(模型);

%生成射频接收机测试数据时间= (0:FrameLength + FilterDelaySamples) * SamplePeriod;%在帧的末尾附加足够的样本以补偿延迟过滤器的%txWaveform = timeseries ([tx;tx (1: FilterDelaySamples + 1),时间);%模拟RF RX的RF Blockset模型set_param(模型,“LoadInitialState”,“关闭”);disp (“模拟LTE帧1……”);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%表示所有剩余帧%生成数据time = ((n-1)*FrameLength+(0:FrameLength) + FilterDelaySamples)*SamplePeriod;txWaveform = timeseries (RepeatFrame、时间);%执行Simulin金宝appk RF Blockset testbenchdisp ([“模拟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),“啊——”);标题(“维生素与RMS %”);包含(的帧数的);

bdclose(模型);明确([模型,“_acc”]);