In this paper, we investigate the potential of reconfigurable intelligent surfaces (RISs) in facilitating passive/device-free three-dimensional (3D) drone localization within existing cellular infrastructure operating at millimeter-wave (mmWave) frequencies and employing multiple antennas at the transceivers. The developed localization system operates in the bi-static mode without requiring direct communication between the drone and the base station. We analyze the theoretical performance limits via Fisher information analysis and Cramér Rao lower bounds (CRLBs). Furthermore, we develop a low-complexity yet effective drone localization algorithm based on coordinate gradient descent and examine the impact of factors such as radar cross section (RCS) of the drone and training overhead on system performance. It is demonstrated that integrating RIS yields significant benefits over its RIS-free counterpart, as evidenced by both theoretical analyses and numerical simulations.
ID статьи : 2404.14879Название : Device-Free 3D Drone Localization in RIS-Assisted mmWave MIMO NetworksАвторы : Jiguang He, Charles Vanwynsberghe, Hui Chen, Chongwen Huang, Aymen FakhreddineКлассификация : eess.SP (обработка сигналов)Дата публикации : апрель 2024 г. (препринт arXiv, обновлено 2 января 2025 г.)Ссылка на статью : https://arxiv.org/abs/2404.14879 В данной работе исследуется потенциал переконфигурируемых интеллектуальных поверхностей (RIS) для реализации пассивной/бесустройственной трёхмерной (3D) локализации беспилотных летательных аппаратов (БПЛА) в инфраструктуре многоантенных сотовых базовых станций в диапазоне миллиметровых волн (mmWave). Разработанная система локализации работает в бистатическом режиме без необходимости прямой связи между БПЛА и базовой станцией. Теоретические пределы производительности анализируются с использованием информационного анализа Фишера и границы Крамера-Рао (CRLB). Кроме того, разработан алгоритм локализации БПЛА с низкой сложностью на основе координатного градиентного спуска, и исследовано влияние таких факторов, как эффективная площадь рассеяния (RCS) БПЛА и затраты на обучение на производительность системы. Как теоретический анализ, так и численное моделирование демонстрируют значительные преимущества интегрированной системы RIS по сравнению с системой без RIS.
Основная проблема, которую решает данная работа, заключается в реализации пассивной трёхмерной локализации летающих БПЛА в существующей сотовой инфраструктуре, то есть в точном определении их пространственного положения без активной связи БПЛА с базовой станцией.
Требования безопасности : с распространением применения БПЛА растёт необходимость в обнаружении и локализации несанкционированных или вредоносных БПЛАПовторное использование инфраструктуры : полное использование существующей инфраструктуры сотовых сетей 5G/6G, избежание развёртывания специализированных радиолокационных системИнтеграция технологий : важный сценарий применения парадигмы интегрированной связи и восприятия (ISAC)Традиционные радиолокационные системы : требуют специального оборудования, высокие затраты и ограниченный диапазон покрытияАктивные системы локализации : требуют связи между целью и базовой станцией, неприменимы к некооперативным целямСуществующие пассивные системы : ограниченная производительность в диапазоне mmWave, недостаточная точность локализацииИспользование RIS для усиления путей распространения сигнала и повышения производительности пассивной локализации Реализация высокоточной трёхмерной локализации в системах mmWave MIMO Разработка практических алгоритмов с низкой сложностью, приближающихся к теоретическим пределам производительности Проектирование системы : предложена архитектура системы пассивной локализации БПЛА на основе RIS-ассистированного mmWave MIMOТеоретический анализ : анализ теоретических пределов производительности системы с использованием матрицы информации Фишера и CRLBРазработка алгоритма : разработан трёхэтапный алгоритм локализации с низкой сложностью на основе координатного градиентного спуска (CGD)Оценка производительности : комплексный анализ влияния ключевых параметров (RCS, затраты на обучение) на производительность системыПроверка преимуществ : теоретическое и экспериментальное подтверждение значительных преимуществ RIS-ассистированной системы по сравнению с традиционными методамиВходные данные : многовременные принятые сигналы на приёмном устройстве пользователя (UE)
Выходные данные : трёхмерные координаты положения БПЛА ( x D , y D , z D ) (x_D, y_D, z_D) ( x D , y D , z D ) Ограничения : отсутствие прямой связи между БПЛА и базовой станцией, система работает в бистатическом режиме
Система состоит из трёх ключевых компонентов:
Базовая станция mmWave (BS) : оснащена равномерной плоской решёткой из M B = M B , y × M B , z M_B = M_{B,y} \times M_{B,z} M B = M B , y × M B , z антеннПереконфигурируемая интеллектуальная поверхность (RIS) : содержит M R = M R , x × M R , y M_R = M_{R,x} \times M_{R,y} M R = M R , x × M R , y отражающих элементовПриёмное устройство пользователя (UE) : оснащено M U M_U M U антеннамиИспользуется модель канала Салеха-Валенсуэлы, учитывающая разреженные характеристики рассеяния в диапазоне mmWave. Ключевые каналы включают:
Канал BS-RIS H 1 ∈ C M R × M B \mathbf{H}_1 \in \mathbb{C}^{M_R \times M_B} H 1 ∈ C M R × M B :
H 1 = e − j 2 π d 1 / λ ρ 1 α x ( θ r , 1 , ϕ r , 1 ) ⊗ α y ( θ r , 1 , ϕ r , 1 ) ( α y ( θ t , 1 , ϕ t , 1 ) ⊗ α z ( ϕ t , 1 ) ) H \mathbf{H}_1 = \frac{e^{-j2\pi d_1/\lambda}}{\sqrt{\rho_1}} \boldsymbol{\alpha}_x(\theta_{r,1}, \phi_{r,1}) \otimes \boldsymbol{\alpha}_y(\theta_{r,1}, \phi_{r,1}) \left(\boldsymbol{\alpha}_y(\theta_{t,1}, \phi_{t,1}) \otimes \boldsymbol{\alpha}_z(\phi_{t,1})\right)^H H 1 = ρ 1 e − j 2 π d 1 / λ α x ( θ r , 1 , ϕ r , 1 ) ⊗ α y ( θ r , 1 , ϕ r , 1 ) ( α y ( θ t , 1 , ϕ t , 1 ) ⊗ α z ( ϕ t , 1 ) ) H Канал BS-БПЛА h 2 \mathbf{h}_2 h 2 , канал RIS-БПЛА h 3 \mathbf{h}_3 h 3 , канал БПЛА-UE h 4 \mathbf{h}_4 h 4 и другиеБазовая станция использует многолучевой метод:
Фиксированный луч f 0 \mathbf{f}_0 f 0 , направленный на RIS Изменяющийся во времени луч f k \mathbf{f}_k f k , направленный в область неба Изменяющийся во времени луч RIS ω k \boldsymbol{\omega}_k ω k , охватывающий целевую область Модель принятого сигнала:
y k = ζ h 4 ( h 3 diag ( ω k ) H 1 F k s k + h 2 F k s k ) + H 5 F k s k + n k \mathbf{y}_k = \zeta\mathbf{h}_4(\mathbf{h}_3\text{diag}(\boldsymbol{\omega}_k)\mathbf{H}_1\mathbf{F}_k\mathbf{s}_k + \mathbf{h}_2\mathbf{F}_k\mathbf{s}_k) + \mathbf{H}_5\mathbf{F}_k\mathbf{s}_k + \mathbf{n}_k y k = ζ h 4 ( h 3 diag ( ω k ) H 1 F k s k + h 2 F k s k ) + H 5 F k s k + n k
где ζ \zeta ζ — эффективная площадь рассеяния (RCS) БПЛА.
Используется метод координатного градиентного спуска (CGD):
Целевая функция:
min f ( h 2 , h 3 , h 4 ) = ∥ Y − P / 2 H ~ Ω ˉ − P / 2 H ^ F ˉ ∥ F 2 \min f(\mathbf{h}_2,\mathbf{h}_3,\mathbf{h}_4) = \left\|\mathbf{Y} - \sqrt{P/2}\tilde{\mathbf{H}}\bar{\boldsymbol{\Omega}} - \sqrt{P/2}\hat{\mathbf{H}}\bar{\mathbf{F}}\right\|_F^2 min f ( h 2 , h 3 , h 4 ) = Y − P /2 H ~ Ω ˉ − P /2 H ^ F ˉ F 2
Итеративное обновление:
h 4 ( k ) = h 4 ( k − 1 ) − η ∂ f ∂ h 4 ∣ h 4 = h 4 ( k − 1 ) \mathbf{h}_4^{(k)} = \mathbf{h}_4^{(k-1)} - \eta \frac{\partial f}{\partial \mathbf{h}_4}\bigg|_{\mathbf{h}_4=\mathbf{h}_4^{(k-1)}} h 4 ( k ) = h 4 ( k − 1 ) − η ∂ h 4 ∂ f h 4 = h 4 ( k − 1 )
Извлечение угловых параметров посредством двумерного поиска:
θ ^ t , 2 , ϕ ^ t , 2 = arg max θ t , 2 , ϕ t , 2 ∣ h ^ 2 ( α y ( θ t , 2 , ϕ t , 2 ) ⊗ α z ( ϕ t , 2 ) ) ∣ \hat{\theta}_{t,2}, \hat{\phi}_{t,2} = \arg\max_{\theta_{t,2},\phi_{t,2}} \left|\hat{\mathbf{h}}_2(\boldsymbol{\alpha}_y(\theta_{t,2}, \phi_{t,2}) \otimes \boldsymbol{\alpha}_z(\phi_{t,2}))\right| θ ^ t , 2 , ϕ ^ t , 2 = arg max θ t , 2 , ϕ t , 2 h ^ 2 ( α y ( θ t , 2 , ϕ t , 2 ) ⊗ α z ( ϕ t , 2 ))
Отображение оценок углов в трёхмерное положение на основе принципа наименьших квадратов:
p ^ D = ( B B + B R + B U ) − 1 ( B B p B + B R p R + B U p U ) \hat{\mathbf{p}}_D = (\mathbf{B}_B + \mathbf{B}_R + \mathbf{B}_U)^{-1}(\mathbf{B}_B\mathbf{p}_B + \mathbf{B}_R\mathbf{p}_R + \mathbf{B}_U\mathbf{p}_U) p ^ D = ( B B + B R + B U ) − 1 ( B B p B + B R p R + B U p U )
Использование многолучевости : одновременное использование двух путей: BS-БПЛА-UE и BS-RIS-БПЛА-UEПроектирование лучей : обеспечение ортогональности лучей базовой станции, снижение помех от прямого путиАлгоритм низкой сложности : метод CGD имеет более низкую вычислительную сложность по сравнению с минимизацией ядерной нормыТеоретические границы : предоставление анализа информации Фишера и CRLB в качестве эталона производительностиКонфигурация антенн : M B = 8 × 8 M_B = 8 \times 8 M B = 8 × 8 , M U = 4 × 4 M_U = 4 \times 4 M U = 4 × 4 , M R = 6 × 6 M_R = 6 \times 6 M R = 6 × 6 Положения узлов :
BS: (0, 0, 26) RIS: (0, 0.5, 25.5) БПЛА: (3, 3, 30) UE: (2, 2, 24) Параметры канала : показатель потерь пути Γ = 2 \Gamma = 2 Γ = 2 , полоса пропускания B = 20 B = 20 B = 20 МГцДисперсия шума : σ 2 = − 101 \sigma^2 = -101 σ 2 = − 101 дБмВ качестве основной метрики оценки используется среднеквадратическая ошибка (RMSE):
RMSE = E { ( p D − p ^ D ) T ( p D − p ^ D ) } \text{RMSE} = \sqrt{E\{(\mathbf{p}_D - \hat{\mathbf{p}}_D)^T(\mathbf{p}_D - \hat{\mathbf{p}}_D)\}} RMSE = E {( p D − p ^ D ) T ( p D − p ^ D )}
С RIS против без RIS : проверка выигрыша в производительности RISТеоретический CRLB : использование в качестве верхней границы производительностиРезультаты моделирования показывают, что в диапазоне SNR от -10 дБ до 10 дБ:
Система с RIS : демонстрирует значительные преимущества в области низкого SNRРазница в производительности : при SNR -10 дБ система RIS снижает RMSE примерно на один порядок по сравнению с системой без RISПриближение к теории : предложенный алгоритм показывает производительность, близкую к теоретической границе CRLBВлияние различного количества временных слотов обучения K:
Большее значение K приводит к более высокой точности локализации Существует точка насыщения производительности, где предельная выгода от увеличения K снижается Значительное влияние эффективной площади рассеяния на производительность:
ζ = 0.5 \zeta = 0.5 ζ = 0.5 : относительно худшая производительностьζ = 1.0 \zeta = 1.0 ζ = 1.0 : эталонная производительностьζ = 2.0 \zeta = 2.0 ζ = 2.0 : значительное улучшение производительностиПреимущество в низком SNR : выигрыш в производительности RIS наиболее значителен при низких отношениях сигнал-шумЧувствительность к RCS : эффективная площадь рассеяния БПЛА является ключевым фактором, влияющим на точность локализацииЭффективность алгоритма : алгоритм CGD достигает производительности, близкой к теоретическому пределу, при низкой вычислительной сложностиТехнологии локализации 5G/6G : методы локализации на основе mmWave и крупномасштабного MIMOИнтеграция радиолокации и восприятия : обнаружение и локализация целей в системах ISACRIS-ассистированная связь : применение переконфигурируемых интеллектуальных поверхностей в беспроводной связиПо сравнению с традиционной радиолокацией: использование существующей сотовой инфраструктуры, снижение затрат на развёртывание По сравнению с активной локализацией: применимость к некооперативным целям, расширение области применения По сравнению с существующими исследованиями RIS: первое применение RIS к сценарию пассивной локализации БПЛА RIS значительно повышает производительность пассивной локализации БПЛА в системах mmWave Предложенный алгоритм CGD достигает отличной производительности при низкой вычислительной сложности Параметры системы (RCS, затраты на обучение) оказывают значительное влияние на производительность Предположения о канале : рассматриваются только пути прямой видимости (LoS), в реальной среде могут быть многолучевые отраженияСтатические сценарии : не учитывается влияние движения БПЛА на точность локализацииИдеальная синхронизация : предполагается идеальная синхронизация по времени и частоте в системеИзвестные положения : предполагается, что положения RIS и UE точно известныОптимизация распределения мощности : исследование оптимальной стратегии распределения мощности между многолучевыми передачами базовой станцииДинамическое отслеживание : расширение на отслеживание траектории движущихся БПЛАОбработка многолучевости : рассмотрение надёжной локализации в условиях распространения без прямой видимостиПрактическое развёртывание : проверка производительности алгоритма в реальных условияхВысокая инновационность : первое применение технологии RIS к пассивной локализации БПЛА, обладает высокой инновационностьюПолнота теории : предоставляет полный анализ информации Фишера и вывод CRLBПрактичность алгоритма : алгоритм CGD имеет низкую сложность, удобен для практической реализацииКомплексный анализ : систематический анализ влияния ключевых параметров на производительностьДостаточная верификация : теоретический анализ и результаты моделирования взаимно подтверждают друг другаИдеализированные предположения : предположение о канале LoS может не соответствовать действительности в реальной средеОграниченные сценарии : в основном применимо к локализации высотных БПЛА в открытых пространствахОграничения экспериментов : отсутствие практических испытаний, только результаты моделированияНедостаточное подавление помех : недостаточное рассмотрение подавления помех в сложной электромагнитной средеАкадемическая ценность : предоставляет новые идеи для исследований применения RIS в приложениях восприятияПрактические перспективы : имеет потенциал применения в умных городах, системах безопасности и мониторингаТехнологический прогресс : способствует развитию и стандартизации технологии ISACСистемы мониторинга БПЛА на крышах городских зданий Сети обнаружения БПЛА вокруг аэропортов Защита воздушного пространства важных объектов Приложения с добавленной стоимостью в сетях 5G/6G Статья ссылается на множество важных работ в смежных областях, включая технологии локализации 5G, системы ISAC и коммуникацию RIS, что обеспечивает прочную теоретическую базу для данной работы.