# Phase plane

Differential equations :
Scope :
Fields
Natural sciencesEngineering
Astronomy
Physics
Chemistry
Biology
Geology
Applied mathematics
Continuum mechanics
Chaos theory
Dynamical systems
Social sciences
Economics
Population dynamics
List of named differential equations :
Classification :
Types
Ordinary
Partial
Differential-algebraic
Integro-differential
Fractional
Linear
Non-linear
By variable type
Dependent and independent variables
Autonomous
Coupled / Decoupled
Exact
Homogeneous / Nonhomogeneous
Features
Order
Operator
Notation :
Relation to processes
Difference (discrete analogue)
Stochastic
Stochastic partial
Delay :
Solution :
Existence and uniqueness
Picard–Lindelöf theorem
Peano existence theorem
Carathéodory's existence theorem
Cauchy–Kowalevski theorem :
General topics
Initial conditions
Boundary values
Dirichlet
Neumann
Robin
Cauchy problem
Wronskian
Phase portrait
Lyapunov / Asymptotic / Exponential stability
Rate of convergence
Series / Integral solutions
Numerical integration
Dirac delta function :
Solution methods
Inspection
Method of characteristics
Euler
Exponential response formula
Finite difference (Crank–Nicolson)
Finite element
Infinite element
Finite volume
Galerkin
Petrov–Galerkin
Green's function
Integrating factor
Integral transforms
Perturbation theory
Runge–Kutta
Separation of variables
Undetermined coefficients
Variation of parameters :
People :
List
Isaac Newton
Gottfried Leibniz
Jacob Bernoulli
Leonhard Euler
Joseph-Louis Lagrange
Józef Maria Hoene-Wroński
Joseph Fourier
Augustin-Louis Cauchy
George Green
Carl David Tolmé Runge
Martin Kutta
Rudolf Lipschitz
Ernst Lindelöf
Émile Picard
Phyllis Nicolson
John Crank :
vte :
In applied mathematics, in particular the context of nonlinear system analysis, a phase plane is a visual display of certain characteristics of certain kinds of differential equations; a coordinate plane with axes being the values of the two state variables, say (x, y), or (q, p) etc. (any pair of variables). It is a two-dimensional case of the general n-dimensional phase space.
The phase plane method refers to graphically determining the existence of limit cycles in the solutions of the differential equation.
The solutions to the differential equation are a family of functions. Graphically, this can be plotted in the phase plane like a two-dimensional vector field. Vectors representing the derivatives of the points with respect to a parameter (say time t), that is (dx/dt, dy/dt), at representative points are drawn. With enough of these arrows in place the system behaviour over the regions of plane in analysis can be visualized and limit cycles can be easily identified.
The entire field is the phase portrait, a particular path taken along a flow line (i.e. a path always tangent to the vectors) is a phase path. The flows in the vector field indicate the time-evolution of the system the differential equation describes.
In this way, phase planes are useful in visualizing the behaviour of physical systems; in particular, of oscillatory systems such as predator-prey models (see Lotka–Volterra equations). In these models the phase paths can "spiral in" towards zero, "spiral out" towards infinity, or reach neutrally stable situations called centres where the path traced out can be either circular, elliptical, or ovoid, or some variant thereof. This is useful in determining if the dynamics are stable or not.
Other examples of oscillatory systems are certain chemical reactions with multiple steps, some of which involve dynamic equilibria rather than reactions that go to completion. In such cases one can model the rise and fall of reactant and product concentration (or mass, or amount of substance) with the correct differential equations and a good understanding of chemical kinetics.

## Example of a linear system

A two-dimensional system of linear differential equations can be written in the form:
    ${\begin{aligned}{\frac {dx}{dt}}&=Ax+By\\{\frac {dy}{dt}}&=Cx+Dy\end{aligned}}$
