144

144 Rozdział 11

harm(2:fN) = 2*harm(2:fN);

% puste miejsca po ujemnych częstotliwościach w spektrum

harm(fN+2:N) =    [];

nh=length(harm);

harml=1;

for i=l:nh

fi=f(i); harmi=harm(i);

if fi == 50 harml=harmi; break;    end

end

fprintf('\n Amplituda harmonicznej 50 Hz wynosi:'); fprintf('\n harml = %8.4fharml);

fprintf('\n Amplitudy harmonicznych odniesiono do harml'); harm=harm/harml;

% Rozkład harmonicznych

% odniesionych do amplitudy harmonicznej 50 Hz plot(f,harm,’b-');

xlabel('f, Hz'); ylabel(’Amplituda’);

title('Amplitudy harmonicznych odniesiono do amplitudy harmonicznej 50 Hz'); grid on; return

Amplitudy harmonicznych odniesiono do amplitudy harmonicznej 50 Hz

Rys. 113. Rozkład harmonicznych uzyskany dla sygnału o przebiegu czasowym 2 s i częstotliwości próbkowania 2000 Hz