2025-11-13T05:04:10.961087

Uncertainty Propagation in Finite Impulse Response Filters: Evaluating the Gaussian Assumption

Couchman, Stanley-Marbell
A common assumption in signal processing is that underlying data numerically conforms to a Gaussian distribution. It is commonly utilized in signal processing to describe unknown additive noise in a system and is often justified by citing the central limit theorem for sums of random variables, although the central limit theorem applies only to sums of independent identically distributed random variables. However, many linear operations in signal processing take the form of weighted sums, which transforms the random variables such that their distributions are no longer identical. One such operation is a finite impulse response (FIR) filter. FIR filters are commonly used in signal processing applications as a pre-processing step. FIR output noise is generally assumed to be Gaussian. This article examines the FIR output response in the presence of uniformly distributed quantization noise. We express the FIR output uncertainty in terms of the input quantization uncertainty and filter coefficients. We show that the output uncertainty cannot be assumed to be Gaussian, but depending on the application a Gaussian estimation may still be useful. Then, we show through detailed numerical simulations that the output uncertainty distribution of the filter can be estimated through its most dominant coefficients.
academic

Propagazione dell'Incertezza nei Filtri a Risposta all'Impulso Finito: Valutazione dell'Assunzione Gaussiana

Informazioni Fondamentali

  • ID Articolo: 2510.11384
  • Titolo: Uncertainty Propagation in Finite Impulse Response Filters: Evaluating the Gaussian Assumption
  • Autori: Jennie Couchman, Phillip Stanley-Marbell
  • Classificazione: eess.SP (Elaborazione dei Segnali)
  • Data di Pubblicazione: 13 ottobre 2025 (preprint arXiv)
  • Link dell'Articolo: https://arxiv.org/abs/2510.11384v1

Riassunto

Questo studio sfida un'assunzione comune nell'elaborazione dei segnali: che i dati sottostanti seguano numericamente una distribuzione gaussiana. Sebbene questa assunzione sia solitamente giustificata facendo riferimento al teorema del limite centrale, il teorema del limite centrale si applica solo alla somma di variabili casuali indipendenti e identicamente distribuite. Tuttavia, molte operazioni lineari nell'elaborazione dei segnali assumono la forma di somme ponderate, il che fa sì che la distribuzione delle variabili casuali non sia più identica. Questo articolo studia specificamente l'operazione dei filtri a risposta all'impulso finito (FIR), esaminando la risposta dell'output FIR in presenza di rumore di quantizzazione uniformemente distribuito, dimostrando che l'incertezza dell'output non può essere assunta come gaussiana, ma l'approssimazione gaussiana potrebbe comunque essere utile a seconda dello scenario applicativo.

Contesto di Ricerca e Motivazione

Problema Centrale

Questo studio mira a risolvere un problema teorico fondamentale nell'elaborazione dei segnali: la distribuzione dell'incertezza dell'output del filtro FIR segue veramente una distribuzione gaussiana?

Importanza del Problema

  1. Accuratezza della Base Teorica: L'assunzione gaussiana è onnipresente nell'elaborazione dei segnali, influenzando la progettazione dei passaggi di elaborazione successivi e la valutazione delle prestazioni
  2. Affidabilità delle Applicazioni Pratiche: Assunzioni di distribuzione errate possono portare a un'analisi imprecisa della propagazione dell'incertezza
  3. Significato Guida per la Pratica Ingegneristica: Fornire una base teorica più accurata per la progettazione di sistemi reali

Limitazioni dei Metodi Esistenti

  1. Uso Errato del Teorema del Limite Centrale: Gli ingegneri spesso applicano erroneamente il TLC a somme ponderate, mentre il TLC si applica solo a variabili indipendenti e identicamente distribuite
  2. Mancanza di Verifica Rigorosa: La ricerca esistente manca di un'analisi matematica rigorosa della distribuzione dell'output del filtro FIR
  3. Disconnessione tra Teoria e Pratica: Le assunzioni teoriche non sono state sufficientemente verificate numericamente

Contributi Fondamentali

  1. Derivazione di Espressioni Matematiche: Prima esplicitazione della forma funzionale dell'incertezza dell'output del filtro FIR in termini di incertezza dell'input e coefficienti del filtro
  2. Prova Teorica: Dimostrazione rigorosa che l'output del filtro FIR non può essere assunto come gaussiano
  3. Metodo di Approssimazione: Proposta di un metodo pratico per stimare la distribuzione dell'incertezza dell'output utilizzando i coefficienti del filtro più significativi
  4. Verifica Sperimentale: Simulazione numerica su larga scala basata su dati EEG reali che verifica i risultati dell'analisi teorica

Spiegazione Dettagliata del Metodo

Definizione del Compito

