3D Gaussian Distribution and RGB Reconstruction (3D Gaussian-RGB) based GNSS Multipath Detection — Part I: Derivation and Analysis of the 3D Gaussian ellipsoid Modelling

Sintesi Tecnica: Distribuzione Gaussiana 3D e Ricostruzione RGB (3D Gaussian-RGB) per la Rilevazione del Multipath GNSS — Parte I: Derivazione e Analisi della Modellazione dell'ellissoide Gaussiano 3D

Definizione del Problema
I sistemi di navigazione satellitare globale (GNSS) affrontano sfide significative in termini di accuratezza e robustezza del posizionamento in ambienti complessi, come i canyon urbani, a causa degli effetti multipath. Sebbene i metodi di rilevazione esistenti utilizzino caratteristiche nel dominio della misura (es. pseudorange, C/N0) o distorsioni nel dominio della correlazione, essi presentano spesso limitazioni. Gli approcci nel dominio della misura possono oscurare le distorsioni del segnale originale a causa dei processi di filtraggio e stima del ricevitore. Al contrario, i metodi esistenti nel dominio della correlazione si affidano frequentemente a correlatori ad alta risoluzione o non riescono a esplorare sufficientemente le relazioni non lineari e le connessioni intrinseche tra le caratteristiche del multipath e i parametri di correlazione. È necessaria una base teorica che bilanci la capacità di rilevazione con la complessità computazionale, tenendo conto delle dipendenze statistiche tra le uscite dei correlatori Early (E), Prompt (P) e Late (L).

Metodologia
Questo studio propone un modello teorico nel dominio della correlazione basato sull'analisi statistica dei valori di correlazione E, P e L dal ramo in quadratura (Q) dei segnali GNSS. Il ramo Q è stato selezionato perché è più sensibile alle distorsioni di codice e fase indotte dal multipath rispetto al ramo in fase (I), specialmente quando il loop di tracking è agganciato.

La metodologia procede attraverso i seguenti passaggi:

1. Modellazione del Segnale: Vengono stabiliti modelli matematici per i segnali LOS (Line-of-Sight), NLOS (Non-Line-of-Sight) e multipath. Il segnale multipath è modellato come una sintesi di componenti LOS e NLOS, introducendo parametri per il ritardo relativo ($\delta\tau_N$), l'attenuazione ($\alpha$), lo spostamento Doppler e la rotazione di fase ($\delta\phi_N$).
2. Normalizzazione e Linearizzazione: I valori di correlazione E/P/L nel ramo Q vengono normalizzati rispetto all'uscita del correlatore Prompt per mitigare le inconsistenze della forza del segnale variabile nel tempo. Viene applicata un'espansione di Taylor multivariata del primo ordine per linearizzare questi valori attorno al vettore del rumore, assumendo che il rumore sia piccolo rispetto alla forza del segnale.
3. Distribuzione Gaussiana 3D Congiunta: Si dimostra che i valori E/P/L linearizzati seguono una distribuzione Gaussiana congiunta 3D. Il vettore media ($\mu_R$) e la matrice di covarianza ($\Sigma_R$) sono derivati, tenendo conto della coerenza del rumore (tramite decomposizione di Cholesky) e del jitter di fase propagato attraverso il loop di tracking.
4. Modellazione dell'Ellissoide di Equiprobabilità: La distribuzione congiunta è rappresentata geometricamente come un ellissoide di equiprobabilità 3D. La decomposizione spettrale della matrice di covarianza fornisce una base ortogonale per l'ellissoide. Il modello caratterizza le condizioni del segnale attraverso tre trasformazioni geometriche:
   • Traslazione: Governata dal vettore media $\mu_R$, che influenza principalmente i segnali multipath a causa delle differenze di fase e ritardo.
   • Scaling (Ridimensionamento): Determinato dagli autovalori della matrice di covarianza, influenzato dalla varianza del rumore e dall'ampiezza del segnale (C/N0).
   • Rotazione: Definita dagli autovettori della matrice di covarianza, dipendente dalla spaziatura ($\Delta$) dei correlatori Early-Late.
5. Derivazione della Soglia di Rilevazione: Utilizzando il criterio di Neyman-Pearson (NP), viene derivata una soglia di rilevazione. Il problema del test d'ipotesi è formulato come la distinzione tra LOS (Ipotesi Nulla $H_0$) e multipath (Ipotesi Alternativa $H_1$). La soglia di rilevazione $\gamma$ corrisponde al confine dell'ellissoide di equiprobabilità per una specificata probabilità di falso allarme ($P_{fa}$).

