Радиотехника и электроника, 2023, T. 68, № 1, стр. 95-102

Сравнение теоретической и реальной пропускной способности PD-NOMA

Я. В. Крюков a*, Д. А. Покаместов a, Е. В. Рогожников a

a Томский государственный университет систем управления и радиоэлектроники
634050 Томск, просп. Ленина, 40, Российская Федерация

* E-mail: kryukov.tusur@gmail.com

Поступила в редакцию 13.04.2022
После доработки 25.05.2022
Принята к публикации 06.06.2022

Полный текст (PDF)

Аннотация

Произведено сравнение теоретической и реальной пропускной способности в системах PD-NOMA. Теоретическая пропускная способность получена с использованием идеальных сигнальных конструкций по теореме Шеннона, а практическая – с учетом использования банка сигнально-кодовых конструкций из стандарта 3GPP LTE. Описан подход для оценки реальной пропускной способности PD-NOMA. Найдена величина погрешности между теоретической и реальной пропускными способностями PD-NOMA.

ВВЕДЕНИЕ

Сотовые системы связи непрерывно развиваются с момента их появления, а начиная с 2020 г. в мировую телекоммуникационную инфраструктуру внедряются сети пятого поколения 5G New Radio (NR). Прогресс беспроводной связи является обязательным условием для развития подавляющего большинства сфер жизнедеятельности людей. Информация стала самым ценным мировым ресурсом, определяющим все экономические взаимоотношения. В свою очередь телекоммуникационные системы – это способ передачи этого ресурса, и именно поэтому наблюдается экспоненциальный рост трафика в беспроводных сетях, благодаря появлению новых телекоммуникационных услуг, развитию межмашинных сетей, применению новых связных технологий, развитию новых направлений в связи.

Одним из ключевых факторов эффективности многоканальной системы связи является множественный доступ (МД). Он отвечает за организацию многоканальной связи, т.е. использование множеством абонентов общего ограниченного ресурса связи системы, осуществляя разделение их каналов в одной или нескольких плоскостях (доменах). В современных беспроводных сетях (например, 5G NR) применяется два подхода для разделения абонентских каналов: 1) ортогональное частотное разделение каналов OFDMA (Orthogonal Frequency-Division Multiple Access) и 2) пространственное разделение каналов SDMA (Space Division Multiple Access). В их основу заложен принцип ортогональности, т.е. максимально возможное исключение межканальной интерференции при формировании группового канала. Такой МД называют ортогональным (Orthogonal Multiple Access – OMA).

В поисковых работах 2012–2013 гг. по тематике 5G рассмотрено применение методов неортогонального множественного доступа (Non-Orthogonal Multiple Access – NOMA) и показана их потенциальная эффективность. В отличии от OMA, группа методов NOMA позволяет увеличить системную пропускную способность за счет контролируемой межканальной интерференции, которая устраняется в точке приема. Одним из наиболее перспективных подходов NOMA является неортогональное мультиплексирование каналов в мощностном домене (Power Domain-NOMA–PD-NOMA).

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

Несмотря на то, что технология PD-NOMA рассматривалась в качестве кандидата для систем связи нового поколения, она не вошла в техническую спецификацию 5G NR. Считается, что главная причина этого – высокая вычислительная сложность обработки PD-NOMA, поэтому ее реализация в сигнальных процессорах массовых мобильных устройств коммерчески невыгодна. Кроме того, большой интерес сегодня представляет технология Massive MIMO (Multiple Input Multiple Output), которая позволяет кратное увеличение пропускной способности системы за счет применения сложных многоэлементных антенных решеток на стороне базовых станций (БС) и не приводит к удорожанию абонентского оборудования. Несмотря на это, PD-NOMA не теряет своей актуальности и продолжает рассматривается в качестве кандидата для будущих систем связи и последующих модернизаций 5G NR.

В открытом доступе находится большое количество работ по тематике PD-NOMA. Впервые эта идея и ее потенциальная эффективность продемонстрирована в работе [1], но оказалась не востребована из-за недостаточного уровня технического прогресса. Интерес вновь возник спустя 40 лет и в работе) [2] была описана современная концепция применения PD-NOMA. Начиная с 2013 г. вышло большое количество работ по PD‑NOMA, однако, до сих пор остается ряд актуальных нерешенных проблем по различным направлениям: объединение абонентов и распределение мощности [3], организация случайного доступа к среде [4], декодирование [5], совместная работа с MIMO [6] и др. Тщательный обзор трендов, подходов и проблем PD-NOMA приведен в работе [7].

В работах [811] детально описана концепция PD-NOMA в контексте применения ее в 5G NR, а так же демонстрируется пропускная способность PD-NOMA и приведено ее равнение ее с OMA. Основной недостаток этих работ заключается в расчете пропускной способности по теореме Шеннона–Хартли, которая подразумевает рад допущений и не позволяет получить высокую точность результатов, близких к практическим. Например, к наиболее важным допущениям относится предположение о том, что для передачи используются идеальные сигнально-кодовые конструкции (СКК). Однако, точная оценка пропускной способности не может быть осуществлена без учета характеристик реальных СКК, а в свободном доступе такие исследования нами не найдены.

