Cette simulation fonctionne sur la base de la solution EDO
$$
\ddot{\theta} + \gamma\dot{\theta} + \dfrac{g}{L}\sin\theta = 0
$$
L'ODE est résolue par la méthode Runge-Kutta du quatrième ordre.
L'accélération due à la gravité est supposée être $g = 9,8$ m/s$^2$ et le facteur d'amortissement comme $\gamma = 0,05$ s$^{-1}$.
De plus, il y a l'inclusion d'une erreur $\alpha$ générée par une variable aléatoire.
C'est-à-dire qu'à chaque fois que l'on calcule la demi-période du pendule, ou sa vitesse, une valeur $\alpha$
est dessiné dans l'intervalle $[-\epsilon, \epsilon]$
et ajouté à la mi-temps ou à la vitesse indiquée sur l'appareil.
$\epsilon$ est l'amplitude du bruit. La valeur de $\alpha$ est composée de deux parties :
erreur aléatoire (~distribution gaussienne autour de la valeur calculée)
et erreur systématique (offset ajouté à la valeur calculée).
En savoir plus: