All revision notes topics

5.17 Phase portraitsIB Maths: Applications and Interpretation HL: Revision notes

Section 1

Coupled linear systems and the phase plane

Two populations, or two connected quantities, can be described by the coupled linear system dxdt=ax+by,dydt=cx+dy,\frac{dx}{dt}=ax+by,\qquad \frac{dy}{dt}=cx+dy, which in matrix form is ddt(xy)=M(xy)\frac{d}{dt}\begin{pmatrix} x \\ y \end{pmatrix}=M\begin{pmatrix} x \\ y \end{pmatrix} with M=(abcd)M=\begin{pmatrix} a & b \\ c & d \end{pmatrix}. A phase portrait is a diagram in the xyxy-plane showing trajectories: the paths (x(t),y(t))(x(t),y(t)) followed as tt increases, with arrows for direction. At each point (x,y)(x,y) the trajectory has direction (dxdtdydt)=M(xy)\begin{pmatrix} \frac{dx}{dt} \\ \frac{dy}{dt} \end{pmatrix}=M\begin{pmatrix} x \\ y \end{pmatrix}. For the systems in this course the origin is the only equilibrium point, where both rates are zero.

Key termstrajectoryphase portraitequilibrium point
Exam tip

To find the direction of travel at a point, substitute it into both equations. At (1,0)(1,0) the direction is (a,c)(a,c).

Section 2

Eigenvalues, eigenvectors and straight-line solutions

Find the eigenvalues λ\lambda of MM from det⁡(M−λI)=0\det(M-\lambda I)=0, then an eigenvector v\mathbf{v} for each from (M−λI)v=0(M-\lambda I)\mathbf{v}=\mathbf{0} (a GDC will do this too). You are only given systems with distinct, non-zero eigenvalues. If the system starts at a point on an eigenvector v\mathbf{v}, the solution is eλtve^{\lambda t}\mathbf{v}. It stays on that line through the origin, moving away if λ>0\lambda>0 and towards the origin if λ<0\lambda<0. These are the straight-line trajectories. Example: M=(4123)M=\begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix} has λ2−7λ+10=0\lambda^2-7\lambda+10=0, so λ=2,5\lambda=2,5 with eigenvectors (1−2)\begin{pmatrix} 1 \\ -2 \end{pmatrix} and (11)\begin{pmatrix} 1 \\ 1 \end{pmatrix}.

Key termseigenvalueeigenvector
Common mistake

Forgetting that an eigenvector can be scaled: (1−2)\begin{pmatrix} 1 \\ -2 \end{pmatrix} and (−24)\begin{pmatrix} -2 \\ 4 \end{pmatrix} give the same line.

Section 3

Exact solutions for real distinct eigenvalues

When the eigenvalues λ1≠λ2\lambda_1\neq\lambda_2 are real, the general solution is (xy)=Aeλ1tv1+Beλ2tv2.\begin{pmatrix} x \\ y \end{pmatrix}=Ae^{\lambda_1t}\mathbf{v}_1+Be^{\lambda_2t}\mathbf{v}_2. Use the initial conditions to find AA and BB. Example: M=(4123)M=\begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}, with x=3x=3 and y=0y=0 at t=0t=0. Then A+B=3A+B=3 and −2A+B=0-2A+B=0, so A=1A=1 and B=2B=2 (using λ1=2\lambda_1=2, v1=(1−2)\mathbf{v}_1=\begin{pmatrix} 1 \\ -2 \end{pmatrix}, λ2=5\lambda_2=5, v2=(11)\mathbf{v}_2=\begin{pmatrix} 1 \\ 1 \end{pmatrix}): x=e2t+2e5t,y=−2e2t+2e5t.x=e^{2t}+2e^{5t},\qquad y=-2e^{2t}+2e^{5t}. Exact solutions are only required in this case. For complex or imaginary eigenvalues you describe the behaviour qualitatively.

Key termsgeneral solution
Common mistake

Mixing up which eigenvector goes with which eigenvalue. Write each pair next to its exponent.

