Fundamentální systém

Fundamentální systém soustavy homogenních lineárních obyčejných diferenciálních rovnic se v matematické analýze nazývá každá báze vektorového prostoru složená z řešení této soustavy.

Fundamentálním systémem soustavy homogenních diferenciálních rovnic nazveme každou množinu vektorových funkcí \{\mathbf{y}_1,\ldots,\mathbf{y}_n\} takovou, že :\mathcal{L} := \{\mathbf{y} \in C^1(\langle a,b\rangle; \mathbb{R}^n)\ |\ \mathbf{y} = \sum_{k=1}^na_k\mathbf{y}_k\ ,\ a_1, \ldots, a_n \in \mathbb{R}\} je množina všech řešení této soustavy.

Znalost fundamentálního systému je předpokladem pro použití metody variace konstant k získání partikuláního řešení nehomogenní lineární diferenciální rovnice prvního řádu i nehomogenních lineárních diferenciálních rovnic vyšších řádů.

Fundamentální systém, (Hlavní) fundamentální matice a Wronskián

Homogenní lineární soustava diferenciálních rovnic prvního řádu

Uvažujme lineární homogenní soustavu diferenciálních rovnic prvního řádu :\mathbf{y}^\prime(x)=A(x)\,\mathbf{y}(x), na nějakém intervalu I \subseteq \mathbb{R}, kde A:\langle a,b\rangle\to\R^{n\times n} je maticové zobrazení. Hledáme řešení této soustavy diferenciálních rovnic v prostoru C^1(\langle a,b\rangle;\R^n) spojitě diferencovatelných funkcí \mathbf{y}:\langle a,b\rangle\to\R^{n}.

Jestliže máme dvě různá řešení této diferenciální rovnice, pak jejich součet a násobek s reálnou konstantou jsou také řešením této rovnice. Prostor řešení je tedy vlastním podprostorem prostoru všech spojitě diferencovatelných funkcí.

Jestliže matice koeficientů A je spojitá maticová funkce, pak lze na ni použít Picardovu-Lindelöfovu větu o existenci a jednoznačnosti. Pak každé řešení diferenciální rovnice je jednoznačně určeno hodnotou \mathbf{y}(a) v počátečních podmínkách a naopak každá počáteční hodnota \ \mathbf{y}(a)=:\mathbf{y}_0\in\R^n této soustavy diferenciálních rovnic určuje jednoznačně řešení. +more Z toho plyne, že prostor řešení je n-dimenzionální.

Definice

Každá báze tohoto n-dimenzionálního prostoru řešení se nazývá fundamentální systém lineární soustavy diferenciálních rovnic. Ve většině případů si volíme bázi soustavy funkcí, které jsou řešeními \{\mathbf{y}_1(x), \ldots, \mathbf{y}_n(x)\}, pro které počáteční hodnota \mathbf{y}_i(a)=\mathbf{e}_i je i-tý kanonický jednotkový vektor.

Pokud \{\mathbf{y}_1, \ldots, \mathbf{y}_n\} je fundamentální systém, pak matici \Phi(x) := (\mathbf{y}_1(x)\ |\ \cdots\ |\ \mathbf{y}_n(x)) \in \mathbb{R}^{n\times n} nazveme fundamentální maticí a její determinant \ \det \Phi(x) budeme nazývat Wronskián. Jestliže \Phi(x_0) pro nějaké x_0 je jednotková matice, pak se \Phi nazývá hlavní fundamentální matice v bodě x_0.

Fundamentální matice \Phi je také řešení homogenní obyčejné (maticové) diferenciální rovnice, konkrétně :\Phi^\prime(x)=A(x)\Phi(x). Prostor řešení původní homogenní soustavy v \mathbb{R}^n pak je \{\mathbf{y} \in C^1(\langle a,b\rangle; \mathbb{R}^n)\ |\ \mathbf{y}(x)=\Phi(x)\cdot \mathbf{c},\ \mathbf{c} \in \mathbb{R}^n\}. +more \Phi je dokonce i hlavní fundamentální matice pro x_0, takže \mathbf{y}(x) := \Phi(x)\mathbf{y}_0 řeší počáteční úlohu pro \mathbf{y}(x_0) = \mathbf{y}_0.

Fundamentální matice \Phi(x)\in\R^{n\times n} je pro každé x\in\langle a,b\rangle invertibilní. Pro Wronskián platí Liouvilleův vzorec.

Homogenní lineární diferenciální rovnice vyššího řádu

Stejně jako v případě rovnic prvního řádu může být prostor řešení lineární soustavy vyššího řádu také vektorový prostor, a proto každou jeho bázi lze považovat za fundamentální systém.

