% przykład użycia fft do obliczena spektrum czyli 
% wykresu amplitudy od częstotliwości 

% - przygotowanie sygnału
Fs=1000; % czestotliwość próbkowania
Ts=1/Fs; % okres
t=0:Ts:1.3; % zadany czas trwania sygnału

% - generacja sygnału - można zadać różne sygnały w celu spraqwdzenia działania  
% tutaj: sygnał zawierający dwie czestotliwości 50 i 75Hz oraz wartość
% stałą równą 1.5
sygnal=3*sin(2*pi*50*t)+2*sin(2*pi*75*t)+1.5;
plot(t,sygnal);

% obliczenia fft
Y=fft(sygnal);
n=(length(Y)-1)
% fft daje wyniki zespolone dlatego jeśli interesuja nas tylko amplituda to
% liczymy moduł
P2 = abs(Y/length(Y)); 
% moduł z fft daje widmo symetryczne, więc potrzebujemy tylko połówki 
P1 = P2(1:n/2+1);
% pierwsza próbka to amplituda o częstotliwości 0 Hz czyli sygnał ciągły, 
% występuje w wyniku tylko raz
% dlatego jedynie resztę należy przemnożyć przez 2 aby dostać poprawną
% wartość amplitudy
P1(2:end-1) = 2*P1(2:end-1);
figure;
% wykres z prawidłowo wyskalowaną osią częstotliwości
plot(0:(Fs/n):(Fs/2),P1)