这个问题可能是由于代码中计算周期的方式不正确引起的。以下是一种正确的代码示例:
import math
# 定义常量
G = 9.8 # 重力加速度
L = 1.0 # 摆长
# 初始化
theta = math.pi / 4
omega = 0.0
t = 0.0
dt = 0.01 # 时间步长
# 计算运动
while t < 10:
alpha = (-G / L) * math.sin(theta)
omega += alpha * dt
theta += omega * dt
# 将角度限制到 [-pi, pi] 之间
while theta > math.pi:
theta -= 2 * math.pi
while theta < -math.pi:
theta += 2 * math.pi
t += dt
# 输出结果
print(theta)
在这个代码示例中,摆锤的运动是按照欧拉方法进行计算的。使用此方法时,我们必须使用足够小的时间步长以确保计算的准确性,同时还要在每次计算后将角度限制在 [-π, π] 之间,以避免出现角度偏移和计算错误。