单摆的数值模拟可以使用python的scipy库中的odeint函数来实现。首先需要确定微分方程,然后使用odeint函数求解。示例代码如下:
from scipy.integrate import odeint
import numpy as np
def pend(y, t, b, c):
theta, omega = y
dydt = [omega, -b*omega - c*np.sin(theta)]
return dydt
b = 0.25
c = 5.0
y0 = [np.pi - 0.1, 0.0]
t = np.linspace