Радиотехника и электроника, 2019, T. 64, № 9, стр. 891-901

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

А. Ю. Беккиев 1*, В. И. Борисов 2**

1 АО “Росэлектроника”
121357 Москва, ул. Верейская, 29, стр. 141, Российская Федерация

2 АО “Концерн "Созвездие”
394018 Воронеж, ул. Плехановская, 14, Российская Федерация

* E-mail: aybekkiev@roselectronics.ru
** E-mail: bvi@sozvezdie.su

Поступила в редакцию 26.09.2018
После доработки 27.11.2018
Принята к публикации 13.12.2018

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

Аннотация

Построена вероятностно-временная модель оценки помехозащищенности систем радиосвязи в условиях радиоэлектронной борьбы. На основе модели получена замкнутая система уравнений, позволяющая оценить помехозащищенность каналов радиосвязи, а также предельные возможности станции ответных помех при подавлении систем радиосвязи с псевдослучайной перестройкой рабочей частоты (ППРЧ). Приведены численные расчеты зависимости скорости переключения частоты сигнала с ППРЧ от вероятности ошибки после подавления.

Ключевые слова: вероятностно-временная модель, системы связи с ППРЧ, вероятность обнаружения

1. ВВЕДЕНИЕ В ПРОБЛЕМУ

Начиная с 80-х годов XX века на развитие систем радиосвязи оказывает существенное влияние радиоэлектронная борьба (РЭБ), в целях которой обеспечивается обнаружение излучений средств радиосвязи и постановка помех (радиоэлектронное подавление) для нарушения функционирования линий радиосвязи. В общем случае процесс РЭБ включает в себя два основных последовательно выполняемых этапа: радиотехническую разведку и радиоподавление. В работе [1] помехозащищенность (ПЗ) радиотехнических систем предложено оценивать вероятностью

(1)
${{P}_{{{\text{пз}}}}} = 1 - {{P}_{{{\text{орг}}}}}{{P}_{{{\text{п}}f}}}{{P}_{{{\text{пд}}}}},$
где Pорг – вероятность того, что радиоподавление вообще будет организовано, Pпf – вероятность того, что спектр помех РЭБ будет попадать в полосу приемного устройства радиотехнической станции (РТС), Pпд – вероятность подавления РТС, т.е. вероятность того, что мощность помехи РЭБ на входе приемного устройства РТС будет превышать значение, достаточное для нарушения работы. Вероятность Pскр в [1] названа скрытностью излучения ((!))
(2)
${{P}_{{{\text{скр}}}}} = 1 - {{P}_{{{\text{п}}f}}},$
а вероятность Pпу – помехоустойчивостью РТС:

(3)
${{P}_{{{\text{пу}}}}} = 1 - {{P}_{{{\text{пд}}}}}.$

Учитывая определения (2) и (3), помехозащищенность выражение для Pпз (1) можно переписать в развернутом виде:

(4)
${{P}_{{{\text{пз}}}}} = 1 - (1 - {{P}_{{{\text{скр}}}}})(1 - {{P}_{{{\text{пу}}}}}){{P}_{{{\text{орг}}}}},$

раскрывающем основные составляющие понятия помехозащищенности.

Из соотношения (4) ясно, что помехозащищенность РТС определяется ее скрытностью, количественной мерой которой является вероятность Pскр, и помехоустойчивостью, количественной мерой которой является вероятность Pпу.

Недостаток существующего указанного выше подхода к оценке помехозащищенности состоит в том, что в нем не учитываются время функционирования РТС и время перехвата передачи информации в РТС станцией радиоэлектронного подавления (РЭП). Последнее можно назвать “реакцией” постановщика помех РЭБ. По этой причине помехозащищенность систем радиосвязи не может рассматриваться в отрыве от современных технических возможностей радиоэлектронной аппаратуры, методов радиотехнической разведки и возможностей постановщика помех. В связи с чем, все вопросы помехоустойчивости должны быть рассмотрены на примерах конкретной структурной схемы функционирования РТС и радиоуправления в условиях РЭБ. Приведем один из возможных обобщенных вариантов такой структурной схемы (рис. 1).

Рис. 1.

Обобщенная структурная схема радиолинии связи в условиях действия помех РЭБ, где ξ(t) – шум, η(t) – помеха СОП.

На рис. 1 обозначено: КИ – кодер для источника, КК – кодер для канала, ДК – декодер для канала, ДИ – декодер для источника, ДД – двоичные данные, τ – время задержки формирования помехи и доставки ее до входа приемника лини радиосвязи. На рис. 2 изображено воздействие помехой на сигнальный символ длительности T.

Рис. 2.

Эпюра воздействия помехой на сигнал длительности T.

Как следует из рисунка при воздействии преднамеренной помехи на сигнальный символ (в частности на скачок частоты псевдослучайной перестройкой рабочей частоты (ППРЧ) при ответной помехе) он поражается лишь частично. Эту особенность необходимо учитывать при оценке помехозащищенности радиолиний связи.

Время реакции постановщиков помех каждые пять лет сокращалось в 100…150 раз и на сегодняшний день составляет единицы миллисекунд. Обнаружитель станции РЭБ способен обнаруживать сигналы линии радиосвязи на интервале длительности одного скачка частоты сигнала с ППРЧ. Это потребовало разработки дополнительных мероприятий по повышению помехозащищенности линий радиосвязи, а также уточнения и дальнейшего развития некоторых положений статистической теории связи. В том числе:

1) разработка вероятностно-временной модели функционирования систем радиосвязи в условиях РЭБ, позволяющей производить оценку помехозащищенности, учитывать энергетические и временные возможности, как систем радиосвязи, так и станции радиоразведки и постановщика помех [2];

2) обоснование повышения помехоустойчивости за счет расширения спектра сигнала [35];

3) синтез оптимального приемника при скачкообразном изменении гауссовской помехи на интервале длительности единичного символа [2].

2. ВЕРОЯТНОСТНО-ВРЕМЕННАЯ МОДЕЛЬ

Вероятностно-временная модель функционирования систем радиосвязи в условиях РЭБ позволяет проводить оценку помехозащищенности, учитывать энергетические и временные возможности как систем радиосвязи, так и станций разведки и постановки помех (см. рис. 1).

