Определение координат источника радиоизлучения Текст научной статьи по специальности «Физика»

УДК 621.372.542

Определение координат источника радиоизлучения

© А. А. Грешилов МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

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

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

Введение. Пеленгация источников радиоизлучения (ИРИ) имеет место в процессе мониторинга радиоэлектронной обстановки. При этом необходимо определять азимутальные, угломестные пеленги ИРИ и амплитуду сигнала, по которым устанавливают координаты источника радиоизлучения. Пеленгатор регистрирует излучение путем записи сигналов на элементы (вибраторы) антенной системы (АС). Выполняя различные действия над сигналами, принимаемыми вибраторами, определяют параметры излучения, а затем координаты источника радиоизлучения. Проблема состоит в том, чтобы понять, какие действия над сигналами окажутся эффективнее. Данной проблеме посвящены работы [1, 2]. В частности, наиболее полно проблема определения координат источника излучения (ИРИ) рассмотрена в [2]. Прием сигналов ИРИ осуществляется в заданной полосе частот AF перемещающимся в пространстве бортовым пеленгатором на летно-подъемном средстве (ЛПС). Измеряются пространственно-информационные параметры обнаруженных сигналов: азимута 9г. и угла места Рг. в системе координат антенной системы с одновременным определением местоположения ЛПС {£ipS, LipS, HipS}, где Др8, LipS и Hp — широта, долгота и высота ЛПС соответственно. Предварительно устанавливают

удаления ИРИ d. от ЛПС и координаты ИРИ Vп. ={X0,Y0,Z0}. в

момент времени t. в левосторонней системе декартовых координат антенной системы пеленгатора. Корректируют координаты ИРИ V у. = {X'0, Y0', Z 0}. с учетом априорно известной ориентации антенной системы бортового пеленгатора относительно ЛПС {kmt, lmt, С ant} путем последовательного умножения значений координат Vпi на соответствующие углам Эйлера матрицы поворота. После чего вычисляют истинные геоцентрические координаты местоположения ИРИ

Vгц. ={X0, Y", Z0|г. с учетом измеренных в момент времени t. пространственных углов ЛПС (крена k ips тангажа /ips курсового угла a lps. и склонения С lps. ), а также координат его местоположения (широты Blps долготы Llps и высоты Hlps .). Склонение Сlps. рассчитывают как разность между путевым |ilps . и курсовым a lps . углами ЛПС. Истинные геоцентрические координаты Vгц. местоположения ИРИ преобразуют в географические координаты Vг. = {B0, L 0, H0}..

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

Этот способ определения пеленгов имеет недостатки:

1) не обеспечивает высокую точность результатов, используя заранее подготовленную таблицу, которая связывает показания на элементах АС с азимутальными и угломестными пеленгами, полученными с помощью внешнего генератора, включая облет ЛПС;

2) степень отличия измеренных параметров Дф l,h изм (fV) от эталонных значений Дф l,h эт (fv), рассчитанных для всех направлений прихода сигнала Д9 к и Ар с и всех fV, вычисляют по формуле [2]

N N / \2

н9,р(fv) = 11 (¿4Ak,с эт(fv)-Дф1,h измf j . (1)

Однако формула (1) не содержит какую-либо информацию о погрешностях указанных величин;

3) отсутствует четкий алгоритм определения погрешностей координат источника излучения;

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

В настоящей работе описывается способ определения координат ИРИ, в котором учитываются указанные недостатки.

Определение азимутальных и угломестных пеленгов. Процедура нахождения координат источника излучения базируется на знании азимутальных и угломестных пеленгов, на основании которых записывают уравнения линий на плоскости или в пространстве. Координаты точки пересечения линий и есть координаты источника излучений. Как показано в работе [2], результаты, полученные в одной системе координат, легко пересчитываются в другие системы координат. Выберем декартову прямоугольную систему координат.

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

А0 = У

(2)

получим

соб 0

f

где A =

sin у2 1 - cos у2 sin у3 1 - cos у3 sin у 1 - cos у

Л

0 =

tg 0( Pi-ф 0) ^

Фс

; 7 =

f P2 - PiCos у2^ Рз - Picos уз

