Digital Beam Pattern Optimisation for the GRAO 32-m Telescope: A Comparative Analysis of FIR Filter Design Methods

Resumen Técnico: Optimización Digital del Patrón de Haz para el Telescopio GRAO de 32 m

Planteamiento del Problema
La utilidad científica del telescopio de 32 m del Observatorio de Radioastronomía de Ghana (GRAO) se ve limitada por imperfecciones residuales en su patrón de haz, incluyendo ondulaciones de lóbulos secundarios cercanos, efectos de difracción estructural y fugas de polarización cruzada. Estos artefactos, derivados de errores de superficie, bloqueo estructural y no idealidades de la alimentación, degradan la resolución angular, la sensibilidad y la precisión polarimétrica. Las estrategias de mitigación tradicionales, como el ajuste mecánico de la alimentación o el bloqueo físico, carecen de flexibilidad una vez que el telescopio está operativo y no pueden adaptarse fácilmente a los requisitos observacionales variables. El artículo identifica la necesidad de un método no invasivo y computacionalmente eficiente para suprimir artefactos de alta frecuencia espacial preservando la resolución limitada por difracción y la fidelidad polarimétrica del telescopio.

Metodología
Los autores proponen un marco que reformula la optimización del haz como un problema de procesamiento digital de señales (DSP) mediante la explotación de la dualidad entre el desplazamiento angular y la frecuencia espacial.

• Base de Datos: Se generaron simulaciones electromagnéticas de alta resolución de la antena GRAO de 32 m utilizando la suite de software GRASP (empleando Óptica Física y Teoría de la Difracción Física). Los datos consisten en una rejilla esférica (201 puntos de elevación, 181 puntos de azimut) que contiene entradas completas de la matriz de Jones compleja ($J_{qh}, J_{qv}, J_{ph}, J_{pv}$) a 5 GHz.
• Mapeo FIR: El patrón de haz angular se trata como una secuencia espacial discreta. Se aplica un filtro de Respuesta Impulsiva Finita (FIR) mediante convolución a estas muestras espaciales. Los coeficientes del filtro ($h[k]$) actúan como pesos espaciales para suprimir las altas frecuencias espaciales (lóbulos secundarios) mientras se preserva el lóbulo principal.
• Estrategias de Diseño: Se adaptaron y compararon dos métodos clásicos de diseño de FIR:
  1. Métodos Basados en Ventanas: Respuestas al impulso derivadas de filtros de paso bajo ideales que fueron truncadas y multiplicadas por funciones de ventana (Hamming, Hann, Blackman, Bartlett, Rectangular y Kaiser).
  2. Algoritmo de Parks–McClellan: Un diseño óptimo de rizado equitativo utilizando la técnica de intercambio de Remez para minimizar la desviación máxima entre la respuesta de frecuencia deseada y la real bajo una norma de Chebyshev.
• Restricciones: Se impusieron restricciones de fase lineal para asegurar un retardo de grupo constante, preservando la fase relativa entre los componentes de polarización, esencial para una calibración polarimétrica precisa.

Contribuciones Clave

1. Mapeo Formal: El artículo establece un mapeo formal de los métodos clásicos de diseño de filtros FIR hacia el dominio angular de un telescopio de plato único, tratando la optimización del haz como una tarea de filtrado espacial.
2. Comparación Cuantitativa: Proporciona una comparación rigurosa de múltiples diseños de filtros espaciales FIR aplicados a los datos de haz simulados del GRAO, evaluándolos frente a métricas relevantes para la astronomía.
3. Nuevas Métricas: Los autores introducen "métricas de filtrado espacial digital" (por ejemplo, $\theta_{50}$, respuesta residual y eficiencia dentro de una ventana angular fija) para caracterizar robustamente los haces filtrados donde las definiciones clásicas (como el HPBW y el nivel de lóbulo secundario máximo) se vuelven ambiguas o engañosas debido al ensanchamiento del lóbulo principal.
4. Vía de Implementación: El estudio describe las compensaciones prácticas entre la longitud del filtro, el costo computacional y la integración en los back-ends astronómicos, demostrando la viabilidad tanto para la calibración offline como para la implementación en tiempo real en FPGA/GPU.

