%配置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;framength = length(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 in dBNocdBm = -98;% Noc in dBm/15kHzNocdBW = NocdBm - 30;dBW/15kHz57中% Noc信噪比= 10^(SNRdB/10);%线性Es/NocEs = SNR*(10^(NocdBW/10));线性Es / RE的百分比FFTLength = 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: frameength +FilterDelaySamples)*SamplePeriod;在帧尾附加足够的样本以弥补延迟过滤器百分比tx波形=时间序列([tx;tx (1: FilterDelaySamples + 1),时间);模拟RF RX的RF块集模型set_param(模型,“LoadInitialState”，“关闭”）;disp (“模拟LTE帧1…”）;sim(模型、时间(结束));在xInitial中保存模型的最终状态，以便下一帧处理xInitial = xFinal;同步接收到的波形Offset = lteDLFrameOffset(rmc，挤压(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波形= timeseries(RepeatFrame,time);%执行Simulin金宝appk RF Blockset测试台disp ([“模拟LTE帧”num2str (n),“……”]);sim(模型、时间(结束));xInitial = xFinal;保存模型状态计算和显示EVM测量evmmeas = hPDSCHEVM(rmc,cec，挤压(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”]);