142

142 Rozdział 11

yh=xh(1:500); % analizowany sygnał [ f,harm]=aharm(yh,fp) ;

fprintf(¹\nkoniec analizy harmonicznych');

return;

Rys. 1 1.2. Odkształcony sinusoidalny sygnał zawierający harmoniczne parzyste i nieparzyste

Przyjęto okres próbkowania dt = 0.0005 s, czyli częstotliwość próbkowania 2000 Hz. Symulowano przebieg czasowy sygnału trwający 2 s oraz 0.25 s. Na rysunku 11.3 pokazano rozkład harmonicznych dla pełnego przebiegu czasowego 2 s, a na rysunku 11.4 - rozkład harmonicznych odpowiadający przebiegowi 0.25 s.

Analizę harmonicznych w odkształconym przebiegu przeprowadzono za pomocą funkcji aharm(). Funkcja aharm() jest zapamiętana w pliku aharm.m. Realizuje ona następujące działania:

•    uzupełnienie zerami sygnału xh, gdy liczba próbek nxh=length(xh) jest mniejsza od częstotliwości próbkowania fp,

•    wybór pierwszych N elementów sygnału, gdy jego długość nxh=Iength(xh) jest większa od częstotliwości próbkowania^,

•    wywołanie transformaty fft,

•    wyznaczenie częstotliwości w poszczególnych punktach transformaty,