User Tools

Site Tools


gibson:teaching:spring-2016:math445:lab11

====== Differences ====== This shows you the differences between two versions of the page.

Link to this comparison view

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.
gibson/teaching/spring-2016/math445/lab11.1461258184.txt.gz · Last modified: 2016/04/21 10:03 by gibson