考虑到有时需要由速度时程计算加速度时程，或者由位移时程计算速度时程，采用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); =peer2acc(direc,filenames{1}); =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');
研究中经常会用到将加速度数据积分成速度和位移，对前几年编写的程序重新进行了改写。为了测试和验证peer2acc函数、dirFolder函数和acc2vd函数，编写了一个简单的测试程序 acc2vd函数主要信息如下： %% Main information % Developed by Decai@hfut 2014/11/14 % Calculate velocity displacement from acceleration through integration %% Description of input output variables % ---input variables % acc: acceleration data % dt: time step % ---output variables % vel: velocity data % dis: displacement data %% Format to use % vel=acc2vd(acc,dt) % or =acc2vd(acc,dt) 函数附件： acc2vd.p 测试程序如下： clear all clc direc='test'; filenames=dirFolder(direc); =peer2acc(direc,filenames{1}); =acc2vd(acc,dt); t=dt:dt:N*dt; subplot(311) plot(t,acc); ylabel('acceleration'); subplot(312) plot(t,vel); ylabel('velocity'); subplot(313) plot(t,dis); xlabel('t') ylabel('displacement'); 运行结果：
PEER地震动记录数据库是目前应用最为广泛的数据库，为了便于利用其中的强震记录，对几年前采用Matlab编写的读取加速度数据函数进行了改写，使用了向量化的编写规则，减少了原有的循环处理，比原来的程序简化了很多。 函数输入记录所在文件夹名和记录文件名，输出一列形式的加速度数据、时间步长和数据点数，输出的加速度数据单位为g。例如：test.AT2强震记录文件存放在文件夹d:\records下，则主程序中采用acc=peer2acc('d:\records','test.AT2')或 =peer2acc(' d:\ records','test.AT2')或 =peer2acc(' d:\ records','test.AT2')，可以得到加速度数据、时间步长和数据点数。 函数说明如下： %% Main information % Developed by Decai@hfut 2014/11/13 % read acceleration data from PEER NGA database -- output acceleration % data in one column %% Description of input output variables % input variables % direc: directory of the input file % filename: input file name % output variables % acc: accelaration history data(Unit: g) % dt: time step % N: the points of acc. history %% Format to use % acc=peer2acc(direc,filename) % or =peer2acc(direc,filename) % or =peer2acc(direc,filename) %% 输入输出变量含义 %输入变量 % direc: 强震记录所在的文件夹名 % filename: 记录文件名 % 输出变量 % acc: 加速度时程数据(单位：g) % dt: 时间步长 % N: 数据点数 函数附件： peer2acc.p
标签: 强震记录