Мы считаем вышеуказанную проблему актуальной, поэтому целью данной работы является оценка практической пропускной способности PD-NOMA, без использования указанного допущения и без учета характеристик реальных СКК.

1. СИСТЕМНАЯ МОДЕЛЬ PD-NOMA

Для исследований PD-NOMA чаще всего используется системная модель нисходящего группового канала с одной приемной и одной передающей антенной. Предполагается, что БС расположена в центре соты и обслуживает двух и более пользователей UE (User Equipment). Задача БС заключается в организации нисходящего группового канала, мультиплексируя UE в мощностном домене. Теоретически, их количество не ограничено, но принято рассматривать наиболее реализуемый на практике сценарий, подразумевающий $K = 2$ пользователей. Увеличение K не дает выгодный прирост пропускной способности системы, но существенно увеличивает вычислительную сложность формирования и обработки сигналов, поэтому сценарии с $K \geqslant 3$ рассматриваются крайне редко.

Обозначим переменными x, y передаваемые информационные символы первого (UE1) и второго (UE2) пользователя. Формирование группового PD-NOMA сигнала $s$ для двух пользователей можно описать с помощью соотношения

(1)
$s = \sqrt {{{p}_{1}}} x + \sqrt {{{p}_{2}}} y,$

где ${{p}_{1}},{{p}_{2}}$ – это весовые мощностные коэффициенты, которые отвечают за разделение общего мощностного ресурса $P$ между UE1 и UE2, при условии $P = {{p}_{1}} + {{p}_{2}}$. Чем меньше отношение сигнал/шум (ОСШ) пользователя, тем больше его весовой мощностной коэффициент, и наоборот. Распределение ${{p}_{1}},{{p}_{2}}$ является одной из приоритетных задач, потому что величина мощностного коэффициента напрямую влияет на помехоустойчивость каждого канала. Поскольку происходит разделение общего доступного мощностного ресурса $P$, то увеличение ${{p}_{1}}$ приводит к уменьшению ${{p}_{2}}$ и, наоборот. При распределении необходимо найти компромисс, который одновременно удовлетворяет обоих пользователей и учитывает системную стратегию распределения радиоресурсов.

В системной модели PD-NOMA удобно считать $P = 1$, воспринимая ${{p}_{1}},{{p}_{2}}$ в качестве долей единицы. Отсортируем пользователей в порядке увеличения их ОСШ и запишем ограничения для мощностных коэффициентов

