====== Differences ====== This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
gibson:teaching:spring-2016:math445:lab11 [2016/04/21 10:02] gibson [Problem 3: nonlinear undamped pendulum] |
gibson:teaching:spring-2016:math445:lab11 [2016/04/26 12:42] (current) vining |
||
---|---|---|---|
Line 56: | Line 56: | ||
Compare the analytic solution $\theta(t) = \theta_0 \cos \omega t$ of the linear undamped pendulum to a numerical solution computed with Matlab's ''ode45'' function. Use constants | Compare the analytic solution $\theta(t) = \theta_0 \cos \omega t$ of the linear undamped pendulum to a numerical solution computed with Matlab's ''ode45'' function. Use constants | ||
- | * $g = 9./8$ (meters per second^2) | + | * $g = 9.8$ (meters per second^2) |
* $\ell = 1.0$ (meters) | * $\ell = 1.0$ (meters) | ||
* $\theta_0 = 0.1$ (radians) | * $\theta_0 = 0.1$ (radians) | ||
Line 117: | Line 117: | ||
---- | ---- | ||
- | ====Problem 3: nonlinear damped pendulum==== | + | ====Problem 4: nonlinear damped pendulum==== |
- | For this lab problem, you are to recreate the previous four plots for the nonlinear damped pendulum, whose equations of motion are | + | For this lab problem, you are to recreate the time series and phase portrait for the nonlinear damped pendulum, whose equations of motion are |
\begin{eqnarray*} | \begin{eqnarray*} | ||
Line 125: | Line 125: | ||
\end{eqnarray*} | \end{eqnarray*} | ||
- | Use parameter values $g=9.8, \ell=1$, $m=1$, and $\alpha = 1$. For the time series plots, initiate the pendulum with $x_1 = \theta=0$ and $x_2 = d\theta/dt = 2$. For the phase portraits, show the range $8 \leq \theta \leq 8$ on the horizontal axis and $-10 \leq d\theta/dt \leq 10$ on the vertical. On top of the quiver plots, show trajectories with initial conditions $\theta=0$ and a variety of $d\theta/dt$ ranging from -10 to 10 in steps of 1. | + | Use parameter values $g=9.8, \ell=1$, $m=1$, and $\alpha = 1$. For the time series plots, initiate the pendulum with $x_1 = \theta=0$ and $x_2 = d\theta/dt = 2$. For the phase portraits, show the range $-8 \leq \theta \leq 8$ on the horizontal axis and $-10 \leq d\theta/dt \leq 10$ on the vertical. On top of the quiver plots, show trajectories with initial conditions $\theta=0$ and a variety of $d\theta/dt$ ranging from -10 to 10 in steps of 1. |
+ | ---- | ||
Turn in your code, your plots, and answer the following questions | Turn in your code, your plots, and answer the following questions | ||
+ | ====Questions (to be answered at the end of your lab)==== | ||
**(a)** Describe the differences that you see in the phase portraits of the nonlinear pendulum compared to the linear pendulum. | **(a)** Describe the differences that you see in the phase portraits of the nonlinear pendulum compared to the linear pendulum. |