Noise-immunity of optimal symbol-by-symbol reception of signals with corrective coding in Galois fields for fading channels
- Authors: Nazarov L.E.1
-
Affiliations:
- Fryazino Branch Kotelnikov Institute of Radioengineering and Electronics RAS
- Issue: Vol 69, No 9 (2024)
- Pages: 887-893
- Section: THEORY AND METHODS OF SIGNAL PROCESSING
- URL: https://journals.rcsi.science/0033-8494/article/view/281994
- DOI: https://doi.org/10.31857/S0033849424090083
- EDN: https://elibrary.ru/HRGMUA
- ID: 281994
Cite item
Full Text
Abstract
A description of the algorithm for optimal symbol-by-symbol reception of signal structures based on correction coding in non-binary Galois fields is given. The results of modeling this algorithm are given in order to study its noise-immunity for transionospheric channels with fading due to scattering on ionospheric irregularities for a number of digital signals with multi-level phase shift keying in combination with a correction code with a parity check in Galois fields. It is shown that for these channels the use of a symbol-by-symbol reception algorithm provides an energy gain of up to 4.5...24 dB in relation to the reception of signals without coding.
Full Text
ВВЕДЕНИЕ
Посимвольный прием цифровых сигналов минимизирует вероятность ошибки на информационный бит Pб в отличие от приема, реализующего правило максимального правдоподобия, минимизирующее вероятность ошибки на дискретное сообщение [1–3]. Разработанные алгоритмы посимвольного приема сигнальных конструкций на основе ряда корректирующих кодов (например, на основе низкоплотностных кодов, турбокодов [3–7]) обеспечивают достижение вероятностных характеристик, близких к предельным характеристикам, определяемых пропускной способностью каналов с аддитивным белым гауссовским шумом (АБГШ) [6, 8–12].
Суть посимвольного приема – принятие решений относительно переданных сигнальных символов на основе вычисленных апостериорных вероятностей [2, 3].
Современной тенденцией при разработке информационных систем является использование цифровых сигналов со сложными «созвездиями», определяющими повышение их информативной емкости [1, 5, 13, 14]. Этот подход дает возможность увеличения скорости передачи информации при ограниченной частотной полосе канала [1, 4, 15]. В работах [9, 12–14, 16–18] приведены алгоритмы оптимального посимвольного приема сигнальных конструкций на основе цифровых сигналов с различными манипуляциями (например, с использованием многоуровневой фазовой манипуляции (ФМ-М)) в сочетании с корректирующими кодами в недвоичных полях Галуа GF (2m), объем которых согласован с объемом 2m соответствующих сигнальных «созвездий». Исследование помехоустойчивости этих алгоритмов приема ряда анализируемых сигнальных конструкций произведено для канала АБГШ [12, 14, 19].
Актуальной является проблема исследования помехоустойчивости алгоритмов оптимального посимвольного приема сигнальных конструкций из этого класса на основе корректирующего кодирования в недвоичных полях Галуа GF (2m) при передаче по каналам с амплитудными замираниями за счет многолучевого распространения [1, 4].
ПОСТАНОВКА ЗАДАЧИ
Пусть – последовательность k информационных символов поля Галуа GF (2m), формируемого по модулю неприводимого многочлена γ (x) степени m [15]. Элементы ai представляются многочленами
с коэффициентами αp (ai) ∈ GF (2) [5, 15].
Кодовое слово , соответствующее вектору корректирующего кода c кодовой скоростью r = k / n, задается соотношением , где H – порождающая матрица кода размером n × k, n – длина кодовых слов [15]. Кодовые символы bi (0 ≤ i ≤ n − 1) в составе однозначно сопоставляются цифровым сигналам с манипуляционным «созвездием» объемом 2m, передаваемым по каналам с помехами и искажениями.
Реализация с выхода демодулятора сигналов поступает на вход решающего устройства, – комплексные отсчеты квадратурных каналов. При использовании оптимального посимвольного приема вычисляются апостериорные вероятности , где β ∈ GF (2m), и принимаются решения относительно символов âi с использованием правила [2, 3]
. (1)
В [2, 12] приведено описание разработанного производительного алгоритма вычисления соотношения (1). Исследование помехоустойчивости этого алгоритма произведено для ряда цифровых сигналов с объемами «созвездий» 2m для канала АБГШ [12, 14, 17]. Показано, что применение анализируемого алгоритма приема сигнальных конструкций на основе цифровых сигналов ФМ-М и корректирующего кода с проверкой на четность в поле Галуа GF (2m) (m = 2, 3, 4) обеспечивает значимый энергетический выигрыш до 2 дБ по отношению к передаче без кодирования [12].
Цель работы – исследование помехоустойчивости разработанного алгоритма оптимального посимвольного приема анализируемых сигнальных конструкций на основе корректирующего кодирования в полях Галуа GF (2m) при многолучевом распространении по трансионосферным каналам (спутниковым ионосферным каналам) с замираниями за счет рассеяния на ионосферных неоднородностях.
АЛГОРИТМ ОПТИМАЛЬНОГО ПОСИМВОЛЬНОГО ПРИЕМА
Цифровые сигналы s (t) характеризуются частотной эффективностью m = log2M (бит/с/Гц), задающей максимальную скорость передачи R = mΔF (бит/с) для канала с частотной полосой ΔF [1, 4].
Сигнальные конструкции на основе s (t) и корректирующих кодов с параметрами (n, k) в полях GF (2m) задаются соотношением [1]
, (2)
где f – несущая частота; U (t) = 1 при iT ≤ t < (i + 1) T, иначе, U (t) = 0; Ai, φi – амплитуда и фаза радиоимпульсов в составе s (t), задающие вид сигнального «созвездия»; T – тактовый интервал.
Радиоимпульсы в составе s (t) формируются, отображая m двоичных информационных символов αp (bi) (p = 0, 1, ..., m − 1) в отсчеты «созвездия» с комплексными огибающими [4]. ФМ-М-сигналы характеризуются постоянством амплитуд Ai = Ac.
Апостериорные вероятности для символов сигнальных конструкций (2) вычисляются с использованием соотношения [2, 10]
, (3)
где – функция правдоподобия;
.
Сложность вычисления (3) определяется требуемым объемом вычислительных операций P1 ≈ 2mk и представляет трудоемкую задачу даже для малых значений m, k.
В работах [2, 12] приведено описание разработанного алгоритма для решения этой задачи с существенно более низкой сложностью реализации. Для пояснения ниже дано его общее описание. Алгоритм включает три этапа.
Первый этап. Вычисляется множество спектральных составляющих Cl (r) для последовательности
, (4)
где i = 0, 1, ..., 2m − 1; l = 0, 1, ..., n − 1 – номера кодовых символов в составе кодового слова; wi (r) – базисные функции Уолша-Адамара с перемежением номеров . Закон перемежения задается аналитически для порождающих многочленов в виде γ (x) = 1+ xk + xm (1 ≤ k ≤ m − 1) либо в виде таблицы, получаемой путем предобработки с целью сопоставления функциям Уолша–Адамара значений , задаваемых моделью канала передачи [2, 12].
Второй этап. Вычисляется множество {Tl (λ)} с использованием Cl (r) и множества кодовых слов R дуального кода CH с параметрами (n, n − k) [15]
. (5)
Обозначение rp : R ∈ CH в (5) определяет операции для последовательности кодовых символов rp в составе кодовых слов R кода CH.
Третий этап. Вычисляют апостериорные вероятности , используя алгоритм быстрого спектрального преобразования в базисе Уолша–Адамара (БПУ) размерностью 2m над {Tl (λ)},
. (6)
Решения относительно символов кодового слова принимаются на основе значений в соответствии с правилом (1).
При исполнении первого этапа (4) полагается, что на временном интервале lT ≤ t < (l + 1) T нормированные отсчеты для ФМ-М-сигналов и канала АБГШ с спектральной плотностью мощности N0 представляют случайные величины со средними
.
и с единичной дисперсией. Здесь Ec – энергии сигнальных символов «созвездия» с огибающими , плотности вероятности при вычислении соотношения (4) задаются как [12]
, (7)
||x|| – евклидова метрика; K – коэффициент нормировки.
Сложность вычисления соотношений (4)–(6) оценивается как P2 ≅ 2m (n – k), для значений n − k << k выполняется условие P2 << P1.
Наиболее простым для реализации является алгоритм посимвольного приема для сигнальных конструкций на основе корректирующего кода с проверкой на четность с параметрами (n, n − 1). В этом случае кодовый вектор кода CH представляет последовательность длительностью k + 1 одинаковых элементов поля GF (2m) и справедливо соотношение P2 ≅ 2m [9, 12, 19].
СТАТИСТИЧЕСКИЕ МОДЕЛИ КАНАЛА С ЗАМИРАНИЕМ
Замирания сигналов приводят к деградации вероятностных характеристик Pб по отношению к распространению в свободном пространстве [2]. Количественные значения соответствующих энергетических потерь определяются статистическими свойствами амплитуды сигналов A как случайного процесса.
Ниже рассматриваются модели замираний сигналов при их многолучевом распространении по трансионосферным каналам [20–24]. При создании и развитии моделей замираний сигналов для данных каналов используются два подхода – на основе аналитического описания процесса распространения сигналов [21, 22] и на основе использования эмпирических соотношений относительно плотности распределения p (A) [24].
Аналитические подходы основаны на решении стохастического нелинейного волнового уравнения относительно электрического поля [22]. Соответствующие решения, полученные с использованием борновского приближения и приближения Рытова, дают возможность оценить вторые статистические моменты функционалов от амплитуды A.
Модели замираний из второго класса связывают параметры эмпирической плотности распределения p (A) с индексом сцинтилляции [21–24]
,
здесь <·> – усреднение по времени, если полагать процесс A эргодическим.
Относительно значений S4 замирания классифицируются как слабые для S4 < 0.3, средние для 0.3 < S4 < 0.6 и сильные для S4 > 0.6 [21].
Вероятность ошибки Pб для ФМ-М-сигналов с амплитудой A имеет вид [1]
. (8)
Здесь p (θ) – плотность распределения фазы, для АБГШ с односторонней спектральной плотностью N0 справедливо соотношение [1]
, (9)
Ec = Eбlog2M – энергия радиоимпульсов; A – энергия на 1 бит.
Для амплитуды A с плотностью распределения p (A) вероятность ошибки Pб с учетом замираний задается соотношением [1]
. (10)
Для слабых замираний плотность распределения p (A) задается логнормальным законом [20, 21]
. (11)
Здесь A0 – амплитуда регулярной сигнальной составляющей; σ2 – средняя мощность многолучевых компонент сигналов на выходе ионосферной линии.
Более общее выражение для задается распределением Накагами [20]
. (12)
Здесь – средняя мощность сигналов; m ≥ 1/ 2 – параметр, задаваемый соотношением
,
; Г (m) – гамма-функция.
Параметры S4 и m для распределения Накагами связаны [20–22]
. (13)
Для m > 1 (для слабых, средних и, отчасти, сильных замираний) распределение Накагами аппроксимируется распределением Релея–Райса [21]
, (14)
где A0 – амплитуда регулярной сигнальной составляющей; I0 (x) – модифицированная функция Бесселя первого рода нулевого порядка.
Вычисление значений Pб с использованием соотношения (10) с целью оценивания энергетических потерь с учетом замираний по отношению к распространению в свободном пространстве выполняется при условии равенства средних энергий
.
Это условие обеспечивается выбором значений A0 и σ2 модели (14), как решения уравнения
. (15)
Значения в (15) для распределения Релея–Райса (14) связаны через коэффициент Райса d [25]
. (16)
Решения системы уравнений (15), (16) относительно имеют вид
.
РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ
Ниже даны результаты моделирования алгоритма посимвольного приема с целью исследования помехоустойчивости сигналов ФМ-4, ФМ-8 и ФМ-16 (коэффициенты частотной эффективности 2, 3 и 4 бит/с/Гц) в сочетании с корректирующим кодом с проверкой на четность в полях Галуа GF (22), GF (23), GF (24) при распространении по трансионосферным каналам с замираниями. Информационный объем передаваемых сообщений составляет 96 битов, т.е. для рассматриваемых сигнальных конструкций количество информационных символов k = 48, 32 и 24, соответствующие кодовые скорости корректирующих кодов равны 48/49, 32/33 и 24/25. Моделирование произведено для условий идеальной синхронизации по частоте, фазе и временным тактам. Приведены оценки вероятностей ошибки Pб и соответствующих энергетических потерь ΔE по отношению к распространению в свободном пространстве.
На рис. 1 приведены вероятности ошибки Pб, вычисленные с использованием соотношений (8), (9), (10) для сигналов ФМ-4, ФМ-8, ФМ-16 при передаче по каналу АБГШ: для вероятности Pб = 10–5 требуются значения Eб / N0, равные 9.5, 13.0 и 17.25 дБ соответственно.
Рис. 1. Вероятности ошибки Pб для сигналов ФМ-4 (кривая 1), ФМ-8 (кривая 2), ФМ-16 (кривая 3) при распространении по каналу АБГШ.
На рис. 2–4 приведены вероятностные кривые Pб для посимвольного приема рассматриваемых сигнальных конструкций – кривые 1 и 2 соответствуют вероятностям Pб без применения и с применением корректирующего кодирования для каналов АБГШ с замиранием с параметрами сцинтилляции трансионосферного канала S4 = 0.3 (см. рис. 2а, 3а, 4а) и S4 = 0.6 (см. рис. 2б, 3б, 4б). Для канала с замиранием с параметром S4 = 0.3 вероятность ошибки Pб = 10–5 при посимвольном приеме сигналов без использования корректирующего кодирования достигается при значениях Eб / N0, равных 16.0, 19.0 и 23.0 дБ соответственно (энергетические потери по отношению к АБГШ каналу без замираний составляют 6.5, 6.0, 5.75 дБ). Для параметра данная вероятность ошибки достигается при значениях 39.5, 43.0 и 45.0 дБ соответственно (энергетические потери по отношению к АБГШ каналу без замираний составляют 30.0, 30.0 и 27.75 дБ).
Рис. 2. Вероятности ошибки при приеме сигнальной конструкции на основе сигналов ФМ-4 для канала с замиранием (параметр S4 = 0.3 (а) и 0.6 (б)): 1 – без кодирования; 2 – с использованием корректирующего кода в поле GF (22).
Рис. 3. Вероятности ошибки при приеме сигнальной конструкции на основе сигналов ФМ-8 для канала с замиранием (параметр S4 = 0.3 (а) и 0.6 (б)): 1 – без кодирования; 2 – с использованием корректирующего кода в поле GF (23).
Рис. 4. Вероятности ошибки при приеме сигнальной конструкции на основе сигналов ФМ-16 для канала с замиранием (параметр S4 = 0.3 (а) и 0.6 (б)): 1 – без кодирования; 2 – с использованием корректирующего кода в поле GF (24).
Для канала с замиранием с параметром S4 = 0.3 вероятность Pб = 10–5 при посимвольном приеме сигналов с использованием корректирующего кода с проверкой на четность достигается при значениях Eб / N0, равных 11.5, 14.25 и 18.0 дБ соответственно (энергетические потери по отношению к АБГШ каналу без замираний составляют 2.0, 1.25 и 0.75 дБ). Для параметра S4 = 0.6 вероятность ошибки Pб = 10–5 достигается при значениях Eб / N0, равных 27.0, 19.0 и 30.0 дБ соответственно (энергетические потери по отношению к АБГШ каналу без замираний составляют 17.5, 6.0 и 12.75 дБ).
В табл. 1 приведены результирующие значения энергетического выигрыша при использовании рассматриваемых цифровых сигналов ФМ-М и корректирующего кодирования в полях Галуа по отношению к передаче без кодирования по трансионосферным каналам с замираниями. Видно, что применение корректирующего кодирования с проверкой на четность в сочетании с алгоритмом оптимального посимвольного приема анализируемых сигнальных конструкций на основе цифровых сигналов ФМ-М для трансионосферных каналов с замираниями обеспечивает значимый энергетический выигрыш до 4.5…24 дБ по отношению к посимвольному приему без корректирующего кодирования.
Таблица 1. Значения энергетического выигрыша (дБ) при использовании сигнальных конструкций на основе сигналов ФМ-М и корректирующего кодирования с кодовой скоростью r в полях Галуа GF (M) по отношению к передаче без кодирования по трансионосферным каналам
Цифровые сигналы | r | S4 = 0.3 | S4 = 0.6 |
ФМ-4 | 48/49 | 4.5 | 12.5 |
ФМ-8 | 32/33 | 4.75 | 24.0 |
ФМ-16 | 24/25 | 5.0 | 15.0 |
ЗАКЛЮЧЕНИЕ
Дано описание алгоритма оптимального посимвольного приема сигнальных конструкций на основе цифровых сигналов и корректирующего кодирования в полях Галуа GF (2m). Основу анализируемого алгоритма составляет БПУ с размерностью 2m. Сложность алгоритма посимвольного приема определяется размерностью дуального кода, что обусловливает перспективность его применения для блоковых корректирующих кодов с низкой избыточностью.
Исследование вероятностных характеристик рассматриваемого алгоритма посимвольного приема произведено путем его моделирования для сигнальных конструкций на основе интенсивно используемых в приложениях цифровых сигналов с многофазовой манипуляцией ФМ-4, ФМ-8 и ФМ-16 и для корректирующих кодов с проверкой на четность с кодовыми скоростями 48/49, 32/33 и 24/25 в полях Галуа GF (22), GF (23), GF (24). Показано, что применение алгоритма посимвольного приема для моделей трансионосферных каналов с замираниями за счет рассеяния на ионосферных неоднородностях обеспечивает значительный энергетический выигрыш, до 4.5…24 дБ, по отношению к приему сигналов без кодирования. Используемые модели каналов с замираниями задавались параметром индекса сцинтилляции S4 = 0.3 (слабые и средние замирания) и S4 = 0.6 (средние и сильные замирания).
ФИНАНСИРОВАНИЕ РАБОТЫ
Работа выполнена в рамках государственного задания ИРЭ им. В.А. Котельникова РАН.
About the authors
L. E. Nazarov
Fryazino Branch Kotelnikov Institute of Radioengineering and Electronics RAS
Author for correspondence.
Email: levnaz2018@mail.ru
Russian Federation, Vvedenskogo Squar., 1 Fryazino, Moscow region, 141190
References
- Proakis J.G., Salehi M. Digital communication. Boston: McGraw-Hill, Higher Education, 2001.
- Bahl L.R., Cocke J., Jelinek F., Raviv J. // IEEE Trans. 1974. V. IT-20. № 3. P. 284.
- Смольянинов В.М., Назаров Л.Е. // РЭ. 1999. Т. 44. № 7. С. 838.
- Скляр Б. Цифровая связь. Теоретические основы и практическое применение. М.: ИД «Вильямс», 2003.
- Johnson S.J. Iterative Error Correction: Turbo, Low-Density Parity-Check and Repeat-Accumulate Codes. Cambridge: Univ. Press, 2010.
- Назаров Л.Е., Головкин И.В. // РЭ. 2010. Т. 55. № 10. С. 1193.
- Терешонок М.В., Кленов Н.В., Лобов Е.М. и др. // РЭ. 2022. Т. 67. № 3. С. 294. https://doi.org/10.31857/S0033849422030160.
- Назаров Л.Е., Батанов В.В. // РЭ. 2022. Т. 67 № 8. С. 782. https://doi.org/ 10.31857/S0033849422080137.
- Ping Li, Chan S., Yeng K.L. // Electronic Lett. 1997. V. 33. № 19. P. 1614.
- Steiner F., Bocherer G., Liva G. // IEEE Comm. Lett. 2018. V. 22. № 11. P. 2210.
- Lin S.-J. // IEEE Trans. 2018. V. СOM-66. № 8. P. 3235.
- Назаров Л.Е. // РЭ. 2023. Т. 68. № 9. С. 873. http://doi.org/ 10.31857/S003384942309019X
- Bourduge J., Poulliat C., Gadat B. // 2023 IEEE Int. Symp. on Information Theory (ISIT). Taipei. 25-30 Jun. N.Y.: IEEE, 2023. P. 2517. https://doi.org/10.1109/ISIT54713.2023.10206851.
- Назаров Л.Е. // Физ. основы приборостроения. 2022. Т. 11. № 3. С. 44. https://doi.org/10.25210/jfop-2203-044049.
- Питерсон У., Уэлдон Э. Коды, исправляющие ошибки. М.: Мир, 1976.
- Назаров Л.Е., Шишкин П.В. // РЭ. 2019. Т. 64. № 9. С. 910.
- Назаров Л.Е., Шишкин П.В. // Журн. радиоэлектроники. 2018. № 12. http://jre.cplire.ru/jre/dec18/10/text.pdf
- Kaipa K. // IEEE Comm. Lett. 2018. V. 22. № 11. P. 2210.
- Yeo S., Park I.-C. // IEEE Trans. 2018. V. IT-64. № 7. P. 5170.
- Crane R.K. // Proc. IEEE. 1977. V. 65. № 2. P. 180.
- Rino C.L. Theory of Scintillation with Applications in Remote Sensing. Hoboken: John Wiley & Sons, 2011.
- Ionospheric Propagation Data and Prediction Methods Required for the Design of Satellite Services and Systems. Recommendation ITU-R P.531-11. Geneva: Int. Telecommun. Union (ITU), 2012. 24 p.
- Назаров Л.Е., Смирнов В.М. // Журн. Радиоэлектроники 2020. № 11. http://jre.cplire.ru/jre/nov20/7/text.pdf https://doi.org/10.30898/1684-1719.2020.11.7.
- Назаров Л.Е., Батанов В.В. // РЭ. 2022. Т. 67. № 11. С. 1133. https://doi.org/10.31857/S0033849422110110
- Тихонов В.И. Статистическая радиотехника. М.: Сов. радио, 1966.
Supplementary files
