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 Articolo : 2404.14879Titolo : Device-Free 3D Drone Localization in RIS-Assisted mmWave MIMO NetworksAutori : Jiguang He, Charles Vanwynsberghe, Hui Chen, Chongwen Huang, Aymen FakhreddineClassificazione : eess.SP (Elaborazione dei Segnali)Data di Pubblicazione : Aprile 2024 (preprint arXiv, aggiornato il 2 gennaio 2025)Link Articolo : https://arxiv.org/abs/2404.14879 Questo articolo esamina il potenziale delle Superfici Intelligenti Riconfigurabili (RIS) per realizzare la localizzazione passiva/senza dispositivo tridimensionale (3D) di droni nell'infrastruttura cellulare multi-antenna a banda millimetrica (mmWave). Il sistema di localizzazione sviluppato opera in modalità bistatica, senza richiedere comunicazione diretta tra il drone e la stazione base. I limiti teorici di prestazione sono analizzati attraverso l'analisi dell'informazione di Fisher e il limite inferiore di Cramér Rao (CRLB). Inoltre, è stato sviluppato un algoritmo di localizzazione dei droni a bassa complessità ma efficace basato sulla discesa del gradiente coordinato, e sono stati investigati i fattori che influenzano le prestazioni del sistema, come la sezione trasversale radar (RCS) del drone e il carico di addestramento. Sia l'analisi teorica che le simulazioni numeriche dimostrano vantaggi significativi del sistema integrato con RIS rispetto ai sistemi senza RIS.
Il problema fondamentale affrontato in questo articolo è: realizzare la localizzazione passiva tridimensionale di droni in volo nell'infrastruttura cellulare esistente, ovvero determinare accuratamente la posizione spaziale del drone senza che esso comunichi attivamente con la stazione base.
Esigenze di Sicurezza : Con la diffusione delle applicazioni dei droni, la necessità di rilevamento e localizzazione di droni non autorizzati o malevoli diventa sempre più urgenteRiutilizzo dell'Infrastruttura : Sfruttare pienamente l'infrastruttura di rete cellulare 5G/6G esistente, evitando la distribuzione aggiuntiva di sistemi radar dedicatiFusione Tecnologica : Scenario di applicazione importante del paradigma di Integrazione di Comunicazione e Sensing (ISAC)Sistemi Radar Tradizionali : Richiedono apparecchiature radar specializzate, con costi elevati e copertura limitataSistemi di Localizzazione Attiva : Richiedono comunicazione tra il bersaglio e la stazione base, non applicabili a bersagli non cooperativiSistemi Passivi Esistenti : Prestazioni limitate nella banda mmWave, precisione di localizzazione insufficienteUtilizzare RIS per potenziare i percorsi di propagazione del segnale e migliorare le prestazioni di localizzazione passiva Realizzare localizzazione tridimensionale ad alta precisione nei sistemi MIMO mmWave Sviluppare algoritmi a bassa complessità pratici, prossimi ai limiti teorici di prestazione Progettazione del Sistema : Propone un'architettura di sistema di localizzazione passiva di droni mmWave MIMO assistita da RISAnalisi Teorica : Analizza i limiti teorici di prestazione del sistema attraverso la matrice di informazione di Fisher e CRLBSviluppo dell'Algoritmo : Sviluppa un algoritmo di localizzazione tridimensionale a tre fasi a bassa complessità basato sulla discesa del gradiente coordinato (CGD)Valutazione delle Prestazioni : Analizza in modo completo l'impatto dei parametri chiave come RCS e carico di addestramento sulle prestazioni del sistemaVerifica dei Vantaggi : Dimostra attraverso teoria e simulazione i vantaggi significativi del sistema assistito da RIS rispetto ai metodi tradizionaliInput : Segnali ricevuti multi-slot presso il dispositivo utente ricevente (UE)
Output : Coordinate tridimensionali della posizione del drone ( x D , y D , z D ) (x_D, y_D, z_D) ( x D , y D , z D ) Vincoli : Nessuna comunicazione diretta tra il drone e la stazione base, il sistema opera in modalità bistatica
Il sistema è costituito da tre componenti chiave:
Stazione Base mmWave (BS) : Equipaggiata con un array planare uniforme di 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 antenneSuperficie Intelligente Riconfigurabile (RIS) : Contiene 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 elementi riflettentiDispositivo Utente Ricevente (UE) : Equipaggiato con M U M_U M U antenneAdotta il modello di canale Saleh-Valenzuela, considerando le caratteristiche di scattering sparse nella banda mmWave. I canali chiave includono:
Canale 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 Canale BS-Drone h 2 \mathbf{h}_2 h 2 , Canale RIS-Drone h 3 \mathbf{h}_3 h 3 , Canale Drone-UE h 4 \mathbf{h}_4 h 4 , e altriLa stazione base adotta un metodo multi-fascio:
Fascio fisso f 0 \mathbf{f}_0 f 0 orientato verso RIS Fascio variabile nel tempo f k \mathbf{f}_k f k orientato verso la regione aerea Fascio variabile nel tempo RIS ω k \boldsymbol{\omega}_k ω k che copre l'area bersaglio Modello del segnale ricevuto:
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
dove ζ \zeta ζ è la sezione trasversale radar (RCS) del drone.
Adotta il metodo della discesa del gradiente coordinato (CGD) per stimare il canale in cascata:
Funzione obiettivo:
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
Aggiornamento iterativo:
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 )
Estrae i parametri angolari attraverso ricerca bidimensionale:
θ ^ 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 ))
Mappa le stime angolari alla posizione tridimensionale basandosi sul principio dei minimi quadrati:
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 )
Sfruttamento dei Percorsi Multipli : Utilizza contemporaneamente i due percorsi BS-Drone-UE e BS-RIS-Drone-UEProgettazione dei Fasci : Assicura l'ortogonalità dei fasci della stazione base, riducendo l'interferenza del percorso direttoAlgoritmo a Bassa Complessità : Il metodo CGD ha complessità computazionale inferiore rispetto alla minimizzazione della norma nucleareLimiti Teorici : Fornisce analisi dell'informazione di Fisher e CRLB come benchmark di prestazioneConfigurazione Antenne : 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 Posizioni dei Nodi :
BS: (0, 0, 26) RIS: (0, 0.5, 25.5) Drone: (3, 3, 30) UE: (2, 2, 24) Parametri del Canale : Indice di perdita di percorso Γ = 2 \Gamma = 2 Γ = 2 , larghezza di banda B = 20 B = 20 B = 20 MHzVarianza del Rumore : σ 2 = − 101 \sigma^2 = -101 σ 2 = − 101 dBmUtilizza l'errore quadratico medio (RMSE) come metrica di valutazione principale:
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 )}
Con RIS vs Senza RIS : Verifica il guadagno di prestazione di RISCRLB Teorico : Serve come limite superiore di prestazioneI risultati della simulazione mostrano, nell'intervallo SNR da -10dB a 10dB:
Sistema con RIS : Mostra vantaggi significativi nella regione a basso SNRDifferenza di Prestazione : A SNR di -10dB, il sistema RIS riduce l'RMSE di circa un ordine di grandezza rispetto al sistema senza RISProssimità Teorica : L'algoritmo proposto ha prestazioni prossime al limite teorico CRLBEffetto di diversi numeri di slot di addestramento K:
Maggiore è il valore di K, maggiore è la precisione di localizzazione Esiste un punto di saturazione delle prestazioni, con rendimenti marginali decrescenti dall'aumento eccessivo di K Impatto significativo della sezione trasversale radar sulle prestazioni:
ζ = 0.5 \zeta = 0.5 ζ = 0.5 : Prestazioni relativamente scarseζ = 1.0 \zeta = 1.0 ζ = 1.0 : Benchmark di prestazione standardζ = 2.0 \zeta = 2.0 ζ = 2.0 : Miglioramento significativo delle prestazioniVantaggio a Basso SNR : Il guadagno di prestazione di RIS è più significativo in condizioni di basso rapporto segnale-rumoreSensibilità RCS : La sezione trasversale radar del drone è un fattore chiave che influenza la precisione di localizzazioneEfficacia dell'Algoritmo : L'algoritmo CGD realizza prestazioni prossime al limite teorico con bassa complessitàTecnologie di Localizzazione 5G/6G : Metodi di localizzazione basati su mmWave e MIMO su larga scalaFusione Radar-Sensing : Rilevamento e localizzazione di bersagli nei sistemi ISACComunicazione Assistita da RIS : Applicazioni di Superfici Intelligenti Riconfigurabili nelle comunicazioni wirelessRispetto ai radar tradizionali: Utilizza l'infrastruttura cellulare esistente, riducendo i costi di distribuzione Rispetto alla localizzazione attiva: Applicabile a bersagli non cooperativi, ampliando la gamma di applicazioni Rispetto alla ricerca RIS esistente: Prima applicazione di RIS allo scenario di localizzazione passiva di droni RIS migliora significativamente le prestazioni di localizzazione passiva di droni nei sistemi mmWave L'algoritmo CGD proposto realizza prestazioni eccellenti con bassa complessità I parametri del sistema (RCS, carico di addestramento) hanno un impatto importante sulle prestazioni Ipotesi di Canale : Considera solo percorsi Line-of-Sight (LoS), l'ambiente reale potrebbe avere multipercorsiScenario Statico : Non considera l'impatto del movimento del drone sulle prestazioni di localizzazioneSincronizzazione Perfetta : Assume sincronizzazione temporale e di frequenza perfetta del sistemaPosizioni Conosciute : Assume che le posizioni di RIS e UE siano note con precisioneOttimizzazione dell'Allocazione di Potenza : Ricerca della strategia di allocazione di potenza ottimale tra i multi-fasci della stazione baseTracciamento Dinamico : Estensione al tracciamento della traiettoria di droni in movimentoElaborazione dei Multipercorsi : Localizzazione robusta in ambienti di propagazione non Line-of-SightDistribuzione Pratica : Verifica delle prestazioni dell'algoritmo in ambienti realiForte Innovatività : Prima applicazione della tecnologia RIS alla localizzazione passiva di droni, con forte innovativitàCompletezza Teorica : Fornisce analisi completa dell'informazione di Fisher e derivazione di CRLBPraticità dell'Algoritmo : L'algoritmo CGD ha bassa complessità, facilitando l'implementazione praticaAnalisi Completa : Analizza sistematicamente l'impatto dei parametri chiave sulle prestazioniVerifica Sufficiente : L'analisi teorica e i risultati di simulazione si verificano reciprocamenteIpotesi Idealizzate : L'ipotesi di canale LoS potrebbe non essere valida negli ambienti realiScenari Limitati : Principalmente applicabile alla localizzazione di droni ad alta quota in spazi apertiLimitazioni Sperimentali : Mancanza di verifica pratica, solo risultati di simulazioneElaborazione dell'Interferenza : Considerazione insufficiente della soppressione dell'interferenza in ambienti elettromagnetici complessiValore Accademico : Fornisce nuove prospettive per la ricerca su RIS nelle applicazioni di sensingProspettive Pratiche : Ha potenziale di applicazione in città intelligenti, sorveglianza e sicurezzaPromozione Tecnologica : Promuove lo sviluppo e la standardizzazione della tecnologia ISACSistemi di monitoraggio di droni sui tetti di edifici urbani Reti di rilevamento di droni intorno agli aeroporti Protezione dello spazio aereo di strutture critiche Applicazioni di servizi a valore aggiunto per reti 5G/6G L'articolo cita numerosi lavori importanti in campi correlati, inclusa la ricerca all'avanguardia in tecnologie di localizzazione 5G, sistemi ISAC e comunicazioni assistite da RIS, fornendo una base teorica solida per questo lavoro.