|
Consider the general system of two first-order ordinary differential equations
Let and denote fixed points with , so
Then expand about so
To first-order, this gives
![d/(dt)[deltax; deltay]=[f_x(x_0,y_0) f_y(x_0,y_0); g_x(x_0,y_0) g_y(x_0,y_0)][deltax; deltay],](/images/equations/LinearStability/NumberedEquation1.gif) |
(9)
|
where the matrix
is called the stability matrix.
In general, given an -dimensional map , let be a fixed point, so that
 |
(10)
|
Expand about the fixed point,
so
 |
(13)
|
The map can be transformed into the principal axis frame by finding the eigenvectors and eigenvalues
of the matrix
 |
(14)
|
so the determinant
 |
(15)
|
The mapping is
![deltax_(princ)^'=[lambda_1 ... 0; | ... |; 0 ... lambda_n].](/images/equations/LinearStability/NumberedEquation6.gif) |
(16)
|
When iterated a large number of times,
only if for all , but
if any . Analysis of the eigenvalues (and eigenvectors)
of therefore characterizes the type of fixed point.
Tabor, M. "Linear Stability Analysis." §1.4 in Chaos and Integrability in Nonlinear Dynamics: An Introduction.
New York: Wiley, pp. 20-31, 1989.
|