Preview

Вестник Концерна ВКО «Алмаз – Антей»

Расширенный поиск

Совместное измерение угловых координат близкорасположенных целей линейной антенной решеткой

https://doi.org/10.38013/2542-0542-2022-3-59-66

Полный текст:

Аннотация

В статье рассматривается метод совместной оценки угловых координат близкорасположенных целей линейной антенной решеткой по сигналам после пространственно-временной обработки. Данный метод позволяет повысить точность траекторного сопровождения в условиях высокой плотности воздушного движения и низкой разрешающей способности по дальности. Полученное выражение совместной оценки угловых координат обобщает формулу максимально правдоподобной оценки угловой координаты одиночной цели. При отношении сигнал/шум выше 11 дБ дисперсия оценки угловой координаты одиночной цели приближается к нижней границе Крамера – Рао. Оценки угловых координат двух близкорасположенных целей с одинаковой амплитудой не имеют систематической ошибки, вызванной интерференционным искажением сигнала.

Для цитирования:


Фетисов С.Е. Совместное измерение угловых координат близкорасположенных целей линейной антенной решеткой. Вестник Концерна ВКО «Алмаз – Антей». 2022;(3):59-66. https://doi.org/10.38013/2542-0542-2022-3-59-66

For citation:


Fetisov S.E. Simultaneous measurement of the angular coordinates of closely located targets by a linear array. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2022;(3):59-66. https://doi.org/10.38013/2542-0542-2022-3-59-66

Введение

Для измерения угловой координаты точечного источника излучения в фазированных антенных решетках (ФАР) часто применяется суммарно-разностный метод [1–3], позволяющий достичь высокую точность измерения при низких вычислительных затратах. Однако суммарно-разностный метод при моноимпульсном измерении не позволяет совместно оценить угловые координаты близкорасположенных целей [2].

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

Известные методы углового разрешения и совместной оценки угловых координат подразделяют на четыре группы [5]: методы линейного предсказания, методы группы «Кейпон», проекционные методы и методы параметрического поиска [6][7]. Методы первых трех групп используют стохастическое описание принимаемого сигнала и требуют оценки ковариационной матрицы, что ограничивает их применение со сжатыми сигналами. В сравнении с методами первых трех групп методы параметрического поиска применимы к временным рядам малой длины [5]. Важно отметить, что методы параметрического поиска по сравнению с другими методами дают наилучшее разрешение и высокую точность оценок при корреляции сигналов от нескольких источников излучения [5].

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

Модель пространственных сигналов после диаграммообразования

Фазовое распределение поля плоской волны, падающей под углом θ на M ненаправленных элементов линейной ФАР, выражается сигнальным вектором [1] [10]

где δ = [δ1, …, δM]T – вектор-столбец действительных множителей. Для эквидистантной линейной ФАР с длиной волны λ и расстоянием d между фазовыми центрами смежных элементов δ = [0, …, M – 1]T × 2πd / λ.

Сумму ΣNi=1ai si) сигналов N источников с комплексными амплитудами a = [a1, a2, …, aN]T и угловыми координатами θ = [θ1, …, θN]T представим векторным произведением S ˣ a. Здесь S – сигнальная матрица, составленная из N сигнальных векторов:

Введем лучеобразующую матрицу [1]

где ε1, …, εL – угловые направления L приемных лучей суммарных ДН, w1, …, wM – весовые коэффициенты амплитудного распределения ФАР.

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

y = W Sa + ψ, (1)

где ψ – вектор-столбец случайных нормально распределенных комплексных значений шума с ковариационной матрицей ψ. Для радиолокатора с цифровым диаграммообразованием – оператор математического ожидания, (∙)H – оператор сопряженного транспонирования, n – пространственно-временной белый гауссов шум приемных каналов ФАР с нулевым средним и дисперсией σ2ш [10]. Действительная и мнимая части случайных значений такого шума считаются статистически независимыми и имеют дисперсию σ2ш/2 [11]. Для радиолокационных станций с аналоговым диаграммообразованием коэффициенты матрицы Ψ в общем случае неизвестны и должны измеряться.

