Esta simulação funciona com base na solução da EDO
$$
\ddot{\theta} + \gamma\dot{\theta} + \dfrac{g}{L}\sin\theta = 0
$$
A EDO é resolvida pelo método de Runge-Kutta de quarta ordem.
A aceleração da gravidade é assumida como $g = 9.8$ m/s$^2$ e o fator de amortecimento como $\gamma = 0.05$ s$^{-1}$.
Há, além disso, a inclusão de um erro $\alpha$ gerado por uma variável aleatória.
Isto é, a cada vez que o meio período do pêndulo é calculado, ou a sua velocidade, um valor $\alpha$
é sorteado dentro do intervalo $[-\epsilon, \epsilon]$
e somado ao meio período ou a velocidade informada no dispositivo.
$\epsilon$ é a amplitude do ruído. O valor de $\alpha$ é composto de duas partes:
erro aleatório (~ distribuição gaussiana em torno do valor calculado)
e erro sistemático (offset somado ao valor calculado).
Leia mais: