5. Analiza częstotliwościowa sygnałów losowych

5.1. Procesy losowe

Sygnał losowy można rozumieć jako proces, w którym wartość chwilowa w danym momencie czasu nie jest z góry określona, lecz opisana w sposób probabilistyczny. Aby to formalnie ująć, każdą próbkę sygnału w chwili  t traktujemy jako zmienną losową.

Zmienna losowa to wielkość, która przyjmuje wartości liczbowe w sposób losowy (stochastyczny), ale z góry określonym prawdopodobieństwem. Jeżeli zamiast jednej zmiennej losowej rozpatrujemy rodzinę zmiennych losowych uporządkowanych w czasie, otrzymujemy proces losowy  X(t) . Każda chwila czasu  t=t_0 odpowiada wówczas innej zmiennej losowej  X(t_0) , a ich zbiór opisuje cały sygnał losowy. W praktyce proces losowy można traktować jako zbiór wielu możliwych realizacji przebiegu sygnału, z których każda jest jedną realizacją tego procesu. Taki opis umożliwia przejście od analizy chwilowych wartości sygnału do badania jego własności statystycznych, uśrednionych w czasie lub po realizacjach.

Klasycznym przykładem sygnałów o charakterze losowym są sygnały biomedyczne, takie jak EEG (elektroencefalograficzne) oraz EMG (elektromiograficzne). EEG rejestruje bioelektryczną aktywność mózgu przy pomocy elektrod umieszczonych na powierzchni skóry głowy. Charakteryzuje się bardzo niską amplitudą i dużą zmiennością, oraz charakterystycznymi rytmami częstotliwościowymi, związanymi z określonymi stanami pracy mózgu (np. rytmy alfa, beta, theta). Z kolei EMG odzwierciedla aktywność elektryczną mięśni szkieletowych, powstającą w wyniku sumowania potencjałów czynnościowych włókien mięśniowych. Jest to sygnał o wyższych amplitudach niż EEG, a jego przebieg zależy od rodzaju i siły skurczu mięśnia. Oba sygnały traktuje się jako procesy losowe, gdyż ich przebieg chwilowy jest trudny do przewidzenia, natomiast możliwe jest badanie ich statystycznych własności.

Opisem probabilistycznym zmiennej losowej  X(t_0) w chwili  t=t_0 jest m.in. dystrybuanta (CDF, ang. Cumulative Distribution Function), którą definiuje się jako:

 F_X(x)=P(X(t=t_0)\leq x)

czyli prawdopodobieństwo, że w chwili  t_0 amplituda sygnału nie przekroczy wartości  x .

Dystrybuanta opisuje, jak rozkładają się wartości amplitud sygnału losowego w danym momencie czasu. Dzięki temu można np. wyznaczać progi detekcji sygnału, badać rozkład amplitud w sygnałach biomedycznych czy porównywać różne stany organizmu (np. sen vs. czuwanie w EEG) w kategoriach zmian rozkładów prawdopodobieństwa.

Dystrybuanta daje pełny opis rozkładu wartości sygnału w wersji skumulowanej – informuje o łącznym prawdopodobieństwie wystąpienia wszystkich wartości mniejszych od wartości  x . W praktyce, wygodniej jest analizować nie tyle wartości skumulowane, co dokładne prawdopodobieństwo pojawienia się amplitud w pobliżu danej wartości. Tę rolę spełnia funkcja gęstości prawdopodobieństwa (PDF, ang. Probability Density Function):

 f_X(x)=\frac{d}{dx}F_X(x)

czyli pochodna dystrybuanty. Funkcja gęstości prawdopodobieństwa pozwala wskazać, które amplitudy sygnału występują najczęściej, a które są rzadkie, co ma kluczowe znaczenie w analizie różnych procesów losowych.

Rys. 5.1. Histogram i funkcja gęstości prawdopodobieństwa dla rozkładu normalnego oraz dystrybuanta
 
Na rys. 5.1 amplitudy sygnału losowego układają się zgodnie z rozkładem normalnym. Po lewej stronie rysunku, histogram próbek wraz z dopasowaną gęstością prawdopodobieństwa (PDF) pokazuje, że najwięcej wartości skupia się wokół zera, a wartości skrajne występują coraz rzadziej. Po prawej stronie rysunku, dystrybuanta (CDF) rośnie od 0 do 1, pokazując skumulowane prawdopodobieństwo – np. dla wartości amplitudy bliskiej 0 prawdopodobieństwo, że sygnał przyjmie mniejszą wartość, wynosi około 0.5. Dopasowanie rozkładu normalnego pozwala więc opisać tło sygnału. W przypadku istotnych odchyleń od kształtu normalnego (np. obecności ogonów, asymetrii czy dodatkowych pików) można wnioskować o występowaniu artefaktów lub nietypowych zdarzeń, które wymagają zastosowania bardziej złożonej diagnostyki.