1

Федеральное государственное бюджетное образовательное учреждение

высшего образования
«МИРЭА – Российский технологический университет»

На правах рукописи

Лелюх Андрей Александрович

ПОВЫШЕНИЕ ПОМЕХОУСТОЙЧИВОСТИ ПРИЕМА СИГНАЛОВ С

КВАДРАТУРНОЙ АМПЛИТУДНОЙ МОДУЛЯЦИЕЙ В СИСТЕМАХ
ШИРОКОПОЛОСНОГО ДОСТУПА ДЛЯ МОБИЛЬНЫХ АБОНЕНТОВ

Специальность 2.2.13 Радиотехника, в том числе системы и устройства

телевидения

Диссертация на соискание ученой степени

кандидата технических наук

Научный руководитель:
д.т.н., профессор Куликов Г.В.

Москва 2022
 2

ВВЕДЕНИЕ .……………………………………………………………..………..5
1. ОБЗОР МЕТОДОВ И СРЕДСТВ ПРЕДОСТАВЛЕНИЯ
ШИРОКОПОЛОСНЫХ ИНФОРМАЦИОННЫХ УСЛУГ ДЛЯ МОБИЛЬНЫХ
АБОНЕНТОВ………………………………………...………………………...….13
1.1. Основные особенности систем спутникового широкополосного доступа для
аэромобильных абонентов ……………………………………………….….15
1.2. Гибридная сеть связи………………………………………..…………….…..17
1.3. Оценка энергетического потенциала спутниковых радиолиний в
авиационной системе связи………………………………………..……….…19
1.4. Основные особенности систем спутникового широкополосного доступа для
абонентов на железнодорожном транспорте ……………………………..22
1.5. Сигналы и помехи в каналах радиосвязи КЦС с наземной базовой сетью 4G-
LTE .. …………………..……………………………………...........................25
1.6. Обзор работ по оценке влияния дестабилизирующих факторов на
помехоустойчивость приема сигналов КАМ...………………….………….30
1.7. Выводы по разделу 1…………………………..……………………….…..….35
2. АНАЛИЗ ПОМЕХОУСТОЙЧИВОСТИ КОГЕРЕНТНОГО ПРИЕМА
СИГНАЛОВ КАМ ПРИ НАЛИЧИИ НЕФЛУКТУАЦИОННЫХ ПОМЕХ……...36
2.1. Методика определения вероятности ошибки когерентного приемника
сигналов КАМ при наличии нефлуктуационных помех..……………………..….36
2.2. Анализ вероятности ошибки когерентного приемника сигналов КАМ при
наличии гармонической помех………………………………………………………40
2.3. Анализ вероятности ошибки когерентного приемника сигналов КАМ при
наличии сканирующей помехи.. ……………………………………………………48
2.4. Анализ вероятности ошибки когерентного приемника сигналов КАМ при
наличии фазоманипулированной помехи……………………………………….…..54
2.5. Анализ вероятности ошибки когерентного приемника сигналов КАМ при
наличии ретранслированной помехи………………………………………….…….65
2.6. Сравнительная оценка влияния различных нефлуктуационных помех……..71
2.7. Выводы по разделу 2……………………………………………………….…….77
 3

3. ВЛИЯНИЕ ПОГРЕШНОСТЕЙ ФОРМИРОВАНИЯ ОПОРНЫХ

КВАДРАТУРНЫХ КОЛЕБАНИЙ НА ПОМЕХОУСТОЙЧИВОСТЬ
КОГЕРЕНТНОГО ПРИЕМА СИГНАЛОВ КАМ……………………………..…… 78
3.1. Исследование влияния погрешностей оценивания фазы несущей на
помехоустойчивость приема сигналов КАМ………………………………….…….78
3.2. Исследование влияния погрешностей тактовой синхронизации на
помехоустойчивость приема сигналов КАМ……………………………………….84
3.3. Исследование влияния амплитудного разбаланса квадратур на
помехоустойчивость приема сигналов КАМ……………………………………….89
3.4. Исследование влияния фазового разбаланса квадратур на помехоустойчивость
приема сигналов КАМ……………………………………………………..………....95
3.5. Сравнение влияния погрешностей синхронизации при приеме сигналов КАМ
и М-ФМ……………………………………………………..……………………..…100
3.6. Выводы по разделу 3 ………………………………….…..………………… 101
4. РАЗРАБОТКА АЛГОРИТМОВ БОРЬБЫ С НЕФЛУКТУАЦИОННЫМИ
ПОМЕХАМИ ПРИ ПРИЕМЕ СИГНАЛОВ КАМ………….……………………..103
4.1. Оптимальный прием сигналов КАМ на фоне гармонической помехи со
случайной начальной фазой……………………………………………………… 103
4.2. Адаптивная фильтрация нефлуктуационных помех при приеме сигналов
КАМ…………………………………………………………….…….……………. 113
4.3. Выводы по разделу 4…………………………………………………………...121
5. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ………………………………..122
5.1. Измерение отношения сигнал/шум в сети LTE на маршруте следования
железнодорожного состава………………………………….……………………..122
5.2. Экспериментальное исследование помехоустойчивости приема сигналов 16-
КАМ в присутствии гармонической помехи……………………….…………….125
5.3. Выводы по разделу 5………………………………………….……………… 130
ЗАКЛЮЧЕНИЕ ………………………………………………..…….……………..131
СПИСОК СОКРАЩЕНИЙ………...…………………………..…………………...134
СПИСОК ЛИТЕРАТУРЫ…………………………………………….……………..135
 4

ПРИЛОЖЕНИЕ 1. Алгоритм и программа моделирования радиоканала с

использованием сигнала 16-КАМ (когерентный прием)…………………….….146
ПРИЛОЖЕНИЕ 2. Программа моделирования радиоканала с использованием
сигнала 16-КАМ (некогерентная обработка помехи)…………………………...155
ПРИЛОЖЕНИЕ 3. Программа моделирования адаптивного фильтра
(фрагмент)………………………………………………………..…….…………... 165
ПРИЛОЖЕНИЕ 4. Акты внедрения и копия патента……………………..……..166
 5

ВВЕДЕНИЕ
В соответствии с Указом Президента РФ от 01.12.2016 № 642 «О стратегии
научно-технологического развития Российской Федерации» в ближайшие 10 - 15
лет приоритетами научно-технологического развития Российской Федерации
обозначены направления, которые позволят получить научные и научно-
технические результаты и создать передовые цифровые, интеллектуальные
производственные технологии, являющиеся основой инновационного развития
внутреннего рынка продуктов и услуг. При этом, развитие информационно-
телекоммуникационных систем является одним из приоритетных направлений
развития науки, технологий и техники в Российской Федерации.
Научные организации и производственные компании ряда зарубежных
стран в последнее время активно разрабатывают системы передачи данных,
обеспечивающие предоставление широкополосных информационных услуг
(включая доступ в Интернет) для мобильных абонентов. Это относится к
системам связи гражданской авиации, железнодорожного транспорта, а также
мобильным объектам иного гражданского назначения. При этом рассматриваются
два основных способа организации связи: с применением спутниковых
группировок и с использованием наземных сотовых сетей.
Основой разработки перспективных систем передачи информации является
использование цифровых методов обработки сигналов. К числу основных
преимуществ таких методов можно отнести высокую пропускную способность
каналов связи при лучших, чем у аналоговых методов, показателях качества, а
также, что немаловажно, простота исполнения приемо-передающей аппаратуры на
цифровой унифицированной элементной базе.
К числу хорошо зарекомендовавших себя и часто используемых во всем мире
в цифровых радиосистемах относятся многопозиционные сигналы с квадратурной
амплитудной модуляцией (КАМ, англ. QAM), которые имеют высокие
энергетические и спектральные свойства и за счет многопозиционности
обеспечивают высокую пропускную способность каналов связи. В частности,
 6

данные сигналы широко используются в беспроводных и сотовых сетях 4G-LTE при

использовании технологии OFDMA.
В монографиях Дж. Прокиса [1], П. Боккера [2], Fuqin Xiong [3] и других
исследователей проанализирована помехоустойчивость приема сигналов с
квадратурной амплитудной модуляцией для случая приема на фоне гауссовского
шума. Современные радиосистемы зачастую работают в более сложной
помеховой обстановке, когда кроме шумовой помехи на входе приемника
присутствуют и нефлуктуационные помехи, вызванные разными причинами
природного и технического характера. Необходимо так же отметить, что эти
помехи могут быть созданы и преднамеренно. Воздействие нефлуктуационных
помех приводит к существенному снижению качества связи, поэтому анализ их
влияния на помехоустойчивость приемников сигналов КАМ и разработка методов
борьбы с ними представляет собой важную научно-техническую задачу. Ее
решению для этих и других видов сигналов посвящены научные работы Тихонова
В.И., Харисова В.Н., Ярлыкова М.С., Савватеева Ю.И., Мартиросова В.Е.,
Рубцова В.Д., Куликова Г.В, Пестрякова А.В., Довбня В.Г., Вольского В.А.,
Выдолбина Ю.К., Савищенко И.В., Яценко С.Ю., Артеменко А.А., Мальцева А.А.
и других.
Все вышесказанное определяет актуальность данной диссертационной
работы, направленной на решение научно-технической задачи повышения качества
связи в радиосистемах, использующих сигналы с квадратурной амплитудной
модуляцией.
Целью работы является повышение помехоустойчивости приема сигналов
с квадратурной амплитудной модуляцией при наличии нефлуктуационных помех
и аппаратурных погрешностей в перспективных системах широкополосного
доступа для мобильных абонентов.
Для достижения поставленной цели необходимо решить следующие
основные задачи:
1. Определение наиболее распространенных видов помех, характерных для
цифровых радиоканалов передачи информации на транспорте.
 7

2. Анализ и сравнение помехоустойчивости демодуляторов дискретных

сигналов с квадратурной амплитудной модуляцией при воздействии
нефлуктуационных помех.
3. Анализ влияния на помехоустойчивость приема сигналов с квадратурной
амплитудной модуляцией погрешностей системы синхронизации и блока
квадратурного преобразования.
4. Синтез и анализ помехоустойчивости алгоритмов борьбы с
нефлуктуационными помехам при приеме сигналов с квадратурной амплитудной
модуляцией.
5. Натурные аппаратурные испытания и компьютерное моделирование
приемников сигналов с квадратурной амплитудной модуляцией при наличии в
радиоканале связи нефлуктуационных помех.
Объектом исследования являются приемники сигналов с квадратурной
амплитудной модуляцией при воздействии нефлуктуационных помех и при
наличии погрешностей работы узлов приемника.
Предметом исследования является помехоустойчивость приемников
сигналов с квадратурной амплитудной модуляцией при воздействии
нефлуктуационных помех и при наличии погрешностей работы узлов приемника.
Методы исследования. В диссертационной работе использовались методы
оптимального приема сигналов, статистической радиотехники, теории
потенциальной помехоустойчивости, теории вероятностей, компьютерного
моделирования и натурные аппаратурные испытания.
Достоверность материалов исследования подтверждается корректным
использованием математического аппарата и результатами как компьютерного
моделирования, так и практической проверки в ходе отладки узлов и испытаний
комплекса цифровой связи (АО «МНИИРС», ООО «ВизКом») на подвижных
объектах, подтверждающими теоретические результаты и хорошо
согласующимися с частными результатами других авторов, известными из
открытых научно-технических источников.
 8

Научная новизна работы:

1. Предложена методика определения вероятности ошибки когерентного
приемника многопозиционных сигналов с квадратурной амплитудной модуляцией
при наличии нефлуктуационных помех.
2. С использованием предложенной методики получены оценки вероятности
битовой ошибки для приемников сигналов с квадратурной амплитудной
модуляцией при наличии гармонической, сканирующей, ретранслированной и
фазоманипулированной помех с разными спектральными и энергетическими
характеристиками.
3. Оценено и проверено компьютерным моделированием влияние
погрешностей формирования опорных квадратурных колебаний на
помехоустойчивость когерентного приема сигналов с квадратурной амплитудной
модуляцией.
4. Разработан алгоритм оптимального приема сигналов с квадратурной
амплитудной модуляцией на фоне гармонической помехи с некогерентной
обработкой помехи и оценена его помехоустойчивость.
5. Предложен алгоритм адаптивной фильтрации нефлуктуационных помех,
использующий информацию о фазовой структуре сигнала с квадратурной
амплитудной модуляцией, и оценена его помехоустойчивость.
На защиту выносятся следующие основные положения диссертационной
работы:
1. Предложенная методика определения вероятности ошибки для
когерентного приема сигналов с квадратурной амплитудной модуляцией на
фоне нефлуктуационных помех позволяет оценить энергетические потери
при разных параметрах помех.
2. Допустимая фазовая погрешность несущей при формировании опорных
колебаний в когерентном демодуляторе сигналов с квадратурной
амплитудной модуляцией в зависимости от позиционности составляет от
10 до 30, погрешность тактовой синхронизации – от 2 до 5%. Допустимый
 9

разбаланс квадратур составляет: по амплитуде – от 3 до 10%, по фазе – от

1,50 до 30.
3. Оптимальный алгоритм приема сигналов с квадратурной амплитудной
модуляцией с некогерентной обработкой гармонической помехи позволяет за
счет использования дополнительной схемы формирования весовых
коэффициентов обеспечить при большой интенсивности помехи энергетический
выигрыш для 4-КАМ около 2 дБ, для 16-КАМ – более 5 дБ.
4. Адаптивный алгоритм подавления гармонической помехи, использующий
информацию о фазовой структуре сигнала с квадратурной амплитудной
модуляцией, позволяет обеспечить энергетический выигрыш для 4-КАМ
1-3 дБ, для 16-КАМ 0,5 дБ.
Практическая значимость работы:
1. Оценена помехоустойчивость приемников сигналов с квадратурной
амплитудной модуляцией при наличии в радиоканале связи различных
нефлуктуационных помех и определены возникающие энергетические потери.
2. Определены допустимые погрешности систем фазовой и тактовой
синхронизации и величины амплитудного и фазового разбаланса квадратур при
когерентном приеме сигналов с квадратурной амплитудной модуляцией
(использовано в АО «МНИИРС»).
3. Разработаны программы моделирования канала передачи информации при
наличии нефлуктуационных помех с использованием сигналов 16-КАМ.
4. Разработаны методы борьбы с нефлуктуационными помехами, в частности,
наиболее опасной гармонической помехой, при приеме сигналов с квадратурной
амплитудной модуляцией, а именно, синтезирован оптимальный приемник
сигналов на фоне гармонической помехи со случайной начальной фазой с
некогерентной обработкой помехи и предложен алгоритм адаптивной фильтрация
нефлуктуационных помех, следящий за фазовой траекторией сигнала КАМ
(патент на полезную модель RU 198305 U1, использовано в ООО «ВизКом»).
Апробация работы. Основные положения диссертации докладывались на
следующих научных конференциях:
 10

1. III Международная заочная научно-практическая конференция «Авиация:

история, современность, перспективы развития», Минск, БГАА, 8–9 ноября
2018 г.
2. IV Международная научно-практическая конференция «Актуальные
проблемы и перспективы развития радиотехнических и
инфокоммуникационных систем» «РАДИОИНФОКОМ-2019», Москва,
ноябрь 2019 г.
3. III Международная научно-практическая конференция «Актуальные
проблемы и перспективы развития авиации», Минск, БГАА, 22–23 мая 2019 г.
4. Международный симпозиум «Надежность и качество», Пенза, 2020.
5. 2020 Systems of Signal Synchronization, Generating and Processing in
Telecommunications (SYNCHROINFO), 1-3 July 2020, Svetlogorsk, Russia.
Публикации. По материалам диссертации опубликовано 19 работ, из них
10 статей – в рецензируемых журналах из Перечня ВАК, 1 статья в издании,
индексируемом в WoS и Scopus, 7 докладов в сборниках трудов российских и
международных конференций с 2018 по 2020 гг., 1 из которых индексирован в
Scopus. Получен 1 патент на полезную модель.
Личный вклад соискателя. Все основные научные результаты диссертации
получены лично автором, либо при его непосредственном участии. Автор
самостоятельно провел все расчеты, он активно участвовал в определении методов
исследования, в обсуждении результатов, в подготовке и проведении натурных
аппаратурных испытаний и компьютерного моделирования.
Исследование проводилось при финансовой поддержке Министерства науки
и высшего образования Российской Федерации (уникальный идентификатор
RFMEFI57418X0190). Результаты диссертационной работы использованы в АО
«МНИИРС» и ООО «ВизКом».
Структура и объем диссертации. Диссертация состоит из введения, 5
разделов, заключения, списка литературы и приложений. Объем работы
составляет 168 страниц, рисунков – 62, таблиц – 7. Список литературы включает
80 наименований.
 11

В первом разделе описаны методы и средства предоставления

широкополосных информационных услуг для мобильных абонентов –
авиапассажиров и пассажиров на железнодорожном транспорте. Приведена
информация об условиях приема радиосигналов в таких каналах связи, проведен
анализ литературных источников, посвященных оценке влияния
дестабилизирующих факторов на помехоустойчивость приема сигналов КАМ.
Во втором разделе приведена методика определения вероятности ошибки
когерентного приемника сигналов с квадратурной амплитудной модуляцией при
наличии нефлуктуационных помех и проведен анализ помехоустойчивости
приема для случаев, когда в радиоканале наряду с шумовой помехой
присутствуют следующие нефлуктуационные помехи: гармоническая,
сканирующая, фазоманипулированная и ретранслированная. По результатам
анализа оценены энергетические потери, вызванные этими помеховыми
воздействиями, и определены наиболее опасные виды помех.
В третьем разделе исследовано влияние погрешностей работы узлов
приемника на помехоустойчивость когерентного приема сигналов КАМ, а
именно, фазовой погрешности при формировании опорных колебаний и
погрешностей тактовой синхронизации, а также амплитудного и фазового
разбаланса квадратур. Это позволило оценить допустимые отклонения
параметров приемника от номинальных значений.
В четвертом разделе разработаны алгоритмы подавления нефлуктуационных
помех, в частности, наиболее опасной гармонической помехи, при приеме сигналов
КАМ. Синтезирован оптимальный приемник сигналов с квадратурной амплитудной
модуляцией на фоне гармонической помехи со случайной начальной фазой с
некогерентной обработкой такой помехи и предложен алгоритм адаптивной
фильтрация нефлуктуационных помех, следящий за фазовой траекторией сигнала
КАМ. Оценен выигрыш в помехоустойчивости предложенных методов борьбы с
нефлуктуационной помехой.
Пятый раздел посвящен экспериментальным исследованиям. Приводятся
результаты замеров отношения сигнал/шум на входе экспериментального модема
 12

в сети LTE на маршруте следования железнодорожного состава Москва – Санкт

Петербург и результаты оценки помехоустойчивости приема сигналов КАМ в
присутствии гармонической помехи на экспериментальном стенде.
В заключении обобщены и приведены результаты анализа
помехоустойчивости приемников сигналов КАМ при воздействии
нефлуктуационных помех, а также при наличии погрешностей работы узлов
приемника, дается оценка эффективности предложенных методов борьбы с
помехами и обозначаются направления дальнейших исследований.
В приложении приведены программы компьютерного моделирования
канала передачи информации с использованием сигналов 16-КАМ, акты об
использовании результатов диссертационной работы и копия документа на объект
интеллектуальной собственности.
 13

1. ОБЗОР МЕТОДОВ И СРЕДСТВ ПРЕДОСТАВЛЕНИЯ

ШИРОКОПОЛОСНЫХ ИНФОРМАЦИОННЫХ УСЛУГ ДЛЯ
МОБИЛЬНЫХ АБОНЕНТОВ

Предоставление широкополосных информационных услуг (включая

доступ в Интернет) для мобильных абонентов, пользующихся самолетами
гражданской авиации, железнодорожным и автомобильным транспортом
осуществляется в мире с применением спутниковых группировок и с
использованием наземных сотовых сетей.
Так в США и Канаде действует широкополосная система прямой связи
«воздух-земля» ATG (Air To Ground) в полосах частот 849-851 МГц/ 894-896
МГц (технология CDMA) и планируется развернуть вторую систему в полосе
частот 2315-2324,5 МГц. В Китае предполагается использовать технологию
CDMA EV-DO или LTE TDD в полосе частот 1785-1805 МГц.
В этих системах предоставление широкополосных услуг на борту
воздушного судна (включая доступ в Интернет) обеспечивается за счет
организации радиолиний между воздушным судном и сетью наземных
базовых станций (рис. 1.1).
Такое решение имеет «плюсы» и «минусы». Преимуществами по
сравнению со связью через спутник можно назвать:
- возможность развертывания в короткое время и низкую стоимость
бортового оборудования;
- малую временную задержку сигнала.
К недостаткам системы можно отнести невозможность использовать ее
над морями и океанами, а также в труднодоступных районах.
 14

Бортовое оборудование DA2GC

Наземные
станции DA2GC

Радио-
линии Транс-
DA2GC портная
сеть

Наземные
станции DA2GC

Рис. 1.1. Организация радиолиний между воздушным судном и сетью

наземных базовых станций.

Примерами зарубежных систем спутниковой связи, реализующими

технологии широкополосного доступа мобильных абонентов, являются
Inmarsat Global Xpress, Intelsat EpicNG, Eutelsat KaSat, ViaSat.
Аэромобильные станции для широкополосного спутникового доступа
производят зарубежные компании Honeywell Aerospace JetWave, ViaSat,
Astronics AeroSat, Gilat и др. При этом пропускная способность одного
спутника класса HTS (HTS – High Throughput Satellite – высокая пропускная
способность) [4-7] составляет десятки-сотни Гбит/с, скорость передачи
информации в прямой линии – десятки-сотни Мбит/с для аэромобильного
терминала и единицы-десятки Мбит/с – для индивидуального авиапассажира.
Необходимо отметить, что в данных системах полностью отсутствует
место отечественного разработчика и, тем более, производителя. Для
развития отечественной отрасли связи и обеспечения независимости
необходимо создание собственных систем и оборудования, работающего на
российские спутниковые группировки.
В настоящее время существуют эксплуатируемые и планируемые к
запуску отечественные спутники связи, которые можно отнести к классу
HTS, например, Экспресс-АМ5, Экспресс-АМ6, Ямал-601.
 15

1.1. Основные особенности систем спутникового широкополосного

доступа для аэромобильных абонентов

Общая архитектура спутниковой системы передачи информации для

предоставления авиапассажирам спутникового широкополосного доступа в
Интернет представлена на рис. 1.2 [8,9].

Рис. 1.2. Архитектура построения комплексов цифровой связи для

летательных аппаратов:
1 – источник питания; 2 – обтекатель и антенна; 3 – модем; 4 – сервер; 5 –
беспроводные точки доступа; 6 – обратная линия; 7 – прямая линия; 8 – канал
«вверх»; 9 – канал «вниз»; 10 – шлюзовая земная станция; 11 – центральный
узел сети; 12 – центр данных; 13 – Интернет; 14 – спутник.

В комплекс цифровой связи для летательных аппаратов входят два

крупных сегмента: среда передачи информации через спутник между
шлюзовой земной станцией и самолетной абонентской станцией и среда
доступа к информации пассажиров самолета.
Первый сегмент состоит из двух линий: прямой – для передачи
информации от земной станции к самолету и обратной – для
противоположного направления передачи.
Второй сегмент на борту самолета служит для организации обмена
 16

информацией между сервером и мобильными устройствами пассажиров; при

этом могут быть использованы технологии Wi-Fi, стандарты сотовой связи
или авиационные стандарты, например, ARINC 429.
Для современного этапа развития и ближайшей перспективы
спутникового широкополосного доступа в интересах аэромобильных
абонентов характерно применение:
- Ku- и Ka-диапазонов радиоволн;
- геостационарных спутников;
- спутниковых ретрансляторов без демодуляции и обработки информации
на борту спутника;
- спутниковых транспондеров с полосой порядка 100 МГц и выше;
- многолучевых бортовых антенн с объединением парциальных
(локальных) зон обслуживания отдельных сфокусированных узких лучей.
- сотового принципа повторного использования частотного ресурса в
парциальных зонах обслуживания.
- «бесшовного» перехода с луча на луч.
- уплотнения каналов, схем модуляции и кодирования в прямой линии,
регламентированных стандартом DVB-S2/S2X [10-12].
Особенностью таких систем связи, связанной с высокой мобильностью
самолета, является то, что для обеспечения стабильной связи на борту
самолета необходимо использовать управляемые антенны, которые могут
отслеживать положение спутника относительно самолета. При этом
достижимые значения эквивалентной изотропно излучаемой мощности
(ЭИИМ) и добротности аэромобильных станций в Ka-диапазоне составляют
не менее 43-48 дБВт и 11-13 дБ/К, соответственно.
На сегодняшний день доступ пассажиров самолётов к широкополосным
мультимедийным услугам осуществляется, в основном, за счёт
международной глобальной системы спутниковой связи «Инмарсат» в L-
диапазоне и Ка-диапазоне. При этом пропускная способность обмена
информации в L-диапазоне очень низкая (декларируемая 492 кбит/с, а
 17

реально меньше), стоимость предоставления услуг очень высокая,

отсутствует возможность просмотра телевизионных каналов и проведения
видеоконференцсвязи, широкого применения возможностей IP-телефонии и
интернета. В Ка-диапазоне скорость обмена информации достаточно
высокая, однако, в системе нет ни одного блока или узла, разработанного
российскими специалистами. Нашей стране, как правило, предлагается
только использовать уже развёрнутые системы связи.
Одним из немногоих предприятий в Российской Федерации, которое
устанавливается на самолётах земные станции спутниковой связи (ЗССС) Ku-
диапазона радиоволн и другое отечественное оборудование для предоставления
пассажирам услуг IP-телефонии, высокоскоростного интернета,
видеоконференцсвязи, приёма телевизионного сигнала высокой чёткости
является АО «Московский научно-исследовательский институт радиосвязи»
(АО «МНИИРС»).

1.2. Гибридная сеть связи

Новым направлением разработок (АО «МНИИРС») является создание

гибридной сети связи (рис. 1.3), основанной на одновременном
использовании систем наземно-воздушной связи в формате LTE и системы
спутниковой связи на базе отечественных космических группировок.

