Verlet Integration
Verlet 积分法是一种用于求解牛顿运动方程的数值方法,被广泛运用于动力学模拟以及视频游戏中。尔莱算法的优点在于:数值稳定性比简单的欧拉方法高很多,并保持了物理系统中的时间可逆性与相空间体积元体积守恒的性质。
基本韦尔莱算法
根据牛顿运动方程有
![](https://i-blog.csdnimg.cn/blog_migrate/5542e18ffdf80f21404421b015fce083.png)
![](https://i-blog.csdnimg.cn/blog_migrate/2dc9fc57d0773d3efbd370241a54ebf2.png)
![](https://i-blog.csdnimg.cn/blog_migrate/3d8097432bd6c9717184567082f730b8.png)
![](https://i-blog.csdnimg.cn/blog_migrate/e4ab26c9b7d6b16788baf90c987fb3eb.png)
![](https://i-blog.csdnimg.cn/blog_migrate/d0430899c461863722ab1d45e8b98581.png)
![](https://i-blog.csdnimg.cn/blog_migrate/b2121dc182bc710eafeda7be1848398d.png)
![](https://i-blog.csdnimg.cn/blog_migrate/c3848cf2d102dd7608becf7567651ad8.png)
![](https://i-blog.csdnimg.cn/blog_migrate/cf07a191c8c08fd2062399bf735add5e.png)
速度表示的韦尔莱算法
一般地,速度表示的韦尔莱算法更为常用,它可以给出同一时间变量下的速度和位置。它实际上与基本韦尔莱算法等价,精度相同。
首先对位置进行泰勒展开
![](https://i-blog.csdnimg.cn/blog_migrate/5b0a047baba982e65374d5b3e5cd587b.png)
![](https://i-blog.csdnimg.cn/blog_migrate/801d076192dae26bf5b091905aab6c11.png)
对两式相减可得
![](https://i-blog.csdnimg.cn/blog_migrate/af456620f9d4b5d8fc3a45704bf44a06.png)
![](https://i-blog.csdnimg.cn/blog_migrate/69038994a1cdde920f7b341e9f46d6ca.png)
![](https://i-blog.csdnimg.cn/blog_migrate/afa35d055696d1547044a8f458a448a9.png)
![](https://i-blog.csdnimg.cn/blog_migrate/f36c1ae482981a6c1fc3ab4977e89e92.png)
![](https://i-blog.csdnimg.cn/blog_migrate/1d805dff31e48adfff0296423288eae8.png)
- 首先通过泰勒展开式
计算得到位置
- 由
和系统的相互作用势条件(如果相互作用仅依赖位置
)可以求得力场
- 由速度表示的韦尔莱公式求出新的速度