OpenSeismatlab可以计算以下强大运动数据处理输出:(1)加速,速度和位移的时间历史,有或没有基线校正,(2)重采样的加速度,速度和位移的时间历史,具有修改的时间步长,(3)峰接地加速,速度和位移(PGA,PGV,PGD),(4)总累积能量及其归一化时间历史,(5)根据Trifunac&Brady(1975)的显着持续时间,(6)arias强度(7)线性弹性响应光谱,(8)恒定延展性无弹性响应谱,(9)傅里叶幅度谱,(10)平均(TM)
乔治帕帕扎阿西罗罗斯(2021年)。Openseismomatlab.(https://www.mathwands.com/matlabcentral/fileexchange/67069-openseSmomatlab),Matlab中央档案交换。检索到。
谢谢你的考虑。我们对这些提案开放。
这可能是包含在GISMO内的一组很好的工具,这些工具是面向对象的,并且具有读取和编写不同地震数据格式的功能。
@ming yin samson lau:
感谢您对此提交的兴趣。如果您有其他问题,请告诉我!
嗨,抱歉似乎我想出了什么问题
我没有下载EL DAT文件并将其放在同一个路径下,以便无法打开文件
谢谢你的工作,你对我的论文提供了很多帮助
嗨,我是Matlab的新用户,因为我正在在地面运动选择和水刺上做我的jornal纸
下载加载项后,我试图在Matlab中使用它,
虽然我在加载地震数据上复制了您的前几个脚本,但它显示了
使用TextScan时出错
文件标识符无效。使用Fopen生成有效的文件标识符。
我做错了什么吗?请帮我
@ Michele Paoletti:代码已更新,以便解决此问题。感谢您的建议,并为延迟回复道歉。
谢谢你的回答乔治。TD被认为是地震波的持续时间。例如,如果我们使用Trifunac持续时间TD的定义是能量从总能量的5%到95%的时间。减少时间持续时间尊重您的公式,我们仅考虑波浪的能量而不是噪声能量(虽然非常低),而是首先存在,然后是地震波。所得到的结果是略低值尊重您的公式。我想知道你的想法。
OpenSeismomatlab在输入地震加速时间历史的整个持续时间内计算arias强度。这是代码段中没有任何限制的原因。TD你是什么意思?
嗨,必须在0到TD的时间间隔内计算Arias强度,而您发布的此行代码没有此限制:
aint2 = cumsum(a。^ 2)* pi * dt /(2 * g);
为什么?
非常感谢Ricardo A.根据您的评论,代码已修复和增强。您现在可以下载最新版本。问候!
2错误:
第188行应该是:地震= Cumtrapz(时间,XGTT);
第190行应该是:Seismic.disp = cumtrapz(时间,地震.vel);
以前,速度乘以额外的“dt”和额外的'dt ^ 2'的位移。谢谢你的工作!
此外,也许您可以添加这两个参数计算:
%平均频率
fi = f(f> 0.25&f <20);
ci = u(f> 0.25&f <20);
FM =((CI(:)'。^ 2)*(fi(:)))/(ci(:)'* ci(:));
地震.FM = FM;
%主要的时期和频率
maxi = 0;
indice = 0;
对于i = 1:长度(Seismic.psa);
如果seismic.psa(i)> maxi
maxi = seismic.psa(i);
indice =地震。Period(i);
别的
maxi = maxi;
indice = indice;
万一
结尾
地震.predspec = maxi;
地震.predperiod = indice;
地震.predfreq = 1 / maxi;
maxi = 0;
indice = 0;