This example demonstrates how to model and test an LTE RF receiver using LTE System Toolbox™ and RF Blockset™.
The figure below shows the main parts of this example. An LTE waveform is generated using the LTE System Toolbox. This waveform is filtered and transmitted through a propagation channel before feeding it to the RF receiver model implemented with RF Blockset. This model is based on commercially available parts. EVM figures are then provided for the output of the RF receiver.
This example is implemented using MATLAB® and Simulink®, which interact at runtime. The functional partition is shown in the figure below
The MATLAB script implements the simulation test-bench, and the Simulink model is the device under test (DUT). LTE frames are streamed between the test-bench and the DUT.
In this section we generate the LTE waveform using the LTE System Toolbox. We use the reference measurement channel (RMC) R.6 as defined in TS 36.101 [ 1 ]. This RMC specifies a 25 resource elements (REs) bandwidth, equivalent to 5 MHz. A 64 QAM modulation is used. All REs are allocated. Additionally, OCNG noise is enabled in unused REs.
Only one frame is generated. This frame will then be repeated a number of times to perform the EVM measurements.
% Configuration TS 36.101 25 REs (5 MHz), 64-QAM, full allocation rmc = lteRMCDL('R.6'); rmc.OCNGPDSCHEnable = 'On'; % Create eNodeB transmission with fixed PDSCH data rng(2); % Fixed random seed (arbitrary) data = randi([0 1], sum(rmc.PDSCH.TrBlkSizes),1); % Generate 1 frame, to be repeated to simulate a total of N frames [tx, ~, info] = lteRMCDLTool(rmc, data); % 1 frame % Calculate the sampling period and the length of the frame. SamplePeriod = 1/info.SamplingRate; FrameLength = length(tx);
This section initializes some of the simulation components:
Band limiting filter: design the filter coefficients, which will be used by the Simulink model. The filter has order 32, with passband frequency equal to 2.25 MHz, and stopband frequency equal to 2.7 MHz.
SNR and signal energy
Number of frames: this is the number of times the generated frame is repeated
Preallocate result vectors
% Band limiting interpolation filter FiltOrd = 32; h = firpm(FiltOrd,[0 2.25e6*2*SamplePeriod 2.7e6*2*SamplePeriod 1],[1 1 0 0]); FilterDelaySamples = FiltOrd/2; % filter group delay % Propagation model SNRdB = 57; % Es/Noc in dB NocdBm = -98; % Noc in dBm/15kHz NocdBW = NocdBm - 30; % Noc in dBW/15kHz57 SNR = 10^(SNRdB/10); % linear Es/Noc Es = SNR*(10^(NocdBW/10)); % linear Es per RE FFTLength = info.Nfft; SymbolPower = Es/double(FFTLength); % Number of simulation frames N>=1 N = 3; % Preallocate vectors for results for N-1 frames % EVM is not measured in the first frame to avoid transient effects evmpeak = zeros(N-1,1); % Preallocation for results evmrms = zeros(N-1,1); % Preallocation for results
This section loads the Simulink model shown below. This model includes the following components:
Reading the LTE waveform and the sampling period from the workspace
Channel model: this includes free space path loss and AWGN
RF receiver, including direct conversion demodulator
ADC and DC offset cancellation
Save results to workspace
% Specify and open Simulink model model = 'RFLTEReceiver'; disp('Starting Simulink'); open_system(model);
Since we are running in accelerated mode we need a temporary directory to store some generated files.
originalDir = pwd; tempDir = tempname; mkdir(tempDir); cd(tempDir);
The RF receiver model includes the elements shown below
The RF demodulator includes the following components and is shown below:
Local oscillator (LO) and phase noise model
Phase shift for I and Q components generator
This section simulates the specified number of frames. This is done in two stages:
Simulate the first frame
Simulate the rest of the frames in a loop
The reason for splitting the processing in these two stages is to simplify the code. During the processing of the first frame we need to take into account the delay of the band limiting filter. This is not the case for subsequent frames, since the filter state is maintained between frames. Therefore, the length of the first frame has to be increased slightly to take into account the delay introduced by the filter.
Simulate First LTE Frame
As mentioned for the first simulated frame we need to increase the length of the signal fed to the Simulink model to compensate for the delay introduced by the filter. Next we launch the simulation of the Simulink model without loading any initial state. After processing the first frame with the Simulink model, its state (
xFinal) is stored and assigned to
xInitial for the next time the model is run.
The output of the Simulink model is stored in the variable
rx, which is available in the workspace. Any delays introduced to this signal are removed after performing synchronization. The EVM is measured on the resulting waveform.
% Generate test data for RF receiver time = (0:FrameLength+FilterDelaySamples)*SamplePeriod; % Append to the end of the frame enough samples to compensate for the delay % of the filter txWaveform = timeseries([tx; tx(1:FilterDelaySamples+1)],time); % Simulate RF Blockset model of RF RX set_param(model, 'LoadInitialState', 'off'); disp('Simulating'); sim(model, time(end)); % Save the final sate of the model in xInitial for next frame processing xInitial = xFinal; % Synchronize to received waveform Offset = lteDLFrameOffset(rmc,squeeze(rx),'TestEVM'); % In this case Offset = FilterDelaySamples therefore the following % frames do not require synchronization
Simulate Successive LTE Frames
Now the rest of the frames can be simulated. First, the model state is set using the value stored in
xInitial at the output of the previous iteration.
% Load state after execution of previous frame. Since we are repeating the % same frame the model state will be the same after every frame execution. set_param(model, 'LoadInitialState', 'on', 'InitialState','xInitial'); % Modify input vector to take into account the delay of the bandlimiting % filter RepeatFrame = [tx(FilterDelaySamples+1:end); tx(1:FilterDelaySamples+1)]; EVMalg.EnablePlotting = 'Off'; cec.PilotAverage = 'TestEVM'; for n = 2:N % for all remaining frames % Generate data time = ( (n-1)*FrameLength+(0:FrameLength) + FilterDelaySamples)*SamplePeriod; txWaveform = timeseries(RepeatFrame,time); % Execute Simulink RF Blockset testbench sim(model, time(end)); xInitial = xFinal; % Save model state % Compute and display EVM measurements evmmeas = hPDSCHEVM(rmc,cec,squeeze(rx),EVMalg); evmpeak(n-1) = evmmeas.Peak; evmrms(n-1) = evmmeas.RMS; end
Low edge EVM, subframe 0: 2.919% High edge EVM, subframe 0: 2.922% Low edge EVM, subframe 1: 2.632% High edge EVM, subframe 1: 2.647% Low edge EVM, subframe 2: 3.016% High edge EVM, subframe 2: 3.023% Low edge EVM, subframe 3: 3.451% High edge EVM, subframe 3: 3.446% Low edge EVM, subframe 4: 2.921% High edge EVM, subframe 4: 2.933% Low edge EVM, subframe 6: 3.054% High edge EVM, subframe 6: 3.048% Low edge EVM, subframe 7: 3.188% High edge EVM, subframe 7: 3.183% Low edge EVM, subframe 8: 3.115% High edge EVM, subframe 8: 3.107% Low edge EVM, subframe 9: 3.065% High edge EVM, subframe 9: 3.066% Averaged low edge EVM, frame 0: 3.049% Averaged high edge EVM, frame 0: 3.050% Averaged EVM frame 0: 3.050% Averaged overall EVM: 3.050% Low edge EVM, subframe 0: 2.931% High edge EVM, subframe 0: 2.936% Low edge EVM, subframe 1: 2.649% High edge EVM, subframe 1: 2.663% Low edge EVM, subframe 2: 3.003% High edge EVM, subframe 2: 3.012% Low edge EVM, subframe 3: 3.430% High edge EVM, subframe 3: 3.424% Low edge EVM, subframe 4: 2.922% High edge EVM, subframe 4: 2.934% Low edge EVM, subframe 6: 3.057% High edge EVM, subframe 6: 3.051% Low edge EVM, subframe 7: 3.189% High edge EVM, subframe 7: 3.183% Low edge EVM, subframe 8: 3.108% High edge EVM, subframe 8: 3.099% Low edge EVM, subframe 9: 3.050% High edge EVM, subframe 9: 3.052% Averaged low edge EVM, frame 0: 3.046% Averaged high edge EVM, frame 0: 3.047% Averaged EVM frame 0: 3.047% Averaged overall EVM: 3.047%
This section plots the measured peak and RMS EVM for each simulated frame.
figure; plot((2:N), 100*evmpeak,'o-') title('EVM peak %'); xlabel('Number of frames'); figure; plot((2:N), 100*evmrms,'o-'); title('EVM RMS %'); xlabel('Number of frames');
Close the Simulink model and remove the generated files.
bdclose(model); clear([model,'_acc']); cd(originalDir); rmdir(tempDir,'s');
This example uses the following helper functions:
3GPP TS 36.101 "User Equipment (UE) radio transmission and reception"