Contributi Principali
Il documento presenta quattro contributi primari:

1. Derivazione della Distribuzione Gaussiana Congiunta 3D: Gli autori derivano la distribuzione Gaussiana congiunta 3D dei valori E-P-L, tenendo esplicitamente conto dell'interdipendenza tra questi valori all'interno della stessa funzione di autocorrelazione (ACF).
2. Modello dell'Ellissoide di Equiprobabilità: Viene stabilito un modello ellissoidale in cui le proprietà geometriche (traslazione, scaling, rotazione) sono analiticamente legate ai parametri del multipath. Ciò consente di quantificare come i parametri del multipath (ritardo, attenuazione, fase) deformino la distribuzione statistica.
3. Soglia di Rivelazione Teorica: Una soglia di rilevazione viene derivata utilizzando il criterio NP, utilizzando l'ellissoide di equiprobabilità come confine decisionale. Ciò fornisce una base teorica per quantificare la probabilità di rilevazione ($P_d$) sotto diverse condizioni di multipath.
4. Validazione tramite Simulazione: I modelli teorici sono validati attraverso simulazioni Monte Carlo, dimostrando l'efficacia dei modelli Gaussiani 3D per la rilevazione del multipath.

Risultati
I risultati delle simulazioni utilizzando un modello semi-analitico confermano le derivazioni teoriche:

• Gaussianità: Per intensità di segnale (C/N0) superiori a 30 dB-Hz, i valori di correlazione E-P-L simulati mostrano una forte gaussianità, allineandosi strettamente con la distribuzione congiunta teorica. Al di sotto di questa soglia, si verificano deviazioni, in particolare negli scenari multipath, a causa dell'influenza non trascurabile dei termini di rumore di ordine superiore.
• Consistenza Geometrica: Gli ellissoidi di equiprobabilità simulati si allineano bene con i modelli teorici per C/N0 > 30 dB-Hz. Nello specifico, i segnali multipath mostrano un effetto di traslazione distinto lungo gli assi X e Z, mentre i segnali LOS e NLOS mostrano una traslazione trascurabile. Le proprietà di scaling e rotazione rimangono coerenti tra i diversi tipi di segnale quando il C/N0 è fisso, poiché sono determinate principalmente dalla spaziatura dei correlatori.
• Prestazioni di Rilevazione: L'analisi della curva ROC (Receiver Operating Characteristic) mostra che le prestazioni di rilevazione migliorano con C/N0 più elevati e minore attenuazione dell'energia ($\alpha$). Le prestazioni sono non monotone rispetto al ritardo relativo del codice ($\delta\tau_N$) e alla rotazione di fase ($\delta\phi_N$), con picchi a valori intermedi (circa 0.6 chip e 0.6$\pi$, rispettivamente). Il metodo presenta punti ciechi di rilevazione per segnali multipath con ritardi o fasi estremamente brevi/lunghi o vicini a 0 o $\pi$.
• Accuratezza della Soglia: La soglia di rilevazione teorica $\gamma$ coincide strettamente con i risultati simulati quando C/N0 > 30 dB-Hz. Le discrepanze aumentano a intensità di segnale inferiori a causa della degradazione della gaussianità della distribuzione.

Significatività e Rivendicazioni
Il documento afferma che questo lavoro fornisce la base statistica per un quadro più ampio di rilevazione del multipath. Stabilendo un rigoroso modello dell'ellissoide Gaussiano 3D, lo studio offre un metodo geometricamente interpretabile per caratterizzare le distorsioni del segnale indotte dal multipath. Gli autori dichiarano che questo modello statistico funge da prerequisito per la Parte II della loro ricerca, dove il vettore E/P/L normalizzato viene proiettato nello spazio colore RGB. Il presente lavoro (Parte I) si concentra sulla derivazione delle soglie teoriche e sull'analisi delle proprietà statistiche, mentre il saggio complementare utilizzerà queste soglie per interpretare le caratteristiche del colore.

Lo studio sottolinea che, sebbene il metodo distingua efficacemente tra segnali LOS e multipath basandosi sulla traslazione della distribuzione, distinguere tra segnali LOS e NLOS rimane complesso sotto simili condizioni di C/N0 a causa della sovrapposizione delle caratteristiche di scaling degli ellissoidi. Gli autori mantengono un ambito modesto, osservando che l'assunzione di una distribuzione Gaussiana a media zero per il correlatore Prompt negli scenari multipath è un'approssimazione giustificata dalla validazione empirica entro margini di errore pratici, nonostante le potenziali perturbazioni della sincronizzazione di fase negli ambienti reali.