Ordinary Differential Equations/Locally linear

We will study autonomous systems $$\mathrm {x} '=\mathrm {f} (\mathrm {x} ),$$ where the components of $\mathrm {f}$ are $C^{1}$ functions so that we are able to Taylor expand them to first order. A system of the form $$\mathrm {x} '=\mathrm {A} \mathrm {x} +\mathrm {g} (\mathrm {x} )$$ is called locally linear around a critical point $\mathrm {x} _{0}$ of $\mathrm {f}$ if $${\frac {\left\|\mathrm {g} (\mathrm {x} )\right\|}{\left\|\mathrm {x} -\mathrm {x} _{0}\right\|}}\to 0{\text{ as }}\mathrm {x} \to \mathrm {x} _{0}.$$

Example presenting the method

We study the damped oscillating pendulum system: $${\frac {dx}{dt}}=y,\,{\frac {dy}{dt}}=-\gamma y-\omega ^{2}\sin(x),$$ where $\gamma$ is called the damping constant and as in the spring problem it is responsible for removing energy.

1. First we find the critical points. From the previous section we have:$$(k\pi ,0){\text{ for any integer }}k.$$
2. Second we Taylor expand the RHS of the system $\mathrm {F} (x,y):={\bigl (}{\begin{smallmatrix}y\\\\-\gamma y-\omega ^{2}\sin(x)\end{smallmatrix}}{\bigr )}$ around arbitrary critical point $(x_{0},y_{0})$: $${\begin{aligned}\mathrm {F} (x,y)&=\mathrm {F} (x_{0},y_{0})+\mathrm {J} _{\mathrm {F} }(x_{0},y_{0}){\binom {x-x_{0}}{y-y_{0}}}+o(\left\|(x-x_{0},y-y_{0})\right\|)\\&={\bigl (}{\begin{smallmatrix}0&1\\&\\-\omega ^{2}\cos(x_{0})&-\gamma \end{smallmatrix}}{\bigr )}{\binom {x-x_{0}}{y-y_{0}}}+o(\left\|(x-x_{0},y-y_{0})\right\|).\end{aligned}}$$
3. Here $\mathrm {J} _{\mathrm {F} }(x_{0},y_{0})$ is the Jacobian matrix at $(x_{0},y_{0})$ which, for function $\mathrm {F} (x,y)={\binom {\mathrm {F} _{1}(x,y)}{\mathrm {F} _{2}(x,y)}}$, is defined as:$${\begin{aligned}J_{\mathrm {F} }(x_{0},y_{0}):={\bigl (}{\begin{smallmatrix}{\frac {\mathrm {d} \mathrm {F} _{1}}{\mathrm {d} x}}(x_{0},y_{0})&{\frac {\mathrm {d} \mathrm {F} _{1}}{\mathrm {d} y}}(x_{0},y_{0})\\&\\{\frac {\mathrm {d} \mathrm {F} _{2}}{\mathrm {d} x}}(x_{0},y_{0})&{\frac {\mathrm {d} \mathrm {F} _{2}}{\mathrm {d} y}}(x_{0},y_{0})\end{smallmatrix}}{\bigr )}\end{aligned}}$$
4. The linearization around $(x_{0},y_{0})=(k\pi ,0)$ for an even integer $k$ is: $${\frac {\mathrm {d} }{\mathrm {d} t}}{\binom {x}{y}}={\bigl (}{\begin{smallmatrix}0&1\\&\\-\omega ^{2}&-\gamma \end{smallmatrix}}{\bigr )}{\binom {x-k\pi }{y}}+o(\left\|(x-k\pi ,y)\right\|).$$
5. The eigenvalues of that matrix are: $$\lambda _{1},\,\lambda _{2}={\frac {-\gamma \pm {\sqrt {\gamma ^{2}-4\omega ^{2}}}}{2}}.$$
6. If $\gamma ^{2}-4\omega ^{2}>0$, then the eigenvalues are real, distinct, and negative. Therefore, the critical points will be stable nodes.We observe that the basins of attractions for each even-integer critical points are well-separated.
7. If $\gamma ^{2}-4\omega ^{2}=0$, then the eigenvalues are repeated, real, and negative. Therefore, the critical points will be stable nodes.
8. If $\gamma ^{2}-4\omega ^{2}<0$, then the eigenvalues are complex with negative real part. Therefore, the critical points will be stable spiral sinks.
9. The linearization around $(x_{0},y_{0})=(k\pi ,0)$ for odd integer $k$ is:$${\frac {\mathrm {d} }{\mathrm {d} t}}{\binom {x}{y}}={\bigl (}{\begin{smallmatrix}0&1\\&\\\omega ^{2}&-\gamma \end{smallmatrix}}{\bigr )}{\binom {x-k\pi }{y}}+o(\left\|(x-k\cdot \pi ,y)\right\|).$$
10. The eigenvalues of that matrix are: $$\lambda _{1},\,\lambda _{2}={\frac {-\gamma \pm {\sqrt {\gamma ^{2}+4\omega ^{2}}}}{2}}.$$
11. Therefore, it has one negative eigenvalue $\lambda _{1}<0$ and one positive eigenvalue $\lambda _{2}>0$, and so the critical points will be unstable saddle points.