Получим общее выражение для вероятности ошибочного различения символов в условиях действия помех РЭБ на конечном времени передачи T = NT0, где T0 – длительность одного символа, N – число символов. В частном случае это может быть скачок частоты медленной (по отношению к скорости передачи) ППРЧ. Для этого обозначим через q0 и q1 соответственно отношение сигнал/шум (ОСШ) без воздействия помехи и при ее воздействии, а соответствующие им вероятности ошибочного различения – через PE0 = = P(q0) и PE1 = P(q1). Средняя вероятность ошибки может быть вычислена по формуле полной вероятности. Обозначим через H0 гипотезу о том, что данный символ (бит) не поражен помехой, а через H1 – о том, что поражен. Из рис. 3 видно, что в случае обнаружения сигнала станцией ответных помех (СОП) (за время τ с вероятностью Pоб(τ)) количество не пораженных помехой символов равно k, а количество пораженных равно − k (поскольку T1 = T − τ = (Nk)T0). Тогда (при условии, что сигнал обнаружен) условная вероятность того, что символ поражен помехой, равна

(5)
$\frac{{N - k}}{N} = \frac{{(N - k){{T}_{0}}}}{{N{{T}_{0}}}} = \frac{{T - \tau }}{T},$
а условная вероятность того, что символ не поражен, равна

(6)
$\frac{k}{N} = \frac{{k{{T}_{0}}}}{{N{{T}_{0}}}} = \frac{\tau }{T}.$
Рис. 3.

Эпюры кодовых символов и сигналов (а, б), шумов (в), запаздывающей помехи (г), дисперсии шумов и помех (д), ОСШ (е) и вероятностей ошибки различения PE0 и PE1 (ж).

Вероятность гипотезы H0 будет равна вероятности суммы событий: сигнал не обнаружен за время τ или обнаружен, но данный символ не поражен помехой. Вероятность такого события равна

(7)
$\begin{gathered} P({{H}_{0}}) = [1 - {{P}_{{{\text{об}}}}}(\tau )] + {{P}_{{{\text{об}}}}}(\tau )\frac{\tau }{T} = \\ = 1 - {{P}_{{{\text{об}}}}}(\tau )\frac{{T - \tau }}{T}. \\ \end{gathered} $

Аналогично, гипотеза H1 равна произведению событий: сигнал обнаружен и данный символ поражен помехой. Вероятность такого события равна

(8)
$P({{H}_{1}}) = {{P}_{{{\text{об}}}}}(\tau )\frac{{T - \tau }}{T}.$

Средняя вероятность ошибки на символ после воздействия помехи РЭБ согласно формуле полной вероятности будет равна

(9)
$\begin{gathered} {{{\bar {P}}}_{E}} = P({{H}_{0}}){{P}_{{E0}}} + P({{H}_{1}}){{P}_{{E1}}} = \\ = {{P}_{{E0}}} + {{P}_{{{\text{об}}}}}(\tau )\frac{{T - \tau }}{T}({{P}_{{E1}}} - {{P}_{{E0}}}). \\ \end{gathered} $

Как следует из формулы (9), средняя вероятность ошибки, в отличие от (4), определяется не только энергетическими характеристиками (отношением сигнал/шум) на входе приемника радиолинии, но и относительным временем поражения символа помехой. Помехоустойчивость радиолинии в этом случае существенно зависит от предельных возможностей станции помех, в частности от времени запаздывания τ и мощности генерируемой помехи.

3. ОЦЕНКА ВРЕМЕНИ ЗАДЕРЖКИ

Задержка за счет разности хода лучей. Время работы без помех τ складывается из следующих составляющих: τ01 – задержка из-за разности хода лучей (передатчик → приемник и передатчик → → станция РЭП → приемник), τ02 – задержка сигнала в тракте приемника обнаружителя. Основной вклад в величину задержки сигнала в тракте приемника вносят времена задержек фильтров основной селекции и время накопления энергии для обнаружения с заданной вероятностью. Проведем оценку данных величин.

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

(10)
${{\tau }_{{01}}} = \frac{1}{c}({{d}_{{J1}}} + {{d}_{{J2}}} - {{d}_{{12}}}),$
где c – скорость света, dJ1 – расстояние между передатчиком и СОП, dJ2 – расстояние между СОП и приемником РТС, d12 – расстояние между передатчиком и приемником (см. рис. 4, где обозначено: ПРД(1) – передатчик РТС, ПРМ(2) – приемник РТС). Цифра 1 обозначает передатчик радиолинии, цифра 2 – приемник, буква J – станцию ответных помех.

Рис. 4.

Геометрия размещения на местности РТС и СОП.

Задержки в аппаратуре радиоразведки СОП. Как отмечено в работах [6, 7], благодаря высоким темпам научно-технического прогресса за последние несколько десятков лет произошло значительное конструктивное и технологическое усовершенствование радиотехнических устройств приема и обработки сигналов, входящих в состав аппаратуры подсистемы радиоразведки (РР) СОП. Использование современной элементной базы, цифровой обработки сигналов, высокопроизводительных сигнальных процессоров и программируемых логических интегральных схем (ПЛИС) позволяет реализовать на практике широкий спектр возможностей по решению задач РР.

В подсистеме РР СОП используются многоканальные обнаружители-пеленгаторы – радиотехнические устройства, с помощью которых решается задача обнаружения и пеленгования источников радиоизлучения (ИРИ). В состав обнаружителя-пеленгатора входят многоканальная антенная система (АС), подключенная к входам многоканального радиоприемного устройства, и модуль цифровой обработки сигналов (ЦОС). Обнаружение и пеленгование ИРИ осуществляется по результатам приема и обработки сигналов ИРИ, принятых в пространственно разнесенных пунктах – приемных элементах АС обнаружителя-пеленгатора.

С помощью обнаружителя-пеленгатора в подсистеме РР станции ответных помех выполняется частотный поиск сигналов ИРИ в спектральной области путем сканирования рабочего диапазона частот. Широкое применение на практике получили моноимпульсные обнаружители-пеленгаторы, с помощью которых осуществляется синхронный прием временных реализаций в текущей частотной полосе мгновенного анализа (ПМА) во всех пространственных каналах, оцифровка реализаций и их быстрое преобразование Фурье (БПФ).

