Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
ℹ Please use `linewidth` instead.
Chap 45 Exercises
\[ \newcommand{\dnorm}{\text{dnorm}} \newcommand{\pnorm}{\text{pnorm}} \newcommand{\recip}{\text{recip}} \]
Exercise 1 Using an angle from -90 to 90 degrees to specify directions, say in what directions each of the eigenvectors is pointed. (Your answer should be good to \(\pm 10 \text{deg}\). Remember, 0 degrees is horizontal.)
Exercise 2 The exercise is based on the interactive applet shown below. Click on the image to open up the applet in another tab of your browser. Then arrange so that the applet’s tab is side by side with this one. That way you can see the questions and figure out the answer at the same time.
We will look at trajectories that go round and round: oscillating solutions to the linear dynamics around a fixed point.
In the applet, click on the left-hand graph to set the \(a\) and \(b\) values in the [ab10] matrix. Clicking on the right-hand graph will start a trajectory (shown in red) from the place where you click.
You can create an oscillatory flow by clicking anywhere in the blue parabolic region in the left-hand graph. Play around with clicking in different spots in the parabolic region. Set the number of time steps for the red trajectory until you can see at least two full revolutions around the origin.
Click near (but not on) the fixed point. Which of these best describes the shape of the trajectory?
```{mcq}
#| label: oscil1-1
#| show_hints: true
1. A tightly wound spiral [ correct hint: You're right ]
2. A loosely wound spiral
3. A rectangular shape
4. Alternating between heading toward the fixed point and heading away from it.
```Leaving the a and b values the same as when you answered the previous question, start the trajectory about half-way from the fixed point. How does the new trajectory compare to the one from the previous problem?
```{mcq}
#| label: oscil1-2
#| show_hints: true
1. A more tightly wound spiral
2. A more loosely wound spiral [ correct hint: Excellent! ]
3. A spiral that goes the other way round
4. A rectangular shape
```Find some $a$ and $b$ that causes the trajectory to go round in a closed loop. (You might not be able to get it to close it exactly, but try the best you can.) Vary $a$ leaving $b$ alone, and vice versa, so you can see how the existence of the closed loop depends on $a$ and $b$. Which of these best describes the requirement for a closed loop?
```{mcq}
#| label: oscil1-3
#| show_hints: true
1. b is close to zero
2. a is close to zero [ correct hint: Correct ]
3. both a and b are close to zero
4. a is bigger than b
```Drawing on your experience from the previous question in vary ing $a$ and $b$, which of these best describes the requirement for a trajectory that spirals **out** from the origin?
```{mcq}
#| label: oscil1-4
#| show_hints: true
1. $a < 0$
2. $0 < a$ [ correct hint: Good job! ]
3. $b < 0$
4. $0 < b$
```Drawing on your experience from the previous question in varying $a$ and $b$, which of these best describes the requirement for a trajectory that spirals **in** toward the origin?
```{mcq}
#| label: oscil1-5
#| show_hints: true
1. $a < 0$ [ correct hint: Excellent! ]
2. $0 < a$
3. $b < 0$
4. $0 < b$
```Now you should be pretty good at navigating the a-b space to make differently-shaped round-and-round trajectories. Focus for the moment on whether the spirals head in to the fixed point (stable) or away from it (unstable). What is the relationship between the eigenvalues and whether the trajectory is stable or not?
```{mcq}
#| label: oscil1-6
#| show_hints: true
1. real part is negative for stable flow [ correct hint: You're right ]
2. real part is positive for stable flow
3. imaginary part is negative for stable flow
4. imaginary part is positive for stable flow
```Explore close to the edges of the parabolic zone in a-b space. What happens to the trajectories compared to when a-b is in the middle of the parabolic zone?
```{mcq}
#| label: oscil1-7
#| show_hints: true
1. the speed of rotation slows down [ correct hint: Yes ]
2. the speed of rotation is unchanged
3. the speed of rotation increases
```Look at the eigenvalues when you move around the a-b space. What aspect of the eigenvalues corresponds to the speed of revolution?
```{mcq}
#| label: oscil1-8
#| show_hints: true
1. the magnitude of the imaginary part [ correct hint: right-o ]
2. the magnitude of the real part
3. whether the imaginary part is positive or negative
4. whether the real part is positive or negative.
```Exercise 3
What are the numerical values of the eigenvalues of the matrix $\left[\begin{array}{rc}-1 & \frac{3}{4} \\ 1 & 0\end{array}\right]$?
```{mcq}
#| label: char1-1
#| show_hints: true
1. $\frac{1}{2}$ and $-\frac{3}{2}$ [ correct hint: Nice. ]
2. $\frac{1}{2}$ and $-\frac{5}{2}$
3. $\frac{1}{2}$ and $+\frac{5}{2}$
4. $-\frac{1}{2}$ and $-\frac{5}{2}$
```What are the numerical values of the eigenvalues of the matrix $\left[\begin{array}{rc}0 & \frac{1}{2} \\ 1 & 0\end{array}\right]$?
```{mcq}
#| label: char1-2
#| show_hints: true
1. $\pm\sqrt{2}/2$ [ correct hint: Good job! ]
2. $\pm\sqrt{2}$
3. $\pm 2$
4. $\pm\frac{1}{2}$
```Exercise 4 We will work with the linear dynamical system
\[\begin{eqnarray} \partial_t x & = -x\\ \partial_t y & = -11 x - 2 y\ .\\ \end{eqnarray}\]
Using integrateODE() find the solution \(x(t)\) for \(0 \leq t \leq 10\) for each of these three initial conditions:
Trajectory 1) \(x_0=1, y_0=-1\)
Trajectory 2) \(x_0=0, y_0=1\)
Trajectory 3) \(x_0=1, y_0=0\)
T1 <- integrateODE( ...., x=1, y=-1, bounds(t=0:10))
T2 <- integrateODE( ...., x=0, y=1, bounds(t=0:10))
T3 <- integrateODE( ...., x=1, y=0, bounds(t=0:10))The three solutions will be the functions T1$x(), T2$x(), and T3$x().
Since the dynamics are linear, there will be a simple relationship between the three solutions, just as there is a simple relationship between three vectors in a 2-dimensional vector space.
Find this simple relationship and make a graph that demonstrate it.
Exercise 5 The only fixed point of the rabbit-fox system is \(\left(r^\star, f^\star\right) = \left(\frac{\alpha}{\beta}, -\frac{\delta}{\gamma}\right)\)
The linearized dynamics around this fixed point are: \[\partial_t r = \left(\alpha + \frac{\beta\gamma}{\delta}\right) [r - r^\star] + \alpha [f - f^\star]\\ \partial_t f = - \delta [r - r^\star] + \left(\delta + \frac{\alpha\beta}{\delta}\right) [f - f^\star]\]
Using the formula for eigenvalues, calculate the eigenvalues symbolically and report on the stability of the fixed point.
Exercise 6 In this exercise, you will be using an app that allows you to set \(a\) and \(b\) in the [ab10] matrix and hunt for eigenvectors.
Open up the app in another browser tab and arrange it side-by-side with this document, so you can see both at the same time.
You used a somewhat similar app when exploring linear, two-dimensional finite-difference equations. Here’s a link, but you don’t need to open that app for this exercise.
This new app does a few things differently:
- In the old app, you had four number-entry boxes to set the values of the [abcd] matrix. Here, because we are using the [ab10] format, you need only set \(a\) and \(b\). You can do this by clicking within the left-hand graph.
- The a-b selector graph is annotated with the type of generic behavior that the differential-equation system will show for any combination of \(a\) and \(b\). This allows you to reason backward from the behavior you are interested in to the corresponding values of \(a\) and \(b\). Since the type of behavior is encoded in the eigenvalues, you are effectively inverting the eigenvalue formula to find \(a\) and \(b\).
- The flow field (right-hand graph) shows stream-lines of continuous-time motion. Imagine you had a bunch of glitter which you cast out over a pool of water with a given flow pattern. If you took a time-exposure photo, you would see each piece of glitter as it moved along with the water in which it is embedded.
- You can click in the flow field to start a trajectory at that point. You can control the time domain—that is, how long is the time exposure of the photo—with the radio buttons at the bottom.
- There is a thin, green line annotating the flow field. This line runs from the fixed point (right in the center) in the direction of the point on which you click to start a trajectory. A vector with the same orientation as the green line is given numerically below the graph. This orientation is also described using an angle measurement in degrees.
- There is no display of the solution, that is, \(x()\) versus \(t\). This is because we want you to understand the geometry of the flow.
First, play around with the app. Choose miscellaneous values of \(a\) and \(b\), observe the flow, and click in the flow field to draw a red trajectory. Notice what happens when you make the number of steps larger or smaller. In unstable directions, you can get a good view of the trajectory by starting it near the fixed point. If you start it at the periphery of the graph, most of the trajectory will be outside the display. And if the number of steps is very small, you may not see it at all because the glyph is so thin near the start.
When coming from a stable direction, click near the periphery of the graph to get a good view, since the red trajectory will naturally be heading in towards the origin. When estimating the orientation of an eigenvector, try to get the green line heading straight down the middle of the red trajectory. This often requires a bit of fiddling, moving the starting point just a tiny amount to find one that is just right.
When you are on an eigenvector, the trajectory will be dead straight on the green line. Any trajectory not starting on an eigenvector will be curved to a greater or lesser extent, depending on the relative sizes of the eigenvalues.
In the following questions, pick the choice closest to the answer you found. The exact choice of \(a\) and \(b\) stemming from a click depends somewhat on the graphics resolution of your display, so you should expect only to be able to get close to the \(a\) and \(b\) values given in the question.
Set $a=-1.0$ and $b = 0.5$ (or as close as you can get) and observe the flow. This is a saddle: stable in one direction, unstable in another. There will be one stable eigenvector, where the flow moves **directly outward** from the origin, and one unstable eigenvector, where the flow moves **directly inward** toward the origin. Which one of these is the orientation (in degrees) of the unstable eigenvector?
```{mcq}
#| label: ef1-1
#| show_hints: true
1. -149.7
2. -63.6
3. -31.7
4. 69.1
5. 73.5
6. 106.6
7. 143.7
8. 149
```Which of these best describes the eigenvalue corresponding to the unstable eigenvector?
```{mcq}
#| label: ef1-a
#| show_hints: true
1. Negative and real
2. Zero and real
3. Positive and real [ correct hint: right-o ]
4. Complex with a negative real part
5. Complex with a positive real part
```Keeping $a=-1.0$ and $b = 0.5$, hunt for the stable eigenvector. Which one of these is the orientation (in degrees) of the stable eigenvector?
```{mcq}
#| label: ef1-2
#| show_hints: true
1. -149.7
2. -63.6
3. -31.7
4. 69.1
5. 73.5
6. 106.6
7. 143.7
8. 149
```Set $a=2$ and $b = -0.5$ and observe the flow. This is a source: all initial conditions lead to motion away from the fixed point. But for this $a$ and $b$, some directions are "faster" than others. Which one of these is the orientation (in degrees) of the fast eigenvector?
```{mcq}
#| label: ef2-1
#| show_hints: true
1. -149.7
2. -63.6
3. -31.7
4. 69.1
5. 73.5
6. 106.6
7. 143.7
8. 149
```Keeping $a=2$ and $b = -0.5$, what is the orientation of the slow eigenvector?
```{mcq}
#| label: ef2-2
#| show_hints: true
1. -149.7
2. -63.6
3. -31.7
4. 69.1
5. 73.5
6. 106.6
7. 143.7
8. 149
```Staying with the flow from $a=2$ and $b = -0.5$, which of these best describes the pair of eigenvalues?
```{mcq}
#| label: ef2-a
#| show_hints: true
1. Both are real and positive [ correct hint: Yes ]
2. Both are real and negative
3. One is complex, the other real
4. Both are complex, one has negative real part
```Set $a=-2$ and $b = -0.5$ and observe the flow. This is a sink: all initial conditions lead to motion toward the fixed point. Which one of these is the orientation (in degrees) of the fast eigenvector?
```{mcq}
#| label: ef3-1
#| show_hints: true
1. -149.7
2. -63.6
3. -31.7
4. 69.1
5. 73.5
6. 106.6
7. 143.7
8. 149
```Keeping $a=-2$ and $b = -0.5$, what is the orientation of the slow eigenvector?
```{mcq}
#| label: ef3-2
#| show_hints: true
1. -149.7
2. -63.6
3. -31.7
4. 69.1
5. 73.5
6. 106.6
7. 143.7
8. 149
```Staying with the flow from $a=2$ and $b = -0.5$, which of these best describes the pair of eigenvalues?
```{mcq}
#| label: ef3-a
#| show_hints: true
1. Both are real and positive [ correct hint: Excellent! ]
2. Both are real and negative
3. One is complex, the other real
4. Both are complex, one has negative real part
5. Both are complex, both have negative real part
```Set $a=-1$ and $b = -1$ and observe the flow. This is a node: the flow is spiral with the fixed point at the center. Is the flow stable or unstable?
```{mcq}
#| label: ef4-1
#| show_hints: true
1. stable [ correct hint: Good job! ]
2. unstable
3. neither stable nor unstable [ hint: Is it getting closer to the fixed point or further away? that is the question that determines stability. ]
```Keeping $a=-1$ and $b = -1$, what can you correctly say about the eigenvectors?
```{mcq}
#| label: ef4-2
#| show_hints: true
1. They have almost the same orientation
2. There is only one eigenvector [ hint: There are almost aways two eigenvectors. ]
3. There are no eigenvectors [ correct hint: right-o This is true in the sense that there is no place where the motion is on a straight line toward or away from the origin. ]
```Staying with the flow from $a=2$ and $b = -0.5$, which of these best describes the pair of eigenvalues?
```{mcq}
#| label: ef4-a
#| show_hints: true
1. Both are real and positive [ correct hint: Nice. ]
2. Both are real and negative
3. One is complex, the other real
4. Both are complex, one has negative real part
5. Both are complex, both have negative real part
```Set $a=0$ and some $b < 0$. The result will be a trajectory that goes around in a closed loop. (You will want to make the number of steps large to see a complete loop.) You may not be able to click on exactly $a=0$, but you should be able to get close enough that the trajectory is a very tightly wound spiral. We are interested here in how many cycles the trajectory covers in a fixed number of time steps, that is, whether the trajectory moves quickly or slowly around the loop. Which of these statements is true about the relationship between $b$ and the speed of the motion?
```{mcq}
#| label: ef5-a
#| show_hints: true
1. There is no relationship between $b$ and speed of oscillation.
2. When $b$ is closer to zero, the speed is higher.
3. When $b$ is further from zero, the speed is higher. [ correct hint: Correct ]
```we will look at the relationship between the eigenvalues and the speed of oscillation. Keep $a=0$ and $b < 0$, which guarantees that the motion will be circular. Vary $b$, each time observing how many loops the trajectory makes (for a fixed number of steps), and what the eigenvalues are. Which of these statements is **not** true about the speed of oscillation and the eigenvalues.
```{mcq}
#| label: ef5-b
#| show_hints: true
1. The larger the magnitude of the imaginary part of the eigenvalues, the faster the oscillation.
2. The two eigenvalues have the same imaginary part, except one is positive and the other negative.
3. The speed of oscillation relates to the *difference* between the imaginary part of the eigenvalues. [ correct hint: Excellent! ]
4. The real part of the eigenvalues is the same for both when there is oscillatory motion.
```Exercise 7 ./Block_5/wolf-tear-bowl.Rmd Does this exist?
No answers yet collected