ATVD地震加速度转化为速度和位移的小程序

求解地震波加速度通过求导得出速度位移的过程:
假设地震波加速度为 a ( t ) a(t) a(t),可通过求解其积分得到速度和位移。首先对加速度进行一次积分得到速度:
v ( t ) = ∫ 0 t a ( τ ) d τ v(t) = \int_0^t a(\tau) d\tau v(t)=∫0t​a(τ)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)=∫0t​v(τ)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​∫0t​a(s)dsdt
该公式将速度的积分替换为了加速度的二重积分。在实际计算中,可以采用数值积分方法,比如梯形法则或者辛普森法则,来计算上式中的积分。








目前对于一些存在基线漂移的情况,无法进行校正。

ATVD地震加速度转化为速度和位移的小程序相关推荐

  1. 时域积分:地震加速度时程曲线转化为速度,位移时程曲线

    采用梯形积分法,可得到速度和位移时程曲线 %%%test%%% fid = fopen('Acc1.AT2'); %%导入加速度数据: Acc = textscan(fid,'%f','headerl ...

  2. mpvue开发小程序所遇问题及h5转化方案

    项目结构 |---build|---pages.js文件目录|---src|---component子组件|---pages|---业务页面|---store,vuex储存|---utils|---请 ...

  3. 如何利用 webp 进行小程序图片加载速度的优化

    导语 最近很长一段时间没有更新博客,一方面是自己最近参与了小程序的开发,另一方面也是自己略有些怠惰,给自己记个过~那么现在既然回到学校那么还是要分享一些知识的. 前一阵子参与微信小程序开发时遇到了一个 ...

  4. matlab加速度转化为位移,科学网—Matlab编写由加速度积分得到速度和位移函数 - 王德才的博文...

    研究中经常会用到将加速度数据积分成速度和位移,对前几年编写的程序重新进行了改写.为了测试和验证peer2acc函数.dirFolder函数和acc2vd函数,编写了一个简单的测试程序 acc2vd函数 ...

  5. matlab加速度转化为位移,matlab数值积分实现加速度、速度、位移的转换(时域频域积分)...

    对此可以用滤波的方法将大的趋势项去掉. 测试的代码如下 % 测试积分对正弦信号的作用 clc clear close all %% 原始正弦信号 ts = 0.001; fs = 1/ts; t = ...

  6. 51单片机基本刷屏测试实验_基于单片机的发动机振动速度、位移和加速度测量方法...

    Single-chip microcomputer-based measuring of engine vibration speed.displacement and acceleration Ab ...

  7. matlab绘制加速度曲线,加速度积分得到速度和位移的问题小结

    在振动信号采集分析中,我们常常使用加速度传感器,是因为加速度传感器具有可靠性高,可测频带宽,结构小巧,抗干扰能力强等优点.但是有些时候,必需要得到速度和位移信号,这个时候就出现了如何通过对加速度积分得 ...

  8. matlab教程lqg函数与收敛速度,科学网—Matlab编写由加速度积分得到速度和位移函数 - 王德才的博文...

    研究中经常会用到将加速度数据积分成速度和位移,对前几年编写的程序重新进行了改写.为了测试和验证peer2acc函数.dirFolder函数和acc2vd函数,编写了一个简单的测试程序 acc2vd函数 ...

  9. Python在振动信号处理中的应用(五):振动加速度信号转换为速度或位移信号

    文章目录 一.概述 二.算法原理 2.1 时域积分 2.1 频域处理 三.Python中相关函数 3.1 scipy.interger介绍 3.2 scipy.fftpack介绍 四.python代码 ...

最新文章

  1. 【不同的Java垃圾回收器的比较】
  2. nginx进程管理之worker进程
  3. 代码夹带是洪水猛兽吗?
  4. 计算机组成原理实验load,计算机组成原理实验报告五
  5. git bash卡顿解决
  6. Python 网站应用从开发到部署
  7. CentOS7恢复rm -rf 误删的xfs系统
  8. 英雄无敌3 Def 文件对应生物列表
  9. 短信系统实现原理(基于redis)
  10. Excel表格无法重命名工作表
  11. 办公室计算机联机,两台电脑怎样联机?
  12. OSG使用OpenGL(以及glad库)绘制自定义图形
  13. android 批量安装apk,如何一键批量安装APP应用apk安装包?
  14. word常见问题_2 添加不同的页码
  15. 【matlab教程】02、拼接矩阵或向量
  16. 跌宕七十年,日本制造业兴衰「启示录」
  17. 成为一名专业Web渗透工程师难吗?
  18. Mark-sweep GC
  19. form layui 同时提交多个对象_layui实现form表单同时提交数据和文件的代码
  20. 校园邮箱免费激活idea等Jetbrains系列旗舰版产品

热门文章

  1. oracle中00604,【案例】Oracle无法启动报错ORA-00604 ORA-01578官方解决办法
  2. 上位机与PLC 通讯源码 上位机与三菱PLC,西门子PLC通讯 同时一起通讯,单独控制,三菱采用官方MX 通讯,支持三菱FX系列,
  3. Base64的编码和解码
  4. centos7分区、挂载、磁盘合并
  5. XSS专栏之常见xss--总结备忘
  6. CRT 开启日志记录功能
  7. sipp: bind_local;watchdog timer trip
  8. 超经典的20道SQL题目(附加解题思路)
  9. 教你如何用手机连接电脑上网
  10. memoized函数