验证
这是为了验证OpenSeismoMatlab对所有可能的选项和应用程序类型都能正常工作
内容
地震数据加载
将使用埃尔森特罗地震的地震加速度时程(埃尔森特罗，1940年，埃尔森特罗终端变电站大楼)
fid = fopen (“elcentro.dat”，“r”）;文本= textscan (fid,“% % f”）;文件关闭(fid);文本t = {1};dt = t - t (1) (2);xgtt = 9.81 *文本{1,2};
没有基线校正的时间历史
sw =“timehist”；baselineSw = false;西南,S1 = OpenSeismoMatlab (dt, xgtt baselineSw);
图(1)情节(S1.time S1.disp)
图(2)情节(S1.time S1.vel)
图(3)图(S1.time S1.acc)
带有基线校正的时间历史
sw =“timehist”；baselineSw = true;S2 = OpenSeismoMatlab (sw, dt, xgtt baselineSw);
图(4)图(S2.time S2.disp)
图(5)图(S2.time S2.vel)
图(6)情节(S2.time S2.acc)
重新采样加速度时间历史从0.02秒到0.01秒。
sw =重新取样的；dti = 0.01;S3 = OpenSeismoMatlab (dt, xgtt西南,[],dti);
图(7)情节(S3.time S3.acc)
PGA
sw =“职业”；S4 = OpenSeismoMatlab (dt、xgtt sw);
S4。PGA
ans = 3.127624200000000
震动
sw =“震动”；S5 = OpenSeismoMatlab (dt、xgtt sw);
S5。震动
ans = 0.007218413820000
PGD
sw =pgd的；S6 = OpenSeismoMatlab (dt、xgtt sw);
S6。PGD
ans = 8.475736406400015 e-05
咏叹调强度和显著持续时间
sw =“咏叹调”；S7 = OpenSeismoMatlab (dt、xgtt sw);
S7。Ecum
ans = 11.251388628141965
图(8)情节(S7.time S7.EcumTH)
S7.t_5_95
Ans = 1.680000000000000 25.500000000000000
S7。道明
ans = 23.840000000000000
S7.arias
ans = 1.645606908074047
线性弹性反应谱和伪谱
sw =“西文”；ksi = 0.05;T = 0.04:0.02:1;S8 = OpenSeismoMatlab (dt, xgtt西南,[],[],ksi, T);
图(9)情节(S8.Period S8.PSa)
图(10)情节(S8.Period S8.PSv)
图(11)情节(S8.Period S8.Sd)
图(12)情节(S8.Period S8.Sv)
图(13)情节(S8.Period S8.Sa)
图(14)情节(S8.Period S8.SievABS)
图(15)情节(S8.Period S8.SievREL)
警告:忽略复数X和/或Y参数的虚部
S8。PredPSa
ans = 8.990993548805100
S8。PredPeriod
ans = 0.500000000000000
等延性响应谱和伪谱
sw =“cd”；ksi = 0.05;T = 0.04:0.02:1;μ= 2;S9 = OpenSeismoMatlab (dt, xgtt西南,[],[],ksi, T,μ);
图(15)情节(S9.Period S9.CDPSa)
图(16)情节(S9.Period S9.CDPSv)
图(17)情节(S9.Period S9.CDSd)
图(18)情节(S9.Period S9.CDSv)
图(19)情节(S9.Period S9.CDSa)
图(20)情节(S9.Period S9.fyK)
图(21)情节(S9.Period S9.muK)
图(22)情节(S9.Period S9.iterK)
傅里叶振幅谱和平均周期
sw =fas的；S10 = OpenSeismoMatlab (dt、xgtt sw);
图(23)情节(S10.freq S10.FAS)
S10。Tm
ans = 0.533389590787339
S10。调频
ans = 3.293755163661824