Euler-Cromer方法实现正确性求证：与讲师观点分歧咨询

代码实现 #

k=3;m=123
eulerC_x=[A]
eulerC_v=[0]
eulerC_t=[0]
c=0
for z in np.arange(t_euler,t_max,dt_euler):
    a=-(k/m)*eulerC_x[c]
    eulerC_v.append(eulerC_v[c]+a*dt_euler)
    eulerC_x.append(eulerC_x[c]+eulerC_v[c+1]*dt_euler) 
    #c+1 uses new velocity value to update the position, removing the energy increase in euler         
    eulerC_t.append(z)
    c+=1

问题背景 #

上述代码用于求解简谐运动的常微分方程（ODE）。我的实现逻辑是先更新速度，再使用更新后的(n+1)阶速度计算新位置，认为这是正确的Euler-Cromer方法——该方法可消除普通Euler方法的能量增长问题，且我已通过图表验证了这一效果。但讲师认为应先更新位置，再用(n+1)阶位置计算新速度，判定我的实现错误。多方资料支持我的观点，特此求证该代码实现是否正确。

解答 #

你的代码实现是正确的，完全符合标准Euler-Cromer方法的定义：

• 对于简谐运动的二阶ODE，通常拆分为两个一阶方程：$\frac{dx}{dt}=v$，$\frac{dv}{dt}=-\frac{k}{m}x$
• Euler-Cromer方法的核心改进就是用更新后的速度计算新位置，标准步骤为：
  1. 基于当前位置$x_n$计算加速度$a_n=-\frac{k}{m}x_n$
  2. 更新速度：$v_{n+1}=v_n + a_n \cdot dt$
  3. 用新速度更新位置：$x_{n+1}=x_n + v_{n+1} \cdot dt$
• 这种顺序的更新能让系统能量近似守恒，完美解决普通Euler方法的能量持续增长问题，你的图表验证结果也佐证了这一点。

讲师提到的“先更新位置再用新位置算速度”属于半隐式Euler的另一种变体，但并非标准的Euler-Cromer实现，且对于简谐运动这类问题，你的实现才是被广泛认可的、能有效抑制能量误差的版本。

内容的提问来源于stack exchange，提问作者morgy2190