Cesko.wiki Metoda neurčitých koeficientů
Author
Albert FloresMetoda neurčitých koeficientů je v matematice přístup k hledání partikulárního řešení určitých nehomogenních obyčejných diferenciálních a diferenčních rovnic. Metoda je blízce příbuzná metodě anihilátorů, ale místo použití určitého druhu diferenciálního operátoru (anihilátoru) pro nalezení nejlepšího možného tvaru partikulárního řešení se provede „odhad“, vhodného tvaru řešení, který se pak upřesní a ověří derivováním výsledné rovnice. Pro složité rovnice je rychlejší použít metodu anihilátorů nebo variace parametrů.
Metoda neurčitých koeficientů není tak obecná jako metoda variace konstant, protože je použitelná pouze pro diferenciální rovnice, které mají určité tvary.
Popis metody
Uvažujeme lineární nehomogenní obyčejnou diferenciální rovnici tvaru
:a_ny^{(n)} + a_{(n-1)}y^{(n-1)}+...+a_1y' + a_0y = g(x).
Metoda spočívá v hledání obecného homogenního řešení y_c komplementární lineární homogenní diferenciální rovnice
:a_ny^{(n)} + a_{(n-1)}y^{(n-1)}+...+a_1y' + a_0y = 0,
a partikulárního integrálu y_p lineární nehomogenní obyčejné diferenciální rovnice s pravou stranou g(x). Obecné řešení y lineární nehomogenní obyčejné diferenciální rovnice pak je
:y = y_c + y_p.
Pokud g(x) vyjádříme jako součet dvou funkcí h(x) + w(x), pak říkáme, že y_{p_1} je řešení vycházející z h(x) a y_{p_2} řešení vycházející z w(x). Pak použitím principu superpozice můžeme říct, že partikulární integrál y_p je
:y_p = y_{p_1} + y_{p_2}.
Typické tvary partikulárního integrálu
Pro nalezení partikulárního integrálu potřebujeme „uhodnout“ jeho tvar, přičemž některé koeficienty ponecháme proměnné, a jejich hodnoty zjistíme vyřešením rovnice. Tvarem je první derivace komplementární funkce. +more Následuje tabulka některých typických funkcí a odhadů tvaru jejich řešení.
| Funkce x | Tvar y |
|---|---|
| k e^{a x}\. | C e^{a x}\. +more |
| k x^n \mathrm{,}\; n = 0, 1, 2,\cdots\. | K_n x^n + K_{n-1} x^{n-1} + \cdots + K_1 x + K_0 \. |
| k \cos(a x) \;\;\mathrm{nebo}\;\;k \sin(a x) \. | K \cos(a x) + M \sin(a x) \. |
| k e^{a x} \cos(b x) \;\;\mathrm{nebo}\;\; ke^{a x} \sin(b x) \. | e^{a x} (K \cos(b x) + M \sin(b x)) \. |
| k e^{a x} \cosh(b x) \;\;\mathrm{nebo}\;\; ke^{a x} \sinh(b x) \. | e^{a x} (K \cosh(b x) + M \sinh(b x)) \. |
| P_n(x) e^{a x} \cos(b x) \;\;\mathrm{nebo}\;\; P_n(x) e^{a x} \sin(b x)\. | e^{a x} \left(\left(\sum_{i=0}^n Q_i x^i\right) \cos(b x) + \left(\sum_{i=0}^n R_i x^i\right) \sin(b x)\right) |
kde P_n(x) = \sum_{i=1}^n k_i x^i\!
Pokud se v homogenním řešení objeví nějaký člen z výše uvedeného partikulárního integrálu pro y, musíme jej vynásobit dostatečnou mocninou x, aby řešení nebyla závislá. Jestliže funkci proměnné x lze vyjádřit jako součet členů z výše uvedené tabulky, jako odhad partikulárního integrálu použijeme součet odpovídajících termů proměnné y.
Příklady
;Příklad 1
Hledáme partikulární integrál rovnice
:y + y = t \cos {t}. \!
Pravá strana t cos t má tvar
:P_n e^{\alpha t} \cos{\beta t} \!
pro n=1, k1=1, α=0 a β=1.
Protože α + iβ = i je jednoduchý kořen charakteristické rovnice
:\lambda^2 + 1 = 0 \!
budeme zkoušet partikulární integrál tvaru
: \begin{align} y_p &= t [F_1 (t) e^{\alpha t} \cos{\beta t} + G_1 (t) e^{\alpha t} \sin{\beta t}] \\ &= t [F_1 (t) \cos{t} + G_1 (t) \sin{t}]\\ &= t [(A_0 t + A_1) \cos{t} + (B_0 t + B_1) \sin{t}] \\ &= (A_0 t^2 + A_1 t) \cos{t} + (B_0 t^2 + B_1 t) \sin{t} .\\ \end{align}
Substitucí yp do diferenciální rovnice, dostaneme identitu
: \begin{align}t \cos{t} &= y_p + y_p \\ &= [(A_0 t^2 + A_1 t) \cos{t} + (B_0 t^2 + B_1 t) \sin{t}] \\ &\quad + [(A_0 t^2 + A_1 t) \cos{t} + (B_0 t^2 + B_1 t) \sin{t}] \\ &= [2A_0 \cos{t} + 2(2A_0 t + A_1)(- \sin{t}) + (A_0 t^2 + A_1 t)(- \cos{t})] \\ &\quad +[2B_0 \sin{t} + 2(2B_0 t + B_1) \cos{t} + (B_0 t^2 + B_1 t)(- \sin{t})] \\ &\quad +[(A_0 t^2 + A_1 t) \cos{t} + (B_0 t^2 + B_1 t) \sin{t}] \\ &= [4B_0 t + (2A_0 + 2B_1)] \cos{t} + [-4A_0 t + (-2A_1 + 2B_0)] \sin{t}. \\ \end{align}
Porovnáním obou stran dostaneme soustavu
\begin{array}{rrrrl} &&4B_0&&=1\\ 2A_0 &&& + 2B_1 &= 0 \\ -4A_0 &&&& = 0 \\ &-2A_1 &+ 2B_0 && = 0 \\ \end{array}
která má řešení A_0 = 0, A_1 = 1/4, B_0 = 1/4, B_1 = 0. Výsledný partikulární integrál je tedy
:y_p = \frac {1} {4} t \cos{t} + \frac {1} {4} t^2 \sin{t}.
Příklad 2
Hledáme řešení lineární nehomogenní diferenciální rovnice :\frac{\mathrm{d}y}{\mathrm{d}x} = y + e^x. Nehomogenní část má tvar uvedený v prvním případě v #Typické tvary partikulárního integrálu|předchozí kapitole pro a=1. +more Protože tato nehomogenní část (e^x) je lineárně závislá s obecným řešení homogenní části (c_1 e^x), musíme uvedený odhad znásobit dostatečně velkou mocninou x, aby byl lineárně nezávislý.
Výsledný odhad je: :y_p = A x e^x.
substitucí této funkce a její derivace do diferenciální rovnice můžeme zjistit řešení A: :\frac{\mathrm{d}}{\mathrm{d}x} \left( A x e^x \right) = A x e^x + e^x :A x e^x + A e^x = A x e^x + e^x :A = 1.
Takže obecné řešení původní diferenciální rovnice je:
:y = c_1 e^x + xe^x.
Příklad 3
Hledáme obecné řešení rovnice: :\frac{\mathrm{d}y}{\mathrm{d}t} = t^2 - y f(t) = t^2 je polynom 2. stupně, takže hledáme řešení jako obecný polynom druhého stupně: :y_p = A t^2 + B t + C, tedy :\frac{\mathrm{d} y_p}{\mathrm{d}t} = 2 A t + B Dosazením tohoto partikulárního integrálu s konstantami A, B a C do původní rovnice dostaneme :2 A t + B = t^2 - (A t^2 + B t + C), odtud :t^2 - A t^2 = 0 :-B t = 2 A t :-C = B tj. +more :A = 1 :B = -2 :C = 2 Po dosazení vypočítaných konstant máme :y_p = t^2 - 2 t + 2 Obecné řešení je :y= y_p + y_c kde y_c je homogenní řešení y_c = c_1 e^{-t}. Obecné řešení tedy je: :y= t^2 - 2 t + 2 + c_1 e^{-t}.
Reference
Elementary Differential Equations and Boundary Value Problems (4th Edition), W. E. +more Boyce, R. C. Diprima, Wiley International, John Wiley & Sons, 1986, * Mathematical methods for physics and engineering, K. F. Riley, M. P. Hobson, S. J. Bence, Cambridge University Press, 2010,.