声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1175|回复: 6

[编程技巧] 数值积分

[复制链接]
发表于 2010-5-24 14:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我有一个随时间变量的加速度信号,现想通过数值积分获得速度信号,该怎么积分呢?
大虾给以指点么
附上加速度数值

timenofirst.txt

113.83 KB, 下载次数: 3

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2010-5-24 15:08 | 显示全部楼层

  1. doc trapz
复制代码

评分

1

查看全部评分

 楼主| 发表于 2010-5-24 23:29 | 显示全部楼层

回复 沙发 rocwoods 的帖子

这组加速度信号积分出的速度应该是一组数,trapz积分处了一个数字,是不是先差值拟合出加速度函数,然后在积分?

>> z=trapz(t,x)

z =

-2.3804e-006

>>
发表于 2010-5-25 07:44 | 显示全部楼层
LZ请help trapz仔细些, 最底下一般不是还有相关函数吗?:lol
或许cumtrapz就是LZ要的!?
个人记得已有许多诸如此类的讨论, 搜索下并看看!:@)
 楼主| 发表于 2010-5-25 13:01 | 显示全部楼层

回复 地板 ChaChing 的帖子

,先去除趋势项后,使用v=cumtrapz(t,a)对上述加度积分得到了速度,但是进一步对速度使用s=cumtrapz(t,v),出现明显和实际位移不符合,能不能直接对a积分两次得到位移?

[ 本帖最后由 dreamstone 于 2010-5-25 19:11 编辑 ]
发表于 2010-6-3 11:25 | 显示全部楼层
数值积分的话,还可以用quad,quadl,quadgk等函数
 楼主| 发表于 2010-6-3 18:35 | 显示全部楼层

回复 6楼 qibbxxt 的帖子

呵呵,谢谢了,我自己通过滤波,去趋势项,自己编了个梯形公式,不过误差还有点大快10%了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 03:29 , Processed in 0.066950 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表