Phase Portraits (DP IB Applications & Interpretation (AI)): Revision Note

Did this video help you?

Phase Portraits

What is a phase portrait for a system of coupled differential equations?

  • Here we are again considering systems of coupled equations that can be represented in the matrix form bold italic x with bold dot on top equals bold italic M bold italic x, where bold italic x with bold dot on top equals open parentheses table row cell x with dot on top end cell row cell y with dot on top end cell end table close parentheses, bold italic M equals open parentheses table row a b row c d end table close parentheses, and bold italic x equals open parentheses table row x row y end table close parentheses

  • A phase portrait is a diagram showing how the values of x and y change over time

    • On a phase portrait we will usually sketch several typical solution trajectories

    • The precise trajectory that the solution for a particular system will travel along is determined by the initial conditions for the system

  • Let lambda subscript 1 and lambda subscript 2 be the eigenvalues of the matrix bold italic M

    • The overall nature of the phase portrait depends in large part on the values of lambda subscript 1 and lambda subscript 2

What does the phase portrait look like when lambda subscript 1 and lambda subscript 2 are real numbers?

  • Recall that for real distinct eigenvalues the solution to a system of the above form is bold italic x equals A straight e to the power of lambda subscript 1 t end exponent bold italic p subscript 1 bold plus B straight e to the power of lambda subscript 2 t end exponent bold italic p subscript bold 2, where lambda subscript 1 and lambda subscript 2 are the eigenvalues of bold italic M and bold italic p subscript 1 and bold italic p subscript 2 are the corresponding eigenvectors

    • AI HL only considers cases where lambda subscript 1 and lambda subscript 2 are distinct (i.e., lambda subscript 1 not equal to lambda subscript 2) and non-zero

  • A phase portrait will always include two ‘eigenvector lines’ through the origin, each one parallel to one of the eigenvectors

    • So if bold italic p subscript 1 equals open parentheses table row 1 row 2 end table close parentheses and bold italic p subscript 2 equals open parentheses table row cell negative 3 end cell row 4 end table close parentheses, for example, then these lines through the origin will have equations y equals 2 x and y equals negative 4 over 3 x, respectively

    • These lines will define two sets of solution trajectories

    • If the eigenvalue corresponding to a line’s eigenvector is positive, then there will be solution trajectories along the line away from the origin in both directions as t increases

    • If the eigenvalue corresponding to a line’s eigenvector is negative, then there will be solution trajectories along the line towards the origin in both directions as t increases

    • No solution trajectory will ever cross an eigenvector line

  • If both eigenvalues are positive then all solution trajectories will be directed away from the origin as t increases

    • In between the ‘eigenvector lines’ the trajectories as they move away from the origin will all curve to become approximately parallel to the line whose eigenvector corresponds to the larger eigenvalue

5-7-1-ib-ai-hl-real-pos-pos-eigenvalues-1
  • If both eigenvalues are negative then all solution trajectories will be directed towards the origin as t increases

    • In between the ‘eigenvector lines’ the trajectories will all curve so that at points further away from the origin they are approximately parallel to the line whose eigenvector corresponds to the more negative eigenvalue

      • They will then converge on the other eigenvalue line as they move in towards the origin

5-7-1-ib-ai-hl-real-neg-neg-eigenvalues-1
  • If one eigenvalue is positive and one eigenvalue is negative then solution trajectories will generally start by heading in towards the origin before curving to head out away again from the origin as t increases

    • In between the ‘eigenvector lines’ the solution trajectories will all move in towards the origin along the direction of the eigenvector line that corresponds to the negative eigenvalue, before curving away and converging on the eigenvector line that corresponds to the positive eigenvalue as they head away from the origin

5-7-1-ib-ai-hl-real-pos-neg-eigenvalues

What does the phase portrait look like when lambda subscript 1 and lambda subscript 2 are imaginary numbers?

  • Here the solution trajectories will all be either circles or ellipses with their centres at the origin

