Алгоритмы DOA-оценки для GNSS-антиджамминга: разреженные решётки, CRB и диаграммообразование

Разбираем, как работают алгоритмы оценки направления прихода (DOA) помех в GNSS, зачем нужны разреженные антенные решётки, как граница Крамера-Рао и пространственный коэффициент корреляции повышают точность, и какие схемы диаграммообразования применяются на практике.

Что такое DOA-оценка для GNSS-антиджамминга?

Обнаружение помехи — первый шаг к её подавлению. Система оценивает направление прихода (Direction of Arrival, DOA) сигнала глушилки с помощью обработки антенных решёток, после чего применяет пространственный фильтр или формирует «нуль» диаграммы направленности, чтобы минимизировать влияние помехи. Полезный сигнал очищается до того, как он попадёт в GNSS-приёмник. Алгоритмы DOA-оценки опираются на функцию мощности P(theta, phi) — пространственный оценщик, зависящий от полярных угловых координат источника помехи.

Точность такой оценки напрямую ограничена разрешающей способностью решётки и числом доступных отсчётов. В реальных условиях GNSS-помеха может быть широкополосной, импульсной или узкополосной, и каждая из этих ситуаций по-разному влияет на спектральные и параметрические методы. Именно поэтому современные исследования сосредоточены на алгоритмах, способных работать при малом числе снимков (snapshots) и низком отношении сигнал/шум.

В 2012 году исследователи Сеульского национального университета (Hong-Ki Song, Sung-Yong Lee, Seung-Hoo Lee, Jung-Beom Kim, Changdon Kee) показали, что корректная пространственная фильтрация на уровне антенной решётки способна сохранить точность позиционирования даже при активном воздействии. Позже коллектив Национального университета обороны Китая (Xianghao Li, Zukun Lu, Feiqiang Chen и соавторы) развил методы, ориентированные на реальные GNSS-сигналы с их специфическими модуляциями.

Типовая цепочка обработки включает четырёхканальный антенный тракт: четырёхэлементная квадратурная фазированная решётка и четырёхканальный когерентный фронтенд. Такая схема распознаёт помеху, реализует пространственный фильтр и восстанавливает сигнал для коммерческого или программно-определяемого приёмника (COTS/SDR). Это делает DOA-оценку не абстрактной теорией, а практическим инструментом защиты навигации.

Как работают алгоритмы DOA-оценки: спектральные и параметрические методы

Все методы оценки DOA помехового сигнала делятся на два больших класса: спектральное оценивание и параметрическое оценивание. Спектральные методы строят оценку на основе функции мощности — например, классического диаграммообразования (метод Бартлетта) или MVDR (Minimum Variance Distortionless Response). Параметрические методы исходят из модели сигнала и оценивают её параметры, что даёт более высокое разрешение при малом числе отсчётов.

Традиционные алгоритмы DOA-оценки часто определяют целевую функцию на основе евклидовой метрики расстояния в пространстве (T Long, 2026). Однако предложенный в той же работе алгоритм демонстрирует превосходную оценку DOA при менее чем 50 снимках и отношении сигнал/шум ниже −5 дБ. Это принципиально важно для GNSS, где помеха может быть слабой по мощности, но достаточной для срыва слежения за сигналом.

Сравнение двух подходов удобно свести в таблицу по ключевым характеристикам.

Выбор между классами — это компромисс между вычислительной сложностью, устойчивостью к нештатным помехам и разрешающей способностью. Жёсткая спецификация модели снижает гибкость и робастность методов при работе с неноминальными помехами (Inside GNSS, 2017). Поэтому на практике часто комбинируют спектральные и параметрические подходы, адаптируя их под конкретную радиочастотную обстановку.

Разреженные решётки и разреженная выборка для подавления помех

Разреженные антенные решётки и разреженная выборка играют заметную роль в антиджамминге GNSS (MG Amin, цитируется 201 раз). Их главное преимущество — снижение аппаратных и вычислительных затрат по сравнению с полной однородной круговой решёткой (Uniform Circular Array, UCA). Вместо всех элементов используется уменьшенный набор, что удешевляет фронтенд и упрощает обработку.

