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