clear
clf

% Här mellan klamrarna klistrar ni in mätdata.
Data = []; 
% Slut på data

tid = Data(:,1)-Data(1,1);  % Första kolumnen är tiden, första elementet är 0 
dt = max(tid)/length(tid);
Yacc = Data(:,5)-mean(Data(:,5)); % Här är accelerationen i y-led, mätt i g
Xacc = Data(:,4)-mean(Data(:,4)); % Här är accelerationen i x-led, mätt i g
figure(1)
plot(tid,Yacc)
title('Acceleration i y-led')
xlabel('Tid (s)')
ylabel('Acceleration i g (1g=9.81m/s2)')
grid on
figure(2)
plot(tid,Xacc)
title('Acceleration i x-led')
xlabel('Tid (s)')
ylabel('Acceleration i g (1g=9.81m/s2)')
grid on

% Här görs fft-analysen------------------------------------
fs = 1/dt;
m = length(tid);
n = pow2(nextpow2(m));
FrekvAmpl_Y = fft(Yacc,n);
FrekvAmpl_X = fft(Xacc,n);
f = (0:n-1)*(fs/n);
lf = length(f);
AmplitudY = FrekvAmpl_Y.*conj(FrekvAmpl_Y)/n;
AmplitudX = FrekvAmpl_X.*conj(FrekvAmpl_X)/n;
% ---------------------------------------------------------

figure(3)
semilogx(f(2:lf/2),AmplitudY(2:lf/2))
title('Frekvensinnehåll i svängningen i y-led')
xlabel('Frekvens (Hz)')
ylabel('Amplitud i y-led')
grid on
figure(4)
semilogx(f(2:lf/2),AmplitudX(2:lf/2))
title('Frekvensinnehåll i svängningen i x-led')
xlabel('Frekvens (Hz)')
ylabel('Amplitud i x-led')
grid on