Практическая реализация пространственно-временной фильтрации спутниковых навигационных сигналов в реальном времени

Полный текст

Введение

Сигналы глобальных спутниковых навигационных систем обладают повышенной уязвимостью к помехам различной природы ввиду своей слабой мощности вблизи поверхности Земли. Также в последние годы целенаправленное излучение помех начинает представлять все большую угрозу для нормального функционирования технических средств, работа которых полагается на спутниковые навигационные системы. Этим обуславливается высокая актуальность задачи построения помехозащищенных приемников спутниковых навигационных сигналов [1; 2].

Одним из часто используемых методов подавления помех является пространственно-временная фильтрация. Этот метод представляет собой обобщение техники пространственной обработки, увеличивающее количество степеней свободы цифрового фильтра путем добавления в процесс его построения расширенной временной составляющей сигнала. Данная методика позволяет более эффективно бороться с присутствующими в сигнале помехами, но при этом кратно увеличивает вычислительную сложность алгоритма фильтрации [3]. Вследствие этого практическая реализация пространственно-временных алгоритмов в реальном времени обычно осуществляется с использованием программируемых логических интегральных схем (ПЛИС) [2].

Целью данной работы являются построение пространственно-временного фильтра, его практическая реализация в реальном времени в ПЛИС в составе помехозащищенного навигационного приемника (ПНП) и проверка качества работы полученной спутниковой навигационной системы.

1. Постановка задачи фильтрации. Построение алгоритма обработки

Задача адаптивной пространственно-временной фильтрации спутниковых навигационных сигналов имеет следующую постановку [4–6]. В горизонтальной плоскости расположена антенная решетка (АР) с $N$ антенными элементами. Над решеткой находятся различные источники полезного сигнала (спутники) и помех. Приходящий на АР входной сигнал представляет собой аддитивную смесь полезного сигнала, помех и шума, причем уровень помехи выше уровня полезного сигнала, который, в свою очередь, выше уровня шума. Требуется построить адаптивный пространственно- временной фильтр, выходной сигнал которого будет представлять собой полезный сигнал, выделенный из смеси с помехой и шумом.

За основу пространственно-временного фильтра была взята следующая процедура пространственной обработки [7]. Пусть $\left(t\right)\in {ℂ}^N$ – входной сигнал фильтра, определенный в дискретные моменты времени $t=\mathrm{1,2,}\dots ,$ $x_1\left(t\right)\in ℂ$ – первый входной канал, $\left(t\right)\in {ℂ}^{N-1}$ – входной сигнал без первого канала. Тогда выходной сигнал $y\left(t\right)\in ℂ$ пространственного фильтра можно построить для всех $t$ как

$y\left(t\right)=x_1\left(t\right)+\left(\left(t\right),\left(t\right)\right),$ (1)

где $\left(\cdot ,\cdot \right)$ – скалярное произведение, $w\left(t\right)=-{\widehat{R}}^{-1}\left(t\right)r\left(t\right)$ – веса фильтра, где $\widehat{R}\left(t\right)=\tilde{x}\left(t\right){\tilde{x}}^{*}\left(t\right)$ – выборочная оценка корреляционной матрицы входного сигнала; верхний индекс * – эрмитово сопряжение; $r\left(t\right)={\overline{x}}_1\left(t\right)\tilde{x}\left(t\right),$ $\overline{z}$ – число, комплексно сопряженное с $z.$

Фильтр (1) можно сделать пространственно-временным, если вместо $\tilde{x}\left(t\right)$ использовать вектор, полученный конкатенацией векторов $\tilde{x}\left(t\right), \tilde{x}(t+1),$$\dots , \tilde{x}(t+T-1)$ для некоторого натурального количества $T$ временных отводов фильтра. Для того чтобы сделать построенную процедуру обработки более подходящей для реализации в ПЛИС, можно вычислять веса фильтра один раз в $K\in \mathbb{N}$ отсчетов с помощью накопленных суммированием матриц ${\widehat{R}}_K$ и $r_K:$

$\begin{array}{l}{w}_K=-{\widehat{R}}_K^{-1}{r}_K,\\ {\widehat{R}}_K=\sum _{t=1}^K\tilde{x}\left(t\right){\tilde{x}}^{*}\left(t\right),\\ r_K=\sum _{t=1}^K{\overline{x}}_1\left(t\right){\tilde{x}}^{*}\left(t\right).\end{array}$ (2)

