【数值积分(欧拉积分、Velet积分与蛙跳积分)】

本文介绍了数值积分的基本概念,包括泰勒级数展开与积分的关系,以及常用的欧拉积分、Velet积分和蛙跳积分。Velet积分通过分离速度和位置的计算提高精度,而蛙跳积分则是对Velet积分的优化,通过计算半时间步长的速度来提升效果。这些方法在实际模拟中对于提高精度有着重要作用。
摘要由CSDN通过智能技术生成

数值积分(欧拉积分、Velet积分与蛙跳积分)

数值积分

我们使用的各种数值积分方法都是对积分做近似,近似的本质类似于泰勒级数展开式。在实际模拟中,不可能使用泰勒级数展开中的无穷项,所以积分方法只用级数的截断类似。更健壮的积分方法是使用级数中的更多项。积分方法对模拟的精度有非常大的影响。最简单的积分有欧拉积分,其次有Velet积分(类似蛙跳积分)、蛙跳积分。理解这些概念首先需了解级数展开与积分。

级数展开与积分

泰勒级数展开
s ( t + h ) = s ( t ) + s ˙ ( t ) h + 1 2 s ¨ ( t ) h 2 + ⋯ \mathbf{s}(t+h)=\mathbf{s}(t)+\dot{\mathbf{s}}(t) h+\frac{1}{2} \ddot{\mathbf{s}}(t) h^{2}+\cdots s(t+h)=s(t)+s˙(t)h+21s¨(t)h2+
欧拉积分只用到了前两项: s ( t + h ) ≈ s ( t ) + s ˙ ( t ) h \mathbf{s}(t+h) \approx \mathbf{s}(t)+\dot{\mathbf{s}}(t) h s(t+h)s(t)+s˙(t)h。因为只用到了一阶导数,所以欧拉积分是一阶方法。
常见的二阶方法有
梯形法 s ( t + h ) = s ( t ) + s ˙ ( t ) + s ˙ ( t + h ) 2 h \mathbf{s}(t+h)=\mathbf{s}(t)+\frac{\dot{\mathbf{s}}(t)+\dot{\mathbf{s}}(t+h)}{2} h s(t+h)=s(t)+2s˙(t)+s˙(t+h)h
中点法 s ( t + h ) = s ( t ) + s ˙ ( t + 0.5 h ) h \mathbf{s}(t+h)=\mathbf{s}(t)+\dot{\mathbf{s}}(t+0.5 h) h s(t+h)=s(t)+s˙(t+0.5h)h

Velet积分

Velet积分假设加速度仅由位置决定,即 a = 1 m f ( x

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值