The Minimum Variance Distortionless Response (MVDR) beamforming technique is widely applied in array systems to mitigate interference. However, applying MVDR to large arrays is computationally challenging; its computational complexity scales cubically with the number of antenna elements. In this paper, we introduce a scalable MVDR beamforming method tailored for massive arrays. Our approach, which is specific to scenarios where the signal of interest is below the noise floor (e.g.,~GPS), leverages the Sherman-Morrison formula, low-rank Singular Value Decomposition (SVD) approximations, and algebraic manipulation. Using our approach, we reduce the computational complexity from cubic to linear in the number of antennas. We evaluate the proposed method through simulations, comparing its computational efficiency and beamforming accuracy with the conventional MVDR approach. Our method significantly reduces the computational load while maintaining high beamforming accuracy for large-scale arrays. This solution holds promise for real-time applications of MVDR beamforming in fields like radar, sonar, and wireless communications, where massive antenna arrays are proliferating.
- ID del Artículo: 2510.14802
- Título: Un Algoritmo de Formación de Haces MVDR Escalable que es Lineal en el Número de Antenas
- Autores: Sanjaya Herath, Armin Gerami, Kevin Wagner, Ramani Duraiswami, Christopher A. Metzler
- Clasificación: eess.SP (Procesamiento de Señales)
- Fecha de Publicación: 16 de octubre de 2025 (preimpresión en arXiv)
- Enlace del Artículo: https://arxiv.org/abs/2510.14802
La técnica de formación de haces de Respuesta Mínima de Varianza sin Distorsión (MVDR) se aplica ampliamente en sistemas de arreglos para suprimir interferencias. Sin embargo, la aplicación de MVDR a arreglos de gran escala presenta desafíos computacionales, con una complejidad que se relaciona cúbicamente con el número de elementos de antena. Este artículo propone un método de formación de haces MVDR escalable para arreglos de gran escala. El método está diseñado específicamente para escenarios donde la señal de interés está por debajo del piso de ruido (como GPS), utilizando la fórmula de Sherman-Morrison, aproximación de descomposición en valores singulares (SVD) de bajo rango y manipulaciones algebraicas. Mediante este enfoque, la complejidad computacional se reduce de una relación cúbica con el número de antenas a una relación lineal. Se evaluó la eficiencia computacional y la precisión de la formación de haces del método propuesto mediante simulación, comparándolo con el método MVDR tradicional. El método reduce significativamente la carga computacional mientras mantiene una alta precisión de formación de haces para arreglos de gran escala.
El problema central que aborda esta investigación es la complejidad computacional de la formación de haces MVDR tradicional en arreglos de antenas de gran escala. Específicamente:
- Cuello de botella de complejidad computacional: El MVDR tradicional requiere calcular la inversa de la matriz de covarianza, con complejidad O(M³), donde M es el número de antenas
- Requisitos de tiempo real: En entornos dinámicos es necesario actualizar frecuentemente la matriz de covarianza, lo que dificulta la implementación en tiempo real
- Tendencia de arreglos de gran escala: En sistemas modernos de radar, sonar y comunicaciones inalámbricas, el tamaño de los arreglos de antenas continúa aumentando (de cientos a miles de antenas)
La importancia de este problema se refleja en:
- Amplitud de aplicaciones: MVDR se aplica ampliamente en detección de objetivos de radar, análisis de escenas acústicas y otros campos
- Necesidades de desarrollo tecnológico: Los arreglos de antenas de gran escala proporcionan alta resolución espacial y capacidad fuerte de supresión de interferencias
- Requisitos de procesamiento en tiempo real: Muchos escenarios de aplicación requieren procesamiento de formación de haces en tiempo real
Los métodos existentes en la literatura presentan las siguientes limitaciones:
- Enfoques algorítmicos: Las aproximaciones de bajo rango basadas en Nyström, descomposición QR, etc., aún tienen complejidad computacional relativamente alta
- Métodos distribuidos: Requieren protocolos de comunicación complejos y mecanismos de sincronización
- Métodos de aprendizaje profundo: Requieren grandes volúmenes de datos de entrenamiento con capacidad de generalización limitada
- Propone un algoritmo MVDR de complejidad lineal: Reduce la complejidad computacional de O(M³) a O(MK²), donde K≪M
- Combina múltiples técnicas matemáticas: Integra ingeniosamente la fórmula de Sherman-Morrison, aproximación SVD de bajo rango y manipulaciones algebraicas
- Optimización para escenarios específicos: Diseñado específicamente para escenarios donde la señal está por debajo del piso de ruido, como aplicaciones GPS
- Mantiene alta precisión de formación de haces: Reduce significativamente la complejidad computacional mientras mantiene un desempeño comparable al MVDR tradicional
- Proporciona marco algorítmico completo: Ofrece procedimientos completos para inicialización, actualización y cálculo de pesos
Considérese un arreglo compuesto por M antenas isotrópicas que reciben una señal de interés desde la dirección objetivo θ₀ e L señales de interferencia desde las direcciones θ₁, θ₂, ..., θₗ. El vector de señal recibido en el tiempo t es:
xt=a(θ0)s0(t)+∑l=1La(θl)sl(t)+sn(t)
donde a(θᵢ) es el vector de dirección, s₀(t) es la señal de interés, sₗ(t) es la señal de interferencia y sₙ(t) es el ruido.
La solución del formador de haces MVDR tradicional es:
w=a(θ0)HR−1a(θ0)R−1a(θ0)
Se utiliza una regla de actualización recursiva con factor de olvido α:
Rn+1=αRn+(1−α)xnxnH
Se utiliza la fórmula de Sherman-Morrison para actualizar la inversa de la matriz de covarianza:
Rn+1−1=αRn−1−α(1−α)α+(1−α)xnHRn−1xnRn−1xnxnHRn−1
Se realiza una aproximación de rango K a la matriz de covarianza:
R≈UKDKUKH
donde U_K ∈ C^{M×K} contiene los primeros K vectores propios y D_K ∈ C^{K×K} contiene los primeros K valores propios.
Mediante la aproximación de bajo rango, la regla de actualización se convierte en:
DK,n+1−1=αDK,n−1−α(1−α)α+(1−α)xnHUKDK,n−1UKHxnDK,n−1UKHxnxnHUKDK,n−1
- Estrategia de reducción de dimensionalidad: Mediante aproximación de bajo rango con K≪M se evita formar explícitamente la matriz de covarianza M×M
- Mecanismo de actualización incremental: Utiliza la fórmula de Sherman-Morrison para lograr actualización con complejidad O(MK²)
- Optimización para escenarios específicos: Para escenarios donde la señal está por debajo del piso de ruido, K típicamente puede establecerse como L+1
- Integración algorítmica: Combina orgánicamente SVD, fórmula de Sherman-Morrison y actualización recursiva
- Configuración del arreglo: Arreglo lineal uniforme (ULA), espaciamiento entre antenas de media longitud de onda
- Número de antenas: M = 50, 75, 100 (experimentos principales), extendido a 500 (pruebas de complejidad)
- Configuración de señales:
- Señal objetivo: distancia de 10 km, velocidad tangencial de 500 m/s
- Señal transmitida: pulso chirp lineal, ancho de banda 300 MHz, duración de pulso 100 ms
- SINR: -10 dB
- Frecuencia de muestreo: 1 MHz
- Parámetros del algoritmo:
- Factor de olvido α = 0.99
- Dimensión de bajo rango K = 10
- Número de instantáneas: 1000
- Ancho del lóbulo principal (MLW): Ancho angular del lóbulo principal del diagrama de radiación
- Nivel de lóbulos secundarios (SLL): Nivel de potencia de los lóbulos secundarios relativo al lóbulo principal
- Profundidad de nulo: Nivel de potencia del nulo relativo al lóbulo principal
- Tiempo de cálculo: Tiempo total para ejecutar 10,000 pasos de tiempo
- Ganancia SINR: Relación entre SINR de salida y SINR de entrada
- MVDR tradicional: Método estándar de inversión de matriz de covarianza
- Comparación de tiempo de ejecución: Pruebas en procesador AMD EPYC 7443 de 24 núcleos
| M | L | MLW (°) | | Profundidad de Nulo (dB) | | SLL (dB) | |
|---|
| | MVDR | Propuesto | MVDR | Propuesto | MVDR | Propuesto |
| 50 | 1 | 2.27 | 2.32 | -49.56 | -42.44 | -13.02 | -9.14 |
| 75 | 2 | 1.36 | 1.33 | -41.02 | -34.84 | -12.75 | -11.05 |
| 100 | 3 | 1.06 | 1.08 | -45.83 | -41.78 | -11.37 | -12.47 |
- MVDR tradicional: Complejidad O(M³), el tiempo de ejecución crece con M³
- Método propuesto: Complejidad O(MK²), el tiempo de ejecución crece linealmente con M
- Mejora de desempeño: Para arreglos con M=500, el tiempo de cálculo se reduce varios órdenes de magnitud
Los resultados experimentales muestran que el método propuesto es comparable al MVDR tradicional en los siguientes aspectos:
- Apuntamiento del lóbulo principal: Apunta correctamente el lóbulo principal hacia la dirección del objetivo
- Formación de nulos: Forma nulos efectivos en las direcciones de las señales de interferencia
- Forma general del diagrama de radiación: Altamente consistente con el MVDR tradicional
Mediante simulación de 100,000 pasos de tiempo se encontró:
- Degradación de SINR: El uso prolongado resulta en una disminución de la ganancia SINR
- Efecto de reinicialización: Después de reinicializar en el paso 50,000, la ganancia SINR se recupera
- Costo de reinicialización: El costo de reinicialización O(M²) sigue siendo inferior al método tradicional O(M³)
- Aproximación de bajo rango de Nyström: Utiliza aproximación de bajo rango para reducir carga computacional y de almacenamiento
- Método de descomposición QR: Seguimiento dinámico de cambios en voz y ruido
- MVDR-SMI: Implementación recursiva utilizando descomposición de Cholesky y transformaciones de Householder
- Algoritmos de paso de mensajes: Implementa formación de haces distribuida mediante comunicación entre nodos locales
- Método ADMM: Distribuye la carga computacional entre múltiples procesadores
- Aprendizaje por refuerzo adversarial profundo: Mejora la capacidad de formación de haces en MIMO de gran escala
- Red neuronal convolucional: Reduce la complejidad de entrenamiento y se utiliza para calibración de arreglos de antenas de gran escala
- Reducción significativa de complejidad computacional: De O(M³) a O(MK²), logrando escalabilidad lineal
- Mantiene alta precisión de formación de haces: Comparable al MVDR tradicional en métricas clave como ancho del lóbulo principal y nivel de lóbulos secundarios
- Aplicable a aplicaciones en tiempo real: Proporciona una solución viable para formación de haces MVDR en tiempo real en arreglos de gran escala
- Restricción de escenarios de aplicación: Diseñado específicamente para escenarios donde la señal está por debajo del piso de ruido, con desempeño deficiente cuando la señal está por encima del piso de ruido
- Degradación de desempeño a largo plazo: Requiere reinicialización periódica para mantener la ganancia SINR
- Costo de inicialización: El cálculo SVD inicial aún requiere complejidad O(M³)
- Extensión a escenarios donde la señal está por encima del piso de ruido
- Métodos de inicialización más eficientes: Como métodos SVD aleatorios
- Estrategia de reinicialización adaptativa: Que dispare automáticamente la reinicialización según la degradación de desempeño
- Fuerte innovación teórica: Combina ingeniosamente múltiples herramientas matemáticas para resolver problemas prácticos
- Alto valor práctico: Resuelve el cuello de botella clave de MVDR en arreglos de gran escala
- Experimentación suficiente: Verifica desde múltiples ángulos incluyendo precisión, complejidad y desempeño a largo plazo
- Completitud algorítmica: Proporciona procedimientos algorítmicos completos y detalles de implementación
- Escenarios de aplicación limitados: Solo aplicable bajo condiciones de señal específicas
- Análisis teórico insuficiente: Carece de garantías teóricas de convergencia y estabilidad
- Orientación en selección de parámetros: Falta orientación teórica para la selección del valor K
- Ausencia de verificación en hardware real: Solo contiene resultados de simulación, sin verificación en plataformas de hardware reales
- Contribución académica: Proporciona nuevas perspectivas de solución para formación de haces en arreglos de gran escala
- Valor de ingeniería: Promete aplicación en campos como radar, sonar y comunicaciones inalámbricas
- Reproducibilidad: Descripción clara del algoritmo, fácil de reproducir y mejorar
- Recepción GPS: Escenario típico donde la señal está por debajo del piso de ruido
- Detección de señales débiles: Aplicaciones que requieren fuerte supresión de interferencias
- Arreglos de antenas de gran escala: Sistemas con cientos a miles de antenas
- Requisitos de procesamiento en tiempo real: Aplicaciones sensibles a la latencia computacional
El artículo cita 21 referencias relacionadas que abarcan teoría fundamental de formación de haces, procesamiento de arreglos de gran escala, algoritmos SVD y otros aspectos, proporcionando una base teórica sólida para la investigación.
Evaluación General: Este es un artículo con importante valor práctico en el campo del procesamiento de señales. Mediante técnicas matemáticas ingeniosas resuelve el cuello de botella computacional de la formación de haces MVDR en arreglos de gran escala. Aunque presenta limitaciones como restricciones en escenarios de aplicación, proporciona una contribución valiosa al desarrollo del campo.