.
Рис. 1.3. Гибридная сеть связи.
 18

В гибридную сеть связи входят:

1) спутниковая радиолиния передачи информации между самолетной
абонентской станцией и спутником;
2) радиолиния передачи информации между самолетной абонентской
станцией и сетью наземных сотовых базовых станций.
Это означает, что частично покрытие будет обеспечиваться спутниками, а
основную роль будут играть наземные базовые станции (ATG, Air-to-Ground)
с направленными вверх антеннами [13]. Именно по такому принципу
работает сеть Gogo в США; главным преимуществом такой технологии
является низкая себестоимость услуги в сравнении со спутниковыми сетями
и высокая скорость передачи данных.
Именно разработка гибридной сети актуальна для России, где имеют
место, как густозаселённые области в Европейской части страны, так и
огромные малозаселённые пространства Сибири, Дальнего Востока и
Северных территорий.
Основным элементом бортовой (самолётной) части гибридной сети
связи будет являться комплекс цифровой связи (КЦС), представляющий
собой совокупность оборудования, работающего с наземной базовой сетью
4G-LTE со скоростью обмена информации до 20 Мбит/с и бортового
оборудования спутниковой связи в Ku-диапазоне со скоростью передачи до 2
Мбит/с.
Комплекс цифровой связи предназначен для обеспечения
широкополосного доступа к мультимедийным услугам авиационным
пассажирам, а также получения экипажем летательного аппарата
информации в интересах управления воздушным движением.
Система предназначена для предоставления следующих услуг:
- доступ в Интернет, в том числе организация Wi-Fi;
- организация IP-телефонии;
- организация видеоконференцсвязи;
- просмотр телевидения;
 19

- прослушивание радиостанций;
- организация GSM-связи;
- обеспечение пассажирам просмотра коллекции видеофильмов,
прослушивание музыки, чтение книг, доступа к компьютерным играм;
- предоставление услуг штатному экипажу (обмен аэронавигационной
информацией, ACARS, CPDLC);
- передача мониторинговой информации о подвижном объекте, работе
его основных систем и оборудования.

1.3. Оценка энергетического потенциала спутниковых радиолиний в

авиационной системе связи

Энергетический потенциал прямой линии (в Гц) определяется

отношением мощности полезного сигнала (в Вт) на входе приемника
аэромобильного терминала к односторонней спектральной плотности
мощности аддитивных помех (в Вт/Гц) [7,13]:

1
 C   C 1  C 1  C 1  C 1  C 
1 
              
 

 , (1.1)
 0  F  0 U  0  D  0  IM  0 CCI  0  ACI 
N  N N N N N

 C   C   C   C   C 
где   ,   ,   ,   ,   – парциальные энергетические
 N 0 U  N 0  D N
 0  IM  N 0 CCI N
 0  ACI

потенциалы (в Гц), определяемые соответственно шумами канала «вверх»,

шумами канала «вниз», интермодуляционными помехами в транспондере,
соканальными помехами вследствие повторного использования частотного
ресурса и междуканальными помехами. В корректно спроектированной
системе основным ограничивающим скорость передачи информации
фактором являются шумы канала «вниз», поэтому далее рассмотрим
соответствующий энергетический потенциал, который определяется
формулой [7]:
 20

 C  EIRPD (G / T ) MS
   ,
N
 0 D kL D

где EIRPD – парциальная ЭИИМ спутникового транспондера в канале «вниз»

для рассматриваемой аэромобильной станции, (G/T)MS – добротность
приемной системы аэромобильной станции, k – постоянная Больцмана, LD –
потери в канале «вниз».
На рис. 1.4 представлены рассчитанные зависимости достижимого
энергетического потенциала канала «вниз» (левая ось ординат) от
парциальной ЭИИМ для следующих условий:
- орбитальная позиция геостационарного спутника 49° в.д.;
- несущая частота 19,45 ГГц (диапазон Ka);
- поляризация круговая левая;
- добротность приемной системы аэромобильной станции 13 дБ/К.

Рис. 1.4. Энергетический потенциал и символьное отношение сигнал/шум

в канале «вниз» в зависимости от парциальной ЭИИМ.
 21

Верхний график соответствует крейсерскому полету самолета. Нижние

графики приведены для сравнительной иллюстрации влияния осадков в Ka-
диапазоне и соответствуют приему спутникового сигнала на земную станцию,
расположенную на территории Сунженского района Республики Ингушетия,
где среднегодовая интенсивность дождей одна из наиболее высоких в России;
при этом были заданы два значения процента потери надежности по дождю Tr
: 0,5% и 0,1%. При расчете учитывались потери в атмосфере, а также
уменьшение добротности приемной системы из-за дождя в соответствии с
методиками, изложенными в рекомендациях ITU-R Р.372, Р.453, Р.531, Р.618,
Р.676, Р.678, Р.834, Р.836, Р.837, Р.838, Р.839, Р.840, Р.841, Р.1511.
Теперь в качестве примера рассмотрим возможность использования
данного энергетического потенциала для обеспечения скорости передачи
информации Rb=50 Мбит/с в канале «вниз» с использованием технологий
стандарта DVB-S2X [7-9]. Выберем из числа возможных в стандарте схему
модуляции-кодирования «256APSK-3/4», обеспечивающую наиболее
высокую спектральную эффективность, равную 5,9 бит/символ [10]. Тогда
для этой схемы требуемое символьное отношение сигнал/шум,
соответствующее предусмотренной стандартом вероятности кадровой
ошибки FER=10-5, должно составлять не менее Es/N0=19,57 дБ. Фактически
достижимое символьное отношение сигнал/шум в канале «вниз» (в разах)
связано с энергетическим потенциалом и информационной скоростью для
используемой схемы модуляции-кодирования следующей зависимостью:
 C 
  (3 / 4) log 2 (256)
 Es   N0  D
   .
 N0  D Rb

Зависимости символьного отношения сигнал/шум в канале «вниз» от

парциальной ЭИИМ соответствуют правой оси ординат на рис. 1.4. Здесь же
пунктирной линией показано требуемое пороговое отношение сигнал/шум. Из
рисунка видно, что превышение этого порогового отношения для
крейсерского полета самолета обеспечивается при парциальной ЭИИМ более
 22

57,62 дБВт. При этом требуемый энергетический запас и, соответственно,

требуемое рабочее значение парциальной ЭИИМ могут быть заданы на
последующем этапе детального анализа всех слагаемых в (1.1). При влиянии
дождя требуемые значения ЭИИМ увеличиваются, причем, чем меньше
допустимый процент потери надежности по дождю, тем большее рабочее
значение ЭИИМ необходимо.
Для повышения энергетической эффективности системы возможно
использование других схем модуляции-кодирования, предусмотренных
стандартом DVB-S2X и требующих меньшего порогового отношения сигнал-
шум, но в обмен на ухудшение спектральной эффективности [10].
На рис. 1.5 в качестве еще одного примера показаны энергетический
потенциал и символьное отношение сигнал-шум в зависимости от
добротности приемной системы при парциальной ЭИИМ 57 дБВт и
использовании схемы модуляции-кодирования «32APSK-7/9» (при прочих
неизменных условиях).
Данная схема обеспечивает спектральную эффективность 3,84
бит/символ и требует символьного отношения сигнал/шум не менее
Es/N0=13,05 дБ.

1.4. Основные особенности систем спутникового широкополосного

доступа для абонентов на железнодорожном транспорте

Существуют проекты (АО «МНИИРС») создания интегрированной

системы связи для обеспечения широкополосного доступа к
мультимедийным услугам для железнодорожного транспорта.
Комплекс цифровой связи (КЦС) предназначен для обеспечения
широкополосного доступа к мультимедийным услугам пассажиров на
подвижных объектах с возможностью одновременной работы через ресурсы
отечественных космических аппаратов (КА), а также наземных пунктов
привязки 3G/4G-LTE «МегаФон», размещённых вдоль маршрута движения.
 23

Рис. 1.5. Энергетический потенциал и символьное отношение сигнал/шум

в канале «вниз» в зависимости от добротности приемной системы.

Система предназначена для предоставления пассажирам всех

вышеприведенных услуг. При этом спутниковый канал связи строится
аналогично авиационной системе, а связь между вагонами может
осуществляться по кабелю или с применением Wi-Fi технологии (рис. 1.6).
Проведенные в АО «МНИИРС» испытания позволили провести
измерения параметров приема спутникового сигнала транспондеров НТВ-
Плюс (гиростабилизированная антенна KVH-Tracvision TV6, апертура 60 см)
со спутника Express AMU-1 и спутника Eutelsat 36B в позиции 36 град. в.д.
Измерения уровня сигнала, отношения сигнал/шум (SNR) проводились
на участках прямой видимости спутника при движении и на стоянке
подвижного состава.
 24

Рис. 1.6. Структура обеспечения широкополосного доступа пассажиров

железнодорожного транспорта.

Таблица 1.1. Результаты измерений параметров при приеме спутникового

сигнала
Точка измерения
выход мультисвитча Terra MR-508
№ Наименование
Частота Стоянка Прямое Обратное
п/п спутника
направление направление
Уровень SNR Уровень SNR Уровень SNR
1 Eutelsat 36B 11785 R -69,61 8,54 -69,7 8,37 -69,64 8,17
2 Eutelsat 36B 11977 R -67,53 8,42 -67,7 8,2 -67,62 7,95
3 Express AMU-1 12284 R -68,23 8,57 -68,51 8,32 -68,4 8,23
4 Express AMU-1 12265 L -65,43 12,35 -65,5 12,1 -65,74 12,1
5 Express AMU-1 12341 L -64,31 11,72 -64,51 11,53 -64,82 11,54

Достигнутая при этих условиях вероятность ошибки на бит составила

10-6.
В разделе 5 (табл. 5.1) приведены результаты замеров отношения
сигнал/шум на входе модема в сети LTE на маршруте следования
железнодорожного состава Москва – Санкт Петербург. Показано, что
величина отношения сигнал/шум может изменяться в диапазоне 0 – 27 дБ.
 25

1.5. Сигналы и помехи в каналах радиосвязи КЦС с наземной базовой

сетью 4G-LTE

Как было указано выше, основным элементом бортовой (самолётной)

части гибридной сети связи и аналогичной части сети связи железнодорожного
состава будет комплекс цифровой связи (КЦС), представляющий собой
совокупность оборудования, работающего с наземной базовой сетью 4G-LTE. В
этом случае при использовании технологии OFDMA [14,15] спектр разбивается
на множество ортогональных поднесущих, общее количество которых
ограничено допустимой шириной спектра сигнала. Каждая из поднесущих
может иметь свой вид модуляции. При этом и в нисходящем канале, и в
восходящем канале используется квадратурная амплитудная модуляция (КАМ,
англ. QAM) разной позиционности – 4-КАМ (QPSK), 16-КAM, 64-КAM (в
зависимости от качества канала – до 256-КАМ).
Сигнал многопозиционной квадратурной амплитудной модуляции на
тактовом интервале, равном длительности канального символа Ts , несущего
информацию о k= log2 M информационных битах, может принимать одно из М
возможных значений:
si (t )  Aср ( Ii cos ω0t  Qi sin ω0t ), t  (0, Ts ], i  0,1,...M 1, (1.2)

где Aср  2 Esср / Т s – средняя амплитуда сигнала ( Esср  Ebср log2 M – энергия

канального символа, средняя по всем информационным комбинациям; Ebср –

средняя энергия, приходящаяся на один бит информации), Ii и Qi –
коэффициенты, определяющие амплитуды квадратурных компонент сигнала;
ω0 – несущая частота.
Примеры сигнальных созвездий сигналов КАМ представлены на рис.
1.7. На них отмечены сигнальные точки, соответствующие канальным
символам, и пунктиром изображены границы областей принятия решения.
При воздействии различных помех происходит смещение сигнальной точки в
 26

«неправильную» область, что вызывает ошибку приема канального символа

и, соответственно, одного или нескольких бит.

10111 10011 00110 00010

1101 1001 0001 0101 10010 10101 10001 00100 00101 00111
01 11
1100 1000 0000 0100 10110 10100 10000 00000 00001 00011

1110 1010 0010 0110 11011 11001 11000 01000 01100 01110

1111 1011 0011 0111 11111 11101 11100 01001 01101 01010
00 10
х
11010 11110 01011 01111
3х

М=4 М=16 М=32

Рис. 1.7. Сигнальные созвездия многопозиционной квадратурной

амплитудной модуляции.
Исследования помехоустойчивости приема сигналов КАМ [1-3] на фоне
белого гауссовского шума показали, что минимально достижимая
вероятность ошибки приема канального символа при оптимальной
когерентной обработке такого сигнала для прямоугольных сигнальных
созвездий при четных k определяется выражением:


3 1
 dt.
 t 2 /2
