PEER地震动记录数据库是目前应用最为广泛的数据库,为了便于利用其中的强震记录,对几年前采用Matlab编写的读取加速度数据函数进行了改写,使用了向量化的编写规则,减少了原有的循环处理,比原来的程序简化了很多。
函数输入记录所在文件夹名和记录文件名,输出一列形式的加速度数据、时间步长和数据点数,输出的加速度数据单位为g。例如:test.AT2强震记录文件存放在文件夹d:records下,则主程序中采用acc=peer2acc('d:records','test.AT2')或[acc,dt]=peer2acc('d:records','test.AT2')或[acc,dt,N]=peer2acc('d:records','test.AT2'),可以得到加速度数据、时间步长和数据点数。
函数说明如下:
%% Main information
% Developed by [email protected] 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 [acc,dt]=peer2acc(direc,filename)
% or [acc,dt,N]=peer2acc(direc,filename)
%% 输入输出变量含义
%输入变量
% direc: 强震记录所在的文件夹名
% filename: 记录文件名
% 输出变量
% acc: 加速度时程数据(单位:g)
% dt: 时间步长
% N: 数据点数
转载本文请联系原作者获取授权,同时请注明本文来自王德才科学网博客。
链接地址:http://blog.sciencenet.cn/blog-708601-843366.html
下一篇:Matlab读取文件夹下指定扩展名的所有文件的文件名