%z,x,y方向的傳感器標定值分別為10.19，10.22，10.74
Fs=1000;

%x=data(1:45000,2)/50/10.22;
%y=data(1:45000,3)/50/10.74;
%z=data(1:45000,4)/50/10.19;

x=data(:,2)/10.22;
y=data(:,3)/10.74;
z=data(:,4)/10.19;

x=data(:,2);
y=data(:,3);
z=data(:,4);


subplot(3,1,1);
[pz,fz] = filtbank(z,Fs); 
hz = bankdisp(pz,fz);
title('座椅Z方向1/3倍頻程分析');
wz=[0.63 0.71 0.8 0.9 1 1 1 1  0.8 0.63 0.5 0.4 0.315 0.25 0.2 0.16 0.125 0.1];


subplot(3,1,2);
[px,fx] = filtbank(x,Fs); 
hx = bankdisp(px,fx);
title('座椅X方向1/3倍頻程分析');
wx=[1 1 0.8 0.63 0.50 0.40 0.315 0.25 0.20 0.16 0.125 0.10 0.08 0.063 0.05 0.04 0.0315 0.025];


subplot(3,1,3);
[py,fy] = filtbank(y,Fs); 
hy = bankdisp(py,fy);
title('座椅Y方向1/3倍頻程分析');
wy=[1 1 0.8 0.63 0.50 0.40 0.315 0.25 0.20 0.16 0.125 0.10 0.08 0.063 0.05 0.04 0.0315 0.025];

wz2=wz.^2;
pz2=pz.^2;
az2=wz2.*pz2;
az2sum=sum(az2);
az=sqrt(az2sum);


wx2=wx.^2;
px2=px.^2;
ax2=wx2.*px2;
ax2sum=sum(ax2);
ax=sqrt(ax2sum);


wy2=wy.^2;
py2=py.^2;
ay2=wy2.*py2;
ay2sum=sum(ay2);
ay=sqrt(ay2sum);


awo=((1.4*ax)^2+(1.4*ay)^2+az^2)^(1/2);