Pe КАМ  4Q( kEbср / N 0 , Q( x)  e (1.3)
M 1 2 x

Здесь Ebср / N0 – среднее отношение сигнал/шум, пересчитанное на один

бит информации.
Графики зависимости вероятности битовой ошибки Peb от Ebср / N0 для
М=4, 16 и 64 показаны на рис. 1.8.
В табл. 1.2 для сравнения эффективности различных методов
многопозиционной манипуляции приведены теоретические значения
отношения сигнал/шум, требуемые для достижения величины Peb=10-10 и
числовой параметр, характеризующий эффективность использования полосы
радиоканала.
 27

Рис. 1.8. Графики зависимости вероятности битовой ошибки от отношения

сигнал/шум.

Таблица 1.2. Сравнение эффективности методов цифровой манипуляции

Метод манипуляции Отношение сигнал/шум, Эффективность
Ebср / N0 , дБ использования полосы
радиоканала, бит/с/Гц
2-ФМ (BPSK) 13,06 1,0
4-ФМ, 4-КAM (QPSK) 13,06 2,0
8-ФМ 16,55 3,0
16-ФМ 21,09 4,0
16-КAM 16,98 4,0
64-КAM 21,40 6,0

Видно, что сигналы КАМ при больших одинаковых значениях

позиционности обладают большей энергетической эффективностью, чем
сигналы с многопозиционной фазовой манипуляцией, но имеют
непостоянную огибающую, что ограничивает их применение в радиоканалах
с нелинейностью.
Для приема сигналов КАМ используются две основные схемы
когерентных демодуляторов – многоканальный демодулятор, принимающий
 28

решение по максимуму правдоподобия (рис. 1.9), и квадратурный

демодулятор (рис. 1.10).


Ts

 x(t )S
0
опM 1 (t )dt
t  Ts
SопM 1 (t)
..
.

Ts

 x(t )S опi (t )dt

t  Ts
Канальный
Sопi (t )
0
Блок выбора символ
x(t ) .. максимума
.


Ts

 x(t )S
0
оп1 (t )dt
t  Ts
Sоп1 (t )


Ts

 x(t )S
0
оп 0 (t )dt
t  Ts
Sоп 0 (t )
Sоп1 (t ) Sопi (t ) SопM 1 (t)
..
.

..
.

Блок формирования Блок тактовой

опорных колебаний синхронизации

Рис. 1.9. Структурная схема многоканального когерентного приемника

сигналов КАМ.

Заметим, что в схеме, приведенной на рис. 1.9, определение максимума

проводится с учетом порога принятия решения, определяемого
полуразностью энергий сравниваемых канальных символов m и i:
Esm  Esi Es ср
 mi   [( I m 2  Qm 2 )  ( I i 2  Qi 2 )] (1.4)
N0 N0
 29


Ts
Xi
 x(t )cos(ω
0
0 t)dt
cos(ω0 t)

Генератор Канальный
x(t ) несущей Yi Решающее символ
частоты Z i  acrtg( )
Xi устройство
π
2
sin(ω0 t)

Ts

 x(t )sin(ω
0
0 t)dt
Yi

Рис. 1.10. Структурная схема квадратурного демодулятора сигналов

КАМ.

В реальных радиоканалах, в том числе каналах систем спутникового

широкополосного доступа для аэромобильных абонентов наряду с шумовой
присутствуют и нефлуктуационные помехи, вызванные разными причинами
– естественными (условия распространения радиоволн) и искусственными
(погрешности аппаратуры, нарушение регламентов связи и т.п.). Их
появление приводит к снижению качества связи.
Нефлуктуационные помехи, чаще всего сигналоподобные, могут иметь
разную структуру, некоторую частотную расстройку относительно несущей
частоты полезного сигнала и уровни, меньшие или соизмеримые с его
уровнем. Воспользуемся моделями помех, описанными в [16]:
- гармоническая помеха:
sп (t )  μA0 cos[(ω0  Δωп )t  φп ] (1.5)
с относительной интенсивностью μ , частотной расстройкой Δωп и
случайной начальной фазой φп , равномерно распределенной на
полуинтервале (-, +]; этой моделью могут быть описаны, например,
излучения гетеродинов или сигналы, более узкополосные, чем полезный
сигнал;
- фазоманипулированная помеха:
 30

sп (t )  μA0a j cos[(ω0  ωп )t  φп ], t  (( j  1)Tп , jTп ], j  1,..., N , (1.6)

где a j  1 – случайный дискретный символ помехи; Tп  T N ; N –

относительная скорость манипуляции помехи; этой моделью могут быть

описаны, например, сигналы соседних радиостанций систем передачи
дискретной информации;
- ретранслированная помеха – отраженный от местных предметов и
задержанный на время τ полезный сигнал:
sп (t )  μsi (t  τ,φп ); (1.7)
- сканирующая помеха – «скользящая» по частоте узкополосная помеха с
девиацией ωд и периодом сканирования Tс :

2ωдt
sп (t )  μA0 cos(ωп (t )t  φп ), ωп (t )  ω0  ωд  . (1.8)
Tс
Этой моделью могут быть описаны, например, сигналы радиостанций с
линейной частотной модуляцией.
В качестве шумовой помехи будем использовать стационарный
случайный процесс типа «белый гауссовский шум» с нулевым средним
значением и корреляционной функцией вида [52]:
N0
R(τ)  δ(τ) , (1.9)
2
где N0 – односторонняя спектральная плотность мощности шума.

1.6. Обзор работ по оценке влияния дестабилизирующих факторов на

помехоустойчивость приема сигналов КАМ

Задачу оценки помехоустойчивости приема многопозиционных

сигналов с квадратурной амплитудной модуляцией решали и решают многие
исследователи. Классические формулы для случая приема на фоне белого
гауссовского шума получены в [1-3]. Подобные ситуации, но с применением
 31

других методик оценки вероятностей ошибок или алгоритмов формирования

и обработки сигналов, рассмотрены в [17-19].
Большой интерес представляют также исследования влияния на
качество приема сигналов КАМ различных дестабилизирующих факторов,
которые условно можно разделить на внутрисистемные и внешние. К первой
группе относятся частотные и амплитудные характеристики каналов и
устройств обработки, погрешности синхронизации и т.п., ко второй –
нефлуктуационные помехи, замирания в радиоканале и др. Рассмотрим их
более подробно.
В [20] получена математическая зависимость эквивалентных
энергетических потерь, возникающих в когерентном приемнике из-за
неравномерности амплитудно-частотных характеристик (АЧХ)
используемых полосовых фильтров. Показано, что прием сигналов с КАМ
повышенной позиционности при одних и тех же неравномерностях АЧХ
сопровождается большими эквивалентными энергетическими потерями. В
[21] эта методика применена к оценке влияния неравномерности АЧХ
согласованного фильтра приемника и неравномерности сигнального
созвездия. Описаны уравнения, определяющие энергетический баланс
радиолинии и суммарные потери, возникающие в тракте. В [22] выявлено,
что помехоустойчивость приема сигналов квадратурной амплитудной
модуляции сильно зависит от фазового рассогласования, величина которого
определяется фазовыми нелинейными искажениями, возникающими в
основном в усилителе мощности передатчика, установившейся фазовой
ошибкой системы восстановления когерентной несущей и фазовыми шумами
генераторов.
Большая часть работ посвящена оценке влияния неидеальности работы
системы синхронизации на качество приема сигналов КАМ, в частности,
ошибок в определении фазы несущей [23,24] и статической и динамической
ошибок функционирования системы тактовой синхронизации [25-27]. В
 32

работе [28] проведен анализ влияния смещения несущей частоты и

временных параметров сигнала.
Влияние дефазирования опорных колебаний на помехоустойчивость
системы с квадратурной модуляцией показано в [29]. Проведенные
теоретические расчеты, построенные математические и компьютерные
модели позволили прогнозировать энергетические и вероятностные
параметры системы. В [30] оценено влияние качества функционирования
синтезаторов колебаний гетеродинов на помехоустойчивость приема
сигналов с квадратурной амплитудной модуляцией высокой позиционности,
построены графики зависимости эквивалентных энергетических потерь от
интегральной относительной нестабильности частоты и спектральной
плотности мощности фазовых шумов колебаний гетеродинов.
Похожая задача, но применительно в точности работы квадратурного
преобразователя приемника решена в [31-33]. В частности, в [31]
представлена экспериментальная оценка нескольких относительно
несложных методов совместной оценки и компенсации I/Q-разбаланса и
сдвига постоянной составляющей. В работе [32] получены результаты для
случая сдвига несущей частоты. В [33] исследуется влияние погрешностей
квадратурного преобразования на качество приема сигнала с модуляцией
КАМ-64, проводится оценка амплитудно-фазового разбаланса и сдвига
постоянной составляющей, наблюдаемых на фоне аддитивного белого шума.
Также проведено сравнение с оптимальным алгоритмом оценки и
компенсации, синтезированным методом нелинейной фильтрации. В [34]
рассмотрено влияние фазового и амплитудного разбаланса квадратур на
вероятность символьной ошибки приема сигналов КАМ.
В [35] получены результаты расчетов помехоустойчивости
посимвольного когерентного приема сигналов квадратурной амплитудной
манипуляции, подвергнувшихся нелинейным (компрессионным) искажениям
в аналоговых трактах цифровых систем передачи информации. Здесь
 33

проведено сравнение влияния компрессионных искажений на

помехоустойчивость приема сигналов 4-КАМ, 16-КАМ и 64-КАМ.
Влияние внешних дестабилизирующих факторов, а именно, внешних
аддитивных структурных и мультипликативных помех, рассмотрено в [36-
41]. В частности, в [36] рассмотрено влияние цифровых помеховых сигналов
на вероятность ошибок приема сигнала с модуляцией 16-КAM и проведено
численное моделирование с целью определения видов помех, которые
вызывают максимальную вероятность ошибок приема. Установлено, что
воздействие помеховых QPSK-сигналов вызывает больше битовых ошибок,
чем воздействие 16-КAM сигналов. Кроме того показано, что на вероятность
битовых ошибок существенно влияет величина частотного сдвига между
информационным и помеховым сигналом. Оценка параметров случайного
сигнала для дальнейшей компенсации его искажений на фоне не только
белого аддитивного шума, но и при наличии полосовой
квазидетерминированной помехи, образованной наложением нескольких
квазигармонических колебаний, разнесенных по частоте, рассматривается в
[37].
В [38,40] выполнена сравнительная оценка помехоустойчивости
приема сигнала КАМ в условиях воздействия гармонической или
фазоманипулированной помех. Установлено, что помехоустойчивость в
условиях воздействия рассмотренных помех снижается с ростом значения
числа позиций при всех возможных сочетаниях замираний сигнала и помехи,
а наибольшая помехоустойчивость его приема в условиях воздействия
гармонической или фазоманипулированной помех наблюдается при наличии
замираний только помехи. Определено, что наиболее неблагоприятной из
рассмотренных является гармоническая помеха, так как при фиксированных
значениях отношения помеха/сигнал она обеспечивает наибольшую
вероятность искажения бита. В [39] на основе математического аппарата
индикаторов событий получены формулы средней вероятности правильного
приема незамирающего M-КАМ сигнала в корреляционном приемно-
 34

решающем устройстве на фоне незамирающей импульсной гармонической

помехи. В работе [41] рассмотрены типичные помехи в каналах связи
цифровых систем передачи с сигналами OFDM/КAM, приведена
классификация помех, определены виды искажений диаграммы созвездия.
Работы [42-44] посвящены исследованиям влияния радиоканалов с
многолучевостью и замираниями. Так в [42] проведен сравнительный анализ
эффективности разнесенного приёма сигналов с квадратурной амплитудной
модуляцией и сигналов с фазовой модуляцией в каналах с дискретной
многолучевостью методом, заключающимся в сравнении изменений
сигнальных расстояний из-за такой многолучёвости. Эффективность
разнесённого приёма исследовалась методом компьютерного моделирования.
Установлено, что в таких каналах сигналы с 16-КАМ имеют значительное
преимущество перед 16-ФМ.
В [43] с использованием модифицированного обобщенного критерия
максимального правдоподобия синтезирован алгоритм приема сигналов
КАМ в каналах связи с замираниями. Для оценки параметров канала
предложено передавать специальную группу пилотных символов, состоящую
из элементов сигнальной последовательности c максимальной энергией. В
[44] рассмотрены алгоритмы оптимального и неоптимального приема КАМ
сигналов в условиях быстрых замираний в радиоканале.
Влияние межсимвольной интерференции показано в [45,46]. В
частности, в работе [45] рассматривается ситуация, когда канал не согласован
с сигналом, а в [46] получены аналитические выражения и графические
зависимости для оценки влияния межсимвольной интерференции на
помехоустойчивость приема сигналов КАМ. Показано, что с ростом
позиционности влияние межсимвольной помехи возрастает пропорционально
кратности модуляции.
 35

1.7. Выводы по разделу 1

При технической реализации методов и средств широкополосного

информационного доступа мобильных абонентов, использующих
авиационный или железнодорожный транспорт, предполагается
использовать, в том числе, и действующие сети мобильной связи технологии
LTE. В этом случае в системах радиосвязи будут использованы
многопозиционные сигналы с квадратурной амплитудной модуляцией.
Исследованием помехоустойчивости приемников сигналов КАМ на
фоне белого гауссовского шума занималось много авторов, и результаты этих
исследований достаточно широко известны. Большой интерес представляют
также исследования влияния на качество приема сигналов КАМ различных
дестабилизирующих факторов – внутрисистемных и внешних, зависящих от
состояния радиоканала. Важность решения этой задачи подтверждается
наличием соответствующих научных публикаций, появившихся в последние
5-10 лет. В них используются различные методики оценки качества связи при
таких условиях, и, к сожалению, не всегда удается адекватно сопоставить
результаты, полученные разными авторами. Представляется целесообразным
получить и сравнить в рамках единой методики параметры
помехоустойчивости приема сигналов КАМ при воздействии
нефлуктуационных помех, а также внутрисистемных погрешностей. Кроме
того, пока мало исследований, посвященных методам борьбы с
нефлуктуационными помехами. В связи со сказанным можно заключить, что
изучение данной темы актуально, и этим вопросам будут посвящены
последующие разделы диссертационной работы.
Результаты данного раздела опубликованы в [5,47-51].
 36

2. АНАЛИЗ ПОМЕХОУСТОЙЧИВОСТИ КОГЕРЕНТНОГО ПРИЕМА

СИГНАЛОВ КАМ ПРИ НАЛИЧИИ НЕФЛУКТУАЦИОННЫХ ПОМЕХ

В данном разделе выполнен анализ помехоустойчивости когерентного

приемника сигналов КАМ (рис. 1.9) при наличии шумовой n(t ) и

нефлуктуационных помех sп (t ) разного вида:

- гармоническая помеха (1.5) с частотной расстройкой ωп

относительно несущей сигнала КАМ;
- сканирующая помеха (1.8) с девиацией частоты ωд ;

- фазоманипулированная помеха (1.6) с частотной расстройкой ωп

относительно несущей сигнала КАМ и длительностью тактового интервала
помехи Tп  Ts ;
- ретранслированная помеха (1.7) с временной задержкой τ  Ts .

2.1. Методика определения вероятности ошибки когерентного

приемника сигналов КАМ при наличии нефлуктуационных помех

Когерентный приемник сигнала (1.2) (рис. 1.9) в присутствии белого

гауссовского шума n(t) содержит М корреляционных каналов и вычисляет М
интегралов свертки
2 Aср Ts
Ji 
N0  x(t )s
0
оп i (t )dt , i  0,1,...M  1 (2.1)

принимаемого колебания

x(t )  si (t )  sп (t )  n(t ) (2.2)

с опорными сигналами
sопi (t )  ( Ii cos ω0t  Qi sin ω0t ) .

Решающее устройство в результате сравнения вычисленных интегралов

(с учетом порога) принимает решение по максимуму правдоподобия в пользу
 37

того или иного канального символа. При этом полагаем, что синхронизация
приемника идеальна.
В этом случае вероятность ошибочного приема m-го канального
символа можно вычислить по формуле
M 1
Pes  1   pi ( J m  J i   mi ) m ,
i 0
(2.3)
m i

где pi ( J m  J i   mi ) m – вероятность того, что выходное значение m-го

коррелятора приемника больше выходного значения любого другого i-го
коррелятора при условии, что передавался m-й символ;
Esm  Esi Es ср
 mi   [( I m 2  Qm 2 )  ( I i 2  Qi 2 )] (2.4)
N0 N0
– порог принятия решения, определяемый полуразностью энергий
сравниваемых канальных символов.
Вероятность битовой ошибки Реb при достаточно большом отношении
сигнал/шум с учетом кодирования Грея можно найти по формуле [1]:

Peb  Pes / log2 M . (2.5)

Предположим сначала, что случайная начальная фаза помехи является

фиксированной величиной. В этом случае распределения случайных
процессов Ji на выходах всех корреляторов демодулятора можно считать
нормальными [52,53] и можно рассчитать условные по этому параметру
статистические характеристики: средние значения mi , их линейные
комбинации mmi и дисперсии Dmi.
Для определения этих характеристик следует подставить выражение
для сигнала и выражение для нефлуктуационной помехи в (2.2) и далее – в
(2.1). При этом учтем, что
N0
 n(t )  0;  n(t1 )n(t2 )  δ(t2  t1 ); ω0Ts  1 .
2
Например, статистические характеристики с индексами m и i
определяются следующим образом.
 38

Средние значения:
2 Aср Ts 2 Aср Ts
mm  J m 
N0   x(t )  s
0
оп m (t )dt 
N0  s
0
m (t )  sп (t )  n(t )  sоп m (t )dt 

Ts Ts
2 Aср 2 Aср 2 Aср Ts

N0  s
0
m (t )  sп (t )  sоп m (t )dt 
N0   n(t )  s
0
оп m (t )dt 
N0  (s
0
m (t )  sп (t ))sоп m (t )dt;

2 Aср Ts 2 Aср Ts
mi  J i 
N0   x(t )  s
0
оп i (t )dt 
N0  s
0
m (t )  sп (t )  n(t )  sопi (t )dt 

2 Aср Ts 2 Aср Ts 2 Aср Ts


N0  s
0
m (t )  sп (t )  sопi (t )dt 
N0   n(t )  s
0
оп i (t )dt 
N0  (s
0
m (t )  sп (t )) sопi (t )dt;

mmi  mm  mi   mi .
Дисперсии:
Dm  J m 2    J m  2  J m 2  mm 2 
2 Ts Ts
4 Aср

N0 2    x(t )x(t )  s
0 0
1 2 оп m (t1 ) sоп m (t2 )dt1dt2  mm 2 

2 Ts Ts
4 Aср

N02   [s
0 0
m (t1 )  sп (t1 )][ sm (t2 )  sп (t2 )]sоп m (t1 ) sоп m (t2 )dt1dt2 

2 Ts Ts
4 Aср

N0 2    n(t )n(t )  s
0 0
1 2 оп m (t1 ) sоп m (t2 )dt1dt2  m0 2 .

Di  J i 2    J i  2  J i 2  mi 2 
2 Ts Ts
4A
 20
N0    x(t )x(t )  s
0 0
1 2 оп i (t1 ) sопi (t2 )dt1dt2  mi 2 

2 Ts Ts
4 A0

N02   [s
0 0
m (t1 )  sп (t1 )][ sm (t2 )  sп (t2 )]sоп i (t1 ) sоп i (t2 )dt1dt2 

2 Ts Ts
4 A0

N02    n(t )n(t )  s
0 0
1 2 оп i (t1 ) sопi (t2 )dt1dt2  mi 2 .

Смешанные моменты:
 39

2 Ts Ts
4 Aср
M mi 
N02    x(t )x(t )  s
0 0
1 2 оп m (t1 ) sопi (t2 )dt1dt2  mm mi 

2 Ts Ts
4 Aср

N02   [s
0 0
m (t1 )  sп (t1 )][ sm (t2 )  sп (t2 )]sоп m (t1 ) sоп i (t2 )dt1dt2 

2 Ts Ts
4 Aср

N02    n(t )n(t )  s
0 0
1 2 оп m (t1 ) sопi (t2 )dt1dt2  mm mi .

Dmi  Dm  Di  2M mi .

Расчеты должны быть проведены по всем возможным комбинациям m

и i. В результате можно определить вероятности, условные по параметрам
помех, входящие в (2.3):

m 1

2
pi ( J m  J i   mi ) m  1  Q( mi ), Q( x)  e t /2 dt. (2.6)
Dmi 2 x

и далее – условные вероятности символьной (2.3) и битовой (2.5) ошибки.

Безусловную вероятность битовой ошибки получаем усреднением
полученных результатов по параметрам помех, например, по случайной
величине φп :
π
1
2π π
Pеb  Peb (φп )dφп .

В данном случае в связи со сложностью аналитического

интегрирования было использовано численное усреднение.
При вычислениях необходимо провести перебор всех возможных
фазовых и амплитудных состояний сигнала КАМ и учесть нормировку по
средней энергии с целью адекватного сравнения полученных результатов с
результатами для других сигналов. Последнее можно сделать,
пронормировав величины коэффициентов I i и Qi . При этом логично
положить, чтобы средние энергии сравниваемых сигналов были одинаковы.
В случае 16-КАМ (рис. 1.7) в соответствии с положением сигнальных
точек имеем три уровня энергии посылок:
 40

( x2  x2 ) ( x2  9x2 ) (9 x 2  9 x 2 )
E1  Ts  Emin , E2  Ts , E3  T  Emax .
2 2 2

Отсюда получим:

x 2T AсрT x  1 A .
2
4 E1  8E2  4 E3
Esср   10  , ср
16 2 2 10

Следовательно, для 16-КАМ при вычислениях в зависимости от выбора

сигнальной точки необходимо использовать корректирующие коэффициенты
1 3
и .
10 10

Подобные расчеты для 32-КАМ дают корректирующие коэффициенты

1 3 5
, и ,
20 20 20

1 3 5 7
а для 64-КАМ: , , и .
42 42 42 42

2.2. Анализ вероятности ошибки когерентного приемника сигналов

КАМ при наличии гармонической помехи
Считаем, что на входе приемника кроме полезного сигнала si (t )
присутствует белый гауссовский шум n(t) и гармоническая помеха

sп (t )  μA0 cos[(ω0  ωп )t  φп ] (2.7)

со случайной начальной фазой φп, частотной расстройкой ∆ωп и

относительной интенсивностью μ.
Вычислим средние значения и дисперсии случайных процессов в схеме
когерентного приемника по методике п. 2.1. Для примера покажем
вычисления для сигналов с индексами m и i.
 41

2 Aср Ts
mm  J m 
N0   x(t )  ( I
0
m cos ω0t  Qm sin ω0t )dt 

2 Ts
2 Aср

N0  (I
0
m cos ω0t  Qm sin ω0t  μ cos[(ω0  ωп )t  φ п ])( I m cos ω0t  Qm sin ω0t )dt 

2 Ts 2 Ts 2 Ts
2 Aср 2 Aср 4 Aср
  I m cos ω0tdt  Q sin ω0tdt  I Qm cos ω0t sin ω0tdt 
2 2 2 2
m m
N0 0
N0 0
N0 0
2 T 2 T
2 Aср Im s 2 Aср Qm s
μ  cos(ωпt  φ п )dt 
2 0
 μ sin(ωпt  φ п )dt.
N0 2 0 N0

Проинтегрировав в указанных пределах, можно получить:

ω t
2 Es ср 2 E sin п
mm  ( I m 2  Qm 2 )  s ср μ 2 ( I cos  Q sin  ),
N0 N0 ωпt m m

ωпt
где параметр    φп , также как и  п , является случайной величиной,
2
равномерно распределенной на полуинтервале (-, +].

2 Aср Ts
mi  J i 
N0   x(t )  ( I cos ω t  Q sin ω t )dt 
0
i 0 i 0

2 Ts
2 Aср

N0  (I
0
m cos ω0t  Qm sin ω0t  μ cos[(ω0  ωп )t  φ п ])( I i cos ω0t  Qi sin ω0t )dt 

2 Ts 2 Ts 2 Ts
2 Aср 2 Aср 2 Aср
  I m I i cos ω0tdt   I mQi cos ω0t sin ω0tdt   I i Qm cos ω0t sin ω0tdt 
2

N0 0
N0 0
N0 0
2 Ts 2 Ts 2 T
2 Aср 2 Aср Ii 2 Aср Qi s
  Q Q sin ω0tdt  μ  cos(ωпt  φ п )dt  μ  sin(ωпt  φ п )dt.
2
m i
N0 0
N0 20 N0 2 0

После интегрирования можно получить:

ωпt
2E 2E sin
mi  s ср ( I m I i  QmQi )  s ср μ 2 ( I cos  Q sin  ).
N0 N0 ωпt i i

2
 42

Дисперсии:

Dm  J m 2    J m  2  J m 2  m0 2 
2 Ts Ts
4 Aср

N0 2    x(t )x(t )  ( I
0 0
1 2 m cos ω0t1  Qm sin ω0t1 )( I m cos ω0t2  Qm sin ω0t2 )dt1dt2  mm 2 

2
4 Aср N 0 Ts Ts

2 0 
 mm  2
2
( I m cos ω0t1  Qm sin ω0t1 )( I m cos ω0t2  Qm sin ω0t2 )dt1dt2  mm 2 
N0 0

2 Es ср
 ( I m 2  Qm 2 ).
N0

Di  J i 2    J i  2  J i 2  mi 2 
2 Ts Ts
4 Aср

N0 2    x(t )x(t )  ( I cos ω t
0 0
1 2 i 0 1  Qi sin ω0t1 )( I i cos ω0t2  Qi sin ω0t2 )dt1dt2  mi 2 

2
4 Aср N 0 Ts Ts

2 0 
 mi  2
2
( I i cos ω0t1  Q0 sin ωi t1 )( I i cos ω0t2  Qi sin ω0t2 )dt1dt2  mi 2 
N0 0

2 Es ср
 ( I i 2  Qi 2 ).
N0

Смешанный момент:

2 Ts Ts
4 Aср
M mi 
N02    x(t )x(t )  ( I
0 0
1 2 m cos ω0t1  Qm sin ω0t1 )( I i cos ω0t2  Qi sin ω0t2 )dt1dt2 

2
4 Aср N 0 Ts Ts

2 0 
mm mi  mm mi  2 ( I m cos ω0t1  Qm sin ωi t1 )( I i cos ω0t2  Qi sin ω0t2 )dt1dt2 
N0 0

2 Es ср
mm mi  ( I m I i  QmQi ).
N0
Далее найдем статистические параметры суммарных процессов
ymi  J m  J i   mi на входе устройства сравнения корреляционных интегралов.

Для примера, среднее значение:

 43

ω t
2 Es ср 2 E sin п
mmi  ( I m 2  Qm 2  I m I i  QmQi ) 
s ср
μ 2 [( I  I ) cos  (Q  Q )sin  ] 
N0 N0 ωпt m i m i

2
E 2E
 s ср [( I m 2  Qm 2 )  ( I i 2  Qi 2 )]  s ср {0,5[( I m  I i ) 2  (Qm  Qi ) 2 ]  (2.8)
N0 N0
ωпt
sin
μ 2 [( I  I ) cos  (Q  Q )sin  ]},
ωпt m i m i

2
дисперсия:

2 Es ср
Dmi  [( I m  I i )2  (Qm  Qi ) 2 ] . (2.9)
N0

Это позволило определить вначале условные параметру η, а затем –

после численного усреднения по диапазону его значений – и безусловные
вероятности ошибки приема канального символа и битовой ошибки (2.5).
При вычислениях проведен перебор всех возможных фазовых и
амплитудных состояний сигнала КАМ и учтена нормировка по средней
энергии (см. п. 2.1).
На рис. 2.1 – 2.4 показаны зависимости вероятности ошибки Peb от
отношения сигнал/шум Ebср / N0  Esср / ( N0 log2 M ) при нулевой расстройке

гармонической помехи пТs, то есть при прицельной помехе. Кривые на

каждом рисунке различаются величиной интенсивности помехи μ. Здесь и
далее на графиках величины Ebср / N0 приведены в дБ.
Анализ графиков позволяет сделать вывод, что с увеличением степени
позиционности М полезного сигнала влияние гармонической помехи
значительно усиливается. Так, если слабая помеха (μ=0,1) для М=4
практически не увеличивает вероятность ошибки Peb, то для М=16, 32 и 64
снижение помехоустойчивости становится очень заметным – при Ebср / N0
>15 дБ вероятность ошибки увеличивается на несколько порядков.
 44

Рис. 2.1. Зависимость вероятности битовой ошибки от отношения

сигнал/шум (дБ) для 4-КАМ в присутствии гармонической помехи.

Рис. 2.2. Зависимость вероятности битовой ошибки от отношения

сигнал/шум (дБ) для 16-КАМ в присутствии гармонической помехи.
 45

Рис. 2.3. Зависимость вероятности битовой ошибки от отношения

сигнал/шум (дБ) для 32-КАМ в присутствии гармонической помехи.

Рис. 2.4. Зависимость вероятности битовой ошибки от отношения

сигнал/шум (дБ) для 64-КАМ в присутствии гармонической помехи.
 46

Прицельные помехи с интенсивностью μ≥0,3 опасны для всех видов

сигналов КАМ с М>4 и практически разрушают прием.
Этот вывод подтверждают и графики, приведенные на рис. 2.5, где для
сравнения даны кривые помехоустойчивости для всех исследованных
сигналов в зависимости от величины μ.

Рис. 2.5. Зависимость вероятности битовой ошибки от интенсивности

гармонической помехи при Ebср/N0=13 дБ.

Был исследован характер изменения помехоустойчивости приема

сигналов КАМ в зависимости от величины расстройки гармонической
помехи пТs. Соответствующие кривые показаны на рис. 2.6.
 47

Рис. 2.6. Зависимость вероятности битовой ошибки от расстройки

гармонической помехи при Ebср/N0=13 дБ, μ=0,5.

Наглядно видно, что во всех случаях помехоустойчивость приема

растет при увеличении расстройки помехи и в пределе стремится к
помехоустойчивости приема без помехи. Точки, в которых вероятность
ошибки достигает своих минимумов, соответствуют нулям огибающей
энергетического спектра полезного сигнала. При этих расстройках
поражающее действие помехи никак не проявляется. Можно также сделать
вывод, что при больших расстройках помехи (п>5/Тs) ее влиянием можно
пренебречь.
 48

2.3. Анализ вероятности ошибки когерентного приемника сигналов

КАМ при наличии сканирующей помехи

Предположим, что наряду с полезным сигналом КАМ на входе

приемника действует мешающий сигнал – сканирующая помеха – с
относительной интенсивностью μ и частотой ωп , изменяющейся вокруг
несущей частоты ω0 сигнала КАМ по пилообразному закону (рис. 2.7):
2ωдt
sп (t )  μA0 cos(ωп (t )t  φп ), ωп (t )  ω0  ωд  , (2.10)
Tс

где ωд – девиация помехи; Tс – период сканирования.

ωп

ω0  ωд

ω0
0 Tc t

ω0  ωд

Рис. 2.7. Закон изменения частоты сканирующей помехи.

Для случая Tc  Ts такую помеху на каждом тактовом интервале

сигнала можно считать гармонической с некоторой расстройкой несущей
частоты ωп относительно величины ω0 . Тогда для расчета вероятности
ошибки Peb можно воспользоваться результатами, полученными в п. 2.2.
Для этого необходимо провести усреднение полученных результатов
Peb (φп , ωп ) по двум параметрам: случайная величина φ п , равномерно
распределенная на интервале (–, ], и величина ωп , равномерно

распределенная на интервале ( ωд , ωд ]:

π ωд
1 1
Pеb   [
2π  π 2ωд 
ωд
Pe (φп , ωп )d ωп ]dφп .
 49

Как следует из п. 2.2, наиболее опасным является случай, когда

диапазон сканирования помехи по частоте совпадает с зоной главного
лепестка спектра сигнала.
На рис. 2.8 показаны зависимости вероятности ошибки Peb от
отношения сигнал/шум Ebср N0 при ωдTs  6 для М =4, 16, 32 и 64.

а)

б)
 50

в)

г)
Рис. 2.8. Зависимости вероятности битовой ошибки от отношения
сигнал/шум при девиации частоты помехи ωдTs  6 :
а – М=4, б – М=16, в – М=32, г – М=64.
 51

Рис. 2.9 иллюстрирует для разных М зависимости вероятности ошибки

Peb от интенсивности сканирующей помехи µ при Ebср N0 =13 дБ и

ωдTs  6.

Рис. 2.9. Зависимость вероятности ошибки Peb от интенсивности

сканирующей помехи μ .

Из рисунков видно, что с увеличением степени позиционности М

влияние сканирующей помехи значительно усиливается. При этом, если при
слабой помехе (μ=0,1) для М=4 и 16 вероятность ошибки практически не
увеличивается и энергетические потери для Peb=10-3 составляют не более 1
дБ, то для М=32 или 64 при такой же помехе снижение помехоустойчивости
становится заметным, а энергетические потери уже достигают 2-5 дБ.
Помехи большой интенсивности (μ≥0,5) опасны для всех видов сигналов
КАМ и полностью разрушают прием.
 52

При увеличении девиации частоты помехи, "накрывающей" зону

главного лепестка спектра сигнала, ее влияние на помехоустойчивость
демодулятора КАМ снижается, оставаясь при этом весьма существенным
(рис. 2.10). Так при увеличении девиации на порядок вероятность ошибки
снижается в несколько раз. Это связано с тем, что основной отрицательный
вклад помехи приходится на ее частотные составляющие, близкие к частоте
несущего колебания полезного сигнала, а при расширении зоны
сканирования их доля пропорционально снижается. При стремлении
девиации частоты помехи к нулю, вероятность ошибки стремится к
величине, полученной для случая прицельной гармонической помехи.

а)
53

б)

в)
 54

г)
Рис. 2.10. Зависимости вероятности битовой ошибки от девиации частоты
сканирующей помехи при Ebср N0 =13 дБ:

а – М=4, б – М=16, в – М=32, г – М=64.

2.4. Анализ вероятности ошибки когерентного приемника сигналов

КАМ при наличии фазоманипулированной помехи

Считаем, что на входе приемника кроме полезного сигнала si (t )

присутствует белый гауссовский шум n(t) и фазоманипулированная помеха с
относительной интенсивностью µ.
Посылку фазоманипулированной помехи длительностью Tп можно
записать в виде:

sп (t )  μAсрa j cos[(ω0  ωп )t  φп ], t  (( j  1)Tп , jTп ], j  1,..., N , (2.11)

где a j  1 – случайный дискретный символ помехи; Tп  Ts N ; N –

относительная скорость манипуляции помехи; ωп – расстройка помехи

 55

относительно центральной частоты спектра сигнала КАМ; φ п – случайная

начальная фаза помехи, равномерно распределенная на полуинтервале (-π, + π].
Временные соотношения между полезным сигналом и
фазоманипулированной помехой показаны на рис. 2.11. В общем случае на
интервале Ts укладывается N тактовых интервалов помехи.

Рис. 2.11. Временные соотношения между полезным сигналом и

фазоманипулированной помехой.

Вычислим средние значения и дисперсии случайных процессов в схеме

когерентного приемника по методике п. 2.1. Для примера покажем
вычисления для сигналов с индексами m и i.
Средние значения
2 Aср Ts
mm  J m 
N0   x(t )  ( I
0
m cosω0t  Qm sin ω0t )dt 

2 Ts
2 Aср

N0  (( I
0
m cosω0t  Qm sin ω0t )  μa j cos[(ω0  ωп )t  φ п ]) 

2 Ts
2 Aср
( I m cosω0t  Qm sin ω0t )dt 
N0  (I
0
m cosω0t  Qm sin ω0t ) 2 dt 

2 Ts /Tп 1 ( j 1)Tп
2 Aср

N0
μ  
j 0
a j cos[(ω0  ωп )t  φп ]( I m cosω0t  Qm sin ω0t )dt 
jTп
 56

2 2
2 Aср I m 2Ts 2 Aср Qm 2Ts
    
N0 2 N0 2
2 Ts /Tп 1 ( j 1)Tп
2 Aср a j I m cos(ω0t  [(ω0  ωп )t  φ п ])  

N0
μ  
j 0
[
2   cos([(ω0  ωп )t  φ п ]  ω0t ) 

jTп

a j Qm
 sin([(ω0  ωп )t  φп ]  ω0t )  sin([(ω0  ωп )t  φп ]  ω0t )  dt 
2
2
Aср Ts
  ( I m 2  Qm 2 ) 
N0
2 Ts /Tп 1 ( j 1)Tп
2 Aср  a j Im aQ 

N0
μ   
2
 cos(ωпt  φп )   j m sin(ωпt  φ п )   dt  ...
2
j 0 jTп  
2 Esср 2 Esср μ sin x N 1
N0 N x 
 ( I  Qm ) 
2 2
a j (I m cos( x(2 j  1)  φ п )  Qm sin( x(2 j  1)  φ п )),
N0 m j 0

ωпTп ωпTs
x  .
2 2N

2 Aср Ts
mi  J i 
N0   x(t )  ( I cosω t  Q sin ω t )dt  ...
0
i 0 i 0

2 Esср 2 Esср μ sin x N 1

N0 N x 
 ( I m I i  QmQi )  a j (I i cos( x(2 j  1)  φп )  Qi sin( x(2 j  1)  φп )).
N0 j 0

С учетом величины порога  mi (2.4):

2 Esср
mmi  {0,5(( I m  I i ) 2  (Qm  Qi ) 2 ) 
N0
(2.12)
N 1 N 1
μ sin x
 [( I m  I i ) a j cos( x(2 j  1)  φп )  (Qm  Qi ) a j sin( x(2 j  1)  φ п )]}.
N x j 0 j 0

Как видно, среднее значение, среди прочего, зависит от случайных

величин: начальной фазы и комбинации информационных символов помехи.
Проведенные расчеты, аналогичные п. 2.1 и 2.2, показали, что величина
дисперсии равна:
2 Esср
Dmi  [( I m  I i ) 2  (Qm  Qi ) 2 ]. (2.13)
N0
 57

Полученные параметры использовались для получения условных

вероятностей (2.6). Для получения безусловной вероятности битовой ошибки
было проведено численное усреднение полученных результатов по


случайным величинам φп и a  a1 ,..., a j ,..., aN : 
π
1 1
Pеb  N
2
  Peb (φп , a)dφп .
a 2π  π

При усреднении перебирались все возможные комбинации канальных

символов и проведена нормировка энергий сигналов.
На рис. 2.12 показаны зависимости вероятности ошибки Peb от
отношения сигнал/шум Ebср / N0 при N=2 и расстройке
фазоманипулированной помехи ωпTs =0. Видно, что с увеличением
позиционности сигнала помехоустойчивость приема довольно резко падает,
и в результате даже помехи малого уровня μ=0,1 вызывают существенные
энергетические потери – единицы дБ при М>16. При μ>0,5 прием сигнала
становится невозможным.
Заметим, что, как и ранее, из-за ограничений применимости формулы
(2.5) результаты при малых отношениях сигнал/шум Ebср / N0 <3дБ и большой
интенсивности помех μ>0,5 представляют собой нижнюю границу
вероятности ошибки и имеют определенную несущественную погрешность.
Выражение (2.12) и графики получены для случая, когда границы
тактовых интервалов сигнала и фазоманипулированной помехи совпадают
(рис. 2.11). Если относительное временное смещение фронтов помехи  п / Т п
не равно нулю, то в начальные и конечные составляющие сумм (2.12)
добавятся корректирующие коэффициенты, пропорциональные величине
этого смещения. Так, для упрощенного случая, когда ωпTs =0,
2 Esср
mmi  {0,5(( I m  I i ) 2  (Qm  Qi ) 2 ) 
N0
(2.14)
μ  N 2

 [( I m  I i )cosφп  (Qm  Qi )sin φ п ][ п a1   a j  (1  п )aN 1 ]}.
N Тп j 0 Тп
58

а)