Section 4

Real eigenvalues: nodes and saddle points

The signs of the real eigenvalues decide the portrait:

  • Both positive: every solution moves away from the origin. The origin is an unstable node (a source).
  • Both negative: every solution moves towards the origin. The origin is a stable node (a sink). Trajectories approach along the eigenvector whose eigenvalue is smaller in magnitude.
  • Opposite signs: the origin is a saddle point. Solutions on the eigenvector line of the negative eigenvalue move in to the origin, but every other solution is turned away and eventually moves off along the eigenvector of the positive eigenvalue. To sketch: draw the eigenvector lines first with arrows, then add curved trajectories that follow them.
Key termssaddle pointstable node
Exam tip

A saddle point is unstable: only solutions that start exactly on one line reach the origin.

Section 5

Complex and imaginary eigenvalues: spirals and ellipses

If the eigenvalues are complex, λ=p±qi\lambda=p\pm qi with q≠0q\neq0, the solutions rotate around the origin:

  • p>0p>0: spiral away from the origin.
  • p<0p<0: spiral towards the origin.
  • p=0p=0 (purely imaginary): closed circles or ellipses around the origin, with neither growth nor decay. The real part pp decides growth or decay, and the imaginary part qq gives the rotation. To find the sense of rotation, test one point: at (1,0)(1,0) the direction is (a,c)(a,c), so c>0c>0 means anticlockwise and c<0c<0 means clockwise. Example: dxdt=−x−4y\frac{dx}{dt}=-x-4y, dydt=x−y\frac{dy}{dt}=x-y has λ=−1±2i\lambda=-1\pm2i, so a spiral towards the origin; at (1,0)(1,0) the direction is (−1,1)(-1,1), anticlockwise.
Key termsspiralimaginary eigenvalue
Common mistake

Saying complex eigenvalues always give spirals out. The sign of the real part decides inwards or outwards.

Section 6

Reading a phase portrait in context

A phase portrait answers the question: what happens to the two quantities in the long term?

  • Stable populations: a sink or an inward spiral means both variables settle at the equilibrium (here the origin).
  • Unstable: a source or an outward spiral means the variables grow without bound.
  • Cycles: imaginary eigenvalues give populations that rise and fall repeatedly.
  • Saddle: the outcome depends on the starting point, and only one special line leads to equilibrium. Always state the conclusion in the context of the question, using the eigenvalues as the reason.
Exam tip

Give the reason, then the conclusion: 'both eigenvalues are negative, so both populations tend to zero'.

That's the notes covered.

Carry on to the next subtopic.

Exam questions on 5.17 Phase portraits

  1. The populations xx and yy (in hundreds) of two species of algae in a pond, tt weeks after they are introduced, satisfy dxdt=4x+y\frac{dx}{dt}=4x+y and dydt=2x+3y\frac{dy}{dt}=2x+3y. A GDC may be used.
    Find an eigenvector of the matrix corresponding to the eigenvalue 22.2 marks
  2. The deviations xx and yy (in thousands) of the populations of two interacting species from their equilibrium values satisfy dxdt=x+2y\frac{dx}{dt}=x+2y and dydt=3x\frac{dy}{dt}=3x. A GDC may be used.
    Find the equations of the two straight-line trajectories that pass through the origin.2 marks
  3. The deviations xx and yy (in hundreds) of two interacting populations from their equilibrium values, tt years after a survey begins, satisfy dxdt=−x−4y\frac{dx}{dt}=-x-4y and dydt=x−y\frac{dy}{dt}=x-y.
    Find the eigenvalues of the matrix (−1−41−1)\begin{pmatrix} -1 & -4 \\ 1 & -1 \end{pmatrix} that represents the system.3 marks
See the full worksheet

Written by the Exaim team, led by Shaun Daswani (Head of Upper Secondary, Improve ME Institute; MSc Financial Mathematics, Imperial College London; BSc, UCL) and Jason Daswani (operational lead, Improve ME Institute; LSE).