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

Unsicherheitsausbreitung in Finite-Impuls-Response-Filtern: Bewertung der Gaußschen Annahme

Grundlegende Informationen

  • Papier-ID: 2510.11384
  • Titel: Uncertainty Propagation in Finite Impulse Response Filters: Evaluating the Gaussian Assumption
  • Autoren: Jennie Couchman, Phillip Stanley-Marbell
  • Klassifizierung: eess.SP (Signalverarbeitung)
  • Veröffentlichungsdatum: 13. Oktober 2025 (arXiv-Preprint)
  • Papierlink: https://arxiv.org/abs/2510.11384v1

Zusammenfassung

Diese Forschung stellt eine häufige Annahme in der Signalverarbeitung in Frage: dass zugrunde liegende Daten numerisch einer Gaußschen Verteilung entsprechen. Obwohl diese Annahme normalerweise durch Bezugnahme auf den zentralen Grenzwertsatz gerechtfertigt wird, gilt dieser Satz nur für Summen unabhängig identisch verteilter Zufallsvariablen. Viele lineare Operationen in der Signalverarbeitung nehmen jedoch die Form gewichteter Summen an, was dazu führt, dass die Verteilung der Zufallsvariablen nicht mehr identisch ist. Dieses Papier untersucht speziell die Operation der Finite-Impuls-Response (FIR)-Filterung und prüft die FIR-Ausgabeantwort unter Vorhandensein von gleichmäßig verteiltem Quantisierungsrauschen. Es wird nachgewiesen, dass die Ausgabeunsicherheit nicht als Gaußsch angenommen werden kann, aber je nach Anwendungsszenario die Gaußsche Schätzung dennoch nützlich sein kann.

Forschungshintergrund und Motivation

Kernproblem

Diese Forschung zielt darauf ab, ein grundlegendes theoretisches Annahmeproblem in der Signalverarbeitung zu lösen: Folgt die Unsicherheitsverteilung der FIR-Filterausgabe wirklich einer Gaußschen Verteilung?

Bedeutung des Problems

  1. Genauigkeit der theoretischen Grundlagen: Die Gaußsche Annahme ist in der Signalverarbeitung allgegenwärtig und beeinflusst die Gestaltung und Leistungsbewertung nachgelagerter Verarbeitungsschritte
  2. Zuverlässigkeit praktischer Anwendungen: Falsche Verteilungsannahmen können zu ungenauen Analysen der Unsicherheitsausbreitung führen
  3. Leitfaden für Ingenieurpraxis: Bereitstellung einer genaueren theoretischen Grundlage für das praktische Systemdesign

Einschränkungen bestehender Methoden

  1. Missbrauch des zentralen Grenzwertsatzes: Ingenieure wenden den CLT häufig fälschlicherweise auf gewichtete Summen an, während der CLT nur für unabhängig identisch verteilte Variablen gilt
  2. Mangel an strenger Validierung: Bestehende Forschung mangelt es an strenger mathematischer Analyse der FIR-Filterausgabeverteilung
  3. Kluft zwischen Theorie und Praxis: Theoretische Annahmen wurden nicht ausreichend numerisch validiert

Kernbeiträge

  1. Ableitung mathematischer Ausdrücke: Erstmalige explizite Ausdruckung der Ausgabeunsicherheit des FIR-Filters als Funktion der Eingabeunsicherheit und Filterkoeffizienten
  2. Theoretischer Beweis: Strenger Beweis, dass die FIR-Filterausgabe nicht als Gaußsch angenommen werden kann
  3. Näherungsmethode: Vorschlag einer praktischen Methode zur Schätzung der Ausgabeunsicherheitsverteilung durch die dominantesten Filterkoeffizienten
  4. Experimentelle Validierung: Großflächige numerische Simulationen basierend auf echten EEG-Daten validieren die theoretischen Analyseergebnisse

Methodische Details

Aufgabendefinition

Untersuchung der Ausbreitungsgesetze der Quantisierungsunsicherheit in FIR-Filtern, speziell Analyse der echten Merkmale der Ausgabeverteilung, wenn die Eingabe gleichmäßig verteiltes Quantisierungsrauschen ist.

Theoretischer Rahmen

Quantisierungsunsicherheitsmodell

Die i-te Stichprobe des Eingangssignals wird als gleichmäßig verteilte Zufallsvariable modelliert: Xi∼U(μi−δ2,μi+δ2)X_i \sim U\left(\mu_i - \frac{\delta}{2}, \mu_i + \frac{\delta}{2}\right)

Die Wahrscheinlichkeitsdichtefunktion ist:

\frac{1}{\delta} & \text{wenn } |x_i - \mu_i| < \frac{\delta}{2} \\ 0 & \text{andernfalls} \end{cases}$$ #### FIR-Filterausgabemodellierung Die FIR-Filterausgabe wird als gewichtete Summe dargestellt: $$Y_n = \sum_{i=0}^{N} b_i X_{n-i}$$ Der Mittelwert und die Varianz der Ausgabe sind jeweils: $$\mu = \sum_{i=0}^{N} b_i \mu_{i-n}$$ $$\sigma^2 = \sum_{i=0}^{N} b_i^2 \frac{\delta^2}{12}$$ #### Ableitung der exakten Verteilung Basierend auf der geschlossenen Lösung von Kamgar-Parsi et al. ist die Wahrscheinlichkeitsdichtefunktion der Ausgabe: $$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)$$ Wobei: - $\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)$ ist die Heaviside-Schrittfunktion ### Technische Innovationen 1. **Näherung dominanter Koeffizienten**: Da die vollständige Berechnung bei großen Filtern nicht praktikabel ist, wird eine Näherung vorgeschlagen, die nur dominante Koeffizienten mit einer Amplitude über 5% des maximalen Koeffizienten verwendet 2. **Verteilungsvergleichsrahmen**: Verwendung der Jensen-Shannon-Distanz zur quantitativen Vergleichung der Unterschiede zwischen echter und theoretischer Verteilung 3. **Statistische Testmethode**: Anwendung des D'Agostino & Pearson-Tests zur Überprüfung der Nicht-Gaußschen Natur ## Experimentelle Einrichtung ### Datensatz - **Datenquelle**: Physionet EEGBCI-Datensatz - **Konfiguration**: 64-Elektroden-EEG-System, Abtastfrequenz 160 Hz, Amplitudenauflösung 1 μV - **Aufgabe**: Motorische Vorstellungsaufgaben (Vorstellung von Links- und Rechtshändigkeit) - **Simulationsumfang**: 5000 wiederholte Simulationen, die 75.735 unterschiedliche Ausgabeverteilungen erzeugen ### Filtereinstellungen - **Typ**: Bandpass-FIR-Filter - **Frequenzband**: 7-35 Hz (entsprechend Mu- und Beta-Frequenzbändern der motorischen Vorstellung) - **Parameter**: - Untere Grenzfrequenz: 7 Hz, Übergangsbandbreite 2 Hz - Obere Grenzfrequenz: 35 Hz, Übergangsbandbreite 8,75 Hz - Fensterfunktion: Hamming-Fenster - Filterlänge: 265 Samples (1,656 Sekunden) ### Bewertungsindikatoren 1. **Jensen-Shannon-Distanz**: Quantifizierung der Unterschiede zwischen Verteilungen 2. **D'Agostino & Pearson-Test**: Überprüfung der Normalität 3. **Kurtosis-Analyse**: Bewertung der Verteilungsformmerkmale 4. **Datenverhältnis innerhalb der Standardabweichung**: Vergleich der Konzentration der Verteilung ### Quantisierungsrausch-Simulation Hinzufügen von gleichmäßig verteiltem Quantisierungsrauschen: $U(-0.5\mu V, 0.5\mu V)$ ## Experimentelle Ergebnisse ### Hauptergebnisse #### Validierung der Nicht-Gaußschen Natur - **Statistische Signifikanz**: 99,7% der Ausgabeverteilungen lehnen die Nullhypothese einer Gaußschen Verteilung auf dem Niveau p < 0,05 ab - **Kurtosis-Merkmale**: Alle Ausgabeverteilungen sind platykurtisch (flachgipflig), mit einer durchschnittlichen Exzess-Kurtosis von -0,2924 #### Verteilungsanpassungsvergleich - **Gewichtete Summenverteilung vs. Gaußsche Verteilung**: - Durchschnittliche JS-Distanz zur gewichteten Summen-Theorieverteilung: 0,03184 - Durchschnittliche JS-Distanz zur Gaußschen Verteilung: 0,03661 - Die gewichtete Summenverteilung liegt deutlich näher an der echten Ausgabeverteilung #### Verteilungsmerkmalsanalyse - **Gaußsche Verteilung**: 38,29% der Werte fallen in Mittelwert ± 0,5 Standardabweichung - **Tatsächliche Ausgabe**: 36,81% der Werte fallen in Mittelwert ± 0,5 Standardabweichung - **Gewichtete Summenverteilung**: 37,00% der Werte fallen in Mittelwert ± 0,5 Standardabweichung ### Wirksamkeit der Näherung dominanter Koeffizienten Durch die Verwendung von nur 19 dominanten Koeffizienten (aus 265 ausgewählt) wurde die vollständige theoretische Verteilung erfolgreich angenähert, was die Machbarkeit der Berechnungsvereinfachungsmethode nachweist. ### Fallstudienanalyse Abbildung 2 zeigt 9 typische Fälle, die verschiedene Elektroden (C3, C4, Cz), verschiedene Zeitpunkte und verschiedene Versuche abdecken und konsistent zeigen: 1. Die tatsächliche Verteilung weicht deutlich von der Gaußschen Annahme ab 2. Die gewichtete Summen-Theorieverteilung passt gut zu den tatsächlichen Daten 3. Die JS-Distanz bestätigt diese Beobachtung quantitativ ## Verwandte Arbeiten ### Erweiterung des zentralen Grenzwertsatzes für gewichtete Summen Das Papier zitiert Forschungen von Cuzick (1995), Weber (2006), Avena & da Costa (2024) und anderen zu den Bedingungen des zentralen Grenzwertsatzes für gewichtete Summen, weist aber darauf hin, dass diese Analysen noch nicht auf FIR-Filter angewendet wurden. ### Filteranwendungen in der EEG-Signalverarbeitung - **Motorische Vorstellungsforschung**: Arbeiten von McFarland et al. (2000), ter Horst et al. (2013), Yu et al. (2022) etablierten die Assoziation zwischen Mu/Beta-Frequenzbändern und motorischer Vorstellung - **BCI-Systeme**: Das BCI2000-System von Schalk et al. (2004) bietet eine Standardplattform für die EEG-Datenerfassung ### Unsicherheitsausbreitungstheorie Basierend auf der Verteilungstheorie gewichteter gleichmäßiger Zufallsvariablen von Kamgar-Parsi et al. (1995), wird diese zum ersten Mal auf das Gebiet der digitalen Signalverarbeitung angewendet. ## Schlussfolgerungen und Diskussion ### Hauptschlussfolgerungen 1. **Theoretischer Durchbruch**: Die Unsicherheitsverteilung der FIR-Filterausgabe ist tatsächlich nicht Gaußsch, sondern eine platykurtische Verteilung mit spezifischer Form 2. **Praktischer Wert**: Obwohl nicht Gaußsch, kann die Gaußsche Näherung in bestimmten Anwendungen noch akzeptabel sein und muss für spezifische Anwendungen bewertet werden 3. **Rechnerische Optimierung**: Die Methode der dominanten Koeffizienten bietet einen praktischen Weg zur Vereinfachung der Berechnung ### Einschränkungen 1. **Spezifisches Rauschmodell**: Die Studie berücksichtigt nur gleichmäßig verteiltes Quantisierungsrauschen 2. **Beschränkung auf lineare Systeme**: Die Analyse gilt nur für lineare FIR-Filter 3. **Bewertung des Anwendungsszenarios**: Die Akzeptabilität der Gaußschen Näherung muss für spezifische Anwendungen bewertet werden ### Zukünftige Richtungen 1. **Erweiterung auf andere Rauschmodelle**: Untersuchung anderer Verteilungstypen der Eingabeunsicherheit 2. **Analyse nichtlinearer Systeme**: Erweiterung auf IIR-Filter und nichtlineare Systeme 3. **Echtzeitrechneroptimierung**: Entwicklung effizienterer Online-Unsicherheitsausbreitungsalgorithmen ## Tiefgreifende Bewertung ### Stärken #### Theoretische Beiträge 1. **Strenge mathematische Analyse**: Erstmalige Bereitstellung eines strengen mathematischen Rahmens für die Unsicherheitsausbreitung in FIR-Filtern 2. **Praktische Rechenmethode**: Die Näherungsmethode der dominanten Koeffizienten löst das Problem der Rechenkomplexität 3. **Umfassende experimentelle Validierung**: Großflächige Simulationen bieten starke empirische Unterstützung #### Methodologische Vorteile 1. **Strenge statistischer Tests**: Verwendung mehrerer statistischer Methoden zur Validierung der Zuverlässigkeit von Schlussfolgerungen 2. **Validierung mit echten Daten**: Validierung basierend auf standardisierten EEG-Datensätzen erhöht die Glaubwürdigkeit der Ergebnisse 3. **Quantitative Analyse**: Jensen-Shannon-Distanz bietet einen objektiven Standard für den Verteilungsvergleich #### Praktischer Wert 1. **Ingenieurliche Leitfäden**: Bietet Signalverarbeitungsingenieuren wichtige theoretische Korrektionen 2. **Rechnerische Machbarkeit**: Die vorgeschlagene Näherungsmethode hat praktischen Anwendungswert ### Mängel #### Theoretische Einschränkungen 1. **Einzelnes Rauschmodell**: Berücksichtigung nur gleichmäßiger Verteilung, während Rauschen in realen Systemen komplexer sein kann 2. **Lineare Annahme**: Keine Berücksichtigung der Auswirkungen nichtlinearer Verarbeitungsschritte 3. **Unabhängigkeitsannahme**: Annahme der gegenseitigen Unabhängigkeit von Eingabesamples, was möglicherweise nicht allen praktischen Situationen entspricht #### Experimentelles Design 1. **Einschränkung des Anwendungsbereichs**: Validierung nur auf EEG-Daten, mangelnde Validierung auf anderen Signaltypen 2. **Einzelner Filtertyp**: Test nur einer spezifischen FIR-Filterkonfiguration 3. **Mangel an Vergleichsmaßstäben**: Keine Vergleiche mit anderen Unsicherheitsausbreitungsmethoden #### Praktische Überlegungen 1. **Rechenkomplexität**: Obwohl eine Näherungsmethode vorgeschlagen wird, können sehr große Systeme immer noch Herausforderungen darstellen 2. **Parameterauswahl**: Der 5%-Schwellenwert für dominante Koeffizienten mangelt es an theoretischer Begründung 3. **Echtzeitfähigkeit**: Keine Diskussion der Machbarkeit in Echtzeitanwendungen ### Auswirkungen #### Akademische Beiträge 1. **Wert der theoretischen Korrektur**: Korrektur einer lange bestehenden falschen Annahme mit wichtigem akademischem Wert 2. **Interdisziplinäre Auswirkungen**: Inspirationen für Signalverarbeitung, Statistik, Metrologie und andere Bereiche 3. **Methodologische Innovation**: Die Methode der dominanten Koeffizienten könnte die Analyse anderer komplexer Systeme inspirieren #### Praktischer Wert 1. **Leitfaden für Ingenieuranwendungen**: Bietet eine genauere theoretische Grundlage für Filterdesign und Unsicherheitsanalyse 2. **Referenz für Standardentwicklung**: Könnte die Entwicklung und Überarbeitung relevanter technischer Standards beeinflussen 3. **Pädagogischer Wert**: Bietet wichtige theoretische Ergänzungen zur Signalverarbeitungsausbildung #### Reproduzierbarkeit 1. **Verwendung von Open-Source-Tools**: Basierend auf Open-Source-Tools wie MNE, leicht zu reproduzieren 2. **Öffentliche Datensätze**: Verwendung öffentlich verfügbarer Standarddatensätze erhöht die Reproduzierbarkeit 3. **Detaillierte Methodenbeschreibung**: Bietet ausreichende Implementierungsdetails ### Anwendungsszenarien #### Direkte Anwendungsbereiche 1. **Biomedizinische Signalverarbeitung**: Filterung von EEG-, ECG- und anderen biologischen Signalen 2. **Kommunikationssysteme**: Kanalfilterung und Entzerrung in der digitalen Kommunikation 3. **Audioverarbeitung**: Rauschunterdrückung und Frequenzauswahl in Audiosignalen #### Erweitertes Anwendungspotenzial 1. **Kontrollsysteme**: Unsicherheitsanalyse von Filterelementen in Rückkopplungsregelungen 2. **Bildverarbeitung**: Unsicherheitsausbreitung bei digitaler Bildfilterung 3. **Sensornetzwerke**: Unsicherheitsverwaltung bei der Datenfusion verteilter Sensoren #### Theoretischer Forschungswert 1. **Stochastische Prozesstheorie**: Neue Erkenntnisse über das Verhalten stochastischer Prozesse in linearen Systemen 2. **Messunsicherheit**: Ergänzung zur Unsicherheitsausbreitungstheorie in der Messwissenschaft 3. **Numerische Analyse**: Referenzen für die Fehlerausbreitungsanalyse in numerischen Berechnungen ## Literaturverzeichnis Das Papier zitiert 12 wichtige Referenzen, die Verteilungstheorie gewichteter Summen, EEG-Signalverarbeitung, statistische Testmethoden und andere Schlüsselbereiche abdecken und eine solide theoretische Grundlage für die Forschung bieten. Besonders bemerkenswert ist die bahnbrechende Arbeit von Kamgar-Parsi et al. von 1995, die den theoretischen Rahmen dieser Forschung begründet. --- **Gesamtbewertung**: Dies ist ein hochqualitatives Forschungspapier mit wichtigem theoretischem Wert und praktischer Bedeutung. Es korrigiert nicht nur ein häufiges Missverständnis in der Signalverarbeitung, sondern bietet auch praktische Lösungen. Obwohl es in Forschungsumfang und experimentellem Design gewisse Einschränkungen gibt, sind seine Kernbeiträge bedeutsam genug, um die Entwicklung verwandter Bereiche positiv zu fördern.