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

Riepilogo Tecnico: Analisi degli Spettri Istantanei di Serie di Impulsi con Applicazione ai Suoni Polmonari

Enunciato del Problema
Il documento affronta il percepito "limite teorico" della risoluzione tempo-frequenza nell'analisi di Fourier. L'autore sostiene che tale limite non è intrinseco alla teoria stessa, ma deriva dalla sua implementazione numerica, in particolare dalla dipendenza dalla Condizione al Contorno Periodica (PBC). Introdotta negli anni '30 (notabilmente da Walker nel 1931) per soddisfare i requisiti degli integrali infiniti nei calcoli di Fourier, la PBC assume implicitamente che la serie temporale analizzata sia infinitamente periodica. Tale assunzione rende l'analisi di Fourier convenzionale inadatta per segnali non periodici, come le serie di impulsi casuali presenti nei suoni polmonari anomali (ad esempio, i crepitii), dove il segnale non si ripete.

Metodologia: Analisi LXC-Fourier
Per superare i limiti della PBC, il documento propone di sostituirla con la Condizione di Estrapolazione Lineare (LXC). Questo metodo consente l'analisi di segnali non periodici estrapolando linearmente la forma d'onda invece di ripeterla. La metodologia si basa sui seguenti componenti:

1. Espansione in Serie AM-FM: Il metodo sviluppa la serie temporale $S(t)$ in una somma di funzioni complesse $e^{H_m(t)}$, dove $H_m(t)$ incorpora sia termini di Modulazione di Ampiezza (AM) sia di Modulazione di Frequenza (FM). Ciò generalizza il concetto di "frequenza istantanea" (van der Pol) aggiungendo le dinamiche di ampiezza.
2. Linearizzazione Locale: Per garantire una soluzione unica per l'espansione AM-FM (un problema precedentemente notato da Daubechies), il documento impiega il concetto di "linearizzazione locale" di Kubo. Il segnale è approssimato attorno a un punto temporale specifico $t_k$ ignorando i termini di ordine superiore, convertendo il problema in un'equazione lineare.
3. Implementazione Numerica:
   • Codifica Predittiva Lineare (LPC): Viene utilizzata una metodologia LPC non standard per risolvere le derivate locali $H'_m(t_k)$, che rappresentano la frequenza istantanea e i tassi di crescita/decadimento dell'ampiezza. Si evita la LPC standard poiché dipende intrinsecamente dalla PBC.
   • Calcolo dell'Ampiezza: Le ampiezze complesse $c_m(t_k)$ sono determinate minimizzando l'errore quadratico tra il segnale osservato e il modello.
   • Spettro Istantaneo: Viene derivato uno spettro istantaneo discreto $F_{disc}(f, t_k)$. A differenza degli spettri di potenza standard, questa formulazione consente la visualizzazione della crescita del segnale (intensità positiva) e del decadimento (intensità negativa).
4. Filtraggio: Il metodo supporta vari filtri (ad esempio, ampiezza, larghezza spettrale, potenza) per isolare componenti specifiche del segnale, come la rimozione del rumore di fondo a bassa frequenza o degli errori computazionali.

Contributi Chiave

• Espansione Teorica: Il documento dimostra che l'analisi LXC-Fourier include teoricamente l'analisi PBC-Fourier convenzionale come caso speciale, posizionando così la LXC come un'espansione naturale e più ampia della teoria di Fourier.
• Miglioramento della Risoluzione: Rimuovendo il vincolo di periodicità, il metodo raggiunge una risoluzione tempo-frequenza superiore. Il documento illustra che mentre la risoluzione PBC è limitata dalla lunghezza della finestra (ad esempio, $1/1.2$ Hz per 1.2 cicli), la LXC può risolvere singoli componenti con alta precisione utilizzando lo stesso numero di campioni.
• Visualizzazione delle Serie di Impulsi: Il metodo consente la decomposizione di serie di impulsi casuali in spettri di impulsi individuali, permettendo l'assemblaggio di uno spettrogramma che visualizza la struttura tempo-frequenza di eventi non periodici.

Risultati: Applicazione ai Suoni Polmonari
Il metodo è stato applicato a un set di dati di suoni polmonari registrati tramite uno stetoscopio digitale (campionamento a 4 kHz, filtro 100–500 Hz), includendo:

1. Crepitii (Impulsi Casuali): L'analisi ha rivelato che i crepitii consistono in una serie di impulsi casuali. Gli spettri istantanei hanno mostrato che ogni impulso decade rapidamente, risultando in spettri ampi. Le frequenze centrali si sono disperse tra 300 e 700 Hz. Notevolmente, alcune frequenze sono apparse al di sopra del limite di 500 Hz dello stetoscopio, suggerendo l'alta intensità degli impulsi.
2. Sibili (Frequenza Costante): L'analisi ha identificato una frequenza costante intorno ai 300 Hz. Una scoperta unica è stata la visualizzazione di fasi alternate di "crescita" (rosso) e "decadimento" (blu) all'interno del suono sibilante, formando una struttura ripetitiva non visibile negli spettrogrammi convenzionali.
3. Suoni Polmonari Normali: Dopo aver filtrato il rumore di fondo a bassa frequenza, non è apparso alcun segnale significativo al di sopra dei 100 Hz, confermando che le frequenze rilevate nei casi anomali erano effettivamente patologiche.

Significato e Affermazioni
Il documento afferma che l'analisi LXC-Fourier visualizza con successo la struttura tempo-frequenza delle serie di impulsi inaccessibili ai metodi convenzionali di Trasformata di Fourier a Finestra Mobile (STFT). Nello specifico, evidenzia che:

• Il "limite teorico" della risoluzione di Fourier è un artefatto dell'assunzione PBC, non una legge fondamentale.
• Il metodo fornisce un'espansione in serie unica per segnali non periodici, consentendo la separazione degli spettri di impulsi individuali.
• La capacità di distinguere tra fasi di crescita e decadimento del segnale offre nuove intuizioni sulla struttura dei suoni polmonari, in particolare sulla natura ripetitiva dei sibili.

Gli autori notano con modestia che, sebbene il metodo riveli componenti ad alta frequenza nei crepitii, l'attuale attrezzatura di registrazione (limitata a 500 Hz) potrebbe essere insufficiente per catturare pienamente il fenomeno, suggerendo la necessità di registrazioni con un intervallo di frequenze più ampio per una futura validazione.