Das Interesse an Wavelets begann etwa um 1980 und wuchs seitdem bis heute (2002) kontinuierlich an. Einige zentrale Ideen der Wavelettheorie existierten auf die eine oder andere Art schon früher in diversen Disziplinen.

Zu Beginn der 1980er Jahre verwendeten Wissenschaftler Wavelets – frz. Ondelettes – als Alternative zur Fourier-Analyse beispielsweise bei der Analyse akustischer oder seismischer Signale. Bei der klassischen Fourierzerlegung einer Funktion ist es nicht möglich, lokale Eigenschaften einer Funktion alleine aufgrund der Kenntnis einiger Fourierkoeffizienten zu analysieren. Demgegenüber hat die Wavelet-Zerlegung einen Vorteil: Die Waveletko-effizienten spiegeln einfach, zuverlässig und präzise die Eigenschaften der zu analysierenden Funktion wider. Anstatt mit den unendlich ausgedehnten Sinus- und Cosinusfunktionen arbeitet man bei der Wavelet-Analyse mit Translationen und Dilatationen einer einzigen Grundfunktion, dem Wavelet (auch mother wavelet genannt).

Vorteilhaft ist, wenn man gut lokalisierte Funktionen verwenden kann. In diesem Fall haben kleine Änderungen im Signal nur kleine Auswirkungen auf die Koeffizienten in der Waveletdar-stellung. Daß solche gut lokalisierten Wavelets, also Funktionen mit kompaktem Träger, existieren, die darüberhinaus auch über eine gewisse Glattheit verfügen, ist nicht von vornherein klar. 1987 wurden Wavelets mit kompaktem Träger und beliebiger Glattheit konstruiert. Dies hat wesentlich zum Erfolg der Waveletmethoden beigetragen.

Historisches

