时域积分：地震加速度时程曲线转化为速度，位移时程曲线

采用梯形积分法，可得到速度和位移时程曲线
%%%test%%%
fid = fopen(‘Acc1.AT2’); %%导入加速度数据；
Acc = textscan(fid,’%f’,‘headerlines’,4);
Acc = Acc{1,1};
fclose(fid);
dt = 0.01; %采样间距
t = 0:dt:(length(Acc)-1)dt;
vn1 = cumtrapz(t, Acc); %梯形积分法
vn = vn1 - repmat(mean(vn1), size(vn1,1), 1);
xn1 = cumtrapz(t, vn);
xn =xn1 - repmat(mean(xn1), size(xn1,1), 1);
figure
plot(t,9810vn,‘k’);
hold on
plot(t,9810*xn,‘r’); %单位cm
legend(‘velocity’,‘displacement’);

此时，会发现积分得到的位移初始值不为0，这和实际情况不符合；
要得到初始值为0的位移，只需要注释掉xn =xn1 - repmat(mean(xn1), size(xn1,1), 1);这一行代码即可
%%%test%%%
fid = fopen(‘Acc1.AT2’); %%导入加速度数据；
Acc = textscan(fid,’%f’,‘headerlines’,4);
Acc = Acc{1,1};
fclose(fid);
dt = 0.01; %采样间距
t = 0:dt:(length(Acc)-1)dt;
vn1 = cumtrapz(t, Acc); %梯形积分法
vn = vn1 - repmat(mean(vn1), size(vn1,1), 1);
xn1 = cumtrapz(t, vn);
%xn =xn1 - repmat(mean(xn1), size(xn1,1), 1);
figure
plot(t,9810vn,‘k’);
hold on
plot(t,9810*xn1,‘r’); %单位cm
legend(‘velocity’,‘displacement’);

运行，即可得到位移和速度曲线：

时域积分：地震加速度时程曲线转化为速度，位移时程曲线相关推荐

ATVD地震加速度转化为速度和位移的小程序
ATVD地震加速度转化为速度和位移的小程序求解地震波加速度通过求导得出速度位移的过程: 假设地震波加速度为 a ( t ) a(t) a(t),可通过求解其积分得到速度和位移.首先对加速度进行一次积 ...
光滑曲线_对第一/二型曲线/曲面积分的小总结
公式第一型曲线积分(Line Integrals): 第二型曲线积分(Line Integrals of Vector Fields): 第一型曲面积分(Surface Integrals): 第二 ...
使用python编写多普勒频移函数，绘制多普勒频移随速度变化的曲线，给出代码并举例，代码以markdown格式给出...
下面是一段使用Python编写多普勒频移函数,绘制多普勒频移随速度变化的曲线的代码:python def DopplerShift(v): # v为物体的速度 f = np.arange(0, 100 ...
基于C语言实现离散时域积分算法
时域积分公式如下: 由于实际应用时传感器采样得到的点都是离散的,因此我们常常需要用到的是离散时域积分算法,对离散点的积分就是求和. matlab实现积分算法: % 积分 clc; t=0:1/100 ...
css贝塞尔曲线多个点_贝塞尔曲线实践
贝塞尔曲线: 贝塞尔曲线本质上是由线段和节点组成的,形象的说节点是可拖动的支点,线段像可伸缩的皮筋.一个常规的曲线往往由4个控制点构成(p0,p1,p2,p3),曲线经过起点(p0)和终点(p1). ...
python科赫曲线画树叶_Python——科赫曲线绘制
基本概念科赫曲线:科赫曲线是一种像雪花的几何曲线,所以又称为雪花曲线,它是de Rham曲线的特例.科赫曲线是出现在海里格·冯·科赫的论文<关于一条连续而无切线,可由初等几何构作的曲线> ...
android贝塞尔曲线,一文解析 Android 贝塞尔曲线
原标题:一文解析 Android 贝塞尔曲线相信很多同学都知道"贝塞尔曲线"这个词,我们在很多地方都能经常看到.利用"贝塞尔曲线"可以做出很多好看的UI效果, ...
java 贝塞尔曲线_贝塞尔曲线：原理、自定义贝塞尔曲线View、使用！！！
一.原理转自:http://www.2cto.com/kf/201401/275838.html Android动画学习Demo(3) 沿着贝塞尔曲线移动的Property Animation Pr ...
计算机图形学曲线生成原理,计算机图形学_曲线及生成.ppt
计算机图形学_曲线及生成华中理工大学计算机学院陆枫 99-7 1999年7月 7.2.1 曲线的表示要求 1)唯一性 2)几何不变性 3)易于定界 4)统一性 5)易于实现光滑连接 6)几何直观 ...

时域积分：地震加速度时程曲线转化为速度，位移时程曲线

时域积分：地震加速度时程曲线转化为速度，位移时程曲线相关推荐

最新文章

热门文章