Один из подходов — переконфигурация разреженной неоднородной решётки: из UCA выбирается сокращённое число элементов. Оптимизационная модель при этом является квадратично-дробной и решается методом полуопределённой релаксации (Semi-Definite Relaxation, SDR). Такой подход позволяет сохранить точность и робастность DOA-оценки при меньшей сложности системы.

Взаимные (coprime) решётки тоже рассматриваются как инструмент защиты от спуфинга — в работе Y Zhao 2022 года копрайм-массивы применяются для обнаружения подмены сигнала. А M Moussa в 2019-м разбирал узкополосные GPS-глушилки и то, как они влияют на приёмники. Всё это наводит на мысль, что разреженные конфигурации — вовсе не вынужденный компромисс, а вполне осознанный инженерный выбор.

Разреженная выборка позволяет дополнительно сократить объём данных без потери ключевой информации о направлении прихода. Однако разреженность накладывает ограничения: при слишком малом числе элементов падает разрешающая способность, и алгоритм может «склеить» близкие источники помех. Поэтому выбор конфигурации всегда балансирует между стоимостью и точностью.

Оптимизация точности: граница Крамера-Рао и пространственный коэффициент корреляции

Минимизация границы Крамера-Рао (Cramér-Rao Bound, CRB) используется как целевая функция для обеспечения точности DOA-оценки. Пространственный коэффициент корреляции (Spatial Correlation Coefficient, SCC) добавляется как ограничение для подавления боковых лепестков (R LANG, 2025, 20 цитирований). Полученная квадратично-дробная модель решается методом полуопределённой релаксации (SDR).

Такой подход позволяет одновременно добиваться высокой точности и снижать уровень ложных пиков в спектре оценки. Боковые лепестки особенно опасны в GNSS: они могут привести к ложному обнаружению помехи в направлении, где её нет, и к подавлению полезного сигнала. Ограничение SCC действует как регуляризатор, удерживающий диаграмму в допустимых рамках.

Работы Rongling LANG, Hao XU, Fei GAO, Zewen TANG, Zhipeng WANG и Amir HUSSAIN (Beihang University и Edinburgh Napier University) развивают это направление, объединяя теоретические границы точности с практическими конструкциями решёток. Их результаты показывают, что совместная оптимизация CRB и SCC даёт устойчивое решение даже при неблагоприятных условиях приёма.

Для инженера это означает, что проектирование антенной системы перестаёт быть «подбором на глаз»: появляется формальный критерий, по которому можно сравнивать конфигурации. Минимизация CRB задаёт нижнюю границу достижимой точности, а SCC гарантирует, что решение не будет разрушено боковыми лепестками.

Сжатое зондирование и методы разреженной реконструкции

Алгоритм CS-SFD оценивает ковариационную матрицу помех через разреженную реконструкцию и применяет диагональную загрузку (diagonal loading) для робастности. Процедура включает получение одиночного снимка (single shot) из пространственно-временной решётки и его векторизацию, дискретизацию направлений и последующую оценку ковариационной матрицы помех.

Диагональная загрузка — важный элемент устойчивости: она стабилизирует обращение матрицы в условиях малого числа отсчётов и коррелированных помех. Без неё разреженная реконструкция может давать нестабильные результаты, особенно когда источников помех несколько.

Методы сжатого зондирования (compressed sensing) в целом позволяют обойтись меньшим числом измерений, чем классические методы, что напрямую снижает требования к вычислительным ресурсам приёмника. Это критично для встраиваемых GNSS-платформ, где энергопотребление и тактовая частота ограничены.

Практическая реализация CS-SFD требует аккуратной настройки параметров регуляризации. Слишком слабая загрузка ведёт к неустойчивости, слишком сильная — к размыванию пиков и потере разрешения. Баланс подбирается под конкретный сценарий помеховой обстановки.

Диаграммообразование и формирование нулей: Бартлетт и MVDR

Диаграммообразование и формирование нулей — практический выход всей цепочки DOA-оценки. Классическое диаграммообразование (метод Бартлетта) и MVDR относятся к спектральному оцениванию. Бартлетт прост в реализации, но имеет широкий главный лепесток и высокий уровень боковых лепестков. MVDR минимизирует мощность выходного сигнала при единичном усилении в заданном направлении, что даёт более глубокие нули в сторону помех.

MVDR особенно эффективен, когда число помех меньше числа антенных элементов. В GNSS-сценариях это обычно выполняется: четыре элемента решётки против одной-двух глушилок. Однако при большом числе источников или при их коррелированности эффективность MVDR падает.

