HangHai


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