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:03] gibson [Problem 3: nonlinear damped 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 time series and phase portrait 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 | ||
| 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. | ||