CC BY
84Радиотехнические средства передачи, приема
и обработки сигналов
УДК 621.396.965
М. Е. Шевченко, Д. О. Задирако, Д. Н. Файзуллина, В. Н. Малышев Санкт-Петербургский государственный электротехнический университет
"ЛЭТИ" им. В. И. Ульянова (Ленина) Н. С. Стенюков, М. С. Шмырин АО «НИИ "Вектор"» (Санкт-Петербург)
Методы и алгоритмы панорамного радиомониторинга при малоэлементных антенных решетках1
Представлены методы панорамного радиомониторинга и разработанные на их основе алгоритмы совместного обнаружения и пеленгования в широкой полосе частот. Приведены структурные схемы и панорамы, иллюстрирующие результаты работы алгоритмов.
Антенная решетка, совместное обнаружение и пеленгование, оценка азимута и угла места, КВ-, УКВ-диапазоны, источники радиоизлучения, перекрытие спектров, фазовое пеленгование
Совместные исследования сотрудников кафедры радиоэлектронных средств (РЭС) Санкт-Петербургского государственного электротехнического университета "ЛЭТИ" им. В. И. Ульянова (Ленина) и АО «НИИ "Вектор"» по совершенствованию и дальнейшему развитию методов ведения радиомониторинга в сложных условиях сигнально-помеховой обстановки проводились в течение многих лет. Под руководством проф. В. А. Богдановича с 1970-х гг. проводились теоретические и экспериментальные исследования в области обнаружения, различения и оценивания сигналов в условиях априорной неопределенности сигнально-поме-ховой обстановки. Полученные результаты отражены в десятках статей и в монографиях [1], [2]. Этой тематике посвящались диссертационные исследования, выпускные работы студентов.
Базовая кафедра средств специальной радиоэлектроники, действующая в АО «НИИ "Вектор» с 1978 г., внедряет новые образовательные технологии обучения через исследования, базирующиеся на привлечении студентов в процессе обучения к выполнению конкретных задач предприятия в таких областях, как разработка средств и методов автоматизации процессов обработки
сигналов; разработка систем управления и связи комплексов радиоэлектронной техники; разработка средств радиомониторинга на базе универсальных программируемых приемоанализирующих цифровых модулей; разработка аппаратно-программных средств пеленгования и местоопределения комплексов радиоэлектронной аппаратуры; разработка аппаратно-программных средств для обработки сигналов со сложной частотно-временной структурой и сигналов систем мобильной связи.
Приобретенный опыт позволил совершенствовать существующие методы ведения радиомониторинга и разработать новые подходы и алгоритмы.
Основными задачами панорамного радиомониторинга являются обнаружение сигналов от источников радиоизлучения (ИРИ), оценка несущей частоты обнаруженных сигналов и пеленгование ИРИ, заключающееся в оценивании направлений прихода сигналов от ИРИ. Направление прихода сигналов от ИРИ характеризуется азимутом и углом места, которые иногда объединяют в понятие угловые координаты (УК) ИРИ. Задачей отслеживания интересующего ИРИ кроме пеленгования (оценивания частоты, азимута и угла места) является определение вида модуляции сигнала и его перехват.
1 При подготовке публикации использовались результаты работ по ОКР "Разработка пассивного когерентного локационного комплекса для охраны важных объектов", выполняемой СПбГЭТУ "ЛЭТИ" по договору с АО «НИИ "Вектор"» в рамках комплексного проекта по созданию высокотехнологичного производства при финансовой поддержке работ по проекту Минобрнауки Российской Федерации (постановление Правительства Российской Федерации от 9 апреля 2010 г. № 218).
© Шевченко М. Е., Задирако Д. О., Файзуллина Д. Н., Малышев В. Н., Стенюков Н. С., Шмырин М. С., 2016
5
Панорамный радиомониторинг проводится в широкой полосе обзора (100 кГц...1 МГц и более) и характеризуется высоким уровнем априорной неопределенности сигнально-помеховой обстановки, выражающейся в следующем:
• произвольном расположении соседних сигналов по частоте, возможном перекрытии их спектров из-за внеполосного излучения, несоблюдения условий электромагнитной совместимости, недостаточного ослабления сигнала от удаленного ИРИ в месте проведения радиомониторинга;
• неизвестном уровне шума в различных частотных диапазонах и его частотной окрашенности;
• произвольной заполненности просматриваемого частотного диапазона;
• произвольной длительности сигналов от ИРИ;
• неизвестном виде модуляции сигналов от ИРИ.
Обнаружение сигналов и оценивание УК ИРИ
при панорамном радиомониторинге удобно проводить методами частотной селекции, позволяющими легко разделить сигналы с различными частотами. Эти операции реализуемы в реальном времени благодаря быстродействующим алгоритмам БПФ.
Наряду с этим в научной литературе для панорамного радиомониторинга широко развиваются методы, основанные на выделении сигнального и шумового подпространств. При таком подходе обнаружение и оценивание азимута и угла места сигналов ИРИ производятся во временной области на основе MUSIC- или ESPRIT-подхода [3], [4]. Эти методы вычислительно более затратны, однако обеспечивают большую точность пеленгования и возможность восстановления квадратурных составляющих сигнала ИРИ даже при перекрытии их спектров.
Алгоритмы обнаружения и пеленгования целесообразно разрабатывать на основе подхода совместного обнаружения и оценивания. Термин "совместный" означает, что обнаружение сигналов ИРИ и оценивание направлений их прихода происходит по одним и тем же наблюдаемым выборкам.
Совмещение процессов обнаружения и оценивания УК ИРИ позволяет повысить скорость обзора частотного диапазона по сравнению с методами, в которых обнаружение в широкой полосе частот и пеленгование на фиксированной частоте производятся различными устройствами или на основе анализа различных принятых данных. Кроме того, совмещение позволяет обнаруживать ИРИ с программной перестройкой рабочей частоты (ППРЧ) и определять их УК, определять УК источников кратковременных сигналов, длительность которых не превосходит времени,
требуемого для настройки пеленгатора на заданную частоту при разделении функций обнаружения в широкой полосе и пеленгования.
При панорамном радиомониторинге неопределенность частот ИРИ требует применения точечных алгоритмов обнаружения и пеленгования. Точечные алгоритмы обнаруживают сигнальные составляющие и формируют оценку УК ИРИ для каждого частотного отсчета (ЧО). В интервальных алгоритмах, в отличие от точечных, обнаружение проводится во всем частотном интервале, для которого формируется единственная оценка [5]. В связи с этим интервальные алгоритмы целесообразно применять при слежении за интересующим ИРИ в узкой полосе. Реализовать их в широкой полосе обзора в реальном времени более сложно.
Наблюдаемыми данными для радиомониторинга являются выборки из комплексных отсчетов квадратурных составляющих процессов, принятых М-канальным радиоприемным устройством (РПУ), подключенным к М элементам антенной решетки (АР). Выборки формируются на выходах линейных трактов приемника (ЛТП) в результате квадратурного преобразования (КП) и АЦП. Предполагается идентичность АЧХ и ФЧХ трактов приема. В реальных условиях полной идентичности не достигается; для ее повышения применяются специальные методы и алгоритмы калибровки (см., например, [6]).
Для каждого т-го канала формируется матрица отсчетов с размерами N х Ь:
ХТт =(хт)п1, т = 0, М-1, I = 1, Ь, п = 1, N. (1)
При наличии в наблюдаемых данных сигналов от d > 0 ИРИ и принятых допущениях п-е элементы матриц записываются в виде
d
(■xtm )„i =Yj\bksk к=11
i + (l -1) N
x exp
J 2%fk'
+ (£m )ni
! + (l -1) N
f
д
h
exP (mk )
m = 0, M -1, l = 1, L,
где d - общее число сигналов в наблюдаемых данных в широкой полосе обзора; Ь^ - амплитуда к-го сигнала ИРИ; Sk - нормированная комплексная огибающая к-го сигнала в нулевом канале РПУ; /к - частота к-го сигнала после преобразования центральной частоты просматриваемой полосы на нулевую частоту; /д - частота
дискретизации; у т£ = / (9^, ) - фазовый сдвиг реализации £-го сигнала в т-м канале РПУ относительно нулевого канала (уо£ = 0) (0£, —
азимут и угол места £-го сигнала); ^т — реализации шума неизвестного уровня.
Использование М-элементных линейных (не круговых) АР позволяет обнаружить М -1 ИРИ и потенциально сформировать оценки М - 2 ИРИ на одной частоте. Увеличение числа антенн М также способствует повышению точности формируемых оценок азимута и угла места, однако требует увеличения количества когерентных каналов приемника или реализации коммутируемых трактов приема.
Увеличение количества антенн, особенно при радиомониторинге в КВ-диапазоне приводит к нежелательному для мобильных средств увеличению габаритов и массы радиоаппаратуры, снижению скорости ее развертывания на местности.
При синтезе алгоритмов совместного обнаружения и пеленгования на основе методов с частотной селекцией возможны две стратегии, различающиеся последовательностью проведения обнаружения и пеленгования ИРИ.
Стратегия первичного обнаружения и вторичного пеленгования заключается в первичном энергетическом обнаружении сигналов с последующим пеленгованием в отсчетах, в которых обнаружены сигнальные составляющие. Согласно
этой стратегии обнаружение основано на сравнении статистики, определяемой амплитудными или энергетическими спектрами принятых процессов Б/ (п) (п — номер ЧО) на интервале наблюдения, с порогом заданного уровня вероятности ложной тревоги (рис. 1, а). Для стабилизации ложной тревоги требуется оценка уровня шума. Оценки азимута 9 и угла места ( формируются только в ЧО, в которых обнаружены сигнальные составляющие. Пример частотно-азимутальной панорамы, сформированной на основании энергетического обнаружения (рис. 1, а) представлен на рис. 1, б.
Стратегия первичного пеленгования и вторичного обнаружения состоит в первичном формировании оценок азимута и угла места ИРИ в каждом ЧО, по которым затем выносится решение о наличии или об отсутствии сигнала от ИРИ. В этой стратегии оценки азимута формируются для каждого ЧО. Признаком сигнальных ЧО является группировка оценок азимута, тогда как оценки, сформированные из шумовых отсчетов, распределены по всей области значений от 0 до 360° (рис. 2, а). Вторичное обнаружение основано на фиксации областей группировки оценок (рис. 2, б). Решающей статистикой может служить дисперсия оценки азимута в отсчете. При таком подходе не требуется формирование оценки уровня шума, но для обнаружения необходимо не менее 10 сформированных оценок азимута.
Б/
4 3 2 1
0
200
400
600
800
9, ...
0; —901
I
-180!
9, ...
90
-90
200
400
600
800
б
Рис. 1
9, ...с
90 0
-90 —180
200
400
600
800
б
Рис. 2
п
а
а
0
п
п
Сравнение указанных стратегий обнаружения и пеленгования ИРИ при отсутствии перекрытия спектров соседних сигналов показало следующее:
1. При одних и тех же данных вторичное обнаружение по оценкам азимута не имеет преимущества в пороговом отношении "сигнал/шум" перед первичным энергетическим обнаружением по амплитудному спектру.
2. Вычисление оценок азимута в каждом ЧО и последующее обнаружение требует больших вычислительных ресурсов и объема хранимых в памяти данных по сравнению с обнаружением по амплитудному спектру и дальнейшим пеленгованием в отсчетах, в которых присутствуют сигнальные составляющие.
3. Стратегия первичного энергетического обнаружения по амплитудному спектру и последующее пеленгование в сигнальных ЧО при одном и том же объеме наблюдаемых данных обеспечивает лучшую точность оценок азимута по сравнению с оценками, сформированными по стратегии первичного пеленгования.
Следует ожидать снижения эффективности вторичного обнаружения по оценкам азимута в условиях многолучевого распространения при пеленговании сильных сигналов из-за появления аномальных оценок азимута [7]. Вторичное обнаружение по оценкам азимута недопустимо при наличии в ЧО составляющих сигнала от нескольких ИРИ, так как в ЧО будет формироваться только одна недостоверная (смещенная) оценка азимута, которая не соответствует ни одному ИРИ.
Характеристики обнаружения и точности пеленгования, а также вычислительные преимущества свидетельствуют о целесообразности практической реализации первичного обнаружения по
амплитудному спектру и вторичного пеленгования в ЧО, содержащих сигналы.
Для определения уровня шума в каналах используется устойчивая при 80 %-й заполненности частотного диапазона квантильная оценка уровня шума по цензурированной выборке отсчетов периодограммы [8].
Совместное обнаружение и пеленгование на основе методов с частотной селекцией выполняется по схеме, приведенной на рис. 3. Матрица
наблюдаемых данных ХТт, т = 0, М — 1, (1) преобразуется в частотную область: Хт = (хт ,
п = 1, N, I = 1, Ь, где
N
(Хт ) = ^(х1т о ) ехР 1=1
п = Г7#, I = Г7,
—] 2пп N '
причем w (0 - функция временного окна.
Ь выборок образуют спектрограмму, из которой формируется накопленный амплитудный спектр данных т-го канала, взаимный спектр данных каналов, оценивается уровень шума, обнаруживаются частотные составляющие сигналов; определяются ЧО, в которых присутствуют составляющие одного сигнала; на основе проверки принадлежности составляющих в соседних ЧО одному сигналу формируются оценки частотных интервалов сигналов и оценки УК |§ к, Р к)
фазовым методом.
Минимальное число антенн, необходимое для однозначного формирования оценок азимута и угла места, равно трем (М = 3) (рис. 4, а, Д - расстояние между элементами АР).
V
V
V
РПУ
ЛТП
ЛТП
ЛТП
ЦОС
КП и АЦП
КП и АЦП
КП и АЦП
Вычислитель X,
БПФ
БПФ
БПФ
0
Оценка уровня шума
X,
Оценка уровня шума
X
м—1
Оценка уровня шума
а
л
а
н
и
и
с
X и & 3 т е
2 ч с
я т
л о
в а X
т ы
с н
о т
с о
е т
и н е с а ч
N в
у
р
а
н
ю
О
и
в о § №
лм вм
3 о н я Я 4
8 ° * п
я
В"
я
4
03
н
а
к р
е в
о р
рП
в т
о а
л н
а и
в «
р р
е о
т о
н к
и Й
Й ы
ы в
н о
т л
о и
т у
с а и
ч в
а о
к н ал
е н
ц и
О и с
Рис. 3
/ \ / \ / ч Л/ \А
А ___Y К1
Л 0 3 а
Рис. 4
N
Л 0 S
б
Алгоритм, реализованный в вычислителе (рис. 3), позволяет обнаружить сигналы от ИРИ и сформировать несмещенные оценки азимута и угла места при отсутствии перекрытия спектра соседних сигналов. Однако при наличии всего лишь трех антенн в АР невозможно сформировать оценки перекрывающихся по спектру ИРИ на одной частоте. Поэтому для исключения формирования недостоверных оценок для ЧО, в которых присутствует несколько сигналов, оценки азимута не формируются.
Переход от равносторонней трехэлементной АР (рис. 4, а) к четырехэлементной квадратной АР (рис. 4, б) (M = 4) в КВ-диапазоне позволяет формировать оценки азимута двух ИРИ, частично или полностью перекрывающихся по спектру. Че-тырехэлементная квадратная АР обладает свойством двойной инвариантности к сдвигу, которое использовано при разработке алгоритма [9]. Двойная инвариантность АР (рис. 4, б) к сдвигу проявляется в следующем:
• фазовые сдвиги k-го сигнала в подрешетке, образованной антеннами 0, 1, преобразуются в фазовые сдвиги подрешетки 3, 2 домножением на множитель exp (j2лЛк sin 9 k cos Pk ):
exP (Л 3k ) exP (j'Y 2k )_
exP (j'Yok ) _ exP (j'Yik )_
где A'k = fk ЛIс - расстояние между элементами АР в длинах волн для k-го сигнала; с - скорость света;
• фазовые сдвиги k-го сигнала в подрешетке 0, 3 преобразуются в фазовые сдвиги подрешетки i, 2 домножением на множитель
exp (j 2^ cos 9 k cos Pk ).
Основные этапы алгоритма сохранили преемственность используемого в настоящее время алгоритма трехэлементной АР Реализация алгоритма для четырехэлементной АР потребовала незначительной доработки программного обеспечения.
• exp (j^k sin 9 k cos Pk ),
Для каждого n-го сигнального ЧО формируется выборка
7 =
^n
(x0n )1 (x1n )1 (x2n )i (x3n )1
(x0n )L
- (x1n )l
- (x2n)l
- (x3n )l
составленная из п-х элементов матрицы (1). Вычисляется корреляционная матрица Я = 2п2'н
( н — символ эрмитова сопряжения), которая раскладывается по собственным векторам
Е = [Е0 Е1 Е2 Е3 ],
Е( = [ Е11 Е21 Е3( ]т , ' = 0 3
( т — символ транспонирования) и собственным числам Л = [^0 ^2 Х3], ^0 2 >^3.
Собственные числа Л сравниваются с порогом, определяемым уровнем вероятности ложной тревоги. Фиксируется количество сигналов от ИРИ, которые присутствуют в п-м ЧО. Если ё < 2, формируется сигнальное подпространство
Е _/Е0, ё = 1, Е \[Е0 Е1 ] = 2.
Если ё = 2, то с помощью матричных преобразований с векторами сигнального подпространства формируются оценки Чхп и Чуп операторов поворота:
ехр (]2пЫп х
¥
¥y
X cos 9 n1 cos pn1) 0
exp ( j2%Л'п x x sin 9 n1 cos Pn1)
0
0
exp ( j2%Л'п x
X cos 9 n 2 cos pn2 ) 0
exp(j2лЛП x X sin 9 n2 cos pn2 )_
Операторы поворота связывают между собой элементы векторов сигнального подпространства:
E31 E32 E01 E02
_ E21 E22 _ _ E11 E12 _
" E11 E12 E01 E02
_ E21 E22 _ _ E31 E32 _
¥
Уп
¥ x
С помощью собственных чисел Л y
нок операторов поворота
Л x оце-
лп
2
n
n
n
^ = тх ^ (л Хn )TХ; = ту ^ (Л уп ) ТУ (Тх, Ту - собственные векторы оценок операторов поворота) вычисляются оценки азимутов
0п = (§п1 §п2 ) = - ^ё [®Б (Л уп )/®Б (Лxn )] и углов места
Рп = (Рп1 Рп2 ) = аГССО^аГВ2 (ЛУп ) + Ш^2 (Лxn ).
Для определения указанных оценок необходимо, чтобы собственные числа оценок операторов поворота VХп и ¥уп соответствовали друг
другу. Соответствие достигается вычислением взаимной корреляции собственных векторов Ту и Тх.
При d = 1 для вычисления оценок азимута при АР по рис. 4, б можно использовать выражения
§п =— аг^ [а^ (узп )/ агИ (ть)] аг8 (?1л ) — аг§ (зп )
или
9п = —45 — аг^
аг8 ( 2п )
Результаты статистического имитационного моделирования показали, что для формирования несмещенных оценок азимута двух сигналов необходимо иметь Ь > 30.
На рис. 5-7 приведены фрагменты амплитудного спектра и частотно-азимутальных панорам реальных записей сигналов, полученных в соответствии со схемой на рис. 3 для М = 4. Конфигурация АР соответствовала рис. 4, б, расстояние между антеннами Д = 5.4 м. Записи сделаны в полосе 1 МГц, частота дискретизации 1280 кГц. Перекрытие спектров двух сигналов создавалось перестройкой частоты местного гетеродина.
На рис. 5, а изображен фрагмент амплитудного спектра записи без перекрытия спектров гетеродина 1 и эфирного сигнала 2, расположенного левее сигнала гетеродина на частоте 8172 кГц. Частотно-азимутальные панорамы в отсутствие перекрытия изображены на рис. 5, б. Для сравнения на рис. 5, в приведена частотно-азимутальная панорама, полученная непосредственно из фазовых сдвигов ЧО без проверки наличия в них составляющих двух сигналов.
На обеих панорамах в отсутствие перекрытия оценки азимута ИРИ местного гетеродина и эфирного сигнала одинаковы и составляют 260° и 58°. На панораме (рис. 5, б) на частоте 8040 кГц выявлено 2 ИРИ с азимутами 250° и 125°. На рис. 5, б на частоте 8040 кГц оценки азимута не локализованы.
Результаты, соответствующие 100 %-му перекрытию спектров гетеродина и эфирного сигнала, приведены на рис. 6, а 30 %-му - на рис. 7.
8/ 80 60 40 20 0
t
и_
8000 8050
250 200 150 100 50 0
_1_
_1_
8000 8050
0, ...'
8/ 8/
80 - 80 -
60 - 60 -
40 - 40 -
20 - Л ll.ll . .1 20 Г
0 0
8100 а
8150 /, кГц
8000 8050
8100 а
8150 /, кГц
8000 8050
_1_
250 200 150 100 50 0
_1_
_1_
_1_
250
200
150
100
50 0
_1_
_1_
8100 б
8150 /, кГц
250 • • • • 1 250
200 • ' 1 200 1
150 —• 150 -
100 • 100 —
50 — • » •• • • 50 —
0 1 | | • 0
8000 8050
и
8100
б
8150 /, кГц
8000 8050
_1_
_1_
250200 ^ 150100-*
500
Ч
_1_
_1_
8100 а
8150 /, кГц
_1_
8100
б
8150 /, кГц
»
_1_
8000 8050 8100 8150 /, кГц 8000 8050 8100 8150 /, кГц 8000 8050 8100 8150 /, кГц
в
Рис. 5
в
Рис. 6
в
Рис. 7
2
1
0
0
0
0
0
Частотно-азимутальные панорамы на рис. 6, б и 7, б получены с проверкой наличия составляющих двух сигналов на одной частоте. На рис. 6, б, 7, б оценка азимута ИРИ гетеродина примерно одинакова: 260 ± 2°, оценка эфирного сигнала 57°. Аномальные оценки, обусловленные интерференцией сигнала гетеродина и эфирного сигнала на частоте 8172 кГц, отсутствуют. На частоте 8040 кГц во всех случаях обнаружены 2 эфирных сигнала и сформированы оценки их азимута.
Для сравнения на рис. 6, в, 7, в изображены частотно-азимутальные панорамы, полученные непосредственно из фазовых сдвигов ЧО без проверки наличия в них составляющих двух сигналов. На указанных панорамах при перекрытии спектров сигнала гетеродина с эфирным сигналом выделяются две области группировки оценок азимута в области перекрытия: 210...280° и 41...85°. Однако, во-первых, эти оценки существенно хуже локализованы, чем на рис. 6, б и 7, б соответственно, и, во-вторых, в каждом такте обработки для одного ЧО формировалась только одна оценка азимута, которая принимала одно произвольное значение в зависимости от уровня замираний эфирного сигнала.
Вблизи частоты 8040 кГц на рис. 6, в, 7, в, как и на рис. 5, в, наблюдаются распределенные в широком диапазоне значений оценки азимута. Существенные флуктуации оценок азимута по области значений характерны при некоторых условиях многолучевого распространения сигнала одного ИРИ либо интерференции сигналов нескольких ИРИ.
Поскольку алгоритм с проверкой перекрытия спектров двух сигналов на частоте 8040 кГц сформировал в каждом такте две оценки, которые хорошо локализованы, можно сделать вывод о наличии сигналов двух ИРИ на частоте 8040 кГц.
Совместное обнаружение и пеленгование при круговой коммутируемой М-элементной АР. В УКВ-диапазоне для радиопеленгации применяются круговые коммутируемые М-элементные АР, конструктивно реализованные в одном изделии, что позволяет не тратить время на развертывание на местности.
Коммутация каналов позволяет использовать двух- или трехканальное РПУ вместо М-каналь-ного. В этом случае один канал является опорным, а другие — коммутируемыми. Поскольку когерентная обработка при коммутации невозможна, пеленг оценивается на основе оценок разностей фаз опорного и коммутируемых каналов.
Моделью наблюдаемых данных также является матрица (1), за исключением того, что в каждый момент времени доступны отсчеты опорного и подключенного к одной из т коммутируемых
антенн канала, т = 1, М -1.
Разности фаз между п-ми ЧО опорного и коммутируемого каналов при круговой АР определяются следующим образом:
2%г
У тп ="7— С08 Рк Х„
2%т
сое| 9к -м I-С089£
т = 1, М -1, где г — радиус АР; Xп — длина волны в п-м ЧО.
При выполнении условия г/Х< 0.25 оценки разностей фаз утп формируются однозначно, поэтому можно использовать оптимальную оценку азимута
" М-1
9п = аг^
^ У тп т=1_
81П
2%т ~М
М-1
^ У тп т=1
С08
2%т ~М
(2)
С ростом частоты условие г/X < 0.25 нарушается. При этом некоторые оценки разностей фаз становятся неоднозначными и оценить азимут ИРИ по (2) невозможно.
При г/X > 0.25 используются оценки, максимизирующие модуль произведения вектора оценок взаимного спектра Уп = {Утп} и ожидаемого
вектора оценок Оп (9, () = {ртп (9, р), т = 1, М, для имеющегося отношения г/Xп :
9п, Р п = агБтах "п (9, Р)п .
9, Р
(3)
Составляющие оценок взаимного спектра определяются следующим образом:
(хт Ь (х0 I
1
у =±у
±тп ^
=1
'п1
а элементы ожидаемого вектора
Отп (9, Р) =
= ехР I- ]'
. 2%г
Хп
С08 Р С081 9- 2Лт I- С08 Р С08 9
1 М 1
Схема совместного обнаружения и пеленгования сигналов на основе данных М-элементной коммутируемой АР УКВ-диапазона приведена на рис. 8.
V
л
АР
i -
Г4*
- <
У
ft о й
V £
/
ЛТП
РПУ
КП и АЦП
ЦОС
ЛТП 4 КП БПФ
! и АЦП
БПФ
Вычислитель
Накопление отсчетов взаимного спектра
Оценка уровня шума по взаимному спектру
Накопление и хранение спектров данных
Обнаружение сигналов в ЧО
Оценка Оценка
фазовых количества
сдвигов в ЧО сигналов в ЧО
Оценка азимута и угла места при наличии не более
Рис. 8
Временные наблюдаемые данные (1) от коммутируемого и опорного каналов преобразуются с помощью БПФ в частотную область:
х0/ =[(х0 )1/ ••• (х0 и
Хт1 = [(хт )ц ••• (хт )N1 ] ,
т = 1, М — 1, I = 17. Далее определяется взаимный спектр
(Хт )/ = (хт )1/ (х0 )1/ • (хт )N/ (х0 )N1.
оценивается уровень шума взаимного спектра ат1 и взаимный спектр нормируется на этот уровень: (Хт )// стт/.
Обнаружение в каждом ЧО выполняется по накопленному суммарному взаимному спектру
X =
1
_fV (m)
Отсчеты накопленного суммарного взаимного спектра X = (Х1...XN) сравниваются с порогом,
установленным в соответствии с заданной вероятностью ложной тревоги, и определяются сигнальные ЧО. В сигнальных ЧО определяется число ИРИ, сигнальные составляющие которых присутствуют в отсчете.
Для круговой АР максимальное число ИРИ, для которых на одной частоте можно сформировать практически несмещенные оценки азимута, равно (М/ 2 — 1) [10].
Сигналы от различных ИРИ приходят с разных направлений, и уровни их спектральных составляющих не коррелированы между собой во времени. Поэтому для оценки числа ИРИ следует сформировать матрицу 12
7 =
^n
(x1 )n1 (Xo)
ni
ni
(x0 )ni (x°) (xM-1 )ni (x0 )
ni
(x0 )ni (xo)
ni
(x1 )nL (x* L (x0 )nL (x0 L
(xM-1 )nL (x0 L
(x0 )nL (x0 L
состоящую из взаимных произведений отсчетов опорного и коммутируемых каналов, нормированных на мощность отсчета опорного канала:
(xm )nl (xo)
nl
(x0 )nl (xo)
m = 1, M, l = 1, L.
nl
Количество значимых сингулярных чисел Zn
или собственных чисел Zn Z„ соответствует числу сигналов, присутствующих в ЧО.
Если это количество равно единице, то для формирования оценок азимута применимы выражения (2) или (3). При наличии более одного сигнала оценки азимута формируются на основе алгоритмов ESPRIT или MUSIC [3], [4], [11] из сингулярных векторов Zn или собственных векторов Z.
Коммутирование приемных трактов ограничивает возможности обработки и снижает помехоустойчивость алгоритмов обнаружения и пеленгования. Использование многоканальных когерентных цифровых приемных устройств УКВ-диапазона, разработанных в СПбГЭТУ "ЛЭТИ" и проходящих натурные испытания, позволит избежать коммутации сигналов элементов АР и реализовать когерентные алгоритмы обработки.
Совместное обнаружение и пеленгование на основе методов выделения сигнального и шумового подпространств. Обнаружение и оце-
нивание производятся во временной области на основе MUSIC- или ESPRIT-подхода. Из многоканальных принятых данных формируются кор-реляционные матрицы наблюдений. Количество собственных чисел или сингулярных чисел такой матрицы, превысивших порог, определенный исходя из заданного уровня вероятности ложной тревоги, указывают число ИРИ, а соответствующие им собственные (сингулярные) векторы образуют сигнальное подпространство. Оставшиеся векторы образуют шумовое подпространство.
При ESPRIT-подходе оценки параметров ИРИ формируются на основе векторов сигнального подпространства, при MUSIC - на основе векторов шумового подпространства.
На основе ESPRIT-подхода разработаны алгоритмы совместного обнаружения и пеленгования при использовании трех- и семиэлементной АР во временной области, подробно описанные в [11].
Алгоритмы, синтезированные на основе выделения сигнального и шумового подпространств, по принципу действия являются алгоритмами совместного обнаружения и пеленгования.
Дискретизация гармонического колебания обусловливает появление фазовых сдвигов между дискретными отсчетами в выборке. Фазовый сдвиг между n-м и (n + i)-м отсчетами, принадлежащими принятому от k-го ИРИ сигналу, в каждой из временных выборок (1) составляет 9k = 2%fk^/д. Формирование оценок частот ИРИ, неперекрывающихся по спектру, основано на определении $k , k = 1, d. Тогда оценки частоты fk =%/д/(2га'X k = Ъ d
Вычисление оценки $k основано на формировании хотя бы одной матрицы, строки которой образованы отсчетами наблюдаемых временных данных канала, сдвинутыми на каждой последующей строке на одну позицию относительно предыдущей [4]. Число строк определяется параметром сдвига V < NL/2, V > d.
Для совместного обнаружения и оценивания частоты и УК ИРИ требуется сформировать такие матрицы для каждого m-го канала:
7T =
(xtm )i (xtm )2 (xtm )2 (xtm )3
(xtm )NL-V+1 (xtm )NL-V
(xtm )v (xtm 1
V+1
(xtm )
NL
m = 0, M -1.
Рис. 9
Для получения оценок азимута и угла места формируется матрица
Жт =Е?0 ••• 2ТМ-1]т ,
образованная матрицами 2Тт, т = 0, М -1. Вычисляется корреляционная матрица Я = Жт Ж^,
которая раскладывается по собственным векторам и числам. Количество значимых собственных чисел, превысивших установленный порог, определенный заданным уровнем вероятности ложной тревоги, определяет общее число ИРИ в рассматриваемой полосе. При уголковой АР, образованной двумя линейными АР, расположенными под углом друг к другу (рис. 9), в случае М = 3 (треугольная АР) на одной частоте можно сформировать оценки только одного ИРИ, при М = 5 — двух, при М = 7 — трех ИРИ.
Оценки азимута, угла места и частоты формируются на основе матричных преобразований векторов сигнального подпространства, состоящего из собственных векторов, которым соответствуют собственные числа, превысившие порог [11].
Разработаны и исследованы по реальным записям сигналов алгоритмы при М = 3 и 5.
На рис. 10, а представлен фрагмент спектра записи реальных сигналов на частоте 18 537 кГц1 в зоне перекрытия спектров сигналов двух ИРИ.
Частотно-азимутальная панорама, полученная алгоритмом на основе сигнального подпространства при М = 7, приведена на рис. 10, б, при М = 3 — на рис. 10, в, а алгоритмом на основе методов с частотной селекцией при М = 3 — на рис. 10, г.
Рис. 10, б и в убедительно демонстрируют наличие двух ИРИ, перекрывающихся по спектру. Оценки азимута и угла места сформированы во всей полосе частот сигналов без разрывов.
Средние значения и среднеквадратические отклонения оценок азимута алгоритма на основе выделения сигнального подпространства при М = 7
1 На рис. 10—12 по оси абсцисс отложено отклонение частоты А/ от центральной частоты 18 537 кГц.
Б/ 600 400 200 0
120 0
—2 —1 0 1
2 А/к, кГц
—3 —1.5
0 1.5
б
А/к, кГц
Рис. 10
составляют для первого ИРИ 72° и 1.3° соответственно, а для второго — 152.8° и 2.8°; при М = 3 для первого ИРИ 72° и 1.7°, а для второго — 148° и 4.7°. Оценка азимута второго ИРИ имеет смещение 4°.
Средние значения и среднеквадратические отклонения оценок, сформированные на основе оценок, полученных алгоритмом с частотной селекцией (рис. 10, г), при М = 3 составляют для первого ИРИ 72° и 2.7° соответственно. В зоне перекрытия оценки обоих ИРИ не сформированы. Если бы в алгоритме отсутствовала проверка наличия составляющих двух сигналов в ЧО, то в зоне перекрытия формировалась бы недостоверная оценка азимута ИРИ 145° и 4.1°, не соответствующая ни одному ИРИ.
Алгоритмы на основе сигнального подпространства позволяют также сформировать квадратурные составляющие сигналов, принадлежащих разным ИРИ, перекрывающимся по частоте. Амплитудные спектры оценок сигналов обоих ИРИ, сформированные алгоритмом при М = 7, приведены на рис. 11, а и б, а при М = 3 — на рис. 12, а и б.
Сравнение рис. 11 и 12 показывает, что обоими алгоритмами сформированы практически одинаковые оценки сигнала от первого ИРИ (азимут 72°). На рис. 12, б спектр сигнала второго ИРИ в области перекрытия со спектром первого ИРИ имеет провал, наличие которого обусловлено плохой разделимостью векторов сигнального подпространства, относящихся к первому и второму ИРИ, при М = 3.
На рис. 11, б спектр оценки второго ИРИ провала не содержит, поскольку М = 7.
Б/ -х10 '
Б/ -х10-~
1.5 1.0
0.5 0
—4—2 0
А/к,
кГц Рис. 11
—4—2 0
А/к, кГц
Б/ -х10-7
Б/ х
х10-6
—4—2 0
А/к, кГц
Рис. 12
—4—2 0
б
А/к, кГц
Результаты экспериментального исследования статистическим имитационным моделированием и по реальным записям сигнала свидетельствуют о том, что алгоритмы совместного обнаружения и пеленгования на основе выделения сигнального подпространства:
• при одинаковой АР обеспечивают более высокую точность пеленгования, чем алгоритмы, синтезированные на основе метода с частотной селекцией;
• позволяют сформировать оценки квадратурных составляющих сигнала ИРИ;
• требуют больших вычислительных ресурсов, чем алгоритмы, синтезированные на основе метода с частотной селекцией.
В настоящей статье представлены основные методы и реализованные на их основе алгоритмы совместного обнаружения и пеленгования, предназначенные для панорамного радиомониторинга. Алгоритм при трехэлементной АР на основе метода с частотной селекцией успешно реализован и применяется на практике.
Анализ баз данных оценок УК ИРИ показывает наличие 10 % аномальных оценок. Появляющиеся аномальные измерения обусловлены в основном многолучевым распространением. Разработка алгоритма пеленгования при многолучевом распространении и реализация в реальном времени являются сложной задачей. Поэтому сейчас активно ведутся исследования по вторичной обработке результатов измерений УК ИРИ.
Проводится реализация алгоритма четырех-элементной АР в реальном времени с учетом перекрытия спектров двух сигналов и алгоритма при коммутируемой круговой АР.
2
2
б
а
а
4
2
2
в
г
а
СПИСОК ЛИТЕРАТУРЫ
1. Богданович В. А., Вострецов А. Г. Теория устойчивого обнаружения, различения и оценивания сигналов. 2-е изд., испр. М.: ФИЗМАТЛИТ, 2004. 320 с.
2. Теория обнаружения сигналов / П. С. Акимов, П. А. Бакут, В. А. Богданович и др.; под ред. А. П. Баку-та. М.: Радио и связь, 1984. 440 с.
3. Roy R., Kailath T. ESPRIT-Estimation of Signal Parameters via Rotational Invariance Techniques // IEEE Trans. on Acoust., Speech, and Sign. Proc. 1989. Vol. ASSP-37, № 7. P. 984-995.
4. Lemma A. N., Van der Veen A. J. Analysis of Joint Angle-Frequency Estimation Using ESPRIT // IEEE Trans. on Sign. Proc. 2003. Vol. SP-51, № 5. P. 1264-1283.
5. Богданович В. А., Стенюков Н. С., Шевченко М. Е. Алгоритмы обнаружения сигналов в системах частотного радиомониторинга // Конф. "Научные, инженерные и производственные проблемы создания технических средств мониторинга электромагнитного поля", СПб., 31 мая-2 июня 2005 г. СПб.: Изд-во СПбИИ РАН "Нестор-История", 2005. С. 33.
6. Lemma A. N., Deprettere E. F., Veen A. J. Experimental analysis of antenna coupling for high-resolution DOA estimation algorithms // Proc. of 2nd IEEE Workshop on Signal Processing Advances in Wireless Communications. Annapolis, 9-12 May 1999. Piscataway: IEEE, 1999. P. 362-365.
7. Гутин В. С., Шевченко М. Е. Возникновение аномальных оценок направлений источников радиоизлучения в условиях многолучевого распространения // 3-я Всерос. конф. "Радиоэлектронные средства получения, обработки и визуализации информации" (РСПОВИ-2013), Смоленск, 26-28 июня 2013 г. / РНТОРЭС им. А. С. Попова. М., 2013. С. 166-169.
8. Шевченко М. Е., Чемаров А. О. Обнаружение и оценивание параметров источников радиоизлучения в широкой полосе обзора. СПб.: Изд-во СПбГЭТУ "ЛЭТИ", 2011. 136 с.
9. Шевченко М. Е., Задирако Д. О. Совместное обнаружение и пеленгование источников радиоизлучения с помощью четырехэлементной антенной решетки при перекрытии спектров соседних сигналов // 5-я Всерос. конф. (с междунар. участием) "Радиоэлектронные средства получения, обработки и визуализации" (РСПОВИ-2015). М., 28-30 окт. 2015 г. / РНТОРЭС им. А. С. Попова. М., 2015. С. 73-77.
10. Шевченко М. Е., Малышев В. Н., Файзуллина Д. Н. Совместное обнаружение и пеленгование с использованием коммутируемой антенной решетки // Изв. вузов России. Радиоэлектроника. 2015. Вып. 5. С. 33-38.
11. Шевченко М. Е. Алгоритмы совместного обнаружения и пеленгования на основе методов сигнальных подпространств. СПб.: Изд-во СПбГЭТУ "ЛЭТИ", 2013. 160 с.
M. E. Shevchenko, D. O., Zadirako, D. N. Faizullina, V. N. Malyshev Saint Petersburg Electrotechnical University "LETI"
N. S. Stenyukov, M. S. Shmyrin JSC «SRI "Vector"» (Saint Petersburg)
Methods and Algorithms of Wide Band Radiomonitoring under Small Element Antenna Arrays
Methods of panoramic radiomonitoring and developed on their basis algorithms for joint detection and direction finding in a wide band of frequencies are presented. The block diagrams and panoramas illustrating the results of the algorithms are given.
Antenna array, joint detection and direction finding, azimuth and elevation angle estimation, KV-VHF-band, radio sources, overlapping spectra, phase direction finding
Статья поступила в редакцию 11 апреля 2016 г.