Каждый отсчет БПФ представляет собой комплексную амплитуду в элементарном частотном канале (ЭЧК), ширина полосы которого обратно пропорциональна длительности временной реализации. Совокупность отсчетов БПФ во всех пространственных каналах обнаружителя-пеленгатора, принадлежащих одному и тому же ЭЧК, характеризует амплитудно-фазовое распределение падающей на АС обнаружителя-пеленгатора радиоволны на частоте данного ЭЧК. Спектральное представление временных реализаций обеспечивает возможность определения спектрального состава радиосигналов (т.е. совокупности отсчетов БПФ, принадлежащим сигналу данного ИРИ) в отсутствии информации о законе их модуляции и способе кодирования информации. Как отмечено в [6], процедура БПФ позволяет “снять” модуляцию сигнала, раскладывая его в базисе ортогональных нормированных функций, и выполнять при ведении РР обработку в условиях имеющейся параметрической и непараметрической априорной неопределенности относительно формы принимаемого радиосигнала. При этом непараметрическая неопределенность относительно формы сигнала преобразуется к параметрической неопределенности для каждого ЭЧК относительно КА напряженности поля волны от ИРИ и направления ее прихода [6].

В рамках задачи обнаружения сигналов в подсистеме РР СОП среди совокупности отсчетов БПФ текущей ПМА определяются “сигнальные” отсчеты, в амплитуде и фазе которых содержится информация об порожденной ИРИ волне. Ввиду того, что ИРИ генерируют модулированные сигналы (а также с учетом влияния весового окна цифровой обработки), сигнал одного и того же ИРИ представлен несколькими элементами разрешения по частоте, что соответствует нескольким отсчетам БПФ. Поэтому по результатам обнаружения с целью повышения ОСШ выполняется группирование (отождествление) отсчетов БПФ по принадлежности к сигналу одного и того же ИРИ во всех пространственных каналах обнаружителя – пеленгатора [6]. Последующее пеленгование и местоопределение каждого обнаруженного ИРИ осуществляется по результатам совместной обработки данных во всех пространственных каналах, содержащих спектр сигнала ИРИ.

С учетом приведенного выше описания основных этапов обработки принимаемых радиосигналов задержка в подсистеме РР современных (перспективных) СОП складывается из задержки t0 сигнала в радиоприемном тракте обнаружителя-пеленгатора и длительности τ02 принимаемой временной реализации. Отметим, что в целях повышения быстродействия цифровой обработки сигналов в аппаратуре подсистемы РР СОП реализуется высокоскоростная конвейерная потоковая обработка данных текущей полосы мгновенного анализа на ПЛИС и сигнальных процессорах. При этом время обработки данных в большинстве случаев на порядок меньше длительности принимаемой временной реализации и в дальнейшем не рассматривается. Длительность временной реализации определяет потенциальные возможности подсистемы РР СОП по обнаружению радиосигналов и выбирается исходя из условия обеспечения требуемых показателей эффективности пространственно-многоканального обнаружения сигнала в спектральной области.

Характеристики пространственно-многоканального обнаружения радиосигналов. Несмотря на отсутствие в подсистеме РР СОП априорной информации о форме (законах амплитудной и фазовой модуляции) сигнала ИРИ, обеспечиваемая многоканальным моноимпульсным обнаружителем-пеленгатором возможность пространственно разнесенного приема сигналов позволяет реализовать обнаружение сигналов в спектральной области на основе следующего отличия в природе сигнальных и помеховых (шумовых) составляющих спектральных отсчетов. Сигнальная составляющая отсчетов характеризует пространственное распределение амплитуды и фазы поля волны ИРИ по раскрыву антенной системы обнаружителя-пеленгатора. Компоненты помеховой составляющей в пространственно разнесенных пунктах приема имеют случайные амплитуды и фазы, не обусловленные падением некоторой радиоволны с плоским волновым фронтом от удаленного ИРИ. Поэтому особенностью задачи обнаружения сигналов ИРИ в подсистеме радиоразведки является обнаружение не произвольной сигнальной составляющей в наблюдаемых данных (НД), а составляющей, обусловленной плоским волновым фронтом волны ИРИ.

Задача обнаружения сигналов ИРИ с помощью пространственно-многоканальной АС обнаружителя-пеленгатора представляет собой бинарную задачу проверки следующих статистических гипотез:

1) гипотеза γ1 о том, что в НД присутствуют сигнальные составляющие от ИРИ, характеризуемого истинным направлением прихода радиоволны α0 = (θ00)T:

(11)
${{\gamma }_{1}}:\,\,\left\{ {{\mathbf{\dot {U}}} = \dot {E}{\kern 1pt} {\mathbf{\dot {H}}}\left( {{{\alpha }_{0}}} \right) + \dot {\xi }} \right\};$

2) гипотеза γ0 о том, что в НД отсутствуют сигнальные составляющие от ИРИ:

(12)
${{\gamma }_{0}}:\,\,\left\{ {{\mathbf{\dot {U}}} = \dot {\xi }} \right\}.$
Здесь ${\mathbf{\dot {U}}}$ – вектор комплексных амплитуд сигнала на выходах многоканальной антенной системы обнаружителя-пеленгатора; ${\mathbf{\dot {H}}}\left( \alpha \right)$ – векторная комплексная диаграмма направленности (ВКДН) АС обнаружителя-пеленгатора [6, 7]; $\dot {E}$ – комплексная амплитуда напряженности электрического поля радиоволны в центре пространственной системы координат, в которых определена диаграмма направленности; $\dot {\xi }$ – комплексный вектор аддитивного гауссовского шума с нулевым средним и произвольной эрмитовой корреляционной матрицей; θ(β) – азимут (угол места) прихода радиоволны, T – знак транспонирования. ВКДН АС обнаружителя-пеленгатора представляет собой N-мерный вектор-функцию, элементы которого являются коэффициентами преобразования комплексной амплитуды (КА) напряженности электрического (либо магнитного) поля радиоволны в КА тока на номинальных нагрузках, подключенных к выходам АС обнаружителя-пеленгатора.

