Un pendule élastique de fréquence propre fixée à 1 Hz est mis en oscillation forcée à l'aide d'un moteur à excentrique. La simulation permet de régler la fréquence f de l'excitation ainsi que le coefficient d'amortissement λ de l'oscillateur.
L'allongement \(x(t)\) du pendule élastique est modélisé via l'équation différentielle
\begin{equation} \ddot x+2\lambda\dot x+{\omega_0}^2 x={\omega_0}^2e(t) \end{equation}où \(e(t)\) est le déplacement de la liaison à l'excitateur, \(\omega_0\) est la pulsation propre du pendule et \(\lambda\) un coefficient d'amortissement pour tenir compte des frottements (que l'on suppose linéaires en vitesse).
Cette équation différentielle est intégrée numériquement via la méthode de Euler[1] avec un pas correspondant à 20ms.