Instantaneous Spectra Analysis of Pulse Series -- Application to Lung Sounds with Abnormalities

Resumen Técnico: Análisis de Espectros Instantáneos de Series de Pulsos con Aplicación a Sonidos Pulmonares

Planteamiento del Problema
El artículo aborda el "límite teórico" percibido de la resolución tiempo-frecuencia en el análisis de Fourier. El autor argumenta que este límite no es inherente a la teoría en sí, sino que surge de su implementación numérica, específicamente de la dependencia de la Condición de Frontera Periódica (PBC). Introducida en la década de 1930 (notablemente por Walker en 1931) para satisfacer los requisitos de integrales infinitas en los cálculos de Fourier, la PBC asume implícitamente que la serie temporal que se analiza es infinitamente periódica. Esta suposición hace que el análisis de Fourier convencional sea inadecuado para señales no periódicas, como las series de pulsos aleatorios encontradas en sonidos pulmonares anormales (por ejemplo, crepitaciones), donde la señal no se repite.

Metodología: Análisis de Fourier LXC
Para superar las limitaciones de la PBC, el artículo propone reemplazarla con la Condición de Extrapolación Lineal (LXC). Este método permite el análisis de señales no periódicas mediante la extrapolación lineal de la forma de onda en lugar de repetirla. La metodología se basa en los siguientes componentes:

1. Expansión en Series AM-FM: El método expande la serie temporal $S(t)$ en una suma de funciones complejas $e^{H_m(t)}$, donde $H_m(t)$ incorpora tanto términos de Modulación de Amplitud (AM) como de Modulación de Frecuencia (FM). Esto generaliza el concepto de "frecuencia instantánea" (van der Pol) añadiendo dinámicas de amplitud.
2. Linealización Local: Para garantizar una solución única para la expansión AM-FM (un problema señalado previamente por Daubechies), el artículo emplea el concepto de Kubo de "linealización local". La señal se aproxima alrededor de un punto temporal específico $t_k$ ignorando términos de orden superior, convirtiendo el problema en una ecuación lineal.
3. Implementación Numérica:
   • Codificación Predictiva Lineal (LPC): Se utiliza un método LPC no estándar para resolver las derivadas locales $H'_m(t_k)$, que representan la frecuencia instantánea y las tasas de crecimiento/decaimiento de la amplitud. Se evita el LPC estándar ya que depende inherentemente de la PBC.
   • Cálculo de Amplitud: Las amplitudes complejas $c_m(t_k)$ se determinan minimizando el error cuadrático entre la señal observada y el modelo.
   • Espectro Instantáneo: Se deriva un espectro instantáneo discreto $F_{disc}(f, t_k)$. A diferencia de los espectros de potencia estándar, esta formulación permite la visualización del crecimiento de la señal (intensidad positiva) y del decaimiento (intensidad negativa).
4. Filtrado: El método admite diversos filtros (por ejemplo, de amplitud, ancho espectral, potencia) para aislar componentes específicos de la señal, como la eliminación de ruido de fondo de baja frecuencia o errores computacionales.

Contribuciones Clave

• Expansión Teórica: El artículo demuestra que el análisis de Fourier LXC incluye teóricamente al análisis de Fourier PBC convencional como un caso especial, posicionando así a la LXC como una expansión natural y más amplia de la teoría de Fourier.
• Mejora de la Resolución: Al eliminar la restricción de periodicidad, el método logra una mayor resolución tiempo-frecuencia. El artículo ilustra que, mientras la resolución de la PBC está restringida por la longitud de la ventana (por ejemplo, $1/1.2$ Hz para 1.2 ciclos), la LXC puede resolver componentes individuales con alta precisión utilizando el mismo número de muestras.
• Visualización de Series de Pulsos: El método permite la descomposición de series de pulsos aleatorios en espectros de pulsos individuales, permitiendo la ensamblaje de un espectrograma que visualiza la estructura tiempo-frecuencia de eventos no periódicos.

Resultados: Aplicación a Sonidos Pulmonares
El método se aplicó a un conjunto de datos de sonidos pulmonares grabados mediante un estetoscopio digital (muestreo de 4 kHz, filtro de 100–500 Hz), incluyendo:

1. Crepitaciones (Pulsos Aleatorios): El análisis reveló que las crepitaciones consisten en una serie de pulsos aleatorios. Los espectros instantáneos mostraron que cada pulso decae rápidamente, resultando en espectros amplios. Las frecuencias centrales se dispersaron entre 300 y 700 Hz. Cabe destacar que algunas frecuencias aparecieron por encima del límite de 500 Hz del estetoscopio, lo que sugiere la alta intensidad de los pulsos.
2. Sibilancias (Frecuencia Constante): El análisis identificó una frecuencia constante alrededor de 300 Hz. Un hallazgo único fue la visualización de fases alternas de "crecimiento" (rojo) y "decaimiento" (azul) dentro del sonido de las sibilancias, formando una estructura repetitiva no visible en los espectrogramas convencionales.
3. Sonidos Pulmonares Normales: Tras filtrar el ruido de fondo de baja frecuencia, no apareció ninguna señal significativa por encima de 100 Hz, confirmando que las frecuencias detectadas en los casos anormales eran efectivamente patológicas.

Significado y Afirmaciones
El artículo afirma que el análisis de Fourier LXC visualiza con éxito la estructura tiempo-frecuencia de series de pulsos que son inaccesibles para los métodos convencionales de Transformada Rápida de Fourier de Corta Duración (STFT). Específicamente, destaca que:

• El "límite teórico" de la resolución de Fourier es un artefacto de la suposición de PBC, no una ley fundamental.
• El método proporciona una expansión de series única para señales no periódicas, permitiendo la separación de espectros de pulsos individuales.
• La capacidad de distinguir entre fases de crecimiento y decaimiento de la señal ofrece nuevas perspectivas sobre la estructura de los sonidos pulmonares, particularmente la naturaleza repetitiva de las sibilancias.

Los autores notan modestamente que, aunque el método revela componentes de alta frecuencia en las crepitaciones, el equipo de grabación actual (limitado a 500 Hz) puede ser insuficiente para capturar completamente el fenómeno, sugiriendo la necesidad de grabaciones con un rango de frecuencia más amplio para una validación futura.