二维平面内点的匀加速运动：Python代码计算与位移公式结果不符

嘿，我来帮你拆解这个问题～

你的代码逻辑和匀加速直线运动的物理模型有个关键偏差：

• 你用循环直接累加当前末速度来计算位移，相当于把每个时间单位（循环的1次迭代）的位移当成了 v_t * Δt（Δt=1）；但匀加速运动中，每个时间段的正确位移应该是平均速度乘以时间，也就是 (v_初 + v_末)/2 * Δt。
• 你的循环执行了10次，每次的vx依次是10*1、10*2...10*10，总位移就是10*(1+2+...+10) = 10*(10*11)/2 = 550——这就是你得到550的原因，而公式½*a*t²是连续时间下的积分结果，和离散循环的错误累加逻辑自然不一样。

方法1：直接用物理公式计算（最简洁） #

既然是匀加速直线运动，直接代入公式即可，完全不需要循环。如果不想用单独的时间变量，直接把数值写进公式里就行：

import math

x = 0.0
y = 0.0
acceleration = 10
angle = 0

# 直接计算总位移并分解到x、y轴
x += 0.5 * acceleration * 10 ** 2 * math.cos(math.radians(angle))
y += -0.5 * acceleration * 10 ** 2 * math.sin(math.radians(angle))

print(x, y)  # 输出：500.0 0.0

方法2：修正循环的位移计算逻辑（如果要模拟运动过程） #

如果一定要用循环来模拟运动，那需要用平均速度计算每个时间段的位移，同时可以用while循环替代time变量来控制迭代次数：

import math

x = 0.0
y = 0.0
acceleration = 10
angle = 0
vx = 0.0
vy = 0.0
current_iteration = 0
max_iterations = 10  # 这里可以换成你需要的终止条件，比如位移阈值

while current_iteration < max_iterations:
    # 记录当前时间段的初速度
    prev_vx, prev_vy = vx, vy
    # 更新末速度
    vx += acceleration * math.cos(math.radians(angle))
    vy += -acceleration * math.sin(math.radians(angle))
    # 用平均速度计算这段时间的位移
    x += (prev_vx + vx) / 2
    y += (prev_vy + vy) / 2
    current_iteration += 1

print(x, y)  # 输出：500.0 0.0

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