%TORNADO g=9.82; load tornado_wdss -ascii % indata : tid(s), höjd(m), gy, gx, gz (m/s^2) tornado_short=tornado_wdss(1:13001,:); tid=tornado_short(:,1)-tornado_short(1,1); hojd =tornado_short(:,2); alat =tornado_short(:,3)/g; along =tornado_short(:,4)/g; avert =tornado_short(:,5)/g; atot=(along.^2+alat.^2+avert.^2).^0.5;
figure(1)
subplot (2,1,1)
plot(tid,avert,tid, atot)
ylabel ('avert/g, atot/g')
axis ([0 140 0 3])
grid on
subplot (2,1,2)
plot(tid,hojd)
ylabel ('h(m)')
xlabel ('t(s)')
grid on
Grafen till höger visar höjd och "vertikal" komponent (i accelerometerns
längdriktning) av g-kraft och absolutbeloppet. Den undre grafen visar
höjddata. Vad kan du utläsa ur grafen?
Nästa graf visar alla tre komponenterna.
figure(2)
subplot (3,1,1)
plot(tid,avert)
ylabel ('avert/g,')
axis ([0 140 0 3])
grid on
subplot (3,1,2)
plot(tid,alat)
ylabel ('alat/g')
grid on
subplot (3,1,3)
plot(tid,along)
ylabel ('along/g')
xlabel ('t(s)')
grid on