Als erste orthonormale Waveletbasis kann die 1910 von A. Haar konstruierte Haar-Basis angesehen werden. Das Haar-Wavelet h für den Hilbertraum L₂(ℝ) ist definiert als \begin{eqnarray}h(x)=\left\{\begin{array}{rl}1 & {\text{f}}{\rm{\ddot {u}}}{\text{r}}\,{x}\in [{0},{1/2})\\ -1 & {\text{f}}{\rm{\ddot {u}}}{\text{r}}\,{x}\in [{1/2},{1})\\ 0 & \text{sonst}.\end{array}\right.\end{eqnarray}

Ausgehend vom mother wavelet h wird eine Orthonormalbasis des L₂(ℝ) durch Translation und Dilatation der Funktion h gewonnen. Eine weitere Orthonormalbasis des \({L}_{2}\) wurde 1923 von Walsh konstruiert. Auch diese kann im Nachhinein als eine typische Basis aus sogenannten Wavelet-Paketen (wavelet packets) interpretiert werden.

Nach diesen beiden klassischen Beispielen gab es zahlreiche Entwicklungen in Mathematik und Ingenieurwissenschaften, die letztlich zur Etablierung der Wavelettheorie beigetragen haben. Wesentliche Impulse kamen aus der Spline-Approximationstheorie, der Signal- und Bildverarbeitung (Laplace-Pyramide, Filtertechniken), sowie der harmonischen Analysis. Der Begriff Wavelet wurde ca. 1982 in Frankreich von dem Geophysiker J. Morlet und dem mathematischen Physiker A. Grossmann eingeführt.

Der französische Mathematiker Yves Meyer ist als einer derjenigen Wissenschaftler anzusehen, die die mathematisch fundierte Grundlegung der Wavelettheorie entscheidend geprägt haben. 1985 konstruierte er eine Waveletbasis mit den in mathematischer Hinsicht interessanten Eigenschaften Glattheit und Orthogonalität. Ebenfalls in dieser Zeit wurde von S. Mallat und Y. Meyer ein wichtiges Fundament und zugleich ein wesentliches Hilfmittel für die Konstruktion von Wavelets geschaffen, die multiresolution analysis oder Multiskalenzer-legung.

Die erste Waveletbasis, bestehend aus orthogonalen, beliebig glatten Funktionen mit kompakten Trägern, wurde schließlich von Ingrid Daube- chies 1987 konstruiert. Für viele Anwendungen ist die Eigenschaft der Daubechies-Wavelets, kompakten Träger zu haben, von besonderem Interesse. Speziell bewirken bei Verwendung gut lokalisierter Wavelets kleine Änderungen im Signal nur kleine Änderungen in wenigen Koeffizienten der Wavelet-darstellung.

erzeugt wird, heißt Wavelet. Man nennt ψ auch mother wavelet.

Dabei heißt \(a\in {\mathbb{R}}\backslash \{0\}\) der Skalenparameter und \(b\in {\mathbb{R}}\) der Verschiebungsparameter. Hat man ein geeignetes Wavelet gewählt, können ähnlich wie mit Fourieranalyse auch mit Hilfe einer Wavelet-zerlegung Funktionen analysiert werden. Skalenwerte a mit |a| >> 1 liefern eine breitere Funktion und dienen der Erfassung langwelliger Anteile der zu analysierenden Funktion, kleine Skalenparameter mit |a| << 1 liefern sehr schmale Wavelets und erfassen lokal präzise hochfrequente Funktionsanteile. Der Vorfaktor |a|^-1/2 dient der Normierung, damit \(\Vert \psi \Vert =1\) gilt.

Es gibt eine Vielzahl unterschiedlicher Wavelets, prinzipiell zugelassen sind alle quadratintegrierba-ren Funktionen \(\psi \in {L}_{2}({\mathbb{R}})\), die die Zulässigkeitsbe-dingung \begin{eqnarray}\displaystyle \mathop{\int}\limits_{{\mathbb{R}}}\frac{|\hat{\psi}(\xi){|}^{2}}{|\xi |}d\xi \lt \infty \end{eqnarray}

erfüllen. Eine wichtige Anwendung von Wavelets ist die Analyse und Approximation von Funktionen bzw. diskreten Signalen. Dazu müssen geeignete Basen {ψ_a,b} (Wavelet-Basis) aus Wavelets konstruiert werden. Für praktische Anwendungen wird häufig eine Diskretisierung der Waveletfunktion ψ vorgenommen. Gängig ist die Festlegung a = 2 und b = 1, in diesem speziellen Fall bildet die Familie \begin{eqnarray}{\psi}_{j,k}:={2}^{j/2}\psi ({2}^{j}.-k),j,k\in{{\mathbb{Z}}} \end{eqnarray}

eine orthonormale Basis des \({L}_{2}({\mathbb{R}})\).

Die Wavelet-Transformierte einer Funktion hängt von der Wahl des Wavelets ψ ab. Es steht eine Vielzahl verschiedener Wavelets zur Verfügung. Mit Hilfe der schnellen Wavelet-Transformation werden einer Funktion ihre Waveletkoeffizienten zugeordnet, die zur Analyse derselben verwendet werden können. Mit einer entsprechenden Rücktransformation wird die Synthese der Funktion aus den Waveletkoeffizi-enten vorgenommen. Bei der Auswahl eines geeigneten Wavelets hat man im Gegensatz zur Fourieranalyse, wo die Basisfunktionen feststehen, viel Freiheit. In der Praxis sind folgende Charakteristika von ψ von Interesse:

• Hat ψ kompakten Träger?
• Bildet die Menge

• Verfügt ψ über eine gewisse Anzahl verschwindender Momente?

Die Anzahl der verschwindenden Momente spielt bei der Kompression von Daten eine Rolle. Hat ein Wavelet eine genügend große Anzahl verschwindender Momente, so sind die Waveletkoeffizienten in glatten Bereichen der zu analysierenden Funktion klein – dort, wo Singularitäten auftreten, sind sie dagegen groß. Dieser Effekt ist für die Datenkompression, bei der kleine Koeffizienten vernachlässigt werden, von Interesse.

Viele Wavelettypen, z. B. die Daubechies-Wavelets, sind allein über ihre Koeffizienten gegeben.

Die Multiskalenzerlegung führt auch zu einem hierarchischen Schema für die Berechnung der Waveletkoeffizienten eines Eingabesignals f. In der Elektrotechnik spricht man von der Zerlegung in Teilbänder (subband coding) mit exakter Rekonstruktion. Das Vorgehen ist in der Abbildung schematisch dargestellt.

Dabei stehen H und \begin{eqnarray}\bar{H}\end{eqnarray} für Faltungen mit dem Filter \begin{eqnarray}{\{{h}_{k}\}}_{k\in {\mathbb{Z}}}\end{eqnarray}, und das Symbol 2 ↓ bezeichnet das sogenannte Downsampling zur Zerlegung des Signals f. Beim Downsampling wird nur jeder zweite Eintrag des Ausgangssignals beibehalten. Für die Rekonstruktion mit Hilfe der Filterfolge \begin{eqnarray}{\{{g}_{k}\}}_{k\in {\mathbb{Z}}}\end{eqnarray} führen wir die Abkürzungen G und \begin{eqnarray}\bar{G}\end{eqnarray} ein. In der Praxis hat \begin{eqnarray}{\{{g}_{k}\}}_{k\in {\mathbb{Z}}}\end{eqnarray} ebenso wie \begin{eqnarray}{\{{h}_{k}\}}_{k\in {\mathbb{Z}}}\end{eqnarray} meist nur wenige von Null verschiedene Einträge. Im Schema wird das Upsampling mit 2 ↑ bezeichnet, dieser Vorgang dient der Rekonstruktion des Signals. Idee dabei ist, die Menge der Werte zu vergrößern, indem vorhandene Werte den geraden Indizes in einer neuen Folge zugeordnet werden, während Folgeglieder mit ungeraden Indizes den Wert 0 erhalten.

Für jede orthonormale Basis aus Wavelets mit kompaktem Träger existieren assoziierte Paare endlicher Filter zur Teilbandzerlegung mit exakter Rekonstruktion.

Wavelet-Analyse hingegen arbeitet mit den skalierten und translatierten Versionen eines einzigen Wavelets ψ. Ein mother wavelet ψ mit kompaktem Träger lebt auf einem endlichen Intervall. Eine Sprungstelle einer Funktion kann analysiert werden, indem nur diejenigen Versionen von ψ betrachtet werden, die sie überlappen. Feinere Details können mit entsprechend fein skalierten Versionen von ψ aufgelöst werden. Die lokale Analyse einer Funktion ist mit Hilfe von nur wenigen Basisfunktionen möglich.

Mehrdimensionale Wavelets erhält man durch Tensorprodukte oder direkt aus mehrdimensionalen Skalierungsfunktionen (z. B. Boxsplines).

Zur numerischen Lösung von Gleichungssystemen mit dünn besetzter Matrix verwendet man gern iterative Verfahren wie zum Beispiel das CG-Verfahren. Die Konvergenz solcher Verfahren hängt häufig von der Konditionszahl der Systemmatrix ab. Eine gleichmäßig beschränkte Konditionszahl ist meist nicht von vornherein gegeben, also ist eine Vorkonditionierung nötig. Diese ist im Wa-veletrahmen unabhängig vom Differentialoperator zu realisieren.

Besondere Stärken von Waveletmethoden ergeben sich auch im Zusammenhang mit adaptiven Verfahren. Ausgehend von der genauen Kenntnis der Basisfunktionen, dem Vorhandensein von verschwindenden Momenten der Wavelets, sowie der Stabilität der entsprechenden Basen können Konvergenzaussagen auch für adaptive Verfahren formuliert werden.

Ein weiteres Anwendungsgebiet von Wavelets in der Numerik ist die Behandlung von Integralgleichungen. Die dabei vorkommenden Systemmatrizen sind in der Regel voll besetzt. Mit Hilfe von Wavelets mit verschwindenden Momenten können die Matrizen ausgedünnt werden, die entsprechenden Gleichungssysteme sind so effizienter zu lösen.

[4] Mallat,S.: A Wavelet Tour of Signal Processing. Academic Press New York, 1998.

[5] Meyer, Y.: Ondelettes et Operateurs. Hermann Editeurs des Sciences et des Arts Paris, 1990.