Моделирование получателя РФ для приема 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) и прием"

Для просмотра документации необходимо авторизоваться на сайте