Resultados

• Rendimiento de Co-Polarización: El filtrado FIR suprimió con éxito las ondulaciones de los lóbulos secundarios cercanos y suavizó el perfil del haz. Sin embargo, como era de esperar con el filtrado de paso bajo, esto resultó en un ensanchamiento del lóbulo principal.
  • Los diseños basados en ventanas (por ejemplo, Blackman, Kaiser) lograron una alta supresión de lóbulos secundarios pero ensancharon el ancho de haz de media potencia (HPBW) aproximadamente un 90–100% (de ~0.109° a ~0.21°).
  • Los diseños de Parks–Mcлеllan ofrecieron compensaciones ajustables. Una configuración de 40 toques proporcionó un resultado equilibrado, reduciendo los lóbulos secundarios pero limitando el ensanchamiento del HPBW a ~80% (a ~0.197°) y manteniendo una mejor estabilidad de fase que los filtros de ventana más largos.
• Rendimiento de Polarización Cruzada: El filtrado mitigó significamente las fluctuaciones de pequeña escala y las distorsiones localizadas en los haces de polarización cruzada ($J_{qv}, J_{ph}$). Mientras que el HPBW clásico permaneció indefinido (0°) debido a la naturaleza de los patrones de polarización cruzada, las métricas digitales mostraron una mejora sistemática en la eficiencia del haz dentro de un radio de 1° (aumentando de ~0.26 a ~0.32) y redujeron la fuga, mejorando la estabilidad polarimétrica.
• Eficiencia Computacional: Las operaciones de convolución fueron computacionalmente ligeras. En el hardware de prueba, un filtro de 65 toques se ejecutó en ~0.29 ms para un solo canal de polarización, lo que sugiere la viabilidad para el procesamiento en tiempo real a tasas de actualización de kilohertz.
• Interpretación de Métricas: El estudio aclara que las degradaciones aparentes en las métricas clásicas (por ejemplo, un aumento en el "Nivel de Lóbulo Secundario Máximo" de -23 dB a -2 dB) son artefactos de la redefinición del límite del lóbulo principal tras el filtrado, no un aumento real de la radiación fuera de eje. Las métricas digitales proporcionan una evaluación más físicamente consistente de la fuga y la eficiencia.

Significancia y Reivindicaciones
El artículo afirma que el filtrado espacial digital ofrece una vía generalizable, no invasiva y computacionalmente eficiente para el control del haz aplicable a instrumentos de plato único y de arreglo de fase (phased-array).

• Impacto Científico: El método mejora la estabilidad de la calibración y la precisión polarimétrica, fortaleciendo directamente la capacidad del telescopio para la Interferometría de Muy Larga Base (VLBI), sondeos de líneas espectrales y temporización de púlsares.
• Ventaja Operativa: A diferencia de la optimización mecánica u óptica, este enfoque no requiere modificaciones de hardware y permite el ajuste dinámico de las características del haz (resolución vs. supresión de lóbulos secundarios) para adaptarse a objetivos científicos específicos.
• Aplicabilidad Amplia: El marco une el modelado electromagnético y el procesamiento de señales, proporcionando una base matemática unificada que puede transferirse a otros telescopios de la Red VLBI Africana y futuros instrumentos.
• Modestia: Los autores señalan que estos resultados se basan en simulaciones electromagnéticas. Afirman que se requiere la validación con mediciones en el cielo (patrones de haz medidos holográficamente) para confirmar el comportamiento real del haz y las ganancias prácticas, con tales mediciones planificadas para el próximo semestre de observación. El método se presenta como un complemento práctico, en lugar de un reemplazo, de la optimización óptica o mecánica convencional.