Веса $w_K$ применяются для вычисления выхода $y\left(t\right)$ для всех $t\in \left[\mathrm{1,}K\right],$ а затем пересчитываются заново для каждого отрезка времени из $K$ последовательных отсчетов. Такая модификация позволяет уменьшить вычислительную нагрузку на систему и увеличить численную устойчивость операции обращения корреляционной матрицы и точность ее выборочного приближения [2].

2. Использованная аппаратура

Построенный пространственно-временной фильтр был реализован в ПНП, устроенном следующим образом. Входными сигналами ПНП служили спутниковые сигналы, принятые на активную АР, представляющую собой объединение в одном корпусе двух восьмиэлементных кольцевых эквидистантных АР с разными радиусами, соответствующими диапазонам L1 и L2. После прихода на АР спутниковые сигналы каждого из диапазонов проходили через каскад «полосовой фильтр (ПФ) – усилитель (У) – ПФ» и сумматор, после чего по кабелю передавались в ПНП. Сигнал, поступивший на вход ПНП, проходил через делитель мощности, после чего результирующие сигналы проходили через каскады «ПФ – У» и смесители, в которых смешивались с гетеродинами, соответствующими каждому из сигнальных диапазонов. Затем сигналы снова проходили через каскады «ПФ – У – ПФ» промежуточной частоты и аналого-цифровые преобразователи (АЦП), после чего подавались на вход ПЛИС AMD Xilinx Kintex UltraScale XCKU035, где выполнялась пространственно-временная обработка. После этого каждый из сигналов попадал в передающий тракт одного из двух трансиверов Analog Devices AD9463, где выполнялась демодуляция сигнала, который затем передавался на вход приемника u-blox ZED-F9P, где происходило решение навигационной задачи. Схема АР и ПНП представлена на рис. 1.

 

Рис. 1. Схема ПНП и АР

Fig. 1. Antenna array and jamming-protected navigational receiver scheme

 

Для данной реализации время накопления $K$ для выборочного приближения корреляционной матрицы было принято равным 1250 отсчетам при обработке сигнала GPS и 5000 отсчетам для сигналов ГЛОНАСС. Количество $T$ временных отводов фильтра было принято равным 5.

В табл. 1 приведены данные об использовании алгоритма фильтрации, реализованного на ПЛИС, различных ее ресурсов. Рассматривались соответствующие значения для таких частей ПЛИС, как таблицы поиска (англ. lookup table, LUT), связанная с ними оперативная память (англ. LUTRAM), триггеры (англ. flip-flop, FF), блочная оперативная память (англ. BRAM) и блоки цифровой обработки сигналов (англ. digital signal processing, DSP).

Таблица 1. Использование алгоритма фильтрации ресурсов ПЛИС

Table 1. FPGA resources used by the filtering algorithm

Ресурс	Доступно	Использовано	Использовано, %
LUT	203128	185216	91,18
LUTRAM	112800	47142	41,79
FF	406256	368358	90,67
BRAM	540	405	75
DSP	1700	1435	84,41

 

Испытания системы приема и обработки спутниковых сигналов проводились применительно к сигналам GPS L1 и ГЛОНАСС L1 и L2. В качестве устройства генерации и излучения некоррелированных широкополосных помех использовалась плата трансивера (ПТ) с системой на кристалле AMD Xilinx Zynq UltraScale+ RFSoC, восемью АЦП и восемью цифро-аналоговыми преобразователями (ЦАП). При работе с сигналами GPS использовались помехи полосой 5 МГц с сигналами ГЛОНАСС – 20 МГц.

3. Эксперименты

Перед проведением испытаний алгоритм фильтрации был протестирован путем компьютерного моделирования с использованием записей сигналов GPS L1 с добавлением различного количества некоррелированных помех. При этом измерялись такие характеристики работы алгоритма фильтрации, как коэффициент подавления (КП) помехи, равный отношению мощностей входного и выходного сигналов [7], выраженному в децибелах, а также отношение сигнал/шум (ОСШ) в выходном сигнале. Цифровая фильтрация с измерением КП была реализована в виде обработки программой на языке C++. Значения ОСШ были получены в результате работы цифрового приемника спутниковых сигналов SoftGNSS [8], реализованного в среде MATLAB. Полученные значения характеристик работы модели приведены в табл. 2 в соответствующем столбце.

 

Таблица 2. Результаты моделирования и испытаний ПНП

Table 2. Modeling and jamming-protected navigational receiver tests results

	КП, дБ				ОСШ, дБ