Как отмечено выше, по результатам пространственно-многоканального обнаружения среди совокупности отсчетов БПФ текущей ПМА определяются “сигнальные” отсчеты, в амплитуде и фазе которых содержится информация об генерируемой ИРИ волне.

Отождествленные к сигналу одного и того же ИРИ спектральные отсчеты образуют спектр обнаруженного сигнала, а их количество R определяется длительностью временной реализации τ (длительностью интервала времени контакта с сигналом) и шириной полосы частот W радиосигнала в соответствии с равенством R = [Wτ], где [⋅] – целая часть числа.

Можно показать [8, 9], что решающая статистика обнаружения радиосигналов в подсистеме РР СОП имеет следующий вид:

(13)
$\Xi = M(\hat {\alpha }),$
где
(14)
$\hat {\alpha } = \arg \mathop {\sup }\limits_\alpha \left[ {M\left( \alpha \right)} \right]$
– максимально правдоподобная оценка азимута и угла места направления прихода радиоволны ИРИ;
(15)
$M\left( \alpha \right) = \frac{{{{{{\mathbf{\dot {H}}}}}^{{\rm H}}}\left( \alpha \right){{{{\mathbf{\dot {K}}}}}^{{ - 1}}}{\mathbf{\dot {W}}}{{{{\mathbf{\dot {K}}}}}^{{ - 1}}}{\mathbf{\dot {H}}}\left( \alpha \right)}}{{{{{{\mathbf{\dot {H}}}}}^{{\rm H}}}\left( \alpha \right){{{{\mathbf{\dot {K}}}}}^{{ - 1}}}{\mathbf{\dot {H}}}\left( \alpha \right)}} - $
пеленгационный рельеф (ПР) [6, 7]; ${\mathbf{\dot {W}}}$ – матрица взаимных энергий измеренных КА, накопленная по R спектральным компонентам сигнала; ${\mathbf{\dot {K}}}$ – матрица пространственной ковариации аддитивного шума, Н – знак эрмитового сопряжения.

Решение о реализации гипотезы γ1 против γ0 принимается в случае выполнения неравенства

(16)
$\Xi \geqslant h,$
где порог обнаружения выбирается в соответствии с критерием Неймана–Пирсона и обеспечивает постоянный уровень ложной тревоги.

Решающая статистика обнаружения (13) представляет собой квадратичную форму, построенную на комплексных гауссовских коррелированных величинах [9]. В соответствии с доказанной в [10] теоремой о распределении положительно определенной квадратичной формы, построенной на комплексных гауссовских коррелированных величинах, решающая статистика (13) подчиняется обобщенному χ2-распределению с 2R степенями свободы, параметром нецентральности [10]

(17)
$\lambda = \left( {\sum\limits_{r = 0}^{R - 1} {{{{\left| {{{{\dot {A}}}_{r}}} \right|}}^{2}}} } \right){{{\mathbf{\dot {H}}}}^{H}}\left( {\hat {\alpha }} \right){{{\mathbf{\dot {K}}}}^{{ - 1}}}{\mathbf{\dot {H}}}\left( {\hat {\alpha }} \right)$
и дисперсиями реальных и мнимых частей гауссовских компонент распределения, равными 1/2. Здесь ${{\dot {A}}_{r}}$ – КА отсчета БПФ, принадлежащего радиосигналу ИРИ, $r = \overline {0,R - 1} .$

Для оценки потенциальных возможностей подсистемы РР СОП рассмотрим случай, когда шумы радиоприемных трактов обнаружителя-пеленгатора можно считать статистически независимыми с одинаковой интенсивностью σ2. Учтем также, что в подсистеме РР современных (перспективных) станций ответных помех широко используются моноимпульсные обнаружители-пеленгаторы с эквидистантной кольцевой антенной решеткой (ЭКАР). В отсутствии взаимных электродинамических влияний антенных элементов для ВКДН ЭКАР справедливо равенство

(18)
${{{\mathbf{\dot {H}}}}^{H}}\left( {\hat {\alpha }} \right){\mathbf{\dot {H}}}\left( {\hat {\alpha }} \right) = {{\left( {{{l}_{0}}} \right)}^{2}}N,$
где l0 – эффективность (“нагруженная” действующая длина) вибратора ЭКАР. Тогда параметр нецентральности решающей статистики обнаружения можно записать как

(19)
$\lambda = \left( {\sum\limits_{r = 0}^{R - 1} {\frac{{{{{\left| {{{{\dot {A}}}_{r}}} \right|}}^{2}}}}{{{{\sigma }^{2}}}}} } \right){{\left( {{{l}_{0}}} \right)}^{2}}N = 2\frac{{\sum\limits_{r = 0}^{R - 1} {{{{\left| {{{{\dot {A}}}_{r}}{{l}_{0}}} \right|}}^{2}}} }}{{2{{\sigma }^{2}}R}}RN.$

Обозначим ${{P}_{s}} = 0.5\sum\nolimits_{r = 0}^{R - 1} {{{{\left| {{{{\dot {A}}}_{r}}{{l}_{0}}} \right|}}^{2}}} $ и Pn = σ2R – мощность радиосигнала и мощность шума соответственно в полосе частот, образованной R спектральными компонентами радиосигнала. С учетом введенных обозначений перепишем правую часть (19) как

(20)
$\lambda = 2\frac{{{{P}_{s}}}}{{{{P}_{n}}}}RN = 2N\left[ {W\tau } \right]\gamma ,$
где γ = Ps/Pn – ОСШ по мощности на входе радиоприемного канала обнаружителя-пеленгатора в полосе частот обнаруженного радиосигнала.

Таким образом, за счет совместной обработки радиосигналов при пространственно-многоканальном обнаружении энергетическое ОСШ на выходе обнаружителя, численно равное параметру нецентральности (19), по сравнению со случаем одноканального энергетического обнаружения увеличивается в N раз. Решающая статистика обнаружения подчиняется обобщенному χ2-распределению с 2[Wτ] степенями свободы и дисперсиями реальных и мнимых частей гауссовских компонент распределения, равными 1/2. В случае наличия радиосигнала распределение является нецентральным с параметром (20).