(2)
$\left\{ {\begin{array}{*{20}{l}} {0 < {{p}_{2}} < {{p}_{1}} < 1,} \\ {{{p}_{1}} + {{p}_{2}} = P.} \end{array}} \right.$

В современных системах для передачи информационных символов используются СКК, состоящие из многопозиционной модуляции и помехоустойчивого кодирования. Чаще всего применяются СКК на базе квадратурной амплитудной модуляции (Quadrature Amplitude Modulation – QAM) и Turbo (или Low-Density Parity-Check – LDPC) помехоустойчивых кодов.

В описываемой системной модели используется квадратная QAM с порядками 4, 16, 64. Обозначим символами ${{Q}_{1}}$ и ${{Q}_{2}}$ порядок QAM-модуляции в канале UE1 и UE2, соответственно. Возможные состояния информационных символов $x$, $y$ определены алфавитами ${{{\mathbf{A}}}_{1}}$, ${{{\mathbf{A}}}_{2}}$ используемых ${{Q}_{1}}$-QAM, ${{Q}_{2}}$-QAM. Созвездие группового сигнала $s$, сформированного с помощью (1), является суперпозицией взвешенных по ${{p}_{1}},{{p}_{2}}$ созвездий ${{Q}_{1}}$-QAM и ${{Q}_{2}}$-QAM. В качестве примера, на рис. 1 представлено созвездие PD-NOMA сигнала с конфигурацией ${{p}_{1}} = 0.9$, ${{Q}_{1}} = 16$ и ${{p}_{2}} = 0.1$, ${{Q}_{2}} = 16$.

Рис. 1.

Сигнальное созвездие группового PD-NOMA сигнала: а – созвездие UE1, б – созвездие UE2, в – созвездие PO‑NOMA.

Сигнал в точке приема ${{r}_{k}}$ пользователя $k$ может быть описан с помощью выражения

(3)
${{r}_{k}} = {{\alpha }_{k}}\left( {\sqrt {{{p}_{1}}} x + \sqrt {{{p}_{2}}} y} \right) + {{n}_{k}},$

где ${{\alpha }_{k}}$ – это канальный коэффициент и ${{n}_{k}}$ – аддитивная помеха в канале пользователя $k$.

Декодирование сигнала $r$ осуществляется с помощью метода последовательного подавления помех SIC (Serial Interference Cancellation) [12]. В приемнике UE1 происходит декодирование ${{r}_{1}}$ и оценка его информационного символа $\hat {x}$ может быть получена с помощью критерия минимального евклидового расстояния, если воспринимать сумму $y$ и ${{n}_{1}}$ в качестве аддитивной помехи

(4)
$\hat {x} = {\text{arg}}\mathop {\min }\limits_{x \in {{A}_{1}}} {{\left| {{{r}_{1}} - \sqrt {{{p}_{1}}} x} \right|}^{2}}.$

В приемнике UE2 декодирование сигнала реализуется в два последовательных этапа. Сначала осуществляется оценка $\hat {x}$ из ${{r}_{2}}$ точно так же, как в (4). Далее, сигнал первого пользователя $x$ регенерируется и вычитается из ${{r}_{2}}$, осуществляя устранение его интерференционной помехи

(5)
$r{\kern 1pt} ' = {{r}_{2}} - \sqrt {{{p}_{1}}} \hat {x}.$

После этого становится возможным получить оценку собственного информационного символа $\hat {y}$ из $r{\kern 1pt} '$, воспринимая в качестве помехи только ${{n}_{2}}$

(6)
$\hat {y} = {\text{arg}}\mathop {\min }\limits_{y \in {{A}_{2}}} {{\left| {r{\kern 1pt} '\,\, - \sqrt {{{p}_{2}}} y} \right|}^{2}}.$

2. РАСЧЕТ ПРОПУСКНОЙ СПОСОБНОСТИ PD-NOMA ПО ШЕННОНУ

Для представленной модели PD-NOMA пропускная способность может быть рассчитана с помощью формулы Шеннона. Обозначим переменными $C_{1}^{T}$ и $C_{2}^{T}$ теоретическую пропускную способность каналов UE1 и UE2 в PD-NOMA, которая может быть получена по формулам

(7)
$\begin{gathered} C_{1}^{T} = B{{\log }_{2}}\left( {1 + \frac{{{{{\left| {{{\alpha }_{1}}} \right|}}^{2}}{{p}_{1}}}}{{{{{\left| {{{\alpha }_{1}}} \right|}}^{2}}{{p}_{2}} + B{{N}_{0}}}}} \right);~ \\ ~C_{2}^{T} = B{{\log }_{2}}\left( {1 + \frac{{{{{\left| {{{\alpha }_{2}}} \right|}}^{2}}{{p}_{2}}}}{{B{{N}_{0}}}}} \right)~, \\ \end{gathered} $

где ${{N}_{0}}$ – спектральная плотности мощности шума в общей полосе пропускания $B$. Предполагается, что значение ${{N}_{0}}$ одинаково для обоих абонентов.

При расчете $C_{1}^{T}$ межканальная помеха мощностью ${{p}_{2}}$ от UE2 воспринимается в качестве нормальной помехи. С другой стороны, при расчете $C_{2}^{T}$ предполагается полное устранение интерференционной помехи мощностью ${{p}_{1}}$, поэтому остается декодировать собственный сигнал мощностью ${{p}_{2}}$ на фоне шумовой составляющей $B{{N}_{0}}$.

Обозначим ОСШ относительно общей мощности P в канале пользователя $k$ переменной ${{W}_{k}}$, тогда

${{W}_{k}} = \frac{{P{{{\left| {{{\alpha }_{k}}} \right|}}^{2}}}}{{B{{N}_{0}}}}.$

Преобразуем предыдущее выражение, учитывая принятую нормировку $P = 1$

(8)
${{W}_{k}} = \frac{{{{{\left| {{{\alpha }_{k}}} \right|}}^{2}}}}{{B{{N}_{0}}}},$

и перепишем (7) с учетом (8)

(9)
$\begin{gathered} C_{1}^{T} = B{{\log }_{2}}\left( {1 + \frac{{{{p}_{1}}{{W}_{1}}}}{{{{p}_{2}}{{W}_{1}} + 1}}} \right);~ \\ C_{2}^{T} = B{{\log }_{2}}\left( {1 + {{p}_{2}}{{W}_{2}}} \right)~. \\ \end{gathered} $

Общая пропускная способность системы в полосе B может быть рассчитана с помощью соотношения

(10)
$C_{\Sigma }^{T} = C_{1}^{T} + C_{2}^{T},$

потому что передача ведется одновременно двумя абонентами в единой полосе.

Основная проблема расчета по (9) заключается в использовании упрощающих допущений:

1) при расчете $C_{1}^{T}$ интерференционная помеха от UE2 мощностью ${{p}_{2}}$ воспринимается как эквивалент нормальной помехи с идентичной мощностью;

2) при расчете $C_{2}^{T}$ не учитывается вероятность ошибки при устранении интерференционной помехи при расчете (4), (5). Предполагается, что интерференционная помеха от UE1 мощностью ${{p}_{1}}$ полностью устраняется;

3) предполагается использование идеальных СКК, демонстрирующих предельную спектральную эффективность по Шеннону;

4) не накладывается ограничение на распределение ${{p}_{1}}$, ${{p}_{2}}$, которое неизбежно возникает при использовании реальных СКК.

