====== Differences ====== This shows you the differences between two versions of the page.
Next revision | Previous revision | ||
gibson:teaching:spring-2016:math445:lecture:pendulum [2016/04/27 13:15] gibson created |
gibson:teaching:spring-2016:math445:lecture:pendulum [2016/04/27 13:24] (current) gibson [Changing a second-order ODE into a system of first-order ODEs] |
||
---|---|---|---|
Line 1: | Line 1: | ||
====== The plane pendulum ====== | ====== The plane pendulum ====== | ||
+ | ==== Mathematical models==== | ||
- | n class we derived four different mathematical models of the plane pendulum, listed here in order of decreasing mathematical complexity. | + | In class we derived four different mathematical models of the plane pendulum, listed here in order of decreasing mathematical complexity. |
** Nonlinear damped pendulum.** This is the most physically realistic model. It includes a linear model of air resistance the $\alpha/m \; d\theta/dt$ term and is accurate for large displacement angles $\theta$. | ** Nonlinear damped pendulum.** This is the most physically realistic model. It includes a linear model of air resistance the $\alpha/m \; d\theta/dt$ term and is accurate for large displacement angles $\theta$. | ||
Line 35: | Line 36: | ||
In class we showed that the linear undamped pendulum has a solution of the form $\theta(t) = \theta_0 \cos \omega t$ where $\omega = \sqrt{g/\ell}$ and $\theta_0$ is an arbitrary initial angular displacement. (Note that $\theta_0$ must be small for the small-angle approximation to be valid!) | In class we showed that the linear undamped pendulum has a solution of the form $\theta(t) = \theta_0 \cos \omega t$ where $\omega = \sqrt{g/\ell}$ and $\theta_0$ is an arbitrary initial angular displacement. (Note that $\theta_0$ must be small for the small-angle approximation to be valid!) | ||
- | We also showed that the 2nd-order equation in the scalar variable $\theta(t)$ can be transformed into a 1st order equation in the vector variable $\vec{x}(t)$, using the substitution | + | |
+ | |||
+ | ---- | ||
+ | |||
+ | ==== Changing a second-order ODE into a system of first-order ODEs ==== | ||
+ | |||
+ | A 2nd-order ordinary differential equation in the scalar variable $\theta(t)$ can be transformed into a 1st order equation in the vector variable $\vec{x}(t)$, using the substitution | ||
\begin{eqnarray*} | \begin{eqnarray*} | ||
Line 41: | Line 48: | ||
\end{eqnarray*} | \end{eqnarray*} | ||
- | If you differentiate the above equation in time and perform a few substitutions between $\theta, d\theta/dt$ and $x_1, x_2$, you get | + | For example, to transform the linear undamped pendulum equation $d^2\theta/dt^2 + (g/\ell) \theta = 0$, |
+ | let $x_1 = \theta$ and $x_2 = d\theta/dt$. Differentiate those two equations to get | ||
+ | $d x_1/dt = d\theta/dt = x_2$ and $d x_2/dt = d^2\theta/dt^2$. Now note that, according to the linear undamped pendulum equation, $d^2\theta/dt^2 = -(g/\ell) \theta = -(g/\ell) x_1$. Putting all this together, we can write | ||
\begin{eqnarray*} | \begin{eqnarray*} |