б)
 59

в)

г)
Рис. 2.12. Зависимость вероятности битовой ошибки от отношения
сигнал/шум (дБ) для КАМ в присутствии фазоманипулированной помехи:
а – 4-КАМ, б – 16-КАМ, в – 32-КАМ, г – 64-КАМ.
 60

На рис. 2.13 в качестве примера для случаев М=16, 64 и N=2, 4 при

Ebср / N0 =13дБ и μ=0,3 показано, что величина вероятности битовой ошибки в

зависимости от смещения  п / Т п изменяется несущественно, особенно при

больших M.

Рис. 2.13. Зависимость вероятности битовой ошибки от временного

смещения фронтов помехи  п / Т п .

Зависимость помехоустойчивости приема сигналов КАМ от

интенсивности помехи иллюстрируют графики, приведенные на рис. 2.14.
Кривые построены для двух случаев: N=2 и N=5. При изменении μ от 0 до 1
вероятность ошибки меняется на несколько порядков. Это подтверждает
вывод о том, что при больших уровнях помех необходимы специальные
средства борьбы с ними. Ситуация осложняется тем, что рассматриваемая
фазоманипулированная помеха является достаточно широкополосной и
степень этой широкополосности зависит от ее относительной скорости
передачи N.
 61

Рис. 2.14. Зависимость вероятности битовой ошибки для КАМ от

интенсивности фазоманипулированной помехи при Ebср/N0=13 дБ.

Из рис. 2.15 видно, как изменяется влияние такой помехи при разных
величинах N. Так, при N=5 действие помехи заметно даже при большой
расстройке фазоманипулированной помехи ωпTs . Волнообразный характер
зависимостей объясняется результатом интерференционного взаимодействия
сигнала и помехи. Так, снижение вероятности ошибки при N=2 на интервалах
вблизи точек ωпTs = 4π и 8π связано с тем, что именно на эти точки
приходятся нули огибающей спектра помехи.
 62

Рис. 2.15. Зависимость вероятности битовой ошибки для КАМ от расстройки

фазоманипулированной помехи при Ebср/N0=13 дБ, μ=0,5.

Дополнительную информацию дает анализ графиков, приведенных на

рис. 2.16, где показаны зависимости вероятности битовой ошибки для КАМ
от относительной скорости передачи фазоманипулированной помехи при
расстройке ωпTs =0. Во всех случаях поражающее действие такой помехи
снижается при увеличении N, что объясняется тем, что помеха становится по
своему характеру шумоподобной. При N=1 ситуация эквивалентна случаю
гармонической помехи (см. п. 2.2).
63

а)

б)
 64

в)

г)
Рис. 2.16. Зависимость вероятности битовой ошибки для КАМ от
относительной скорости передачи фазоманипулированной помехи при
Ebср/N0=13 дБ: а – 4-КАМ, б – 16-КАМ, в – 32-КАМ, г – 64-КАМ.
 65

2.5. Анализ вероятности ошибки когерентного приемника сигналов КАМ

при наличии ретранслированной помехи

Считаем, что на входе приемника кроме полезного сигнала si (t )

присутствует белый гауссовский шум n(t) и ретранслированная помеха с
относительной интенсивностью µ, например, задержанный на время   Ts
сигнал, отраженный от земной поверхности или местных предметов:
sп (t )  μsi (t   ,φп ) . (2.15)
Представление о временных соотношениях между полезным сигналом
КАМ и ретранслированной помехой дает рис. 2.17.

Рис. 2.17. Временные соотношения между полезным сигналом и

ретранслированной помехой.

Для анализа помехоустойчивости приемника КАМ сигнала аналогично

пп. 2.2-2.4 рассчитаем статистические характеристики распределений,
условные по случайной начальной фазе φ п :
 66

2 Esср 
mmi  (0.5 ( I m  I i ) 2  (Qm  Qi ) 2   μ( I m I j  Q m Q j  I i I j  Q i Q j ) cos φ 
N0 Ts
 
 μ( I mQ j  I j Qm  I i Q j  I j Qi ) sin φ  (1  )μ( I m 2  Qm 2  I m I i  QmQi ) cos φ  (2.16)
Ts Ts

(1  )μ( I mQi  I i Qm ) sin φ),
Ts

2 Es ср
Dmi  [( I m  I i )2  (Qm  Qi )2 ] , (2.17)
N0

где φ  ω0  φп .

Полученные характеристики позволили определить условные

вероятности, входящие в (2.3), и условную вероятность битовой ошибки. Для
получения безусловной вероятности битовой ошибки численным методом
проведено усреднение результатов по случайной величине φ :

1
Pеb 
2 P

eb (φ)dφ

и выполнена нормировка энергий сигналов с помощью коэффициентов I i и Qi

(см. п. 2.1).
Естественно, что из-за суперпозиции прямой и задержанной посылок с
разными фазами условия приема полезного сигнала могут как существенно
ухудшаться, так и несколько улучшаться. В среднем это все же приводит к
увеличению вероятности ошибки.
На рис. 2.18 приведены зависимости усредненной по начальной фазе
помехи вероятности битовой ошибки Реb от относительной задержки помехи
τ/Ts для интенсивности помехи µ=0,5 при отношении сигнал/шум Ebср / N0 =13
дБ. Увеличение вероятности ошибки в областях τ/Ts=0 и τ/Ts=1 объясняется
тем, что при усреднении по начальной фазе помехи доля «опасных»
ситуаций, когда сигнал и помеха противофазны, и степень снижения уровня
сигнала за счет суперпозиции колебаний, здесь выше, чем для области
τ/Ts=0,5.
 67

Рис. 2.18. Зависимость усредненной вероятности ошибки Реb от задержки

ретранслированной помехи τ/Ts:
1 - 4-КАМ, 2 - 16-КАМ, 3 - 32-КАМ, 4 - 64-КАМ.

На рис. 2.19 представлены зависимости вероятности битовой ошибки

для всех исследованных сигналов от величины μ при отношении Ebср / N0 =13 дБ
и относительной задержке помехи τ/Ts =0,5. Из их сравнения следует, что
ретранслированные помехи большой интенсивности (μ≥0,3) опасны для всех
видов сигналов КАМ и практически разрушают прием.
Наиболее показательной является зависимость вероятности ошибки Реb
от отношения сигнал/шум на входе приемника Ebср / N0 . На рис. 2.20 приведены
рассчитанные зависимости для М = 4, 16, 32 и 64 при разных интенсивностях
помехи и τ/Ts =0,5. Здесь же для сравнения линией 1 приведен график
помехоустойчивости рассматриваемого приемника по отношению только к
белому гауссовскому шуму.
 68

Рис. 2.19. Зависимость вероятности ошибки Реb от интенсивности

ретранслированной помехи µ:
1 - 4-КАМ, 2 - 16-КАМ, 3 - 32-КАМ, 4 - 64-КАМ.

а)
69

б)

в)
 70

г)

Рис. 2.20. Зависимость вероятности ошибки Реb от отношения сигнал/шум

Ebср / N 0 (дБ) для сигналов 4-КАМ (а), 16-КАМ (б), 32-КАМ (в), 64-КАМ (г):

1 – µ =0, 2 – µ =0,1, 3 – µ =0,3, 4 – µ =0,5, 5 – µ =0,9.

С увеличением степени позиционности М полезного сигнала влияние

ретранслированной помехи значительно усиливается. Анализ показывает, что
присутствие на входе приемника ретранслированной помехи с
интенсивностью μ=0,1 при Реb=10-3 для 4-КАМ эквивалентно ухудшению
отношения сигнал/шум всего на 0,2 дБ, для 16-КАМ это ухудшение
составляет уже 0,8 дБ, для 32-КАМ – 2 дБ, а для 64-КАМ – более 3 дБ. При
увеличении интенсивности помехи энергетические проигрыши существенно
возрастают.
 71

2.6. Сравнительная оценка влияния различных нефлуктуационных

помех
Для сравнения результатов расчетов вероятностей ошибки при разных
видах нефлуктуационных помех на рис. 2.21 показаны сводные зависимости
вероятности ошибки Peb от отношения сигнал/шум Ebср N0 при одинаковой

интенсивности мешающих воздействий μ= 0,1 для М=4 и М=32. Рассмотрен

случай, когда все помехи являются прицельными, т.е. центральные частоты
их спектров совпадают с несущей частотой ω0 сигнала КАМ. Для
фазоманипулированной помехи принято N=2, для ретранслированной – τ/Ts
=0,5, для сканирующей – ωдTs  6 .

а)
 72

б)
Рис. 2.21. Зависимость вероятности битовой ошибки Peb от отношения
сигнал/шум Ebср N0 для М=4 (а) и М=32 (б).

Наибольшее снижение помехоустойчивости наблюдается в случае

гармонической помехи: для М=32 энергетические потери по сравнению со
случаем отсутствия помехи составляют около 2 дБ при вероятности ошибки
10-2. В остальных случаях проигрыши меньше и составляют следующие
значения: фазоманипулированная – 1,2 дБ; ретранслированная – 0,6 дБ,
сканирующая – 0,5 дБ. Для М=4 при вероятности ошибки 10-6 проигрыши
меньшие и составляют, соответственно, 1,2 дБ; 0,65 дБ; 0,5 дБ; 0,4 дБ.
При возрастании интенсивности помех ситуация сильно ухудшается.
Это подтверждают графики, приведенные на рис. 2.22, где показаны кривые
помехоустойчивости для всех вышерассмотренных нефлуктуационных помех
в зависимости от интенсивности помех μ при М=4 и М=32. Видно, что
возрастание величины μ от 0 до 1 вызывает ухудшение вероятности ошибки
на несколько порядков.
 73

Рис. 2.22. Зависимость вероятности битовой ошибки Peb от интенсивности

помех μ для М=4 и М=32.

Показанное соотношение влияния помех соблюдается не всегда и

может меняться, так как зависит от множества параметров помех. Заметим
также, что влияние помех значительно снижается при их большой частотной
расстройке относительно центральной частоты спектра полезного сигнала.
По графикам раздела 2 рассчитаны величины энергетических
проигрышей при приеме сигналов КАМ на фоне исследованных
нефлуктуационных помех, результаты анализа приведены в табл. 2.1 и 2.2.

Таблица 2.1. Величины энергетического проигрыша при приеме сигналов

КАМ на фоне нефлуктуационных помех малой интенсивности.
Тип помехи Энергетический проигрыш (дБ) при вероятности ошибки 10-3
(интенсивность µ = 0,1) М=4 М = 16 М = 32 М = 64

Гармоническая 0,3 1,5 3,1 7

Ретранслированная 0,22 1 2 4
Фазоманипулированная 0,22 1,3 2,5 6
Сканирующая 0,15 0,8 1,8 3,8
 74

Таблица 2.2. Величины энергетического проигрыша при приеме сигналов

КАМ на фоне нефлуктуационных помех большой интенсивности.
Тип помехи Энергетический проигрыш (дБ) при заданной вероятности
(интенсивность µ = 0,5) ошибки
М=4 М = 16 М = 32 М = 64
10-3 10-3 10-3 10-3
Гармоническая >10 Проигрыш Проигрыш Проигрыш
Ретранслированная 7 критически критически критически
Фазоманипулированная 7,4 велик велик велик
Сканирующая 6
В заключении проведем сравнение помехоустойчивости когерентных
приемников сигналов КАМ и сигналов с многопозиционной фазовой
манипуляцией (М-ФМ) [54-57] одинаковой позиционности и с одинаковой
средней энергией в присутствии рассмотренных нефлуктуационных помех.
На рис. 2.23 показаны сводные графики для вероятности ошибки Реb при
μ=0.1 и следующих параметрах помех: фазоманипулированная – при
относительной скорости манипуляции помехи N=3, ретранслированная – при
относительной задержке /Тs=0,5, сканирующая – при дTs =6. Рассмотрены
случаи, когда центральные частоты спектров помех совпадают с несущей
частотой ω0 полезного сигнала.

а)
 75

б)

в)
Рис. 2.23. Зависимости вероятности битовой ошибки от отношения
сигнал/шум для КАМ (штриховые линии) и М-ФМ (сплошные линии): 1 –
без помехи; при нефлуктуационных помехах: 2 –сканирующая; 3 –
фазоманипулированная; 4 – ретранслированная; 5 – гармоническая.
а – 4-КАМ и 4-ФМ; б – 16-КАМ и 16-ФМ; в – 32-КАМ и 32-ФМ.
 76

Из сравнения этих графиков следует, что в радиоканалах со сложной

помеховой обстановкой предпочтение следует отдать сигналам КАМ.
Энергетический выигрыш при приеме по сравнению со случаем М-ФМ во
всех случаях составляет более 3 дБ при одинаковой позиционности сигналов.
На рис. 2.24 показаны зависимости вероятности битовой ошибки от
интенсивности помех для КАМ и М-ФМ при М=16.

Рис. 2.24. Зависимость вероятности битовой ошибки от интенсивности

помех при М=16, Еb/N0= Еbср/N0=13 dB для КАМ (штриховые линии) и М-ФМ
(сплошные линии): 1 – сканирующая; 2 – фазоманипулированная; 3 –
ретранслированная; 4 – гармоническая.

Видно, что при малом уровне помех μ<0,3 и одинаковой

позиционности сигналов когерентные приемники КАМ обладают большей
помехоустойчивостью, чем аналогичные приемники М-ФМ. Вероятность
ошибки в первом случае на один-два порядка меньше. В обоих случаях при
μ>0,3 качественный прием сигналов становится невозможным без
использования специальных мер по борьбе с этими помехами.
 77

2.7. Выводы по разделу 2

Полученные в разделе 2 результаты позволяют:

1. Оценить вероятность битовой ошибки когерентного приемника
сигналов КАМ при наличии в канале связи наряду с шумовой помехой
нефлуктуационных сигналоподобных помех разного вида:
гармонической, сканирующей, фазоманипулированной и
ретранслированной.
2. Определить виды нефлуктуационных помех, оказывающих наибольшее
влияние на помехоустойчивость сигналов КАМ.
Наиболее опасной с точки зрения снижения помехоустойчивости
приема сигналов с квадратурной амплитудной модуляцией является
прицельная гармоническая помеха. Далее опасность нефлуктуационных
помех убывает следующим образом: фазоманипулированная,
ретранслированная, сканирующая.
При разработке методов и устройств приема сигналов КАМ в сложной
помеховой обстановке необходимо предусмотреть синтез алгоритмов
подавления таких помех. При этом в первую очередь необходимо обратить
внимание на синтез алгоритмов, направленных на подавление узкополосных
(гармонических) помех, обладающих наибольшим мешающим эффектом.
Материалы раздела 2 опубликованы в [49,51,54,58-62,80].
Подобная методика расчета вероятности битовой ошибки при наличии
нефлуктуационных помех использована научной группой кафедры при
анализе помехоустойчивости приема сигналов М-ФМ в работах [54-57,63].
Существенным отличием приведенной методики от ранее использованной
является тот факт, что пороги сравнения корреляционных интегралов в
приемнике сигналов КАМ не являются нулевыми и, кроме того, в данном
случае проводится усреднение по всем комбинациям канальных символов.
Характерно, что в частном случае при М = 4 результаты расчетов в
обоих случаях совпадают.
 78

3. ВЛИЯНИЕ ПОГРЕШНОСТЕЙ ФОРМИРОВАНИЯ ОПОРНЫХ

КВАДРАТУРНЫХ КОЛЕБАНИЙ НА ПОМЕХОУСТОЙЧИВОСТЬ
КОГЕРЕНТНОГО ПРИЕМА СИГНАЛОВ КАМ

Важной составной частью схем когерентных демодуляторов сигналов

КАМ является блок формирования опорных квадратурных колебаний и блок
тактовой синхронизации. Это относится как к схеме приемника,
построенного по квадратурной схеме (рис. 1.10), так и к схеме
многоканального когерентного приемника с решением по максимуму
правдоподобия (рис. 1.9). Данные блоки строятся на основе классических
систем ФАПЧ. Однако, любая система ФАПЧ имеет погрешности, которые
приводят к неточностям формирования опорных колебаний и тактовых
сигналов. Эти погрешности, в свою очередь, вызывают ошибки при
демодуляции принимаемого сигнала, которые могут приводить к
возрастанию вероятности ошибки. В данном разделе анализируется влияние
различных погрешностей формирования опорных квадратурных колебаний и
тактовых сигналов на помехоустойчивость когерентного приема сигналов
КАМ.

3.1. Исследование влияния погрешностей оценивания фазы несущей на

помехоустойчивость приема сигналов КАМ

Рассмотрим работу когерентного приемника сигнала (1.2),

реализующего алгоритм по максимуму правдоподобия, и оценим влияние
погрешностей фазы несущей на величину вероятности ошибки. Для этого
используем методику анализа помехоустойчивости, приведенную в п. 2.1,
применительно к случаю приема сигнала на фоне только белого шума (1.9). В
этом случае
x(t )  si (t )  n(t ).
 79

Приемник вычисляет М интегралов свертки Ji принимаемого процесса

со всеми копиями сигнала – опорными колебаниями:

2 Aср Ts
Ji 
N0  x(t )s
0
опi (t )dt.

Неточности оценивания фазы несущей, приводящие к фазовым

погрешностям формирования опорных колебаний, изменяют статистические
параметры распределений на выходах всех корреляторов приемника –
средние значения mmi и дисперсии Dmi – и, соответственно, влияют на
величины вероятностей символьной и битовой ошибки.
При наличии неточности работы системы ФАПЧ по фазе несущей
частоты все опорные сигналы будут иметь некоторый фазовый сдвиг φ
(рис. 3.1):
sопi (t )  Ii cos(ω0t  φ)  Qi sin(ω0t  φ), i  0,1,...M  1.

Рис. 3.1. Пример искажения сигнального созвездия 16-КАМ при наличии

ошибки оценивания фазы несущего колебания.

Вычислим средние значения и дисперсии (в качестве примера расчет

показан для сигнала с индексом «0»):
 80

Ts
2 Aср
m0  J 0 
N0   x(t )  ( I
0
0 cos(ω0t  φ)  Q0 sin(ω0t  φ))dt 

Ts
2 Aср

N0   s (t )  n(t )  ( I
0
0 0 cos(ω0t  φ)  Q0 sin(ω0t  φ))dt 
2
Aср Ts ( I 0 2  Q0 2 )
 cos φ.
N0
Ts
2 Aср
mi  J i 
N0   x(t )  ( I
0
i cos(ω0t  φ)  Qi sin(ω0t  φ))dt 
2
Aср Ts
 (( I 0 I i  Q0Qi ) cos φ  ( I 0Qi  I i Q0 ) sin φ).
N0

2 Esср
(( I 0  Q0  I 0 I i  Q0Qi ) cos φ  ( I 0Qi  I iQ0 )sin φ). (3.1)
2 2
m0  mi 
N0

D0  J 0 2    J 0  2  J 0 2   m0 2 
2 Ts Ts
4 Aср
    x(t1 ) x(t2 )  ( I 0 cos(ω0t1  φ)  Q0 sin(ω0t1  φ)) 
N02 0 0

( I 0 cos(ω0t2  φ)  Q0 sin(ω0t2  φ))dt1dt2  m0 2 

2 Ts Ts
4 Aср N 0
 m0 2    ( I 0 cos(ω0t1  φ)  Q0 sin(ω0t1  φ)) 
N02 2 0 0

2 Esср
( I 0 cos(ω0t2  φ)  Q0 sin(ω0t2  φ))dt1dt2  m0 2  ( I 0 2  Q0 2 ).
N0
2 Esср
Di  J i 2    J i 2  J i 2  mi 2  ( I i2  Qi2 ).
N0
2 Ts Ts
4 Aср
M 0i 
N02    x(t )x(t )  ( I
0 0
1 2 0 cos(ω0t1  φ)  Q0 sin(ω0t1  φ)) 

2 Esср
( I i cos(ω0t2  φ)  Qi sin(ω0t2  φ))dt1dt2  m0 mi  ( I 0 I i  Q0Qi ).
N0
2 Esср 2 Esср 4 Esср
Dmi  D0  Di  2M 0i  ( I 02  Q02 )  ( I i2  Qi2 )  ( I 0 I i  Q0Qi ) 
N0 N0 N0
2 Esср
(3.2)
 (( I 0  I i )  (Q0  Qi ) ).
2 2

N0
 81

Подставим выражения для среднего и дисперсии в (2.6), (2.3), (2.4) и

получим зависимости вероятности битовой ошибки от параметров: отношения
сигнал/шум Ebср/N0 (дБ) и фазового сдвига Δφ. При вычислениях необходимо
провести перебор всех возможных фазовых и амплитудных состояний сигнала
КАМ и учесть нормировку по средней энергии, как в п. 2.1.
Графики зависимостей вероятности битовой ошибки от отношения
сигнал/шум для М = 4, 16, 32, 64 приведены на рис. 3.2, а от величины фазового
сдвига Δφ – на рис. 3.3.

а)
82

б)

в)
 83

г)
Рис. 3.2. Зависимости вероятности битовой ошибки от отношения
сигнал/шум (дБ) при наличии фазовой ошибки формирования опорных
колебаний:
а – 4-КАМ; б – 16-КАМ; в – 32-КАМ; г – 64-КАМ.

Рис. 3.3. Зависимости вероятности битовой ошибки от фазового сдвига Δφ.

 84

Проведенное моделирование работы когерентного приемника при

фазовой погрешности формирования опорных колебаний (программа
приведена в Приложении П1) подтвердило правильность теоретических
результатов. На рис. 3.2 дополнительно показаны точки, полученные в
результате компьютерного эксперимента.
Из графиков видно, что неточность оценивания фазы несущей и,
соответственно, фазовая погрешность формирования опорных колебаний,
может сильно снижать помехоустойчивость когерентного приема сигналов
КАМ, и с увеличением позиционности сигналов это влияние усиливается.
Так, фазовая погрешность Δφ = 0,15 радиан (8,50) для Peb = 10-2 при М = 4
вызывает дополнительные энергетические потери около 1 дБ, при М = 16 –
около 3 дБ, при М = 32 – около 5 дБ, а при М = 64 – около 9 дБ.

3.2. Исследование влияния погрешностей тактовой синхронизации на

помехоустойчивость приема сигналов КАМ

Оценим влияние погрешностей тактовой синхронизации на величину

вероятности ошибки когерентного приемника сигнала КАМ, реализующего
алгоритм по максимуму правдоподобия. Для этого, как и ранее, рассчитаем
статистические параметры распределений на выходах всех корреляторов
приемника.
Предположим, что тактовая частота в приемнике соответствует
номинальной, но все тактовые моменты смещены во времени на одинаковую
величину ξ (статическое смещение) (рис. 3.4).

Sk Si Sj
t
Ts ξ Ts ξ Ts

Рис. 3.4. Временное рассогласование при статистическом смещении тактовых

моментов.
 85

В этом случае в вычислении интегралов свертки будет использоваться

реализация сигнала КАМ:
 Aср ( I i cosω0t  Qi sin ω0t ); t  (ξ, Ts )
s(t )  
 Aср ( I j cosω0t  Q j sin ω0t ) j ); t [Ts , Ts  ξ)
и опорные колебания:
sопi (t )  Aср ( I i cosω0t  Qi sin ω0t ), t  (ξ,Ts  ξ).

Интегралы свертки будут иметь вид:

2 Aср Ts ξ
Ii 
N0 
ξ
x(t )sопi (t )dt.

На выходах корреляторов средние значения m0 и mi определяются

следующим образом:
2 Aср Ts ξ 2 Aср Ts
m0  J 0 
N0 
ξ
 x(t ) ( I 0 cos ω0t  Q0 sin ω0t )dt 
N0   s (t )  (I
ξ
0 0 cos ω0t  Q0 sin ω0t )dt 

Ts  ξ
2 Aср 2 Es ср ξ ξ

N0 
Ts
 s j (t )  ( I 0 cos ω0t  Q0 sin ω0t )dt 
N0
(( I 0 2  Q0 2 )(1 
Ts
)  ( I 0 I j  Q0Q j ) ).
Ts
2 Aср Ts ξ 2 Aср Ts
mi  I i 
N0   x(t )  ( I cos ω t  Q sin ω t )dt 
ξ
i 0 i 0
N0   s (t )  (I cos ω t  Q sin ω t )dt 
ξ
0 i 0 i 0

2 Aср Ts ξ 2 Es ср ξ ξ

N0 
Ts
 s j (t )  ( I i cos ω0t  Qi sin ω0t )dt 
N0
(( I 0 I i  Q0Qi )(1 
Ts
)  ( I i I j  QiQ j ) ).
Ts

2 Es ср 
m0  mi  {( I 0 2  Q0 2  I 0 I i  Q0Qi )  [( I 0  I j )( I 0  I i )  (Q0  Q j )(Q0  Qi )}. (3.3)
N0 Ts

Дисперсии и смешанные моменты определяются по формулам:

2 Es ср
D0  J 0 2    J 0 2  J 0 2  m0 2  ( I 0 2  Q0 2 ),
N0

2 Es ср
Di  J i 2    J i 2  J i 2  mi 2  ( I i 2  Qi 2 ),
N0

2 Es ср
M 0i  ( I 0 I i  Q0Qi ).
N0
 86

2 Es ср
Dmi  D0  Di  2M 0i  (( I 0  I i )2  (Q0  Qi )2 ). (3.4)
N0
Подставив полученные выражения в (2.6), (2.3), (2.4), получим

зависимости вероятности битовой ошибки от параметров: отношения

сигнал/шум Ebср/N0 (дБ) и величины относительного смещения тактовых

моментов ξ/Ts. При вычислениях, как и ранее, проведен перебор всех

возможных фазовых и амплитудных состояний сигнала КАМ и учтена

нормировка по средней энергии.

Соответствующие графики для М = 4, 16, 32, 64 приведены на рис. 3.5 и

3.6.

а)
87

б)

в)
 88

