考虑到有时需要由速度时程计算加速度时程,或者由位移时程计算速度时程,采用Matlab编写了由速度时程计算加速度时程函数,采用了五点中心差分方法。
函数说明如下:
%% Main information
% Developed by Decai@hfut 2014/11/15
% Calculate acceleration & displacement from velocity using centered five
% point derivative method
%% Description of input & output variables
% ---input variables
% vel: velocity data
% dt: time step
% ---output variables
% acc: acceleration data
%% Format to use
% acc=v2a(vel,dt)
函数附件:v2a.p
测试程序:
clear all
clc
direc='test';
filenames=dirFolder(direc);
[acc,dt,N]=peer2acc(direc,filenames{1});
[vel,dis]=acc2vd(acc,dt);
t=dt:dt:N*dt;
acc0=v2a(vel,dt);
plot(t,acc,'r');
hold on
plot(t,acc0,'k');
xlabel('t')
ylabel('acceleration');
转载本文请联系原作者获取授权,同时请注明本文来自王德才科学网博客。
链接地址:http://blog.sciencenet.cn/blog-708601-843648.html
上一篇:Matlab编写由加速度积分得到速度和位移函数