|
|
| Line 1: |
Line 1: |
| '''Linear dynamical systems''' are [[dynamical systems]] whose evaluation functions are [[linear]]. While dynamical systems in general do not have [[closed-form expression|closed-form solutions]], linear dynamical systems can be solved exactly, and they have a rich set of mathematical properties. Linear systems can also be used to understand the qualitative behavior of general dynamical systems, by calculating the equilibrium points of the system and approximating it as a linear system around each such point.
| | Greetings! I am Marvella and I really feel comfy when people use the full title. I am a meter reader but I plan on changing it. The thing she adores most is body building and now she is attempting to make cash with it. Her family life in Minnesota.<br><br>Also visit my web page ... [http://fastrolls.com/index.php?do=/profile-72113/info/ at home std testing] |
| | |
| ==Introduction==
| |
| | |
| In a linear dynamical system, the variation of a state vector
| |
| (an <math>N</math>-dimensional [[vector space|vector]] denoted <math>\mathbf{x}</math>) equals a constant matrix
| |
| (denoted <math>\mathbf{A}</math>) multiplied by
| |
| <math>\mathbf{x}</math>. This variation can take two forms: either
| |
| as a [[flow (mathematics)|flow]], in which <math>\mathbf{x}</math> varies
| |
| continuously with time
| |
| | |
| :<math>
| |
| \frac{d}{dt} \mathbf{x}(t) = \mathbf{A} \cdot \mathbf{x}(t)
| |
| </math>
| |
| | |
| or as a mapping, in which
| |
| <math>\mathbf{x}</math> varies in [[discrete time|discrete]] steps
| |
| | |
| :<math>
| |
| \mathbf{x}_{m+1} = \mathbf{A} \cdot \mathbf{x}_{m}
| |
| </math>
| |
| | |
| These equations are linear in the following sense: if
| |
| <math>\mathbf{x}(t)</math> and <math>\mathbf{y}(t)</math>
| |
| are two valid solutions, then so is any [[linear combination]]
| |
| of the two solutions, e.g.,
| |
| <math>\mathbf{z}(t) \ \stackrel{\mathrm{def}}{=}\ \alpha \mathbf{x}(t) + \beta \mathbf{y}(t)</math>
| |
| where <math>\alpha</math> and <math>\beta</math>
| |
| are any two [[scalar (mathematics)|scalars]]. The matrix <math>\mathbf{A}</math>
| |
| need not be [[Symmetry in mathematics#Symmetry in linear algebra|symmetric]].
| |
| | |
| Linear dynamical systems can be solved exactly, in contrast to most nonlinear ones. Occasionally, a nonlinear system can be solved exactly by a change of variables to a linear system. Moreover, the solutions of (almost) any nonlinear system can be well-approximated by an equivalent linear system near its [[fixed point (mathematics)|fixed points]]. Hence, understanding linear systems and their solutions is a crucial first step to understanding the more complex nonlinear systems.
| |
| | |
| ==Solution of linear dynamical systems==
| |
| | |
| If the initial vector <math>\mathbf{x}_{0} \ \stackrel{\mathrm{def}}{=}\ \mathbf{x}(t=0)</math>
| |
| is aligned with a [[right eigenvector]] <math>\mathbf{r}_{k}</math> of
| |
| the [[matrix (mathematics)|matrix]] <math>\mathbf{A}</math>, the dynamics are simple
| |
| | |
| :<math>
| |
| \frac{d}{dt} \mathbf{x}(t) =
| |
| \mathbf{A} \cdot \mathbf{r}_{k} = \lambda_{k} \mathbf{r}_{k}
| |
| </math>
| |
| | |
| where <math>\lambda_{k}</math> is the corresponding [[eigenvalue]];
| |
| the solution of this equation is
| |
| :<math>
| |
| \mathbf{x}(t) =
| |
| \mathbf{r}_{k} e^{\lambda_{k} t}
| |
| </math>
| |
| as may be confirmed by substitution.
| |
| | |
| If <math>\mathbf{A}</math> is [[diagonalizable matrix|diagonalizable]], then any vector in an <math>N</math>-dimensional space can be represented by a linear combination of the right and [[left eigenvector]]s (denoted <math>\mathbf{l}_{k}</math>) of the matrix <math>\mathbf{A}</math>.
| |
| | |
| :<math>
| |
| \mathbf{x}_{0} =
| |
| \sum_{k=1}^{N}
| |
| \left( \mathbf{l}_{k} \cdot \mathbf{x}_{0} \right)
| |
| \mathbf{r}_{k}
| |
| </math>
| |
| | |
| Therefore, the general solution for <math>\mathbf{x}(t)</math> is
| |
| a linear combination of the individual solutions for the right
| |
| eigenvectors
| |
| :<math>
| |
| \mathbf{x}(t) =
| |
| \sum_{k=1}^{n}
| |
| \left( \mathbf{l}_{k} \cdot \mathbf{x}_{0} \right)
| |
| \mathbf{r}_{k} e^{\lambda_{k} t}
| |
| </math>
| |
| | |
| Similar considerations apply to the discrete mappings.
| |
| | |
| == Classification in two dimensions ==
| |
| [[Image:LinDynSysTraceDet.jpg|right|thumb|400px|Classification of 2D fixed point according to the trace and the determinant of the Jacobian matrix.]]
| |
| The roots of the [[characteristic polynomial]] det('''A''' - λ'''I''') are the eigenvalues of '''A'''. The sign and relation of these roots, <math>\lambda_n</math>, to each other may be used to determine the stability of the dynamical system
| |
| :<math> | |
| \frac{d}{dt} \mathbf{x}(t) = \mathbf{A} \mathbf{x}(t).
| |
| </math>
| |
| For a 2-dimensional system, the characteristic polynomial is of the form <math>\lambda^2-\tau\lambda+\Delta=0</math> where <math>\tau</math> is the [[trace (linear algebra)|trace]] and <math>\Delta</math> is the [[determinant]] of '''A'''. Thus the two roots are in the form:
| |
| :<math>\lambda_1=\frac{\tau+\sqrt{\tau^2-4\Delta}}{2}</math>
| |
| :<math>\lambda_2=\frac{\tau-\sqrt{\tau^2-4\Delta}}{2}</math>
| |
| Note also that <math>\Delta=\lambda_1\lambda_2</math> and <math>\tau=\lambda_1+\lambda_2</math>. Thus if <math>\Delta<0</math> then the eigenvalues are of opposite sign, and the fixed point is a saddle. If <math>\Delta>0</math> then the eigenvalues are of the same sign. Therefore if <math>\tau>0</math> both are positive and the point is unstable, and if <math>\tau<0</math> then both are negative and the point is stable. The [[discriminant]] will tell you if the point is nodal or spiral (i.e. if the eigenvalues are real or complex).
| |
| | |
| <!--
| |
| ==Linear systems and higher-order differential equations==
| |
| | |
| ==Linearizing nonlinear systems==
| |
| -->
| |
| | |
| ==See also==
| |
| | |
| * [[Linear system]]
| |
| * [[Dynamic systems]]
| |
| * [[List of dynamical system topics]]
| |
| <!--
| |
| ==External links==
| |
| -->
| |
| | |
| [[Category:Dynamical systems]]
| |
Greetings! I am Marvella and I really feel comfy when people use the full title. I am a meter reader but I plan on changing it. The thing she adores most is body building and now she is attempting to make cash with it. Her family life in Minnesota.
Also visit my web page ... at home std testing