R中的不定积分

Question

我期待计算方程的不定积分.

我有来自加速计的数据通过一个可视化的C程序馈入R,并从那里它是足够简单拿出一个方程来表示加速曲线.这一切都很好,但我也需要计算冲击速度.从好醇"天高中我的理解,我的加速曲线的不定积分将产生的速度方程.

我知道这是很容易与执行数值积分integrate()功能,还有什么是可以与之相比的不定积分？

Answer 1

library(Ryacas)
x <- Sym("x")
Integrate(sin(x), x)

给

expression(-cos(x))

另一种方式:

yacas("Integrate(x)Sin(x)")

你可以在这里找到功能参考

Answer 2

如果您提到的 NA 在指示没有加速度输入的意义上提供了信息，那么它们应该被零替换。假设您的数据acc.vec和设备以以下速率记录rec_per_sec：

acc.vec[is.na(ac.vec)] <- 0
vel.vec <- cumsum(acc.vec)/recs_per_sec

我不认为构建最佳拟合曲线会提高这种情况下的准确性。绘制速度与时间的关系图：

plot(1:length(acc.vec)/recs_per_sec, vel.vec, 
       xlab="Seconds", ylab="Integrated Acceleration = Velocity")