Точные аналитические выражения для вероятностей ложной тревоги Pf и правильного обнаружения Pd радиосигнала могут быть выражены в виде разложения по центральному χ2-распределению, как в [9], либо через обобщенную Q-функцию Маркума, как в [4]

(21)
${{P}_{f}} = \int\limits_h^\infty {{{W}_{\Xi }}(\left. z \right|{{\gamma }_{0}})dz} = {{Q}_{R}}\left( {0,\sqrt h } \right),$
(22)
${{P}_{d}} = 1 - \int\limits_0^h {{{W}_{\Xi }}(\left. z \right|{{\gamma }_{1}})dz} = {{Q}_{R}}\left( {\sqrt \lambda ,\sqrt h } \right),$
где ${{W}_{\Xi }}(\left. z \right|{{\gamma }_{i}})$ – одномерная плотность вероятности статистики Ξ при реализации гипотезы γi, i = 0, 1:
${{W}_{\Xi }}(\left. z \right|{{\gamma }_{0}}) = \frac{{{{z}^{{R - 1}}}\exp \left( { - z} \right)}}{{\Gamma \left( R \right)}},$
$\begin{gathered} {{W}_{\Xi }}(\left. z \right|{{\gamma }_{1}}) = \frac{1}{2}{{\left( {\frac{z}{\lambda }} \right)}^{{(R - {{1)} \mathord{\left/ {\vphantom {{1)} 2}} \right. \kern-0em} 2}}}}\exp [ - (z + {{\lambda )} \mathord{\left/ {\vphantom {{\lambda )} 2}} \right. \kern-0em} 2}] \times \\ \times \,\,{{I}_{{R - 1}}}(\sqrt {z\lambda } ),\,\,\,\,\,z \geqslant 0, \\ \end{gathered} $
IR(⋅) – модифицированная функция Бесселя 1‑го рода порядка R [11], QR(a,b) – обобщенная Q‑функция Маркума [4]:

${{Q}_{R}}(a,b) = \frac{1}{{{{a}^{{R - 1}}}}}\int\limits_b^\infty {{{x}^{R}}\exp \left( { - \frac{{{{x}^{2}} + {{a}^{2}}}}{2}} \right){{I}_{{R - 1}}}(ax)dx} .$

При фиксированных вероятностях ложной тревоги и правильного обнаружения минимальная длительность принимаемой временной реализации (длительность контакта с сигналом) определяется из выражений (21) и (22). При выполнении условия [Wτ] $ \gg $ 1 решающая статистика обнаружения может быть аппроксимирована гауссовским распределением [11]. В этом случае вероятности ложной тревоги Pf и правильного обнаружения Pd определяются формулами

(23)
${{P}_{f}} = {{Q}_{R}}\left( {0,\sqrt h } \right) \approx Q\left( {\frac{{h - 2R}}{{2\sqrt R }}} \right),$
(24)
${{P}_{d}} = {{Q}_{R}}\left( {\sqrt \lambda ,\sqrt h } \right) \approx Q\left( {\frac{{h - 2R - \lambda }}{{2\sqrt {R + \lambda } }}} \right),$
где Q(x) – Q-функция Гаусса или дополнительная функция к гауссовскому интегралу вероятностей:

$Q(x) = \frac{1}{{\sqrt {2\pi } }}\int\limits_x^\infty {\exp ({{ - {{t}^{2}}} \mathord{\left/ {\vphantom {{ - {{t}^{2}}} 2}} \right. \kern-0em} 2})dt} .$

Если придать фиксированные значения вероятностям (23) и (24), то в результате применения к правым и левым частям данных выражений функции, обратной к Q(⋅), получим

(25)
$h - 2R = 2\sqrt R {{Q}^{{ - 1}}}\left( {{{P}_{f}}} \right),$
(26)
$h - 2R - \lambda = 2\sqrt {R + \lambda } {{Q}^{{ - 1}}}\left( {P_{d}^{{}}} \right).$

Совместное решение уравнений (25) и (26) относительно параметра R с учетом соотношения (20) приводит к выражению

(27)
$R = \frac{{{{{\left( {{{Q}^{{ - 1}}}\left( {{{P}_{f}}} \right) - \sqrt {1 + 2\gamma N} {{Q}^{{ - 1}}}\left( {{{P}_{d}}} \right)} \right)}}^{2}}}}{{{{\gamma }^{2}}{{N}^{2}}}}.$

Поскольку R = [Wτ], то необходимая для обеспечения заданных вероятностей ложной тревоги и правильного обнаружения сигнала длительность временной выборки в подсистеме РР СОП определяется по формуле

(28)
${{\tau }_{{02}}} = \frac{{{{{\left( {{{Q}^{{ - 1}}}\left( {{{P}_{f}}} \right) - \sqrt {1 + 2\gamma N} {{Q}^{{ - 1}}}\left( {{{P}_{d}}} \right)} \right)}}^{2}}}}{{{{\gamma }^{2}}{{N}^{2}}W}}.$

На рис. 5а приведены зависимости безразмерной величины τ02W от входного ОСШ γ (в децибелах) для различных значений вероятностей обнаружения Pd: 0.8 (кривая 1), 0.9 (кривая 2), 0.99 (кривая 3), 0.999 (кривая 4) при вероятности ложной тревоги Pf = 10−7 и N = 1. Аналогичные зависимости приведены на рис. 5б для N = 5.

Рис. 5.

Зависимость параметра ${{\tau }_{{01}}}W$ от отношения сигнал–шум на входе радиоприемного канала обнаружителя-пеленгатора.

Как следует из приведенных рисунков и формулы (28), необходимое время накопления обратно пропорционально квадрату числа пространственных каналов для малого ОСШ и обратно пропорционально для большого ОСШ.

Таким образом, число пространственно разнесенных каналов является определяющим фактором, влияющим на время реакции СОП при обнаружении слабых сигналов.

Задержка сигнала в радиоприемном тракте обнаружителя-пеленгатора. Радиоприемный тракт обнаружителя-пеленгатора использует многоканальную систему, состоящую из набора полосовых фильтров (ПФ). Каждый фильтр представляет собой каскад из N1 звеньев второго порядка с передаточными функциями Hi(s), i = 1, 2, …, N1. Результирующая передаточная функция всего фильтра равна