г)
Рис. 3.5. Зависимости вероятности битовой ошибки от отношения
сигнал/шум (дБ) при наличии смещения тактовых моментов:
а – 4-КАМ; б – 16-КАМ; в – 32-КАМ; г – 64-КАМ.

Рис. 3.6. Зависимости вероятности битовой ошибки от величины смещения

тактовых моментов ξ/Ts.
 89

Проведенное моделирование работы когерентного приемника при

погрешностях работы системы тактовой синхронизации (программа
приведена в Приложении П1) подтвердило правильность теоретических
результатов. На рис. 3.5 дополнительно показаны точки, полученные в
результате компьютерного эксперимента.
Из приведенных графиков видно, что смещение тактовых моментов
также сильно влияет на помехоустойчивость когерентного приема сигналов
КАМ. Величина относительного смещения всего в 3% (ξ/Ts = 0,03) для Peb =
10-2 при М = 4 вызывает энергетические потери около 0,4 дБ, при М = 16 –
около 0,6 дБ, при М = 32 – около 0,8 дБ, а при М = 64 – более 1 дБ.

В табл. 3.1 приведены сводные данные по допустимым погрешностям

синхронизации при приеме сигналов КАМ.

Таблица 3.1. Допустимые погрешности оценивания фазы несущей и

тактовой синхронизации (при энергетических потерях 0,5 дБ).

М Допустимая фазовая Допустимая погрешность тактовой

погрешность Δφ синхронизации ξ/Ts
4 3o (π/60) 0,05
16 2o (π/90) 0,04
32 1,5o (π/120) 0,03
64 1o (π/180) 0,02

Рассмотрим влияние амплитудного и фазового разбаланса квадратур,

возникающего в блоке формирования опорных колебаний.

3.3. Исследование влияния амплитудного разбаланса квадратур на

помехоустойчивость приема сигналов КАМ

При наличии амплитудного разбаланса квадратур опорные сигналы в

схеме когерентного приемника будут отличаться от идеальных. Учтем этот
факт введением коэффициента амплитудного разбаланса а:
sопi (t )  Ii cosω0t  aQi sin ω0t , i  0,1,...M  1.
 90

Вычислим средние значения и дисперсии (в качестве примера расчет

показан для сигнала с индексом «0»):
Ts
2 Aоп
m0  J 0 
N0   x(t )  ( I
0
0 cos ω0t  aQ0 sin ω0t )dt 

Ts
2 Aоп
N 0 0
  s0 (t )  n(t )  ( I 0 cos ω0t  aQ0 sin ω0t )dt 

Aср Aоп 2 Aср Aоп 2 2 Esср

 I 0 Ts  aQ0 Ts  ( I 0 2  aQ0 2 ).
N0 N0 N0

Здесь и далее для простоты записей принято Aср  Aоп .

Ts
2 Aоп
mi  J i 
N0   x(t )  ( I
0
i cos ω0t  aQi sin ω0t )dt 

Aср Aоп Aср Aоп 2 Esср 2 Esср

 I 0 I iTs  aQ0QiTs  I 0 Ii  aQ0Qi 
N0 N0 N0 N0
2 Esср
 ( I 0 I i  aQ0Qi ).
N0

2 Esср
m0  mi  ( I 0 ( I 0  I i )  aQ0 (Q0  Qi )). (3.5)
N0

D0  J 0 2    J 0  2  J 0 2   m0 2 
Ts Ts
2
4 Aоп
    x(t1 ) x(t2 )  ( I 0 cos ω0t1  aQ0 sin ω0t1 ) 
N02 0 0

( I 0 cos ω0t2  aQ0 sin ω0t2 )dt1dt2  m0 2 

Ts Ts
2
4 Aоп N 0
 m0 2    ( I 0 cos ω0t1  aQ0 sin ω0t1 ) 
N02 2 0 0

2 Esср
( I 0 cos ω0t2  aQ0 sin ω0t2 )dt1dt2  m0 2  ( I 0 2  a 2Q0 2 ).
N0
2 Esср
Di  J i 2    J i 2  J i 2  mi 2  ( I i2  a 2Qi2 ).
N0
2 Ts Ts
4A
M 0i  оп2
N0    x(t )x(t )  ( I
0 0
1 2 0 cos ω0t1  aQ0 sin ω0t1 ) 

2 Esср
( I i cos ω0t2  aQi sin ω0t2 )dt1dt2  m0 mi  ( I 0 I i  a 2Q0Qi ).
N0
 91

2 Esср 2 Esср 4 Esср

D0i  D0  Di  2M 0i  ( I 0 2  a 2Q0 2 )  ( I i2  a 2Qi2 )  ( I 0 I i  a 2Q0Qi ) 
N0 N0 N0
(3.6)
2 Esср
 (( I 0  I i )  a (Q0  Qi ) ).
2 2 2

N0

Подставив полученные выражения в (2.6), (2.3), (2.4), получим

зависимости вероятности битовой ошибки от параметров: коэффициента
амплитудного разбаланса а и отношения сигнал/шум Ebср/N0 (дБ). Заметим,
что при вычислениях необходимо провести перебор всех возможных
состояний сигнала КАМ и учесть нормировку по средней энергии.
Соответствующие графики для М = 4, 16, 32, 64 приведены на рис. 3.7 и
3.8.

а)
92

б)

в)
 93

г)
Рис. 3.7. Зависимости вероятности битовой ошибки коэффициента
амплитудного разбаланса а:
а – 4-КАМ; б – 16-КАМ; в – 32-КАМ; г – 64-КАМ.

а)
94

б)

в)
 95

г)
Рис. 3.8. Зависимости вероятности битовой ошибки от отношения
сигнал/шум (дБ) при наличии амплитудного разбаланса квадратур:
а – 4-КАМ; б – 16-КАМ; в – 32-КАМ; г – 64-КАМ.

Полученные результаты позволяют утверждать, что при М=16, 32 и 64

отклонение амплитуды одного из опорных квадратурных колебаний от
номинального значения (а=1) приводит к резкому снижению
помехоустойчивости. Допустимым отклонением можно считать величину 3-
5%, при этом значение Peb при Ebср/N0 = 13 дБ возрастает в 2 раза, а
соответствующий энергетический проигрыш в отношении сигнал/шум
составляет от 0,2 до 0,5 дБ.

3.4. Исследование влияния фазового разбаланса квадратур на

помехоустойчивость приема сигналов КАМ

При наличии фазового разбаланса опорный сигнал в одной из квадратур

будет иметь некоторый дополнительный фазовый сдвиг φ :
sопi (t )  Ii cosω0t  Qi sin(ω0t  φ), i  0,1,...M  1.
 96

Вычислим средние значения и дисперсии:

T
2A s
m0  J 0  оп   x(t )  ( I 0 cos ω0t  Q0 sin(ω0t  φ))dt 
N0 0
T
2A s
 оп   s0 (t )  n(t )  ( I 0 cos ω0t  Q0 sin(ω0t  φ))dt 
N0 0
2 Esср
 ( I 0 2  Q0 2 cos φ  I 0Q0 sin φ).
N0

T
2 Aоп s
N 0 0
mi  J i   x(t )  ( I i cos ω0t  Qi sin(ω0t  φ))dt 

2 Esср
 ( I 0 I i  Q0Qi cos φ  I 0Qi sin φ).
N0

2 Esср
m0  mi  ( I 0 ( I 0  I i )  Q0 (Q0  Qi ) cos φ  I 0 (Q0  Qi ) sin φ). (3.7)
N0

2 Esср 2 Esср 4 Esср

D0  Di  2 M 0i  ( I 02  Q02 )  ( I i2  Qi2 )  ( I 0 I i  Q0Qi ) 
N0 N0 N0
(3.8)
2 Esср
 (( I 0  I i )  (Q0  Qi ) ).
2 2

N0

Подставив полученные выражения в (2.6), (2.3), (2.4), получим

зависимости вероятности битовой ошибки от параметров: фазового сдвига φ
и отношения сигнал/шум Ebcp/N0 (дБ). При вычислениях, как и ранее,
проведен перебор всех возможных состояний сигнала КАМ и учтена
нормировка по средней энергии.
Соответствующие графики для М = 4, 16, 32, 64 приведены на рис. 3.9 и
3.10.
 97

Рис. 3.9. Зависимости вероятности битовой ошибки от величины фазового

разбаланса квадратур φ.

а)
98

б)

в)
 99

г)
Рис. 3.10. Зависимости вероятности битовой ошибки от отношения
сигнал/шум (дБ) при наличии фазового разбаланса квадратур:
а – 4-КАМ; б – 16-КАМ; в – 32-КАМ; г – 64-КАМ.

Видно, что фазовый разбаланс квадратур сильно влияет на

помехоустойчивость когерентного приема сигналов КАМ. Так фазовая
погрешность в 0,05 рад (3 градуса) эквивалентна при Peb =10-3 ухудшению
отношения сигнал/шум 0,25 дБ для 4-КАМ, 0,5 дБ для 16-КАМ, 0,6 дБ для 32-
КАМ и около 1 дБ для 64-КАМ. При большей фазовой погрешности
энергетические проигрыши многократно возрастают, особенно при больших М.
В табл. 3.2 приведены сводные данные по величине допустимого
разбаланса квадратур при приеме сигналов КАМ.

Таблица 3.2. Допустимый амплитудный и фазовый разбаланс квадратур

(при энергетических потерях 0,5 дБ).

М Допустимый амплитудный Допустимый фазовый разбаланс

разбаланс
4 - 3o (π/60)
16 10% 3o (π/60)
32 5% ~2o (π/90)
64 ~3% ~1,5o (π/120)
 100

3.5. Сравнение влияния погрешностей синхронизации при приеме

сигналов КАМ и М-ФМ

Подобная задача с использованием аналогичного математического

аппарата применительно к сигналам с многопозиционной фазовой
манипуляцией решена в [64-66]. Проведем сравнение полученных
результатов.
На рис. 3.11 приведены сравнительные зависимости вероятности
битовой ошибки Peb от величины фазовой погрешности формирования
опорных колебаний при Ebср/N0=13 дБ для КАМ и М-ФМ одинаковой
позиционности и с одинаковой средней энергией.

Рис. 3.11. Зависимость вероятности битовой ошибки Peb от величины фазовой

погрешности Δφ для КАМ (сплошные линии) и М-ФМ (штриховые линии).

Сравнительные зависимости вероятности битовой ошибки Peb от

величины погрешности тактовой синхронизации τ/Ts при Ebср/N0=13 дБ для
КАМ и М-ФМ одинаковой позиционности и с одинаковой средней энергией
показаны на рис. 3.12.
 101

Рис. 3.12. Зависимость вероятности битовой ошибки Peb от величины

погрешности тактовой синхронизации τ/Ts для КАМ (сплошные линии) и М-
ФМ (штриховые линии).

Видно, что погрешности оценивания фазы опорных колебаний и

тактовой синхронизации и для сигналов М-ФМ, и для сигналов КАМ
похожим образом влияют на помехоустойчивость когерентного приема.
Более точные расчеты показывают, что для приемника М-ФМ это влияние
выражено сильнее.

3.6. Выводы по разделу 3

Проведенное исследование влияния фазовых погрешностей при

формировании опорных колебаний и погрешностей тактовой синхронизации
при когерентном приеме сигналов КАМ позволяет сделать следующие
выводы:
1. Указанные погрешности систем фазовой и тактовой синхронизации
могут сильно снижать помехоустойчивость приема, и с увеличением
позиционности сигналов это влияние усиливается.
 102

2. Если допустить величину энергетических потерь при приеме, равную

0,5 дБ за счет каждой из этих погрешностей, то для Peb = 10-2
допустимая фазовая погрешность Δφ должна быть от ~30 при М = 4
до ~10 при М = 64, а допустимая погрешность тактовой
синхронизации ξ/Ts, соответственно, от ~5% при М = 4 до ~2% при М
= 64 (табл. 3.1).
3. Для обеспечения более жестких требований к величине потерь и
величине Peb требования к указанным погрешностям существенно
возрастают.

Исследование влияния амплитудного и фазового разбаланса квадратур

при когерентном приеме сигналов КАМ позволяет сделать следующие
выводы:
1. Амплитудный разбаланс квадратур при М > 16 ведет в
существенному снижению помехоустойчивости приема сигналов
КАМ. Допустимым отклонением амплитуды (при энергетических
потерях 0,5 дБ) в этом случае можно считать величину 3-5%. При
М=4 амплитудный разбаланс в широком диапазоне значений
практически не сказывается на помехоустойчивости (табл. 3.2).
2. Фазовый разбаланс квадратур сильно влияет на помехоустойчивость
когерентного приема сигналов КАМ. С увеличением позиционности
сигналов это влияние также усиливается. Если для 4-КАМ и 16-КАМ
можно считать допустимой фазовую погрешность (при
энергетических потерях 0,5 дБ) в 0,05 рад (3 градуса), то 64-КАМ эта
погрешность уже приводит в существенным энергетическим потерям
(около 1 дБ). При большей фазовой погрешности энергетические
проигрыши многократно возрастают, особенно при больших М (табл.
3.2).

Материалы раздела 3 опубликованы в [67, 77, 78].

 103

4. РАЗРАБОТКА АЛГОРИТМОВ БОРЬБЫ С

НЕФЛУКТУАЦИОННЫМИ ПОМЕХАМИ ПРИ ПРИЕМЕ СИГНАЛОВ
КАМ

Проведенный в разделе 2 теоретический анализ вероятности битовой

ошибки при приеме сигналов КАМ, показал, что качество приема
значительно снижается, если в канале связи кроме шумовой помехи
присутствуют еще и нефлуктуационные помехи разных видов. Особенно это
касается случаев, когда такая помеха является узкополосной
(гармонической), совпадающей по частоте с частотой полезного сигнала.
Такое снижение объясняется тем, что алгоритмы работы демодулятора
оптимизированы для приема сигналов на фоне белого гауссовского шума и
не учитывают характеристики нефлуктуационных помех [68]. В связи с этим
в данном разделе решается задача разработки и анализа алгоритмов приема,
позволяющих осуществить оптимальную обработку или подавление
нефлуктуационных помех.

4.1. Оптимальный прием сигналов КАМ на фоне гармонической помехи

со случайной начальной фазой

Проведем синтез и анализ оптимального алгоритма приема

многопозиционных сигналов КАМ с некогерентной обработкой
гармонической помехи.
Для этого рассмотрим оптимальный прием сигнала КАМ (1.2) на фоне
белого гауссовского шума n(t) с односторонней спектральной плотностью
мощности N0 и гармонической помехи
sп (t ,φп )  μAср cos(ωпt  φп ), (4.1)

с относительной интенсивностью µ, частотой ωп  2πfп , близкой к частоте

полезного сигнала, и со случайной начальной фазой φп. При этом полагаем,
что начальная фаза полезного сигнала известна точно, а распределение
 104

случайной величины φп является равномерным в диапазоне (-π,π]. В этом

случае можно говорить о некогерентной обработке помехи.
Процесс, поступающий на вход приемника:
x(t )  si (t )  sп (t , φп )  n(t )

представляет собой аддитивную смесь помехи, шума и сигнала КАМ (1.2).

Для простоты дальнейших записей примем начальную фазу сигнала КАМ
равной нулю.
Введем обозначение для суммы сигнала и помехи:
sсп (t , Ci ,φп )  si (t )  sп (t ,φп ). (4.2)
С учетом этого совместную апостериорную вероятность канального
символа Ci и фазы φп можно определить следующим образом [52,53]:
p ps [ sсп (t , Ci , φп )]  p ps (Ci , φп )  Kp pr (Ci ) p pr (φ п ) 
T Ts
1 s 2 2
 exp[ 
N0 0
sсп (t , Ci , φп )dt 
N0  x(t )s
0
сп (t , Ci , φ п )dt ].

Здесь p pr (Ci ) – априорная вероятность канального символа Ci,

p pr (φп )  1/ 2π , K – коэффициент, включающий в себя все составляющие,

не содержащие информации о полезном сигнале и нефлуктуационной

помехе.
Преобразуем первый интеграл под экспонентой:
Ts Ts Ts Ts

s (t , Ci , φп )dt   si (t , Ci )dt   sп (t , φп )dt  2  si (t , Ci ) sп (t , φ п )dt 

2 2 2
сп
0 0 0 0

 Es  Eп  2 Ri (φп ),
где Es , Eп – энергии посылок сигнала и помехи длительности Тs, Ri (φп ) –
коэффициент корреляции между сигналом и помехой.
Следовательно, искомая апостериорная вероятность равна:
 105

T
1 1 2 s
p ps (Ci , φп )  Kp pr (Ci ) exp[ ( Ei  Eп  2 Ri (φп ))]exp[  x(t ) sсп (t , Ci ,φ п )dt ] 
2π N0 N0 0
T
1 2 2 s
 K1 exp[ Ri (φп )]exp[  x(t ) sсп (t , Ci ,φ п )dt ].
2π N0 N0 0
Здесь в коэффициент K1 включены члены, не зависящие от значения

символа Сi и от фазы φ п .
Для получения апостериорной вероятности дискретного символа Сi
необходимо усреднить полученную величину по фазе помехи φ п . Получим:
π T
1 2 2 s
2π π
p ps (Ci )  K1 (exp[ Ri (φп )]exp[  x(t )sсп (t , Ci ,φп )dt ])dφп .
N0 N0 0
Введем обозначения для интегралов, определяющих величины
корреляции между процессом x(t) и копиями сигналов КАМ:
T
2 s
N 0 0
J0  x(t ) si (t , Ci  0)dt;

….. (4.3)
T
2 s
N 0 0
J M 1  x(t ) si (t , Ci  M  1)dt.

Аналогично введем обозначения для интегралов, определяющих

величины корреляции между образцом принимаемой помехи,
сформированным путем вычитания копий сигнала из принимаемого процесса
x(t), и копией помехи sп (t , φп ) :
T
2 s
N0 0
a0  [ x(t )  si (t , Ci  0)]sп (t , φп )dt ,

….. (4.4)
T
2 s
N0 0
aM 1  [ x(t )  si (t , Ci  M  1)]sп (t , φп ).
 106

Тогда с учетом (4.1), (4.2)-(4.4) алгоритм принятия решения о значении

канального символа Сi можно записать следующим образом с
использованием модифицированных функций Бесселя I0 (·):

 
π
1
Ci  max  p ps (Ci )  max exp( J i )  exp(ai )dφп   max exp( J i ) I 0 (U i ) ,
 2π  π  (4.5)

p ps (Ci  0)  p ps (C j  0)
j i  , i, j  0,..., M 1.
Последнее выражение можно записать по-другому:


exp( J i ) I 0 (U i )  exp( J j ) I 0 (U j )
j i  (4.6)

или


J i  ln( I 0 (U i ))  J j  ln( I 0 (U j ))
j i , i, j  0,..., M 1.
При сравнении вычисляемых величин для нахождении максимума
max{·} необходимо учитывать, что пороги принятия решения не нулевые, а
равны полуразности энергий соответствующих посылок сигнала (2.4).
Аргументы функций Бесселя вычисляются через квадратурные
составляющие величин (4.4):

U i  X i 2  Yi 2 ;

2μAср Ts
Xi 
N0  [ x(t )  s (t , C )]cos ω tdt;
0
i i п

2μAср Ts (4.7)
Yi 
N0  [ x(t )  s (t , C )]sin ω tdt.
0
i i п

Структурная схема приемника, реализующая эти процедуры, приведена

на рис. 4.1. В части оценки канального символа она в основном совпадает с
решающим правилом для приемника сигналов КАМ на фоне только белого
гауссовского шума. Но кроме схемы вычисления корреляционных
интегралов 8 дополнительно имеются схемы вычисления весовых
коэффициентов (в виде бесселевых функций), величина которых
определяется степенью корреляции помехового колебания в принимаемой
 107

смеси x(t) с копией помехи в приемнике. Для примера на рис. 4.1 показан
только одна дополнительная схема 9, остальные строятся аналогично.

Рис. 4.1. Структурная схема оптимального приемника сигналов КАМ с

некогерентной обработкой гармонической помехи: 1 – фазовращатель на 900,
2 − интегратор, 3 – блок вычисления экспоненты, 4 – блок возведения в
квадрат, 5 – блок вычисления модуля, 6 − формирователь весовых
коэффициентов, 7 − блок выбора максимума (с учетом порога), 8 − схема
оценки канального символа, 9 – схема формирования весовых
коэффициентов, 10 – генератор копии помехи, 11 – блок формирования
опорных колебаний.

Для оценки эффективности схемы некогерентной обработки

гармонической помехи было проведено моделирование приема для сигналов
4-КАМ и 16-КАМ. Определялась вероятность битовой ошибки Peb в
зависимости от параметров помехи и настроек приемника.
 108

А. Зависимость вероятности битовой ошибки от расстройки помехи.

Моделирование оптимального приемника сигналов КАМ показало, что

эффективность работы схемы некогерентной обработки помехи зависит от
расстройки ∆fп частоты помехи относительно центральной частоты спектра
полезного сигнала. Рис. 4.2 иллюстрирует зависимость вероятности битовой
ошибки Peb от расстройки помехи ∆fпTs при µ = 0,5 и отношении сигнал/шум
Ebср / N0 =7дБ. Видно, наибольший выигрыш по помехоустойчивости

наблюдается при расстройке ∆fпTs = 0,5. Так, для 4-КАМ вероятность ошибки
снижается в 2,5 раза, для 16-КАМ – в 2,8 раза. При прицельной помехе (∆fпTs
= 0) схема оказывается неэффективной, а при ∆fпTs = 1 выигрыш отсутствует
в связи с тем, что частота помехи совпадает с частотой первого нуля спектра
сигнала, и такая помеха не влияет на помехоустойчивость приема.
Рассмотрим более подробно случай, когда частота помехи не совпадает
с центральной частотой спектра полезного сигнала, и, как следует из рис. 4.2,
схема некогерентной обработки помехи дает выигрыш в
помехоустойчивости.

а)
 109

б)
Рис. 4.2. Зависимости вероятности битовой ошибки от расстройки помехи с
обработкой (1) и без обработки помехи (2): а – 4-КАМ, б – 16-КАМ.

Б. Зависимость вероятности битовой ошибки от отношения сигнал/шум.

На рис. 4.3 показаны зависимость вероятности битовой ошибки Peb от

отношения сигнал/шум (дБ) при ∆fпTs=0,5. Штриховые линии соответствуют
случаю включенной схемы некогерентной обработки помехи. Для сравнения
сплошными линиями даны кривые, полученные без обработки помехи.
Наглядно виден положительный эффект такой обработки – вероятность
битовой ошибки значительно снижается. Особенно это заметно при больших
значениях интенсивности помехи μ=0,5 и 0,9. Так при μ=0,5 и Peb=10-2 для 4-
КАМ включение схемы некогерентной обработки помехи дает
энергетический выигрыш около 2 дБ, для 16-КАМ – более 5 дБ.
 110

а)

б)
Рис. 4.3. Зависимости вероятности битовой ошибки от отношения
сигнал/шум (дБ) при µ=0 (1), µ=0,1 (2), µ=0,5 (3), µ=0,9 (4): а – 4-КАМ, б –
16-КАМ.

В. Влияние неточности установки в приемнике частоты копии помехи.

Алгоритм (4.6) предполагает знание в приемнике частоты помехи, как

образца. Естественно, что это предположение часто не выполняется.
Рассмотрим, как влияет отклонение частоты копии помехи Δfппр,
установленной в приемнике, на величину вероятности ошибки Peb. На рис.
4.4 представлены графики, иллюстрирующие эту зависимость для разных
 111

сигналов при Ebср/N0 = 7 дБ. Здесь по оси абсцисс отложена величина ΔfппрTs,
нормированная к символьной скорости передачи полезного сигнала.
Естественно, что обработка помехи является наиболее эффективной,
когда в приемнике точно известна частота помехи ΔfппрTs=0. Небольшие
неточности установки в приемнике частоты копии помехи ΔfппрTs=±0,1
незначительно снижают помехоустойчивость приема сигнала, причем это
снижение тем сильнее, чем больше интенсивность помехи, но при μ=0,9 и
расстройке ΔfппрTs =±0,5 вероятность битовой ошибки может увеличиваться
уже на порядок.

а)

б)
Рис. 4.4. Влияние неточности установки в приемнике частоты копии помехи
при µ=0,5 (1), µ=0,9 (2): а – 4-КАМ, б – 16-КАМ.
 112

Г. Влияние неточности установки в приемнике уровня копии помехи.

Алгоритм (4.6) предполагает знание в приемнике также и уровня помехи.

Рассмотрим, как влияет отклонение уровня копии помехи Δµппр,
установленного в приемнике, на величину вероятности битовой ошибки Peb. На
рис. 4.5 представлены графики, иллюстрирующие эту зависимость для разных
сигналов при Ebср/N0 = 7 дБ и номинальной интенсивности помехи μ=0,9.
Из рис. 4.5 видно, что величина вероятности ошибки минимальна при
отклонении интенсивности ∆µппр= 0. Небольшие неточности установки в
приемнике уровня копии помехи ∆µппр= ±0,1 незначительно снижают
помехоустойчивость приема сигнала, но при ∆µппр= ±0,4 наблюдается
увеличение вероятности битовой ошибки в 3–5 раз.

Рис 4.5. Влияние неточности установки в приемнике уровня копии помехи

для 4-КАМ (1), 16-КАМ (2).
По результатам анализа полученных данных можно сделать следующие
выводы:

1. Оптимальная некогерентная обработка гармонической помехи

повышает помехоустойчивость приема многопозиционных сигналов
КАМ за счет использования дополнительной схемы коррекции порогов
в блоке принятия решения приемника.
2. Эффективность работы схемы зависит от расстройки частоты
гармонической помехи относительно центральной частоты спектра
 113

полезного сигнала. Наибольший выигрыш в помехоустойчивости

наблюдается при расстройке ∆fпTs = 0.5.
3. Выигрыш в помехоустойчивости приема тем выше, чем больше
интенсивность гармонической помехи.
4. Схема некогерентной обработка гармонической помехи сохраняет
работоспособность и при небольших неточностях установки частоты
(ΔfппрTs=±0,1) и уровня (∆µппр=±0,1) копии помехи в приемнике.
Аналогичная методика синтеза алгоритма компенсации влияния
гармонической помехи была использована в научной группе кафедры для
сигналов М-ФМ. При М = 4 полученные результаты в обоих случаях совпали.

