УДК 621.396 ГРНТИ 47.49.02
СИНТЕЗ И АНАЛИЗ АЛГОРИТМОВ ТРЁХМЕРНОГО ОПРЕДЕЛЕНИЯ КООРДИНАТ НАЗЕМНЫХ ЦЕЛЕЙ ПО РЕЗУЛЬТАТАМ БОРТОВОГО ПЕЛЕНГОВАНИЯ
В.А. УФАЕВ, доктор технических наук
ВУНЦВВС «ВВА имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж) М.П. БЕЛЯЕВ, кандидат технических наук
ВУНЦ ВВС «ВВА имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж)
С применением критериев максимального правдоподобия и минимума расстояния до линии пеленга с учётом рельефа местности синтезированы алгоритмы трёхмерного определения координат наземных целей по результатам бортового пеленгования. Методами Крамара и линеаризации получены формулы расчёта потенциальной точности, моделированием подтверждена их адекватность.
Ключевые слова: бортовое пеленгование, определение координат, синтез и анализ алгоритмов, рельеф местности, потенциальные точности, имитационное статистическое моделирование.
SYNTHESIS AND ANALYSIS OF GROUND TARGETS THREE-DIMENSIONAL DETERMINATION COORDINATE ALGORITHMS ACCORDING TO THE RESULTS OF THE ONBOARD DIRECTION FINDING
V.A. UFAEV, Doctor of Technical Sciences
MESC AF «N.E. Zhukovsky and Y.A. Gagarin Air Force Academy» (Voronezh)
M.P. BELYAEV, Candidate of Technical Sciences
MESC AF «N.E. Zhukovsky and Y.A. Gagarin Air Force Academy» (Voronezh)
Using the maximum likelihood criteria and the minimum distance to the bearing line, taking into account the terrain, algorithms for the ground targets three-dimensional coordinates determination based on the results of airborne direction finding were synthesized. Using the methods of Cramar and linearization, formulas for calculating potential accuracy have been obtained, and their adequacy has been confirmed by modeling.
Keywords: onboard direction finding, coordinate determination, synthesis and analysis algorithms, terrain, potential accuracy, simulation statistical modeling.
Введение. В работах, посвящённых местоопределению наземных источников радиоизлучения с борта летательного аппарата, как правило, ограничиваются определением координат излучателя на плоскости, горизонтальных координат, и по результатам одномерного по азимуту пеленгования [1, 2].
Потребность трёхмерного в пространстве местоопределения возникает в связи с неопределённостью высоты места излучения, обусловленной перепадами высот местности при расположении излучателя вблизи земной поверхности или высоты поднятия над ней. Здесь более широкие возможности реализуемы с применением способов двухмерного пеленгования. Необходим синтез соответствующих алгоритмов и анализ их потенциальной точности. При этом может быть использован методический аппарат работ [3-5] с учётом специфики новой задачи.
Цель статьи - синтез и анализ алгоритмов трёхмерного определения координат наземных целей по результатам пеленгования с борта летательного аппарата.
Задача синтеза состоит в формировании правила обработки результатов пеленгования, обеспечивающего минимум ошибок местоопределения по заданному критерию.
Определим следующие общие условия её решения. Положение источника излучения зададим в земной (топоцентрической) системе координат, на горизонтальной плоскости: х - абсцисса, у - ордината, и по аппликате-высоте к . Пространственное положение бортового пеленгатора Xп, Уп, Нп, где п = 0,...,N — 1 - номер позиции пеленгования при общем количестве
N, считаем известным.
Далее рассмотрим варианты решения задачи.
Алгоритм А максимального правдоподобия совместной оценки горизонтальных координат и высоты места излучения.
Направление на излучатель зададим посредством вектора двухмерного пеленга, направленного на излучатель из места положения летательного аппарата, и определим его в виде составляющих по осям координат:
In У,h ) =
x - X
dn (x, y, h У
, Rn (x, У, h ) =
У - Y
n
dn (x У, h У
Vn (x, y, h ) =
h - H
_n
dn (x У, h)
(1)
где dn(x,У,h)^^¡x-Xn)r7{y-Yñ)r7(h-Hñ)2
- наклонная дальность до цели.
Эти составляющие соответственно равны: COS в' sin6, COS в' COS 6, sin в, где - П — 6 < П - азимут цели, отсчёт положительных значений от оси ординат У по часовой стрелке, — п/2 — в — П/2 - угол места, отсчёт положительных значений от горизонтальной плоскости вверх.
Измерение вектора пеленга обеспечивается, например фазовым методом с помощью объёмной антенной решётки [6] или способами амплитудного двухмерного пеленгования
[7]. Результаты измерений, отметим их скобкой сверху, In, Rn, Vn являются независимыми нормально распределёнными случайными величинами с дисперсией с2 и математическим ожиданием (1).
Исходя из числа измерений и количества неизвестных параметров число позиций пеленгования должно быть не менее двух N > 2 .
В соответствии с [3, 4] образуем функционал правдоподобия в виде суммы квадратов разностей измеренных и расчётных значений составляющих векторов пеленгов нормированной на удвоенную дисперсию измерений:
f (x,У,h) = С • £(/ - In(x,У,h))2 + (Rn - Rn(x,У,h))2 + (Vn - Vn(x,У,h))2 . (2)
2C n
Этот функционал представляет собой показатель степени экспоненты функции правдоподобия. Оценка координат излучателя есть положение его минимума в области местонахождения источника:
(x,У,h)= argmin f (x,у,h).
(3)
x, y ,h
Потенциальную точность оценки координат по алгоритму (3) найдём методом Кра-мара [8]. Для этого определим информационную матрицу:
а хх а ху а и хИ
А = а ху а уу а и уИ , (4)
а ь хИ а уИ аии _
с элементами в виде вторых производных функционала правдоподобия (2) по координатным осям в точке оценочных координат и измерениях равных истинным значениям. Вследствие последнего условия в процессе дифференцирования следует учитывать равенство нулю разностей в функционале (2), что приводит к комбинациям первых производных, например для варианта ху получим:
а =
хУ
ЛДх,у,И) =±У —Я дЯ I АI V
дхду а п дх ду дх ду дх ду
(5)
д Я (х - X)•(у - Г) д Я ¿2 (х,у,И)-(у - 7 )2
где —Яп =----^-, -Я ~--
дх
—Я=-дИ п
д_ дИ
I =-
—Т=-дИ п
й 1(х, у И) '
(И - Н )•(у - г) д п /
й 1 (х y, И) дх
(И - Нп).(х - 7 ^) д -п) V Т
йп3 (x, ^ И) , V дх
й пy, И)-(И - г)2 и /
¿1 (х y, И) (х,у,и)-(х-х )2
д_ дИ
(x, y, И)
<3 (х y, И) д „ д
ду
дИ
Я
—1=д I
ду ду
Указание зависимости составляющих вектора пеленга Яп, 1п ,Тп от координат (х,у,И) опущено для упрощения записи.
Вследствие однотипности формул для других комбинаций параметров ограничимся приведенным фрагментом (5).
После определения информационной матрицы (4) получим обратную ей матрицу ковариации ошибок:
В = А-
(6)
с элементами Ъ^, где /,] = 0,1, 2 - номера её строк и столбцов.
В последующем сравнительную оценку эффективности алгоритмов будем выполнять по двум показателям: линейной ошибке определения координат в горизонтальной плоскости, как среднего квадратичного дисперсий по осям х и у , и средней квадратиче-ской ошибке определения высоты места излучения:
а =4 Ъо,о + Ъ1
■у] Ъ2,2
(7)
Согласно формуле (3) в алгоритме А необходим поиск минимума трёхмерного функционала. В алгебраической форме это не удаётся, а привлечение численных методов
п
п
требует значительных вычислительных затрат. От этого недостатка свободен следующий алгоритм.
Алгоритм А ' оценки координат по критерию минимума расстояний до линии пеленга. Соответствующий трёхмерный функционал определяется как сумма квадратов расстояний до линий пеленга в координатных плоскостях:
f (x, y, h )=Xfe (x, y, h)+ L\ n (y, h) + L\ n (x, h ))•
(8)
где wn - весовые коэффициенты.
Расстояние до линии пеленга в горизонтальной и вертикальных координатных плоскостях равно:
Lln ^ h )=■
(x - x„ )• rn-(y - yn )• in
№ +12 + V2 '
L2 n (y, , L3 n (x, h)=(x - ^ )V -(h - Hn ) í
(9)
^K+íñ^fi
VRñr+í^+Vñ
Минимизация функционала (8) с составляющими (9) достигается алгебраически, путём его дифференцирования, приравнивания результатов нулю, составления и решения системы линейных уравнений с определением пространственных координат цели по формулам:
x " ci c 2 c3 " -1 "C4 _
y = c3 C 5 C 6 C7
h - c3 -C 6 C8 _C 9 _
(10)
где С3 =-£Vn • I • W, с4 = y(lm(zMR„2 + Vn2)-(Re(z)• Rn + (H)-Vn)• í)• W„,
n n
c6 = -УVn • Rn • wn' c7 = Z(Re(z)^(í„2 + V2)-(]m(z,)• ín + (H)• V„)• Rn)• Wn,
n n
c8 = -V(R 2 +12)• W , c9 = У(lm(z )• V • I + Re(z )• V • R -H •(R2 +12))•
8 ^^ \ n n / n ' 9 ^^ \ V n / n n V n / n n n \ n n //
w
Промежуточные суммы ck вычисляются рекуррентно, добавлением результатов преобразования вновь поступивших измерений к ранее полученным.
В функционал (8) введены весовые коэффициенты Wn. Это обусловлено неравноценностью вклада в оценку координат измерений различных пунктов пеленгования. Действительно, с применением метод а линеаризации функций случайных аргументов [9], в
данном случае аргументы Rn, ín ,Vn, можно показать, что математическое ожидание суммы квадратов расстояний до линии пеленга пропорционально квадрату наклонной дальности до источника от пеленгатора. Исходя из этого, установим обратные квадрату дальности весовые коэффициенты. Наклонная дальность равна (Hп - h)/sinв « Hn / Vn. В приближении принято Hn >> h, что обычно выполняется из условий обеспечения безопасности полёта, а значение синуса угла места заменено его оценочным значением исхо-
дя из определения V (х,у,И) в формуле (1). С учётом приближённого равенства весовые
коэффициенты определяются результатами только текущих измерений: wn = /Н2п .
Алгоритм Б максимального правдоподобия оценки координат с учётом рельефа местно-
сти.
Во многих практических случаях излучатель расположен вблизи земной поверхности и высотой его поднятия над Землёй можно пренебречь.
Допустим, что известна функция рельефа местности как зависимость высоты от горизонтальных координат И(х, у). Отсчёт высоты осуществляют от уровня мирового океана. Необходимые при этом данные содержатся в цифровых картах местности, см. например [10]. При больших высотах полёта и дальностях до источника излучения следует дополнительно учитывать сферичность земной поверхности. Для этого используем уравнение сферы и переопределим функцию рельефа:
И(х, у)=: И(х, у)+ г ■
1 1 +
2 2 х2 + у2
(11)
где Г - средний радиус Земли.
При известной функции рельефа местности И(х, у) преобразования (1)-(3) упрощаются с переходом в двухмерные:
йп(х,у)=4(х - хп)2 + (у - хп)2 + (И(х,у)- Нп)2 , (12)
1п^у) = т-А, Я^у)=ТГ^\, V(x,у)= ^уу)-Нп, (13)
й пу) йп(x, у) йп (х у)
/(х, у)=^С-Т(1п - 1п (х, у ))2 +(Яп - Яп (х, у))2 +(Тп - V (х, у))2, (14)
2 ■ С п
х
, у = а^шт /(х, у).
(15)
После определения горизонтальных координат (15) высота места излучения определяется тривиально по заданной функции рельефа местности:
И = И(х, у).
(16)
В связи с уменьшением количества неизвестных параметров, ограничение по числу позиций пеленгования ослабляется до предельного N > 1.
Применительно к функционалу правдоподобия (14) информационная матрица становится двухмерной:
А =
а а
хх ху
аа
ху уу
(17)
2
г
х, у
Элементы матрицы определяются аналогично (5), необходимые при этом производные находятся по формулам:
—я=-
^ п
дх
д Л
(х -X) + (к(х,у)-Нп)- — к(х,у) -(у -7)
дх )
¿1(x, У )
д
¿П(х,у)- (у-7 ) + (к(х,у)-Нп)■— Н(х,у) -(у -7)
А я =_к_Ф )
ду п ¿1(х у) '
( д Л
д (х,у)- (х-X ) + (к(х,у)-Нп)■— к(х,у) -(х-Хп)
д/ =_к_дх )_
дх " ¿Зп (х, у)
—/ =
^ п
ду
д
(у -7 ) + (к(х,у)-Нп)-— к(х,у) -(х -Хп)
_ду )_
а I(x, у)
—К =
^ п
дх
д ( д Л
(х,у)— к(х,у)- (х-X) + (к(х,у)-Нп)— к(х,у) -(к(х,у)-Нп)
Эх
Эх
аI(x, у)
—Г =
^ п
ду
д ( д Л
¿2(х,у)— к(х,у)- (у-Уп)+(к(х,у)-Нп)— к(х,у) -(к(х,у)-Нп)
ду
ду
(18)
^ п3 (x, у )
Обратная матрице (17) матрица ковариации ошибок В, аналогична (6), но двухмерная. Дисперсия оценок по осям координат, коэффициент корреляции и линейная ошибка определения координат в горизонтальной плоскости определяются через её элементы:
в = Л р =--. Ь°л , с =.ЖТо~. (19)
у 1,1 5 ' ху \т-\т-\ г \ х у ^^
вх = ь°,°,
в - в
х у
Формулу расчёта средней квадратической ошибки определения высоты с учётом функциональной взаимосвязи (16) получим методом линеаризации [9], определяя производные рельефа местности по осям абсцисс и ординат в точке истинного положения цели:
в
1Г к( х у)
кдх
+ в
к(х, у )! + 2^- к(х, у к(х, у )рурх ■ в у . (2°) ду ) дх ду
При одномерном по азимуту пеленговании для расчёта потенциальной точности в (12)-(14) и далее в (17)-(20) компоненты связанные с высотой к(х,у)-Нп, Уп(х,у) и
Уп - Уп(х,у) исключают (обнуляют).
2
Алгоритм Б' оценки координат по критерию минимума расстояний до линии пеленга с учётом рельефа местности.
Казалось бы нет препятствий выполнить синтез по принципу, принятому для алгоритма А', да ещё с упрощением по причине снижения размерности задачи. Но гладко было на бумаге, да забыли про овраги, и в данном случае это не идиома, а по сути, буквально.
Что ж, по аналогии с (8), но памятуя зависимость к(х, у), запишем соответствующий рассматриваемому варианту двухмерный функционал:
I (х, у) _ X (¿2п (х, у) + Ь\п (у) + Ь\п (х)) • ^
(21)
Здесь расстояние до линии пеленга в горизонтальной и вертикальных плоскостях равно:
Ап (х у )=
(х - X )• Я -(у - У )• I
\ п / п \у П / п
л/Я
+12 + V2
ь (у)_ (у - Уп )• Vn-(h(x, у)-Я )• Я (х )= (х— хп )• V-(к(х, у)-Я )•/
2п / I ——9 3п\/
(22)
4 Я 2 +1 2 + V 2
V п п п
Л1Я 2 +1 2 + V 2
V " " "
С учётом квадратичной формы функционала (20) с компонентами (21) минимизация его в алгебраическом виде возможна лишь при линейной аппроксимации земной поверхности к(х,у), что собственно и принято в известных способах [11-13].
Чтобы снять данное ограничение, используем метод последовательных приближений.
На первом этапе, даже не прибегая к линейной аппроксимации, что возможно, определим среднюю высоту местности в рабочей зоне. Для этой средней высоты минимизацией функционала (21) с компонентами (22) получим первичную оценку горизонтальных координат и затем по (16) высоты.
На втором этапе повторим процесс применительно к первичной оценке высоты места излучения.
Преимуществом при этом является применение единого расчётного алгоритма координат источника на плоскости на обоих этапах оценивания.
Определяя положение минимума функционала (20) получим:
х
(с4 + к • с6)• с2 -С •(с5 + к • с7) у_ С •(с5 + к • с7)-С-(с4 + к • с6) с2
С1 • С2
^ у
С1 • С2
34
С 2
. (23)
Величина к есть средняя высота местности или первичная оценка места излучения на первом и втором этапе соответственно.
Исходные величины Ск при этом есть результат функционального преобразования измерений и суммирования по формулам:
c =T(R2 + V2 )■ w', c2 =y(/2 + V 2 )• w', c3 =-YR ■ I- w',
1 ^^ \ n n / n ? 2 ^^ \ n n / n ? 3 ^^ n n n *
n n n
c4 = Z((R 2 + V 2 )■ X -(y ■ R + H ■ v )■ I )■ w,
4 ^^ \\ n n / n \ n n n n / n / n ?
n
=У((/2 + V2)■ Y - (X ■ I + H ■ ?)■ R )■ w',
5 ^^ W n n / n \ n n n n / n / n 5
n
c6 =YV ■ I ■ w', c7 =YV ■ R ■ w,
6 ^^ n n n 5 7 ^^ n n n 5
(24)
где w' = w U R
^ n n \
+12 + V2
При большом количестве засечек пеленга определение величин ск целесообразно проводить рекуррентно, путём накопления преобразованных результатов измерений по (24). В этом случае не требуется регистрации всех предыдущих результатов измерений, достаточно преобразовать вновь полученные результаты и добавить к уже имеющимся. Существенно, что на втором и последующих этапах, далее принято 4 этапа, полный перерасчёт по (24) не выполняется,
достаточно переопределить к и выполнить преобразования согласно (23).
Формулы (23), (24) применимы и для расчёта координат по результатам одномерного пеленгования, при этом число позиций пеленгования не менее двух N > 2 и источник не должен
находиться на линии базы. В этом случае следует установить к = 0, V = 0, Wn = 1. Небольшое увеличение точности даёт повторный расчёт при весовых коэффициентах обратно пропор-
циональных w
квадрату
1/((x - X. )2 +(y - Y, )2).
оценочной горизонтальном дальности
до
источника
Функционалы по критерию минимума расстояний до линии пеленга трёхмерный (8) и двухмерный (21) не являются функционалами правдоподобия, непосредственно к ним метод оценки потенциальной точности Крамара не применим. Погрешности соответствующих оценок координат по алгоритмам ускоренного расчёта А', Б' определим путём имитационного статистического моделирования в сравнении с расчётными потенциальными точностями алгоритмов максимального правдоподобия.
Принят вариант равноточного пеленгования со средней квадратической ошибкой 1 град с помощью объёмной решётки минимального состава из четырёх антенн [5]. Рассматривалась система пеленгаторов, размещённых на одинаковой высоте Hn = 2000 м, эквидистантно на окружности радиусом 5000 м в точках с координатами Xn = 5000 ■ sin(2п ■ N-1 ■ n), Yn = 5000 ■ COs(2n ■ N-1 ■ n). Источник излучения располагался на оси ординат на удалении от цента до 10 км с дискретностью 50 м. Таким образом, линия положения источника проходит через центр системы и опорный пункт пеленгования.
Имитировалась холмистая местность с изменением высоты по закону:
h(x, y )= 0,25 ■ h
Г г
1 + cos
w
2 ■ n ■ U ■ x 20000
лл
JJ
1 + cos
2 ■ n ■ U ■ y 20000
ЛЛЛ
(25)
JJJ
где ктах = 500 м - максимальная высота, и - параметр рельефа (число холмов по координатным осям).
Варианты рельефа с одной и тремя возвышенностями по каждой из осей координат показаны на рисунке 1.
X, км -10
Рисунок 1 - Рельеф местности: а - (7=1, б - (7=3
По совокупности 5 10 статистических экспериментов в каждой точке оценивалась средняя квадратичная погрешность сг определения горизонтальных координат и средняя квадратическая ошибка определения высоты места излучения ск для системы из N=3 пеленгаторов на среднепересечённой местности (=3. Результаты показаны на рис.2-4 точками, сплошной линией даны расчетные зависимости для потенциальной точности, в нижней части рисунков штрих-пунктиром дано сечение рельефа местности по линии положения источника в удобном для восприятия масштабе.
-1 10 -8000 -6000 -4000 -2000
Рисунок 2 - Погрешности определения координат по алгоритму Б' с одномерным пеленгованием
Рисунок 3 - Погрешности определения координат по алгоритмам Б' с двухмерным пеленгованием
Рисунок 4 - Погрешности определения координат по алгоритму А'
Результаты расчётов и моделирования в основном хорошо совпадают. Максимальные различия до 15% наблюдаются рис.4б на дальностях более 6000 м, где высок и общий уровень погрешности. Причинами этого являются квазиоптимальный характер алгоритмов и погрешности расчётов, основанных на параболической максимизации функции правдоподобия и применении метода линеаризации функций случайных аргументов.
Из других результатов отметим следующее.
Наивысшие точности определения координат достигаются при двухмерном пеленговании с учётом рельефа местности по алгоритму Б', до 5,5 раз относительно варианта Б' одномерного пеленгования. Алгоритм оценки А' координат в горизонтальной плоскости занимает промежуточное положение, но существенно уступает им по точности при определении высоты места излучения.
Выводы. Полученные результаты могут быть использованы при ведении и оценке эффективности воздушной радиотехнической разведки, в бортовых системах самонаведения на наземные цели, в том числе при автоматической посадке летательного аппарата по излучению приводного радиомаяка.
и
СПИСОК ЛИТЕРАТУРЫ
1. Мельников Ю.П. Воздушная радиотехническая разведка (методы оценки эффективности). М.: Радиотехника, 2005. 253 с.
2. Вассенков А.В., Гузенко О.Б., Дикарев А.С., Скобелкин В.А. Способ определения координат источника радиоизлучений при амплитудно-фазовой пеленгации с борта летательного аппарата. Патент РФ №2432580, 2011.
3. Рабцун А.В. Оценка максимального правдоподобия координат источников сигналов при многопозиционном пеленговании. М.: Радиотехника, 1987. №9. С. 9-10.
4. Уфаев В.А., Уфаев А.В. Уфаев Д.В. Алгоритмы ускоренного расчёта координат в угломерных системах // Антенны, 2013. Вып. 6 (193). С. 31-38.
5. Уфаев В.А. Определение местоположения наземных целей по результатам двухмерного пеленгования с лётно-подъёмных средств //Антенны, 2015. Вып. 5(216). С. 58-64.
6. Уфаев А.В., Уфаев В.А. Синтез и исследование алгоритмов двухмерного пеленгования с применением объёмной антенной решётки //Антенны, 2013. Вып. 5(192). С. 53-58.
7. Козьмин В.А., Уфаев В.А., Уфаев А.В. Способы амплитудного двухмерного пеленгования //Антенны, 2017. Вып. 6 (238). С. 38-44.
8. Фалькович С.Е. Оценка параметров сигнала. М.: Сов. радио, 1970. 352 с.
9. Вентцель Е.С. Теория вероятностей. М.: Высш. шк., 2002. 420 с.
10. Инструментальные средства визуального программирования ГИС-приложений в среде ОС Windows. Геоинформационная система "Карта 2005". г. Ногинск: КБ Панорама, 2005. www.gisinfo.ru (дата обращения: 21.06.2018).
11. Вассенков А.В., Гузенко О.Б., Дикарев А.С., Скобелкин В.А. Способ определения координат источника радиоизлучений при амплитудно-фазовой пеленгации с борта летательного аппарата. Патент РФ №2432580, 2011.
12. Виноградов А.Д., Дмитриев И.С. Способы и характеристики местоопределения наземного источника радиоизлучения при радиопеленговании с борта летательного аппарата //Антенны, 2016. Вып. 6(226). С. 4-44.
13. Березин А.В., Богданов Ю.Н., Васенков А.В., Виноградов А.Д., Дмитриев И.С., Попов С.А. Способ определения координат наземного источника радиоизлучений при радиопеленговании с борта летательного аппарата. Патент РФ №2610150, 2017.
REFERENCES
1. Mel'nikov YU.P. Vozdushnaya radiotekhnicheskaya razvedka (metody otsenki ehffektivnosti). M.: Radiotekhnika, 2005. 253 p. (in Russian).
2. Vassenkov A.V., Guzenko O.B., Dikarev A.S., Skobelkin V.A. Sposob opredeleniya koordinat istochnika radioizlucheniy pri amplitudno-fazovoy pelengatsii s borta letatel'nogo apparata. Patent RF №2432580. 2011. (in Russian).
3.Rabtsun A.V. Otsenka maksimal'nogo pravdopodobiya koordinat istochnikov signalov pri mnogopozitsionnom pelengovanii. M.: Radiotekhnika, 1987. №9. P. 9-10. (in Russian).
4. Ufaev V.A., Ufaev A.V. Ufaev D.V. Algoritmy uskorennogo raschyota koordinat v uglomernyh sistemah // Antenny, 2013. Vyp. 6 (193). P. 31-38. (in Russian).
5. Ufaev V.A. Opredelenie mestopolozheniya nazemnyh tseley po rezul'tatam dvuhmernogo pelengovaniya s lyotno-pod"yomnyh sredstv //Antenny, 2015. Vyp. 5(216). P. 58-64. (in Russian).
6. Ufaev A.V., Ufaev V.A. Sintez i issledovanie algoritmov dvuhmernogo pelengovaniya s primeneniem ob"yomnoy antennoy reshyotki //Antenny, 2013. Vyp. 5(192). P. 53-58. (in Russian).
7. Koz'min V.A., Ufaev V.A., Ufaev A.V. Sposoby amplitudnogo dvuhmernogo pelengovaniya //Antenny, 2017. Vyp. 6 (238). P. 38-44. (in Russian).
8. Fal'kovich S.E. Otsenka parametrov signala. M.: Sov. radio, 1970. 352 p. (in Russian).
9. Venttsel' E.S. Teoriya veroyatnostey. M.: Vyssh. shk., 2002. 420 p. (in Russian).
10. Instrumental'nye sredstva vizual'nogo programmirovaniya GIS-prilozheniy v srede OS Windows. Geoinformatsionnaya sistema "Karta 2005". g. Noginsk: KB Panorama, 2005. www.gisinfo.ru (data obrashheniya: 21.06.2018). (in Russian).
11.Vassenkov A.V., Guzenko O.B., Dikarev A.S., Skobelkin V.A. Sposob opredeleniya koordinat istochnika radioizlucheniy pri amplitudno-fazovoy pelengatsii s borta letatel'nogo apparata. Patent RF №2432580, 2011. (in Russian).
12. Vinogradov A.D., Dmitriev I.S. Sposoby i harakteristiki mestoopredeleniya nazemnogo istochnika radioizlucheniya pri radiopelengovanii s borta letatel'nogo apparata //Antenny, 2016. Vyp. 6(226). P. 4-44. (in Russian).
13.Berezin A.V., Bogdanov YU.N., Vasenkov A.V., Vinogradov A.D., Dmitriev I.S., Popov S.A. Sposob opredeleniya koordinat nazemnogo istochnika radioizlucheniy pri radiopelengovanii s borta letatel'nogo apparata. Patent RF №2610150, 2017. (in Russian).
© Уфаев В.А., Беляев М.П., 2018
Уфаев Владимир Анатольевич, доктор технических наук, ведущий научный сотрудник 23 отдела научно-исследовательского 2 управления научно-исследовательского научно-исследовательского центра (проблем применения, обеспечения и управления авиацией Военно-воздушных сил), Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А, [email protected].
Беляев Максим Павлович, кандидат технических наук, заместитель начальника 23 отдела научно-исследовательского 2 управления научно-исследовательского научно-исследовательского центра (проблем применения, обеспечения и управления авиацией Военно-воздушных сил), Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж), Россия, 394064, г. Воронеж, ул. Старых Большевиков, 54А, [email protected].