myLearning

Topic 4 of 6

Step 20: Differential equations

Systems and phase portraits

Step 20 · Differential equations topics

Write linked differential equations as x′ = Ax, solve them with eigenvalues and eigenvectors, and read the long-run behavior from a phase portrait.

Systems as matrix equations

Two quantities that change together, like two populations or two tanks, give a system: x′ = ax + by and y′ = cx + dy. Write it as one matrix equation, x′ = Ax.

(xy)′=(abcd)(xy)\begin{pmatrix} x \\ y \end{pmatrix}' = \begin{pmatrix} a & b \\ c & d \end{pmatrix}\begin{pmatrix} x \\ y \end{pmatrix}

Example: Writing a system

Write x′ = 3x + y, y′ = x + 3y as a matrix equation.

  1. The coefficients form the rows.

Answerx′ = Ax with A = (3113)\begin{pmatrix} 3 & 1 \\ 1 & 3 \end{pmatrix}

Example: A second-order equation as a system

Write y″ + 3y′ + 2y = 0 as a system.

  1. Let x₁ = y and x₂ = y′.
  2. Then x₁′ = x₂ and x₂′ = −2x₁ − 3x₂.

Answerx′ = Ax with A = (01−2−3)\begin{pmatrix} 0 & 1 \\ -2 & -3 \end{pmatrix}

Eigenvalue solutions

If Av = λv, then x=eλtv\mathbf{x} = e^{\lambda t}\mathbf{v} solves x′ = Ax: it moves along the line through v, growing or shrinking. With two independent eigenvectors, the general solution combines two of these.

x(t)=c1eλ1tv1+c2eλ2tv2\mathbf{x}(t) = c_1e^{\lambda_1 t}\mathbf{v}_1 + c_2e^{\lambda_2 t}\mathbf{v}_2

Example: Two real eigenvalues

Solve x′ = Ax with A = (3113)\begin{pmatrix} 3 & 1 \\ 1 & 3 \end{pmatrix}.

  1. The eigenvalues are 4 and 2, with eigenvectors ⟨1, 1⟩ and ⟨1, −1⟩.
  2. Each gives a straight-line solution.

Answerx=c1e4t⟨1,1⟩+c2e2t⟨1,−1⟩\mathbf{x} = c_1e^{4t}\langle 1, 1\rangle + c_2e^{2t}\langle 1, -1\rangle

Example: Complex eigenvalues

What do solutions of x′ = Ax look like for A = (0−110)\begin{pmatrix} 0 & -1 \\ 1 & 0 \end{pmatrix}?

  1. The eigenvalues are ±i.
  2. x = cos t, y = sin t is one solution: it runs round the unit circle.

AnswerCircles around the origin

Phase portraits

A phase portrait draws solution paths in the xy-plane, with the direction of motion. The eigenvalues of A decide the picture around the origin, where the system is at rest.

EigenvaluesPictureStable?
Both negativeNode: all paths flow inYes
Both positiveNode: all paths flow outNo
Opposite signsSaddleNo
Complex, real part negativeSpiral inYes
Complex, real part positiveSpiral outNo
Purely imaginaryCenter: closed loopsNeither
A saddle, for x′ = x and y′ = −y: paths come in along the y-axis and leave along the x-axis, bending away from the origin.

Example: Classifying

Classify the origin for x′ = Ax with A = (−200−5)\begin{pmatrix} -2 & 0 \\ 0 & -5 \end{pmatrix}.

  1. The eigenvalues are −2 and −5, both negative.

AnswerA stable node: every path flows into the origin

Example: A spiral

Classify the origin for A = (−1−22−1)\begin{pmatrix} -1 & -2 \\ 2 & -1 \end{pmatrix}.

  1. The eigenvalues are −1 ± 2i.
  2. Their real part, −1, is negative.

AnswerA stable spiral: paths wind into the origin

Try it yourself

Work each one on paper first, then open the answer to check.

  1. Write x′ = 2x − y, y′ = 4x + y as x′ = Ax.

    Show answer

    (2−141)\begin{pmatrix} 2 & -1 \\ 4 & 1 \end{pmatrix}

  2. What are the eigenvalues of (200−3)\begin{pmatrix} 2 & 0 \\ 0 & -3 \end{pmatrix}, and what kind of point is the origin?

    Show answer

    2 and −3: a saddle

  3. The eigenvalues are −1 and −4. Is the origin stable?

    Show answer

    Yes: a stable node

  4. The eigenvalues are 2 ± 3i. What does the phase portrait show?

    Show answer

    An unstable spiral

  5. Av = 5v. Which solution of x′ = Ax does v give?

    Show answer

    e5tve^{5t}\mathbf{v}

  6. The eigenvalues are ±4i. What do the paths look like?

    Show answer

    Closed loops around the origin: a center

What next?