NonLinear Lecture 1.2
NonLinear Lecture 1.2
Plane
Lecture 1.2
Paul White
Outline
• “Quotes” and preamble.
• State-space
• Phase plots
• Trajectories and phase portraits
• Equilibrium points
• Solutions of linear equations
2
Phase Plane
• Phase plane methods explore the behaviour of non-linear
systems in a graphical manner.
• The phase plane is a plot of x1 against x2.
• A solution of the equations defines x1(t) and x2(t), each point in
time specifies a point in the phase plane and as time evolves
that point moves along a curve within the phase plane.
• Curves that are drawn represent a solution – these curves are
called trajectories (a.k.a. phase paths or integral curves).
• By studying the nature of the trajectories we seek to
understand the behaviour of the system.
3
Example
Solution of mass-spring-damper
system: x1 is displacement, x2 is
velocity, initial condition (1,0).
4
Phase Portrait
• A plot of a range of trajectories for a particular equation in
some sense characterises all the solutions for that system.
– Note that the trajectories are curves with arrows indicating the direction
of travel of the solution.
• The trajectories are plotted for a range of different initial
conditions.
• The phase portrait is a valuable tool when trying to
“understand” the behaviour of a particular system.
• We seek methods to sketch phase portraits without explicitly
solving the system of equations.
5
Many similar solutions
plotted for different initial
conditions on the circle
x12+x22=1.
6
Finding Trajectories
• The trajectories for a system can in principle be computed
directly.
dx1
= f1 ( x1 , x2 )
dt dx1 f1 ( x1 , x2 )
⇒ =
dx2 dx2 f 2 ( x1 , x2 )
= f 2 ( x1 , x2 )
dt
• Solving this equation defines x1 as a function of x2, which
defines a trajectory.
• However, this equation is rarely any more tractable than the
original system of equations.
• Phase plane methods seek to plot trajectories using without
resorting to solving the above equation directly.
7
Example
• Undamped oscillator:
x + ω2 x = 0
x1 = x2
x2 x =
= x1 x ⇒
x2 = −ω2 x1
∫ ∫
dx1 x
=− 22 ⇒ ω2 x1dx1 =− x2 dx2
dx 2 ω x1
⇒ ω2 x12 + x2 2 =C
dx1 f1 ( x1 , x2 )
=
dx2 f 2 ( x1 , x2 )