УДК 621.396.67
В. И. Порсев, А. И. Гелесев, А. Г. Красько Угловое сверхразрешение сигналов с использованием «виртуальных» антенных решеток
Проведен анализ опубликованных материалов по «виртуальным» антенным решеткам, выявлены ограничения по их использованию в радиолокационных средствах в условиях априорной неопределенности угловых положений источников сигналов. Показана возможность повышения разрешающей способности по угловым координатам в условиях априорной неопределенности угловых положений источников сигналов с использованием «виртуальной» антенной решетки при типовых отношениях сигнал/шум, применяемых при обработке радиолокационных сигналов. Приведены результаты цифрового моделирования сигналов, подтверждающие аналитические выкладки.
Ключевые слова: угловое сверхразрешение сигналов, «виртуальные» антенные решетки, экстраполяция пространственной структуры.
По мере развития информационных радиоэлектронных систем (РЭС) требования к качеству их основных тактико-технических характеристик постоянно повышаются [1, 2]. Особую важность в процессе функционирования РЭС приобретает достоверность определения состава групповых целей, которая зависит прежде всего от разрешающей способности РЭС по угловым координатам, скорости и дальности [3].
До конца 1950-х гг. в различных областях естественных наук в качестве границы Ц разрешающей способности был принят предел, установленный Релеем в 1888 г. Данный Ц предел обусловливается характеристиками Т реальной аппаратной функции системы ,Р(0), £ определяющими разрешающую способность < по данному параметру 0 [4]. « Величина углового разрешения лимити-
^ рована реальными массогабаритными характе-» ристиками антенной системы. В связи с этим для повышения достоверности определения со-щ става групповых целей проводятся исследова-о ния сверхрелеевского углового разрешения сиг-о налов, превышающего границы релеевского. * В настоящее время основные задачи
¡Е сверхреелевского разрешения относят к новому классу некорректно поставленных задач [4].
Для анализа эффективности разрабаты-^ ваемых методов сверхрелеевского разрешения предложены различные оценки, из которых ис-5 пользуются в основном три показателя [2, 4-7]:
2 _
(П
Г! © Порсев В. И., Гелесев А. И., Красько А. Г., 2019
1) первый определяет относительную величину превышения критерия Релея по анализируемому параметру;
2) второй - ошибку найденных значений параметров;
3) третий - отношение сигнал/шум (ОСШ), при котором достигается заданная степень превышения критерия Релея по анализируемому параметру.
По сравнению с релеевским разрешением предлагаемые методы сверхрелеевской обработки, как новый класс некорректно поставленных задач, менее устойчивы к различным мешающим факторам. В связи с этим первый и третий показатели часто объединяют в один, в результате относительную величину превышения критерия Релея ставят в зависимость от ОСШ [2, 4-7].
Анализ эффективности известных методов углового сверхразрешения, например, алгоритмов Берга, Кейпона, MUSIC и других показывает, что пара равномощных источников, разнесенных на полуширину Д0/2 и треть Д0/3 диаграммы направленности (ДН), разрешаются, если ОСШ для каждого из источников составит 17.. .22 дБ и 25.. .32 дБ соответственно [2, 5, 6].
Помимо обязательного наличия большого ОСШ, существенно превышающего отношение сигнал/шум, достигаемое в радиолокационных системах (РЛС), известным методам углового сверхразрешения присущи и другие недостатки. Так, при использовании данных методов часто отмечается не только
плохая обусловленность корреляционных матриц сигналов, приводящая к неустойчивости вычислительной процедуры, но и потребность априорного знания количества пеленгуемых источников. Последнее нередко приводит к смещенным оценкам угловых координат.
Для преодоления этих недостатков предлагаются различные варианты формирования «виртуальных» антенных решеток (ВАР) [8-19].
Сигналы «виртуальных» элементов ВАР формируют в границах «виртуальной» апертуры, которая находится вне апертуры реальной антенной решетки (РАР), посредством пролонгации сигналов существующих АЭ различными способами [8-19]. Однако аналитически формируемые сигналы должны быть адекватны реальным сигналам, которые могли бы быть получены на входах настоящей апертуры, если бы она была равной «виртуальной».
Из этого следует, что использование сигналов «виртуальных» элементов в ВАР эквивалентно увеличению апертуры реальной антенной системы на величину «виртуальной» апертуры.
При этом следует заметить, что ОСШ на выходах элементов РАР, требуемое для формирования известных ВАР [8-14], намного больше ОСШ, достижимого на выходах антенных элементов в РЛС.
В результате этого известные ВАР применяются при пеленгации, поскольку большое ОСШ на выходе РАР, требуемое для реализации существующих алгоритмов углового сверхрелеевского разрешения, при пеленгации на самом деле существует [8-14].
ОСШ на выходах антенных элементов РАР РЛС намного меньше [1-3], чем требуется для реализации известных ВАР [8-17], поэтому сформировать ВАР на основе антенной системы РЛС, используя известные методы экстраполяции сигналов [8-17], невозможно.
Для того чтобы преодолеть этот недостаток формирования ВАР, в работах [15-17] влияние случайной составляющей на значение принимаемого сигнала предлагается уменьшить за счет коллективной обработки сигналов [15, 17], что должно обеспечить ОСШ на выходах антенных элементов РАР не менее 10.. .14 дБ.
Однако приведенные численные результаты исследований в работах [15-17] показывают, что предлагаемый метод формирования ВАР эффективен лишь в частном случае, когда положение источников сигналов симметрично относительно ожидаемого углового направления приема сигналов 0от.
В случае произвольного расположения источников сигналов в пространстве предлагаемый метод формирования ВАР приводит к большим ошибкам определения их угловых координат.
Для лучшего понимания предлагаемой обработки и упрощения получаемых выражений, но без потери общности решения задачи, вместо фазированных антенных решеток рассмотрим реальные (РЛ) и «виртуальные» (ВЛ) линейные антенные решетки (ЛАР) изотропных АЭ.
Сигнал 0т), принимаемый с ожидаемого углового направления приема 0от и представляющий собой аддитивную смесь сигналов от М-точечных объектов, находящихся в одном разрешаемом объеме на угловых направлениях 0т, может быть записан [2, 3] так:
S(t, Öm ) = Z Sm (t-Tm, 0„ ), m = 1, M. (1)
m=1
Здесь ат(1 - тт, 0т) - принимаемый сигнал от т-го объекта;
тт = Ят /с - время запаздывания прихода фазового фронта сигнала от т-го объекта на первый АЭ;
Ят - модуль радиус-вектора, соединяющий первую точку приема антенной плоского фронта волны от т-го точечного объекта и сам объект.
В большинстве случаев пространственно-временная обработка реализуется на разделяемых (факторизуемых) этапах обработки, которая возможна, если обрабатываемые сигналы узкополосны в пространственно-временном смысле [2, 3, 20-22].
С учетом условия факторизации сигнал от т-го объекта, принимаемого п-м АЭ ЛАР, можно записать в виде [2, 3, 20-22]:
$тп(1 - Тт - Ттп) ~ - Тт)е ^ S0mn}, (2)
та
х
(U
ч
та 0-
та
О
О.
£
ф ц
(Ч
о сч
Ol
<
I
(0 те
0 ü СО те
1
о.
ф
£
и
V CQ
сч
■Clin
с?
сч
■Clin сч
(П (П
где £
£0тп = е -/фп(0т) - комплексные временная и пространственная составляющие узкополосного, в пространственно-временном смысле сигнала;
п = ;
ттп = [(п - 1)й&т0т]/с - время запаздывания прихода плоского фазового фронта сигнала от т-го объекта относительно первого АЭ ЛАР на п-й антенный элемент РЛ с шагом АЭ ё;
£т(? - тт) - комплексная огибающая сигнала;
/0 - несущая частота;
Фп(0т) = [-2п(п - 1)^т0т]/Х;
X - длина волны.
Известно, что комплексная решающая статистика (КРС) 2т принятия решения об обнаружении полезного детерминированного сигнала от одного т-го объекта на фоне аддитивного эквивалентного пространственно-временного белого шума может быть записана [2, 3, 20-22] так:
^т = ^ К01 Бо0,т = (3.1)
2УГ (т ® I, )(1 ® б^) = (3.2)
= CT
= CT-2Yr (le® SI)( ®l). (3.3)
Здесь Y = Setm + N - вектор-столбец комплексной огибающей принимаемого сигнала от m-го объекта;
Setm = Sem ® S0f, S0etm = S0em ® S0f - вектор-столбцы комплексных огибающих принимаемого и ожидаемого полезных сигналов от m-го объекта с пространственными Sem =
= IIS II S = IIS II S = e jyn(emn) S = ||<Jemn|b ^oem _n°oemn||, <Jemn & i ^odmn
= j^, n и временной S* = Ы, Sotk=
= So(tk), k = 1, K структурами, получаемыми при пространственной и временной дискретизации;
0om - угловое направление да-й источник сигнала;
N = IIN„II, n = 1, N - блочный вектор-столбец пространственно-временного белого шума;
Nn = IIWkII, Nnk = Nn(tk), k = 1K - вектор-столбец временной структуры шума n-го пространственного канала.
Корреляционная функция пространственно-временного белого шума при равных
мощностях внутренних шумов АЭ ЛАР К^ в выражении (3.1) определяется следующим равенством [20]:
K et - Ke 1
Kt -ct2Ie 1
It,
(4)
где Kt = M |n] = CT2It - корреляционная матрица временной структуры белого шума; M[*] - математическое ожидание; а2 = M \iNnkN*nk J - дисперсия шума; Ie - пространственная структура корреляционной матрицы Ket, характеризующая независимость шумов различных пространственных каналов;
Ie, It - единичные матрицы размерностью N и K, соответственно;
® - операция кронекеровского произведения матриц;
*, T - верхние индексы комплексного сопряжения и транспонирования матриц.
Соотношение (3.2) определяет алгоритм, в котором сначала осуществлена пространственная обработка, а затем временная. Выражение (3.3) - это алгоритм, в котором первоначально выполнена временная обработка после каждого АЭ, а затем - пространственная [2, 3, 20-22].
При приеме сигналов от М источников сигналов длительностью тс с шириной спектра AFC, не разрешаемых по дальности (xm max -- Tm min) << Тс и частоте (F max - F min) << AFc, используя выражение (3.1), сигнальную составляющую КРС Zm можно представить [2, 3, 20-22] так:
M [Zm ] = ст-2 M
M
KS.
m=1
0m ^ Sot + N) (S,
oet
® Sot )
M
= qO N Kp(0o
m =1
m = 1, M.
г, 0om ),
(5)
Здесь Ыр(0т, 0от) = ^вт^Мт = N Smс(яNd,ЛЛ)
е-/(пёт(#-1)»т/х - функция углового рассогласования пространственных структур сигналов (ФР), принимаемых с истинного 0т и ожидаемого 0от направлений прихода сигнала;
ёт = dcos0m - проекция шага РЛ на плоский фазовый фронт волны, приходящей под ожидаемым углом 0от;
- |0m - 0от| < А©;
q20 = s„t S"ot / g2 - ОСШ m-го объекта после временной согласованной обработки в одном АЭ.
Из выражения (4) следует, что нормированная ФР (НФР) p(0m, 0om) совпадает с известным выражением нормированной ДН ЛАР при равномерном амплитудном и линейном фазовом распределениях.
Для наглядности принципа предлагаемого углового сверхразрешения в условиях априорной неопределенности угловых положений источников сигналов ограничимся моделированием обработки, определяемой выражением (3.3).
В качестве примера была выбрана РЛ с 64 АЭ, ширина ДН которой составляет А© ~ 1,2°. На основе этой РЛ была сформирована «виртуальная» ЛАР с N = 64 реальными и N = 128 «виртуальными» АЭ, сумма которых (N + N) = 192.
С целью сравнения углового разрешения были рассмотрены «виртуальная» и ЛАР с равным числом АЭ - 192, при априорно неопределенном пространственном расположении источников сигналов, не разрешаемых РЛ с 64 АЭ.
Были произвольно выбраны несимметричные угловые положения первого 0i = 11,7°, второго 02 = 12,1° и третьего 03 = 12,6° неразрешаемых источников сигналов, которые смещены относительно ожидаемого углового направления прихода сигнала 0от = 12°.
На рис. 1 приведен результат моделирования с помощью трех НФР диаграммы направленности РЛ с 64 АЭ по мощности без шумов, описывающих прием сигналов от трех неразрешаемых источников в зависимости от угловой расстройки 0 = 0om - 0m. Кроме того, для понимания принципов не только известной обработки (3), но и предлагаемой, были смоделированы распределения в антенных элементах РЛ реальной составляющей пространственной структуры принимаемых сигналов
M
(РРС) Re X (S6mnSo6mn ) и значений ее фаз (РФС)
m=1
M _
arg E(SemnS0emn ). П = 1.N.
m=1
Результаты моделирования РРС и РФС в отсутствие шумов при трех принимаемых не разрешаемых сигналах РЛ с 64 АЭ приведены на рис. 2.
о ^ -------^^
-2,5 -2,0 -1,5 -1,0 -0,5 0 0,5 1,0 1,5 2,0 9
Рис. 1. ДН РЛ с 64 АЭ при обработке трех неразрешаемых источников сигналов
Рис. 2. РРС ( и ) и РФС ( * ) РЛ с 64 АЭ при обработке трех неразрешаемых источников сигналов
На рис. 3 приведены результаты моделирования РРС и РФС без шумов в РЛ со 192 АЭ при фазировании, соответственно, на первый 0от = 01, второй 0от = 02 и третий 0от = 03 источники сигналов.
Известно [14], что в теории прогнозирования любой прогноз будет являться близким к действительности лишь в том случае, если на интервале прогноза будут действовать те же закономерности, которые были на интервале предыстории.
Между тем сравнение РРС реальной ЛАР " с 64 АЭ (см. рис. 2) и РРС реальной ЛАР со £ 192 АЭ (см. рис. 3) показывает, что даже при Б отсутствии шума закономерности РРС, наблю- ч даемые в РЛ с 64 АЭ, не встречаются в РРС ^ после 65-го АЭ, в РЛ со 192 АЭ. I
Задача существенно усложняется при ^ наличии шума. ОСШ при моделировании по- | сле временной обработки в одном простран- го
о сч
01
<
I
(0 та
0 ^
СО та
1
о.
ф
£
и ф
со
сч ■ci-io
с?
сч ■ci-io сч
(П (П
24 48 72 96 120 144 168 п в
24 48 72 96 120 144 168 п
б
Рис. 3. РРС ( д. ) и РФС ( * ) без шумов в РЛ со 192 АЭ при фазировании на разрешаемые источники сигналов на:
а - 01; б - 02; в - 0з
ственном канале РЛ с 64 АЭ было выбрано равным -4,5 дБ. Выбор такого ОСШ после АЭ обусловлен областью стандартных показателей качества радиолокационного обнаружения, которые определяются ОСШ после пространственно-временной обработки. Так, ОСШ после пространственно-временной обработки, опре-
деляемое областью стандартных показателей качества обнаружения, было выбрано равным 13,3 дБ. При этом ОСШ, необходимое для получения стандартных показателей качества радиолокационного обнаружения, на антенных элементах РЛ, составляет величину, намного меньшую, чем требуется для формирования ранее предлагавшихся ВАР [14, 15, 17], а после всей обработки - величину, меньшую, чем для реализации известных алгоритмов углового сверхразрешения, например, Берга, Кейпона, MUSIC и др. [2, 5, 6].
Для достоверности анализа влияния шумов на РРС реальных ЛАР при разрешении разных источников сигналов моделировалась одна и та же реализация пространственно-временного шума.
На рис. 4 приведены РРС и РФС для РЛ со 192 АЭ при приеме сигналов с шумами, когда РЛ сфазирована, соответственно, на первый 0om = 0i, на второй 0от = 02 и на третий 0от = 03 разрешаемые источники сигналов.
Из результатов моделирования следует, что, поскольку шум превалирует, РРС реальных ЛАР со 192 АЭ отличаются друг от друга незначительно (см. рис. 4). Такая нивелировка РРС, несмотря на то что зависимости РРС для разных ожидаемых угловых направлениях приема сигналов без шумов различаются существенно (см. рис. 3), подтверждает низкую вероятность получения точной экстраполяции сигналов известными методами при малом ОСШ на антенных элементах РЛ [17].
Таким образом, при типовых показателях качества обнаружения сигналов радиолокационными средствами в условиях априорной неопределенности угловых положений их источников решение задачи углового сверхразрешения с использованием «виртуальных» антенных решеток представляет не только научный, но и практический интерес.
Приведенные результаты моделирования РРС подтверждают известный вывод о том, что в условиях априорной неопределенности угловых положений источников сигналов, когда ОСШ на антенных элементах меньше (много меньше) единицы, задача экстраполяции пространственной структуры совокупного сигнала (1), сформированного М точечными
0 24 48 72 96 120 144 168 п
0 24 48 72 96 120 144 168 п в
2----
1
О -1
-3
Ш
Ad®
шщ;
I
МЯШИЦНЙ
АЭ, аналогично реальному, может быть записана так:
о = о ]фпв(0т) — 0 -/2я(N+nв-1)rfsmemД,
Поскольку в ряде случаев исследователи стремятся, чтобы апертура ВЛ в G раз превышала апертуру РЛ, то количество АЭ в ВЛ составит
N = СМ,
где G - количество итераций сигналов РЛ, необходимых для формирования требуемого количества Мв.
Тогда ожидаемый фазовый набег сигнала от т-го объекта в пв-м антенном элементе ВЛ фпвй(0т) при ^-й итерации определится соотношением
фпвй(0т) = - 1(0т) + фnвg (0т).
Здесь фnвg(0m) = -2л(пвй - 1)й&т0т]Д - фазовый набег в пвё-м антенном элементе ВЛ при g-й
итерации;_
_nвg = ^ - Щ +1 - ^ - ,(g + Щ - gN = = 1, N - локальный номер антенного элемента ВЛ на g-й итерации;
g = Щ
фм^ - 1(0т) = - 1)Жвт0т/Х - груп-
повой фазовый набег всех антенных элементов ВЛ для g-й итерации.
Поскольку экстраполяция принимаемого сигнала Sefm может осуществляться только в присутствии аддитивного шума, то экстраполируется принимаемая реализация Y:
YG = У^Х = Se0mG ® Т, g = 1Ю
-4
0 24 48 72 96 120 144 168 п
б
Рис. 4. РРС ( л. ) и РФС ( * ) с шумами в РЛ со 192 АЭ при фазировании на разрешаемые источники сигналов:
а - на первый; б - на второй; в - на третий
объектами, находящимися в одном разрешаемом объеме, известными методами не может быть решена.
При приеме сигнала из дальней зоны фазовая структура его пространственной составляющей для любого пв-го «виртуального»
II = lle
-2n(g-1)NdsmQoml\
Здесь 80Ота = || $0
g = 1. о.
Тогда после G итераций принимаемых сигналов РЛ комплексная решающая статистика 2та принятия решения об обнаружении разрешаемого полезного детерминированного сигнала от т-го объекта на фоне аддитивного эквивалентного пространственно-временного белого шума N аналогично выражениям (3.1)-(3.3) может быть записана так:
^тО = ХО КвЮ SoвtmG = (5.1)
= (s;emG ® It )(10S0t^ = (5.2)
та
х
(U
ч
та Q.
та
О О.
Ё V
ц
(Ч
а
о сч
01
<
I
со те
0 ^
CÛ те
1
о.
ф
£
и ф
CÛ
сч ■clin
с?
сч ■clin сч
(П (П
= а-2\Т ( ® Б0,т )(^0бте ®1), (5.3)
где Ке«_=М[Ко ]; g = 1,0;
N0 = |^0от^|| = 80отО® К, ^оетО ^0от£ Бо0т||;
SoеtmG = БоетО ® Зогт - ожидаемый полезный сигнал с учетом О итераций;
1е0 - единичная матрица размерностью
N + #в.
Тогда, принимая во внимание выражение (5.1), сигнальную составляющую КРС 2т0 можно представить в виде формулы
M[ZmG ] = S0tmGK01GSMmG +
M
^ S QtrG K01G SoQtmG, (6)
r=1,r Ф m
где S0tmG = S0mG ® Sotm S0mG = ||S0omG ^©m11
g = lG;
S0trG = S0mrG ® Sotr S0mrG = S0mG ^ C S0r _ S0r ® Sotr;
Ser = 11S 0rn||;
S _j2n(n-l)dsin 0r/X.
S0rn _ e ;
J0tr>
г Ф т, г = 1, М.
Получить аналитическое выражение оценки истинного углового направления на т-й источник сигнала 9т не представляется возможным, поэтому минимизация расстройки \т = (0то - 0т) ожидаемого углового направления приема 0то ВЛ и истинного направления на т-й источник 0т сигнала будет достигнута при максимизации решающей статистики [3, 21]:
тах((М[^то(Ото,-)]), - = И
°то1
Количество точек I = Д0/Д0тов определится требуемой точностью и возможным шагом дискретизации ожидаемого углового направления прихода сигнала Д0тов в ВЛ в пределах ширины Д0 ДН РЛ.
В результате этого максимальное значение ДН ВЛ, сфазированной на ожидаемое угловое направление 0тог-, обеспечит выделение сигналов с углового направления 0т.
Для подтверждения возможности обработки с помощью ВЛ, предложенной авторами
данной статьи в патенте [18], было проведено математическое моделирование, ограниченное двумя итерациями сигналов РЛ. Для наглядности углового сверхразрешения в условиях априорной неопределенности угловых положений источников сигналов вначале осуществлялась временная согласованная обработка сигналов в каждом пространственном канале, а затем - пространственная с помощью ВЛ.
На рис. 5 представлены результаты моделирования без шумов и с пространственно-временными шумами РРС и РФС «виртуальной» ЛАР, формируемой двумя итерациями, с N = 64 реальными и N = 128 «виртуальными» АЭ при максимизации решающей статистики max
®moi
(ReM[ZmG(9mOT)]), i = i = l, I при фазировании, соответственно, на первый 0о1, второй 0о2 и третий 9о3 источники сигналов.
Моделирование идеального случая, когда шумы отсутствуют, приведено с целью пояснения предлагаемого принципа формирования ВЛ.
На рис. 6 приведены результаты моделирования диаграммы направленности по мощности ВЛ со (N + N5) = 192 АЭ в присутствии пространственно-временного шума, разрешающей, соответственно, первый 0о1, второй 0о2 и третий 0о3 источники сигналов, которые не разрешаются РЛ с N = 64 АЭ (см. рис. 1), при максимизации решающей статистики
max(ReM[.Z^(Qmoi)]), i = U.
^moi
Как и ранее, ОСШ при моделировании после временной обработки в одном пространственном канале РЛ с 64 АЭ было равно -4,8 дБ, а на выходе РЛ составило 13,3 дБ.
Подведем итоги моделирования.
1. Наличие превалирующего шума в каналах приема радиолокационных систем исключает возможность формирования ВАР путем экстраполяции сигналов с выходов антенных элементов РАР известными методами.
2. Для углового сверхразрешения радиолокатором при ОСШ, типичном для обнаружения целей, целесообразно реализовать способ, предложенный авторами данной статьи в патенте [18]. При его использовании после требуемого
0 24 48 72 96 120 144 168 п
а
0 24 48 72 96 120 144 168 п
б
2 1 0 -1 -2 -3 -4
аиШКИИЯИЯЙ!^!
О 24 48 72 96 120 144 168 п в
О 24 48 72 96 120 144 168 и г
О 24 48 72 96 120 144 168 п
д
О 24 48 72 96 120 144 168 п
е
Рис. 5. РРС ( ) и РФС ( ± ) без шумов и с шумами «виртуальной» ЛАР с 64 реальными и со 128 «виртуальными» АЭ при разрешении первого источника (а, б), второго - (в, г), третьего - (д, е)
те
х
(U
ч
те Q.
те
о
о.
£
V
ц
(Ч
"@!>??"!>? @> ??@!
о см
Ol
<
I
to та S
0 ü CÛ та
1 О.
<D
â
О <D CQ
CM ■clin 9
CM ■clin
CM
w w
2,0 1,8 1,6 1,4 1,2 1,0 0,8 0,6 0,4 0,2 0
2,0 1,8 1,6 1,4 1,2 1,0 0,8 0,6 0,4 0,2
1
/ 1 л!
/л // ¥ ! /X
Я.// „ ' 1/ 1 ' X. m
,5 -2,0 -1,5 -1,0 -0,5 0 0,5 1,0 1,5 2,0 0 а
— 2,5 -2 — ,0 -1 ,5 -1 ,0 -0 ,5 в 0 (1 5 1 ,0 1 ,5 2 0 0
1
г
1
1
\ К
л, / \1 •
\ Ii у V
W ï/4 А
-2,5 -2,0 -1,5 -1,0 -0,5 0 0,5 1,0 1,5 2,0 6 б
Рис. 6. ДН ВЛ с двойной итерацией РЛ с 64 реальными АЭ при фазировании на первый (а), второй (б) и третий (в) источники сигналов
количества итерации пространственной структуры сигналов РЛ с учетом ожидаемых угловых направлений их приема фазирование ВЛ на соответствующий источник сигналов определится максимизацией решающей статистики.
3. Предлагаемый метод углового сверхразрешения в условиях априорной неопределенности угловых положений источников сигналов
обеспечил разрешение ВЛ сигналов с ОСШ. После временной обработки в одном пространственном канале РЛ с 64 АЭ оно составило -4,8 дБ, что намного меньше результатов, полученных по известным методам экстраполяции сигналов. При этом ОСШ на выходе РЛ равно 13,3 дБ, что меньше, чем требуется для реализации известных алгоритмов после пространственно-временной обработки. Смоделированное значение ОСШ находится в области типовых показателей качества обнаружения радиолокационными средствами.
4. Предлагаемый метод формирования ВАР позволяет разрешать в основном лепестке ДН источники сигналов, расположенные произвольно относительно ожидаемого углового направления приема сигналов (см. рис. 6), которые не разрешаются РАР (см. рис. 2).
5. При реализации предлагаемого метода и алгоритма углового сверхразрешения в условиях априорной неопределенности угловых положений источников сигналов выявлено наличие следующих признаков, позволяющих отнести данную задачу к новому классу некорректно поставленных задач:
• необходимость разработки новых приемов решения задачи углового сверхразрешения не только в теории антенн и пространственно-временного обнаружения сигналов, но и при построении алгоритмов пространственной обработки сигналов антенными решетками;
• формирование ВАР возможно при использовании новых закономерностей, ранее не применявшихся в классической теории не только антенн, но и пространственно-временного обнаружения сигналов;
• потребность в создании новых алгоритмов углового сверхразрешения при использовании нестандартной априорной информации о пространственной структуре источников сигналов ВАР.
6. Реализована двойная итерация РРС реальной ЛАР с 64 АЭ, учитывающая на каждом шаге изменения ожидаемых фазовых структур итерируемых сигналов, необходимых для формировании ВЛ (см. рис. 5, а, в, д). Установлено, что структуры итерируемых сигналов с разными ожидаемыми угловыми направлениями их приема 0о1, 0о2 и 0о3 отличаются друг от друга,
несмотря на одни и те же исходные сигналы РЛ с 64 АЭ (см. рис. 2).
7. РРС и РФС реальной (см. рис. 3, а, б, в) и «виртуальной» (см. рис. 5, а, в, д) ЛАР с одним количеством АЭ существенно отличаются друг от друга при наличии и отсутствии шумов.
8. Наличие нескольких разрешаемых целей приводит к некоторому смещению максимума ДН ВЛ (см. рис. 6). Это вызвано предлагаемым методом решения нового класса некорректно поставленных задач.
Таким образом, предложенный метод обеспечивает в радиолокационных средствах при отношении сигнал/шум, определяемом стандартными показателями качества радиолокационного обнаружения, угловое сверхразрешение «виртуальными» антенными решетками источников сигналов при их произвольном угловом положении. Список литературы
1. Ширман Я. Д. Разрешение и сжатие. М.: Советское радио, 1974. 360 с.
2. Радиоэлектронные системы: Основы построения и теория. Изд. 2-е, перераб. и доп. / Я. Д. Ширман, С. Т. Багдасарян, А. С. Малярен-ко и др.; под ред. Я. Д. Ширмана. М.: Радиотехника, 2007. 510 с.
3. Обработка сигналов в радиотехнических системах / А. Д. Долматов, А. А. Елисеев, А. П. Лу-кошкин и др.; под ред. А. П. Лукошкина. Л.: Изд-во Ленингр. ун-та, 1987. 400 с.
4. Агеев А. Л., Антонов Т. В. О новом классе некорректно поставленных задач // Известия УрГУ. 2008. № 58. С. 24-42.
5. Ратынский М. В. Адаптация и сверхразрешение в антенных решетках. М.: Радио и связь, 2003. 200 с.
6. Варюхин В. А. Основы теории многоканального анализа. Киев: Наук. думка, 2015. 168 с.
7. Климов С. А. Метод повышения разрешающей способности радиолокационных систем при цифровой обработке сигналов // Журнал радиоэлектроники. 2013. № 1. С. 1-25.
8. Buhren M., Pesavento M., Bohme J. F. A new approach to array interpolation by generation of artificial shift invariances: interpolated ESPRIT. Proceedings IEEE Int. Conf. Acoustics, Speech, and Signal Proces. (ICASSP). 2003. Vol. 5. Pp. 362-370.
9. Bajwa W. U., Gedalyahu K., Eldar Y. C.
Identification of Parametric Under spread Linear Systems and Super-Resolution Radar // IEEE Transactions on Signal Processing. 2011. Vol. 52. No. 5. Pp. 2548-2561.
10. Ашихмин С. А, Першин П. В, Федоров С. М. Улучшение разрешения двух источников с помощью формирования «виртуальной» антенной решетки // Воронеж: Вестник Воронежского института МВД России. 2018. № 3. С. 64-69.
11. Herman M. A., Strohmer T. High-resolution radar via compressed sensing // IEEE Trans. Signal Processing. 2009. Vol. 57. № 6. Pp. 2275-2284.
12. Hyberg P. Antenna Array Mapping for DOA Estimation in Radio Signal Reconnaissance. PhD thesis. Royal Institute of Technology (KTH), Stockholm, Sweden. 2005. DOI: 10.11648/j. ajea.20150301.11
13. Belloni F, Richter A, Koivunen V. Reducing Excess Variance in Beamspace Methods for Uniform Circular Array. In Proceedings of the IEEE Workshop on Statistical Signal Processing (SSP), Bordeaux, France. July 17-20, 2005. Pp.940-943.
14. Орлов И. Я, Евсеев А. П., Вьюгин П. Н, Пучков А. В. Анализ проблем создания «сверхразрешающей» цифровой антенной решетки с помощью алгоритма полигармонической экстраполяции // Труды 24-й Междунар. на-учн.-техн. конф. «Радиолокация, навигация, связь». Воронеж, 2018. С. 1-12.
15. Лаговский Б. А. Восстановление изображения групповой цели цифровыми антенными решетками // Антенны. 2011. Вып. 2 (165). С. 40-46.
16. Лаговский Б. А., Шумов И. Ю. Восстановление двумерных изображений источников излучения со сверхразрешением // Антенны. 2013. Вып. 4 (191). С. 60-65.
17. Лаговский Б. А. Сверхразрешение на основе синтеза апертуры цифровыми антенными решетками // Антенны. 2013. Вып. 6 (193). С. 9-16. —
18. Гелесев А. И., Зайцев В. Е, Красько А. Г., I Порсев В. И. Пат. 2642883 РФ. Способ углово- | го сверхразрешения цифровыми антенными ° решетками // Бюл. 2018. № 4. 15 с. ^
19. Порсев В. И, Гелесев А. И, Ворошилин Е. П. те
и др. Выделение парциальных траекторий с |
о
применением синтезируемой виртуальной апертуры антенны РЛС // Вестник воздушно-косми- g ческой обороны. 2016. Вып. 3 (11). С. 87-92. ^
20. Пространственно-временная обработка сигналов / И. Я. Кремер, А. И. Кремер, В. М. Петров и др.; под ред. И. Я. Кремера. М.: Радио и связь, 1987. 224 с.
21. Теория обнаружения сигналов / П. С. Акимов, П. А. Бакут, В. А. Богданович и др.; под
ред. П. А. Бакута. М.: Радио и связь, 1984. 440 с.
22. Гелесев А. И. Радиолокационные и радионавигационные системы. Введение в теорию. М.: Изд-во филиала Воен. Акад. РВСН, 1999. 128 с. Поступила 18.10.19
Порсев Валерий Иосифович - доктор технических наук, профессор, заместитель генерального директора по научной работе Акционерного общества «Всероссийский научно-исследовательский институт радиотехники», г. Москва.
Область научных интересов: радиолокация, математическое и имитационное моделирование.
Гелесев Александр Иванович - доктор технических наук, профессор, начальник сектора Акционерного общества «Всероссийский научно-исследовательский институт радиотехники», г. Москва. Область научных интересов: пространственно-временная обработка сигналов.
Красько Александр Григорьевич - инженер 1-й категории Акционерного общества «Всероссийский научно-исследовательский институт радиотехники», г. Москва.
Область научных интересов: антенные решетки, помехозащищенность, математическое моделирование антенных систем.
о см
OI
<1
I
(1 те
Angular superresolution of signals using virtual antenna arrays
We analysed existing publications concerning virtual antenna arrays and determined the limitations of using them in radar systems for the case of prior uncertainty regarding angular positions of signal sources. The paper shows that it is possible to increase angular coordinate resolution for the case of prior uncertainty regarding angular positions of signal sources by employing a virtual antenna array at typical signal-to-noise ratios used in radar signal processing. We provide results of simulating the signals numerically, which confirm our analytical calculations.
Keywords: angular superresolution of signals, virtual antenna arrays, spatial structure extrapolation.
t" Porsev Valeriy Iosifovich - Doctor of Engineering Sciences, Deputy General Director in Research Work, Joint stock company "All-Russian Scientific and Research Institute of Radio Engineering", Moscow. Science research interests: radiolocation, mathematical and physical simulation.
Gelesev Aleksandr Ivanovich - Doctor of Engineering Sciences, Professor, Head of Sector, Joint stock company "All-Russian Scientific and Research Institute of Radio Engineering", Moscow. Science research interests: spatio-temporal signal processing.
Krasko Aleksandr Grigorevich - Engineer of the 1st rank, Joint stock company "All-Russian Scientific and Research O Institute of Radio Engineering", Moscow.
™ Science research interests: antenna arrays, interference immunity, mathematical simulation of antenna systems.