CDS 140a Winter 2016 Homework 1
D. MacMartin and J. Doyle  Issued: 4 Jan 2016 
CDS 140, Winter 2016
{{#if: cds140wi16_hw1.pdf  (PDF) }} 
Due: 13 Jan 2015 at 5:00 pm In class or to box across 107 STL (?) 
Note: In the upper left hand corner of the second page of your homework set, please put the number of hours that you spent on this homework set (including reading).
 Perko, Section 1.2, Exercise 6: Let the $n \times n$ matrix $A$ have real, distinct eigenvalues. Let $\phi(t, x)$ the the solution of the initial value problem \begin{align} \dot x &= A x &\qquad x(0) &= x_0. \end{align} Show that for each fixed $t \in {\mathbb R}$, \begin{align} \lim_{y_0 \to x_0} \phi(t, y_0) = \phi(t, x_0). \end{align} This shows that the solution $\phi(t, x_0)$ is a continuous function of the initial condition.
 (Based on Perko, Section 1.3, Exercise 5, 6)
(a) For each matrix below, find the eigenvalues and eigenvectors of $A$ and $e^A$:(i) $\begin{bmatrix} a & 0 \\ 0 & b \end{bmatrix}$ (ii) $\begin{bmatrix} 1 & 0 \\ a & 1 \end{bmatrix}$ (iii) $\begin{bmatrix} 5 & 6 \\ 3 & 4 \end{bmatrix}$ (remember to show the steps required for these [simple!] computations, don't just plug in values from MATLAB or Mathematica; see notes at the bottom of the page).
(b) Show that if $x$ is the eigenvector of $A$ corresponding to the eigenvalue $\lambda$, then $x$ is also an eigenvector of $e^A$ corresponding to the eigenvalue $e^\lambda$.
 (Based on Perko, Section 1.4, Exercise 6) Let $A:{\mathbb R}^n \to {\mathbb R}^n$ be a linear transformation that leaves a subspace $E \subset {\mathbb R}^n$ invariant (i.e., for all $x \in E$, $A x \in E$). Show that if $x(t)$ is the solution of the initial value problem \begin{align} \dot x &= A x &\qquad x(0) &= x_0 \end{align} with $x_0 \in E$, then $x(t) \in E$ for all $t \in {\mathbb R}$.
 Perko, Section 1.6, Exercise 2: Solve the initial value problem \begin{align} \dot x &= A x &\qquad x(0) &= x_0 \end{align} with \begin{align} A &= \begin{bmatrix} 0 & 2 & 0 \\ 1 & 2 & 0 \\ 0 & 0 & 2 \end{bmatrix}. \end{align} Determine the stable and unstable subspaces and sketch the (3D) phase portrait. (Hint: see Figure 1 in Section 1.6 for an example of a 3D phase portrait.)
 Perko, Section 1.9, problem 5, parts (c), (d2): Let <amsmath>A</amsmath> be an <amsmath>n \times n</amsmath> nonsingular matrix and let <amsmath>x(t)</amsmath> be the solution of the initial value problem (1) with <amsmath>x(0) = x_0</amsmath>. Show that
(c) If <amsmath>x_0 \in E^c</amsmath>, <amsmath>x_0 \neq 0</amsmath> and <amsmath>A</amsmath> is semisimple, then there are postive constants <amsmath>m</amsmath> and <amsmath>M</amsmath> such that for all <amsmath>t \in R</amsmath>, <amsmath>m \leq x(t) \leq M</amsmath>; Note: in the book, Perko defines <amsmath>\sim</amsmath> to mean "set subtraction". So <amsmath>E \sim \{0\}</amsmath> in the book is the set <amsmath>E</amsmath> minus the point 0.

Consider the system
<amsmath> \frac{dx}{dt} = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix} x + \begin{bmatrix} 0 \\ 1 \end{bmatrix}u, \qquad y = \begin{bmatrix} 1 & 0 \end{bmatrix} x.
</amsmath>(a) Show that the unforced system (<amsmath>u = 0</amsmath>) is stable but not asymptotically stable.
(b) Given <amsmath>x(0) = x_0</amsmath> and <amsmath>u(t) = \cos(\omega*t)</amsmath>, solve for the output <amsmath>y(t)</amsmath>. Show that when <amsmath>\omega=1</amsmath> the output is unbounded.
Notes:
 The problems are transcribed above in case you don't have access to Perko. However, in the case of discrepancy, you should use Perko (third edition) as the definitive source of the problem statement.
 There are a number of problems that can be solved using MATLAB, Mathematica or a similar program. If you just give the answer with no explanation (or say "via MATLAB"), the TAs will take off points. Instead, you should show how the solutions can be worked out by hand, along the lines of what is done in the text book. It is fine to check everything with MATLAB or your favorite software tool.
 For numerical calculations, it is OK to use MATLAB to invert a matrix. But you should not use it to compute the matrix exponential and just put down the answer. Instead, show how to get the matrix exponential into a form in which the calculation can be done by hand and then carry out the computation.
 For phase portraits, you should generate the diagram by hand and make sure to label any important features. Describe why the portrait looks as it does based on the relevant properties of the dynamical system (eg, eigenvalues of the A matrix).
 For the final exam, you will not be allowed to use MATLAB, Mathematica or similar programs, so make sure you understand what you are computing and drawing!