clc;
close all
clear all;
fp=48000; % częstotliwośc próbkowania
load DM_lab04_drgania_swobodne_dane % wczytanie danych pomiarowych
t=yt(:,3); % czas jest w trzeciej kolumnie macierzy danych pomiarowych
% zmierzone drgania są w pierwszej i drugiej kolumnie macierzy danych pomiarowych (róznią sie dokładnością) - tutaj wybrano drugą
y=yt(:,1)/(0.408/100);  % skalownie danych z 1 kanału czujnika na jednoski układu SI ([m/s^2])
y=yt(:,2)/(0.811/100);  % skalownie danych z 2 kanału czujnika na jednoski układu SI ([m/s^2]); 
%WYZNACZENIE WIDMA PRZYSPIESZENIA DRGAŃ (na podstawie szybkiej transformaty Fouriera - FFT)
L=2^19; % długosc transformaty Fouriera (to liczba częstotliwości reprezentowanych w widmie)
Y=abs(fft(y,L))/(L/2); % wyznaczenie modułu widma (abs) i wyskalowanie waplitud do rzeczywistych wartości (przez dzielenie przez (L/2))
f=0:fp/L:fp/2-fp/L; % wygenerowanie wartości osi częstotliwości

figure

subplot(2,1,1)
plot(t,y)
xlabel('czas [s]');
ylabel('Amplituda przyspieszenia y(t) [m/s^2]');
xlim([0 .4]) % SKALOWANIE ZAKRESU OSI CZASU
subplot(2,1,2)
plot(f,Y(1:L/2))
xlabel('Częstotliwość [Hz]');
ylabel('Amplituda przyspieszenia Y(f) [m/s^2]');
xlim([0 4000]) % SKALOWANIE ZAKRESU OSI CZĘSTOTLIWOŚCI