Pro definici fundamentální matice skalární lineární diferenciální rovnice řádu n :y^{(n)}(x) = \sum_{k=0}^{n-1}a_k(x)y^{(k)}(x) nejdříve uvažujeme diferenciální rovnici odpovídající soustavě n rovnic prvního řádu :\mathbf{Y}'(x) = A(x)\mathbf{Y}(x), kde A(x) := \begin{pmatrix} 0&1&&0\\ &\ddots&\ddots&\\ &&\ddots&1\\ a_0(x)&a_1(x)&\cdots&a_{n-1}(x)\\ \end{pmatrix} . [Vztah je následující: y(x) je řešením skalární rovnice řádu n právě tehdy když \mathbf{Y}(x) := \begin{pmatrix} y(x)\\y'(x)\\ \vdots\\ y^{(n-1)}(x)\\ \end{pmatrix} je řešením výše uvedené soustavy prvního řádu. +more].

Fundamentální maticí rovnice :y^{(n)}(x) = \sum_{k=0}^{n-1}a_k(x)y^{(k)}(x) nazýváme každou fundamentální matici \Phi soustavy prvního řádu :\mathbf{Y}'(x) = A(x)\mathbf{Y}(x)\ . Přirozeně \Phi nazýváme hlavní fundamentální maticí v x_0, jestliže \Phi(x_0) je jednotková matice. +more Determinant \det\Phi se nazývá Wronskián.

Pro redukci rovnice na soustavu prvního řádu: jestliže \{y_1,\ldots,y_n\} je fundamentální systém, pak :\Phi(x) := \begin{pmatrix} y_1(x)&\cdots&y_n(x)\\ y_1'(x)&\cdots&y_n'(x)\\ \vdots&\cdots&\vdots\\ y_1^{(n-1)}(x)&\cdots&y_n^{(n-1)}(x)\\ \end{pmatrix} je fundamentální matice.

Konstrukce fundamentálního systému

Zkonstruovat fundamentální systém je v obecném případě obtížné. Přesné postupy jsou známy pouze pro diferenciální rovnice se speciální strukturou, jako jsou skalární diferenciální rovnice prvního řádu, soustavy diferenciálních rovnic prvního řádu s konstantními koeficienty, diferenciální rovnice vyšších řádů s konstantními koeficienty a Eulerova rovnice. +more Pokud je známé řešení homogenní diferenciální rovnice vyššího řádu, lze snížit řád diferenciální rovnice použitím postupu pro redukci řádu objeveného Jean le Rond d'Alembertem.

Lineární diferenciální rovnice prvního řádu

Nechť A je primitivní funkce k funkci a, potom :\ \{y(x) = e^{A(x)}\} je fundamentální systém rovnice y'(x) = a(x)y(x).

Lineární soustava diferenciálních rovnic prvního řádu s konstantními koeficienty

Pro soustavu lineárních diferenciálních rovnic s konstantními koeficienty :\ \mathbf{y}'(x) = A\cdot \mathbf{y}(x)\ ,\ A \in \mathbb{R}^{n\times n} nejdříve nalezneme Jordanovu normální formu J matice A a příslušnou Jordanovu bázi B = \{\mathbf{b}_1,\ldots,\mathbf{b}_n\}. Pokud \lambda je komplexní vlastní hodnota s vlastními vektory báze \mathbf{c}_1,\ldots,\mathbf{c}_k, můžeme v Jordanově bázi vybrat bázové vektory tak, že \overline{\mathbf{c}_1},\ldots,\overline{\mathbf{c}_k} se stanou bázovými vektory \overline{\lambda}.

Nyní provedeme pro každý řetězec hlavních vektorů zvlášť: Pokud \mathbf{v}_1, \ldots, \mathbf{v}_k \in B je (úplný) řetězec hlavních vektorů pro vlastní číslo \lambda, tj. :\ (A-\lambda I)\mathbf{v}_{i+1} = \mathbf{v}_i, pak dodávají do fundamentálního systému k řešení (hlavního vektoru) :y_1(x) = e^{\lambda x}\mathbf{v}_1\ ,\ y_2(x) = e^{\lambda x}\left(\frac{x^1}{1. +more}\mathbf{v}_1 + \mathbf{v}_2\right)\ ,\ y_3(x) = e^{\lambda x}\left(\frac{x^2}{2. }\mathbf{v}_1 + \frac{x^1}{1. }\mathbf{v}_2 + \mathbf{v}_3\right)\ ,\ \ldots\ , obecně :y_i(x) = e^{\lambda x} \sum_{j=1}^i{\frac {x^{i-j}}{(i-j). }\mathbf{v}_j}\ ,\ i = 1, \ldots, k Probráním všech řetězců hlavních vektorů sestrojíme (obecně komplexní) fundamentální systém.

Lineární diferenciální rovnice vyššího řádu s konstantními koeficienty

