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

Системы управления,связи и безопасности №1. 2019

Systems of Control, Communication and Security ISSN 2410-9916

УДК 621.396.663

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

Дворников С. В., Конюховский В. С., Симонов А. Н.

Постановка задачи: увеличение количества радиоизлучений приводит к интерференции попадающих в полосу пропускания радиоприемного устройства сигналов различных источников. Устойчивые к интерференции известные методы пеленгования используют структурно-статистическую избыточность для пространственного разделения источников радиоизлучений (ИРИ) и оценивания пеленгов. Это затрудняет их применение в условиях массогабаритных и временных ограничений. Целью работы является обеспечение интерференционной устойчивости узкобазисных малоэлементных амплитудных пеленгаторов за счет присущей им пространственной избирательности. Предлагается, не применяя частотной селекции, использовать спектральные различия радиоизлучений, выделяя пространственные параметры из интерференционной смеси сигналов в приемных каналах пеленгатора. Используемые методы. Решение задачи автоматического амплитудного пеленгования интерферирующих радиоизлучений основано на отыскании огибающей семейства кривых, в качестве которых используется годограф суммарного поля. Для формализации задачи измерения пеленгов рассмотрены два альтернативных подхода: на основе метода главных компонент и на основе кластерного анализа. Новизна представленного решения заключается в использовании гипотезы о том, что годограф суммарного поля касается огибающего его параллелограмма в точках, наиболее удаленных от центра изображения на индикаторе. Также к элементам новизны следует отнести использование скачка производной от функции, описывающей стороны параллелограмма, в качестве признака для нахождения его вершин. Результат: использование представленного решения по автоматическому пеленгованию интерферирующих радиоизлучений позволяет применить узкобазисные малоэлементные амплитудные пеленгаторы в сложной сигнально-помеховой обстановке. Проведенное моделирование показало возможность обнаружения многосигнальной ситуации и измерения пеленга на два ИРИ, при этом удается добиться углового разрешения в один градус при отношении сигнал/шум 10 дБ за время, соответствующее одному периоду колебаний разностной частоты. Практическая значимость: представленное решение предлагается реализовать в средствах радиомониторинга, функционирующих в насыщенной сигнально-помеховой обстановке, когда налагаются дополнительные ограничения на их массогабаритные и временные показатели.

Ключевые слова: амплитудное пеленгование, автоматическое пеленгование, радиопеленгование, интерференция радиоизлучений.

Введение

Пеленгование источников радиоизлучений является одной из важнейших составляющих радиомониторинга [1-5] наряду с обнаружением, измерением параметров и распознаванием радиосигналов, перехватом сообщений и определением местоположения ИРИ. От эффективности пеленгования во многом зависят точность и достоверность добываемой системами радиомониторинга информации. Широкое применение кодирования, маскирования и шифрования

Библиографическая ссылка на статью:

Дворников C. В., Конюховский В. С., Симонов А. Н. Автоматическое амплитудное пеленгование интерферирующих радиоизлучений // Системы управления, связи и безопасности. 2019. № 1. С. 67-81. DOI: 1024411/2410-9916-2019-10104. Reference for citation:

Dvornikov S. V., Konyukhovsky V. S., Simonov A. N. Automatic Amplitude Direction Finding of Interfering Signals. Systems of Control, Communication and Security, 2019, no. 1, pp. 67-81. DOI: 1024411/2410-9916-2019-10104 (in Russian).

DOI: 1024411/2410-9916-2019-10104

Системы управления,связи и безопасности №1. 2019

Systems of Control, Communication and Security ISSN 2410-9916

передаваемой информации приводит к тому, что нередко координатная и пеленговая информация остаются единственными достоверными и надежными сведениями об источниках радиомониторинга [6].

Функционирование средств пеленгования в настоящее время происходит в условиях сложной сигнально-помеховой обстановки, связанной с лавинообразным увеличением количества излучающих радиоэлектронных средств и преднамеренных помех [7, 8]. Это приводит к тому, что в полосу пропускания радиоприемного устройства пеленгатора попадает несколько сигналов от различных источников радиоизлучений, приводя к их интерференции. В средствах радиосвязи в этой ситуации возникает задача помехозащиты, когда необходимо выделить полезный сигнал абонента и исключить мешающие сигналы от других ИРИ. В средствах радиомониторинга при пеленговании, зачастую, все принимаемые радиоизлучения являются полезными, то есть получены от источников, представляющих интерес с точки зрения добывания информации. Тогда возникает задача одновременного пеленгования нескольких интерферирующих радиоизлучений.

