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

Resumen Técnico: Distribución Gaussiana 3D y Reconstrucción RGB (3D Gaussian-RGB) basada en la detección de multipath de GNSS — Parte I: Derivación y Análisis del Modelado del Elipsoide Gaussiano 3D

Planteamiento del Problema
Los Sistemas Globales de Navegación por Satélite (GNSS) enfrentan desafíos significativos en la precisión y robustez del posicionamiento en entornos complejos, como cañones urbanos, debido a los efectos de multipath (multitrayectoria). Aunque los métodos de detección existentes utilizan características en el dominio de la medición (por ejemplo, pseudodistancia, C/N0) o distorsiones en el dominio de la correlación, estos suelen presentar limitaciones. Los enfoques en el dominio de la medición pueden ocultar las distorsiones originales de la señal debido al filtrado del receptor y a los procesos de estimación. Por el contrario, los métodos existentes en el dominio de la correlación dependen frecuentemente de correladores de alta resolución o no logran explorar suficientemente las relaciones no lineales y las conexiones intrínsecas entre las características de multipath y los parámetros de correlación. Existe la necesidad de un modelo teórico que equilibre la capacidad de detección con la complejidad computacional, considerando al mismo tiempo las dependencias estadísticas entre las salidas de los correladores Early (E), Prompt (P) y Late (L).

Metodología
Este estudio propone un modelo teórico de correlación basado en el análisis estadístico de los valores de correlación E, P y L de la rama de cuadratura (Q) de las señales GNSS. Se selecciona la rama Q porque es más sensible a las distorsiones de código y fase inducidas por el multipath que la rama de fase (I), particularmente cuando el bucle de seguimiento está enganchado.

La metodología procede a través de los siguientes pasos:

1. Modelado de la Señal: Se establecen modelos matemáticos para las señales LOS (Línea de Vista), NLOS (Sin Línea de Vista) y multipath. La señal de multipath se modela como una síntesis de los componentes LOS y NLOS, introduciendo parámetros para el retardo relativo ($\delta\tau_N$), desvanecimiento ($\alpha$), desplazamiento Doppler y rotación de fase ($\delta\phi_N$).
2. Normalización y Linealización: Los valores de correlación E/P/L en la rama Q se normalizan respecto a la salida del correlador Prompt para mitigar las inconsistencias de la fuerza de la señal variables en el tiempo. Se aplica una expansión de Taylor multivariante de primer orden para linealizar estos valores alrededor del vector de ruido, asumiendo que el ruido es pequeño en relación con la fuerza de la señal.
3. Distribución Gaussiana Conjunta 3D: Se demuestra que los valores E/P/L linealizados siguen una distribución Gaussiana conjunta 3D. Se derivan el vector de media ($\mu_R$) y la matriz de covarianza ($\Sigma_R$), contabilizando la coherencia del ruido (vía descomposición de Cholesky) y el jitter de fase propagado a través del bucle de seguimiento.
4. Modelado del Elipsoide de Equiprobabilidad: La distribución conjunta se representa geométricamente como un elipsoide de equiprobabilidad 3D. La descomposición espectral de la matriz de covarianza proporciona una base ortogonal para el elipsoide. El modelo caracteriza las condiciones de la señal a través de tres transformaciones geométricas:
   • Traslación: Gobernada por el vector de media $\mu_R$, afectando principalmente a las señales de multipath debido a las diferencias de fase y retardo.
   • Escalamiento: Determinado por los autovalores de la matriz de covarianza, influenciado por la varianza del ruido y la amplitud de la señal (C/N0).
   • Rotación: Definida por los autovectores de la matriz de covarianza, dependiente del espaciamiento de los correladores Early-Late ($\Delta$).
5. Derivación del Umbral de Detección: Utilizando el criterio de Neyman-Pearson (NP), se deriva un límite de detección. El problema de prueba de hipótesis se enmarca en la distinción entre LOS (Hipótesis Nula $H_0$) y multipath (Hipótesis Alternativa $H_1$). El umbral de detección $\gamma$ corresponde al límite del elipsoide de equiprobabilidad para una probabilidad de falsa alarma ($P_{fa}$) especificada.

Contribuciones Clave
El artículo realiza cuatro contribuciones principales:

1. Derivación de la Distribución Gaussiana Conjunta 3D: Los autores derivan la distribución Gaussiana conjunta 3D de los valores E-P-L, contabilizando explícitamente la interdependencia entre estos valores dentro de la misma función de autocorrelación (ACF).
2. Modelo de Elipsoide de Equiprobabilidad: Se establece un modelo de elipsoide donde las propiedades geométricas (traslación, escalamiento, rotación) están vinculadas analíticamente a los parámetros de multipath. Esto permite cuantificar cómo los parámetros de multipath (retardo, desvanecimiento, fase) deforman la distribución estadística.
3. Umbral de Detección Teórico: Se deriva un umbral de detección utilizando el criterio NP, utilizando el elipsoide de equiprobidad como frontera de decisión. Esto proporciona una base teórica para cuantificar la probabilidad de detección ($P_d$) bajo diversas condiciones de multipath.
4. Validación mediante Simulación: Los modelos teóricos se validan mediante simulaciones de Monte Carlo, demostrando la efectividad de los modelos Gaussianos 3D para la detección de multipath.

Resultados
Los resultados de la simulación utilizando un modelo semi-analítico confirman las derivaciones teóricas:

• Gaussianidad: Para fuerzas de señal (C/N0) superiores a 30 dB-Hz, los valores de correlación E-P-L simulados exhiben una fuerte gaussianidad, alineándose estrechamente con la distribución conjunta teórica. Por debajo de este umbral, ocurren desviaciones, particularmente en escenarios de multipath, debido a la influencia no despreciable de términos de ruido de orden superior.
• Consistencia Geométrica: Los elipsoides de equiprobabilidad simulados se alinean bien con los modelos teóricos cuando C/N0 > 30 dB-Hz. Específicamente, las señales de multipath exhiben un efecto de traslación distintivo a lo largo de los ejes X y Z, mientras que las señales LOS y NLOS muestran una traslación insignificante. Las propiedades de escalamiento y rotación permanecen consistentes entre los tipos de señal cuando el C/N0 es fijo, ya que están determinadas primordialmente por el espaciamiento del correlador.
• Desempeño de Detección: El análisis de la curva ROC muestra que el desempeño de la detección mejora con un mayor C/N0 y un menor desvanecimiento de energía ($\alpha$). El desempeño es no monotónico con respecto al retardo de código relativo ($\delta\tau_N$) y la rotación de fase ($\delta\phi_N$), alcanzando su punto máximo en valores intermedios (aprox. 0.6 chips y 0.6$\pi$, respectivamente). El método presenta puntos ciegos de detección para señales de multipath con retardos o fases extremadamente cortos/largos o fases cercanas a 0 o $\pi$.
• Precisión del Umbral: El umbral de detección teórico $\gamma$ coincide estrechamente con los resultados simulados cuando C/N0 > 30 dB-Hz. Las discrepancias aumentan en niveles de fuerza de señal más bajos debido a la degradación de la gaussianidad de la distribución.

Significancia y Reivindicaciones
El artículo afirma que este trabajo proporciona la base estadística para un marco de detección de multipath más amplio. Al establecer un modelo riguroso de elipsoide Gaussiano 3D, el estudio ofrece un método geométricamente interpretable para caracterizar las distorsiones de la señal inducidas por el multipath. Los autores declaran que este modelo estadístico sirve como prerrequisito para la Parte II de su investigación, donde el vector E/P/L normalizado se proyecta en el espacio de color RGB. El presente trabajo (Parte I) se centra en la derivación de los umbrales teóricos y el análisis de las propiedades estadísticas, mientras que el artículo complementario utilizará estos umbrales para interpretar las características de color.

El estudio enfatiza que, si bien el método distingue eficazmente las señales LOS de las de multipath basándose en la traslación de la distribución, distinguir entre señales LOS y NLOS sigue siendo un desafío bajo condiciones de C/N0 similares debido al solapamiento de las características de escalamiento de los elipsoides. Los autores mantienen un alcance modesto, señalando que la suposición de Gauss de media cero para el correlador Prompt en escenarios de multipath es una aproximación justificada por la validación empírica dentro de márgenes de error prácticos, a pesar de las posibles perturbaciones de sincronización de fase en entornos del mundo real.