Wavelet: Definition, egenskaber og anvendelser i signalbehandling

Lær wavelets: definition, matematiske egenskaber og praktiske anvendelser i signalbehandling — fra teori til wavelettransformation og effektiv tids‑frekvens‑analyse.

Forfatter: Leandro Alegsa

wavelet er en matematisk funktion, der bruges til at beskrive en funktion eller et signal ved hjælp af enklere basisfunktioner med begrænset varighed i tid (eller rum) og frekvens. Mange signalbehandlingsopgaver formuleres ved en wavelettransformation, hvor signalet betragtes gennem en "linse" med forstørrelse bestemt af wavelettens skala. Derved afsløres kun de oplysninger, der er bestemt af den valgte wavelets form, hvilket gør wavelets særligt velegnede til at finde lokale og tidsbegrænsede træk i signaler.

Det engelske udtryk "wavelet" blev introduceret i begyndelsen af 1980'erne af de franske fysikere Jean Morlet og Alex Grossman. De brugte det franske ord "ondelette" (som betyder "lille bølge"). Senere blev dette ord overført til engelsk ved at oversætte "onde" til "wave", hvilket resulterede i betegnelsen "wavelet".

Matematisk er en wavelet en (kompleks eller reel) funktion fra Hilbert-rummet ψ ∈ L 2 ( R ) {\displaystyle \psi \in L^{2}(\mathbb {R} )}{\displaystyle \psi \in L^{2}(\mathbb {R} )}. For praktiske anvendelser bør ψ opfylde visse betingelser:

Energibetingelse

Waveletten skal have begrænset energi, dvs. være kvadratsummerbar:

∫ - ∞ ∞ | ψ ( t ) | 2 d t < ∞ {\displaystyle \int _{-\infty }^{\infty }|\psi (t)|^{2}dt<\infty } {\displaystyle \int _{-\infty }^{\infty }|\psi (t)|^{2}dt<\infty }

Tilladelighed (admissibility)

For at kunne genskabe et signal fra dets continuous wavelet transform (CWT) kræves en admissibility-betingelse. Den kan skrives som

∫ 0 ∞ | ψ ^ ( ω ) | 2 ω d ω < ∞ {\displaystyle \int _{0}^{\infty }{{|{{\hat {\psi }}(\omega )|^{2}}} \over {\omega }}d\omega <\infty } {\displaystyle \int _{0}^{\infty }{{|{\hat {\psi }}(\omega )|^{2}} \over {\omega }}d\omega <\infty }, hvor ψ ^ {\displaystyle {\hat {\psi }}}{\displaystyle {\hat {\psi }}} er en Fouriertransformation af ψ {\displaystyle \psi \,} {\displaystyle \psi \,}

Integranden ovenfor definerer ofte en konstant Cψ kaldet admissibility- eller konstanten for waveletten:

Cψ = ∫0 |Ĥψ(ω)|² / ω dω < ∞.

Nulmiddelværdi

Admissibility-betingelsen medfører, at waveletten har nulmiddelværdi (ingen DC-komponent):

∫ - ∞ ∞ ∞ ψ ( t ) d t = 0 {\displaystyle \int _{-\infty }^{\infty }\psi (t)dt=0} {\displaystyle \int _{-\infty }^{\infty }\psi (t)dt=0}

Moderwavelet og skalering/translation

Funktionen ψ {\displaystyle \psi \,}{\displaystyle \psi \,} kaldes moderwavelet. Dens oversatte (forskudte) og skalerede normaliserede versioner er givet ved

