Opis softwaru falkowego

Software do falkowej dekompozycji i rekonstrukcji sygnalow z uzyciem falki Daubechies D5

Zamieszczamy tu trzy funkcje Matlab'owe (wvdecomp.m, wvreconstr.m i robfalkef.m) oraz plik pomocniczy (falkfskal.mat) potrzebny to poprawnego dzialania dwoch z tych funkcji.

Do wszystkich trzech funkcji dziala szczegolowy help np. >> help wvdecomp

WVDECOMP.M

Pierwsza funkcja jest WVDECOMP.M - wymaga ona, by w workspace'sie znajdowala sie macierz z sygnalami (w kolumnach!), ktore chcemy analizowac. W wyniku "samotlumaczaca sie" macierz komorkowa, ze wspolczynnikami falkowymi na kolejnych poziomach (poda tez zakresy czestosci). Jesli chcemy uzyc N poziomow, to dlugosc sygnalu w punktach musi dzielic sie przez 2^N (i miec 3*2^N zapasu na obu koncach, bo beda tam zaburzenia przy rekonstrukcji). Funkcja wymaga, by w bierzacym katalogu (lub w katalogu z aktywna sciezka dostepu) znajdowal sie plik FALFSKAL.MAT (potrzebne jest to do wyswietlenia rysunkow)

FALFSKAL.MAT

plik z falkami i funkcjami skalujacymi juz zapamiatanymi (zeby nie wyliczac, bo zajume to sporo czasu). Kozystaja z niego WVDECOMP.M i ROBFALKEF.M

Otrzymane wspolczynniki mozna dowolnymi metodami analizowac...

Kazdy wspolczynnik mozemy interpretowac jako iloczyn skalarny (rzut) sygnalu na poszczegolne falki czy f. skalujace (z tym, ze wspolczynniki otrzymane przez liczenie iloczynow skalarnych bylyby mniej dokladne).

ROBFALKEF.M

Funkcja, ktora w wyniku daje konkretne falki lub funkcje skalujace (np. do rysowania i ogladania na tle sygnalu). Zobacz "help robrfalkef". Indeksy J,K zgodne sa z indeksami J,K w macierzy komorkowej utworzonej funkcja WVDECOMP.M - dzieki czemu wyliczone wspolczynniki mozna wygodnie interpretowac.

Mozemy z naszych sygnalow wyrzucic sobie (lub zostawic) pewne czestosci (przefiltrowac). Robimy to wyzerowujac wspolczynniki falkowe poziomow odpowiadajacych danym czestosciom. Mozemy tez "odszumic" nasze sygnaly zostawiajac tylko te wspolczynniki, ktore uznamy na wazne (patrzac np. na srednie amplitudy jakie kolejne wspolczynniki osiagaja przy naszych sygnalach).

WVRECONSTR.M

Po filtrowaniu lub odszumianiu sygnaly mozna zrekonstruowac skryptem WVRECONSTR.M Wyniki rekonstrukcji nie powinny miec bledu wzglednego wiekszego niz 1/1000 poza koncowkami sygnalow. Koncowki, to poczatkowe i koncowe 3*2^N punktow sygnalow (gdzie N, to ilosc poziomow uzytych do dekompozycji).

PRZYKLADOWA LITERATURA (zajrzyj tez na strone naszej pracowni)

Wypych M., Kublik E., Wojdyłło P., Wróbel A. 2003. Sorting functional classes of evoked potentials by wavelets. Neuroinformatics 1: 193-202.