Метод оценки угловых координат

Апостериорная плотность вероятности (ПВ) неизвестных параметров N целей при приеме реализации y определяется согласно формуле Байеса:

При отсутствии априорной ПВ p(a,θ) и выборе дельта-функции риска неизвестные параметры находятся в максимуме функции правдоподобия p(y | a,θ) или, что равносильно, в максимуме логарифма отношения правдоподобия (ЛОП) [1] [2]. Запишем ЛОП l = ln(p(y | a,θ)/p(y)) принятого сигнала y с учетом (1), указанных свойств шума и вида совместной ПВ нормально распределенных комплексных случайных величин [11]:

(2)

Раскрыв скобки в выражении (2), получим

или

(3)

где

Найдем экстремум ЛОП (3), приравняв к нулю производную l по a:

Используя формулы дифференцирования скаляра по вектору [12]

где – оператор комплексного сопряжения, после преобразований получим

(4)

Из уравнения (4) следует, что ЛОП достигает максимума при

(5)

Подставив (5) в (3), получим ЛОП, значение которого определяется только количеством целей и их угловыми координатами

(6)

где

Важно отметить, что для N = 1 выражение (6) упрощается и совпадает с известным выражением максимально правдоподобной оценки угловой координаты одиночной цели [13]:

На рисунке 1а показан график функции ЛОП, полученный при моделировании сигналов одиночной цели. Функция меняет свое значение в зависимости от угловой координаты и достигает своего максимума в точке с абсциссой, приближенной к истинному значению угла (1° на рисунке).

График функции ЛОП при моделировании сигналов двух целей приводится в виде линий уровня на рис. 1б. Ввиду равнозначности параметров θ1 и θ2 график симметричен относительно диагональной прямой1 θ1 = θ2 .

Подбор значений параметров θ1, ..., θN для обеспечения максимума (6) при известном количестве целей может происходить последовательно сначала в грубом приближении, потом в виде уточнения. Для грубой оценки угловых координат N целей с разделением углового сектора обзора на K – 1 равных интервалов потребуется выполнить векторно-матричных умножений. Так, например, при угловом секторе обзора 60° и шаге проверки 2° (K = 31) при количестве целей 1, 2, 3 и 4 потребуется выполнить 62, 930, 8990 и 62930 векторно-матричных умножений соответственно. С использованием аппаратного ускорения указанное число матрично-векторных умножений при N < 4 на сегодняшний день не является значительным. Однако при N = 4 проверка гипотез может оказаться затратной как по времени вычисления, так и по занимаемой памяти для хранения предварительно рассчитанных матриц C. Процедура уточнения после грубой оценки угловых координат занимает меньшее время и может проводиться методом половинного деления тестируемого интервала, что требует 3Nn проверок (6), где n – число уточнений, а основание 3 соответствует трем гипотезам: значение угловой координаты цели меньше, больше, либо неизменно на текущем шаге. Использование итеративного уточнения элементов матрицы B–1 методом Ньютона позволит избежать необходимости прямого обращения матрицы B на этапе уточнения оценок. Заметим, что матрица B обратима и функция (6) имеет максимум при N < min (L, M), что ограничивает количество одновременно оцениваемых угловых координат.

Точность оценки угловых координат

Эффективность несмещенных оценок определяется близостью их дисперсии к нижней границе Крамера – Рао (НГКР) [2]. Поскольку выполняется условие дифференцируемости по параметрам a и θ, то оценки асимптотически не смещены и можно ограничиться исследованием их дисперсии [2].

Выразим НГКР для дисперсии оценки угловой координаты одиночной цели по сигналу суммарных ДН. Следует отметить, что известное из [1] выражение потенциальной точности для M-канальной эквидистантной ФАР не учитывает специфику оценки по сигналам суммарных ДН и поэтому в данном исследовании не используется.

Произведем замену u = sin(θ1), A = |a1|, ϕ = arg(a1) и определим коэффициенты матрицы Фишера [2]