ψ a , b ( t ) = 1 a ψ ( t - b a ) {\displaystyle \psi _{a,b}(t)={\frac {1}{\sqrt {a}}}}\psi \left({{{t-b} \over {a}}\right)} {\displaystyle \psi _{a,b}(t)={\frac {1}{\sqrt {a}}}\psi \left({{t-b} \over {a}\right)}

Den oprindelige moderwavelet har parametrene a = 1 og b = 0. Parameteren b styrer translation (forskydning) og a styrer dilatation (skala):

  • Små værdier af a svarer til høje frekvenser (fin tidsopløsning).
  • Store værdier af a svarer til lave frekvenser (god frekvensopløsning).

Egenskaber, der er nyttige i praksis

  • Tids- og frekvenslokalisering: Wavelets er lokaliserede i både tid og frekvens, hvilket gør dem velegnede til at analysere ikke-stationære signaler og transiente fænomener.
  • Vanishing moments: En wavelet siges at have N vanishing moments, hvis ∫ t^k ψ(t) dt = 0 for k = 0, …, N−1. Flere vanishing moments betyder bedre evne til at repræsentere polynomielle forløb sparselt og detektere kanter/transienter.
  • Kompakt støtte (for nogle wavelets): Kort varighed i tid giver effektiv beregning og god lokal analyse.
  • Regulær og symmetri: Jevnhed (smoothness) og symmetri/antisymmetri påvirker fasereaktioner og genskabelsesnøjagtighed.
  • Ortogonale og biortogonale familier: Nogle waveletfamilier giver ortogonale baser (bruges i DWT) eller biortogonale systemer (bruges for eksempel i billedkompression).

Wavelettransformationer

To centrale typer transformationer bruges ofte:

  • Continuous Wavelet Transform (CWT): For et signal f(t) er CWT defineret ved

W_f(a,b) = ∫_{-∞}^{∞} f(t) (1/√a) overline{ψ((t−b)/a)} dt,

hvor overline{ψ} er kompleks konjugeret. CWT er redundante og giver en kontinuerlig tids-skala repræsentation, velegnet til tids-frekvens-analyse og visualisering (scalogram).

  • Rekonstruktion fra CWT sker via en inverse formel med admissibility-konstanten C_ψ:

f(t) = (1/C_ψ) ∫_{0}^{∞} ∫_{-∞}^{∞} W_f(a,b) (1/√a) ψ((t−b)/a) db da / a^2.

  • Discrete Wavelet Transform (DWT): En diskret version opnås ved at vælge et gitter i (a,b), ofte dyadisk: a = 2^j, b = k 2^j. DWT implementeres effektivt via filterbank-strukturer (højpas/lavpas) og multiresolution analysis (MRA) og anvendes bredt i kompression og hurtig signalbehandling.

Almindelige waveletfamilier

  • Haar – den simpleste wavelet; diskret og ortogonal, bruges ofte til introduktion og hurtig kompression.
  • Daubechies – kompakt understøttede ortogonale wavelets med varierende antal vanishing moments; god balance mellem tid/frekvens og kompressionsegenskaber.
  • Symlets og Coiflets – varianter af Daubechies med forbedret symmetri eller ekstra momenter.
  • Morlet – kompleks wavelet (en sinusbølge modulert af en Gauss), nyttig til kontinuerlig tids-frekvens-analyse.
  • Mexican hat (Ricker) – anden afledte af Gauss; reelt og god til detektering af singulære begivenheder.
  • Meyer – glat og med god frekvensadskillelse.

Anvendelser i signalbehandling

  • Støjreduktion (denoising): Ved at transformere signalet og anvende thresholding på koefficienterne kan støj fjernes effektivt uden at fjerne vigtige signalstrukturer.
  • Datakompression: Waveletbaserede metoder (fx JPEG2000) udnytter, at mange signaler har en sparsom repræsentation i waveletbasen.
  • Funktionsekstraktion og diagnosesystemer: Wavelets er gode til at identificere transienter, kanter og rytmiske træk – nyttigt i ECG-analyse, maskindiagnostik og seismologi.
  • Tids-frekvens-analyse: CWT giver detaljeret opdeling af, hvordan frekvensindhold varierer over tid.
  • Billedbehandling: Kantsøgning, teksturanalyse, komprimering og restaurering af billeder.
  • Feature-udtrækning til mønstergenkendelse og maskinlæring på signaler og billeder.

Valg af wavelet i praksis

Valget af moderwavelet afhænger af signalets egenskaber og formålet med analysen. Overvej:

  • Antal vanishing moments (hvor godt polynomier undertrykkes).
  • Supportens længde (kompakt støtte giver lokal analyse).
  • Symmetri eller linear fase (vigtigt i billedbehandling for at undgå faseforvrængning).
  • Regulær/jevn (smoothness) for at undgå artefakter ved genskabelse.
  • Computationsomkostning (kort filterlængde er ofte hurtigere).

Praktiske bemærkninger

  • Grænseeffekter: Ved diskret transform håndteres kanter med spejling, nul-padding eller periodisering — valget påvirker resultatet.
  • Redundans vs. effektivitet: CWT er redundant og velegnet til analyse; DWT er kompakt og hurtig til behandling og kompression.
  • Implementering: Wavelet-biblioteker findes i mange værktøjer (fx MATLAB, Python's PyWavelets), som håndterer filterbanker, DWT og CWT.

Waveletmetoder kombinerer lokal tidsanalyse med multiskala-frekvensinformation og er derfor en af de mest fleksible teknikker i moderne signalbehandling — fra støjreduktion og kompression til detektion af kortvarige fænomener og avanceret tids-frekvens-analyse.

Morlet waveletZoom
Morlet wavelet

Spørgsmål og svar

Spørgsmål: Hvad er en wavelet?


A: En wavelet er en matematisk funktion, der bruges til at nedskrive en funktion eller et signal i form af andre funktioner, der er mere enkle at studere. Den kan ses under linsen med en forstørrelse, der er givet af wavelettens skala, så vi kun kan se den information, der er bestemt af dens form.

Spørgsmål: Hvem har indført udtrykket "wavelet"?


A: Det engelske udtryk "wavelet" blev indført i begyndelsen af 1980'erne af de franske fysikere Jean Morlet og Alex Grossman, som brugte det franske ord "ondelette" (som betyder "lille bølge"). Senere blev dette ord overført til engelsk ved at oversætte "onde" til "wave", hvilket gav os "wavelet".

Spørgsmål: Hvad skal en wavelet opfylde for at kunne anvendes i praksis?


Svar: For at kunne anvendes i praksis skal en wavelet have en begrænset energi og opfylde en betingelse for at kunne anvendes. Ifølge denne betingelse skal den have en middelværdi på nul og også opfylde et integral over frekvensen, som er mindre end uendeligt.

Spørgsmål: Hvad forstås ved translation og dilatation, når der henvises til wavelets?


Svar: Translation henviser til forskydning eller flytning af moderwaveletten langs tidsaksen, mens dilatation henviser til skalering eller strækning/indskrænkning af moderwaveletten langs tidsaksen. Disse to parametre (translation og dilatation) beskrives ved henholdsvis b og a.

Spørgsmål: Hvad betyder det, at en wavelet har nul middelværdi?


Svar: Nul middelværdi betyder, at når man integrerer over alle værdier af t fra negativ uendelig til positiv uendelig, skal summen være lig med 0, dvs. ∫-∞∞∞ψ(t)dt=0 . Dette krav følger som nævnt ovenfor af selve betingelserne for tilladelighed.

Spørgsmål: Hvordan defineres moderwavelets?


A: Moderwavelets defineres som normaliserede versioner af oversatte (forskudte) og udvidede (skalerede) versioner af originale moderwavelets, som har parametrene "a" = 1 og "b" = 0 .


Søge
AlegsaOnline.com - 2020 / 2025 - License CC3