验证OpenSeismoMatlab的巴特沃斯高通滤波器

内容

参考

Boore, d . m . (2005)。垫和过滤器:处理强震动数据。美国地震学会公报》,95 (2),745 - 750。

描述

核实上述引用的图1,1940年埃尔森特罗模拟记录。过滤,过滤后的位移和速度加速度。过滤是通过使用两个传递的四阶高通(即频率低于截止频率减毒)巴特沃斯滤波器的截止频率,如图1所示。

地震运动

地震数据加载

eqmotions = {“Imperial_Valley_El_Centro_9_EW”};data =负载([eqmotions {1},“.dat”]);t =数据(:1);dt = t - t (1) (2);xgtt =数据(:,2);

应用巴特沃斯高通滤波器

开关

sw =“butterworthhigh”;

巴特沃斯滤波器的顺序

边境= 4;

截止频率

方法= 0.1;

应用OpenSeismoMatlab

S1 = OpenSeismoMatlab (dt, xgtt、西南边境,方法;

过滤后的加速度

cxgtt = S1.acc;

情节的时间初始地面运动的历史

%初始化图图()举行在情节(t, 0(大小(t)),“k”,“线宽”,1)%绘制加速时间初始地面运动的历史情节(t, xgtt“k”,“线宽”,1)%完成图持有从网格在包含(的时间(秒))ylabel (Acc。(厘米/ sec2)”)

获得位移和速度时间的历史

开关

sw =“timehist”;

不要使用基线校正

baselineSw = false;

OpenSeismoMatlab适用于最初的地面运动

S2 = OpenSeismoMatlab (sw, dt, xgtt baselineSw);

OpenSeismoMatlab应用于过滤地面运动

S3 = OpenSeismoMatlab (sw, dt, cxgtt baselineSw);

绘制位移时间历程

%初始化图图()举行在情节(S3.time 0(大小(S3.time)),“k”,“线宽”,1)%画最初的地面运动的位移随时间的变化p1 =情节(S2.time S2.disp,“k”,“线宽”1);%绘制位移时间过滤后的地面运动的历史p2 =情节(S3.time S3.disp,“r”,“线宽”1);%完成图持有从网格在传奇((p1, p2) {“初始”,“过滤”})包含(的时间(秒))ylabel (“Disp (cm)。”)

绘制速度时间的历史

%初始化图图()举行在情节(S3.time 0(大小(S3.time)),“k”,“线宽”,1)%画出速度时间初始地面运动的历史p1 =情节(S2.time S2.vel,“k”,“线宽”1);%画出速度时间过滤后的地面运动的历史p2 =情节(S3.time S3.vel,“r”,“线宽”1);%完成图持有从网格在传奇((p1, p2) {“初始”,“过滤”})包含(的时间(秒))ylabel (“韦尔(厘米/秒)。)

版权

版权(c) 2018 - 2022年由乔治Papazafeiropoulos

发表与MATLAB®R2017b