141

Analiza harmonicznych 141

Sygnał ten został wygenerowany według następującego wzoru:

xh = -72 sin(27i50r) + 0.9-72 sin(27tlOOr) + 0.8-72 sin(27tl 500 +

+ 0.7 -72 sin(2rt200r) + 0.6-72 sin(27t250f) + 0.5-72 sin(27i300r) +

+ 0.4-72 sin(2rc350r) + 0.3-72 sin(27t400r) + 0.2-72 sin(27t450r) +    (11.14)

+ 0.1-72 sin(2n500O + 0.09yfl sin(2jt550f) + 0.08V2 sin(27t600f) +

+ 0.07-72 sin(27t650?) + 0.06-72 sin(2Ti700f) + 0.05-72 sin(27i750r).

Sygnał zawiera kolejne harmoniczne od 50 Hz do 750 Hz o malejącej amplitudzie.

W celu wyszukania harmonicznych w sygnale opracowano funkcję ahtest, która generuje odkształcony sygnał, wykonuje jego wykres czasowy, przeprowadza analizę harmonicznych za pomocą funkcji fft dla różnej długości wektora sygnału:

function ahtest

I testowanie transformaty fft tmax=2; dt=0.0005;

t=0:dt:tmax; % przedział czasu analizy fp=l/dt; % częstotliwość próbkowania p2=sqrt(2)

% odkształcony sygnał sinusoidalny xh= 1.00*p2*sin(2*pi*l*50*t)+...

0.90*p2*sin(2*pi*2*50*t)+...

0.80*p2*sin(2*pi*3*50*t) +...

0.70*p2*sin(2*pi*4*50*t) + . . . .

0.60*p2*sin(2*pi*5*50*t) + . . .

0.50*p2*sin(2*pi*6*50*t)+. . .

0.40*p2*sin(2*pi*7*50*t) + . . .

0.30*p2*sin(2*pi*8*50*t) + . . .

0.20*p2*sin(2*pi*9*50*t) + . . .

0.10*p2*sin(2*pi*10*50*t) + . . .

0.09*p2* sin(2*pi*11*50*t) + . . .

0,08*p2*sin(2*pi*12*50*t) + . . .

0.07*p2*sin(2*pi*13*50*t) + . . .

0.0 6 *p2*sin(2*pi*14*50*t) + . . .

0.05*p2*sin(2*pi* 15 * 50*t) ; trys = t(1:100) ; xhrys=xh(1:100) ; plot(trys , xhrys , 'b—’);

title('Odkształcony sygnał sinusoidalny'); xlabel('t, s'); ylabel(¹xh'); grid on; fprintf('\n Dalej? Enter'); pause;

% ANALIZA HARMONICZNYCH

[f, harm]=aharm(xh,fp); % analiza harmonicznych fprintf('\nDalej? Wciśnij ENTER'); pause;

% powtórzenie analizy dla początkowej części sygnału % liczba danych mniejsza do częstotliwości próbkowania