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).
				        