3. РАСЧЕТ РЕАЛЬНОЙ ПРОПУСКНОЙ СПОСОБНОСТИ PD-NOMA

Рассмотрим собственный подход, позволяющий рассчитать реальную пропускную способность системы PD-NOMA, не прибегая к упомянутым выше допущениям.

3.1. Выбор СКК в PD-NOMA

Реальная пропускная способность ${{C}^{R}}$ в пользовательском канале напрямую определяется спектральной эффективностью $\gamma $ используемой СКК, которая выбрана для передачи данных этого пользователя. В современных системах используется принцип адаптивного выбора СКК, позволяющий приспособиться к условиям канала распространения и выбирать СКК с наибольшей спектральной эффективностью, сохраняя при этом требуемый уровень вероятности битовых ошибок при приеме. Как правило, для каждой конкретной системы связи набор доступных СКК формируется исходя из предполагаемых эксплуатационных сценариев и рабочего диапазона ОСШ. Все доступные СКК описаны в стандарте конкретной системы связи, а каждая из них определена типом модуляции, скоростью и методом помехоустойчивого кодирования.

В качестве примера можно рассмотреть набор СКК из стандарта LTE, представленный в [13] (табл. 1 ). Набор включает в себя ${{M}_{o}} = 15$ разных СКК, пронумерованных индексом $m \in 1, \ldots ,{{M}_{o}}$. Каждая конструкция определена следующими параметрами: $Q$ – порядок QAM-модуляции, $R$ – скорость помехоустойчивого Turbo кодирования, $Z$ – пороговый уровень ОСШ, $\gamma $ спектральная эффективность. Пороговый уровень $Z$ показывает минимальное значение ОСШ в канале, при котором для передачи может быть выбрана данная СКК.

Используя данные указанной таблицы адаптивный выбор индекса СКК в системе OMA для пользователя $k$ может быть описан с помощью простого выражения

(11)
${{m}_{k}} = {\text{max}}\left\{ {i,{{Z}_{i}} \leqslant {{W}_{k}}} \right\},\,\,\,\,i = 1,2, \ldots ,{{M}_{o}},$

где ${{W}_{k}}$ – текущее значение ОСШ в канале пользователя $k$. Поскольку межканальная интерференция отсутствует, то СКК для каждого пользователя выбирается независимо от остальных.

Как говорилось ранее, в отличии от OMA, в PD-NOMA осуществляется передача информационных символов сразу двух пользователей в едином частотно-временном ресурсном сегменте. Для пользователей могут быть выбраны разные СКК с различными мощностными коэффициентами, увеличивая или уменьшая взаимную интерференционную помеху. Как следствие, увеличение $C_{1}^{R}$ приводит к уменьшению $C_{2}^{R}$, и, наоборот. Поэтому должен быть реализован совместный адаптивный выбор пары СКК с индексами ${{m}_{1}},{{m}_{2}}$ и расчет мощностных коэффициентов сразу для двух абонентов UE1 и UE2, одновременно учитывая характеристики каждого из них и достигая компромисса.

Распределение ресурсов связи является важной задачей в любой многоканальной системе. Как правило, в OMA осуществляется распределение ресурса только в частотно-временном домене, однако, в PD-NOMA добавляется дополнительный мощностной домен. В свою очередь, распределение ресурса в мощностном домене напрямую влияет на совместный выбор СКК. Рассмотрим эту проблему более подробно.

Существует несколько подходов к расчету мощностных коэффициентов. В работе [14] предложено использовать адаптивный к каналу фазовый поворот созвездия QAM и алгоритм расчета мощностных коэффициентов с помощью треугольника Паскаля. Недостаток работы заключается в том, что рассматривается использование только QPSK и 64-QAM, а алгоритм расчета коэффициентов, по оценке самих авторов, имеет высокую вычислительную сложность.

В работе [15] предложен подход, согласно которому в первую очередь определяется СКК и минимальный мощностной коэффициент для UE1, а далее, по остаточному принципу для UE2. Однако, в этой работе интерференционная помеха представляется в качестве аддитивного белого гауссовcкого шума (АБГШ), что является допущением и увеличивает погрешность результатов.

В работе [16] предложена карта парного выбора СКК для двух пользователей в зависимости от их ОСШ, а также описано преимущество и результативность этого подхода. Недостатком работы является малое количество рассматриваемых СКК (8 типов).

Мы применяем близкий к [16] подход, используя СКК из [13] (табл. 1 ) для формирования таблицы парных СКК для PD-NOMA и расширяя количество доступных СКК по сравнению с [16]. Кроме того, мы опираемся на СКК из реального стандарта связи и формируем парные СКК с помощью собственного подхода, основываясь на критерии максимизации пропускной способности системы. Предполагается, что для декодирования PD-NOMA сигнала применяется SIC на уровне символов модуляции (без исправления ошибок с помощью помехоустойчивого кодирования при устранении сигнала UE1 в приемнике UE2).