Существует и альтернатива — Constrained Adaptive Monopulse (CAMP), применяемый в схеме с поддержкой системы опорной ориентации (Attitude Reference System, ARS). CAMP позволяет измерять рыскание (yaw) напрямую, без итерационного поиска, что существенно снижает вычислительную нагрузку.

Выбор метода диаграммообразования зависит от приоритетов: максимальная простота, максимальная глубина подавления или минимальная задержка обработки. На практике часто используется комбинация — грубая оценка спектральными методами и точная настройка параметрическими.

Схема с поддержкой ARS: сокращение вычислений и отслеживание ориентации

Схема с поддержкой ARS сокращает время вычислений примерно в 18 раз по сравнению с высокооптимизированным одномерным поиском максимального правдоподобия, сохраняя высокую точность оценки ориентации при четырёх глушилках (ION 2026 Pacific PNT Meeting, 13–16 апреля 2026, Гонолулу, Гавайи, стр. 833–844). Авторы работы — Chuanrui Wang, Xiaowei Cui, Gang Liu, Mingquan Lu.

Ключевая идея: трёхмерная оценка ориентации сводится к одномерному отслеживанию рыскания. Локально изменение угла рыскания эквивалентно одномерному сдвигу DOA, что позволяет измерять рыскание напрямую с помощью CAMP без итерационного поиска. Это радикально упрощает вычисления.

Дополнительно применяется мультиспутниковое объединение на основе обратно-дисперсионного взвешивания (inverse-variance weighting), которое подавляет выбросы, вызванные помехами или плохой геометрией созвездия. Такая схема устойчива к отдельным сбоям измерений.

Для высокоточных применений, таких как железнодорожные системы ERTMS, важны вероятности пропуска обнаружения (Pmd) и частота отказов (FR). Исследования Cosimo Stallo и коллег (Radiolabs, University of Roma 3, University of Rome Tor Vergata, Hitachi Rail STS, RFI, European Space Agency, Qascom) показывают, как DOA-оценка встраивается в контур управления безопасностью.

Аппаратные спецификации и параметры типовой системы

Типовая система антиджамминга включает четыре антенны и четыре встроенных усилителя. Поддерживаются трёхдиапазонные GPS L1/L2/L5, GLONASS G1/G2/G5, BeiDou B1/B2, Galileo E1/E5, а также L-диапазонные сервисы коррекции. Характеризация решётки проводится в сферической ближнепольной системе StarLab 18 ГГц.

Моделирование выполняется при менее чем 50 снимках и отношении сигнал/шум ниже −5 дБ. Время вычислений сокращается примерно в 18 раз по сравнению с одномерным поиском максимального правдоподобия. Схема с ARS выдерживает до четырёх глушилок при отслеживании ориентации.

Сводка ключевых параметров приведена в таблице.

Отдельного внимания заслуживает временная антиджамминговая фильтрация (TDAJ). Такие фильтры могут вызывать искажение сигнала и смещение дальномерных измерений. Для компенсации применяется предварительный фильтр оптимального порядка, адаптивно корректирующий искажения. Это важно для высокоточного позиционирования, где ошибка в несколько наносекунд уже критична.

Компромиссы, ограничения и вопросы безопасности

Разреженные решётки снижают аппаратные и вычислительные затраты по сравнению с полной UCA. Однако жёсткая спецификация модели уменьшает гибкость и робастность методов при работе с неноминальными помехами (Inside GNSS, 2017). Аналоговая эквализация обладает низкой точностью и ограниченной применимостью, а цифровая калибровка требует высокой вычислительной сложности.

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

С точки зрения регулирования: глушилки запрещены в большинстве стран, но при этом дёшевы и легкодоступны. Простые устройства способны нарушить работу GNSS-сервисов на расстоянии нескольких километров. Инциденты с GNSS-помехами фиксировались в разных странах, включая активность со стороны Северной Кореи.

Для критических инфраструктур, таких как ERTMS, требования безопасности включают вероятность пропуска обнаружения (Pmd) и частоту отказов (FR). Это означает, что алгоритмы DOA-оценки должны не только работать точно, но и быть сертифицируемыми в рамках формальных процедур обеспечения безопасности.