Studiare le leggi di propagazione dell'incertezza di quantizzazione nei filtri FIR, analizzando specificamente le caratteristiche reali della distribuzione dell'output quando l'input è rumore di quantizzazione uniformemente distribuito.

Quadro Teorico

Modello di Incertezza di Quantizzazione

L'i-esimo campione del segnale di input è modellato come una variabile casuale uniformemente distribuita: Xi∼U(μi−δ2,μi+δ2)X_i \sim U\left(\mu_i - \frac{\delta}{2}, \mu_i + \frac{\delta}{2}\right)

La sua funzione di densità di probabilità è:

\frac{1}{\delta} & \text{se } |x_i - \mu_i| < \frac{\delta}{2} \\ 0 & \text{altrimenti} \end{cases}$$ #### Modellazione dell'Output del Filtro FIR L'output del filtro FIR è rappresentato come una somma ponderata: $$Y_n = \sum_{i=0}^{N} b_i X_{n-i}$$ La media e la varianza dell'output sono rispettivamente: $$\mu = \sum_{i=0}^{N} b_i \mu_{i-n}$$ $$\sigma^2 = \sum_{i=0}^{N} b_i^2 \frac{\delta^2}{12}$$ #### Derivazione della Distribuzione Esatta Basata sulla soluzione in forma chiusa di Kamgar-Parsi et al., la funzione di densità di probabilità dell'output è: $$f_{Y_n}(y_n) = \frac{(-1)^{N+1}}{N! \tilde{b}} \sum_{s_0=\pm1} \cdots \sum_{s_N=\pm1} \tilde{s} \times (y_n - \mu - S)^N \times \Theta(y_n - \mu - S)$$ Dove: - $\tilde{b} = \prod_{i=0}^{N} b_i \delta$ - $\tilde{s} = \prod_{i=0}^{N} s_i$ - $S = \sum_{i=0}^{N} s_i b_i \frac{\delta}{2}$ - $\Theta(t)$ è la funzione gradino di Heaviside ### Punti di Innovazione Tecnica 1. **Approssimazione mediante Coefficienti Dominanti**: Poiché il calcolo completo non è praticabile per filtri di ordine elevato, viene proposto di utilizzare solo i coefficienti dominanti la cui ampiezza supera il 5% del coefficiente massimo 2. **Quadro di Confronto delle Distribuzioni**: Utilizzo della distanza Jensen-Shannon per confrontare quantitativamente le differenze tra la distribuzione reale e quella teorica 3. **Metodo di Test Statistico**: Utilizzo del test di D'Agostino & Pearson per verificare la non-gaussianità ## Configurazione Sperimentale ### Dataset - **Fonte dei Dati**: Dataset EEGBCI di Physionet - **Configurazione**: Sistema EEG a 64 elettrodi, frequenza di campionamento 160 Hz, risoluzione di ampiezza 1 μV - **Compito**: Compito di immaginazione motoria (immaginazione di presa della mano sinistra e destra) - **Scala di Simulazione**: 5000 simulazioni ripetute, generando 75.735 distribuzioni di output diverse ### Configurazione del Filtro - **Tipo**: Filtro FIR passa-banda - **Banda di Frequenza**: 7-35 Hz (corrispondente alle bande mu e beta dell'immaginazione motoria) - **Parametri**: - Frequenza di taglio inferiore: 7 Hz, larghezza banda di transizione 2 Hz - Frequenza di taglio superiore: 35 Hz, larghezza banda di transizione 8,75 Hz - Funzione finestra: Finestra di Hamming - Lunghezza del filtro: 265 campioni (1,656 secondi) ### Metriche di Valutazione 1. **Distanza Jensen-Shannon**: Quantificazione della differenza tra distribuzioni 2. **Test di D'Agostino & Pearson**: Verifica della normalità 3. **Analisi della Curtosi**: Valutazione delle caratteristiche della forma della distribuzione 4. **Proporzione di Dati entro Deviazioni Standard**: Confronto del grado di concentrazione della distribuzione ### Simulazione del Rumore di Quantizzazione Aggiunta di rumore di quantizzazione uniformemente distribuito: $U(-0.5\mu V, 0.5\mu V)$ ## Risultati Sperimentali ### Risultati Principali #### Verifica della Non-Gaussianità - **Significatività Statistica**: Il 99,7% delle distribuzioni di output rifiuta l'ipotesi nulla di distribuzione gaussiana al livello p<0,05 - **Caratteristiche di Curtosi**: Tutte le distribuzioni di output sono distribuzioni platicurtiche, con curtosi in eccesso media di -0,2924 #### Confronto dell'Adattamento della Distribuzione - **Distribuzione di Somma Ponderata vs Distribuzione Gaussiana**: - Distanza JS media dalla distribuzione teorica di somma ponderata: 0,03184 - Distanza JS media dalla distribuzione gaussiana: 0,03661 - La distribuzione di somma ponderata è chiaramente più vicina alla distribuzione di output reale #### Analisi delle Caratteristiche della Distribuzione - **Distribuzione Gaussiana**: 38,29% dei valori rientra in media ± 0,5 deviazioni standard - **Output Reale**: 36,81% dei valori rientra in media ± 0,5 deviazioni standard - **Distribuzione di Somma Ponderata**: 37,00% dei valori rientra in media ± 0,5 deviazioni standard ### Efficacia dell'Approssimazione mediante Coefficienti Dominanti Utilizzando solo 19 coefficienti dominanti (selezionati da 265), è stata ottenuta con successo un'approssimazione della distribuzione teorica completa, dimostrando la fattibilità del metodo di semplificazione computazionale. ### Analisi di Casi Specifici La Figura 2 presenta 9 casi tipici, coprendo diversi elettrodi (C3, C4, Cz), diversi punti temporali e diversi trial, mostrando coerentemente: 1. La distribuzione reale si discosta chiaramente dall'assunzione gaussiana 2. La distribuzione teorica di somma ponderata si adatta bene ai dati reali 3. La distanza JS quantifica oggettivamente questa osservazione ## Lavori Correlati ### Estensione del Teorema del Limite Centrale per Somme Ponderate L'articolo cita ricerche di Cuzick (1995), Weber (2006), Avena & da Costa (2024) e altri sulle condizioni del teorema del limite centrale per somme ponderate, ma sottolinea che queste analisi non sono ancora state applicate ai filtri FIR. ### Applicazioni di Filtraggio nell'Elaborazione dei Segnali EEG - **Ricerca sull'Immaginazione Motoria**: I lavori di McFarland et al. (2000), ter Horst et al. (2013), Yu et al. (2022) hanno stabilito l'associazione tra le bande mu/beta e l'immaginazione motoria - **Sistemi BCI**: Il sistema BCI2000 di Schalk et al. (2004) fornisce una piattaforma standard per l'acquisizione di dati EEG ### Teoria della Propagazione dell'Incertezza Basata sulla teoria della distribuzione della somma di variabili casuali uniformi ponderate di Kamgar-Parsi et al. (1995), applicata per la prima volta al campo dell'elaborazione dei segnali digitali. ## Conclusioni e Discussione ### Conclusioni Principali 1. **Scoperta Teorica**: La distribuzione dell'incertezza dell'output del filtro FIR non è effettivamente gaussiana, ma una distribuzione platicurtica con una forma specifica 2. **Valore Pratico**: Sebbene non gaussiana, l'approssimazione gaussiana rimane accettabile in alcune applicazioni e deve essere valutata in base all'applicazione specifica 3. **Ottimizzazione Computazionale**: Il metodo dei coefficienti dominanti fornisce un percorso pratico di semplificazione computazionale ### Limitazioni 1. **Modello di Rumore Specifico**: Lo studio considera solo il rumore di quantizzazione uniformemente distribuito 2. **Limitazione ai Sistemi Lineari**: L'analisi si applica solo ai filtri FIR lineari 3. **Valutazione dello Scenario Applicativo**: È necessario valutare l'accettabilità dell'approssimazione gaussiana per applicazioni specifiche ### Direzioni Future 1. **Estensione ad Altri Modelli di Rumore**: Ricerca dell'incertezza di input di altri tipi di distribuzione 2. **Analisi di Sistemi Non Lineari**: Estensione ai filtri IIR e ai sistemi non lineari 3. **Ottimizzazione del Calcolo in Tempo Reale**: Sviluppo di algoritmi di propagazione dell'incertezza online più efficienti ## Valutazione Approfondita ### Punti di Forza #### Contributi Teorici 1. **Analisi Matematica Rigorosa**: Prima fornitura di un quadro matematico rigoroso per la propagazione dell'incertezza nei filtri FIR 2. **Metodo di Calcolo Pratico**: Il metodo di approssimazione mediante coefficienti dominanti risolve il problema della complessità computazionale 3. **Verifica Sperimentale Completa**: La simulazione su larga scala fornisce un forte supporto empirico #### Vantaggi Metodologici 1. **Rigore dei Test Statistici**: Utilizzo di molteplici metodi statistici per verificare l'affidabilità delle conclusioni 2. **Verifica su Dati Reali**: La verifica basata su dataset EEG standard aumenta la credibilità dei risultati 3. **Analisi Quantitativa**: La distanza Jensen-Shannon fornisce uno standard oggettivo per il confronto delle distribuzioni #### Valore Pratico 1. **Significato Guida per l'Ingegneria**: Fornisce una correzione teorica importante per gli ingegneri dell'elaborazione dei segnali 2. **Fattibilità Computazionale**: Il metodo di approssimazione proposto ha valore di applicazione pratica ### Insufficienze #### Limitazioni Teoriche 1. **Modello di Rumore Singolare**: Considera solo la distribuzione uniforme, mentre il rumore nei sistemi reali potrebbe essere più complesso 2. **Assunzione di Linearità**: Non affronta gli effetti dei componenti di elaborazione non lineare 3. **Assunzione di Indipendenza**: Assume che i campioni di input siano mutuamente indipendenti, il che potrebbe non corrispondere a alcune situazioni reali #### Progettazione Sperimentale 1. **Limitazione del Dominio Applicativo**: Verifica solo su dati EEG, mancanza di verifica su altri tipi di segnali 2. **Tipo di Filtro Singolare**: Test solo su una configurazione specifica del filtro FIR 3. **Mancanza di Benchmark Comparativi**: Nessun confronto con altri metodi di propagazione dell'incertezza #### Considerazioni di Praticità 1. **Complessità Computazionale**: Sebbene sia proposto un metodo di approssimazione, potrebbe affrontare sfide per sistemi su larga scala 2. **Scelta dei Parametri**: La soglia del 5% per i coefficienti dominanti manca di giustificazione teorica 3. **Applicabilità in Tempo Reale**: Non è stata discussa la fattibilità dell'applicazione in sistemi in tempo reale ### Impatto #### Contributi Accademici 1. **Valore di Correzione Teorica**: La correzione di un'assunzione errata di lunga data ha un importante valore accademico 2. **Impatto Interdisciplinare**: Ha implicazioni significative per l'elaborazione dei segnali, la statistica, la metrologia e altri campi 3. **Innovazione Metodologica**: Il metodo di approssimazione mediante coefficienti dominanti potrebbe ispirare l'analisi di altri sistemi complessi #### Valore Pratico 1. **Guida per l'Applicazione Ingegneristica**: Fornisce una base teorica più accurata per la progettazione di filtri e l'analisi dell'incertezza 2. **Riferimento per la Definizione di Standard**: Potrebbe influenzare la formulazione e la revisione di standard tecnici correlati 3. **Valore Educativo**: Fornisce un importante supplemento teorico all'insegnamento dell'elaborazione dei segnali #### Riproducibilità 1. **Utilizzo di Strumenti Open Source**: Basato su strumenti open source come MNE, facilitando la riproduzione 2. **Dataset Pubblici**: L'utilizzo di dataset standard pubblici aumenta la riproducibilità 3. **Descrizione Dettagliata del Metodo**: Fornisce dettagli di implementazione sufficienti ### Scenari di Applicazione #### Campi di Applicazione Diretta 1. **Elaborazione di Segnali Biomedici**: Filtraggio di segnali biologici come EEG, ECG 2. **Sistemi di Comunicazione**: Filtraggio di canale e equalizzazione nelle comunicazioni digitali 3. **Elaborazione Audio**: Riduzione del rumore e selezione della frequenza nei segnali audio #### Potenziale di Applicazione Estesa 1. **Sistemi di Controllo**: Analisi dell'incertezza nel componente di filtraggio del controllo retroazionato 2. **Elaborazione di Immagini**: Propagazione dell'incertezza nel filtraggio di immagini digitali 3. **Reti di Sensori**: Gestione dell'incertezza nella fusione di dati di sensori distribuiti #### Valore di Ricerca Teorica 1. **Teoria dei Processi Stocastici**: Fornisce nuove intuizioni sul comportamento dei processi stocastici nei sistemi lineari 2. **Scienza della Misurazione**: Fornisce un supplemento alla teoria della propagazione dell'incertezza nella scienza della misurazione 3. **Analisi Numerica**: Fornisce un riferimento per l'analisi della propagazione degli errori nel calcolo numerico ## Bibliografia L'articolo cita 12 importanti riferimenti bibliografici, coprendo la teoria della distribuzione di somme ponderate, l'elaborazione dei segnali EEG, i metodi di test statistici e altri campi chiave, fornendo una solida base teorica per la ricerca. Particolarmente degno di nota è il lavoro pioneristico di Kamgar-Parsi et al. del 1995, che ha gettato le fondamenta per il quadro teorico di questa ricerca. --- **Valutazione Complessiva**: Questo è un articolo di ricerca di alta qualità con importante valore teorico e significato pratico. Non solo corregge un errore comune nel campo dell'elaborazione dei segnali, ma fornisce anche soluzioni pratiche. Sebbene presenti alcune limitazioni nella portata della ricerca e nella progettazione sperimentale, i suoi contributi fondamentali sono sufficientemente importanti da avere un impatto positivo significativo sullo sviluppo dei campi correlati.