which can be organized into a matrix equation:
    ${\begin{aligned}&{\frac {d}{dt}}{\begin{bmatrix}x\\y\\\end{bmatrix}}={\begin{bmatrix}A&B\\C&D\\\end{bmatrix}}{\begin{bmatrix}x\\y\\\end{bmatrix}}\\&{\frac {d\mathbf {v} }{dt}}=\mathbf {A} \mathbf {v} .\end{aligned}}$
where A is the 2 × 2 coefficient matrix above, and v = (x, y) is a coordinate vector of two independent variables.
Such systems may be solved analytically, for this case by integrating:
${\frac {dy}{dx}}={\frac {Cx+Dy}{Ax+By}}$
although the solutions are implicit functions in x and y, and are difficult to interpret.

### Solving using eigenvalues

More commonly they are solved with the coefficients of the right hand side written in matrix form using eigenvalues λ, given by the determinant:
    $\det \left({\begin{bmatrix}A&B\\C&D\\\end{bmatrix}}-\lambda \mathbf {I} \right)=0$
and eigenvectors:
    ${\begin{bmatrix}A&B\\C&D\\\end{bmatrix}}\mathbf {x} =\lambda \mathbf {x}$
The eigenvalues represent the powers of the exponential components and the eigenvectors are coefficients. If the solutions are written in algebraic form, they express the fundamental multiplicative factor of the exponential term. Due to the nonuniqueness of eigenvectors, every solution arrived at in this way has undetermined constants c$_1$, c$_2$, …, c$_n$.
The general solution is:
    $x={\begin{bmatrix}k_{1}\\k_{2}\end{bmatrix}}c_{1}e^{\lambda _{1}t}+{\begin{bmatrix}k_{3}\\k_{4}\end{bmatrix}}c_{2}e^{\lambda _{2}t}.$
where λ$_1$ and λ$_2$ are the eigenvalues, and (k$_1$, k$_2$), (k$_3$, k$_4$) are the basic eigenvectors. The constants c$_1$ and c$_2$ account for the nonuniqueness of eigenvectors and are not solvable unless an initial condition is given for the system.
The above determinant leads to the characteristic polynomial:
    $\lambda ^{2}-(A+D)\lambda +(AD-BC)=0$
which is just a quadratic equation of the form:
    $\lambda ^{2}-p\lambda +q=0$
where $p=A+D=\mathrm {tr} (\mathbf {A} )\,,$ ("tr" denotes trace) and $q=AD-BC=\det(\mathbf {A} )\,.$
The explicit solution of the eigenvalues are then given by the quadratic formula:
    $\lambda ={\frac {1}{2}}(p\pm {\sqrt {\Delta }})\,$
where $\Delta =p^{2}-4q\,.$

### Eigenvectors and nodes

The eigenvectors and nodes determine the profile of the phase paths, providing a pictorial interpretation of the solution to the dynamical system, as shown next.
The phase plane is then first set-up by drawing straight lines representing the two eigenvectors (which represent stable situations where the system either converges towards those lines or diverges away from them). Then the phase plane is plotted by using full lines instead of direction field dashes. The signs of the eigenvalues indicate the phase plane's behaviour:
* If the signs are opposite, the intersection of the eigenvectors is a saddle point.
* If the signs are both positive, the eigenvectors represent stable situations that the system diverges away from, and the intersection is an unstable node.
* If the signs are both negative, the eigenvectors represent stable situations that the system converges towards, and the intersection is a stable node.
The above can be visualized by recalling the behaviour of exponential terms in differential equation solutions.

### Repeated eigenvalues

This example covers only the case for real, separate eigenvalues. Real, repeated eigenvalues require solving the coefficient matrix with an unknown vector and the first eigenvector to generate the second solution of a two-by-two system. However, if the matrix is symmetric, it is possible to use the orthogonal eigenvector to generate the second solution.

### Complex eigenvalues

Complex eigenvalues and eigenvectors generate solutions in the form of sines and cosines as well as exponentials. One of the simplicities in this situation is that only one of the eigenvalues and one of the eigenvectors is needed to generate the full solution set for the system.