%配置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 2.7e6*2*SamplePeriod 1]，[1 1 0 0]);FilterDelaySamples = FiltOrd/2;%滤波器组延迟%传播模型SNRdB = 57;数据库中% Es/NocNocdBm = -98;% Noc在dBm/15kHzNocdBW = NocdBm - 30;% Noc在dBW/15kHz57SNR = 10^(SNRdB/10);线性Es/NocEs = SNR*(10^(NocdBW/10));%线性Es / REFFTLength = info.Nfft;SymbolPower = Es/double(FFTLength);%模拟帧数N>=1N = 3;为N-1帧的结果预分配向量% EVM不在第一帧测量，以避免瞬态效应evmpeak = 0 (n -1,1);%结果预分配evmrms = 0 (n -1,1);%结果预分配

指定并打开Simulink模型金宝app模型=“RFLTEReceiverModel”；disp (“开始仿真软件”金宝app）;open_system(模型);

生成射频接收器的测试数据time = (0:FrameLength+FilterDelaySamples)*SamplePeriod;在帧的末尾添加足够的样本以弥补延迟。过滤器的%tx波形=时间序列([tx;tx (1: FilterDelaySamples + 1),时间);模拟RF RX的RF块集模型set_param(模型,“LoadInitialState”，“关闭”）;disp (“模拟LTE帧1…”）;sim(模型、时间(结束));将模型的最终状态保存在xInitial中以备下一个帧处理xInitial = xFinal;%同步接收波形Offset = lteDLFrameOffset(rmc,squeeze(rx)，“TestEVM”）;在这种情况下，Offset = FilterDelaySamples因此如下%帧不需要同步

执行前一帧后的加载状态。因为我们在重复%相同帧模型状态在每帧执行后都是相同的。set_param(模型,“LoadInitialState”，“上”，“InitialState”，“xInitial”）;%修改输入矢量，以考虑带宽限制的延迟%过滤器RepeatFrame = [tx(FilterDelaySamples+1:end);tx (1: FilterDelaySamples + 1)];EVMalg。EnablePlotting =“关闭”；cec。PilotAverage =“TestEVM”；为n = 2: n%为所有剩余帧%生成数据time = ((n-1)*FrameLength+(0:FrameLength) + FilterDelaySamples)*SamplePeriod;tx波形=时间序列(RepeatFrame,time);执行Simulink金宝app RF Blockset测试台disp ([“模拟LTE帧”num2str (n),“……”]);sim(模型、时间(结束));xInitial = xFinal;%保存模型状态%计算和显示EVM测量evmmeas = hPDSCHEVM(rmc,cec,squeeze(rx)，EVMalg);evmpeak(n-1) = evmmeasure . peak;evmrms(n-1) = evmmeasure . rms;结束

图;情节(100 * evmpeak (2: N),“啊——”)标题('EVM峰值%'）;包含('帧数'）;图;情节(100 * evmrms (2: N),“啊——”）;标题(' evm RMS %'）;包含('帧数'）;

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