(29)
$H(s) = \prod\limits_{i = 1}^{{{N}_{1}}} {{{H}_{i}}(s)} .$

Частотную характеристику, которая получается из (29) при s = jω, можно выразить через амплитудно-частотную $\left| {H\left( {i\omega } \right)} \right|$ и фазо-частотную $\varphi (\omega )$ характеристики следующим образом:

(30)
$\begin{gathered} H(j\omega ) = \left| {H\left( {i\omega } \right)} \right|\exp \left[ {j\varphi (\omega )} \right] = \\ = \prod\limits_{i = 1}^{{{N}_{1}}} {\left| {{{H}_{i}}(j\omega )} \right|\exp [j{{\varphi }_{i}}(\omega )]\,} . \\ \end{gathered} $

Групповое время задержки tф (время запаздывания сигнала) всего фильтра равно

(31)
$\begin{gathered} {{t}_{{\text{ф}}}} = - \frac{d}{{d\omega }}\varphi (\omega ) = - \sum\limits_{i = 1}^N {\frac{{d{{\varphi }_{i}}(\omega )}}{{d\omega }}} = \\ = {{t}_{{{\text{ф}}1}}} + {{t}_{{{\text{ф}}2}}} + \ldots + {{t}_{{{\text{ф}}{{N}_{1}}}}}, \\ \end{gathered} $
где ${{t}_{{{\text{ф}}i}}} = {{ - d{{\varphi }_{i}}(\omega )} \mathord{\left/ {\vphantom {{ - d{{\varphi }_{i}}(\omega )} {d\omega }}} \right. \kern-0em} {d\omega }}$ – групповое время задержки i-го звена. Таким образом, групповое время задержки фильтра равно сумме времен задержки звеньев.

Известно [12], что частотная характеристика типового звена второго порядка полосно-пропускающего фильтра описывается соотношением

(32)
${{H}_{i}}(j\omega ) = \frac{{({{{{\omega }_{0}}} \mathord{\left/ {\vphantom {{{{\omega }_{0}}} Q}} \right. \kern-0em} Q})j\omega }}{{{\kern 1pt} - {{\omega }^{2}} + ({{{{\omega }_{0}}} \mathord{\left/ {\vphantom {{{{\omega }_{0}}} Q}} \right. \kern-0em} Q})j\omega + \omega _{0}^{2}}} = \frac{{ja}}{{b + ja}},$
где ω0 – центральная частота, Q = ω0/Δω – добротность контура, $\Delta = {{\omega }_{В}} - {{\omega }_{Н}},$ ${{\omega }_{В}},\;{{\omega }_{Н}}$ – верхняя и нижняя частоты по уровню 3 дБ, a = (ω0/Q)ω, $b = \omega _{0}^{2} - {{\omega }^{2}}.$

(33)
${{\varphi }_{i}}(\omega ) = {\text{arctg}}\frac{{{{H}_{{\operatorname{Im} }}}(j\omega )}}{{{{H}_{{\operatorname{Re} }}}(j\omega )}} = {\text{arctg}}\frac{b}{a}.$

После дифференцирования (33) и соответствующих преобразований получаем

(34)
${{\tau }_{i}}(\omega ) = - \frac{{d{{\varphi }_{i}}(\omega )}}{{d\omega }} = \frac{{{{\omega }_{0}}}}{Q}\frac{{\omega _{0}^{2} + {{\omega }^{2}}}}{{{{\omega }^{2}}{{{({{{{\omega }_{0}}} \mathord{\left/ {\vphantom {{{{\omega }_{0}}} Q}} \right. \kern-0em} Q}{\kern 1pt} )}}^{2}} + {{{(\omega _{0}^{2} - {{\omega }^{2}})}}^{2}}}}.$

Анализ зависимости групповой задержки по частоте показывает, что максимальная задержка происходит на частоте ω = ω0, не зависит от центральной частоты и равна

(35)
${{\tau }_{i}}({{\omega }_{0}}) = 2\frac{Q}{{{{\omega }_{0}}}} = \frac{2}{{\Delta \omega }}.$

Групповое время задержки всего каскада из N1 звеньев второго порядка ограничено величиной

(36)
${{t}_{0}} = {{t}_{{\text{ф}}}} = \frac{{2{{N}_{1}}}}{{\Delta \omega }}.$

4. ПРИМЕНЕНИЕ ВЕРОЯТНОСТНО-ВРЕМЕННОЙ МОДЕЛИ К ОБОСНОВАНИЮ ТРЕБОВАНИЙ К СИСТЕМАМ СВЯЗИ С ППРЧ

Естественной реакцией РТС на угрозу ответных помех со стороны СОП в системах связи с ППРЧ является увеличение скорости перестройки рабочей частоты. Однако в этом случае возникают следующие проблемы [13]: 1) увеличение стоимости синтезатора частот при ухудшении его надежности, 2) усложнение устройств синхронизации, 3) дополнительные затраты времени при переключении с частоты на частоту.

Для оценки необходимой скорости переключения ППРЧ предположим, что фиксировано расстояние между передатчиком и приемником d12. Определим вначале положение СОП, обеспечивающее заданный коэффициент перекрытия сигнала ответной помехой ρ:

(37)
$\rho = \frac{{T - \tau }}{T} = \frac{{T - {{\tau }_{{01}}} - ({{\tau }_{{02}}} + {{t}_{0}})}}{T},$
откуда с учетом (10) следует
(38)
$\frac{1}{c}({{d}_{{J1}}} + {{d}_{{J2}}} - {{d}_{{12}}}) + ({{\tau }_{{02}}} + {{t}_{0}}) = (1 - \rho )T$
или

(39)
${{d}_{J}}_{1} + {{d}_{J}}_{2} = {{d}_{{12}}} + [(1 - r)T - ({{t}_{{02}}} + {{t}_{0}})]c.$

Уравнение (39) позволяет определить область размещения СОП, изнутри которой с заданным коэффициентом перекрытия ρ может подавляться сигнал на входе приемника РТС. Эта область представляет собой эллипс, в фокусах которого находятся передатчик и приемник РТС, а полуоси определяются следующими выражениями

