Difference between revisions of "CDS 140a Winter 2015 Homework 3"
(4 intermediate revisions by the same user not shown)  
Line 8:  Line 8:  
 pdf = cds140wi15_hw3.pdf   pdf = cds140wi15_hw3.pdf  
}} __MATHJAX__  }} __MATHJAX__  
{{warningThis problem set is under construction. This banner will be removed when the problems are finalized.}}  <! {{warningThis problem set is under construction. This banner will be removed when the problems are finalized.}} >  
'''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  '''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  
Line 15:  Line 15:  
<ol>  <ol>  
<! 2014 TA comments: A lot of people also were confused about solving for the unstable manifold by simply letting xdot =  f(x) on #1. >  <! 2014 TA comments: A lot of people also were confused about solving for the unstable manifold by simply letting xdot =  f(x) on #1. >  
<li>'''Perko, Section 2.7, problem 2'''  <li>'''Perko, Section 2.7, problem 2, 3'''<br>  
Find the first three successive approximations $u^{(1)}(t,a)$, $u^{(2)}(t,a)$, and $u^{(3)}(t,a)$ for  (a) Find the first three successive approximations $u^{(1)}(t,a)$, $u^{(2)}(t,a)$, and $u^{(3)}(t,a)$ for  
<center><amsmath>  <center><amsmath>  
\aligned  \aligned  
Line 24:  Line 24:  
</amsmath></center>  </amsmath></center>  
and use $u^{(3)}(t,a)$ to approximate $S$ near the origin. Also approximate the unstable manifold $U$ near the origin for this system. Note that $u^{(2)}(t,a)=u^{(3)}(t,a)$ and therefore $u^{(j+1)}(t,a)=u^{(j)}(t,a)$ for $j\geq 2$. Thus $u(t,a)=u^{(2)}(t,a)$ which gives the exact function defining $S$.  and use $u^{(3)}(t,a)$ to approximate $S$ near the origin. Also approximate the unstable manifold $U$ near the origin for this system. Note that $u^{(2)}(t,a)=u^{(3)}(t,a)$ and therefore $u^{(j+1)}(t,a)=u^{(j)}(t,a)$ for $j\geq 2$. Thus $u(t,a)=u^{(2)}(t,a)$ which gives the exact function defining $S$.  
(b) Sketch $S$, $U$, $E^s$ and $E^u$.  
</li>  </li>  
Line 42:  Line 44:  
* Note: The result holds for ${\Bbb R^n}$, but you only need to consider the case $n = 2$ and $k = 1$ (although you are free to prove the more general case if you prefer).  * Note: The result holds for ${\Bbb R^n}$, but you only need to consider the case $n = 2$ and $k = 1$ (although you are free to prove the more general case if you prefer).  
* Hint: One way to show $S$ is an invariant manifold in $\mathbb R^2$ is to show that the normal vector (orthogonal to the tangent to $S$ at $(x,h(x))$) is orthogonal to the vector field $(f,g)$ at that point. (It is sufficient to prove the result for $\mathbb R^2$.)  * Hint: One way to show $S$ is an invariant manifold in $\mathbb R^2$ is to show that the normal vector (orthogonal to the tangent to $S$ at $(x,h(x))$) is orthogonal to the vector field $(f,g)$ at that point. (It is sufficient to prove the result for $\mathbb R^2$.)  
</li>  </li>  
Line 77:  Line 58:  
is transformed into a linear system using this transformation, i.e., if $y = H(x)$, show that $\dot y = A x$ for an appropriate $A$. Use this transformation to compute and sketch the stable and unstable manifolds for the nonlinear system.  is transformed into a linear system using this transformation, i.e., if $y = H(x)$, show that $\dot y = A x$ for an appropriate $A$. Use this transformation to compute and sketch the stable and unstable manifolds for the nonlinear system.  
* Note: since the nonlinear system can be transformed into a linear one, it follows that $S = H^{1}(E^s)$ and $U = H^{1}(E^u)$.  * Note: since the nonlinear system can be transformed into a linear one, it follows that $S = H^{1}(E^s)$ and $U = H^{1}(E^u)$.  
</li>  </li>  
Line 114:  Line 70:  
\endaligned  \endaligned  
</amsmath></center>  </amsmath></center>  
and show that the successive approximations $\Phi^(k)\rightarrow\Phi$ and $\Psi^(k)\rightarrow\Psi$ as $k\rightarrow\infty$ for all $x=(y_1,y_2,z)\in{\mathbb R}^3$. Define $H_0=(\Phi,\Psi)^T$ and use this homeomorphism to find  and show that the successive approximations $\Phi^{(k)}\rightarrow\Phi$ and $\Psi^{(k)}\rightarrow\Psi$ as $k\rightarrow\infty$ for all $x=(y_1,y_2,z)\in{\mathbb R}^3$. Define $H_0=(\Phi,\Psi)^T$ and use this homeomorphism to find  
<center><amsmath>  <center><amsmath>  
H (x) =\int_0^1 e^{As} H_0 \bigl(\phi_s(x) \bigr) ds,  H (x) =\int_0^1 e^{As} H_0 \bigl(\phi_s(x) \bigr) ds, 
Latest revision as of 06:15, 22 January 2015
R. Murray  Issued: 19 Jan 2015 
CDS 140, Winter 2015  Due: 28 Jan 2015 at 12:30 pm In class or to box across 107 STL 
__MATHJAX__
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 2.7, problem 2, 3
(a) Find the first three successive approximations $u^{(1)}(t,a)$, $u^{(2)}(t,a)$, and $u^{(3)}(t,a)$ for<amsmath> \aligned \dot{x}_1&=x_1,\\ \dot{x}_2&=x_2+x_1^2 \endaligned
</amsmath>and use $u^{(3)}(t,a)$ to approximate $S$ near the origin. Also approximate the unstable manifold $U$ near the origin for this system. Note that $u^{(2)}(t,a)=u^{(3)}(t,a)$ and therefore $u^{(j+1)}(t,a)=u^{(j)}(t,a)$ for $j\geq 2$. Thus $u(t,a)=u^{(2)}(t,a)$ which gives the exact function defining $S$.
(b) Sketch $S$, $U$, $E^s$ and $E^u$.
 Consider a dynamical system with $x = (u, v) \in {\Bbb R}^n$. For the case $n = 2$, prove that if
