clear all g=9.82; load hanghai_wdss -ascii % indata : tid(s), höjd(m), ay, ax, az %t h g_lag g_vert g_long (g) % z - fram/bak, x=vertikal, y=sidled (i västen) %hanghai=hanghai_wdss; %(1:1936,:); hanghai=hanghai_wdss; %(1:1936,:); tid=hanghai(:,1)-hanghai(1,1); hojd =hanghai(:,2); alat =hanghai(:,3)/g; along =hanghai(:,5)/g; avert =hanghai(:,4)/g; figure(1) subplot (2,1,1) plot(tid,avert) ylabel ('a_{vert}/g') %axis([0 60 -1 4]) grid on subplot(2,1,2) plot(tid,hojd) %axis([0 60 -15 25]) grid on ylabel ('hojd (m)') xlabel('t(s)') figure (2) subplot(3,1,1) plot(tid,avert,tid, along) ylabel ('a_{long}/g, a_{vert}/g') %axis([0 60 -1 4]) grid on subplot(3,1,2) plot(tid,along) %axis([0 15 -1 1]) grid on ylabel ('a_{long}/g') xlabel('t(s)') subplot(3,1,3) plot(tid,alat) %axis([45 60 -1 1]) grid on ylabel ('a_{lat}/g') xlabel('t(s)') figure (3) atot=(along.^2+alat.^2+avert.^2).^0.5; subplot(2,1,1) plot(tid,atot,tid,avert) %axis([0 60 -0.5 3.5]) grid on ylabel ('g_{vert}, g_{tot}') xlabel('t(s)') figure (1) subplot (2,1,1) plot(tid,atot) ylabel ('a_{tot}/g') %axis([0 60 -1 4]) axis ([0 130 0 2.5]) grid on subplot(2,1,2) axis ([0 130 -10 15]) grid on