Using Euler-Cromer to find the motion of a simple pendulum?

Question

Rachel 2012-4-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/34444-using-euler-cromer-to-find-the-motion-of-a-simple-pendulum

My professor wants us to find the number of periods a simple, meter-long pendulum will have in a day if given an initial angle of 50 degrees using the Euler-Cromer method. My problem is that I can't come up with al the equations I need.

I know that I need: dθ/dt = f(t, ω) and dω/dt = g(t, θ).

I managed to get dω/dt = -9.8*sin(θ). I can't figure out what dθ/dt should be, though. Once I have that, I understand how to use for loops to find θ and ω for time steps and if statements to fin the number of periods throughout the day (even if it is inefficient), I just can't seem to get started. Can anybody give me some help as to what dθ/dt is? Remember, the starting angle is 50 degrees, so the small angle approximation cannot be used.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt Kindig 2012-4-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/34444-using-euler-cromer-to-find-the-motion-of-a-simple-pendulum#answer_43193

Can you use the ODE solvers built-in to Matlab (something like 'ode45')? If so, that would be the easiest way, and would avoid looping completely.

If the objective of the assignment is to implement your own ODE solver, then I would look up the definition for Runge-Kutta integration on Wikipedia or similar.

By the way, d(theta)/dt is just w (omega).

2 个评论
显示无隐藏无

Rachel 2012-4-4

We are supposed to implement our own ODE solver, and I know about the Runge-Kutta method, but we are told explicitly in the problem to use the Euler-Cromer.

I tried using dθ/dt(i) = ω(i), but I end up with a pendulum that swings in circles, not an arc of 100 degrees total.

Matt Kindig 2012-4-4