Известно много эффективных методов пеленгования, которые достаточно хорошо изучены и успешно реализованы на практике [1-3, 9, 10]. Исторически первыми были амплитудные пеленгаторы от рамочных поворотных до широкобазисных стационарных. С развитием радиоэлектронных технологий, в частности с появлением высокоточных, высокостабильных опорных генераторов, более распространенными стали фазовые пеленгаторы. Частотные пеленгаторы основаны на использовании эффекта Доплера, возникающего при движении приемного антенного элемента относительно источника радиоизлучения. Кроме классических пеленгаторов в современных средствах радиомониторинга все чаще используют корреляционные методы пеленгования, опирающиеся на использование многоэлементных антенных систем и многоканальных приемных трактов. Следует также обратить внимание на перспективные собственно-структурные методы, обладающие повышенным пространственным разрешением, основанные на использовании статистических свойств радиосигналов и шумов.

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

DOI: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

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

Постановка задачи

Одним из наиболее известных является классический амплитудный пеленгатор [12], построенный по схеме Эдкока (рис. 1).

Рис. 1. Схема узкобазисного амплитудного пеленгатора

Информацию о направлении на ИРИ получают, используя амплитудные уравнения двух ортогональных узко базисных ориентированных по сторонам света пеленгацнонных пар антенн (ППА), амплитуды разностных ЭДС которых можно представить в виде:

2л, ^

sin В cos [3,

:Дзв =2£sin — ¿>sin6cosß

V А,

К.,

. (2п »in -

£дсю = 2Е sin|^—ôcos 0COS ßj = Еп

cosöcosß,

где £Дзв, £Дсю - амплитуды разностных ЭДС в ППА «запад - восток» и «север - юг» соответственно; Е - амплитуда ЭДС в антенном элементе; Ь - половина базы антенной системы; Етях = АЕЬ я/А, - максимальное значение амплитуд разностных ЭДС; 8, |3 - азимут и угол места на ИРИ соответственно.

Из сравнения выражений (1) и (2) получают: 6 = аг аё(ЕА„/ЕАсю).

В описанном пеленгаторе сигналы с каждого антенного элемента (после усиления в УКс, УКю, УКз, УКв) объединяются на разность, а в качестве индикатора пеленга используется электронно-лучевая трубка. Разностные ЭДС с ППА «север - юг» и «запад - восток» подаются соответственно на вертикальные и горизонтальные отклоняющие пластины. Электронно-лучевая трубка (ЭЛТ) помимо индикации выполняет роль устройства для определения функции агс^, поэтому на экране положение светящейся линии соответствует азимуту на ИРИ.

DPI: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

В [9, 10] доказана возможность и показан способ одновременного пеленгования нескольких ИРИ, частоты которых находятся в полосе пропускания приемных каналов описанного пеленгатора.

Если в полосу пропускания (2Дсоп) приемных каналов пеленгатора (см. рпс. 1) попадают два синусоидальных колебания, соизмеримых по амплитуде, с

со2 — со.

частотами со, и со2, причем со2 со, = Í 2 < 2Дсоп и J—--1« 1, где

со

со = 0,5(со, + со2), тогда суммарные напряжения на выходе вертикального и горизонтального каналов будут равны:

ив (í) = UBl sin го,/ + Ub2 sin co2f, wr (í) = UTl sin со,/ + Uv2 sin co2/, где UB] = [/, eos 9,, UTl = Ul sin 8,, Ub2 = U2 cosG2, Ur2 = U2 sin92, 6, Ф 92.

Напряжения двух сигналов в каждом канале могут быть представлены как биения в виде колебаний одной частоты со,, но с переменными во времени (с частотой Í2) амплитудой и фазой, а отклонения луча на экране ЭЛТ по осям Оу и Ох описываются выражениями

у = L£/B(/)sin со,/, х = L£/r(/)sin(co,/ - Ф(/)),

где UB(t) = ^U2Bl + U2b2 + 2 UBlUB2 eos Q/ ; Ur (t) = ^ + U;2 + 2UrlUr2 eos Q/ ;

^./ч ^/ч ^/ч sin Q / A , . U _ sin Q /

Ф(/) = Фг(/)-Фв(/); Фв(/) = агс1ё---— ; Фг(/) = агс1ё й

UBl + Ub2 cos П/ " UY] + Ui2 cos Q/

L - чувствительность ЭЛТ.

Уравнение движения луча в координатах хОу будет:

--2———-со8Ф

ÚUl{t) ÜU г (/ )UB (/) l}U;it) U

