clear all
close all

% USTAWIENIA POMIARU
l_sek_pom = 5;    % czas pomiaru [s]

fs        = 48000; 
lBit      = 24; 
lKan      = 1;    % tylko kanał 1
kal_kan_1 = 0.408/100;
ID        = -1;   % domyślny mikrofon

% NAGRYWANIE SYGNAŁU
recObj = audiorecorder(fs, lBit, lKan, ID);
recordblocking(recObj, l_sek_pom);
y = getaudiodata(recObj);

% OSIE CZASU I PRZESKALOWANIE
t = (0:length(y)-1)'/fs;          % czas [s]
a = y(:,1) / kal_kan_1;           % przyspieszenie [m/s^2]

% CAŁKOWANIE NUMERYCZNE -> prędkość i przemieszczenie
v = cumtrapz(t, a);               % prędkość [m/s]
v = detrend(v);
d = cumtrapz(t, v);               % przemieszczenie [m]
d = detrend(d,12);

% RYSOWANIE WYKRESÓW
figure('NumberTitle','off','Name','Przyspieszenie, prędkość i przemieszczenie','Position',[20 50 1000 600]);

subplot(3,1,1);
plot(t, a);
grid on;
title('Przyspieszenie (kanał 1)');
xlabel('t [s]');
ylabel('a [m/s^2]');

subplot(3,1,2);
plot(t, v);
grid on;
title('Prędkość');
xlabel('t [s]');
ylabel('v [m/s]');

subplot(3,1,3);
plot(t, d);
grid on;
title('Przemieszczenie');
xlabel('t [s]');
ylabel('d [m]');