A Directivity-Dependent Rician K-Factor Model for Indoor Industrial Channels

Resumen Técnico: Un Modelo de Factor K de Rician Dependiente de la Directividad para Canales Industriales Interiores

Planteamiento del Problema
El despliegue de enlaces de banda E (71–86 GHz) en grandes entornos industriales parcialmente reverberantes (por ejemplo, almacenes, salas de ensamblaje) enfrenta desafíos significativos debido a las altas dispersiones de retardo. Estos entornos actúan como cavidades reverberantes, produciendo componentes de multitrayectoria que degradan severamente el rendimiento de banda ancha. Si bien estudios empíricos en diversas frecuencias (de 28 GHz a 304 GHz) han confirmado que las antenas directivas suprimen las componentes de multitrayectoria de llegada tardía y reducen la dispersión de retardo de la raíz cuadrática media (RMS) ($\sigma_\tau$), la literatura existente carece de un modelo basado en física y de forma cerrada. Los enfoques actuales dependen de ajustes empíricos que no pueden predecir $\sigma_\tau$ ni el retardo excesivo medio ($\bar{\tau}$) a partir de parámetros medibles de la sala o de la ganancia de la antena. Existe una necesidad específica de un modelo unificado que mapee la ganancia total del sistema ($G$) hacia las estadísticas de retardo utilizando un único escalar físicamente interpretable para cuantificar la desviación respecto a una cavidad isotrópica ideal.

Metodología
Los autores derivan un modelo basado en física que vincula la directividad de la antena con el factor K de Rician y las estadísticas de retardo mediante los siguientes pasos:

1. Derivación del Factor K: El factor K de Rician (relación entre la potencia de la línea de visión y la potencia dispersa) se modela como una función de la ganancia total del sistema ($G = G_{Tx} + G_{Rx}$). Los autores introducen un factor de reverberancia ($a_K \in [0, 1]$) para tener en cuenta el crecimiento sublineal de la potencia dispersa con la ganancia, causado por reflexiones especulares en campos anisotrópicos. La relación se define como $K(G) = K_0 \cdot 10^{a_K G/10}$, donde $K_0$ es el factor K base para un enlace omnidireccional.
2. Identidad General: Se deriva una identidad matemática general que conecta $\sigma_\tau$, $\bar{\tau}$ y $K$ para cualquier Perfil de Retardo de Potencia (PDP) de dispersión arbitrario. Esta identidad revela que el par $(\bar{\tau}, \sigma_\tau)$ traza un semicírculo parametrizado por $K$, con un radio determinado exclusivamente por el coeficiente de variación ($CV_s$) del perfil de dispersión.
3. Modelo de Dispersión Exponencial: Los autores aplican este marco al caso especial motivado físicamente de un PDP de decaimiento exponencial (modelo Sabine–Rice), donde $CV_s = 1$. Esto produce expresiones explícitas de forma cerrada para $\bar{\tau}(G)$ y $\sigma_\tau(G)$.
4. Validación mediante Simulación: El modelo se valida utilizando simulaciones de Rayos de Disparo y Rebote (SBR) en MATLAB. El escenario es el salón industrial FILL en Austria (aproximadamente 57.300 m³, recinto de hormigón) a 75 GHz. Las simulaciones se ejecutaron para 100 colocaciones aleatorias de enlaces Tx/Rx utilizando cuatro configuraciones de antena: omnidireccional (0 dBi) y tres configuraciones de pares de bocinas (ganancia total de 20, 30 y 40 dBi).
5. Métodos de Extracción: Se utilizaron dos métodos independientes para extraer los factores K de los datos de simulación:
   • Método A: Extracción basada en momentos utilizando la identidad general derivada.
   • Método B: Integración directa de potencia del primer rayo de llegada frente a la suma de todos los rayos subsiguientes.

Resultados Clave

• Supresión de Retardo: El aumento de la ganancia de la antena de 0 dBi a 40 dBi suprimió la potencia de multitrayectoria en más de 40 dB dentro de 100 ns. El retardo excesivo medio ($\bar{\tau}$) disminuyó de 104,1 ns a 9,2 ns, y la dispersión de retardo RMS ($\sigma_\tau$) disminuyó de 125,2 ns a 38,2 ns.
• Validación del Modelo: El modelo derivado de tres parámetros $\{\tau_s, K_0, a_K\}$ ajustó con éxito los datos de simulación. Para el salón FILL, los parámetros extraídos fueron $\tau_s \approx 127$ ns, $K_0 = 0,32$ y $a_K = 0,40$ (con un $R^2 = 0,91$ para el ajuste del factor K).
• Comportamiento Asintótico: Los resultados confirman la predicción teórica de que, para factores K altos, $\sigma_\tau$ decae a la mitad de la tasa de $\bar{\tau}$ con respecto a la ganancia de la antena (en dB). Específicamente, $\bar{\tau} \propto 10^{-a_K G/10}$ mientras que $\sigma_\tau \propto 10^{-a_K G/20}$.
• Análisis de Discrepancias: La extracción de K basada en momentos (Método A) arrojó valores de K más altos que el método de relación de potencia (Método B) a ganancias directivas. El artículo atribuye esto a que el PDP de dispersión se desvía del modelo exponencial ($CV_s > 1$) bajo condiciones directivas, donde la dispersión difusa se suprime más que las reflexiones especulares fuertes.

Significado e Implicaciones de Diseño
El artículo afirma proporcionar el primer modelo compacto y de forma cerrada que vincula la ganancia de la antena con la dispersión de retardo en canales industriales reverberantes. Su principal significado radica en habilitar reglas de diseño para la especificación de la ganancia de la antena:

• El modelo permite a los ingenieros calcular la ganancia total mínima del sistema requerida para cumplir objetivos específicos de dispersión de retardo.
• Para el salón FILL simulado, lograr $\sigma_\tau \leq 50$ ns requiere aproximadamente 40 dBi de ganancia total, mientras que lograr $\sigma_\tau \leq 5$ ns requeriría aproximadamente 90 dBi.
• Los autores señalan que el objetivo de 5 ns es efectivamente inalcanzable solo con antenas pasivas, destacando la necesidad de ecualización complementaria o formación de haces extrema.
• Las expresiones derivadas facilitan la optimización conjunta de la ganancia de la antena y la sobrecarga del ecualizador, ya que la reducción en $\sigma_\tau$ se traduce directamente en menos taps de ecualizador y menor costo computacional en módems de portadora única.

El artículo concluye que, aunque el modelo está validado para un salón de hormigón vacío específico, el marco es general. El trabajo futuro sugerido por los autores incluye la validación experimental en salas amuebladas, la caracterización de la dependencia de la ganancia del coeficiente de variación de dispersión ($CV_s$) y la integración del modelo con modelos de canal de fábrica estandarizados (por ejemplo, 3GPP).