Получено уравнение спирали с непрерывно изменяющимися во времени конфигурацией и наклоном звеньев. Каждое звено спирали предетавляет собой фигуру, близкую к эллипсу, частота изменения формы и наклона большой оси которого равна Q = co2-co, (Í2 = 2nF , F = /2 - /,). Огибающая спирали будет представлять собой параллелограмм со сторонами - линиями пеленгов на первый и второй источники радиоизлучений (рис. 2 а, б).

Доказательство того, что линия (3) вписана в параллелограмм, в обшем виде требует совместного решения систем уравнений, получаемых из (3):

y(jt, Q/) = 0 и

d

x(y,Dt) = 0 и av(v-£1/) 0.

a

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

DPI: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

индикатор ЭЛТ приводит к усложнению изображения на экране, при трех ИРИ изображение имеет вид параллелепипеда, представленного в аксонометрической проекции, причем по трем его ребрам определяются азимуты всех трех пеленгуемых источников (см. рис. 2 в).

а б в

Рис. 2. Изображение на экране индикатора: а - поворотный визир для отсчета пеленгов; б - изображение при двух ИРИ; в - изображение при трех ИРИ

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

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

Идея, положенная в основу метода, состоит в том, что каждая точка т, в которой годограф поля касается описывающего параллелограмма АВСО, является наиболее удаленной от центра по сравнению с другими точками, имеющими тот же полярный угол ф (рис. 3).

Для расчета азимутов теоретически достаточно для каждого ИРИ выбрать пару точек, расположенных на соответствующей стороне параллелограмма, и оценить угол наклона прямой, проходящей через них. Например, для нахождения азимута на первый ИРИ можно выбрать соответствующую ему сторону ВС

Р01: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

(или DA) (см. рис. 3), для которой, выделив две произвольные точки mi = (r¡,(р.) п т. - (г.,фу), рассчитать угол ее наклона

r¡ sin ф  - f'j sin ф7

9, = arctg-

t) СОБф. - г - cos ф -

Рис. 3. Геометрическая основа метода

В современных средствах пеленгования, построенных на аппаратно-программной платформе, все преобразования сигналов выполняются путем цифровой обработки отсчетов, полученных после дискретизации сигналов. Поэтому положение точки т с декартовыми (х, у) или полярными (г, ф) координатами при движении по отображаемой спирали будет также меняться дискретно, ее временная выборка будет представлять собой массив из N цифровых дискретных отсчетов М = {т.}, г = 1, N, а изображение на индикаторе можно представить совокупностью N точек т1 = (х., у1) (рис. 4 а).

y¡

в »V

» •

А

• т,

ТП, Й • <

-7

* »

О

im

д: —»•

•• V *. D

i-

л' • и 1 * *. .*.v: .% •

».. i

A i" "

• • J

• " L

I 1

Г

I •

—TL

Фл

фо Q фс ф/ ф,

Фй 7t

ф

а б

Рис. 4. Графическое представление массива данных: а - для декартовых координат; б - для полярных координат

После преобразования в полярные координаты (г, ф) изображение полученного массива М изменится (см. рис. 4 б), однако на нем по прежнему можно выделить положение соответствующих точек т и сторон параллелограмма

DOI: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

АВСИ, при этом поиск максимального удаления гтах для заданного ф становится процедурой обработки полученного массива.

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

Дискретность отсчетов приводит к тому, что выбранные точки массива М, соответствующие максимальному удалению гтах, смещены от сторон параллелограмма АВСО (рис. 5). По этой причине азимуты на ИРИ, рассчитанные по формуле (4), будут иметь методическую ошибку.

.........'

A i-"-

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

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

Рис. 6. Устранение методической ошибки: а - увеличение выборки сигналов; б - принцип гистограммной обработки

Р01: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

