CC BY
73_ДОКЛАДЫ АН ВШ РФ_
2020_октябрь-декабрь_№ 4 (49)
- ТЕХНИЧЕСКИЕ НАУКИ -
УДК 621.391, 621.396,519.254
ПРОСТРАНСТВЕННО-ВРЕМЕННАЯ ОБРАБОТКА ШИРОКОПОЛОСНЫХ СИГНАЛОВ НА ОСНОВЕ МОДЕЛИ ЛИНЕЙНОГО ПРЕДСКАЗАНИЯ*
Д.Н. Зима, Д.О. Соколова, А.А. Спектор
Новосибирский государственный технический университет
Основные разработки в области радиолокационных систем обзора направлены на исследование их функционирования с учетом влияния различных факторов, как, например, воздействия помех. Наиболее важным представляется обнаружение широкополосного сигнала, позволяющего увеличить разрешающую способность по дальности и скорости. При этом возникают проблемы подавления широкополосных помех существующими способами. В работе развиваются методы обработки широкополосных сигналов в условиях действия активных помех применительно к использованию в системах с различными вариантами пространственно-временных антенных элементов, в частности на примере линейных антенных решеток. Подход основан на представлении сигналов и помех, регистрируемых цифровой антенной решеткой, в виде многомерных пространственно-временных процессов, т. е. функций пространственной и временной координат. Это обусловлено как пространственным распределением элементов антенной решетки, так и пространственным распределением помех. Байесовский обнаружитель сигнала является оптимальным алгоритмом и обладает наилучшими характеристиками, но при этом крайне сложна его практическая реализация, осуществляемая в области пространственно-временных координат. Исследуемые алгоритмы обработки строятся на основе модели линейного предсказания, т. е. путем использования модели марковского случайного процесса для описания помехи на пространственно-распределенных антенных элементах. Особое внимание уделяется разработке алгоритмов, реализуемых при ограниченных вычислительных ресурсах и работе в реальном времени, что представляет собой проблему статистических методов обработки сигналов.
Ключевые слова: пространственно-временной сигнал, широкополосный сигнал, активные помехи, марковский случайный процесс, байесовская обработка, радиолокационные системы, цифровая антенная решетка.
Б01: 10.17212/1727-2769-2020-4-17-31
Введение
В настоящее время интенсивно развиваются методы широкополосной радиолокации. Широкополосные сигналы используют для обзора радиолокационной обстановки, в системах охраны, в системах передачи информации [1]. Это обусловлено необходимостью повышения разрешающей способности радиолокационных систем (РЛС) для измерения координат, разрешения целей по дальности и скорости, распознавания пространственно-распределенных целей. К широкополосным сигналам можно отнести и шумовые сигналы [2]. Наиболее эффективным способом борьбы с активными помехами является применение в РЛС цифровых антенных решеток (ЦАР). ЦАР способны ослабить действие помех путем формирования глубоких провалов в диаграмме направленности в направлении на постановщика помех [3, 4].
*
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта №19-37-90069.
© 2020 Д.Н. Зима, Д.О. Соколова, А .А. Спектор
В настоящее время обычно используется адаптивная обработка сигналов в РЛС с ЦАР. Если расположение источника помех и источника полезного сигнала различно, то применение адаптивных цифровых антенных решеток обеспечивает высокую эффективность подавления помех, которую дополнительно можно повысить путем применения других методов обработки, например, корреляционной обработки [5, 6]. При этом адаптивная обработка не эффективна при борьбе с помехами, приходящими с направлений, близких к направлению прихода полезного сигнала. Основная проблема подавления широкополосных помех заключается в зависимости положения нулей синтезируемой диаграммы как от направления прихода сигнала, так и от частоты, что приводит к невозможности сформировать адаптивной антенной решеткой достаточно глубокие нули диаграммы направленности антенны, соответствующие всем частотным точкам широкополосной помехи [8, 9].
Задача исследования состоит в определении особенностей пространственно-временного представления широкополосного сигнала и влияния на это активных помех. За счет этого можно упростить реализацию алгоритмов обработки сигналов в радиолокационных системах, повысить точность радиотехнических систем.
1. Моделирование пространственно-временных процессов на антенной
решетке
Обнаружение широкополосного сигнала пространственно-распределенными элементами не всегда эффективно, поскольку только в некоторых случаях можно разделить временную и пространственную фильтрацию [9]. Поэтому возникла идея исследования сигнала, регистрируемого на антенной решетке, в качестве многомерного. Рассмотрим процесс моделирования пространственно-временного процесса на примере линейной антенной решетки (ЛАР). В этом случае сигнал будет зависеть от одной пространственной и одной временной переменных.
Допустим, ЛАР имеет N пространственно-распределенных антенных элементов. Временной сигнал и помеха в общем случае могут быть описаны общем выражением вида
u(t) = U(t)cos (ra0t + y(t) -ф),
где модулирующие процессы U (t) и y(t) определяются видом полезного сигнала; ®о - несущая частота; ф - начальная фаза. В общем случае имеем сигналы на выходах антенных элементов:
Uj = Uj (t), i = 1, N.
Пространственный сигнал образуется в любой фиксированный момент времени t совокупностью отсчетов на всех элементах ЛАР. Фиксируя время t, имеем вектор:
II ,|T
U = Ui, U2, ..., Un| .
Считаем фронт волны плоским, тогда сигналы на различных элементах ЛАР отличаются временными сдвигами, зависящими от направления прихода волны 6 . Углу 6 соответствуют задержки между двумя соседними элементами ЛАР Т0 :
d sin (8)
т0 =-• (1)
c
Пространственно-временной сигнал можно представить в виде
Щ (t) = U (t - (i - 1)то) cos (rao (i -1)^0 + V (t - (i - 1)to )),
а при фиксированных 8 и t он становится пространственным сигналом, т. е. колебательной функцией номера антенного элемента i. Обозначим параметр ®Н = юото, который определяет среднюю нормированную частоту пространственных колебаний:
2 d sin (8)
юн = 2п—:-• (2)
Л
Для случая эквидистантной линейной антенной решетки при d = Л /2 выражение (2) принимает вид
гаН = п sin (8).
Без учета флюктуаций, вызванных модуляциями, пространственный сигнал соответствует:
ui = cos(н (i -1)) = cos (п sin(8n )(i -1)) . (3)
Из (3) нетрудно заметить, что при изменении угла 8 изменяется пространственная частота колебаний помехи. При 8 = ±90° частота максимальна. Пространственные сигналы при направлении прихода помех 8 = 90° и 8 = -90° отличаются фазой. При угле 8 = ±90° частота пространственных колебаний максимальна и равна 1Л , пространственный период равен длине волны Л. Очевидно, размещение элементов на расстоянии d = Л /2 соответствует при 8 = ±90° взятию двух отсчетов на интервале одного периода, что, по теореме Котельникова, соответствует минимальной допустимой частоте дискретизации для восстановления непрерывной функции. Таким образом, в каждый момент времени на ЛАР имеется гармонический пространственный сигнал, что справедливо при переходе от дискретной пространственной переменной к непрерывной. Тогда при последовательном попадании на линейку пространственно-распределенных элементов различных участков широкополосного сигнала в дискретном спектре пространственного сигнала видно только одну частоту, независимо от наличия модуляции. Данное свойство удобно при наличии нескольких сигналов помехи с разными направлениями прихода. Данный характер сохраняется при рассмотрении пространственно-временного колебания.
При замене пространственной дискретной переменной i на непрерывную координату х = d(i -1) можно представить (3) в виде
u(t, х) = U | t -dт0 |cos
х 1Í х
t -1 т° Мt -1 т°
(4)
Очевидно, что двумерный сигнал u (t, х) обладает свойством
u(t, х) = u(t - ktX),
(5)
где kt = sin 9 / c - коэффициент преобразования временной и пространственной координат. Тогда, используя (4) и (5), с помощью двумерного преобразования
Фурье спектр пространственно-временного сигнала можно представить в виде • + /2 _ . _ .
U(®t, ®x) = j j ы(t _ ktx)e~j™tte~j™xxdtdx =
-да_ xm/2
+ xm/2 , ,
= j u(tx)e~dtx j e_.(®tkt+®x)*dx =
_ xm /2
sin[(®tkt +®X)Xm / 2]
= U (®t) Xm
(®tkt +®х))Xm /2
где и - частоты, соответствующие временной и пространственной переменным и называемые временная и пространственная частоты. Качественный вид двумерного спектра пространственно-временного сигнала с линейной частотной модуляцией показан на рис. 1. Для различных углов прихода колебания положение активной области спектра различно, что приведено на контурных графиках амплитудного спектра пространственно-временного сигнала на рис. 2.
Рис. 1 - Спектр пространственно-временного сигнала с линейной антенной решетки
Fig. 1 - Spectrum of spatiotemporal signal from linear antenna array
Временная частота, Гц Временная частота, Гц
Рис. 2 - Контурный график спектра пространственно-временного сигнала при угле прихода колебания: а - е = 0°; б - е = 45° Fig. 2 - Contour plot of the spectrum of spatiotemporal signal at the angle of arrival of the oscillation: a - е = 0°; b - е = 45°
2. Понятия узкополосного и широкополосного сигнала
в пространственном смысле
Широкополосными, или сложными, в радиотехнике принято называть сигналы, у которых произведение эффективной ширины спектра на длительность много больше единицы. Большинство используемых сигналов в радиолокации узкополосные, но интерес к широкополосным и сверхширокополосным сигналам нарастает. В РЛС с цифровыми антенными решетками, использующих широкополосные сигналы, недостаточным оказывается деление сигналов на широкополосные и узкополосные, и происходит дальнейшая классификация на широкополосные и узкополосные сигналы в пространственном смысле [9]. Сигнал считается узкополосным в пространственном смысле, если значения его модулирующих функций в любой момент времени на крайних элементах антенной решетки при любом направлении прихода волны практически одинаковы. В иных случаях сигнал рассматривают как широкополосный в пространственном смысле.
Использование широкополосных сигналов в пространственном смысле приводит к изменению процессов, связанных с формированием сигналов и их обработкой. Задача обнаружения широкополосного в пространственно-временном смысле сигнала пространственно-распределенными элементами представляется трудоемкой, поскольку только в некоторых случаях можно разделить временную и пространственную фильтрацию [9]. Только если сигнал относится к числу узкополосных в пространственном смысле, то при его формировании, управлении характеристиками направленности ЦАР, построении алгоритмов обработки возможно разделение операций временных и пространственных преобразований. Для случая широкополосных сигналов в пространственном смысле необходимо разработать эффективный подход к обработке таких сигналов.
3. Байесовская пространственно-временная обработка во временной области
Рассмотрим применение оптимального байесовского критерия обнаружения сигнала для случая обработки пространственно-временного сигнала в областях временной и пространственной переменных. При фиксированном времени / имеем вектор
1|Т
и = И1, И2, ..., и
где иI = иI (/), и = и (/).
Вектор и из-за наличия помех является случайным с гауссовским распределением, а совместная плотность распределения вероятностей при отсутствии и наличии сигнала имеет вид соответственно:
^с(и) =-
1
Wl(U) =
(72л) («Лег |я|) 1
ехр •!-1 и Т Я-1и I/2 Ч 2
((2л) (ег|я|)
I/2 ехр 2(и - Я)Т^-1(и - Я) \.
Обнаружитель работает по методу максимального правдоподобия:
I(и) = ехр {-2(и - Я)ТЯ-1(и - Я) + 2иТЯ-1и| =
w0(U) [ 2 2
= ехр
ЯТ я-1и + и Т Я-1Я - ЯТЯ-1Я
Решающая статистика является корреляционной статистикой и имеет вид
г (и) = итя-1я .
Таким образом, для применения такой обработки необходимо находить обратную корреляционную матрицу Я пространственно-распределенной помехи вида
иПг V) = Е} (1) + иШг V) .
Общий вид корреляционной матрицы:
Я =
Г11 г12
\г11 12
'11
ГП
а ее элементы:
г12 = ищ «)ищ (Г) = Е}=1 Оп/]) (т¥2}) + а25(/1, ¿2),
где £>П} - дисперсия }-й помехи на элементах антенны; т^} - временной сдвиг }-й помехи в ¿1 -м и ¿2 -м элементах ЦАР:
^'>=f j).
8c; , 8^j - векторы групповых направлений на цель и j -ю помеху; (т) -нормированная АКФ j -й помехи:
sin(rc-Af т) я-4/дт '
Adhh (8П) = (i 2 "¿1)dsin (8П). Временная задержка между элементами ¿2 и ¿j, согласно выражению (1):
= Ad12 (8п ) = (¿2 - i1) d sin(8n)
TV2 = = .
12 c c
Применение оптимального байесовского критерия обнаружения сигнала требует обращения корреляционных матриц помехи порядка, соответствующего квадрату произведения числа приемных пространственно-распределенных элементов kN и числа моментов времени кт . Это требует трудоемких расчетов, приводит к появлению ошибок, а при малых углах прихода помехи корреляционная матрица близка к вырожденной матрице.
4. Пространственно-временная обработка на основе модели линейного
предсказания
Проблемы, возникающие при байесовской обработке пространственно-временных сигналов во временной области, можно решить, если использовать модель марковского случайного процесса для описания помехи на пространственно-распределенных антенных элементах.
Обработка строится на модели авторегрессии, т. е. линейного предсказания [10], минимального порядка:
L
u¿ =Z atui-i +v. (6)
l=1
Для квазигармонических процессов это второй порядок. Поэтому при фиксированном t и при единственном источнике помех с углом 8п используется модель
u¿ = ajm^j + a2Ui-2 + v¿,
где v¿ - независимый гауссовский шум предсказания. Для определения этих коэффициентов предсказания используется метод моментов, который основан на использовании соотношений, связывающих искомые параметры с моментами наблюдаемых процессов [11].
В качестве показателя соответствия помехи марковской модели (рис. 3) приведена временная реализация помехи до и после применения обработки. Из рисунка видно, что частота пересечения сигналом нулевого уровня увеличилась, что может служить признаком уменьшения корреляции между соседними отсчетами. На рис. 4 видно, что АКФ результирующего информационного шума стремится к дельтаобразному виду.
Áj (т)
i ¡ A лЛ X AL\ А Л Л / \ / \
/ 4l. A ilyv
г Ii mi v r-\ И у/V v V \J \ / мАл/ V ,W / V V\ f '■■ i Vi V i 5 0 v v
О 10 20 30 40 50 60 70
.......Исходная помеха
— Помеха после декорреляции
Рис. 3 - Временная реализация помехи до и после проведения процедуры декорреляции
Fig. 3 - Time realization of the interference before and after the decorrelation procedure
л =: А л :
а A A 1
\/ И 1 / V У
* У
— Исходная помеха
— Помеха после декорреляции
Рис. 4 - АКФ помехового сигнала до и после проведения декорреляции Fig. 4 - ACF of the interference signal before and after decorrelation
Если сигнал принимается с известного направления на фоне белого шума, то его обработка может быть всегда разделена на пространственную и временную [9, 12]. Тогда в данном случае совместное распределение отсчетов помехи на всех элементах линейной антенной решетки для модели с двумя коэффициентами предсказания на всех I элементах решетки:
N
W0(U) = W0(Mi, U2, ..., Ui ) = W0(Mi, M2)n^(Mi|Mi-1, U-2 ) ,
i=3
где w0(u1, u2) =
1
2ст„ 12
rexp
'a V" '12
плотность вероятности перехода, при этом:
2 2 u - 2v12U]U2 + U2
2°2 (1 - r22)
; n(
Ui Ui-1,
'l—2 ) =
w0(U31) w0(U 21)
Щ(и3г) =■
1
-ехр
-1 и
™0(и 2/0 =
(2*)2У 1
2лл/ Л
ехр
1
~ и 2,ТК2-1и 2/
(7)
(8)
Укороченные векторы сигнала помехи в выражениях (7) и (8) имеют вид:
II ||Т
и2/ = \и,, И,--^ ,
Т
и3/ = ||И/, И/-1,Иг-2|| , а корреляционные матрицы укороченных векторов:
'11 '12 '13
Я2 = '11 '12 , Л = '21 '22 '23
'21 '22
'31 '32 '33
С учетом введенных обозначений, совместное распределение отсчетов помехи на всех N элементах решетки:
^о(и) = ^о(и 22)П^Ти34-
,=3 ^0(и2,0
Рассмотрим гипотезу наличия полезного сигнала:
и'1(и) = ^0(и - £),
где £ - вектор сигнала на решетке в момент времени t. Тогда, отношение правдоподобия:
1(и) = ^1(и) = М'0(и22 -£22) пП-^0(и3/ -£3,0 -0(и2,0
^0(и) ^0(и22) ,=3 ^0(и3,0 ^0(и2, -£2,0
(9)
В (8) сигнальные укороченные вектора имеют вид:
$21 =1 к, 5г -1| ,
£3г =||кг, кг-1, кг-21| .
С учетом введенные обозначений рассмотрим в (9) раздельно дроби:
^0(и3 , - £3,0 =
3 =■
^0(и3,)
= ехр
иТЛ3 1£3г - "2 £3ТЛ3 1£3г
12, = "0(и 2,-£2,0 = ехр
^0(и 2,)
и2гЛ2 1£2г £2г Л2 1£2г
N
I (и) = 122 П-^ = ехр
3 121
и 22 Л21£22-2 £22 Я2 1£22
N г
;ехр < ^ [и31,Л3-1£3г -^Л-^ -
4(
(Я3 1£3г -£21Л2
Т
Решающая статистика 1п(1 (и)) = Ь(и):
Ь(и) = и2Т2Л2 1£22 -1 £22Я2 ^£22 +
N
+Х
г=3
иТЯ3 1£3г - и2гЯ2 1£2г 1Я3 1£3г - £2г'Я2 1£2г )
(10)
В (10) присутствуют слагаемые одинаковые при обеих гипотезах и неинформативные, таким образом они могут быть отнесены к порогу принятия решения, поэтому решающая статистика имеет вид:
N
I (и) = и22Я-1 £22 + X и3.Л3-1£3, - иТ^2,
1=3
(11)
Перейти от пространственной к пространственно-временной обработке сигналов линейной антенной решетки можно путем введения понятия кадра. Кадр -совокупность отсчетов в данный момент времени t. Предполагаем отдельные кадры независимыми.
Т
Имеем вектор и ^) =
^) И_^) )
. В данном случае рассматривается дис-
кретное время t = 1,1, являющееся номером кадра. Решающая статистика строится на основе отношения правдоподобия:
I (и) =
wо(U)
(12)
где
и = ||и(1)Т, и(2)Т,..., и^)Т||Т.
(13)
В предположении о независимости векторов (12), из (13) имеем
I (и) = П I (0(и ()), t=1
где I ^ )(и(t)) определяется формулой, в которой все входящие векторы приобретают верхний индекс t:
Решающая статистика для I кадров имеет корреляционный вид в соответствии с (11):
I
г (и) = £ г(г\и(г)) =
г=1
Таким образом, происходит временное накопление вышеописанной пространственной обработки. Следует отметить, что обработка сигнала при наличии нескольких помех требует большего количества коэффициентов предсказания в (6).
5. Исследование алгоритма пространственно-временной обработки
В работе предлагается подход к обработке, основанный на представлении сигналов и помех, поступающих на элементы ЦАР, как функций временной и пространственных координат. Это обусловлено пространственным распределением приемных антенных элементов антенной решетки и пространственного распределения помех. В качестве полезного сигнала используется широкополосный импульсный сигнал с линейной частотной модуляцией, а в качестве помех рассматриваются активные шумовые помехи, прошедшие входные устройства, которые осуществляют частотную селекцию и формируют спектр помехо-вого колебания.
Оценим выигрыш отношения сигнал/помеха, получаемый при обработке для фиксированного момента времени, т.е. пространственно-временной сигнал трансформирован в пространственный.
Отношение сигнал/помеха на входе пространственного фильтра:
2 = ^ 9вых 2 ,
где £ - амплитуда входного сигнала; ст2 - дисперсия помехи на каждом элементе антенны.
Для байесовской обработки пространственно-временного сигнала во временной отношение сигнал/помеха принимает вид
Р
2 С вых ,„т „-1 „ч
?в2ых = р—=¿а1 к ),
РПвых
где £ = 5,..., |Т = |1,1,..., 1||Т - вектор полезного сигнала, 1 = ||1,1,..., 1||Т -
вектор единичных элементов, Я - корреляционная матрица пространственно-распределенной помехи.
Определим отношение сигнал/помеха на выходе пространственного фильтра для модели с линейным предсказанием (марковская модель второго порядка):
г (и) = о, ^ (и )2 = Рпвых = г (и )2 (и )Т
С учетом описанного в предыдущем разделе и на основе выражения (11) находятся 2(и) и отдельные компоненты 2 (и) .
При теоретическом исследовании критерием качества предложенного метода являлось улучшение отношения сигнал/помеха после пространственной обработки к входному отношению сигнал/помеха. На рис. 5 отображено улучшение отношения сигнал/помеха для байесовской обработки пространственного сигнала и для обработки на основе модели линейного предсказания в зависимости от направления прихода одного помехового колебания и при фиксированном
направлении на полезный сигнал (0С = 0°).
Угол прихода сигнала помехи, град.
— Байесовский алгоритм обработки во временной области — Алгоритм обработки на основе линейного предсказания
Рис. 5. - Зависимость улучшения отношения сигнал/помеха от направления
прихода помехи Fig. 5 - Dependence of the improvement in the signal-to-noise ratio on the direction of arrival of the interference
Из рис. 5 видно, что алгоритмы обработки пространственных сигналов на основе модели линейного предсказания позволяют избежать недостатков, присущих оптимальному байесовскому алгоритму обработки во временной области без
ухудшения выигрыша по отношению сигнал/помеха. При переходе к пространственно-временной обработке на основе марковской модели улучшение отношения сигнал/помеха приближается к значениям, полученным при оптимальной байесовской обработке. Дальнейшее исследование показало, что при увеличении количества коэффициентов предсказания и объема выборки помехи качество обработки увеличивается, но это приводит к значительным вычислительным затратам.
Заключение
В работе предложен и развит подход к синтезу алгоритма обработки сигналов, регистрируемых антенной решеткой, направленные на повышение их быстродействия. На основе теории многомерной обработки сигналов предложен подход к моделированию сигналов и помех, поступающих на элементы антенной решетки, как многомерных пространственно-временных процессов. Для различных углов прихода колебания положение активной области спектра различно, что может стать основой для разделения сигналов и в том числе - для борьбы с помехами. Использование байесовской обработки пространственно-временных сигналов во временной области является оптимальной, но требует трудоемких расчетов, приводит к появлению ошибок, а при малых углах прихода помехи корреляционная матрица близка к вырожденной матрице, что решается путем использования модели марковского случайного процесса для описания помехи на пространственно-распределенных антенных элементах.
ЛИТЕРАТУРА
1. Садомовский А.С. Использование широкополосных сигналов в системах радиосвязи // Радиоэлектронная техника. - Ульяновск, 2011. - № 1 (4). - С. 161-165.
2. Lukin K.A. Radar design using Noise/Randome waveforms // 2006 International Radar Symposium, IRS-2006. - Krakow, 2006. - P. 1-4. - DOI: 10.1109/IRS.2006.4338071.
3. Активные фазированные антенные решетки / под ред. Д.И. Воскресенского,
A.И. Канащенкова. - М.: Радиотехника, 2004. - 488 с.
4. Нечаев Ю.Б., Борисов Д.Н., Пешков И.В. Автокалибровочный алгоритм компенсации амлитудно-фазовых ошибок в каналах цифровой антенной решетки // Вестник ВГУ. Серия: Физика. Математика. - 2011. - № 1. - С. 59-69.
5. Адаптивные антенные решетки. В 2 ч. Ч. 1 / В.А. Григорьев, С.С. Щесняк,
B.Л. Гулюшин и др.; под общ. ред. В.А. Григорьева. - СПб.: Университет ИТМО, 2016. - 179 с.
6. Адаптивные антенные решетки. В 2 ч. Ч. 2. / В.А. Григорьев, С.С. Щесняк, В.Л. Гулюшин и др.; под общ. ред. В.А. Григорьева. - СПб.: Университет ИТМО, 2016. - 118 с.
7. Ратынский М.В. Адаптация и сверхразрешение в антенных решетках. - М.: Радио и связь, 2003. - 200 с.
8. Principles of modern radar. Vol. 3. Radar applications / ed. by W.L. Melvin, J.A. Scheer. -Edison: SciTech Publishing, 2014. - 820 р.
9. Обработка сигналов в радиотехнических системах / под ред. А.П. Лукошкина. - Л.: Изд-во ЛГУ, 1987. - 400 с.
10. Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи. - М.: Советское радио, 1973. - 232 с.
11. Соколова Д.О. Непараметрическое обнаружение и классификация в сейсмических системах охраны: дис. ... канд. техн. наук. - Новосибирск, 2013. - 147 с.
12. Якубов В.П. Статистическая радиофизика. - Томск: НТЛ, 2006. - 132 с.
SPATIOTEMPORAL PROCESSING OF BROADBAND SIGNALS BASED ON THE LINEAR PREDICTION MODEL
Zima D.N., Sokolova D.O., Spector A.A.
Novosibirsk State Technical University, Novosibirsk, Russia
The main developments in the field of radar surveillance systems are aimed at studying their functioning, taking into account the influence of various factors, such as the effect of interference. The most important thing seems to be the detection of a broadband signal, which makes it possible to increase the range and speed resolution. This raises the problem of suppressing broadband interference with existing methods. The paper develops methods for processing broadband signals in the presence of active interference as applied to the use in systems with various variants of spatio-temporal antenna elements, in particular on the example of linear antenna arrays. The approach is based on the representation of signals and interference recorded by a digital antenna array in the form of multidimensional spatiotemporal processes, i.e. functions of spatial and temporal coordinates. This is due to both the spatial distribution of the antenna array elements and the spatial distribution of interference. Bayesian signal detector is the optimal algorithm and has the best characteristics, but at the same time its practical implementation is extremely difficult, carried out in the field of spatiotemporal coordinates. The investigated processing algorithms are based on the linear prediction model, i.e. by using the model of a Markov random process to describe interference on spatially distributed antenna elements. Particular attention is paid to the development of algorithms that can be implemented with limited computing resources and work in real time, which is a problem of statistical methods of signal processing.
Keywords: spatiotemporal signal, broadband signal, active interference, Markov random process, Bayesian processing, radar systems, digital antenna array. DOI: 10.17212/1727-2769-2020-4-17-31
REFERENCES
1. Sadomovskii A.S. Ispol'zovanie shirokopolosnykh signalov v sistemakh radiosvyazi [The use of broadband signals in radio communication systems]. Radioelektronnaya tekhnika [Radioe-lectronic equipment]. Ul'yanovsk, 2011, no. 1 (4), pp. 161-165.
2. Lukin K.A. Radar design using Noise/Randome waveforms. 2006International Radar Symposium, IRS-2006, Krakow, 2006, pp. 1-4. DOI: 10.1109/IRS.2006.4338071.
3. Voskresenskii D.I., Kanashchenkov A.I., eds. Aktivnye fazirovannye antennye reshetki [Active phased antenna arrays]. Moscow, Radiotekhnika Publ., 2004. 488 p.
4. Nechaev Yu.B., Borisov D.N., Peshkov I.V. Avtokalibrovochnyi algoritm kompensatsii am-litudno-fazovykh oshibok v kanalakh tsifrovoi antennoi reshetki [Autocalibration algorithm of compensation of gain and phase imperfections in channels of smart antenna]. Vestnik Vo-ronezhskogo gosudarstvennogo universiteta. Seriya: Fizika. Matematika = Proceedings of Voronezh State University. Series: Physics. Mathematics, 2011, no. 1, pp. 59-69.
5. Grigor'ev V.A., ed. Adaptivnye antennye reshetki. V 2 ch. Ch. 1 [Adaptive antenna arrays. In 2 pt. Pt. 1]. St. Petersburg, ITMO University Publ., 2016. 179 p.
6. Grigor'ev V.A., ed. Adaptivnye antennye reshetki. V 2 ch. Ch. 2 [Adaptive antenna arrays. In 2 pt. Pt. 2]. St. Petersburg, ITMO University Publ., 2016. 118 p.
7. Ratynskii M.V. Adaptatsiya i sverkhrazreshenie v antennykh reshetkakh [Adaptation and superresolution in antenna arrays]. Moscow, Radio i svyaz' Publ., 2003. 200 p.
8. Melvin W.L., Scheer J.A., eds. Principles of modern radar. Vol. 3. Radar applications. Edison, SciTech Publishing, 2014. 820 p.
9. Lukoshkin A.P., ed. Obrabotka signalov v radiotekhnicheskikh sistemakh [Signal processing in radio-technical systems]. Leningrad, Leningrad University Publ., 1987.400 p.
10. Kazakov V.A. Vvedenie v teoriyu markovskikh protsessov i nekotorye radiotekhnicheskie zadachi [Introduction to the theory of Markov processes and some radio engineering problems]. Moscow, Sovetskoe radio Publ., 1973. 232 p.
11. Sokolova D.O. Neparametricheskoe obnaruzhenie i klassifikatsiya v seismicheskikh siste-makh okhrany. Diss. kand. tekhn. nauk [Nonparametric detection and classification in seismic security systems. PhD eng. sci. diss.]. Novosibirsk, 2013. 147 p.
12. Yakubov V.P. Statisticheskaya radiofizika [Statistical radiophysics]. Tomsk, NTL Publ., 2006. 132 p.
СВЕДЕНИЯ ОБ АВТОРАХ
Зима Дарья Николаевна - родилась в 1994 году, аспирант кафедры теоретических основ радиотехники Новосибирского государственного технического университета. Область научных интересов: статистическая обработка сигналов. Опубликовано 12 научных работ. (Адрес: 630073, Россия, г. Новосибирск, пр. Карла Маркса, 20. E-mail: d.zima@corp.nstu.ru).
Zima Daria Nikolaevna (b. 1994) - a postgraduate student at the department of radio engineering fundamentals in the Novosibirsk State Technical University. Her research interests are currently focused on statistical signal processing. She is the author of 12 scientific papers. (Address: 20, Karl Marx Av., Novosibirsk, 630073, Russia. E-mail: d.zima@corp.nstu.ru).
Соколова Дарья Олеговна - родилась в 1984 году, канд. техн. наук, доцент кафедры теоретических основ радиотехники Новосибирского государственного технического университета. Область научных интересов: статистическая обработка сигналов. Опубликовано 35 научных работ. (Адрес: 630073, Россия, г. Новосибирск, пр. Карла Маркса, 20. E-mail: d.sokolova@corp.nstu.ru).
Sokolova Darya Olegovna (b. 1984) - Candidate of Sciences (Eng.), associate professor at the departmentof radio engineering in the Novosibirsk State Technical University. Her research interests are currently focused on statistical signal processing. She is the author of 35 scientific papers. (Address: 20, Karl Marx Av., Novosibirsk, 630073, Russia. E-mail: d.sokolova@corp.nstu.ru).
Спектор Александр Аншелевич - родился в 1946 году, д-р техн. наук, профессор, заведующий" кафедрой" теоретических основ радиотехники Новосибирского государственного технического университета. Область научных интересов: локационные, навигационные и инфо-коммуникационные технологии, цифровые статистические методы обработки сигналов и изображений. Опубликовано около 200 научных работ. (Адрес: 630073, Россия, Новосибирск, пр. Карла Маркса, 20. E-mail: spectoraa@mail.ru).
Spector Alexander Anshelevich (b. 1946) - Doctor of Sciences (Eng.), professor, head of the department of radio engineering fundamentals in the Novosibirsk State Technical University. His research interests are currently focused on location, navigation, and information communication technologies, digital statistical methods of signal and image processing. He is the author of about 200 scientific papers (Address: 20, Karl Marx Av., Novosibirsk, 630073, Russia. E-mail: spectoraa@mail.ru).
Статья поступила 21 ноября 2020 г.
Received November 21, 2020
To Reference:
Zima D.N., Sokolova D.O., Spector A.A. Prostranstvenno-vremennaya obrabotka shirokopo-losnykh signalov na osnove modeli lineinogo predskazaniya [Spatiotemporal processing of broadband signals based on the linear prediction model]. Doklady Akademii nauk vysshei shkoly Rossiiskoi Federatsii = Proceedings of the Russian higher school Academy of sciences, 2020, no. 4 (49), pp. 17-31. DOI: 10.17212/1727-2769-2020-4-17-31.