(40)
$\begin{gathered} a = \frac{1}{2}\{ [(1 - \rho )T - ({{\tau }_{{02}}} + {{t}_{0}})]c + {{d}_{{12}}}\} = \frac{1}{{2\varepsilon }}{{d}_{{12}}}, \\ b = \sqrt {{{a}^{2}} - {{d_{{12}}^{2}} \mathord{\left/ {\vphantom {{d_{{12}}^{2}} 4}} \right. \kern-0em} 4}} , \\ \end{gathered} $
где ε – эксцентриситет эллипса.

Пусть требуется создать СОП из произвольной внутренней точки эллипса с фокусным расстоянием d12 и эксцентриситетом ε такую ответную помеху, которая делает вероятность приема сигнала в РТС больше некоторой предельно допустимой величины Pдоп. В соответствии с вероятностно-временной моделью при расположении СОП внутри эллипса, в фокусах которого находятся передатчик и приемник РТС, должно выполняться неравенство

(41)
${{\bar {P}}_{E}} = {{P}_{{E0}}} + \frac{{T - \tau }}{T}{{P}_{{{\text{об}}}}}({{\tau }_{{02}}})({{P}_{{E1}}} - {{P}_{{E0}}}) \geqslant {{P}_{{{\text{доп}}}}}.$

Учитывая, что $\rho = {{(T - \tau )} \mathord{\left/ {\vphantom {{(T - \tau )} T}} \right. \kern-0em} T}$ из неравенства (41) определяется необходимый коэффициент перекрытия сигнала ответной помехой

(42)
$1 \geqslant \rho \geqslant \frac{{{{P}_{{{\text{доп}}}}} - {{P}_{{E0}}}}}{{{{P}_{{{\text{об}}}}}({{\tau }_{{02}}})({{P}_{{E1}}} - {{P}_{{E0}}})}}.$

Из неравенства (42) вытекает ограничение на вероятность обнаружения сигнала

(43)
${{P}_{{{\text{об}}}}}({{\tau }_{{02}}}) \geqslant \frac{{{{P}_{{{\text{доп}}}}} - {{P}_{{E0}}}}}{{({{P}_{{E1}}} - {{P}_{{E0}}})}}.$

После подстановки (42) в (40) получим условия существования эллипса с заданными вероятностными характеристиками подавления

(44)
$\begin{gathered} T\left[ {1 - \frac{{{{P}_{{{\text{доп}}}}} - {{P}_{{E0}}}}}{{{{P}_{{{\text{об}}}}}(\tau )({{P}_{{E1}}} - {{P}_{{E0}}})}}} \right] - ({{\tau }_{{02}}} + {{t}_{0}}) - \\ - \,\,\frac{{1 - \varepsilon }}{\varepsilon }\frac{{{{d}_{{12}}}}}{c} > 0. \\ \end{gathered} $

Исходя из условий существования эллипса с заданными свойствами, получается ряд ограничений на параметры. Рассмотрим ограничения на скорость переключения частоты сигналов с ППРЧ, которую обозначим через ${{R}_{h}} = {1 \mathord{\left/ {\vphantom {1 {{{T}_{h}}}}} \right. \kern-0em} {{{T}_{h}}}}$ (число скачков в секунду). Длительность скачка частот Th складывается из двух составляющих – времени переключения Tsw и времени передачи данных Tdw [13]:

(45)
${{T}_{h}} = {{T}_{{dw}}} + {{T}_{{sw}}}.$

Для военных систем связи типа SINCGARS доля времени, в течение которого передаются данные, составляет [13]

(46)
$\frac{{{{T}_{{dw}}}}}{{{{T}_{h}}}} \approx 0.8,$

т.е. временной параметр T в (44) связан с длительностью скачка частоты Th соотношением:

(47)
$T = {{T}_{{dw}}} = 0.8{{T}_{h}},$

откуда, с учетом (42), получается верхняя граница для скорости переключения ППРЧ, при которой возможно подавление РТС до требуемой вероятности ошибки

(48)
$\begin{gathered} {{R}_{h}} \leqslant \frac{{0.8}}{{\left( {{{\tau }_{{02}}} + {{t}_{0}} + \frac{{1 - \varepsilon }}{\varepsilon }\frac{{{{d}_{{12}}}}}{c}} \right)}} \times \\ \times \,\,\left( {1 - \frac{{{{P}_{{{\text{доп}}}}} - {{P}_{{E0}}}}}{{{{P}_{{{\text{об}}}}}{\text{(}}{{\tau }_{{\text{0}}}})[{{P}_{{E1}}} - {{P}_{{E0}}}]{\kern 1pt} }}} \right). \\ \end{gathered} $

Из анализа (48) следует, что при заданных вероятностях ошибки на бит в РТС до подавления PE0 и допустимой вероятности ошибки Pдоп остается свободный параметр PE1, который определяется мощностью ответной помехи. Действительно, если пренебречь затуханиями при распространении сигнала, то для бинарной частотной модуляции вероятность ошибки равна

(49)
${{P}_{{E1}}} = \frac{1}{2}\exp \left\{ { - \frac{1}{2}\frac{{{{P}_{s}}}}{{\sigma _{n}^{2} + \sigma _{j}^{2}}}} \right\},$
где Ps – мощность сигнала, $\sigma _{n}^{2},$ $\sigma _{j}^{2}$ – мощности шума и ответной помехи соответственно. После логарифмирования и соответствующих преобразований из (49) получаем

(50)
$\frac{{\sigma _{j}^{2}}}{{\sigma _{n}^{2}}} = \frac{{\ln 2{{P}_{{E0}}}}}{{\ln 2{{P}_{{E1}}}}} - 1.$

Полученное соотношение показывает, во сколько раз мощность ответной помехи должна превосходить мощность аддитивного гауссовского шума, чтобы обеспечить заданные вероятности PE0 и PE1. Из (48) также видно, что $\sigma _{j}^{2}$ → ∞ при PE1 → 1/2.

