ATVD地震加速度转化为速度和位移的小程序
ATVD地震加速度转化为速度和位移的小程序
求解地震波加速度通过求导得出速度位移的过程:
假设地震波加速度为 a ( t ) a(t) a(t),可通过求解其积分得到速度和位移。首先对加速度进行一次积分得到速度:
v ( t ) = ∫ 0 t a ( τ ) d τ v(t) = \int_0^t a(\tau) d\tau v(t)=∫0ta(τ)dτ
其中, t t t 是当前时间, τ \tau τ 是积分变量。该式表示在 0 0 0 时刻的速度为 0 0 0,在 t t t 时刻的速度为积分后的结果。然后,可对速度再次积分得到位移:
x ( t ) = ∫ 0 t v ( τ ) d τ = ∫ 0 t ∫ 0 τ a ( s ) d s d τ x(t) = \int_0^t v(\tau) d\tau = \int_0^t \int_0^\tau a(s) ds d\tau x(t)=∫0tv(τ)dτ=∫0t∫0τa(s)dsdτ
这个式子表示在 0 0 0 时刻的位移也是 0 0 0,在 t t t 时刻的位移为速度的积分值。注意,对速度积分的上限是 τ \tau τ,因为速度是加速度的积分,而加速度本身是一个函数,需要在积分中引入一个变量 τ \tau τ,作为积分的上限。
上述公式可以进一步简化为:
x ( t ) = 1 2 ∫ 0 t ∫ 0 t a ( s ) d s d t x(t) = \frac{1}{2} \int_0^t \int_0^t a(s) ds dt x(t)=21∫0t∫0ta(s)dsdt
该公式将速度的积分替换为了加速度的二重积分。在实际计算中,可以采用数值积分方法,比如梯形法则或者辛普森法则,来计算上式中的积分。
目前对于一些存在基线漂移的情况,无法进行校正。
ATVD地震加速度转化为速度和位移的小程序相关推荐
- 时域积分:地震加速度时程曲线转化为速度,位移时程曲线
采用梯形积分法,可得到速度和位移时程曲线 %%%test%%% fid = fopen('Acc1.AT2'); %%导入加速度数据: Acc = textscan(fid,'%f','headerl ...
- mpvue开发小程序所遇问题及h5转化方案
项目结构 |---build|---pages.js文件目录|---src|---component子组件|---pages|---业务页面|---store,vuex储存|---utils|---请 ...
- 如何利用 webp 进行小程序图片加载速度的优化
导语 最近很长一段时间没有更新博客,一方面是自己最近参与了小程序的开发,另一方面也是自己略有些怠惰,给自己记个过~那么现在既然回到学校那么还是要分享一些知识的. 前一阵子参与微信小程序开发时遇到了一个 ...
- matlab加速度转化为位移,科学网—Matlab编写由加速度积分得到速度和位移函数 - 王德才的博文...
研究中经常会用到将加速度数据积分成速度和位移,对前几年编写的程序重新进行了改写.为了测试和验证peer2acc函数.dirFolder函数和acc2vd函数,编写了一个简单的测试程序 acc2vd函数 ...
- matlab加速度转化为位移,matlab数值积分实现加速度、速度、位移的转换(时域频域积分)...
对此可以用滤波的方法将大的趋势项去掉. 测试的代码如下 % 测试积分对正弦信号的作用 clc clear close all %% 原始正弦信号 ts = 0.001; fs = 1/ts; t = ...
- 51单片机基本刷屏测试实验_基于单片机的发动机振动速度、位移和加速度测量方法...
Single-chip microcomputer-based measuring of engine vibration speed.displacement and acceleration Ab ...
- matlab绘制加速度曲线,加速度积分得到速度和位移的问题小结
在振动信号采集分析中,我们常常使用加速度传感器,是因为加速度传感器具有可靠性高,可测频带宽,结构小巧,抗干扰能力强等优点.但是有些时候,必需要得到速度和位移信号,这个时候就出现了如何通过对加速度积分得 ...
- matlab教程lqg函数与收敛速度,科学网—Matlab编写由加速度积分得到速度和位移函数 - 王德才的博文...
研究中经常会用到将加速度数据积分成速度和位移,对前几年编写的程序重新进行了改写.为了测试和验证peer2acc函数.dirFolder函数和acc2vd函数,编写了一个简单的测试程序 acc2vd函数 ...
- Python在振动信号处理中的应用(五):振动加速度信号转换为速度或位移信号
文章目录 一.概述 二.算法原理 2.1 时域积分 2.1 频域处理 三.Python中相关函数 3.1 scipy.interger介绍 3.2 scipy.fftpack介绍 四.python代码 ...
最新文章
- 【不同的Java垃圾回收器的比较】
- nginx进程管理之worker进程
- 代码夹带是洪水猛兽吗?
- 计算机组成原理实验load,计算机组成原理实验报告五
- git bash卡顿解决
- Python 网站应用从开发到部署
- CentOS7恢复rm -rf 误删的xfs系统
- 英雄无敌3 Def 文件对应生物列表
- 短信系统实现原理(基于redis)
- Excel表格无法重命名工作表
- 办公室计算机联机,两台电脑怎样联机?
- OSG使用OpenGL(以及glad库)绘制自定义图形
- android 批量安装apk,如何一键批量安装APP应用apk安装包?
- word常见问题_2 添加不同的页码
- 【matlab教程】02、拼接矩阵或向量
- 跌宕七十年,日本制造业兴衰「启示录」
- 成为一名专业Web渗透工程师难吗?
- Mark-sweep GC
- form layui 同时提交多个对象_layui实现form表单同时提交数据和文件的代码
- 校园邮箱免费激活idea等Jetbrains系列旗舰版产品