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讀取檔案夾下指定擴充名的所有檔案的檔案名