Моделирование получателя РФ для приема LTE

Этот пример демонстрирует, как смоделировать и протестировать LTE получатель РФ с помощью LTE Toolbox™ и RF Blockset™.

Пример требует LTE Toolbox.

Образцовое описание

Фигура ниже показов системная модель для этого примера. Форма волны LTE сгенерирована, отфильтрована, передана через канал распространения и подана в модель получателя RF Blockset. Модель может быть собрана с помощью коммерчески доступных частей. Измерения EVM выполняются на получателе РФ вывод.

Этот пример реализован с помощью MATLAB® и Simulink®, которые взаимодействуют во времени выполнения. Функциональный раздел:

Испытательный стенд измерения реализован со скриптом MATLAB с помощью модели Simulink в качестве устройства под тестом (DUT). Кадры LTE передаются потоком между испытательным стендом и DUT.

Сгенерируйте форму волны LTE

В этом разделе мы генерируем форму волны LTE с помощью LTE Toolbox. Мы используем ссылочный канал измерения (RMC) R.6, как задано в TS 36.101 [1]. Этот RMC задает 25 ресурсов элементы (REs) пропускная способность, эквивалентная 5 МГц. Используются 64 модуляции QAM. Все REs выделяются. Кроме того, шум OCNG включен в неиспользованном REs.

Только один кадр сгенерирован. Этот кадр будет затем повторен неоднократно, чтобы выполнить измерения EVM.

% 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);

Инициализируйте компоненты симуляции

Этот раздел инициализирует некоторые компоненты симуляции:

  • Полоса, ограничивающая фильтр: разработайте коэффициенты фильтра, которые будут использоваться моделью Simulink. Фильтр имеет порядок 32 с частотой полосы пропускания, равной 2,25 МГц и частотой полосы задерживания, равной 2,7 МГц.

  • NumberOfFrames : это - число раз, сгенерированный кадр повторяется

  • Предварительно выделите итоговые векторы

% 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

% 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

Загрузите испытательный стенд RF Blockset

Этот раздел загружает модель Simulink, показанную ниже. Модель выполняет следующие функции:

  • Чтение формы волны LTE и периода выборки из рабочей области

  • Фильтрация Bandlimiting

  • Сохранение результатов к рабочей области

Модель включает следующие компоненты:

  • Модель Channel: включает потерю пути к свободному пространству

  • Получатель РФ: включает прямой демодулятор преобразования

  • ADC и DC смещают отмену

% Specify and open Simulink model
model = 'simrfV2_lte_receiver';
open_system(model);

Модель получателя РФ

Модель получателя РФ включает элементы, показанные ниже

Эта модель была первоначально сгенерирована с помощью приложения RF Budget Analyzer, и позже изменена, чтобы включать дополнительные нарушения. Начальный Бюджетный Анализ РФ Получателя РФ

Моделируйте кадры

Этот раздел моделирует конкретное количество кадров. Это сделано на двух этапах:

  • Моделируйте первый кадр

  • Моделируйте остальную часть кадров в цикле

Причина разделения обработки на этих двух этапах состоит в том, чтобы упростить код. Во время обработки первого кадра мы должны учесть задержку полосы, ограничивающей фильтр. Дело обстоит не так для последующих кадров, поскольку состояние фильтра сохраняется между кадрами. Поэтому длина первого кадра должна быть увеличена немного, чтобы учесть задержку, введенную фильтром.

Моделируйте первый кадр LTE

Как упомянуто для первого моделируемого кадра мы должны увеличить длину сигнала, питаемого модель Simulink, чтобы компенсировать задержку, введенную фильтром. Затем мы запускаем симуляцию модели Simulink, не загружая начального состояния. Модель собирается моделировать в режиме Accelerator, чтобы уменьшить время выполнения. После обработки первого кадра с моделью Simulink ее состояние (xFinal) хранится и присвоенный xInitial в следующий раз, когда модель запущена.

Вывод модели Simulink хранится в переменной rx, которая доступна в рабочей области. Любые задержки, введенные этому сигналу, удалены после выполнения синхронизации. EVM измеряется на получившейся форме волны.

% 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');
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

Моделируйте последовательные кадры LTE

Теперь остальная часть кадров может быть моделирована. Во-первых, образцовое состояние установлено с помощью значения, сохраненного в xInitial при выводе предыдущей итерации.

% 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 = simrfV2_lte_receiver_evm_cal(rmc,cec,squeeze(rx),EVMalg);
    evmpeak(n-1) = evmmeas.Peak;
    evmrms(n-1) = evmmeas.RMS;
end
Low edge EVM, subframe 0: 3.165%
High edge EVM, subframe 0: 3.173%
Low edge EVM, subframe 1: 3.024%
High edge EVM, subframe 1: 3.036%
Low edge EVM, subframe 2: 3.004%
High edge EVM, subframe 2: 2.992%
Low edge EVM, subframe 3: 2.964%
High edge EVM, subframe 3: 2.927%
Low edge EVM, subframe 4: 3.003%
High edge EVM, subframe 4: 3.012%
Low edge EVM, subframe 6: 3.001%
High edge EVM, subframe 6: 2.983%
Low edge EVM, subframe 7: 3.005%
High edge EVM, subframe 7: 3.002%
Low edge EVM, subframe 8: 3.019%
High edge EVM, subframe 8: 3.002%
Low edge EVM, subframe 9: 3.050%
High edge EVM, subframe 9: 3.041%
Averaged low edge EVM, frame 0: 3.024%
Averaged high edge EVM, frame 0: 3.017%
Averaged EVM frame 0: 3.024%
Averaged overall EVM: 3.024%
Low edge EVM, subframe 0: 3.165%
High edge EVM, subframe 0: 3.173%
Low edge EVM, subframe 1: 3.024%
High edge EVM, subframe 1: 3.036%
Low edge EVM, subframe 2: 3.004%
High edge EVM, subframe 2: 2.992%
Low edge EVM, subframe 3: 2.964%
High edge EVM, subframe 3: 2.927%
Low edge EVM, subframe 4: 3.003%
High edge EVM, subframe 4: 3.012%
Low edge EVM, subframe 6: 3.001%
High edge EVM, subframe 6: 2.983%
Low edge EVM, subframe 7: 3.005%
High edge EVM, subframe 7: 3.003%
Low edge EVM, subframe 8: 3.019%
High edge EVM, subframe 8: 3.002%
Low edge EVM, subframe 9: 3.050%
High edge EVM, subframe 9: 3.041%
Averaged low edge EVM, frame 0: 3.024%
Averaged high edge EVM, frame 0: 3.017%
Averaged EVM frame 0: 3.024%
Averaged overall EVM: 3.024%

Визуализируйте измеренный EVM

Этот раздел строит измеренный пик и RMS EVM для каждого моделируемого кадра.

hf(1) = figure;
plot((2:N), 100*evmpeak,'o-')
title('EVM peak %');
xlabel('Number of frames');
hf(2) = figure;
plot((2:N), 100*evmrms,'o-');
title('EVM RMS %');
xlabel('Number of frames');

Чистка

Закройте модель Simulink и удалите сгенерированные файлы.

bdclose(model);
close(hf)
clear 'hf'

Приложение

Этот пример использует следующую функцию помощника:

Выбранная библиография

  1. 3GPP TS 36.101 "Передача радио оборудования пользователя (UE) и прием"