Результатом указанной п гстограммной обработки будет массив Мт = {тк}, к = \,К, К « N, содержащий координаты отобранных точек т* Ф*) (Р"с. 7).

......

» » ' «С

в. .

......

D

а : .

1111111 I И I I I I 1111111111 II1111

11111 111111111111 И I I I I 111111 I—к

о

Дф

Дф

а б

Рис. 7. Результаты гистограммной обработки: а - для декартовых координат; б - для полярных координат.

Далее необходимо осуществить поиск вершин параллелограмма. Для этого предлагается воспользоваться признаком скачка производной от функции г(ф) на границах отрезков А В, ВС, СО и О А, т. е критерием того, что при аргументах <рА, ф3, фс н фс (см. рис. 4 б) производная г'(ф) резко изменяется.

Такое изменение соответствует скачкообразному увеличению в вершинах параллелограмма абсолютного значения второй производной г по ф :

г"(ф)>С, фе(фл,ф3,фс,фд),

г"(ф)<С, ф^(фл,ф3,фс,фд),

где С - значение порога, определяемое эмпирически.

При цифровой обработке производные следует записать в дискретной форме:

>-;(ф)=

гДФ) =

гк - Гк-1

ФА +1 - ФА Фа-ФА-1

Результатом поиска по критерию (5) будут граничные значения фк, соответствующие вершинам параллелограмма АВСИ, после чего имеется возможность разделить все точки (гк,ф^) на группы М1 = {т^ }, / = 1.4, к1 = 1, К1, где А', - количество точек массива МТ для сторон АВ, ВС, СО, ОА соответственно.

Для расчета азимутов на каждый ИРИ предлагается использовать метод главных компонент [13], для чего все точки пгк (гк,ф^ ) массива МТ преобразуются в декартовы координаты пгк = (л^, ук ): ^ =Г,С05ф,, у, =^шфг

Р01: 1024411/2410-9916-2019-10104

Системы управления, связи и безопасности_№1. 2019

Systems of Control, Communication and Security ISSN 2410-9916

Разброс координат (xk, yk ) относительно каждой стороны параллелограмма приблизительно описывается эллипсом рассеяния [13], угловое положение большой оси которого равно углу наклона соответствующей стороны. Чтобы найти эти параметры необходимо для каждой группы точек М1 составить ковариационную матрицу:

R f Dxi СО\{х,У}^

