Esta simulación funciona en base a la solución EDO.
$$
\ddot{\theta} + \gamma\dot{\theta} + \dfrac{g}{L}\sin\theta = 0
$$
La EDO se resuelve mediante el método de Runge-Kutta de cuarto orden.
La aceleración debida a la gravedad se supone como $g = 9,8$ m/s$^2$ y el factor de amortiguación como $\gamma = 0,05$ s$^{-1}$.
Además, se incluye un error $\alpha$ generado por una variable aleatoria.
Es decir, cada vez que se calcula el medio periodo del péndulo, o su velocidad, se obtiene un valor $\alpha$
se dibuja dentro del intervalo $[-\epsilon, \epsilon]$
y se suma al tiempo medio o a la velocidad reportada en el dispositivo.
$\epsilon$ es la amplitud del ruido. El valor de $\alpha$ se compone de dos partes:
error aleatorio (~distribución gaussiana alrededor del valor calculado)
y error sistemático (compensación sumada al valor calculado).
Lea mas: