Optimal Calibration of the Endpoint-corrected Hilbert Transform

Resumen Técnico: Calibración Óptima de la Transformada de Hilbert con Corrección de Extremo (Endpoint-corrected)

Planteamiento del Problema
La estimación precisa y de baja latencia de la fase instantánea de las oscilaciones es crítica para la detección y actuación en bucle cerrado, particularmente en sistemas de neuroestimulación con seguimiento de fase (phase-locked) y sistemas de control en tiempo real. El enfoque estándar consiste en construir una señal analítica mediante la Transformada de Hilbert (HT) sobre una ventana deslizante. Sin embargo, las implementaciones en tiempo discreto utilizando la Transformada Discreta de Fourier (DFT) asumen periodicidad dentro de la ventana. Cuando la señal no es perfectamente periódica a través de los límites de la ventana, esta suposición crea discontinuidades que se manifiestan como oscilaciones de Gibbs (Gibbs ringing), causando distorsiones significativas de amplitud y fase en los bordes de la ventana. Dado que los sistemas en tiempo real dependen de la muestra más reciente (el extremo o endpoint), estos artefactos de frontera corrompen directamente la estimación crítica de la fase.

Aunque se introdujo la Transformada de Hilbert con corrección de extremo (ecHT) para mitigar estos problemas mediante la aplicación de un filtro de banda estrecha causal al espectro analítico, sus distorsiones sistemáticas en el extremo carecían de un análisis fundado y de forma cerrada. Las evaluaciones previas fueron mayoritariamente empíricas, dejando sin definir la naturaleza del sesgo residual y la estrategia óptima de corrección.

Metodología
Este estudio proporciona una caracterización analítica de la ecHT basada en la DFT. Los autores derivan analíticamente el operador de extremo para una señal de tono único de banda estrecha, descomponiendo la salida de la ecHT en dos componentes distintos:

1. Una ganancia compleja determinista ($G_+$): Un término de frecuencia positiva que induce un sesgo de amplitud y fase calibrable.
2. Un término de fuga residual ($G_-$): Un término de frecuencia negativa que escapa a la máscara analítica y al filtro, creando una ondulación dependiente de la fase que varía con la fase instantánea.

Basándose en esta descomposición, el artículo deriva una calibración escalar óptima de Error Cuadrático Medio (MSE). Los autores demuestran que, para cualquier implementación de ecHT, existe un escalar complejo único $C_{opt}$ que minimiza el error entre la salida calibrada y la señal analítica ideal. Este escalar se deriva analíticamente como:
$$C_{opt} = \frac{E[\hat{Z}^* Z]}{E[|\hat{Z}|^2]}$$
donde $Z$ es el extremo verdadero y $\hat{Z}$ es la estimación de la ecHT. Geométricamente, esto representa una proyección ortogonal (un filtro de Wiener escalar) que extrae la información máximamente recuperable de forma lineal. Los autores demuestran que $C_{opt}$ depende únicamente de la ganancia de frecuencia positiva y de la fuga de frecuencia negativa, las cuales pueden precomputarse para cualquier longitud de ventana, ancho de banda, orden y frecuencia central sin necesidad de ajustes basados en datos.

Contribuciones Claves

1. Descomposición Analítica: El artículo establece que los errores de la ecHT consisten en un sesgo sistemático corregible (determinado por $G_+$) y un piso de varianza irreducible (determinado por $G_-$).
2. Calibración Óptima (cecHT): Los autores introducen la "ecHT calibrada" (cecHT), que aplica el escalar derivado $C_{opt}$ para eliminar el sesgo de fase sistemático. Esto resulta en un error de fase de media casi nula bajo condiciones estables.
3. Reglas de Diseño: El estudio proporciona límites explícitos para los errores de fase y amplitud en el extremo y relaciona los parámetros de diseño prácticos (longitud de ventana, ancho de banda/orden del filtro y desajuste de la frecuencia central) con el sesgo residual mediante un retardo de grupo en el extremo.
4. Eficiencia Computacional: La calibración añade una carga despreciable (una única multiplicación compleja) al flujo de trabajo existente de la ecHT basado en FFT de $O(L \log L)$, preservando la compatibilidad con el tiempo real.

Resultados
La eficacia de la cecHT fue validada mediante simulaciones y datos del mundo real:

• Simulaciones: En simulaciones de tono único, la cecHT redujo consistentemente el error de fase medio de aproximadamente 8.53° (ecHT no calibrada) a 0.40°, con una reducción del error máximo de 11.70° a 1.23°. La calibración hizo que el algoritmo fuera virtualmente independiente de la elección de parámetros (ancho de banda, orden del filtro, deriva de frecuencia), permitiendo ajustar los parámetros para otras restricciones (por ejemplo, robustez al ruido) sin sacrificar la precisión de fase.
• Ondas Alfa de EEG: Aplicada al conjunto de datos de Ondas Alfa de EEG, la cecHT redujo el error de fase medio de 4.9° a 0.0° ($p < 10^{-6}$) mientras preservaba el Valor de Bloqueo de Fase (PLV) y reducía significativamente el Índice de Retraso de Fase (PLI), indicando la eliminación de una asimetría de fase.
• Temblor Esencial: En el análisis de datos de temblor, la calibración redujo el error de fase medio de -15.9° a 5.2°. Al combinarlo con actualizaciones periódicas de la frecuencia central ($f_0$) para rastrear la deriva espectral, el error medio se redujo aún más a 0.0°, demostrando que los errores residuales eran principalmente inducidos por la deriva y no intrínsecos al método de calibración.

Significancia y Reivindicaciones
El artículo sostiene que la cecHT ofrece un "reemplazo práctico directo" para los flujos de trabajo de ecHT existentes, proporcionando un estimador estadísticamente robusto con reglas de diseño claras. Al separar el sesgo corregible de la fuga intrínseca, el marco ofrece tanto una herramienta práctica para una mejora inmediata como un lenguaje diagnóstico para comprender los límites de la estimación de fase.

Los autores enfatizan que, si bien la calibración elimina el sesgo sistemático codificado en la respuesta de fase del filtro, no puede eliminar el error de piso irreducible causado por la fuga espectral y la no estacionariedad. Por lo tanto, el método logra estimaciones sin sesgo bajo condiciones estables, donde el jitter restante refleja ya sea las limitaciones del estimador o la inestabilidad genuina de la señal. El estudio posiciona a la cecHT como una solución intermedia entre los filtros causales simples y los modelos adaptativos complejos, adecuada para el procesamiento de señales neurales en tiempo real, siempre que se gestionen la estabilidad de la frecuencia y la robustez ante artefactos.