5-7-1-ib-ai-hl-imaginary-eigenvalues
  • You can determine the direction (clockwise or anticlockwise) and the shape (circular or elliptical) of the trajectories by considering the values of fraction numerator d x over denominator d t end fraction and fraction numerator d y over denominator d t end fraction for points on the coordinate axes

    • For example, consider the system bold italic x with bold dot on top equals open parentheses table row 1 cell negative 2 end cell row 1 cell negative 1 end cell end table close parentheses bold italic x

      • The eigenvalues of open parentheses table row 1 cell negative 2 end cell row 1 cell negative 1 end cell end table close parentheses are straight i and negative straight i, so the trajectories will be elliptical or circular

    • When ­x equals 1 and y equals 0, fraction numerator d x over denominator d t end fraction equals 1 left parenthesis 1 right parenthesis minus 2 left parenthesis 0 right parenthesis equals 1 and fraction numerator d y over denominator d t end fraction equals 1 left parenthesis 1 right parenthesis minus 1 left parenthesis 0 right parenthesis equals 1

      • This shows that from a point on the positive x-axis the solution trajectory will be moving ‘to the right and up’ in the direction of the vector open parentheses table row 1 row 1 end table close parentheses

    • When ­x equals 0 and y equals 1, fraction numerator d x over denominator d t end fraction equals 1 left parenthesis 0 right parenthesis minus 2 left parenthesis 1 right parenthesis equals negative 2 and fraction numerator d y over denominator d t end fraction equals 1 left parenthesis 0 right parenthesis minus 1 left parenthesis 1 right parenthesis equals negative 1

      • This shows that from a point on the positive y-axis the solution trajectory will be moving ‘to the left and down’ in the direction of the vector open parentheses table row cell negative 2 end cell row cell negative 1 end cell end table close parentheses

    • The directions of the trajectories at those points tell us that the directions of the trajectories will be anticlockwise

    • They also tell us that the trajectories will be ellipses

      • For circular trajectories, the direction of the trajectories when they cross a coordinate axis will be perpendicular to that coordinate axis

What does the phase portrait look like when lambda subscript 1 and lambda subscript 2 are complex numbers?

  • In this case lambda subscript 1 and lambda subscript 2 will be complex conjugates of the form a plus-or-minus b straight i, where a and b are non-zero real numbers

    • If a equals 0, b not equal to 0, then we have the imaginary eigenvalues case above

  • Here the solution trajectories will all be spirals

    • If the real part of the eigenvalues is positive (i.e., if a greater than 0), then the trajectories will spiral away from the origin

    • If the real part of the eigenvalues is negative (i.e., if a less than 0), then the trajectories will spiral towards the origin

5-7-1-ib-ai-hl-complex-eigenvalues
  • You can determine the direction (clockwise or anticlockwise) of the trajectories  by considering the values of fraction numerator d x over denominator d t end fraction and fraction numerator d y over denominator d t end fraction for points on the coordinate axes

    • For example, consider the system bold italic x with dot on top equals open parentheses table row 1 5 row cell negative 2 end cell 3 end table close parentheses bold italic x

      • The eigenvalues of open parentheses table row 1 cell negative 5 end cell row 2 3 end table close parentheses are 2 plus 3 straight i and 2 minus 3 straight i, so the trajectories will be spirals

      • Because the real part of the eigenvalues left parenthesis 2 right parenthesis is positive, the trajectories will spiral away from the origin

    • When ­x equals 1 and y equals 0, fraction numerator d x over denominator d t end fraction equals 1 left parenthesis 1 right parenthesis plus 5 left parenthesis 0 right parenthesis equals 1 and fraction numerator d y over denominator d t end fraction equals negative 2 left parenthesis 1 right parenthesis plus 3 left parenthesis 0 right parenthesis equals negative 2

      • This shows that from a point on the positive x-axis the solution trajectory will be moving ‘to the right and down’ in the direction of the vector open parentheses table row 1 row cell negative 2 end cell end table close parentheses

    • When ­x equals 0 and y equals 1, fraction numerator d x over denominator d t end fraction equals 1 left parenthesis 0 right parenthesis plus 5 left parenthesis 1 right parenthesis equals 5 and fraction numerator d t over denominator d t end fraction equals negative 2 left parenthesis 0 right parenthesis plus 3 left parenthesis 1 right parenthesis equals 3

      • This shows that from a point on the positive y-axis the solution trajectory will be moving ‘to the right and up’ in the direction of the vector open parentheses table row 5 row 3 end table close parentheses

    • The directions of the trajectories at those points tell us that the directions of the trajectory spirals will be clockwise

Worked Example

Consider the system of coupled differential equations

 fraction numerator d x over denominator d t end fraction equals negative 2 x plus 2 y
fraction numerator d y over denominator d t end fraction equals x minus 3 y

Given that negative 1 and negative 4 are the eigenvalues of the matrix open parentheses table row cell negative 2 end cell 2 row 1 cell negative 3 end cell end table close parentheses, with corresponding eigenvectors open parentheses table row 2 row 1 end table close parentheses and open parentheses table row cell negative 1 end cell row 1 end table close parentheses, draw a phase portrait for the solutions of the system.

5-7-1-ib-ai-hl-phase-portraits-we-solution

You've read 0 of your 5 free revision notes this week

Unlock more, it's free!

Join the 100,000+ Students that ❤️ Save My Exams

the (exam) results speak for themselves:

Did this page help you?

Roger B

Author: Roger B

Expertise: Maths Content Creator

Roger's teaching experience stretches all the way back to 1992, and in that time he has taught students at all levels between Year 7 and university undergraduate. Having conducted and published postgraduate research into the mathematical theory behind quantum computing, he is more than confident in dealing with mathematics at any level the exam boards might throw at you.