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('位移时程')