Диагональные элементы матрицы Ф–1 определяют нижнюю границу дисперсии оценки неизвестных параметров одиночной цели: синуса угловой координаты , амплитуды и начальной фазы прихода сигнала ϕ̂ . Ограничимся рассмотрением первого элемента на диагонали матрицы Ф–1:

(7)

где U = WH Ψ–1W – матрица M × M; ( ) – оператор произведенияАдамара.

Дисперсия оценки угловой координаты одиночной цели связана с дисперсией оценки синуса угловой координаты приблизительным равенством Далее знак приблизительного равенства опустим, пренебрегая несимметричностью распределения случайной величины оценки на больших углах приема сигнала θ1 относительно нормали ФАР. Для ФАР с цифровым диаграммообразованием выражение (7) с учетом сказанного преобразуется:

(8)

где Q = WH (WWH)–1 W.

Величина  в выражении (8) определяет отношение сигнал/шум (ОСШ) в приемном канале ФАР. Величина ОСШ в децибелах после пространственно-временной обработки

где g – коэффициент увеличения ОСШ в децибелах при наблюдении цели в максимуме пространственно-временной функции неопределенности [1].

На графике рисунка 2 сплошной кривой показана зависимость нижней границы среднеквадратического отклонения (СКО) оценки угловой координаты одиночной цели от значения ОСШ qdB. Отмеченные точки в окрестности этой кривой были получены статистической обработкой результатов многократной оценки с использованием выражения (6).

Приведенные на графике рисунка 2 данные рассчитывались при моделировании сигналов ФАР, содержащей 32 ненаправленных элемента с расстоянием λ/2 между фазовыми центрами смежных элементов и спадающим к краям ФАР амплитудным распределением. Согласно графику, значения СКО оценки угловой координаты одиночной цели приближаются к нижней границе Крамера – Рао при ОСШ более 11 дБ. При меньших значениях наблюдается отклонение от предела, что объясняется ростом числа аномальных измерений [1].

При оценке координат двух близкорасположенных целей с одинаковой амплитудой проверялось отсутствие систематического смещения оценок. Для этого анализировалась статистика оценок при ОСШ qdB 12 дБ и 15 дБ при угловом расстоянии между целями в 1,0 и 1,5 ширины ДН по уровню половинной мощности. Полученные гистограммы оценок приводятся на рисунке 3.

Наименее точный результат наблюдался при угловом расстоянии в одну ширину ДН и ОСШ 12 дБ – в 5% случаев сигналы двух целей не различались, что соответствует на графике рисунка 3 локальному возрастанию плотности оценки в районе 0°. Такой исход наблюдался при синфазном приходе сигнала двух целей и существенном влиянии шума. При повышении ОСШ до 15 дБ такой исход наблюдался в два раза реже, а при угловом расстоянии между целями в 1,5 ширины ДН практически исчезал.

При ОСШ больше 11 дБ и угловом расстоянии между целями более 1 ширины ДН плотность распределения оценки угловой координаты каждой цели имеет симметричный унимодальный вид. Математическое ожидание оценки близко соответствует ее моде, что говорит о незначительной ошибке смещения при взаимно мешающем влиянии двух целей, сигналы которых моделировались со случайной начальной фазой. Отдельные эксперименты также показывают незначительное отклонение гистограмм от симметричной формы при синфазном либо противофазном приходе сигналов.

Заключение

В статье предложен метод совместной оценки угловых координат близкорасположенных целей по сигналам суммарных диаграммах направленности линейной ФАР. Полученное выражение совместной оценки обобщает формулу максимально правдоподобной оценки угловой координаты одиночной цели.

Отличительной особенностью метода является его применимость с временными рядами малой длины, что позволяет использовать метод после пространственно-временной обработки и обнаружения. Результаты статистического моделирования по временному ряду единичной длины подтвердили, что случайная ошибка оценки угловой координаты одиночной цели незначительно отличается от нижней границы Крамера – Рао при отношении сигнал/ шум выше 11 дБ.

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

