Numerical instability of spherical pendulum
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Bas Siebers
el 25 de Feb. de 2015
Comentada: Mischa Kim
el 26 de Feb. de 2015
Hi,
I am trying to simulate a spherical pendulum. The equation of motion of the spherical pendulum are:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/147306/image.png)
So far, I was able to simulate the equation of motion with a ode45 solver. However I experiencing numerical instabilities when the phi angle approach zero.
Does any one have an idea to get rid off these numerical instabilities?
Thank you in advance,
Bas
0 comentarios
Respuesta aceptada
Mischa Kim
el 25 de Feb. de 2015
Bas, the plus sign in your equation does look a bit strange. Shouldn't that be a minus instead?
8 comentarios
Mischa Kim
el 26 de Feb. de 2015
Angular momentum is conserved (constant), not azimuth rate. The equation for the angular momentum is
m*l^2*sin(phi)^2*thetadot = const.
Therefore, when phi decreases, thetadot increases, just like described above.
Más respuestas (0)
Ver también
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!