4.2. Адаптивная фильтрация нефлуктуационных помех при приеме

сигналов КАМ
В п. 4.1 был рассмотрен оптимальный алгоритм компенсации влияния
нефлуктуационной гармонической помехи на прием сигналов КАМ, при
реализации которого считаются известными (или известными с
определенной точностью) частота и амплитуда помехи. На практике такие
условия крайне редки. В случае, когда априорных сведений о помехах мало,
эффективным способом борьбы с ними является использование адаптивных
фильтров (АФ) – нерекурсивных цифровых фильтров с регулируемыми
весовыми коэффициентами (ВК).
Структурная схема такого фильтра представлена на рис. 4.6.

Рис. 4.6. Структура нерекурсивного цифрового фильтра.

 114

Элемент выходной последовательности такого фильтра в пространстве

комплексных огибающих можно записать следующим образом:
. . T .
y i  Wi Xi ,
. . . .
где Xi  [ xi , xi 1 ,... xi  N 1 ] – вектор последовательности входных отсчетов
. . i . i . i
фильтра; Wi  [ w0 ,... wk , wN 1 ] – вектор ВК в i-й момент времени; т – знак
т

транспонирования; N – длина фильтра.

АФ имеют разные алгоритмы настройки вектора ВК, при этом методы
адаптации разделяются на два типа [69-71]:
- адаптация с эталонной моделью, когда для формирования на
приемной стороне эталонной реакции системы при передаче сигнала
используется дополнительная обучающая последовательность;
- адаптация без эталонной модели, когда дополнительная обучающая
последовательность не требуется, а вектор ВК подстраивается в процессе
приема сигнала по результату анализа реакции фильтра на соответствие
выбранному критерию.
В качестве такого критерия могут быть использованы характерные
свойства полезного сигнала, отличающие его от помехи. При этом в
принимаемой смеси «сигнал + помеха» будет максимизироваться сигнальная
составляющая и подавляться помеховая.
В [72] предложены схема и два алгоритма работы АФ, совместно
использующие априорную информацию о двух параметрах полезного
сигнала – постоянной огибающей и фазовой структуре (форме фазового
импульса). Первый алгоритм, изначально в отдельности предложенный в [73]
осуществляет подстройку вектора ВК в соответствии с уравнением:
. . . . . *
Wi 1  Wi  d  (| y i |2 1) y i Xi , (4.8)
где d – коэффициент, определяющий степень инерционности и устойчивость
процесса адаптации; (·)* – знак комплексного сопряжения.
 115

.
При этом минимизируется целевой функционал M{ (| y i |2 1) }, где M {}
– знак математического усреднения.
Такой алгоритм применим к сигналам с постоянной огибающей, к
числу которых сигнал КАМ относится только в частном случае М=4. Этот
случай подробно рассмотрен в [74].
Работа второго алгоритма [75] основана на использовании известных
данных о форме фазового импульса, а именно, о модуле производной
информационной составляющей фазы (наклоне фазового импульса) сигнала.
В случае сигнала КАМ этот параметр равен нулю.
Данный алгоритм контролирует форму (наклон) фазового импульса
принимаемого сигнала в соответствии с уравнением:
. . . . .
Wi 1  Wi  b | y i ' | y i X*i , (4.9)

где b – коэффициент, определяющий степень инерционности и устойчивость

этого процесса адаптации.
В результате работы АФ происходит выравнивание фазового импульса
и, как следствие, уменьшение влияния помехи и повышение точности
фильтрации.
Рассмотрим результаты моделирования алгоритма (4.9).
Моделирование проводилось для сигналов 4-КАМ и 16-КАМ в условиях
наличия шумовой и гармонической помех. Фрагмент программы,
описывающий работу АФ, приведен в Приложении 3.

А. Проверка работоспособности алгоритма.

При моделировании рассмотрен случай, когда частота помехи

совпадает с центральной частотой спектра сигнала, что согласно
исследованиям [38,58] является наихудшей ситуацией с точки зрения
помехоустойчивости приема. При проверке работоспособности алгоритма
оценивалась амплитудно-частотная характеристика (АЧХ) адаптивного
фильтра и изменение ее формы во времени. АЧХ получалась из импульсной
 116

характеристики фильтра – вектора ВК – с помощью процедуры быстрого

преобразования Фурье. Длина фильтра полагалась равной N = 256, так как
предварительная проверка показала, что при значениях N > 200 зависимость
качества фильтрации от N стабилизируется.
Адаптивный фильтр в модели приемника был включен после
квадратурного преобразователя – переноса спектра на нулевую частоту, и его
АЧХ, наблюдаемая в процессе адаптации при указанной выше частоте
помехи, приведена на рис. 4.7.
Как видно, фильтр осуществляет режекцию конвертированной смеси
сигнала и гармонической помехи в узкой полосе в районе нулевой частоты.
Кроме того, наблюдается дополнительный провал АЧХ в более
высокочастотной области, связанный, видимо, с настройкой фильтра на
дополнительные комбинационные частоты, возникающие при оцифровке
сигнала и помехи.

Рис. 4.7. АЧХ адаптивного фильтра.

 117

Спектр S(ω) смеси «сигнал + гармоническая помеха» при µ = 0,9 в

конце процесса адаптации до адаптивного фильтра и после него показан на
рис. 4.8. Видно, что в данном случае помеха подавляется на ~ 8 дБ.

а б
Рис. 4.8. Спектр смеси «сигнал + гармоническая помеха» до (а) и после (б)
адаптивного фильтра.

Динамику процесса адаптации фильтра при разных величинах b

демонстрирует изменение во времени коэффициента передачи фильтра K0 на
нулевой частоте, характеризующего степень подавления гармонической
помехи (рис. 4.9).

Рис. 4.9. Зависимость коэффициента передачи фильтра K0 на нулевой частоте

от времени адаптации (n – количество принятых информационных битов):
4-КАМ (штриховые линии) и 16-КАМ (сплошные линии).
 118

Можно наблюдать разную скорость адаптации фильтра для случаев 4-

КАМ и 16-КАМ, а также тот факт, что при 16-КАМ коэффициент передачи
фильтра K0 в конце адаптации не опускается ниже 0,1, в то время, как при 4-
КАМ он стремится к нулю. В результате анализа этих графиков для модели
была выбрана величина коэффициента b = 10-5.

Б. Зависимость вероятности битовой ошибки от интенсивности

гармонической помехи.
Зависимости вероятности битовой ошибки Peb от интенсивности
гармонической помехи µ при Ebср/N0 = 7 дБ приведены на рис. 4.10. Из
полученных графиков видно, что в случае 4-КАМ с возрастанием
интенсивности помехи возрастает и эффективность ее фильтрации, но для 16-
КАМ адаптивный фильтр улучшает помехоустойчивость приема лишь в
очень малой степени. При малых величинах µ < 0,2 наблюдается даже
ухудшение помехоустойчивости за счет частичного подавления фильтром
полезного сигнала. Кроме этого для сигнала 16-КАМ в этом диапазоне µ
наблюдалась неустойчивость процесса адаптации.

а)
 119

б)
Рис. 4.10. Зависимость вероятности битовой ошибки Peb от
интенсивности гармонической помехи µ при Ebср/N0 = 7 дБ:
а) М = 4, б) М = 16.
Штриховые линии – без АФ; сплошные линии – с АФ (алгоритм (4.9)).

В. Зависимость вероятности битовой ошибки от отношения сигнал/шум.

Зависимости вероятности битовой ошибки Peb от отношения

сигнал/шум Ebср/N0 при наличии гармонической помехи с относительной
интенсивностью µ приведены на рис. 4.11.
Кривые для М = 4 совпадают с результатами, полученными в [76] для
сигнала 4-ФМ, что закономерно. В этом случае использование адаптивного
фильтра существенно повышает помехоустойчивость приема сигналов.
Энергетический выигрыш при µ = 0,3 составляет около 2 дБ, при µ = 0,5
доходит до 4 дБ и увеличивается с возрастанием µ.
При М = 16 исследуемый алгоритм работы АФ, хотя и сохраняет
работоспособность, но его эффективность низка – энергетический выигрыш при
µ = 0,5 составляет ориентировочно от 0,1 до 0,5 дБ. Это объясняется тем, что
сигналы КАМ с М ≥ 16 имеют амплитудную модуляцию, и скачки амплитуды,
возникающие в тактовые моменты, приводят к резким изменениям комплексной
огибающей и ее производной, что нарушает процесс адаптации (4.9).
 120

а)

б)
Рис. 4.11. Зависимости вероятности битовой ошибки Peb от отношения
сигнал/шум (дБ) для: а) М = 4, б) М = 16.
Штриховые линии – без АФ; сплошные линии – с АФ (алгоритм (4.9)).

Аналогичный адаптивный алгоритм компенсации гармонической

 121

помехи был использован в научной группе кафедры для сигналов М-ФМ.

При М = 4 полученные результаты в обоих случаях совпали.

4.3. Выводы по разделу 4

В разделе 4 проведена разработка алгоритмов подавления
нефлуктуационных помех, в частности, наиболее опасной гармонической
помехи, при приеме сигналов КАМ, а именно, синтезирован и
проанализирован оптимальный приемник сигналов КАМ на фоне
гармонической помехи со случайной начальной фазой с некогерентной
обработкой помехи (4.6), (4.7), предложен и проанализирован алгоритм
адаптивной фильтрация нефлуктуационных помех (4.9), следящий за фазовой
траекторией сигнала КАМ.
Схемы имеют разную сложность и эффективность. Так, первый
алгоритм является более сложным в реализации и требует знания на
приемной стороне некоторых параметров гармонической помехи – частоты и
амплитуды, что на практике часто невыполнимо. Кроме того, алгоритм
оказался неэффективным при «прицельной» помехе, частота которой
совпадает с несущей частотой полезного сигнала. В этой связи более
предпочтителен второй алгоритм, который адаптивно подстраивает
характеристики фильтра под параметры гармонической помехи и
осуществляет ее подавление. Так, моделированием установлено, что при
наличии «прицельной» помехи для М = 4 энергетический выигрыш может
достигать единиц дБ, а для М = 16 – величины 0,5 дБ.

Материалы раздела 4 опубликованы в [72, 79].

 122

5. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ

5.1. Измерение отношения сигнал/шум в сети LTE на маршруте

следования железнодорожного состава

Маршрут: Москва - Санкт Петербург

Состав измерительного стенда (рис. 5.1):
- Industrial Cellular Router R3000-4L (Robustel) – 1шт.
- антенна (штатная) AST – 2 шт.
- ноутбук Lenovo с программным обеспечением – 1 шт.

а)
 123

б)

Рис. 5.1. Состав измерительного стенда (а) и схема измерений (б).

Методика измерения отношения сигнал/шум и результаты

1. Измерение уровня шума Рш (dBm) на входе модема (проводится в

экранированной безэховой камере).

Измеренное значение: Рш = -113 dBm.
2. Измерение уровня сигнал+шум Рс+ш (dBm) на входе модема
(проводится на маршруте следования железнодорожного состава).
Измеренные значения: вторая и пятая строки в табл. 5.1.
3. Расчет отношения сигнал/шум q = Рс/Рш (dB) на входе модема по

формуле:
P
q  с  10lg(10 /10  1),
Pш

где  (dB)  Pс+ш  Pш .

Рассчитанные значения: третья и шестая строки в табл. 5.1.

Для наглядности результаты измерений отношения сигнал/шум по

маршруту представлены в виде диаграммы (рис. 5.2).
 (dB)
(dB)
(dBm)

Место
Место

(dBm )
Уровень
Уровень

маршрута
маршрута

Отношение
Отношение

сигнал/шум
сигнал/шум

сигнал+шум
сигнал+шум
Нет сети LTE Лес 23 -90 Москва-Октябрьская
18,9 -94 Лес (Новгородская обл.) 14,9 -98 ст. Лихоборы
9,5 -103 ст. Окуловка 0 -110 ст. Фирсановка
Нет сети LTE Лес 22 -91 ст. Крюково
21 -92 Лес Нет сети LTE Лес
3,3 -108 ст. Торбино Нет сети LTE Лес
Нет сети LTE Лес Нет сети LTE Лес
Нет сети LTE Лес 11,7 -101 Лес
0 -110 разъезд н.п. Бурга 8,4 -104 Лес
124

Нет сети LTE Лес 23 -90 Лес (Московская обл.)

27 -86 Въезд ст. М. Вишера 0 -110 Лес (Тверская обл.)
0 -110 Малая Вишера (вокзал) 26 -87 г. Тверь
25 -88 Выезд ст. М. Вишера Нет сети LTE Лес
Нет сети LTE Лес 14,9 -98 Лес
29 -84 ст. Чудово Нет сети LTE Лес
Нет сети LTE Лес 24 -89 ст. Спирово
23 -90 Лес (Ленинградская обл.) 11,7 -101 г. Вышний Волочек
11,7 -101 Разъезд Нет сети LTE Лес
11,7 -101 ст. Колпино 17,9 -95 Окраина г. Бологое
27 -86 ст. Обухово 4,7 -107 г. Бологое (вокзал)
Таблица 5.1. Результаты измерений отношения сигнал/шум.

16,9 -96 С. Петербург-Московский

 125

35

30

25

20

15

10

0
Маршрут Москва - С. Петербург

Рис. 5.2. Результаты измерений отношения сигнал/шум (дБ) по маршруту.

5.2. Экспериментальное исследование помехоустойчивости приема
сигналов 16-КАМ в присутствии гармонической помехи

Структурная схема экспериментального стенда приведена на рис. 5.3.

Рабочее место включает следующие приборы:
 векторный генератор сигналов N5166B CXG Keysight;
 аналоговый генератор сигналов N5183B MXG Keysight;
 анализатор сигналов (спектра) N9021B MXA Keysight.

Генератор сигналов векторный

N5166B CXG

Σ
Анализатор сигналов
Генератор сигналов аналоговый (спектра) N9021B MXA
N5183B MXG с измерительным
приложением N9054EM0E
Генератор сигналов аналоговый
N5183B MXG
Σ
Рис. 5.3. Структурная схема экспериментального стенда.
 126

Векторный генератор сигналов служит для формирования

информационного радиосигнала со следующими характеристиками:
 модуляция – 16-КАM;
 фильтр RRC с коэффициентом скругления 0,35;
 несущая частота – 3000 МГц;
 частота следования канальных символов – 10 МГц;
 информационная последовательность – ПСП длиной 10240 символов с
циклическим повторением.
Для формирования цифрового информационного радиосигнала
использовалось ПО PathWave N7608APPC.
Аналоговый генератор сигналов служит для формирования
гармонической помехи. В эксперименте использовались следующие значения
частоты гармонической помехи:
 Fп1 = 3000,01 МГц (примерно равна частоте несущей информационного
сигнала);
 Fп2 = 3004,99 МГц (отстоит от несущей информационного сигнала
примерно на половину частоты следования символов (низкое
подавление фильтра RRC));
 Fп3 = 3006,51 МГц (частота на спаде частотной характеристики фильтра
RRC).
Анализатор спектра выполняет роль приемника информационного
радиосигнала, прием осуществлялся с использованием измерительного
приложения для анализа сигналов с цифровыми видами модуляции
N9054EM0E Vector Modulation Analysis (VMA).
Для объединения сигналов генераторов используется сумматор
мощности.
В качестве аддитивного белого гауссовского шума в эксперименте
используются собственные шумы анализатора спектра.
 127

Проведенные измерения
1. Измерена вероятность битовой ошибки от отношения сигнал/шум без
гармонической помехи.
2. Проведена проверка возможности приема в присутствии помехи. Для
частоты помехи Fп1 прием оказался возможен только при
относительном уровне помехи не более -17 дБ относительно
информационного сигнала, для Fп2 – не более -14 дБ, а для Fп3 – не
более +2 дБ.
3. Измерена вероятность битовой ошибки от отношения сигнал/шум в
присутствии гармонической помехи для трех случаев частоты помехи.
Результаты измерений приведены на рис. 5.4.

1,0E-02
Вероятность биитовой ошибки

1,0E-03

1,0E-04

1,0E-05

1,0E-06
15 16 17 18 19 20 21 22
Отношение сигнал/шум, дБ

Без помехи Fп1 Fп2 Fп3

Рис. 5.4. Результаты эксперимента.

Спектральные характеристики и сигнальные созвездия, полученные с

помощью измерительного приложения VMA, приведены на рис. 5.5.
128

а)

б)
 129

в)

г)
Рис. 5.5. Спектральные характеристики и сигнальные созвездия сигнала 16-КАМ:
а – отношение сигнал/шум 25 дБ, без помехи; б – отношение сигнал/шум 25
дБ, помеха Fп3 с уровнем 0 дБ; в – отношение сигнал/шум 25 дБ, помеха Fп2 с
уровнем -20 дБ; г – отношение сигнал/шум 25 дБ, помеха Fп1 с уровнем -20 дБ.
 130

5.3. Выводы по разделу 5

Проведенные экспериментальные испытания позволяют сделать

следующие выводы:
1. По п. 5.1: отношение сигнал/шум (на входе модема) в сети LTE на
маршруте следования железнодорожного состава Москва – Санкт
Петербург изменяется в диапазоне 0 – 27 дБ. С учетом этого диапазона
отношений сигнал/шум в теоретической части работы проводилась оценка
помехоустойчивости приемников сигналов КАМ.
2. По п. 5.2: прием сигнала КАМ становится невозможным при большом
уровне прицельной гармонической помехи (при относительном уровне
помехи более -17 дБ); при увеличении расстройки гармонической помехи
вероятность битовой ошибки снижается на 1-2 порядка. Это подтверждает
выводы, сделанные в п. 2.2 теоретической части работы.
 131

ЗАКЛЮЧЕНИЕ

В соответствии с поставленной целью в результате теоретических и

экспериментальных исследований получены следующие основные
результаты:
1. Определены наиболее типичные виды нефлуктуационных помех,
которые могут присутствовать в радиоканалах систем широкополосного
доступа для мобильных абонентов. К ним относятся гармоническая,
сканирующая, фазоманипулированная и ретранслированная помехи.
2. С использованием методов статистической радиотехники проведен
анализ влияния нефлуктуационных помех на помехоустойчивость
когерентного приемника сигналов с квадратурной амплитудной модуляцией.
3. Показано, что наиболее опасной с точки зрения снижения
помехоустойчивости приема сигналов с квадратурной амплитудной
модуляцией является прицельная гармоническая помеха. Далее опасность
нефлуктуационных помех убывает следующим образом:
фазоманипулированная, ретранслированная, сканирующая. Так, для М=32
при интенсивности гармонической помехи μ = 0,1 энергетические потери по
сравнению со случаем отсутствия помехи составляют около 2 дБ (при
вероятности ошибки 10-2). При воздействии других помех проигрыши
составляют следующие значения: фазоманипулированная – 1,2 дБ;
ретранслированная – 0,6 дБ, сканирующая – 0,5 дБ. Для М=4 (при
вероятности ошибки 10-6) проигрыши меньшие и составляют,
соответственно, 1,2 дБ; 0,65 дБ; 0,5 дБ; 0,4 дБ.
При интенсивности помех μ>0,3 качественный прием сигналов
становится невозможным без использования специальных мер по борьбе с
этими помехами. В связи с этим при разработке методов и устройств приема
сигналов КАМ в сложной помеховой обстановке необходимо предусмотреть
синтез алгоритмов подавления таких помех. При этом в первую очередь
необходимо обратить внимание на синтез алгоритмов, направленных на
 132

подавление узкополосных (гармонических) помех, обладающих наибольшим

мешающим эффектом.
4. Проведено сравнение помехоустойчивости когерентных
приемников сигналов КАМ и сигналов с многопозиционной фазовой
манипуляцией (М-ФМ) одинаковой позиционности и с одинаковой средней
энергией в присутствии рассмотренных нефлуктуационных помех, которое
позволяет утверждать, что в радиоканалах со сложной помеховой
обстановкой предпочтение следует отдать сигналам КАМ. Энергетический
выигрыш при приеме по сравнению со случаем М-ФМ во всех случаях
составляет более 3 дБ при одинаковой позиционности сигналов.
5. Проведен анализ влияния фазовых погрешностей при формировании
опорных колебаний и погрешностей тактовой синхронизации при
когерентном приеме сигналов КАМ, который позволил сделать следующий
вывод:
- при допустимой величине энергетических потерь при приеме, равной
0,5 дБ за счет каждой из этих погрешностей, максимальная фазовая
погрешность должна быть от ~30 при М = 4 до ~10 при М = 64, а допустимая
погрешность тактовой синхронизации – от ~5% при М = 4 до ~2% при М =
64; для обеспечения более жестких требований к величине потерь и величине
Peb требования к указанным погрешностям существенно возрастают.
6. Проведен анализ влияния амплитудного и фазового разбаланса
квадратур при когерентном приеме сигналов КАМ, который позволил
сделать следующие выводы:
- амплитудный разбаланс квадратур при М≥16 ведет в существенному
снижению помехоустойчивости приема сигналов КАМ. Допустимым
отклонением амплитуды в этом случае можно считать величину 5%;
- фазовый разбаланс квадратур сильно влияет на помехоустойчивость
когерентного приема сигналов КАМ. С увеличением позиционности
сигналов это влияние также усиливается. Если для 4-КАМ и 16-КАМ можно
считать допустимой фазовую погрешность в 0,05 рад (30), то 64-КАМ эта
 133

погрешность уже приводит в существенным энергетическим потерям (около

1 дБ). При большей фазовой погрешности энергетические проигрыши
многократно возрастают, особенно при больших М.
7. Синтезирован оптимальный приемник сигналов КАМ на фоне
гармонической помехи со случайной начальной фазой с некогерентной
обработкой помехи. Моделированием установлено, что при μ=0,5 для 4-КАМ
включение схемы некогерентной обработки помехи дает энергетический
выигрыш около 2 дБ, для 16-КАМ – более 5 дБ. Наибольший выигрыш в
помехоустойчивости наблюдается при расстройке ∆fпTs = 0,5.
8. Предложен алгоритм адаптивной фильтрация нефлуктуационных
помех, следящий за фазовой траекторией сигнала КАМ. Моделированием
установлено, что энергетический выигрыш при µ = 0,5 для 4-КАМ доходит
до 4 дБ, для 16-КАМ – до 0,5 дБ.
9. Разработано программное обеспечение для моделирования канала
передачи информации с использованием сигналов КАМ.
10. Проведены натурные аппаратурные исследования
помехоустойчивости приема сигналов 16-КАМ в присутствии гармонической
помехи, результаты которых подтвердили теоретические выводы.

Результаты работы использованы в АО «МНИИРС», ООО «ВизКом» и

могут быть использованы на других предприятиях радиоэлектронной
отрасли, занимающихся разработкой систем передачи дискретных
сообщений.
Направлением дальнейших исследований является совершенствование
методов и алгоритмов для борьбы с нефлуктуационными помехами.
 134

СПИСОК СОКРАЩЕНИЙ

АФ – адаптивный фильтр
АЧХ – амплитудно-частотная характеристика
ВК – весовой коэффициент
ЗССС – земная станция спутниковой связи
КАМ – квадратурная амплитудная манипуляция
КЦС – комплекс цифровой связи
М-ФМ – многопозиционная фазовая манипуляция
ОФМ – относительная фазовая манипуляция
ОСШ – отношение сигнал/шум
ЧМ – частотная манипуляция
ФАПЧ – фазовая автоподстройка частоты
ФМ – фазовая манипуляция
ЭИИМ – эквивалентная изотропная излучаемая мощность
АPSK – amplitude phase shift keying
ATG – air-to-ground
BER – bit error rate
BPSK – binary phase-shift keying
CDMA – code division multiple access
DVB – digital video broadcasting
GSM – global system for mobile communications
LTE – long-term evolution
M-PSK – multiple phase-shift keying
OFDMA – orthogonal frequency-division multiple access
PSK – phase shift-keying
SNR – signal-to-noise ratio
QAM – quadrature amplitude modulation
QPSK – quadrature phase-shift keying
 135

СПИСОК ЛИТЕРАТУРЫ
1. Прокис Дж. Цифровая связь. Пер. с англ. / Под ред. Д.Д. Кловского. –
М.: Радио и связь. 2000. – 800 с.
2. Боккер П. Передача данных. Том 1. (Техника связи в системах
телеобработки данных). Том I. Основы: Пер. с нем./Под ред. Д. Д. Кловского.
– М.: Связь, 1980. – 264 с.
3. Fuqin Xiong. Digital modulation techniques. Second edition. Artech House
Telecommunications Library. Artech House. 2006. – 1039 p.
4. Официальный сайт компании Inmarsat, PLC, URL:
http://www.inmarsat.com (дата обращения: 12.12.2020).
5. Официальный сайт компании Intelsat, URL: http://www.intelsat.com
(дата обращения: 12.12.2020).
6. Официальный сайт компании Eutelsat, URL: http://www.eutelsat.com
(дата обращения: 12.12.2020).
7. Swinford R., Grau B. High Throughput Satellites. Delivering future
capacity needs. – 2015. [Электронный ресурс] URL:
https://www.adlittle.com/sites/default/files/viewpoints/ADL_High_Throughput_Sa
tellites-Viewpoint.pdf (дата обращения: 12.12.2020).
8. Нестеров А.В., Гавриков Н.С., Лелюх А.А. Структура построения
комплексов цифровой связи для летательных аппаратов// Информационно-
измерительные и управляющие системы. 2018. №7. С. 32-37.
9. Булыгина Ю.В., Нестеров А.В. Разработка комплекса цифровой связи
для обеспечения широкополосного доступа к мультимедийным услугам
пассажиров на подвижных объектах// Успехи современной
радиоэлектроники. 2017. №1. С. 17-22.
10. Digital Video Broadcasting (DVB); Implementation guidelines for the
second generation system for Broadcasting, Interactive Services, News Gathering
and other broadband satellite applications; Part 2 - S2 Extensions (DVB-S2X). -
Digital Video Broadcasting (DVB): DVB Document A171-2, March 2015, 183 p.
 136

URL: https://dvb.org/?standard=dvb-s2x-implementation-guidelines (дата

обращения: 12.12.2020).
11. DVB. [Электронный ресурс] URL: https://www.dvb.org/standards/dvb-
s2x (дата обращения: 12.12.2020).
12. Minoli D. Innovations in satellite communications and satellite technology
the industry implications of DVB-S2X, high throughput satellites, Ultra HD,
M2M, and IP. – USA: WILEY, 2015. - 441 p.
13. High Throughput Ku-band for Aero Applications / Ch. McLain, S. Panthi,
J. Hetrick. – 30-th International Communications Satellite Systems Conference
(ICSSC) 24–27 September 2012, Ottawa, Canada.
14. Бакулин М.Г., Крейнделин В.Б., Шлома А.М., Шумов
А.П. Технология OFDM. Учебное пособие для вузов. — М.: Горячая линия -
Телеком, 2015. — 360 с.
15. Samuel C. Yang. OFDMA System Analysis and Design (Mobile
Communications)/ – Publ.: Norwood, MA: Artech House, 2010 – 300 p.
16. Куликов Г.В. Методы помехоустойчивого приема модулированных
сигналов с непрерывной фазой в каналах связи с нефлуктуационными
помехами: дис. … докт. тех. наук: 05.12.04 / Куликов Геннадий
Валентинович; МИРЭА. – М., 2003. – 356 с.
17. Алехин В.А., Шеболков В.В. Помехоустойчивость сигналов с
Квадратурной амплитудной модуляцией//Известия ЮФУ. Технические
науки. 2009. № 1 (90). С. 7-14.
18. Мартиросов В.Е., Рамирес Агилар Х.А. Помехоустойчивость
когерентного приема сигналов QAM//Электросвязь. 2007. № 5. С. 44-48.
19. Выболдин Ю.К. Сравнительная оценка помехоустойчивости приема
сигналов с КАМ// В сборнике: Научная сессия ГУАП. сборник докладов.
Санкт-Петербургский государственный университет аэрокосмического
приборостроения. 2016. С. 26-31.
20. Бабанин И.Г., Мухин И.Е. Оценка влияния неравномерности
амплитудно-частотной характеристики полосового фильтра на
 137

эквивалентные энергетические потери в системах связи, использующих

сигналы с квадратурной амплитудной модуляцией//Телекоммуникации.
2012. № 11. С. 19-23.
21. Бабанин И.Г., Михайлова Н.Ю., Николаенко А.И., Коптев Д.С.
Процедура оценки влияния неидеальности параметров оптимального
демодулирующего устройства системы приема сигналов с квадратурной
амплитудной модуляцией различной позиционности и коэффициента
неравномерности сигнального созвездия на эквивалентные энергетические
потери//Радиотехника. 2019. Т. 83. № 10 (16). С. 60-69.
22. Яценко С.Ю. Анализ влияния искажений и шумов на
помехоустойчивость приема сигналов квадратурной амплитудной модуляции
//Журнал радиоэлектроники. 2016. № 5. С. 3.
23. Савищенкo Н.В. Помехоустойчивость когерентного приема
многопозионных сигналов КАМ и ФМ при неидеальной синхронизации//
Информационно-управляющие системы. 2010. № 1 (44). С. 53-63.
24. Артеменко А.А., Мальцев А.А., Рубцов А.Е. Влияние неточности
оценивания фазы несущей на вероятность битовых ошибок в М-КАМ
системах передачи данных//Вестник Нижегородского университета им. Н.И.
Лобачевского. 2007. № 2. С. 81-87.
25. Довбня В.Г., Севрюков А.Е. Влияние качества функционирования
системы тактовой синхронизации на помехоустойчивость приема КАМ-
сигналов//Телекоммуникации. 2016. № 4. С. 43-47.
26. Довбня В.Г., Севрюков А.Е. Влияние качества функционирования
системы тактовой синхронизации на помехоустойчивость приема КАМ-
сигналов//Проектирование и технология электронных средств. 2015. № 4. С.
24-28.
27. Залевский А.А., Рогожников Е.В. Влияние временной
синхронизации на вероятность битовой ошибки в системе связи с QAM-
модуляцией//Электронные средства и системы управления. Материалы
 138

докладов Международной научно-практической конференции. 2016. № 1-1.

С. 51-53.
28. Saeedi-Sourck H., Sadri S., Wu Y., Bergmans J.W.M., Farhang-
Boroujeny B. Sensitivity analysis of offset QAM multicarrier systems to residual
carrier frequency and timing offsets//Signal Processing. 2011. Т. 91. № 7. С. 1604-
1612.
29. Гуревич В.Э., Девяткин Д.В., Харьков О.И. Исследование реальной
помехоустойчивости системы радиодоступа с квадратурной модуляцией при
дефазировании опорных колебаний//В сборнике: Актуальные проблемы
инфотелекоммуникаций в науке и образовании. Сборник научных статей: в 4
томах. 2017. С. 154-159.
30. Довбня В.Г., Коптев Д.С. Влияние качества функционирования
гетеродинов на помехоустойчивость приема сигналов с квадратурной
амплитудной модуляцией //Радиотехника. 2020. Т. 84. № 9 (17). С. 40-48.
31. Хасьянова Е.Р., Седов М.О. Экспериментальное исследование
методов оценки и компенсации погрешностей квадратурного преобразования
сигналов с модуляцией М-КАМ//Системы синхронизации, формирования и
обработки сигналов. 2017. Т. 8. № 1. С. 119-122.
32. Пестряков А.В., Поборчая Н.Е., Хасьянова Е.Р. Оценка влияния
сдвига несущей частоты на эффективность работы алгоритма компенсации
погрешностей квадратурного преобразования//В сборнике: 26-я
Международная Крымская конференция «СВЧ-техника и
телекоммуникационные технологии» (КрыМиКо'2016). материалы
конференции: в 13 томах. 2016. С. 413-417.
33. Пестряков А.В., Хасьянова Е.Р. Оценка и компенсация
погрешностей квадратурного преобразования при приеме сигнала с
модуляцией КАМ-64 с применением метода статистического усреднения //
В сборнике: Технологии информационного общества. X Международная
отраслевая научно-техническая конференция: сборник трудов. 2016. С. 196.
 139

34. Park J., Park S., Cho P.D., Lee H. Ho. Effect of phase/quadrature error
and 1/Q gain imbalance for QAM symbol error probability//
IEICE Transactions on Communications. 2007. Т. E90-B. № 11. С. 3287.
35. Мартиросов В.Е. Влияние нелинейных (компрессионных)
искажений на помехоустойчивость приема сигналов QAM//Радиотехника.
2008. № 9. С. 4-11.
36. Башкиров А.В., Крисилов А.В., Машин В.В., Муратов А.В.,
Хорошайлова М.В. Анализ вероятностей символьных ошибок для 16-QAM-
сигналов в присутствии QPSK-сигналов и 16-QAM-сигналов//Радиотехника.
2020. Т. 84. № 6 (12). С. 5-9
37. Поборчая Н.Е., Смердова Е.О. Вариационный алгоритм
компенсации искажений КАМ сигнала на фоне квазидетерминированной
полосовой помехи и аддитивного белого шума//
Системы синхронизации, формирования и обработки сигналов. 2014. Т.
5. № 4. С. 141-147.
38. Ложкин К.Ю., Петров А.В., Прожеторко С.С. Аналитические
зависимости средней вероятности искажения бита M-КАМ сигнала на фоне
гармонической или фазоманипулированной помех//
Электромагнитные волны и электронные системы. 2018. Т. 23. № 5. С. 32-41.
39. Ложкин К.Ю., Миронов В.А., Прожеторко С.С., Токарева Н.Ю.
Эффективность приема М-КАМ-сигнала в условиях воздействия импульсной
гармонической помехи// Радиотехника. 2018. № 11. С. 52-57.
40. Ложкин К.Ю., Петров А.В., Миронов В.А., Михалёв В.В.,
Прожеторко С.С. Аналитические зависимости средней вероятности
искажения бита М-КАМ-сигнала на фоне гармонической или
фазоманипулированной помех с учетом замираний//Радиотехника. 2020. Т.
84. № 4 (8). С. 27-35.
41. Геранин Р.В., Киселёв В.В., Львов А.А., Руденко А.Е., Светлов М.С.
Классификация и оценивание помех при использовании QAM
 140

сигналов//Вестник Саратовского государственного технического

университета. 2010. Т. 4. № 2 (50). С. 10-12.
42. Сидельников Г.М. Сравнительный анализ эффективности
разнесенного приема сигналов с квадратурной амплитудной и с фазовой
модуляцией в канале с дискретной многолучёвостью //Вестник Поволжского
государственного технологического университета. Серия: Радиотехнические
и инфокоммуникационные системы . 2020. № 2 (46). С. 18-30.
43. Выболдин Ю.К., Зайченко К.В. Алгоритм приема сигналов с
квадратурной амплитудной модуляцией в каналах связи с замираниями//
Радиотехника. 2018. № 9. С. 193-199.
44. Zakharov Y.V., Zhang J., Baronkin V.M. Optimal and mismatched detection
of QAM signals in fast fading channels with imperfect channel estimation//IEEE
Transactions on Wireless Communications. 2009. Т. 8. № 2. С. 617-621.
45. Поборчая Н.Е. Анализ алгоритмов компенсации искажений сигнала
КАМ-64, наблюдаемого на фоне аддитивного белого шума и МСИ//
Системы синхронизации, формирования и обработки сигналов. 2015. Т.
6. № 4. С. 173-175.
46. Довбня В.Г., Вольский В.А. Помехоустойчивость приема
многопозиционных сигналов с квадратурной амплитудной манипуляцией в
канале с межсимвольной интерференцией//Известия Юго-Западного
государственного университета. Серия: Управление, вычислительная
техника, информатика. Медицинское приборостроение. 2013. № 2. С. 17-20.
47. Kulikov G., Nesterov A., Leluh A. Simulation of the Radio
Communication Channel with Unmanned Aerial Vehicle in Urban Conditions//
Journal of Electrical and Electronic Engineering, Vol. 6, Is. 4, August 2018, pp.
104-110, doi: 10.11648/j.jeee.20180604.11.
48. Нестеров А.В., Лелюх А.А., Гавриков Н.С., Куликов Г.В. Оценка
энергетических характеристик радиолиний перспективных систем
спутникового широкополосного доступа для аэромобильных абонентов//
 141

Журнал радиоэлектроники [электронный журнал]. 2019. №1. Режим доступа:

http://jre.cplire.ru/jre/jan19/11/text.pdf. DOI 10.30898/1684-1719.2019.1.11.
49. Куликов Г.В., Нестеров А.В., Лелюх А.А., Нгуен Ван Зунг. О
помехоустойчивости приема сигналов с многопозиционной манипуляцией в
перспективных системах спутникового широкополосного доступа для
аэромобильных абонентов//В сб. материалов III международной заочной
научно-практической конференции «Авиация: История, Современность,
Перспективы развития», Минск, БГАА, 8–9 ноября 2018 г., с. 84-86.
50. Нестеров А.В., Лелюх А.А., Гавриков Н.С., Куликов Г.В. Основные
тенденции и перспективы развития спутникового широкополосного доступа
для аэромобильных абонентов // В сб. материалов III международной заочной
научно-практической конференции «Авиация: История, Современность,
Перспективы развития», Минск, БГАА, 8–9 ноября 2018 г., с. 97-98.
51. Куликов Г.В., Лелюх А.А. Влияние структурных помех на качество
связи в перспективных системах спутникового широкополосного доступа для
аэромобильных абонентов// Актуальные проблемы и перспективы развития
авиации: сборник материалов III Международной научно-практической
конференции БГАА, Минск, 22–23 мая 2019 г. / ред. А. А. Жукова [и др.]; под
научн. ред. А. А. Шегидевича. – Минск: БГАА, 2019.
52. Тихонов В.И. Статистическая радиотехника. - 2-е изд., перераб. и
доп. - М. : Радио и связь, 1982. - 624 с. : ил.
53. Тихонов В.И. Оптимальный прием сигналов. М.: Радио и связь, 1983.
54. Куликов Г.В., Нгуен Ван Зунг, Нестеров А.В., ЛелюхА.А.
Помехоустойчивость приема сигналов с многопозиционной фазовой
манипуляцией в присутствии гармонической помехи // Наукоемкие
технологии. 2018. №11. С. 32–38.
55. Куликов Г.В., Нгуен Ван Зунг. Анализ помехоустойчивости приема
сигналов с многопозиционной фазовой манипуляцией при воздействии
сканирующей помехи // Российский технологический журнал, 2018, Т. 6, №
6, с. 5–12. DOI: 10.32362/2500-316X-2018-6-6-5-12.
 142

56. Куликов Г.В., Нгуен Ван Зунг, До Чунг Тиен. Влияние

фазоманипулированной помехи на помехоустойчивость корреляционного
демодулятора сигналов с многопозиционной фазовой манипуляцией //
Российский технологический журнал [электронный журнал]. 2019. Т. 7. № 2.
С. 18–28. Режим доступа: https://rtj.mirea.ru/upload/medialibrary/4c2/RTZH_
2_2019_18_28.pdf.
57. Нгуен Ван Зунг. Помехоустойчивость корреляционного приемника
сигналов с многопозиционной фазовой манипуляцией при наличии
ретранслированной помехи // Журнал радиоэлектроники [электронный
журнал]. 2019. №3. Режим доступа: http://jre.cplire.ru/jre/mar19/4/text.pdf. DOI
10.30898/1684-1719.2019.3.4.
58. Куликов Г.В., Нестеров А.В., Лелюх А.А. Помехоустойчивость
приема сигналов с квадратурной амплитудной манипуляцией в присутствии
гармонической помехи // Журнал радиоэлектроники [электронный журнал].
2018. № 11. Режим доступа: http://jre.cplire.ru/jre/nov18/9/text.pdf. DOI
10.30898/1684-1719.2018.11.9.
59. Куликов Г.В., Лелюх А.А. Анализ помехоустойчивости приема
сигналов с квадратурной амплитудной манипуляцией при воздействии
сканирующей помехи // Журнал радиоэлектроники [электронный журнал].
2019. № 2. Режим доступа: http://jre.cplire.ru/jre/feb19/5/text.pdf.
60. Куликов Г.В., Лелюх А.А., Баталов Е.В., Кузеленков П.И.
Помехоустойчивость приема сигналов с квадратурной амплитудной
манипуляцией в присутствии фазоманипулированной помехи // Журнал
радиоэлектроники [электронный журнал]. 2019. № 7. Режим доступа:
http://jre.cplire.ru/jre/jul19/10/text.pdf DOI 10.30898/1684-1719.2019.7.10
61. Куликов Г.В., Лелюх А.А., Баталов Е.В., Граченко Е.Н.
Помехоустойчивость приема сигналов с квадратурной амплитудной
манипуляцией в канале связи с нефлуктуационными помехами// В кн.: Сб.
науч. трудов IV Международной НПК "Радиоинфоком-2019", с. 100-103, – М.:
МИРЭА, ноябрь 2019 г.
 143

62. Куликов Г.В., Лелюх А.А., Граченко Е.Н. Помехоустойчивость

когерентного приемника сигналов с квадратурной амплитудной
манипуляцией при наличии ретранслированной помехи// Радиотехника и
электроника, 2020, т. 65, №8, с. 804-808. DOI: 10.31857/S0033849420070074.
63. Нгуен Ван Зунг. Помехоустойчивость приема дискретных сигналов
с многопозиционной фазовой манипуляцией при наличии нефлуктуационных
помехам: дис. … канд. тех. наук: 05.12.04 / Нгуен Ван Зунг; РТУ МИРЭА. –
М. 2020. – 179 с.
64. Куликов Г.В., Ван Зунг Нгуен. Влияние погрешностей
синхронизации на помехоустойчивость когерентного приема сигналов М-
ФМ// Российский технологический журнал [электронный журнал].
2019;7(5):47-61. https://doi.org/10.32362/2500-316X-2019-7-5-47-61.
65. Куликов Г.В., Ван Зунг Нгуен. Влияние неточности тактовой
синхронизации на помехоустойчивость когерентного приема сигналов М-
ФМ// В кн.: Сб. науч. трудов IV Международной НПК "Радиоинфоком-2019",
с. 109-113, – М.: МИРЭА, ноябрь 2019 г.
66. Куликов Г.В., Ван Зунг Нгуен. Влияние неточности оценивания
частоты и фазы несущей на помехоустойчивость когерентного приема
сигналов М-ФМ// В кн.: Сб. науч. трудов IV Международной НПК
"Радиоинфоком-2019", с. 114-122, – М.: МИРЭА, ноябрь 2019 г.
67. Kulikov G.V., Nguen Van Dung, Kulagin V.P., Lelyukh A.A. Influence
of synchronization errors on the noise immunity of receiving multi-position M-
PSK and M-QAM signals// 2020 Systems of Signal Synchronization, Generating
and Processing in Telecommunications (SYNCHROINFO), Svetlogorsk, Russia,
2020, pp. 1-5, doi: 10.1109/SYNCHROINFO49631.2020.9166034.
68. Бывшев М.Е. и др. Оптимальный прием сигналов на фоне помех и
шумов. / Под ред. Ю.И. Савватеева. М.: Радиотехника, 2011.
69. Курицын С.А. Методы адаптивной обработки сигналов передачи
данных. – М.: Радио и связь, 1988. – 144 с.: ил.
 144

70. Уидроу Б., Стирнз С. Д. Адаптивная обработка сигналов. Пер. с

англ. под ред. Шахгильдяна В. В. – М.: Радио и связь, 1989. – 440 с.
71. Джиган В.И. Адаптивная фильтрация сигналов: теория и
алгоритмы. – М.: Техносфера, 2013. – 530 с.
72. Куликов Г.В., Кулагин В.П., До Чунг Тиен, Лелюх А.А.
Адаптивный цифровой фильтр для подавления нефлуктуационных помех.
Патент на полезную модель RU 198305 U1. Дата регистрации: 30.06.2020.
Номер заявки: 2020108179. Опубликовано: 30.06.2020 Бюл. № 19
73. Treichler J.R., Agee B.G. A new approach to multipath correction of
constant modulus signals. IEEE Trans. 1983. V. ASSP-31. № 2. P. 459–472.
74. Куликов Г.В., До Чунг Тиен. Эффективность адаптивного фильтра с
алгоритмом слежения за огибающей при приеме сигналов с
многопозиционной ФМ на фоне нефлуктуационных помех// Российский
технологический журнал. 2020;8(5):34-43. https://doi.org/10.32362/2500-316X-
2020-8-5-34-43.
75. Куликов Г.В. Методы помехоустойчивого приема модулированных
сигналов с непрерывной фазой в каналах связи с нефлуктуационными
помехами: дис. … докт. тех. наук: 05.12.04 / Куликов Геннадий
Валентинович; МИРЭА. – М., 2003. – 356 с.
76. Куликов Г.В., До Чунг Тиен. Эффективность фазового алгоритма
адаптивной фильтрации при приеме сигналов с многопозиционной фазовой
манипуляцией// Журнал радиоэлектроники. 2020. № 4. Режим доступа:
http://jre.cplire.ru/jre/apr20/9/text.pdf. DOI 10.30898/1684-1719.2020.4.9.
77. Куликов Г.В., Лелюх А.А. Влияние амплитудного и фазового
разбаланса квадратур на помехоустойчивость когерентного приема сигналов
с квадратурной амплитудной модуляцией// Российский технологический
журнал. 2021;9(1):29-37.
78. Куликов Г.В., Лелюх А.А., Баталов Е.В. Влияние погрешностей
фазовой и тактовой синхронизации на помехоустойчивость когерентного
 145

приема сигналов с квадратурной амплитудной модуляцией// Российский

технологический журнал. 2021;9(2):35-43.
79. Куликов Г.В., Лелюх А.А., Граченко Е.Н. Адаптивная фильтрация
гармонической помехи при приеме сигналов с квадратурной амплитудной
модуляцией// Журнал радиоэлектроники [электронный журнал]. 2021. №8.
80. Куликов Г.В., Лелюх А.А., Нгуен Ван Зунг, Кулагин В.П.
Помехоустойчивость приема многопозиционных сигналов М-ФМ и М-КАМ
в канале связи с нефлуктуационными помехами// Труды международного
симпозиума «Надежность и качество», 2020, Пенза, т. 2, с. 333-336.
 146

Приложение 1. Алгоритм и программа моделирования радиоканала с

использованием сигнала 16-КАМ (когерентный прием)

П1.1. Описание алгоритма работы программы моделирования

радиоканала с использованием сигнала КАМ

Алгоритм работы программы (рис. П.1) содержит следующие

процедуры:

Ввод исходных данных. Задание основных параметров используемых

моделей сигналов и помех. Задание значения рассинхронизации фазы.

1. Модулятор КАМ сигнала: формирование информационной

последовательности (генератор случайных чисел), формирование КАМ
сигнала квадратурным методом.

2. Канал связи: формирование шумовой, гармонической и

ретранслированной помех с заданными параметрами, формирование смеси
сигнала и помех.

3. Демодулятор КАМ сигнала: обработка смеси сигнала и помех в

квадратурном демодуляторе, формирование принятой последовательности
данных.

4. Вычисление вероятности ошибки: сравнение переданного и принятого

массивов данных, вычисление вероятности ошибки.
 147

Начало

Ввод исходных данных

Формирование информационной последовательности I и Q

Формирование сигнала синфазного Формирование сигнала

канала квадратурного канала
ichanel=I.*cos(2*pi*fc*t+dfi) qchanel=-Q.*sin(2*pi*fc*t+dfi)

Формирование КАМ сигнала

S = ichanel + qchanel

Формирование шумовой помехи

n(t)

Формирование гармонической и ретранслированной помехи

S2(m1, t), S22(m2, t)

Формирование сигнала на выходе канала связи

x(t) = S(t) + n(t) + S2(m1, t) + SS2(m2, t)

Вычисление корреляционных интегралов

min(I)=|I-X| min(Q)=|Q-Y|

error_bit = ctr1/N;
4. Вычисление
вероятности
ошибки

Конец

Рис. П.1. Структура алгоритма моделирования.

 148

П1.2. Программа моделирования радиоканала с использованием

сигналов 16-КАМ

Данная программа моделирует информационный радиоканал, в

котором задается полезный сигнал 16-КАМ, гауссовский шум,
гармоническая и ретранслированная помехи с заданными параметрами, а
также предусмотрена возможность задания фазы рассинхронизации каналов.
Осуществляется обработка смеси сигналов и помех в демодуляторе,
сравнение переданного и принятого сигналов и вычисление вероятности
битовой ошибки. Значения битовой ошибки записываются в текстовый файл
в зависимости от отношения сигнал/шум при фиксированном значении
рассинхронизации фазы каналов.

Последовательность операций:

1. Запуск приложения Matlab или Scilab.

2. Создание нового файла, подготовка памяти (строка программы 1).

3. Ввод исходных данных (строки программы 2-15):

- создание вспомогательного текстового файла data.txt;

- амплитуда сигнала 16-КАМ A0;

- длина информационной последовательности N;

- частота несущего колебания fc;

- частота информационных символов fz;

- период дискретизации td;

- количество отсчетов на бит M;

- амплитудные значения I для сигнальных точек 16-КАМ Iamp;

- амплитудные значения Q для сигнальных точек 16-КАМ Qamp;

- фаза рассинхронизации dfi;

- интенсивность гармонической помехи m1;

 149

- интенсивность ретранслированной помехи m2;

- счетчик для набора вектора ошибок ce;

- время задержки ретранслированного сигнала относительно длительности

канального символа p.

4. Формирование 16-КАМ сигнала (строки программы 16-76):

- формирование случайной информационной последовательности b;

- формирование информационной последовательности I и Q каналов;

- формирование сигнала синфазного канала: ichanel=I.*cos(2*pi*fc*t+dfi);

- формирование сигнала квадратурного канала: qchanel=-

Q.*sin(2*pi*fc*t+dfi);

- формирование 16-КАМ сигнала: S = ichanel + qchanel;

5. Модель канала связи (строки программы 77-92):

- формирование гауссовского шума noise;

- формирование гармонической помехи S2;

- формирование ретранслированной помехи S22;

- формирование на выходе суммы сигнала 16-КАМ + гармоническая помеха

+ ретранслированная помеха + гауссовский шум S1 = S + S2 + S22 +noise.

6. Демодуляция 16-КАМ сигнала (строки программы 93-206):

- вычисление корреляционных интегралов синфазного канала SS2M и

квадратурного канала SS3M;

- определение ближайшего значения интеграла по оси I – min(abs(Iamp -

SS2M(n)));

- определение ближайшего значения интеграла по оси Q – min(abs(Qamp -

SS3M(n)));

- принятие решения об информационном символе d;

 150

7. Вычисление вероятности битовой ошибки (строки программы 207-219):

- вычисление вероятности битовой ошибки error_bit;

- запись в текстовый файл массива ошибок fprintf(fileID,"%4f\n",verbit(ce)).

Листинг программы
clc;clear all; close all;
% Ввод исходных данных
fileID = fopen('data.txt','w'); %создание вспомогательного
текстового файла data.txt
A0 = 1; % амплитуда сигнала
N = 108000; % длина информационной последовательности
fc = 10; % частота несущего колебания
fz = 1; % частота информационных символов
fd = 20*fc; % частота дискретизации
td = 1/fd; % период дискретизации
M = fd/fz; % количество отсчетов на бит
Iamp = [126.4911 379.4733 -126.4911 -379.4733]; % амплитудные
значения I для сигнальных точек 16-КАМ
Qamp = [126.4911 379.4733 -126.4911 -379.4733]; % амплитудные
значения Q для сигнальных точек 16-КАМ
dfi = 0; % фаза рассинхронизации
m1 = 0; % интенсивность гармонической помехи
m2 = 0; % интенсивность ретранслированной помехи
ce = 1; % счетчик для массива ошибок, чтобы строить график
p = 0; % время задержки относительно длительности канального
символа
% Модулятор 16-КАМ сигнала
b = randi([0 1], 1, N);% формирование случайной информационной
последовательности
% Деление на I канал и Q канал с домножением на корректирующие
коэффициенты
for n = 1 : length(b)/4
if b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 0 && b(4*n)
== 1
I(n) = -3/sqrt(10);
Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 1
I(n) = -1/sqrt(10);
Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 1
I(n) = 1/sqrt(10);
Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 0 &&
b(4*n) == 1
I(n) = 3/sqrt(10);
 151

Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = -3/sqrt(10);
Q(n) = 1/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = -1/sqrt(10);
Q(n) = 1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = 1/sqrt(10);
Q(n) = 1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = 3/sqrt(10);
Q(n) = 1/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 1 &&
b(4*n) == 0
I(n) = -3/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 0
I(n) = -1/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 0
I(n) = 1/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 1 &&
b(4*n) == 0
I(n) = 3/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = -3/sqrt(10);
Q(n) = -3/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = -1/sqrt(10);
Q(n) = -3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = 1/sqrt(10);
Q(n) = -3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = 3/sqrt(10);
Q(n) = -3/sqrt(10);
end
end
I = repmat(I, 4*M, 1); % копируем последовательность I по
строкам
 152

