动力学方程的数值解（动力学方程+行星运动轨道）

动力学方程的数值解（动力学方程+行星运动轨道）

---

Zhuofei(fregot@icloud.com) julia version:1.5.0

---

计算动力学方程

动力学方程的数值解

例如，动力学方程如下：

$-kx=m\frac{dv}{dt}$

如果假设 $k/m=1$, 则:

$-x=\frac{dv}{dt}$

假设初始状态为: $t=0$, $x=1,v=0$.

令 $\epsilon =0.01s$, 可以得到:

$x(t+\epsilon )=x(t)+\epsilon \ast v(t)$

$v(t+\epsilon )=v(t)+\epsilon \ast a(t)$

又因为:

$a(t)=\frac{dv}{dt}$