Inženýrská matematika | Robert Mařík © 2007-2012 |
Poznámka 3.1 (klasifikace stacionárních bodů). Podle chování trajektorií v okolí stacionárních bodů rozdělujeme tyto stacionární body do několika navzájem disjunktních skupin. Nechť \( \displaystyle (x^{*},y^{*})\) je singulárním bodem systému (1.1).
Definice 3.1 (Jacobiho matice). Matice
\[ J(x,y) = \left (\array{ { \partial\, f \over \partial\, x} (x,y)& { \partial\, f \over \partial\, y} (x,y) \cr { \partial\, g \over \partial\, x} (x,y)& { \partial\, g \over \partial\, y} (x,y) \cr } \right ) \] |
se nazývá Jacobiho matice soustavy (1.1).
Definice 3.2. Charakteristickou rovnicí matice \( \displaystyle A\) rozumíme kvadratickou rovnici \( \displaystyle \det (A -\lambda I) = 0\) s neznámou \( \displaystyle \lambda \), tj. charakteristickou rovnicí matice \( \displaystyle A = \left (\array{ a& b\cr c& d} \right )\) je rovnice
\[ \lambda ^{2} - (a + d)\lambda + ad - bc = 0. \] |
Kořeny této rovnice (reálné nebo komplexní) nazýváme vlastní čísla matice \( \displaystyle A\).
Věta 3.1 (klasifikace stacionárních bodů pomocí vlastních čísel Jacobiho matice). Předpokládejme, že nula není vlastním číslem Jacobiho matice pro dvourozměrný autonomní systém (1.1).
Vlastní hodnoty, \( \displaystyle \lambda _{1,2}\in \mathbb{R}\) | typ stac. bodu | průběh trajektorií |
\( \displaystyle \lambda _{1}\geq \lambda _{2} > 0\) | nestabilní uzel | |
\( \displaystyle \lambda _{1} > 0 >\lambda _{2}\) | sedlo | |
\( \displaystyle 0 >\lambda _{1}\geq \lambda _{2}\) | stabilní uzel | |
Vlastní hodnoty, \( \displaystyle \lambda _{1,2}\not \in \mathbb{R}\) | typ stac. bodu | průběh trajektorií |
\( \displaystyle ℜ(\lambda _{1,2}) > 0\) | nestabilní ohnisko | |
\( \displaystyle ℜ(\lambda _{1,2}) < 0\) | stabilní ohnisko | |
\( \displaystyle ℜ(\lambda _{1,2}) = 0\) | ohnisko nebo bod rotace | nebo kterákoliv z předchozích možností |
Poznámka 3.2. Zjednodušeně řečeno, kdykoliv se mezi vlastními hodnotami Jacobiho matice v bodě \( \displaystyle S\) objeví vlastní hodnota se zápornou reálnou částí, existuje trajektorie konvergující do bodu \( \displaystyle S\). Pokud má některé vlastní hodnota kladnou reálnou část, existuje trajektorie vycházející z bodu \( \displaystyle S\). Pokud mají vlastní hodnoty nenulovou imaginární část, dochází v okolí bodu \( \displaystyle S\) k oscilacím.
Jiná možnost, jak určit typ stacionárních bodů, je obsažena v následující větě. V této větě \( \displaystyle D\) značí determinant Jacobiho matice v bodě \( \displaystyle (x^{*},y^{*})\), tj. a \( \displaystyle \Delta \) stopu1 Jacobiho matice v tomto bodě, tj.
\[ \begin{align*} D & =\det J(x^{*},y^{*}) ={ \partial\, f \over \partial\, x} (x^{*},y^{*}){ \partial\, g \over \partial\, y} (x^{*},y^{*}) -{ \partial\, f \over \partial\, y} (x^{*},y^{*}){ \partial\, g \over \partial\, x} (x^{*},y^{*}), & & \\\Delta & =\mathop{\mathrm{Tr}} J(x^{*},y^{*}) ={ \partial\, f \over \partial\, x} (x^{*},y^{*}) +{ \partial\, g \over \partial\, y} (x^{*},y^{*}). & & \\\end{align*}\]Věta 3.2. Nechť \( \displaystyle (x^{*},y^{*})\) je stacionární bod soustavy (1.1) a \( \displaystyle J(x^{*},y^{*})\) hodnota Jacobiho matice v tomto bodě. Pomocí čísel \( \displaystyle D\) a \( \displaystyle \Delta \) lze rozhodnout o kvalitě stacionárního bodu \( \displaystyle (x^{*},y^{*})\) podle následující tabulky.
\( \displaystyle D^{} < 0\) | sedlo | ||
\( \displaystyle D > 0\) | \( \displaystyle \Delta > 0\) | \( \displaystyle \Delta ^{2}\geq 4D\) | nestabilní uzel |
\( \displaystyle \Delta ^{2} < 4D\) | nestabilní ohnisko | ||
\( \displaystyle \Delta < 0\) | \( \displaystyle \Delta ^{2}\geq 4D\) | stabilní uzel | |
\( \displaystyle \Delta ^{2} < 4D\) | stabilní ohnisko | ||
\( \displaystyle \Delta ^{} = 0\) | bod rotace nebo ohnisko | ||
Robert Mařík, Ústav matematiky, Lesnická a dřevařská fakulta Mendelovy univerzity v Brně | © 2007-2012 |