Количество помех	GPS (модель)	GPS	ГЛОНАСС L1	ГЛОНАСС L2	GPS (модель)	GPS	ГЛОНАСС L1	ГЛОНАСС L2
1	70	54	44	52	37	39	41	42
2	70	51	38	50	36	40	38	41
3	69	50	34	46	39	41	33	42
4	69	48	29	40	43	38	32	40
5	67	40	26	36	38	36	32	37
6	67	38	22	30	41	35	30	35
7	64	35	20	27	38	33	29	34

 

После моделирования были проведены испытания ПНП для сигналов всех типов и диапазонов с измерением тех же характеристик работы. Эксперименты по приему и обработке спутниковых сигналов проводились на крыше здания, под открытым небом. На расстоянии 3 м по кругу от ПНП с АР было размещено 7 излучающих антенн, равноудаленных друг от друга и подключенных к выходам ЦАП ПТ. Схема установки для испытаний изображена на рис. 2.

 

Рис. 2. Схема установки для испытаний. Обозначения: 1 – АР, 2 – ПНП, 3 – ПТ, 4 – излучающие антенны

Fig. 2. Tests setting scheme. Designations: 1 – antenna array, 2 – jamming-protected navigational receiver, 3 – transceiver circuit board, 4 – radiating antennae

 

Рассматривалось несколько сигнально-помеховых ситуаций, отличавшихся количеством излучаемых помех (от 1 до 7). Во всех случаях после фильтрации через достаточно короткое время приемник входил в режим Fixed, то есть получал навигационное решение. Результаты испытаний приведены в табл. 2 в соответствующих столбцах.

Для каждой сигнально-помеховой ситуации проводилось по три эксперимента, в которых обрабатывались выборки длиной ${10}^5,$ ${10}^6$ и ${10}^7$ отсчетов. Если задать желаемую точность как $\epsilon ={10}^{-5},$ то нетрудно убедиться, что выполняется неравенство

$\frac{\epsilon \sqrt{N_e}}{\sigma }\le \tau \left(\alpha ,\nu \right),$

где $N_e=3$ – количество экспериментов; $\sigma =\mathrm{5,4745}\cdot {10}^6$ – стандартное отклонение набора длин выборок; $\tau =\mathrm{4,303}$ – значение распределения Стьюдента при $\alpha =\mathrm{0,95,}$ $\nu =N_e-1.$ Это означает, что данный набор длин выборок удовлетворяет критерию Стьюдента и является достаточным для получения достоверных результатов экспериментов [9].

Из результатов моделирования и экспериментов с ПНП можно сделать следующие выводы. В модельном случае построенный алгоритм пространственно-временной обработки позволяет получить достаточно большие значения КП и ОСШ, что говорит о хорошем качестве работы алгоритма в упрощенных модельных условиях. Значения КП, полученные при испытаниях ПНП, заметно ниже модельных данных, однако значения ОСШ при этом получились лишь незначительно ниже. В совокупности с фактом получения навигационного решения во всех экспериментах это означает, что построенный ПНП справляется с поставленной задачей. Также можно отметить, что динамика изменения значений КП в зависимости от количества помех в реальном и модельном случаях гораздо более схожая, чем аналогичная динамика изменения значений ОСШ.

Заключение

В данной работе было проведено построение адаптивного алгоритма пространственно-временной обработки спутниковых навигационных сигналов и описание его практической реализации в ПЛИС в составе ПНП, а также приведены результаты моделирования работы полученного алгоритма и испытаний ПНП в различных сигнально-помеховых ситуациях. На основе результатов экспериментов были сделаны выводы о работоспособности и эффективности построенного алгоритма фильтрации как в модельном варианте, так и в составе реального ПНП, что позволяет рекомендовать использование данного алгоритма при создании систем помехозащищенной спутниковой навигации в будущем.

Об авторах

Евгений Иванович Глушанков

Санкт-Петербургский государственный университет телекоммуникаций имени профессора М.А. Бонч-Бруевича

Email: glushankov57@gmail.com
ORCID iD: 0000-0003-4148-3208

доктор технических наук, профессор кафедры радиосистем и обработки сигналов

Россия, Санкт-Петербург

Владимир Игоревич Царик

ООО «Эйртэго»

Автор, ответственный за переписку.
Email: wladimirzarik@mail.ru
ORCID iD: 0000-0003-3428-9976

магистр по направлению «Прикладная математика и информатика» (Санкт-Петербургский государственный университет, 2020 г.), соискатель ученой степени кандидата технических наук, ведущий инженер 

Россия, Санкт-Петербург

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

Дополнительные файлы