1^соу{у,х}г Dyt У где Dxj, Dyt - дисперсии, a cov{x, y}f, cov{y,Jc}( - ко вар нации, соответственно, координат (х, у) для каждой стороны I.

Далее вычисляются собственные значения {А^Д^ и собственные векторы {v,,v2}/ матриц R,. Будем считать, что собственные значения в {/„,

расположены в порядке убывания. Информацию об азимуте на каждый ИРИ извлекают из углового положения / - й стороны, которую получают, используя собственные векторы, соответствующие большему собственному значению:

01=arctg^, / = 1,3; vn,

9., = arctg1^5-, 1 = 2,4,

где vn = {xn,ynf - первый собственный вектор матрицы R, ; х l,yl - координаты собственного вектора. Индекс "1" в vn указывает на собственный вектор,

соответствующий максимальному собственному значению, параметр / устанавливает выбор соответствующей стороны параллелограмма.

Кроме метода главных компонент можно предложить альтернативный метод обработки массива Afr, основанный на разбиении множества точек тк ={хк,ук) на кластеры. Для этого, используя множество [тк}, сформируем

множество векторов А = {ак}, ак = щтк+1, к = 1,К—1, такое, что

а*

( Y -X ^ Ai+1 Лк

\Уш~Ук

к /

Векторы аА можно рассматривать как векторы производных по сторонам параллелограмма, которые образуют с осью Ох углы

Лу*

ук = arctg

Axt

Массив {уА} можно разбить на четыре кластера, соответствующие сторонам параллелограмма (рис. 8). Для этого сначала находят минимальное утш и максимальное утах значения в массиве и рассчитывают среднее Чо»г +7пип) между ними. Значение у а¥Г позволяет разделить массив

{уА} на два кластера. Каждый кластер задается множеством индексов элементов, содержащихся в нем. Аналогично производится разбиение каждого из по-

DOI: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

лученных кластеров. В результате индексы к[ = 1, К,, I- 1,4, соответствующие сторонам параллелограмма, попадают в четыре различных кластера.

Ук' L .

DA -л 1 1 1 1 1 О |Фс фл 71 Ф

1 1 1 1 1 1 фо 1 1 1 1 * 1 1 1 !

! >» 1 1 1 1 1 ВС 1 \ К ! * ч 1 1 1 itt и«****

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

Для расчета азимутов на каждый из двух ИРИ следует усреднить значения ук, сгруппированные по параллельным сторонам:

(

9, = 0,5

^_

К,

+

±3_

К,

е, = о,5

ХУ* Хуа

h

К,

+

Кл

Результаты моделирования

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

°есР >

где ст0ср = л/0,5(ад[ + ); Д0 = |0,-02|; а9], а62 и 0,, ()2 -соответственно СКО

пеленгования и азимуты первого и второго ИРИ.

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

а СКО пеленгования находить, нормируя его относительно большего собственного значения X,

"12

= arctg1

"12

"Il

где G0/ - СКО пеленгования.

Р01: 1024411/2410-9916-2019-10104

Systems of Control, Communication and Security

ISSN 2410-9916

Результаты моделирования для двух вариантов длительности т выборок сигналов (один и десять периодов разностной частоты Г) при различных отношениях сигнал/шум (ОСШ) д представлены на рис. 9.

А о0ср, град j аеср, град

if ъ t 1 IF Ъ 1 1 1 1 гчг

ъ i ъ

ъ 1 ч ч

ч t / а ч ъ /б ч

V 1 ч ч

ч

8 * ! 1 т / в

% / /

а ' »

4

Рис. 9. Зависимость разрешающей способности от ОСШ: а - оценка по методу главных компонент; б - оценка по методу кластерного анализа; в - оценка на основе отношения собственных значений

Анализ полученных зависимостей (рис. 9) показал следующее.

1. Более высокой разрешающей способностью обладает метод главных компонент. Этот результат является ожидаемым, так как данный метод основан на собственных векторах ковариационных матриц К,, то есть наилучшим образом использует все имеющиеся в сигналах корреляционные свойства.

2. Для метода главных компонент отмечено парадоксальное превышение точности пеленгования в условиях интерференции (рис. 9) по сравнению с од-посишальной ситуацией [12]. Объяснение этого кроется в том, что даже выборка длиной всего в один период разностной частоты (т = 1//7) все же содержит определенный набор статистических данных. Использование этого набора и приводит к повышению точности по сравнению с оценкой для односигнальной обстановки, полученной по мгновенным отсчетам.

3. Применение метода кластерного анализа возможно только при значительных ОСШ (д>20 дБ). Причиной этого, предположительно, является неполное использование статистических связей рассчитанных значений ук, так как при формировании массива {уА} используются только соседние отсчеты выборки {аА.}.

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

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

Р01: 1024411/2410-9916-2019-10104

Системы управления,связи и безопасности №1. 2019

Systems of Control, Communication and Security ISSN 2410-9916

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

Выводы

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

Литература

1. Рембовский А. М., Ашихмин А. В., Козьмин В. А. Радиомониторинг -задачи, методы, средства / под ред. А. М. Рембовского. - М.: Горячая линия -Телеком, 2010. - 624 с.

2. Липатников В. А., Царик О. В. Методы радиоконтроля. Теория и практика. Монография. - СПб.: ГНИИ «Нацразвитие», 2018. - 608 с.

3. Poisel R. A. Electronic Warfare Target Location Methods. - Artech House, Norwood MA, 2005. - 272 p.

4. Радзиевский В. Г., Сирота А. А. Теоретические основы радиоэлектронной разведки. - М.: Радиотехника, 2004. - 432 с.

5. Макаренко С. И. Информационное противоборство и радиоэлектронная борьба в сетецентрических войнах начала XXI века. Монография. - СПб.: Наукоемкие технологии, 2017. - 546 с.

6. Кудрявцев А. М., Смирнов А. А., Федянин А. В. Алгоритм "трассовой" обработки данных радиомониторинга // Научно-технические ведомости Санкт-Петербургского государственного политехнического университета. Информатика. Телекоммуникации. Управление. 2010. № 1 (93). С. 38-42.

7. Челышев В. Д., Якимовец В. В. Зарубежные радиоэлектронные системы наземного и спутникового мобильного радиосервиса. - СПб.: ВАС, 2012. - 388 с.

8. Агиевич С. Н., Волощук Э. В., Круглов С. А., Луценко С. А., Пономарев А. А. Способ радиоподавления спутниковых каналов управления // Патент на изобретение RUS 2677261, опубл. 16.01.2019, бюл. № 2.

9. Мезин В. К. Автоматические радиопеленгаторы. - М.: Сов. радио, 1969. - 216 с.

10. Кукес И. С., Старик М. Е. Основы радиопеленгации. - М.: Сов. радио, 1964. - 640 с.

11. Комарович В. Ф., Никитченко В. В. Методы пространственной обработки радиосигналов. - Л.: ВАС, 1989. - 278 с.

12. Симонов А. Н., Волков Р. В., Дворников С. В. Основы построения и функционирования угломерных систем координатометрии источников радиоизлучений: учебное пособие / под ред. А. Н. Симонова. - СПб.: ВАС, 2017. - 248 с.

13. Налимов В. В. Теория эксперимента. - М.: Наука, 1971. - 208 с.

DOI: 1024411/2410-9916-2019-10104

Системы управления,связи и безопасности №1. 2019

Systems of Control, Communication and Security ISSN 2410-9916

References

1. Rembovsky A. M., Ashikhmin A. V., Kozmin В. А. Radiomonitoring -zadachi, metody, sredstva [Radiomonitoring - problems, methods, tools]. Moscow, Goryachaya Liniya - Telekom Publ., 2010. 624 p. (in Russian).

2. Lipatnikov V. A., TSarik O. V. Metody radiokontrolya. Teoriya i praktika. Monografiia [Methods of radiomonitoring. Theory and practice. Monograph]. Saint Petersburg, Humanitarian National Research Institute «Natsrazvitie» Publ., 2018. 608 p. (in Russian).

3. Poisel R. A. Electronic warfare target location methods. Artech House, Norwood MA, 2005. 272 p.

4. Radzievsky V. G., Sirota A. A. Teoreticheskie osnovy radioehlektronnoj razvedki [Theoretical basis of radioelectronic reconnaissance]. Moscow, Radiotekhnika Publ., 2004. 432 p. (in Russian).

5. Makarenko S. I. Informatsionnoe protivoborstvo i radioelektronnaia borba v setetsentricheskikh voinakh nachala XXI veka. Monografiia [Information warfare and electronic warfare to network-centric wars of the early XXI century. Monograph]. Saint Petersburg, Naukoemkie Tekhnologii Publ., 2017. 546 p. (in Russian).

6. Kudryavtsev A. M., Smirnov A. A., Fedyanin A. V. Algorithm of "route" processing of radiomonitoring data. Nauchno-tekhnicheskie vedomosti Sankt-Peterburgskogo gosudarstvennogo politekhnicheskogo universiteta. Informatika. Telekommunikatsii. Upravlenie, 2010, no. 1 (93), pp. 38-42 (in Russian).

7. Chelyshev V. D., Yakimovets V. V. Zarubezhnye radioehlektronnye sistemy nazemnogo i sputnikovogo mobil'nogo radioservisa [Foreign ground and satellite mobile radioservice radioelectronic systems]. Saint Petersburg, Military Telecommunications Academy Publ., 2012. 388 p. (in Russian).

8. Agievich S. N., Voloshchuk E. V., Kruglov S. A., Lutsenko S. A., Ponomarev A. A. Method of radiosuppression of satellite control channels. Patent Russia, no. RUS 2677261, Publish. 16.01.2019, bul. no. 2 (in Russian).

9. Mezin V. K. Avtomaticheskie radiopelengatory [Automatic direction finder]. Moscow, Sovetskoe Radio Publ., 1969. 216 p. (in Russian).

10. Kukes I. S., Starik M. E. Osnovy radiopelengatsii [Basics of radio direction finding]. Moscow, Sovetskoe Radio Publ., 1964. 640 p. (in Russian).

11. Komarovich V. F., Nikitchenko V. V. Metody prostranstvennoj obrabotki radiosignalov [Methods of spatial processing of radio signals]. Leningrad, Military Telecommunications Academy Publ., 1989. 278 p. (in Russian).

12. Simonov A. N., Volkov R. V., Dvornikov S. V. Osnovy postroeniya i funktsionirovaniya uglomernykh sistem koordinatometrii istochnikov radioizluchenij [Basics of construction and functioning of triangulation emitters location systems]. Saint Petersburg, Military Telecommunications Academy Publ., 2017. 248 p. (in Russian).

13. Nalimov V. V. Teoriya ehksperimenta [Theory of experiment]. Moscow, Nauka Publ., 1971. 208 p. (in Russian).

Статья поступила 5 ноября 2018 г.

DOI: 1024411/2410-9916-2019-10104

Системы управления,связи и безопасности №1. 2019

Systems of Control, Communication and Security ISSN 2410-9916

Информация об авторах

Дворников Сергей Викторович - доктор технических наук, профессор. Профессор кафедры радиосвязи. Военная академия связи имени Маршала Советского Союза С.М. Буденного. Область научных интересов: радиомониторинг, формирование и анализ сигналов, координатометрия источников радиоизлучений, помехозащищенность радиолиний. E-mail: practicdsv@yandex.ru

Конюховский Владимир Самуилович - кандидат физико-математических наук, профессор. Профессор кафедры математики. Военная академия связи имени Маршала Советского Союза С.М. Буденного. Область научных интересов: аналитическая геометрия, векторная алгебра, цифровая обработка сигналов и изображений, компьютерное моделирование. E-mail: konyukhovsky@yandex.ru

Симонов Алексей Николаевич - кандидат технических наук, доцент. Докторант. Военная академия связи имени Маршала Советского Союза С.М. Буденного. Область научных интересов: радиомониторинг, анализ и обработка сигналов, пеленгование и координатометрия источников радиоизлучений. E-mail: sanmailbox@yandex.ru

Адрес: 194064, Россия, г. Санкт-Петербург, Тихорецкий пр., д. 3.

Automatic Amplitude Direction Finding of Interfering Signals

S. V. Dvornikov, V. S. Konyukhovsky, A. N. Simonov

Problem statement. Increasing number of radio emissions leads to interference of signals of different sources falling into the bandwidth of the radio receiver. Sustainable to the interference of the known methods of direction finding use structural and statistical redundancy for spatial separation of emitters and the estimation of the bearing. This makes it difficult to use them in conditions of mass and time constraints. The aim of the work is to ensure the interference stability of narrow-band low-element amplitude direction finders due to their inherent spatial variability. It is proposed to use the spectral differences of radio emissions, highlighting the spatial parameters of the interference mixture of signals in the receiving channels of the direction finder, without using frequency selection. Methods. The solution of the problem of automatic amplitude bearing of interfering radio emissions is based on finding the envelope of the family of curves, as which the total field hodograph is used. Two alternative approaches are proposed for formalization of bearing measurement: on the basis of principal components method and on the basis of cluster analysis. The novelty of the presented solution lies in the use of the idea that the hodograph of the total field concerns the envelope of its parallelogram at the points most distant from the center of the image on the indicator. Also, the elements of novelty include the use of a jump derived from the function describing the sides of the parallelogram, as a sign for finding its vertices. Result. The use of the proposed solution for automatic direction finding of interfering radio emissions makes it possible to apply narrow-band small-element amplitude direction finders in a complex signal-noise environment. The simulation showed the possibility of detecting a multisignal situation and measuring the bearing on two emitters, while it is possible to achieve an angular resolution of one degree at a signal/noise ratio of 10 dB in a time corresponding to one period of oscillations of the difference frequency. Practical relevance. The presented solution is suggested to use in radiomonitoring equipment, operating in the saturated signal-to-interference about the situation when you have additional restrictions on their dimensions and weight.

Key words: amplitude direction finding, automatic direction finding, radio direction finding, interference of radio signals.

DOI: 1024411/2410-9916-2019-10104

Системы управления,связи и безопасности №1. 2019

Systems of Control, Communication and Security ISSN 2410-9916

Information about Authors

Sergey Viktorovich Dvornikov - Dr. habil. of Engineering Sciences, Professor. Professor at the Department of radio communications. Military Telecommunications Academy named after Marshal of the Soviet Union S. Budyonny. Field of research: the radiomonitoring, formation and analysis of signals, emitters location, noise immunity of radio links. E-mail: practicdsv@yandex.ru

Vladimir Samuilovich Konyukhovsky - Ph.D. of Physico-mathematical Sciences, Professor. Professor at the Department of mathematics. Military Telecommunications Academy named after Marshal of the Soviet Union S. Budyonny. Field of research: analytical geometry, vector algebra, digital signal and image processing, computer simulation. E-mail: konyukhovsky@yandex.ru

Alexey Nikolaevich Simonov - Ph.D. of Engineering Sciences, Associate Professor. Doctoral Candidate. Military Telecommunications Academy named after Marshal of the Soviet Union S. Budyonny. Field of research: radiomonitoring, signal analysis and processing, direction finding and location of emitters. E-mail: sanmail-box@yandex.ru

Address: Russia, 194064, Saint-Petersburg, Tihoreckiy avv., 3.

DOI: 1024411/2410-9916-2019-10104