Пример. Рассмотрим систему связи с ППРЧ в выделенной полосе 30 МГц, ширина полосы сигнала 3.75 кГц, расстояние между передатчиком и приемником 10 км, вероятность ошибки на бит в РТС до подавления PE0 = 10−5, что соответствует входному ОСШ 13.35 дБ. Пренебрегая затуханием сигнала при распространении, считаем, что и на входе обнаружителя СОП отношение сигнал-шум также 13.35 дБ, заданы вероятности обнаружения и ложной тревоги Pd = 0.99, Pf = 10−5. В соответствии с (28) время накопления энергии для обнаружения с заданными вероятностями для N = 1 равно τ02 = 2.21 × 10−4 c, для N = 5 аналогичное время равно τ02 = 3.38 × 10−5 c, и для N = 7 τ02 = = 2.33 × 10−5 c.

Аналоговая часть обнаружителя пеленгатора состоит из каскадов входных полосовых фильтров второго порядка с полосой согласованной с шириной спектра обнаруживаемого сигнала. В этом случае в соответствии с (36), например, для 4-каскадного фильтра имеем задержку сигнала t0 = 3.4 × 10–4 с. Суммарная задержка составит ${{t}_{0}} + {{\tau }_{{02}}} = 3.4 \times {{10}^{{ - 4}}}$ + $2.26 \times {{10}^{{ - 4}}}$ = $5.66 \times {{10}^{{ - 4}}}{\text{с}}$ для N = 1 и ${{t}_{0}} + {{\tau }_{{02}}} = 3.4 \times {{10}^{{ - 4}}}$ + $3.38 \times {{10}^{{ - 5}}}$ = $ = 3.73 \times {{10}^{{ - 4}}}\,\,{\text{с}}$ для N = 5. Если в аналоговой части применяются фильтры с шириной полосы, значительно превосходящих полосу сигнала, то задержка главным образом определяется временем накопления энергии для обнаружения с заданными вероятностями.

На рис. 6а представлено семейство кривых, определяющих верхнюю границу доступных для подавления числа скачков частоты (в сек), в зависимости от PE1 при допустимой вероятности Pдоп = 0.05 и N = 1 для различных значений эксцентриситета ε: 0.1 (кривая 1), 0.5 (кривая 2), 0.99 (кривая 3). Аналогичные кривые для N = 5 с учетом задержки в аналоговых цепях приведены на рис. 6б. При построении кривых на рис. 6в и 6г не учитывались задержки в аналоговых входных цепях, а задержка определялась только необходимым временем накопления энергии для обнаружения с заданными вероятностями. Рисунок 6в соответствует N = 1, а рис. 6г построен для N = 5.

Рис. 6.

Зависимость скорости переключения ППРЧ от вероятности ошибки после подавления.

Кривые на рис. 6в и 6г соответствуют предельным границам необходимой скорости переключения ППРЧ.

Таким образом, СОП, находясь за пределами эллипса с заданным эксцентриситетом, при указанных скоростях переключения ППРЧ не сможет подавить РТС до требуемой вероятности ошибки. Случай ε = 0.99 соответствует ситуации, когда СОП практически может находится на отрезке прямой между передатчиком и приемником РТС. В этом случае требуется максимальная скорость переключения. Дальнейшее уменьшение $\varepsilon $ увеличивает границы области, запрещенной для нахождения СОП. В связи с этим необходимая скорость переключения ППРЧ уменьшается.

На рис. 7 изображен график относительной мощности ответной помехи, обеспечивающей заданные вероятности PE0 и PE1. Как видно из графика, с ростом PE1 (т.е., с увеличением мощности ответных помех) при Pдоп = 0.05 граница для скачков частоты резко возрастает, и, начиная с PE1 = 0.1 (что соответствует увеличению мощности помех на 7.5 дБ), далее начинается незначительное плавное увеличение.

Рис. 7.

Зависимость относительной мощности помех от вероятности подавления.

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

ВЫВОДЫ

1. Построена вероятностно-временная модель оценки помехозащищенности систем радиосвязи в условиях радиоэлектронной борьбы.

2. На основе модели получена замкнутая система уравнений, позволяющая оценить предельные возможности станции ответных помех при подавлении систем радиосвязи с ППРЧ.

3. Станция ответных помех способна подавлять системы радиосвязи с ППРЧ до определенной скорости переключения рабочей частоты. При этом граница скорости переключения зависит не только от временных параметров (таких как время геометрической разности хода электромагнитной волны, время задержки сигнала в фильтрах селекции, время интегрирования в обнаружителе), но и от энергетического параметра (мощности ответной помехи).

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

  1. Типугин В.Н., Вейцель В.А. Радиоуправление. М.: Сов. радио, 1962.

  2. Борисов В.И. Помехозащищенность систем радиосвязи. Основы, теории и принципы реализации. М.: Наука, 2009.

  3. Борисов В.И., Зинчук В.М. Помехозащищенность систем радиосвязи. Вероятностно-временной подход. М.: Радио и связь, 1999.

  4. Борисов В.И., Зинчук В.М., Лимарев А.Е. и др. Помехозащищенность систем радиосвязи с расширением спектра сигналов методом псевдослучайной перестройки рабочей частоты. М.: Радио и связь, 2000.

  5. Борисов В.И., Зинчук В.М., Лимарев А.Е., Шестопалов В.И. Помехозащищенность систем радиосвязи с расширением спектра сигналов прямой модуляцией псевдослучайной последовательностью. М.: РадиоСофт, 2011.

  6. Артемов М.Л., Афанасьев О.В., Сличенко М.П. // Радиотехника. 2014. № 11. С. 11.

  7. Артемов М.Л., Сличенко М.П. // Антенны. 2018. № 5. С. 31.

  8. Артемова Е.С., Сличенко М.П. // Антенны. 2018. № 5. С. 47.

  9. Артемов М.Л., Борисов С.Г., Сличенко М.П. // Радиотехника. 2014. № 11. С. 11.

  10. Cличeнкo M.П. // PЭ. 2014. T. 59. № 5. C. 473.

  11. Справочник по специальным функциям c формулами, графиками и математическими таблицами / Под ред. Абрамовица М., Стиган И. М.: Наука, 1979.

  12. Джонсон Д., Джонсон Дж., Мур Г. Справочник по активным фильтрам. М.: Энергоатомиздат, 1983.

  13. Torrieri D.J. // IEEE J. on Selected Areas in Commun. 1989. V. 7. № 4. P. 569.

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