При исследовании статистических свойств оценок считалось, что количество наблюдаемых целей известно. Если количество целей неизвестно, то оно может определяться последовательной проверкой гипотезы о наличии мешающего сигнала [9] с учетом условия N < min (L, M). На практике количество оцениваемых угловых координат ограничено сложностью вычислительного процесса. Исходя из рассчитанного числа векторно-матричных умножений, применение метода с использованием распространенных средств аппаратного ускорения возможно для совместной оценки угловых координат до четырех точечных целей.

1. Поиск оценок для экономии вычислительного ресурса следует проводить в одной полуплоскости относительно диагонали.

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

1. Ширман Я. Д., Манжос В. Н. Теория и техника обработки радиолокационной информации на фоне помех. – М.: Радио и связь, 1981. – 416 с.

2. Радиотехнические системы: Учеб. для вузов по спец. «Радиотехника» / Ю. П. Гришин, В. П. Ипатов, Ю. М. Казаринов и др.; Под ред. Ю. М. Казаринова. – М.: Высш. шк., 1990. – 496 с.

3. S. M. Sherman and D. K. Barton, Monopulse Principles and Techniques, Artech House Publishers, 2011.

4. Порсев В. И., Николаев А. П., Кривоножко И. С. Многоканальное накопление радиолокационных сигналов, отраженных от высокоскоростных целей, движущихся с ускорением. Вестник Концерна ВКО «Алмаз – Антей». 2018;(1):23–34.

5. Nickel U. Angular superresolution with phased array radar: a review of algorithms and operational constraints. “IEE Proc.”, Pt.F, 1987, Vol. 134, No. 1, pp. 53–59.

6. Ратынский М. В. Адаптация и сверхразрешение в антенных решетках. М.: Радио и связь, 2003, 200 с.

7. Ермолаев В. Т. Методы оценивания параметров источников сигналов и помех, принимаемых антенной решеткой / В. Т. Ермолаев, А. Г. Флаксман; Федер. агентство по образованию, Нижегород. гос. ун-т им. Н. И. Лобачевского. – Н. Новгород : [б. и.], 2007, – 98 с.

8. J. Bohme. Source-parameter estimation by approximate maximum likelihood and nonlinear regression. IEEE Journal of Oceanic Engineering, 1985, vol. 10, no. 3, pp. 206–212.

9. Nickel U. Superresolution using an active antenna array. IEE Conf. Publ., 1982, vol. 216, pp. 87–91.

10. Кремер И. Я. Пространственно-временная обработка сигналов / И. Я. Кремер, А. И. Кремер, В. М. Петров и др.; Под ред. И. Я. Кремера. – М.: Радио и связь, 1984. – 224 с.

11. N. R. Goodman. «Statistical Analysis Based on a Certain Multivariate Complex Gaussian Distribution (An Introduction).» Ann. Math. Statist. 34 (1) 152–177, March, 1963.

12. Kay, S. M., Fundamentals of Statistical Signal Processing: Estimation Theory, PTR Prentice Hall, 1993.

13. Журавлев А. К., Лукошкин А. П., Поддубный С. С. Обработка сигналов в адаптивных антенных решетках. – Л.: Изд-во Ленинградского ун-та, 1983, 240 с.


Об авторе

С. Е. Фетисов
Акционерное общество «Всероссийский научно-исследовательский институт радиотехники»
Россия

Фетисов Сергей Евгеньевич – ведущий инженер-программист отдела 025.
Область научных интересов: радиолокация, цифровая обработка сигналов.

Москва



Рецензия

Для цитирования:


Фетисов С.Е. Совместное измерение угловых координат близкорасположенных целей линейной антенной решеткой. Вестник Концерна ВКО «Алмаз – Антей». 2022;(3):59-66. https://doi.org/10.38013/2542-0542-2022-3-59-66

For citation:


Fetisov S.E. Simultaneous measurement of the angular coordinates of closely located targets by a linear array. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2022;(3):59-66. https://doi.org/10.38013/2542-0542-2022-3-59-66

Просмотров: 86


Creative Commons License
Контент доступен под лицензией Creative Commons Attribution 4.0 License.


ISSN 2542-0542 (Print)