Практические шаги внедрения DOA-оценки

Первый шаг — получить одиночный снимок из пространственно-временной решётки и векторизовать его. Затем дискретизировать направления и оценить ковариационную матрицу помех через разреженную реконструкцию. После этого применить диагональную загрузку для обеспечения робастности.

Для переконфигурации разреженной решётки: выбрать сокращённый набор элементов, минимизировать CRB, ограничить SCC и решить задачу через SDR. Для отслеживания рыскания с поддержкой ARS: использовать крен и тангаж от ARS, свести задачу к одномерному рысканию, применить CAMP и объединить мультиспутниковые измерения.

Важно контролировать качество калибровки антенного тракта. Рассогласование амплитуд и фаз между каналами напрямую ухудшает точность DOA и может свести на нет преимущества даже хорошего алгоритма. Регулярная проверка в ближнепольной системе помогает поддерживать характеристики решётки.

Наконец, следует помнить о компромиссе между сложностью и устойчивостью. Чем больше априорных предположений заложено в модель, тем лучше результат в номинальных условиях и тем хуже — при нештатных помехах. Гибридные схемы, сочетающие спектральные и параметрические методы, обычно дают наиболее сбалансированное решение.

Часто задаваемые вопросы

<h3>Как DOA-оценка поддерживает GNSS-антиджамминг?</h3><p>Обнаружение помехи — первый шаг к её подавлению. Система оценивает направление прихода сигнала глушилки с помощью обработки антенной решётки, после чего применяет пространственный фильтр или формирует нуль диаграммы, чтобы минимизировать влияние помехи. Полезный сигнал очищается до попадания в GNSS-приёмник, что сохраняет возможность позиционирования.</p>

<h3>Какие основные категории методов DOA-оценки существуют?</h3><p>Методы оценки DOA помехового сигнала делятся на спектральное и параметрическое оценивание. Диаграммообразование и формирование нулей, включая классическое диаграммообразование (Бартлетт) и MVDR, относятся к спектральному оцениванию. Параметрические методы опираются на модель сигнала и дают более высокое разрешение при малом числе отсчётов.</p>

<h3>Зачем использовать разреженные решётки для DOA-оценки помех GNSS?</h3><p>Разреженные решётки и разреженная выборка снижают аппаратные и вычислительные затраты. Из однородной круговой решётки можно выбрать сокращённый набор элементов, снизив сложность при сохранении точности и робастности DOA-оценки. Это особенно важно для встраиваемых GNSS-платформ с ограниченными ресурсами.</p>

<h3>Какую роль играет граница Крамера-Рао в DOA-оценке?</h3><p>Минимизация границы Крамера-Рао (CRB) используется как целевая функция для обеспечения точности DOA-оценки. Пространственный коэффициент корреляции (SCC) добавляется как ограничение для подавления боковых лепестков, а полученная квадратично-дробная модель решается методом полуопределённой релаксации (SDR). Это даёт устойчивое решение при неблагоприятных условиях приёма.</p>

Часто задаваемые вопросы

Как DOA-оценка поддерживает GNSS-антиджамминг?

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

Какие основные категории методов DOA-оценки существуют?

Методы оценки DOA помехового сигнала делятся на спектральное и параметрическое оценивание. Диаграммообразование и формирование нулей, включая классическое диаграммообразование (Бартлетт) и MVDR, относятся к спектральному оцениванию. Параметрические методы опираются на модель сигнала и дают более высокое разрешение при малом числе отсчётов.

Зачем использовать разреженные решётки для DOA-оценки помех GNSS?

Разреженные решётки и разреженная выборка снижают аппаратные и вычислительные затраты. Из однородной круговой решётки можно выбрать сокращённый набор элементов, снизив сложность при сохранении точности и робастности DOA-оценки. Это особенно важно для встраиваемых GNSS-платформ с ограниченными ресурсами.

Какую роль играет граница Крамера-Рао в DOA-оценке?

Минимизация границы Крамера-Рао (CRB) используется как целевая функция для обеспечения точности DOA-оценки. Пространственный коэффициент корреляции (SCC) добавляется как ограничение для подавления боковых лепестков, а полученная квадратично-дробная модель решается методом полуопределённой релаксации (SDR). Это даёт устойчивое решение при неблагоприятных условиях приёма.