Q = repmat(Q, 4*M, 1); % копируем последовательность Q по

строкам
I = I(:)'; % вытягиваем по столбцам
Q = Q(:)'; % вытягиваем по столбцам
t = 0 : td : (N * M - 1) * td; % массив времени
ichanel = I.* cos(2 * pi * fc * t+dfi); % сигнал I - канал
qchanel = - Q .* sin(2 * pi * fc * t+dfi); % сигнал Q - канал
S = ichanel + qchanel; % итоговый сигнал I + Q
% Канал связи
for SNR = 0:13 % Цикличное изменение отношения сигнал/шум
% Формирование помех
% Гауссовский шум
sigma = (A0/2) .* sqrt((fz*fd)/(10 .^ (SNR/10)));
noise = sigma .* randn(1,length(t));
% Гармоническая помеха со случайной фазой
S2 = m1 .* A0 .* sin(2*pi*fc*t + (2*pi/360) *
fix((1:numel(S))/60000));
% Ретранслированная помеха со случайной фазой
SS2 = m2 .* (I .* cos(2 * pi * fc * t + (2*pi/360) * fix((1 :
numel(S))/ 60000)) - Q .* sin(2 * pi * fc * t + (2*pi/360) *
fix((1 : numel(S))/ 60000)));
S22 = [zeros(1, p * M) SS2(1 : end - p * M)];
% Сумма полезного сигнала и помех
S1 = S + S2 + S22 + noise;
% Демодулятор 16-КАМ сигнала (по амплитуде)
% Создание массивов сигналов
S2 = cos(2*pi*fc*t);
S3 = -sin(2*pi*fc*t);
S1m = reshape(S1,4*M,N/4)';
S2m = reshape(S2,4*M,N/4)';
S3m = reshape(S3,4*M,N/4)';
SS2M = zeros(1,N/4);
SS3M = zeros(1,N/4);
index = SS2M;
d = zeros(1,N);
% Вычисление корреляционных интегралов и принятие решения об
информационных символах
for n = 1 : N/4
for ii = 1 : 4*M
SS2M(n) = SS2M(n) + S1m(n,ii).*S2m(n,ii);
SS3M(n) = SS3M(n) + S1m(n,ii).*S3m(n,ii);
end
[~,isw(n)] = min(abs(Iamp - SS2M(n)));
switch isw(n)
case 1
[~,qsw(n)] = min(abs(Qamp - SS3M(n)));
switch qsw(n)
case 1
d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 0;
d(4*n) = 0;
case 2
 153

d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 0;
d(4*n) = 1;
case 3
d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 1;
d(4*n) = 0;
case 4
d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 1;
d(4*n) = 1;
end
case 2
[~,qsw(n)] = min(abs(Qamp - SS3M(n)));
switch qsw(n)
case 1
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 0;
case 2
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 1;
case 3
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 1;
d(4*n) = 0;
case 4
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 1;
d(4*n) = 1;
end
case 3
[~,qsw(n)] = min(abs(Qamp - SS3M(n)));
switch qsw(n)
case 1
d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 0;
d(4*n) = 0;
case 2
d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 0;
d(4*n) = 1;
case 3
 154

d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 1;
d(4*n) = 0;
case 4
d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 1;
d(4*n) = 1;
end
case 4
[~,qsw(n)] = min(abs(Qamp - SS3M(n)));
switch qsw(n)
case 1
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 0;
case 2
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 1;
case 3
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 1;
d(4*n) = 0;
case 4
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 1;
d(4*n) = 1;
end
end
end
% Вычисление вероятности битовой ошибки
ctr1 = 0; % счетчик ошибок
for e = 1 : N
if b(e) - d(e) ~= 0
ctr1 = ctr1 + 1;
end
end
error_bit = ctr1/N;
verbit(ce) = error_bit; % массив из значений вероятностей
битовых ошибок
fprintf(fileID,"%4f\n",verbit(ce)); % запись в текстовый файл
массива ошибок
ce = ce+1;
ctr1 = 0;
end
fclose(fileID);
 155

Приложение 2. Программа моделирования радиоканала с

использованием сигнала 16-КАМ (некогерентная обработка помехи)
clc; clear;
// Ввод исходных данных
A0 = 1;
N = 1080;
fc = 10;
fz = 1;
fd = 20*fc;
td = 1/fd;
M = fd/fz;
SNR = 7; //осш в дБ
m1 = 0.9;
dm1=0.5; //отклонение m1 в приемнике

// Модулятор 16-КАМ сигнала

b = grand(1, N, 'uin', 0, 1);
// деление на I канал и Q канал
for n = 1 : length(b)/4
if b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 0 && b(4*n)
== 1
I(n) = -3/sqrt(10);
Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 1
I(n) = -1/sqrt(10);
Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 1
I(n) = 1/sqrt(10);
Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 0 &&
b(4*n) == 1
I(n) = 3/sqrt(10);
Q(n) = 3/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = -3/sqrt(10);
Q(n) = 1/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = -1/sqrt(10);
Q(n) = 1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = 1/sqrt(10);
Q(n) = 1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 0 &&
b(4*n) == 0
I(n) = 3/sqrt(10);
Q(n) = 1/sqrt(10);
 156

elseif b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 1 &&

b(4*n) == 0
I(n) = -3/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 0
I(n) = -1/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 0
I(n) = 1/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 1 &&
b(4*n) == 0
I(n) = 3/sqrt(10);
Q(n) = -1/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 1 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = -3/sqrt(10);
Q(n) = -3/sqrt(10);
elseif b(4*n-3) == 1 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = -1/sqrt(10);
Q(n) = -3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 0 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = 1/sqrt(10);
Q(n) = -3/sqrt(10);
elseif b(4*n-3) == 0 && b(4*n-2) == 1 && b(4*n-1) == 1 &&
b(4*n) == 1
I(n) = 3/sqrt(10);
Q(n) = -3/sqrt(10);
end
end
I = repmat(I, 4*M, 1);
Q = repmat(Q, 4*M, 1);
I = I(:)';
Q = Q(:)';
t = 0 : td : (N*M-1)*td;
i = fix(t/(td*4*M)) + 1;
ichanel = I(i).* cos(2*%pi*fc*t);
qchanel = - Q(i) .* sin(2*%pi*fc*t);
S1 = ichanel + qchanel;
// Канал связи
// Формирование помех
// Гауссовский шум
sigma = (A0/2) .* sqrt((fz*fd)/(10 .^ (SNR/10)));
noise = sigma .* rand(1, length(t), 'normal');
N0=2;
// Гармоническая помеха со случайной фазой
S2 = m1 .* A0 .* sin(2*%pi*(fc+0.125)*t + (2*%pi/360) .*
fix((1:length(S1))/600));
 157

// Сумма 16KAМ сигнала + гармоническая помеха + гауссовский шум

S3 = S1 + S2 + noise;
// Демодулятор 16-КАМ сигнала
// Формирование опорных сигналов

//квадратурные копии помехи

P1=cos(2*%pi*(fc+0.125)*t);
P2=sin(2*%pi*(fc+0.125)*t);
// Формирование опорных сигналов
S01=(-3/sqrt(10))*cos(2*%pi*fc*t)-(3/sqrt(10))*sin(2*%pi*fc*t);
S02=(-1/sqrt(10))*cos(2*%pi*fc*t)-(3/sqrt(10))*sin(2*%pi*fc*t);
S03=(1/sqrt(10))*cos(2*%pi*fc*t)-(3/sqrt(10))*sin(2*%pi*fc*t);
S04=(3/sqrt(10))*cos(2*%pi*fc*t)-(3/sqrt(10))*sin(2*%pi*fc*t);
S05=(-3/sqrt(10))*cos(2*%pi*fc*t)-(1/sqrt(10))*sin(2*%pi*fc*t);
S06=(-1/sqrt(10))*cos(2*%pi*fc*t)-(1/sqrt(10))*sin(2*%pi*fc*t);
S07=(1/sqrt(10))*cos(2*%pi*fc*t)-(1/sqrt(10))*sin(2*%pi*fc*t);
S08=(3/sqrt(10))*cos(2*%pi*fc*t)-(1/sqrt(10))*sin(2*%pi*fc*t);
S09=(-3/sqrt(10))*cos(2*%pi*fc*t)-(-1/sqrt(10))*sin(2*%pi*fc*t);
S010=(-1/sqrt(10))*cos(2*%pi*fc*t)-(-
1/sqrt(10))*sin(2*%pi*fc*t);
S011=(1/sqrt(10))*cos(2*%pi*fc*t)-(-1/sqrt(10))*sin(2*%pi*fc*t);
S012=(3/sqrt(10))*cos(2*%pi*fc*t)-(-1/sqrt(10))*sin(2*%pi*fc*t);
S013=(-3/sqrt(10))*cos(2*%pi*fc*t)-(-
3/sqrt(10))*sin(2*%pi*fc*t);
S014=(-1/sqrt(10))*cos(2*%pi*fc*t)-(-
3/sqrt(10))*sin(2*%pi*fc*t);
S015=(1/sqrt(10))*cos(2*%pi*fc*t)-(-3/sqrt(10))*sin(2*%pi*fc*t);
S016=(3/sqrt(10))*cos(2*%pi*fc*t)-(-3/sqrt(10))*sin(2*%pi*fc*t);
// Формирование опорных сигналов для смещения всех квадрантов
созвездия в центр
Ssm1=(-2/sqrt(10))*cos(2*%pi*fc*t)+(2/sqrt(10))*sin(2*%pi*fc*t);
Ssm2=(2/sqrt(10))*cos(2*%pi*fc*t)+(2/sqrt(10))*sin(2*%pi*fc*t);
Ssm3=(2/sqrt(10))*cos(2*%pi*fc*t)-(2/sqrt(10))*sin(2*%pi*fc*t);
Ssm4=(-2/sqrt(10))*cos(2*%pi*fc*t)-(2/sqrt(10))*sin(2*%pi*fc*t);
// Формирование опорных сигналов разностей сигнал-опорное
колебание= оценка помехи
S01m=S3-S01;
S02m=S3-S02;
S03m=S3-S03;
S04m=S3-S04;
S05m=S3-S05;
S06m=S3-S06;
S07m=S3-S07;
S08m=S3-S08;
S09m=S3-S09;
S010m=S3-S010;
S011m=S3-S011;
S012m=S3-S012;
S013m=S3-S013;
S014m=S3-S014;
S015m=S3-S015;
S016m=S3-S016;
// Формирование матриц всех сигналов
 158

S31 = matrix(S3, 4*M, N/4)';

S106 = matrix(S06, 4*M, N/4)';

S107 = matrix(S07, 4*M, N/4)';
S1010 = matrix(S010, 4*M, N/4)';
S1011 = matrix(S011, 4*M, N/4)';

S1sm1 = matrix(Ssm1, 4*M, N/4)';

S1sm2 = matrix(Ssm2, 4*M, N/4)';
S1sm3 = matrix(Ssm3, 4*M, N/4)';
S1sm4 = matrix(Ssm4, 4*M, N/4)';

S101m = matrix(S01m, 4*M, N/4)';

S102m = matrix(S02m, 4*M, N/4)';
S103m = matrix(S03m, 4*M, N/4)';
S104m = matrix(S04m, 4*M, N/4)';
S105m = matrix(S05m, 4*M, N/4)';
S106m = matrix(S06m, 4*M, N/4)';
S107m = matrix(S07m, 4*M, N/4)';
S108m = matrix(S08m, 4*M, N/4)';
S109m = matrix(S09m, 4*M, N/4)';
S1010m = matrix(S010m, 4*M, N/4)';
S1011m = matrix(S011m, 4*M, N/4)';
S1012m = matrix(S012m, 4*M, N/4)';
S1013m = matrix(S013m, 4*M, N/4)';
S1014m = matrix(S014m, 4*M, N/4)';
S1015m = matrix(S015m, 4*M, N/4)';
S1016m = matrix(S016m, 4*M, N/4)';

P11 = matrix(P1, 4*M, N/4)';

P22 = matrix(P2, 4*M, N/4)';
// обнуление массивов для интеграторов
SS44=zeros(1,N/4);
U6=SS44; SS06i=SS44;
U7=SS44; SS07i=SS44;
U10=SS44; SS010i=SS44;
U11=SS44; SS011i=SS44;

X6=SS44; Y6=SS44;
X7=SS44; Y7=SS44;
X10=SS44; Y10=SS44;
X11=SS44; Y11=SS44;

tong6=zeros(1,N/4);
tong7=zeros(1,N/4);
tong10=zeros(1,N/4);
tong11=zeros(1,N/4);

d=zeros(1,N);

for n=1:N/4
SS07i(n)=0; X7(n)=0; Y7(n)=0;
SS06i(n)=0; X6(n)=0; Y6(n)=0;
 159

SS010i(n)=0; X10(n)=0; Y10(n)=0;

SS011i(n)=0; X11(n)=0; Y11(n)=0;

for ii=1:4*M
//вычисление основных интегралов
SS07i(n)=SS07i(n)+(2/N0)*(S31(n,ii).*S107(n,ii));
SS06i(n)=SS06i(n)+(2/N0)*(S31(n,ii).*S106(n,ii));
SS010i(n)=SS010i(n)+(2/N0)*(S31(n,ii).*S1010(n,ii));
SS011i(n)=SS011i(n)+(2/N0)*(S31(n,ii).*S1011(n,ii));
//вычисление корректирующих слагаемых в алгоритме
X7(n)=X7(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S107m(n,ii));
Y7(n)=Y7(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S107m(n,ii));
X6(n)=X6(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S106m(n,ii));
Y6(n)=Y6(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S106m(n,ii));
X10(n)=X10(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1010m(n,ii));
Y10(n)=Y10(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1010m(n,ii));
X11(n)=X11(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1011m(n,ii));
Y11(n)=Y11(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1011m(n,ii));
end

U7(n)=sqrt(X7(n).^2+Y7(n).^2);
U6(n)=sqrt(X6(n).^2+Y6(n).^2);
U10(n)=sqrt(X10(n).^2+Y10(n).^2);
U11(n)=sqrt(X11(n).^2+Y11(n).^2);

tong7(n)=(SS07i(n))+log(besseli(0,U7(n)));
tong6(n)=(SS06i(n))+log(besseli(0,U6(n)));
tong10(n)=(SS010i(n))+log(besseli(0,U10(n)));
tong11(n)=(SS011i(n))+log(besseli(0,U11(n)));

//определение квадранта принимаемого сигнала

[value, index(n)] = max([tong7(n) tong6(n) tong10(n)
tong11(n)]);
switch index(n)
case 1 //квадрант 1
cas(n)=1;
case 2 //квадрант 2
cas(n)=2;
case 3 //квадрант 3
cas(n)=3;
case 4 //квадрант 4
cas(n)=4;
end

if cas(n)==1
//работа в 1 квадранте с переносом созвездия в центр
SS07i(n)=0; X7(n)=0; Y7(n)=0;
SS06i(n)=0; X6(n)=0; Y6(n)=0;
SS010i(n)=0; X10(n)=0; Y10(n)=0;
SS011i(n)=0; X11(n)=0; Y11(n)=0;

for ii=1:4*M
SS07i(n)=SS07i(n)+(2/N0)*(S31(n,ii)+S1sm1(n,ii))*S107(n,ii);
 160

SS06i(n)=SS06i(n)+(2/N0)*(S31(n,ii)+S1sm1(n,ii))*S106(n,ii);
SS010i(n)=SS010i(n)+(2/N0)*(S31(n,ii)+S1sm1(n,ii))*S1010(n,ii);
SS011i(n)=SS011i(n)+(2/N0)*(S31(n,ii)+S1sm1(n,ii))*S1011(n,ii);

X7(n)=X7(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S104m(n,ii));
Y7(n)=Y7(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S104m(n,ii));
X6(n)=X6(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S103m(n,ii));
Y6(n)=Y6(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S103m(n,ii));
X10(n)=X10(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S107m(n,ii));
Y10(n)=Y10(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S107m(n,ii));
X11(n)=X11(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S108m(n,ii));
Y11(n)=Y11(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S108m(n,ii));
end

U7(n)=sqrt(X7(n).^2+Y7(n).^2);
U6(n)=sqrt(X6(n).^2+Y6(n).^2);
U10(n)=sqrt(X10(n).^2+Y10(n).^2);
U11(n)=sqrt(X11(n).^2+Y11(n).^2);

tong7(n)=(SS07i(n))+log(besseli(0,U7(n)));
tong6(n)=(SS06i(n))+log(besseli(0,U6(n)));
tong10(n)=(SS010i(n))+log(besseli(0,U10(n)));
tong11(n)=(SS011i(n))+log(besseli(0,U11(n)));

[value, index1(n)] = max([tong7(n) tong6(n) tong10(n)

tong11(n)]);
switch index1(n)
case 1 //квадрант 1-1
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 1;
case 2 //квадрант 1-2
d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 0;
d(4*n) = 1;
case 3 //квадрант 1-3
d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 0;
d(4*n) = 0;
case 4 //квадрант 1-4
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 0;

end
elseif cas(n)==2
//работа во 2 квадранте с переносом созвездия в центр
SS07i(n)=0; X7(n)=0; Y7(n)=0;
 161

SS06i(n)=0; X6(n)=0; Y6(n)=0;

SS010i(n)=0; X10(n)=0; Y10(n)=0;
SS011i(n)=0; X11(n)=0; Y11(n)=0;

for ii=1:4*M
SS07i(n)=SS07i(n)+(2/N0)*(S31(n,ii)+S1sm2(n,ii))*S107(n,ii);
SS06i(n)=SS06i(n)+(2/N0)*(S31(n,ii)+S1sm2(n,ii))*S106(n,ii);
SS010i(n)=SS010i(n)+(2/N0)*(S31(n,ii)+S1sm2(n,ii))*S1010(n,ii);
SS011i(n)=SS011i(n)+(2/N0)*(S31(n,ii)+S1sm2(n,ii))*S1011(n,ii);

X7(n)=X7(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S102m(n,ii));
Y7(n)=Y7(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S102m(n,ii));
X6(n)=X6(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S101m(n,ii));
Y6(n)=Y6(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S101m(n,ii));
X10(n)=X10(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S105m(n,ii));
Y10(n)=Y10(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S105m(n,ii));
X11(n)=X11(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S106m(n,ii));
Y11(n)=Y11(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S106m(n,ii));
end

U7(n)=sqrt(X7(n).^2+Y7(n).^2);
U6(n)=sqrt(X6(n).^2+Y6(n).^2);
U10(n)=sqrt(X10(n).^2+Y10(n).^2);
U11(n)=sqrt(X11(n).^2+Y11(n).^2);

tong7(n)=(SS07i(n))+log(besseli(0,U7(n)));
tong6(n)=(SS06i(n))+log(besseli(0,U6(n)));
tong10(n)=(SS010i(n))+log(besseli(0,U10(n)));
tong11(n)=(SS011i(n))+log(besseli(0,U11(n)));

[value, index2(n)] = max([tong7(n) tong6(n) tong10(n)

tong11(n)]);
switch index2(n)
case 1 //квадрант 2-1
d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 0;
d(4*n) = 1;
case 2 //квадрант 2-2
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 1;
case 3 //квадрант 2-3
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 0;
d(4*n) = 0;
case 4 //квадрант 2-4
d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 0;
 162

d(4*n) = 0;
end

elseif cas(n)==3
//работа в 3 квадранте с переносом созвездия в центр
SS07i(n)=0; X7(n)=0; Y7(n)=0;
SS06i(n)=0; X6(n)=0; Y6(n)=0;
SS010i(n)=0; X10(n)=0; Y10(n)=0;
SS011i(n)=0; X11(n)=0; Y11(n)=0;

for ii=1:4*M
SS07i(n)=SS07i(n)+(2/N0)*(S31(n,ii)+S1sm3(n,ii))*S107(n,ii);
SS06i(n)=SS06i(n)+(2/N0)*(S31(n,ii)+S1sm3(n,ii))*S106(n,ii);
SS010i(n)=SS010i(n)+(2/N0)*(S31(n,ii)+S1sm3(n,ii))*S1010(n,ii);
SS011i(n)=SS011i(n)+(2/N0)*(S31(n,ii)+S1sm3(n,ii))*S1011(n,ii);

X7(n)=X7(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1010m(n,ii));
Y7(n)=Y7(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1010m(n,ii));
X6(n)=X6(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S109m(n,ii));
Y6(n)=Y6(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S109m(n,ii));
X10(n)=X10(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1013m(n,ii));
Y10(n)=Y10(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1013m(n,ii));
X11(n)=X11(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1014m(n,ii));
Y11(n)=Y11(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1014m(n,ii));
end

U7(n)=sqrt(X7(n).^2+Y7(n).^2);
U6(n)=sqrt(X6(n).^2+Y6(n).^2);
U10(n)=sqrt(X10(n).^2+Y10(n).^2);
U11(n)=sqrt(X11(n).^2+Y11(n).^2);

tong7(n)=(SS07i(n))+log(besseli(0,U7(n)));
tong6(n)=(SS06i(n))+log(besseli(0,U6(n)));
tong10(n)=(SS010i(n))+log(besseli(0,U10(n)));
tong11(n)=(SS011i(n))+log(besseli(0,U11(n)));

[value, index3(n)] = max([tong7(n) tong6(n) tong10(n)

tong11(n)]);
switch index3(n)
case 1 //квадрант 3-1
d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 1;
d(4*n) = 0;
case 2 //квадрант 3-2
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 1;
d(4*n) = 0;
case 3 //квадрант 3-3
d(4*n-3) = 1;
d(4*n-2) = 1;
d(4*n-1) = 1;
 163

d(4*n) = 1;
case 4 //квадрант 3-4
d(4*n-3) = 1;
d(4*n-2) = 0;
d(4*n-1) = 1;
d(4*n) = 1;
end

elseif cas(n)==4
//работа в 4 квадранте с переносом созвездия в центр
SS07i(n)=0; X7(n)=0; Y7(n)=0;
SS06i(n)=0; X6(n)=0; Y6(n)=0;
SS010i(n)=0; X10(n)=0; Y10(n)=0;
SS011i(n)=0; X11(n)=0; Y11(n)=0;

for ii=1:4*M
SS07i(n)=SS07i(n)+(2/N0)*(S31(n,ii)+S1sm4(n,ii))*S107(n,ii);
SS06i(n)=SS06i(n)+(2/N0)*(S31(n,ii)+S1sm4(n,ii))*S106(n,ii);
SS010i(n)=SS010i(n)+(2/N0)*(S31(n,ii)+S1sm4(n,ii))*S1010(n,ii);
SS011i(n)=SS011i(n)+(2/N0)*(S31(n,ii)+S1sm4(n,ii))*S1011(n,ii);

X7(n)=X7(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1012m(n,ii));
Y7(n)=Y7(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1012m(n,ii));
X6(n)=X6(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1011m(n,ii));
Y6(n)=Y6(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1011m(n,ii));
X10(n)=X10(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1015m(n,ii));
Y10(n)=Y10(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1015m(n,ii));
X11(n)=X11(n)+(2*(m1+dm1)*A0/N0).*(P11(n,ii).*S1016m(n,ii));
Y11(n)=Y11(n)+(2*(m1+dm1)*A0/N0).*(P22(n,ii).*S1016m(n,ii));
end

U7(n)=sqrt(X7(n).^2+Y7(n).^2);
U6(n)=sqrt(X6(n).^2+Y6(n).^2);
U10(n)=sqrt(X10(n).^2+Y10(n).^2);
U11(n)=sqrt(X11(n).^2+Y11(n).^2);

tong7(n)=(SS07i(n))+log(besseli(0,U7(n)));
tong6(n)=(SS06i(n))+log(besseli(0,U6(n)));
tong10(n)=(SS010i(n))+log(besseli(0,U10(n)));
tong11(n)=(SS011i(n))+log(besseli(0,U11(n)));

[value, index4(n)] = max([tong7(n) tong6(n) tong10(n)

tong11(n)]);
switch index4(n)
case 1 //квадрант 4-1
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 1;
d(4*n) = 0;
case 2 //квадрант 4-2
d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 1;
 164

d(4*n) = 0;
case 3 //квадрант 4-3
d(4*n-3) = 0;
d(4*n-2) = 0;
d(4*n-1) = 1;
d(4*n) = 1;
case 4 //квадрант 4-4
d(4*n-3) = 0;
d(4*n-2) = 1;
d(4*n-1) = 1;
d(4*n) = 1;
end
end
end
error=sum(b ~= d) / N;

disp(error);
 165

Приложение 3. Программа моделирования адаптивного фильтра

(фрагмент)

…
// сумма: сигнал KAМ + гармоническая помеха + гауссовский шум
S3 = S1 + S2 + noise;

//формирование квадратур
S4=cos(2*%pi*fc*t);
S5=-sin(2*%pi*fc*t);
preI=S3.*S4;
preQ=S3.*S5;
X=preI+%i.*preQ;

p=10^(-6);
C=zeros(1,256);
C(1)=0.01;
W=zeros(1,256);

y=zeros(1,N*M);

for nn=2:N*M
for k=1:256
if nn>=k
W(k)=X(nn-k+1);
end
end
y(nn)=C*(W.');
C=C-p.*(abs((y(nn)-y(nn-1)))).*y(nn).*conj(W);
end

//формирование выходного сигнала АФ из квадратур

I2=real(y);
Q2=imag(y);
S33=I2.* cos(2*%pi*fc*t)- Q2.* sin(2*%pi*fc*t);
166
167
168