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

Resumo Técnico: Análise de Espectros Instantâneos de Série de Pulsos com Aplicação em Sons Pulmonares

Enunciação do Problema
O artigo aborda o "limite teórico" percebido da resolução tempo-frequência na análise de Fourier. O autor argumenta que esse limite não é inerente à própria teoria, mas decorre de sua implementação numérica, especificamente da dependência da Condição de Fronteira Periódica (PBC). Introduzida na década de 1930 (notadamente por Walker em 1931) para satisfazer os requisitos de integrais infinitas nos cálculos de Fourier, a PBC assume implicitamente que a série temporal sendo analisada é infinitamente periódica. Essa suposição torna a análise de Fourier convencional inadequada para sinais não periódicos, como séries de pulsos aleatórios encontrados em sons pulmonares anormais (por exemplo, estertores), onde o sinal não se repete.

Metodologia: Análise de Fourier-LXC
Para superar as limitações da PBC, o artigo propõe substituí-la pela Condição de Extrapolação Linear (LXC). Este método permite a análise de sinais não periódicos ao extrapolar linearmente a forma de onda, em vez de repeti-la. A metodologia é construída sobre os seguintes componentes:

1. Expansão em Série AM-FM: O método expande a série temporal $S(t)$ em uma soma de funções complexas $e^{H_m(t)}$, onde $H_m(t)$ incorpora termos de Modulação de Amplitude (AM) e Modulação de Frequência (FM). Isso generaliza o conceito de "frequência instantânea" (van der Pol) ao adicionar dinâmicas de amplitude.
2. Linearização Local: Para garantir uma solução única para a expansão AM-FM (um problema anteriormente notado por Daubechies), o artigo emprega o conceito de "linearização local" de Kubo. O sinal é aproximado em torno de um ponto temporal específico $t_k$, ignorando termos de ordem superior, convertendo o problema em uma equação linear.
3. Implementação Numérica:
   • Codificação Preditiva Linear (LPC): Um método LPC não padrão é utilizado para resolver as derivadas locais $H'_m(t_k)$, que representam a frequência instantânea e as taxas de crescimento/decadência da amplitude. A LPC padrão é evitada, pois depende inerentemente da PBC.
   • Cálculo de Amplitude: As amplitudes complexas $c_m(t_k)$ são determinadas minimizando o erro quadrático entre o sinal observado e o modelo.
   • Espectro Instantâneo: Um espectro instantâneo discreto $F_{disc}(f, t_k)$ é derivado. Diferentemente dos espectros de potência padrão, esta formulação permite a visualização do crescimento do sinal (intensidade positiva) e do decaimento (intensidade negativa).
4. Filtragem: O método suporta vários filtros (por exemplo, amplitude, largura espectral, potência) para isolar componentes específicos do sinal, como a remoção de ruído de fundo de baixa frequência ou erros computacionais.

Principais Contribuições

• Expansão Teórica: O artigo demonstra que a análise de Fourier-LXC inclui teoricamente a análise de Fourier-PBC convencional como um caso especial, posicionando assim a LXC como uma expansão natural e mais ampla da teoria de Fourier.
• Melhoria da Resolução: Ao remover a restrição de periodicidade, o método alcança maior resolução tempo-frequência. O artigo ilustra que, enquanto a resolução da PBC é restrita pelo comprimento da janela (por exemplo, $1/1,2$ Hz para 1,2 ciclos), a LXC pode resolver componentes individuais com alta precisão usando o mesmo número de amostras.
• Visualização de Série de Pulsos: O método permite a decomposição de séries de pulsos aleatórios em espectros de pulsos individuais, permitindo a montagem de um espectrograma que visualiza a estrutura tempo-frequência de eventos não periódicos.

Resultados: Aplicação em Sons Pulmonares
O método foi aplicado a um conjunto de dados de sons pulmonares gravidos por um estetoscópio digital (amostragem de 4 kHz, filtro de 100–500 Hz), incluindo:

1. Estertores (Pulsos Aleatórios): A análise revelou que os estertores consistem em uma série de pulsos aleatórios. Os espectros instantâneos mostraram que cada pulso decai rapidamente, resultando em espectros largos. As frequências centrais dispersaram-se entre 300 e 700 Hz. Notavelmente, algumas frequências apareceram acima do limite de 500 Hz do estetoscópio, sugerindo a alta intensidade dos pulsos.
2. Sibilos (Frequência Constante): A análise identificou uma frequência constante em torno de 300 Hz. Uma descoberta única foi a visualização de fases alternadas de "crescimento" (vermelho) e "decaimento" (azul) dentro do som do sibilos, formando uma estrutura repetitiva não visível em espectrogramas convencionais.
3. Sons Pulmonares Normais: Após filtrar o ruído de fundo de baixa frequência, nenhum sinal significativo apareceu acima de 100 Hz, confirmando que as frequências detectadas nos casos anormais eram de fato patológicas.

Significado e Alegações
O artigo afirma que a análise de Fourier-LXC visualiza com sucesso a estrutura tempo-frequência de séries de pulsos que são inacessíveis aos métodos convencionais de Transformada de Fourier de Curta Janela (STFT). Especificamente, destaca-se que:

• O "limite teórico" da resolução de Fourier é um artefato da suposição da PBC, não uma lei fundamental.
• O método fornece uma expansão de série única para sinais não periódicos, permitindo a separação de espectros de pulsos individuais.
• A capacidade de distinguir entre fases de crescimento e decaimento do sinal oferece novas perspectivas sobre a estrutura dos sons pulmonares, particularmente a natureza repetitiva dos sibilos.

Os autores observam modestamente que, embora o método revele componentes de alta frequência nos estertores, o equipamento de gravação atual (limitado a 500 Hz) pode ser insuficiente para capturar totalmente o fenômeno, sugerindo a necessidade de gravações com uma faixa de frequência mais ampla para validação futura.