你确定你应该只考虑加速度大于一些绝对阈值?我(不是earthquake-experienced)看起来可能你应该开始集成加速度超过阈值时,直到地震再次安静下来,事情解决(一旦块滑动也会这样做,直到地球停止移动,摩擦停止)?
集成加速度信号的获得累积位移
6视图(30天)
显示旧的评论
大家好！
我必须及时整合两次地震的加速度图获得的累计位移刚性的图块。特别是我的信号只分析的阈值设置块的加速度(得+)运动。同时,位移只需要评估,只要地面和刚体之间的相对速度不为零。我已经看到cumtrapz命令可以使用,但是,我不能解决我的问题。谁能帮我吗?
这就是我应该重现:
这是我的实际代码…我认为有一些问题与漂移的因素。
clc;清晰的所有;关闭所有;
data =负载(“accelerogramma.txt”);% carico我达缇女士
时间=数据(:1);%在secondi节奏
T =长度(时间);
加速度=数据(:,2);% accelerazione m / / s2
图
情节(时间、加速度);xlim (max(时间)[0]);% grafico输入
唉= 0.09;
* 0 + ay ay_plot =时间;
持有在
情节(时间,ay_plot); xlim (max(时间)[0]);
情节(时间、零(T)); xlim (max(时间)[0]);
dt = abs((时间(1,1)- (2,1)));
速度= cumtrapz(时间、加速度);
位移= cumtrapz(时间、速度);
图
次要情节(1,1)
情节(时间、加速度、时间、零(T)); xlim (max(时间)[0]);
次要情节(1、2)
情节(时间、速度、时间、零(T)); xlim (max(时间)[0]);
次要情节(3、1,3)
情节(时间、位移、时间、零(T)); xlim (max(时间)[0]);
acc_rel = acceleration-ay;
vel_rel = cumtrapz(时间,acc_rel);
为i = 1: T
如果vel_rel (i) < 0
vel_rel (i) = 0;
其他的
结束
结束
图
vel_rel情节(时间)
持有在
acc_rel情节(时间)
disp_rel = cumtrapz(时间,vel_rel);
图
情节(时间,disp_rel);
接受的答案
马修•诺伊
2021年10月25日
你好费德里科•
检查我的代码版本:
clc;清晰的所有;关闭所有;
data =负载(“accelerogramma.txt”);% carico我达缇女士
时间=数据(:1);%在secondi节奏
样品长度=(时间);
加速度=数据(:,2);% accelerazione m / / s2
图
情节(时间、加速度);xlim (max(时间)[0]);% grafico输入
唉= 0.09;
* 0 + ay ay_plot =时间;
持有在
情节(时间,ay_plot); xlim (max(时间)[0]);
情节(时间、零(样本));xlim (max(时间)[0]);
dt =意味着(diff(时间));
印第安纳州=找到(加速度> ay);
印第安纳州=印第安纳州(1);
acc_r = 0(样本)- ay;%相对加速度
速度= 0(1、样本);
= 1;
为ci =印第安纳州:样本
acc_r (ci) =加速度(ci) - ay;%相对加速度
速度(ci) =速度(ci-1) + * dt / 2 * (acc_r (ci) + acc_r (ci-1));%迭代trapz集成
如果速度(ci) < 0% = 0 desactivate一体化进程
= 0;
速度(ci) = 0;
结束
如果acc_r (ci) > 0% = 1激活一体化进程
= 1;
结束
结束
位移= cumtrapz(时间、速度);
图
次要情节(1,1)
情节(时间、加速度、时间、零(样本));xlim (max(时间)[0]);
次要情节(1、2)
情节(时间、速度、时间、零(样本));xlim (max(时间)[0]);
次要情节(3、1,3)
情节(时间、位移、时间、零(样本));xlim (max(时间)[0]);