С помощью аналитических расчетов и математического моделирования разработан банк из ${{M}_{n}} = 185$ парных СКК для системы PD-NOMA, которые приведены в [13] (табл. 2 ). С помощью предложенного банка можно осуществить адаптивный к ОСШ обоих пользователей совместный выбор двух СКК, которые могут быть использованы при мультиплексировании PD-NOMA и не приведут к превышению порогового уровня вероятности битовых ошибок.

Используя данные упомянутой таблицы из [13], адаптивный выбор совместной пары индексов СКК ${{m}_{1}},{{m}_{2}}$ для PD-NOMA можно записать виде

(12)
$\left\{ {{{m}_{1}},{{m}_{2}}} \right\} = {\text{max}}\left( {\left. i \right|\left\{ {\begin{array}{*{20}{c}} {{{T}_{{i,1}}} \leqslant {{W}_{1}},} \\ {{{T}_{{i,2}}} \leqslant {{W}_{2}},} \end{array}~~} \right.i = 1,2, \ldots ,{{M}_{n}}} \right),$

где ${{m}_{1}}~$ – номер СКК (из табл. 1 работы [13]) для UE1 с его ОСШ ${{W}_{1}}$, данные которого располагаются на неортогональном слое с долей мощности ${{p}_{1}}$; ${{m}_{2}}$ – номер СКК для UE2 с его ОСШ ${{W}_{2}}$, данные которого располагаются на неортогональном слое с долей мощности ${{p}_{2}}$; ${{T}_{{i,1}}}$ и ${{T}_{{i,2}}}$ – соответствующие пороговые значения ОСШ для выбора i-й пары. Согласно (12) осуществляется выбор пары ССК, которая обеспечивает максимальную суммарную спектральную эффективность единицы частотно-временного ресурса при заданных ${{W}_{1}},{{W}_{2}}$.

Реальная пропускная способность пользователей в полосе B определяется спектральной эффективностью ${{\gamma }_{{{{m}_{1}}}}}$,$~{{\gamma }_{{{{m}_{2}}}}}$ выбранных для них СКК с помощью (12)

(13)
$C_{1}^{R} = B{{\gamma }_{{{{m}_{1}}}}};\,\,\,\,C_{2}^{R} = B{{\gamma }_{{{{m}_{2}}}}},$

а общая пропускная способность системы $C_{\Sigma }^{R}~$ рассчитывается по (10).

Для примера, рассмотрим выбор пары СКК в PD-NOMA для сценария ${{W}_{1}} = 10$ и ${{W}_{2}} = 30$ дБ. Согласно (12) будет выбрана пара $i = 143$, соответствующая ${{m}_{1}} = 6$ и ${{m}_{2}} = 15$, а суммарная спектральная эффективность составит γ6 + γ15 = $ = 6.7305$ бит/с/Гц. При этом гарантировано, что можно рассчитать такие коэффициенты ${{p}_{1}},~{{p}_{2}}$, учитывая (2), которые обеспечат вероятность ошибки на блок (BLER, Block Error Rate) менее 10% при обработке группового сигнала методом SIC.

3.2. Ограничение на соотношение мощностных коэффициентов

Использование QAM сигналов накладывает ограничение на соотношение мощностных коэффициентов, нарушение которого приводит к выходу созвездия UE2 за границы области детектирования символов на созвездии UE1 и существенному росту вероятности ошибок декодирования $\hat {x}$ в (4). В системной модели PD-NOMA по Шеннону это ограничение не учитывается, потому что использование QAM не рассматривается.

Получим ограничение на соотношение мощностных коэффициентов. Минимальное евклидово расстояние между соседними точками на взвешенном по мощности сигнальном созвездии ${{Q}_{k}}$-QAM (рис. 1) для пользователя $k$ можно рассчитать с помощью

(14)
${{D}_{k}} = 2\sqrt {\frac{{{{p}_{k}}}}{{\frac{2}{3}\left( {{{Q}_{k}} - 1} \right)}}} .$

Расстояние между крайними точками созвездия UE2 на квадратурной (или синфазной) оси рассчитывается с помощью $\left( {\sqrt {{{Q}_{2}}} - 1} \right){{D}_{2}}$, а ширина области детектирования символов на созвездии UE1 равна ${{D}_{1}}$. Тогда, чтобы избежать выхода созвездия UE2 за границы области детектирования символов на созвездии UE1, геометрическое ограничение можно записать виде

(15)
$\left( {\sqrt {{{Q}_{2}}} - 1} \right){{D}_{2}} < {{D}_{1}}.$

Обозначим переменными $p_{1}^{{{\text{min}}}}$ и $p_{1}^{{{\text{max}}}}$ минимальный и максимальный граничный мощностной коэффициент UE1 при известных ${{Q}_{1}}$, ${{Q}_{2}}$. Величина $p_{1}^{{{\text{min}}}}$ может быть получена из (15) с помощью (14), а $p_{1}^{{{\text{max}}}} = 1$ из (2). Запишем ограничения для $p_{1}^{{{\text{min}}}}$ и $p_{1}^{{{\text{max}}}}$

(16)
$\left\{ {\begin{array}{*{20}{l}} {p_{1}^{{{\text{min}}}} = \frac{{\left( {\sqrt {{{Q}_{2}}} - 1} \right)\left( {{{Q}_{1}} - 1} \right)}}{{{{Q}_{1}}\left( {\sqrt {{{Q}_{2}}} - 1} \right) + 2}},} \\ {p_{1}^{{{\text{max}}}} = 1.} \end{array}} \right.~$

Учитывая (2) и (16), запишем граничные коэффициенты для $p_{2}^{{{\text{min}}}}$ и $p_{2}^{{{\text{max}}}}$

(17)
$\left\{ {\begin{array}{*{20}{l}} {p_{2}^{{{\text{min}}}} = 0,} \\ {p_{2}^{{{\text{max}}}} = 1 - p_{1}^{{{\text{min}}}}.} \end{array}} \right.~$

Используя (16), (17), общее ограничение на распределение ${{p}_{1}},~{{p}_{2}}~$ будет

(18)
$\left\{ {\begin{array}{*{20}{l}} {\frac{{\left( {\sqrt {{{Q}_{2}}} - 1} \right)\left( {{{Q}_{1}} - 1} \right)}}{{{{Q}_{1}}\left( {\sqrt {{{Q}_{2}}} - 1} \right) + 2}}~ < {{p}_{1}} \leqslant 1,} \\ {0 \leqslant {{p}_{2}} < 1 - \frac{{\left( {\sqrt {{{Q}_{2}}} - 1} \right)\left( {{{Q}_{1}} - 1} \right)}}{{{{Q}_{1}}\left( {\sqrt {{{Q}_{2}}} - 1} \right) + 2}}.} \end{array}} \right.$

В случае ${{p}_{1}} = 1$, ${{p}_{2}} = 0~$ система переходит в режим OMA, т.к. весь доступный мощностной ресурс отдается для передачи данных UE1.

4. СРАВНЕНИЕ ТЕОРЕТИЧЕСКОЙ И РЕАЛЬНОЙ ПРОПУСКНОЙ СПОСОБНОСТИ PD-NOMA

Приведем сравнение теоретической $C_{\Sigma }^{T}$ и реальной $C_{\Sigma }^{R}$ пропускной способности системы PD-NOMA в зависимости от ОСШ пользователей ${{W}_{1}},{{W}_{2}}$ в диапазоне от –10 до 40 дБ с шагом ${{\Delta }}W = 0.5$ дБ. Выбор диапазона обусловлен тем, что он встречается на практике при эксплуатации реальных систем связи. Полоса пропускания системы нормирована $B = 1$.

Максимальное значение $C_{k}^{R}$ ограничено спектральной эффективностью самой высокоскоростной СКК с индексом $m = 15$, т.е. ${\text{max}}\left( {C_{k}^{R}} \right) = {{\gamma }_{{15}}}$. Для справедливого сравнения максимальное значение $C_{k}^{T}$ искусственно ограничено, ${\text{max}}\left( {C_{k}^{T}} \right) = {{\gamma }_{{15}}}$.

Подбор весовых мощностных коэффициентов реализовали с помощью грубого перебора значения ${{p}_{1}}$ с шагом ${{\Delta }}{{p}_{1}} = 0.005$, используя (18). Для этого применяется разработанная ранее имитационная математическая модель системы PD-NOMA с каналом АБГШ, которая позволяет рассчитать BLER на стороне каждого пользователя при заданной конфигурации (${{Q}_{1}},{{Q}_{2}}$ и ${{p}_{1}},{{p}_{2}}$) мультиплексирования в мощностном домене. Эта модель более подробно описана в [13]. Таким образом, для каждой совместной пары СКК из табл. 2 указанной работы, методом подбора найдены такие значения ${{p}_{1}},{{p}_{2}}$, при которых BLER обоих абонентов не превышает 10% при заданных ${{W}_{1}},{{W}_{2}}$.

На рис. 2 представлена трехмерная зависимость реальной пропускной способности системы $C_{\Sigma }^{R}\left( {{{W}_{1}},{{W}_{2}}} \right)$ от ОСШ в канале UE1 и UE2. Опишем полученный результат:

Рис. 2.

Зависимость реальной (а) и теоретической (б) пропускной способности системы PD-NOMA от ОСШ обоих пользователей.

1) при высоком значении ОСШ в каналах обоих пользователей (${{W}_{1}},{{W}_{2}} > 30$ дБ) БС выбирает высокоскоростные СКК для обоих пользователей. В наилучшем сценарии можно получить двойной выигрыш в пропускной способности по сравнению с OMA, потому что одновременно двум пользователям будет доступна для передачи СКК с максимальной спектральной эффективностью;

2) высокое ОСШ в канале одного из пользователей (${{W}_{{1|2}}} > 30$ дБ). В большинстве работ PD-NOMA описывается подобный сценарий, рассматривающий мультиплексирование каналов пользователей с низким (на окраине соты) и высоким (в центре соты) ОСШ. В этом случае БС назначает одному пользователю высокоскоростную СКК, а второй пользователь ведет передачу с помощью помехоустойчивой, но низкоскоростной СКК;

3) при низком ОСШ обоих пользователей (${{W}_{1}},{{W}_{2}} < 15$ дБ) БС назначает им низкоскоростные СКК, а PD-NOMA предоставляет меньшую спектральную эффективность, чем OMA. Это связано с тем, что суммарное влияние введённой межканальной помехи и АБГШ не позволяет производить декодирование PD-NOMA методом SIC с требуемым уровнем ошибок.

На рис. 3 представлена зависимость теоретической пропускной способности $C_{\Sigma }^{T}\left( {{{W}_{1}},{{W}_{2}}} \right)$, полученная с помощью формулы (9) с использованием идентичных коэффициентов ${{p}_{1}},{{p}_{2}}$. Видно существенное различие между $C_{\Sigma }^{T}$ и $C_{\Sigma }^{R}$, что демонстрирует относительно невысокую точность оценки пропускной способности PD-NOMA по (9).

Рис. 3.

Зависимости $C_{\Sigma }^{T}\left( {{{W}_{2}}} \right)$ (штриховая кривая) и $C_{\Sigma }^{R}\left( {{{W}_{2}}} \right)$ (сплошная кривая) при ${{W}_{1}} = 0$ дБ для PD-NOMA (а) и ОМА (б).

Чтобы наглядно оценить величину погрешности между $C_{\Sigma }^{T}$ и $C_{\Sigma }^{R}$ рассмотрим их поведение на двумерной плоскости при фиксированном значении ОСШ одного из пользователей. На рис. 3а приведены $C_{\Sigma }^{T}\left( {{{W}_{2}}} \right)$ и $C_{\Sigma }^{R}\left( {{{W}_{2}}} \right)$ при фиксированном значении ${{W}_{1}} = 0$ дБ.

Обозначим разность в ${{W}_{2}}$ при $C_{\Sigma }^{T} = ~C_{\Sigma }^{R}$ переменной $\Delta {{W}_{n}}$. Величина $\Delta {{W}_{n}}$ показывает разность в ОСШ, которая должна быть скомпенсирована для достижения $C_{\Sigma }^{T} = ~C_{\Sigma }^{R}$. Иными словами, это различие в помехоустойчивости между теоретической и реальной системами. Из данных, приведенных на рисунке видно, что значение $\Delta {{W}_{n}}$ в измеряемом диапазоне значений составляет, в среднем, $\Delta {{W}_{n}} \approx 8{\kern 1pt} - {\kern 1pt} 10$ дБ. При ${{W}_{2}} > 33$ дБ достигается ${\text{max}}\left( {{{C}_{2}}} \right) = ~{{\gamma }_{{15}}}$, а значение ${{C}_{\Sigma }}$ ограничено низким ${{C}_{1}}$ из-за малого ${{W}_{1}}$ (область обведена на рисунке).

Далее, получим оценку погрешности между $C_{\Sigma }^{T}$ и $C_{\Sigma }^{R}~$ для систем OMA в идентичных условиях. На рис. 3б представлены зависимости $C_{\Sigma }^{T}\left( {{{W}_{2}}} \right)$ и $CС_{\Sigma }^{R}\left( {{{W}_{2}}} \right)$ при ${{W}_{1}} = 0$ дБ. В рамках мультиплексирования OMA оба пользователя используют для передачи всю доступную мощность P, однако поровну разделяют частотный ресурс B, поэтому ${{B}_{1}} = {{B}_{2}} = 0.5.$ Выбор СКК происходит согласно (11), а расчет $C_{\Sigma }^{R}$ с помощью (13). Разность в ${{W}_{2}}$ между $C_{\Sigma }^{T}\left( {{{W}_{2}}} \right)$ и $C_{\Sigma }^{R}\left( {{{W}_{2}}} \right)$ для OMA обозначена переменной ${{\Delta }}{{W}_{o}}$. Из данных, приведенных на рисунке видно, что, в среднем, ${{\Delta }}{{W}_{o}} \approx 3{\kern 1pt} - {\kern 1pt} 5$ дБ. Для демонстрируемого сценария величина ${{\Delta }}{{W}_{o}}$ меньше, чем $\Delta {{W}_{n}}$ приблизительно в два раза.

Сравнивая зависимости PD-NOMA (рис. 3а) с OMA (рис. 3б) можно сделать вывод о том, что теоретический расчет пропускной способности PD‑NOMA обладает меньшей точностью по сравнению с OMA, из-за большей погрешности ${{\Delta }}{{W}_{n}} > {{\Delta }}{{W}_{o}}$.

Ранее говорилось, что в большинстве работ сравнение пропускной способности систем PD‑NOMA и OMA осуществляется исключительно с помощью теоремы Шеннона. Однако нами показано, что теоретическое сравнение является недостоверным, а расчетная пропускная способность будет существенно отличаться.

Недостоверность сравнения наблюдается, например, при значении ${{W}_{1}} = 10$ дБ. При теоретическом расчете $C_{\Sigma }^{T} = 4$ для PD-NOMA и $C_{\Sigma }^{T} = 3$ для OMA, соответственно PD-NOMA обладает 25% теоретическим выигрышем. Однако, реальная пропускная способность системы $C_{\Sigma }^{R} = 0.95$ для PD-NOMA и $C_{\Sigma }^{R} = 1.5$ для OMA, что соответствует 35% проигрышу PD-NOMA.

Отметим, что при иной конфигурации мультиплексирования могут быть получены результаты сравнения, отличные от наших. Однако, при любой конфигурации погрешность PD-NOMA окажется больше $\Delta {{W}_{n}} > \Delta {{W}_{o}}$ из-за допущений, описанных в разд. 3.

ЗАКЛЮЧЕНИЕ

Основной вывод работы заключается в том, что полученная другими исследователями расчетная пропускная способность PD-NOMA является недостижимой в реальных условиях и требует более точной оценки. Для этого предложена набор парных СКК, позволяющий выполнить адаптивный выбор конфигурации PD-NOMA мультиплексирования двух пользователей в канале с АБГШ.

Приведено сравнение теоретической и реальной пропускной способности системы PD-NOMA. В первом случае для расчетов используется теорема Шеннона, которая применяется в большинстве научных работ. Реальная пропускная способность получена с учетом СКК из стандарта 3GPP LTE.

Показано, что реальная пропускная способность PD-NOMA существенно отличается от теоретической из-за использования упрощающих допущений, а погрешность может является основной причиной некорректных результатов в работах, посвященных сравнению PD-NOMA с OMA.

Продемонстрировано, что результат сравнения теоретической пропускной способности систем PD-NOMA и OMA может не соответствовать результату сравнения реальной пропускной способности.

Авторы заявляют об отсутствии конфликта интересов.

Список литературы

  1. Cover T. // IEEE Trans. 1972. V. IF-18. № 1. P. 2.

  2. Benjebbour A., Saito Y., Kishiyama Y. et al. // Proc. 2013 IEEE Int. Symp. on Intelligent Signal Processing and Communication Systems. Naha. 12–15 Nov. N.Y.: IEEE, 2013. P. 770.

  3. Liu X., Wang J., Zhao N. et al. // IEEE Wireless Commun. Lett. 2019. V. 8. № 3. P. 965.

  4. Zhang J., Tao X., Wu H. et al. // IEEE Internet of Things J. 2020. V. 7. № 7. P. 6369.

  5. Ding Z., Schober R., Poor H. // IEEE Commun. Lett. 2020. V. 24. № 11. P. 2373.

  6. Wang J., Li Y., Ji C. et al. // IEEE Trans. 2020. V. COM-68. № 4. P. 2293.

  7. Dai L., Wang B., Ding Z. et al. // IEEE Commun. Surveys & Tutorials. 2018. V. 20. № 3. P. 2294.

  8. Saito Y., Benjebbour A., Kishiyama Y., Nakamura T. // Proc. 2013 IEEE 24th Annual Int. Symp. on Personal, Indoor and Mobile Radio Commun. (PIMRC). London. 8–11 Sept. N.Y.: IEEE, 2013. P. 611.

  9. Ding Z., Lei X., Karagiannidis K. et al. // IEEE J. Selected Areas Commun. 2017. V. 35. № 10. P. 2181.

  10. Benjebbour A., Saito K., Li A. et al. // 2015 Int. Conf. on Wireless Networks and Mobile Commun.(WINCOM). Marrakesh. 20–23 Oct. N.Y.: IEEE, 2015. Article № 7381343.

  11. Liu C.H., Liang D.C. // IEEE Trans. 2018. V. WC-17. № 5. P. 3524.

  12. Higuchi K., Benjebbour A. // IEICE Trans. 2015. V. COM-98. № 3. P 403.

  13. Kryukov Y., Pokamestov D., Abenov R. et al. // J. Phys: Conf. Series. 2021. V. 2134. № 1. P. 12023.

  14. Abdel Moniem M., Gasser S., El-Mahallawy M. et al. // Applied Sciences. 2019. V. 9. № 15. P. 3022.

  15. Kim S., Kim H., Hong D. // Proc. 2018 29th Annual Int. Symp. on Personal, Indoor and Mobile Radio Commun. (RIMRC). Bologna. 9–12 Sept. N.Y.: IEEE, 2018. Article № 8580995.

  16. Hsieh H.-Y., Yang M.-J., Wang C.-H. // Proc. 2016 IEEE 27th Annual Int. Symp. on Personal, Indoor, and Mobile Radio Commun. (RIMRC). Valencia. 4–8 Sept. N.Y.: IEEE, 2016. Article № 7794796.

Дополнительные материалы отсутствуют.