n J

v Pn- Picos у

n J

уm — угол между m-м вибратором и направлением отсчета; 0 — азимутальный пеленг; Р — угломестный пеленг; ф0 — начальная фаза сигнала;

R

pm = cos ^0-ут ^cos Р + Ф0

или Рт = cos(0-ут)cosР + Фс; т = 1, ..., n;

PX „ ф0К Л .

где pm =-; ф0 =-; у1 = 0 — начало отсчета; к — длина волны

2кЯ 2kR

сигналов ИРИ; R — радиус антенной системы.

Сразу получаем оценку начальной фазы сигнала ф0, затем определяем оценку азимутального пеленга 0 из найденного значения tg 0(P -ф0) и оценку угломестного пеленга Р . Поскольку получены аналитические формулы для вычисления начальной фазы сигнала ф0, азимутального пеленга 0 , а затем угломестного пеленга Р , то для них достаточно просто вычислить дисперсии, как для функции случайного аргумента [4-6].

В работе [3] приведено и другое решение, в котором аналитически исключается начальная фаза сигнала. Обозначим

(

cos у2 -1-

Р2 - Р

(

sin у

Рз - Р1

Р2 - Р1

I

Lcos у 3 +

Л

Р2 - Р

Р3 - Ри \

= A;

Ръ - Р1

А и В — известные константы, т. е.

sin у3

= Б:

(3)

(4)

J

tg 0 =

A Б

При известном 0 получим

cos Р =

Р2 - Р1

cos(0-у2) -cos 0

(5)

(6)

Пеленги 0 и в определяются по каждому элементу АС. Из полученного ряда значений находим среднее значение пеленгов, их дисперсии и коэффициенты корреляции.

Следует отметить, что операции, имеющие место в формулах (2) и (3), (4)—(6), не представляют большой вычислительной сложности, для их выполнения требуется немного времени, а ошибка в определении пеленгов уменьшается, поскольку в предлагаемом алгоритме учитывают или исключают начальную фазу сигнала ф0, влияющую

на значение пеленгов.

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

Для вычисления дисперсии найденных оценок применяется следующая формула [4-6]:

D[ f ( x) ] = S f^D (*').

В исследуемом случае в качестве f (x) выступают формулы для cos Р, tg 0, ф0. В качестве x, выступают все другие переменные,

P -ф

входящие в формулу, например: cosР = —--. Тогда f (x) = cosР;

cos 0

P = xj; ф0 = x2; cos0 = x3;

D (cos Р) = ]Г( ^ )2 D (x,).

Определение координат источника излучения путем объединения всей информации по пеленгам. Итак, известен набор пеленгов и их погрешностей от разных источников (пеленгаторов разного вида), а также координаты регистраторов сигналов и погрешности этих координат. Азимутальный и угломестный пеленги определяют координаты направляющего вектора прямой в пространстве, проходящей через точку с известными координатами (регистратор сигналов) и через точку с неизвестными координатами x и y (источник излучения). Каноническое уравнение прямой в пространстве, проходящей через точку M1 (x0, y0, z0) параллельно вектору

S = li + mj + nk , имеет вид

x — x0 = y — y0 = z — Z0 l m n

Уравнение этой же прямой можно записать как пересечение двух плоскостей:

"х-хо = ^-Уо , I т

У - Уо = 2 - 2о т п

Уравнение прямой на плоскости (например, плоскости ХУ) имеет

вид

х - х0 _ У - Уо

l m

Для прямой на плоскости введем l = cos a; m = sin а. Тогда х sin а - у cos а = х0 sin а - у0 cos а

или

х tg а - y = b, (7)

где b = хо tg а - yo.

В данном уравнении прямой (7) две случайные величины — tg а и

b. При известных дисперсиях а2(х0); а2(у0); а2(а) получим дисперсию

2 2 2 2 2 2 D(b) = D(xo tg а - yo) = (xo/cos а) а (а) + tg а а (хо) + а (уо);

D(tg а) = (1/cos4а) а2(а).

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

х tg а - y = bi, i = 1, ..., N,

в которой необходимо определить координаты (x, y) источника излучений.

Метод наименьших квадратов применить нельзя, так как он предназначен только для случая, когда в левой части уравнения отсутствуют случайные величины. Если запишем уравнение один раз с правой частью bi, а второй раз — с tg аг-, то получим две прямые, выходящие из точки расположения пеленгатора под разными углами. Истинное направление будет находиться между ними. В данном случае применим методы конфлюэнтного анализа [4] для построения линии ортогональной регрессии, в которой учитываются погрешности всех исходных данных и которая наиболее близка к истинному направлению. Получим следующий функционал [4]:

F = 1/2 (bi - х tg «i + y ))[D (bt) + x2D( tg «i) точка минимума которого определяет точечные оценки координат

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

^ = 0;

дх

=0.

дУ

Ковариационную матрицу полученных точечных оценок х и у (матрицу М) вычисляем при найденных точечных оценках:

Г д2 Г д2Г V

М = дх2 дхду

д2 Г д2 Г

ч дудх д2 у )

Уравнение прямой линии в пространстве

х - х0 = у - у0.

1 т

у - у0 г - ¿0

т

п

эквивалентно системе уравнений плоскостей

т(х - хо) -1 (у -уо) + 02 = 0;

0х + п(у - Уо) - т(г - го) = 0.

Система для определения координат (х, у, г) источника излучения в данном случае будет иметь 2Ы уравнений с тремя неизвестными (х, у, г). Она содержит четыре случайные величины:

ъ\ъ1г х о2(ъ2г х т / 1г, п / т,

где Ъ1 г = (тг/1г)Х0 - У0, Ъ2г = (Пг/тг)у0г - ¿0,.

Очевидно, что эта система распадается на две системы уравнений: на плоскости ХУ и на плоскости YZ.

Функционал конфлюэнтного анализа в данном случае имеет следующий вид [4]:

Г = 1

г=1

(Ъ1г - (т, / I)х + у)2 _ + _ (Ъ2г - (п / тг)у - г)2

о2(ЪЬ.) + х2о2(т / 1г) с2(Ъ21 ) + у2о2(п /тг) Точечные оценки координат ИРИ вычисляем из системы уравнений

^ = о,

дх

^ = о,

ду

^ _ 0.

дг

Ковариационная матрица М точечных оценок х, у и г

М =

д2 ^ д2 ^ д2 ^

дх2 дхду дхдг

д2 ^ д2 ^ д2 ^

дудх ду2 дудг

д2 ^ д2 ^ д2 ^

дгдх дгду дг2

По ковариационным матрицам строим на плоскости эллипс рассеяния, а в пространстве — эллипсоид рассеяния значений координат ИРИ.

Оценки координат ИРИ и их погрешности рассчитывают по явным формулам, на что не требуется много времени.

Реализация алгоритма. На рисунке приведена блок-схема алгоритма определения координат ИРИ. Алгоритм выполняется в следующей последовательности.

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

2. Пеленгующие устройства передают комплексные амплитуды с выходов элементов АС в блок логарифмирования, куда вводят и функцию комплексной огибающей выходов элементов АС для каждой АС.

3. После логарифмирования данные поступают в блок сравнения действительных и мнимых частей, в котором определяется амплитуда зарегистрированного сигнала и формируются уравнения из равенства мнимых частей.

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

5. В вычислителе 1 при полученных точечных оценках пеленгов и начальной фазы сигнала рассчитывается ковариационная матрица оценок пеленгов и начальной фазы сигнала.

6. Данные о пеленгах и их погрешностях поступают в базу пеленгов и их погрешностей, куда также поступают подобные данные с других пеленгующих устройств.

7. Выбираем вид прямой — на плоскости или в пространстве. Согласно этому формируются системы уравнений для прямой на плоскости или в пространстве.

8. Соответственно в вычислителе 2 или 3 определяются точечные оценки координат ИРИ на плоскости или в пространстве.

9. При известных точечных оценках координат ИРИ в вычислителях 4 или 5 определяются ковариационные матрицы оценок на плоскости или в пространстве.

10. По ковариационным матрицам строятся эллипс рассеяния для оценки координат ИРИ на плоскости или эллипсоид рассеяния для оценки координат ИРИ в пространстве.

Блок-схема алгоритма определения координат ИРИ:

1 — пеленгующие устройства; 2 — комплексные амплитуды с выходов элементов АС; 3 — комплексная огибающая выходов элементов АС; 4 — блок логарифмирования; 5 — блок сравнения действительных и мнимых частей; 6 — амплитуда зарегистрированного сигнала; 7 — формирование уравнений из равенства мнимых частей; 8 — решатель систем алгебраических уравнений 9 — точечные оценки пеленгов и начальной фазы сигнала; 10 — вычислитель 1 (вычисление ковариационной матрицы оценок пеленгов и начальной фазы сигнала); 11 — база пеленгов и их погрешностей; 12 — прямая в пространстве; 13 — формирование системы уравнений для прямой на плоскости; 14 — формирование системы уравнений для прямой в пространстве; 15 — вычислитель 2 (точечные оценки координат ИРИ на плоскости); 16 — вычислитель 3 (точечные оценки координат ИРИ в пространстве); 17 — вычислитель 4 (вычисление ковариационной матрицы оценок координат ИРИ на плоскости); 18 — вычислитель 5 (вычисление ковариационной матрицы оценок координат ИРИ в пространстве); 19 — эллипс рассеяния для оценки координат ИРИ; 20 — эллипсоид рассеяния для оценки координат ИРИ; 21 — блок вывода результатов; 22 — выход

Формируется полная информация о координатах ИРИ, позволяющая принимать решения о дальнейших действиях.

Пример. При двух измерениях пеленгов от одного ИРИ получены уравнения двух прямых в пространстве:

(10 - х)/(-5) = (1 - у)/4 = (-3 - 7)/В

и

(1 - х)/4 = (6 - у)/(-1) = (15 - 7)/(-10).

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

40 - 4х = (-5) + 5у; 2 - 2у = (-3) - z

и

(-1) + х = 24 - 4у; 60 - 10у = 15 - z.

Нетрудно проверить, что прямые пересекаются в точке (5, 5, 5).

ЛИТЕРАТУРА

[1] Брызгалов А.П., Блаиров С.С., Хныкин А.В., Фальков Э.Я. Способ определения координат источника радиоизлучения. Пат. РФ № 2339966, G0155/00, 2005.

[2] Иванов Ю.В., Наумов А.С., Саяпин В.Н., Смирнов П.Л., Соломатин А.И., Терентьев А.В., Царик О.В., Шепилов А.М. Способ и устройство определения координат источника радиоизлучения. Пат. РФ № 2419106, G01S13/46, 2006.

[3] Грешилов А.А. Определение амплитуды, азимутов и начальной фазы сигналов. Инженерный журнал: наука и инновации, 2013, вып. № 12. URL: http://engjournal.ru/catalog/appmath/hidden/1156.html.

[4] Грешилов А.А. Математические методы принятия решений. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2006, 584 с.

[5] Грешилов А.А. Некорректные задачи цифровой обработки информации и сигналов. Изд. 2-е, доп. Москва, Логос, 2009, 360 с.

[6] Грешилов А.А., Плохута П.А. Способ определения азимутальных и уг-ломестных пеленгов источников радиоизлучения с повышенным быстродействием. Пат. РФ № 2380720, МПК G 01 S5/04, 27.01.2010.

Статья поступила в редакцию 05.07.2013 г.

Ссылку на эту статью просим оформлять следующим образом: Грешилов А.А. Определение координат источника радиоизлучения. Инженерный журнал: наука и инновации, 2013, вып. 12

URL: http://engjoumal.ru/catalog/appmath/hidden/1156.html

Грешилов Анатолий Антонович родился в 1939 г., окончил. Московский инженерно-физический институт в 1964 г. Д-р техн. наук, профессор кафедры «Высшая математика» МГТУ им. Н.Э. Баумана. Автор более 150 научных работ, в том числе более 30 монографий, 30 авторских свидетельств и патентов в области разработки математических методов учета неопределенности исходной информации в задачах математической физики, распознавания образов, прогнозирования и в других технических приложениях. e-mail: agresh@mail.ru