加速度基线校正问题探讨 - 图文(2)

v1(i)=v1(i-1)+a1(i-1)*Delta_t/2+a1(i)*Delta_t/2;

x1(i)=x1(i-1)+v1(i-1)*Delta_t+a1(i-1)*Delta_t^2/3+a1(i)*Delta_t^2/6; end

for i=2:N

v2(i)=v2(i-1)+a2(i-1)*Delta_t/2+a2(i)*Delta_t/2;

x2(i)=x2(i-1)+v2(i-1)*Delta_t+a2(i-1)*Delta_t^2/3+a2(i)*Delta_t^2/6; end

plot(t,a0,t,a1,t,a2);grid on,set(gca, 'XTick', [0:2:56]),set(get(gca,'Children'),'linewidth',1.25),

legend('原始加速度','去平均趋势','去线性趋势'),xlabel('时间(s)'),ylabel('加速度（g）'),title('加速度时程')

figure;plot(t,v0,t,v1,t,v2);grid on,set(gca, 'XTick', [0:2:56]),set(get(gca,'Children'),'linewidth',1.5),

legend('原始加速度','去平均趋势','去线性趋势'), xlabel('时间(s)'),ylabel('速度（×10m/s）'),title('速度时程')

figure;plot(t,x0,t,x1,t,x2);grid on,set(gca, 'XTick', [0:2:56]),set(get(gca,'Children'),'linewidth',1.5),

legend('原始位移','去平均趋势','去线性趋势'),xlabel('时间(s)'),ylabel('位移（×10m）'),title('位移时程')