<amsmath> \aligned \dot{u}&=f(u,v),\qquad u\in\mathbb{R}^k\\ \dot{v}&=g(u,v),\qquad v\in\mathbb{R}^{nk} \endaligned
</amsmath>then the manifold $S=\{(u,v)\in\mathbb R^k\times\mathbb R^{nk}v=h(u)\}$ is an invariant manifold of the system if
<amsmath> g(u,h(u))=Dh(u)f(u,h(u))
</amsmath>Use this result to compute the stable manifold for the system of problem 1 above using the Taylor series for $h(x)$ to define $S=\{(x_1,x_2)x_2=h(x)\}$ and matching coefficients to solve for $h(x)$.
 Note: The result holds for ${\Bbb R^n}$, but you only need to consider the case $n = 2$ and $k = 1$ (although you are free to prove the more general case if you prefer).
 Hint: One way to show $S$ is an invariant manifold in $\mathbb R^2$ is to show that the normal vector (orthogonal to the tangent to $S$ at $(x,h(x))$) is orthogonal to the vector field $(f,g)$ at that point. (It is sufficient to prove the result for $\mathbb R^2$.)
 Perko, Section 2.6, problem 3
Show that the continuous map $H:{\Bbb R}^3 \to {\Bbb R}^3$ defined by
<amsmath> \begin{bmatrix} y_1 \\ y_2 \\ y_3 \end{bmatrix} = H(x) = \begin{bmatrix} x_1 \\ x_2 + x_1^2 \\ x_3 + x_1^2 / 3 \end{bmatrix}
</amsmath>has a continuous inverse $H^{1}:{\Bbb R}^3 \to {\Bbb R}^3$ and that the nonlinear system
<amsmath> \frac{dx}{dt} = \begin{bmatrix} x_1 \\ x_2 + x_1^2 \\ x_3 + x_1^2 \end{bmatrix}
</amsmath>is transformed into a linear system using this transformation, i.e., if $y = H(x)$, show that $\dot y = A x$ for an appropriate $A$. Use this transformation to compute and sketch the stable and unstable manifolds for the nonlinear system.
 Note: since the nonlinear system can be transformed into a linear one, it follows that $S = H^{1}(E^s)$ and $U = H^{1}(E^u)$.

Perko, Section 2.8, Problem 1
Solve the system
<amsmath> \aligned \dot{y}_1&=y_1\\ \dot{y}_2&=y_2+z^2\\ \dot{z}&=z \endaligned
</amsmath>and show that the successive approximations $\Phi^{(k)}\rightarrow\Phi$ and $\Psi^{(k)}\rightarrow\Psi$ as $k\rightarrow\infty$ for all $x=(y_1,y_2,z)\in{\mathbb R}^3$. Define $H_0=(\Phi,\Psi)^T$ and use this homeomorphism to find
<amsmath> H (x) =\int_0^1 e^{As} H_0 \bigl(\phi_s(x) \bigr) ds,
</amsmath>where $A$ is the linearization of the nonlinear dynamics at the origin and $\phi_t(x)$ is the flow of the full system. Use the homemorphism $H$ to find the stable and unstable manifolds
<amsmath> W^s(0)=H^{1}(E^s)\quad {\mathrm{and}}\quad W^u(0)=H^{1}(E^u)
</amsmath>for this system.
 Note: this problem involves some simple but somewhat tedious computations. If you know how to use Mathematica or a similar program, you may wish to carry out the computations for the successive approximations using that software. However, make sure to show the results at each step of the calculation.
 Hint: You should find
<amsmath> H(y_1,y_2,z) = \begin{bmatrix} y_1 \\ y_2z^2/3 \\ z \end{bmatrix}, \qquad \aligned W^s(0) &=\{x\in{\mathbb{R}}z=0\}, \\ W^u(0) &=\{x\in{\mathbb{R}}y_1=0,y_2=z^2/3\}. \endaligned
</amsmath>
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 as the definitive source of the problem statement.