<पी>
यह सिमुलेशन ईडीओ समाधान के आधार पर काम करता है
$$
\ddot{\theta} + \gamma\dot{\theta} + \dfrac{g}{L}\sin\theta = 0
$$
<पी>
ODE को चौथे क्रम रंज-कुट्टा विधि द्वारा हल किया जाता है।
गुरुत्वाकर्षण के कारण त्वरण $g = 9.8$ m/s$^2$ और अवमंदन कारक $\गामा = 0.05$ s$^{-1}$ माना जाता है।
<पी>
इसके अलावा, एक यादृच्छिक चर द्वारा उत्पन्न त्रुटि $\alpha$ का समावेश भी है।
अर्थात्, हर बार पेंडुलम की आधी अवधि, या उसकी गति की गणना करने पर, एक मान $\alpha$ होता है
अंतराल $[-\epsilon, \epsilon]$ के भीतर खींचा जाता है
और आधे समय या डिवाइस पर रिपोर्ट की गई गति में जोड़ा गया।
$\epsilon$ शोर का आयाम है। $\alpha$ का मान दो भागों से बना है:
यादृच्छिक त्रुटि (~परिकलित मूल्य के आसपास गाऊसी वितरण)
और व्यवस्थित त्रुटि (गणना मूल्य में ऑफसेट जोड़ा गया)।
<पी>
और पढ़ें:
-
S. Thornton, J. Marion
Dinâmica clássica de partículas e sistemas
(Cengage Learning, 2011)