Fundamentální systém pro skalární lineární diferenciální rovnici n-tého řádu s konstantními koeficienty : y^{(n)}(x) - \sum_{k=0}^{n-1}a_ky^{(k)}(x) = 0\ ,\ a_1, \ldots, a_n \in \mathbb{R} získáme řešením charakteristické rovnice P(\lambda) = 0 s charakteristickým polynomem :P(\lambda) := \lambda^n - \sum_{k=0}^{n-1}a_k\lambda^k. Jestliže \lambda_1,\ldots,\lambda_k jsou (po dvou různé) kořeny polynomu P s násobnostmi \mu_1,\ldots,\mu_k, pak (komplexní) fundamentální systém \mu_i lineárně nezávislých řešení pro kořen \lambda_i je :y_{i,1}(x) = e^{\lambda_i x}\ ,\ y_{i,2}(x) = xe^{\lambda_i x}\ ,\ \ldots\ ,\ y_{i,\mu_i}(x) = x^{\mu_i-1}e^{\lambda_i x}.

[To vysvětluje způsob vyjadřování: pomocí výše uvedené transformace získáme skalární rovnici n-tého řádu ze soustavy diferenciálních rovnic prvního řádu, takže máme matici koeficientů jako charakteristický polynom přesně toto, jak bylo zde ukázáno.]

Reálný fundamentální systém

Výše uvedeným postupem získáme n lineárně nezávislých řešení, která mohou obsahovat komplexní čísla. Díky tomu, že pokud má charakteristické rovnice s reálnými koeficienty komplexní kořeny, jsou to dvojice komplexně sdružených čísel, využitím linearity diferenciální rovnice a Eulerova vzorce e^{ix} = \cos x + i\sin x lze komplexní fundamentální systém převést na reálný tak, že každou dvojici komplexně sdružených řešení y(x), \overline{y}(x) z (komplexního) fundamentálního systému převedeme na reálná řešení {\rm Re\;}y(x), {\rm Im\;}y(x).

Periodická soustava diferenciálních rovnic prvního řádu

Pro soustavu :\ \mathbf{y}'(x) = A(x)\mathbf{y}(x) s \omega-periodickou spojitou maticí koeficientů A: \mathbb{R} \rightarrow \mathbb{R}^{m \times m} nelze explicitně zkonstruovat fundamentální systém - ale díky Floquetově teorii lze zjistit jaká bude struktura fundamentální matice této soustavy.

Příklady

Lineární soustava diferenciálních rovnic prvního řádu s konstantními koeficienty

Řešíme následující soustavu diferenciálních rovnic :\mathbf{y}'(x) = A\cdot \mathbf{y}(x)\ ,\ A := \begin{pmatrix}3&-1&1\\2&0&1\\1&-1&2\\\end{pmatrix}. Matice A má jednoduché vlastní číslo 1 a dvojité vlastní číslo 2. +more Prostor vlastních vektorů je E(A,1) = \langle\begin{pmatrix}0\\1\\1\\\end{pmatrix}\rangle\ ,\ E(A,2) = \langle\begin{pmatrix}1\\1\\0\\\end{pmatrix}\rangle. Pro hlavní řetěz vektorů pro vlastní číslo 2 je stále potřeba :\textrm{Kern}(A-2I)^2 = \langle\begin{pmatrix}1\\1\\0\\\end{pmatrix}, \begin{pmatrix}0\\0\\1\\\end{pmatrix}\rangle. Zvolíme například :\mathbf{v}_2 := \begin{pmatrix}0\\0\\2\\\end{pmatrix} \in \textrm{Kern}(A-2I)^2 \setminus \textrm{Kern}(A-2I). Pak musíme jako hlavní vektor v první fázi zvolit \mathbf{v}_1 := (A-2I)\mathbf{v}_2 = \begin{pmatrix}2\\2\\0\\\end{pmatrix}. Tím dostaneme fundamentální systém \{\mathbf{y}_1, \mathbf{y}_2, \mathbf{y}_3\}, kde :\mathbf{y}_1(x) := e^x\cdot\begin{pmatrix}0\\1\\1\\\end{pmatrix}\ ,\ \mathbf{y}_2(x) := e^{2x}\cdot\begin{pmatrix}2\\2\\0\\\end{pmatrix}\ ,\ \mathbf{y}_3(x) := e^{2x}\cdot\left[\begin{pmatrix}2x\\2x\\0\\\end{pmatrix} + \begin{pmatrix}0\\0\\2\\\end{pmatrix}\right]\ .

Lineární diferenciální rovnice vyššího řádu s konstantními koeficienty

Uvažujme diferenciální rovnici :y^{(4)}(x) - y(x) = 0. Tato rovnice má charakteristický polynom \lambda^4-1, který má čtyři kořeny 1, -1, i, -i. +more Odtud hned získáme komplexní fundamentální systém :\{e^x, e^{-x}, e^{ix}, e^{-ix}\}, kterému odpovídá reálný fundamentální systém :\{e^x, e^{-x}, \sin x, \cos x\}.

Reference

Literatura

Carmen Chicone: Ordinary Differential Equations with Applications. 2. +more Auflage. Texts in Applied Mathematics 34, Springer-Verlag 2006, . * Harro Heuser: Gewöhnliche Differentialgleichungen. Teubner, 1995. S. 250. *.

fundamentální systém