Ec. DiferencialesSystems of Differential Equations

Systems with Complex Eigenvalues

Real and imaginary parts give spiraling solutions.

Sin empezar Ir a la práctica

Esta lección aún no está traducida. Por ahora se muestra en inglés; la versión en español está en camino.

When $A$ has complex eigenvalues $\alpha \pm \beta i$, the eigenvector method still works, but it produces complex solutions. As with complex roots of one equation, their real and imaginary parts are the real solutions, and they spiral.

Idea clave

If $A\mathbf{v} = \lambda\mathbf{v}$ with $\lambda = \alpha + \beta i$ and $\mathbf{v} = \mathbf{a} + i\mathbf{b}$, then $$\mathbf{u}_1 = e^{\alpha t}(\mathbf{a}\cos\beta t - \mathbf{b}\sin\beta t),$$ $$\mathbf{u}_2 = e^{\alpha t}(\mathbf{a}\sin\beta t + \mathbf{b}\cos\beta t)$$ are real solutions, and $\mathbf{x} = c_1\mathbf{u}_1 + c_2\mathbf{u}_2$. Trajectories spiral in if $\alpha < 0$, out if $\alpha > 0$, and close up if $\alpha = 0$.

Por qué funciona

$e^{\lambda t}\mathbf{v}$ $= e^{\alpha t}(\cos\beta t + i\sin\beta t)(\mathbf{a}$ $+ i\mathbf{b})$. Multiplying out, the real part is $\mathbf{u}_1$ and the imaginary part is $\mathbf{u}_2$. Because $A$ is real, the equation $(\mathbf{u}_1 + i\mathbf{u}_2)' = A(\mathbf{u}_1 + i\mathbf{u}_2)$ splits into a real part and an imaginary part that must each balance, so $\mathbf{u}_1$ and $\mathbf{u}_2$ are both solutions. The direction of rotation is easy to read without any of this: at the point $(1, 0)$ the velocity is the first column of $A$, so if $a_{21} > 0$ the motion turns upward there, counterclockwise.

Ejemplo resuelto: Eigenvalues and rotation

Describe the trajectories of $\mathbf{x}' = \begin{bmatrix} -1 & -2 \\ 2 & -1 \end{bmatrix}\mathbf{x}$.
  1. $\tau = -2$, $\Delta = 5$: $\lambda^{2} + 2\lambda + 5 = 0$, so $\lambda = -1 \pm 2i$.
  2. $\alpha = -1 < 0$: they spiral in.
  3. At $(1, 0)$ the velocity is $(-1, 2)$: moving up, so counterclockwise.

Respuesta: spirals in, counterclockwise

Ejemplo resuelto: Solve with an initial condition

Solve the same system with $\mathbf{x}(0) = \begin{bmatrix} 1 \\ 0 \end{bmatrix}$.
  1. For $\lambda = -1 + 2i$: $(A - \lambda I)\mathbf{v} = 0$ has first row $-2i v_1 - 2v_2 = 0$, so $\mathbf{v} = \begin{bmatrix} 1 \\ -i \end{bmatrix}$: $\mathbf{a} = \begin{bmatrix} 1 \\ 0 \end{bmatrix}$, $\mathbf{b} = \begin{bmatrix} 0 \\ -1 \end{bmatrix}$.
  2. $\mathbf{u}_1 = e^{-t}\begin{bmatrix} \cos 2t \\ \sin 2t \end{bmatrix}$ and $\mathbf{u}_2 = e^{-t}\begin{bmatrix} \sin 2t \\ -\cos 2t \end{bmatrix}$.
  3. At $t = 0$: $c_1\begin{bmatrix} 1 \\ 0 \end{bmatrix} + c_2\begin{bmatrix} 0 \\ -1 \end{bmatrix} = \begin{bmatrix} 1 \\ 0 \end{bmatrix}$, so $c_1 = 1$, $c_2 = 0$.

Respuesta: $x_1 = e^{-t}\cos 2t$, $x_2 = e^{-t}\sin 2t$: a point circling counterclockwise as it shrinks toward the origin

Ejemplo resuelto: A center

Describe $\mathbf{x}' = \begin{bmatrix} 0 & 4 \\ -1 & 0 \end{bmatrix}\mathbf{x}$.
  1. $\tau = 0$, $\Delta = 4$: $\lambda = \pm 2i$, so $\alpha = 0$.
  2. At $(1, 0)$ the velocity is $(0, -1)$: clockwise.

Respuesta: closed loops (ellipses), traversed clockwise with period $\pi$

Error común

Judging stability from $\beta$ or from the eigenvector. Only the real part $\alpha$ (half the trace) decides spiral in or out; $\beta$ sets how fast it turns. And the direction of rotation is not in the eigenvalue's sign of $\beta$: check the velocity at one point.
The spiral sink of the examples. Press Run to send a particle along the trajectory; the strip below shows $x_1(t)$ and $x_2(t)$ as decaying oscillations a quarter period apart.

Práctica

Apruebas con 3 respuestas correctas, las 2 últimas seguidas. Cada respuesta correcta da 3 XP.

Cargando la práctica…

Se apoya en

Lleva a

Este es uno de los peldaños más altos de Ecuaciones Diferenciales.

Fuentes y lecturas

  • Paul's Online Math Notes, Differential Equations: Complex Eigenvalues
  • MIT OpenCourseWare 18.03SC Differential Equations (Fall 2011), Unit IV: First-order Systems