Obyčejná diferenciální rovnice

Obyčejné diferenciální rovnice jsou matematické rovnice, které obsahují neznámou funkci jedné nezávislé proměnné a její derivace. Název „obyčejné“ se používá jako protiklad k termínu parciální diferenciální rovnice, ve kterých se vyskytuje více než jedna nezávislá proměnná.

Nejjednodušší třídou obyčejných diferenciálních rovnic jsou lineární diferenciální rovnice, které mají strukturu množiny řešení pevnou, danou linearitou.

Řešení diferenciálních rovnic často nelze vyjádřit v analytickém tvaru pomocí elementárních funkcí. Někdy je možno tato řešení nalézt ve tvaru řady nebo integrálu. Pro řešení těchto diferenciálních rovnic se používají grafické nebo numerické metody, které řešení aproximují a mohou dávat užitečné informace, například o kvalitativním chování řešení.

Pozadí

Dráhu dělové střely lze popsat křivkou určenou obyčejnou diferenciální rovnicí, která je odvozena z druhého Newtonova zákona.

Obyčejné diferenciální rovnice se přirozeně objevují v nejrůznějších souvislostech v matematice i dalších vědách (přírodních i společenských), nejčastěji prostřednictvím rychlostí změn vyjádřených derivacemi. Pokud vyjádříme vztahy mezi funkcemi, jejich diferenciály a derivacemi pomocí rovnic, dostaneme diferenciální rovnice popisující změnu, vývoj nebo dynamiku jevů. Někdy se diferenciální rovnice objeví při práci s veličinami definovanými jako rychlost změny jiných veličin (časová derivace) nebo jejich gradienty (spády).

Diferenciální rovnice se používají v geometrii a analytické mechanice, ve fyzice a astronomii (nebeská mechanika), geologii (modelování počasí, strukturální geologie), chemii (dynamika reakcí)[1], biologii (infekční nemoci, genetické změny, systémová biologie), ekologii a modelování populace (konkurence populací, model dravce a kořisti), ekonomice (vývoj kurzů akcií, úrokových sazeb a tržní rovnováha změn ceny).

Diferenciální rovnice prvního řádu mají názornou geometrickou aplikaci díky směrovému poli, kdy je křivce znázorňující řešení rovnice v každém bodě předepsán sklon, závisející na souřadnících bodu, kterým křivka prochází.

Studiem diferenciálních rovnic se zabývalo mnoho matematiků, mezi jinými Isaac Newton, Gottfried Wilhelm Leibniz, Bernoulliovi, Jacopo Riccati, Alexis Claude Clairaut, Jean le Rond d'Alembert a Leonhard Euler, kteří přispěli k rozvoji tohoto oboru.

Jednoduchých příkladem diferenciální rovnice je Newtonův druhý pohybový zákon – vztah mezi polohou x hmotného bodu, na který působí síla F, a časem t vede k diferenciální rovnici

pro pohyb hmotného bodu s konstantní hmotností m. Pokud F závisí na poloze x(t) částice v čase t, neznámá funkce x(t) se objevuje na obou stranách diferenciální rovnice, jak je zřejmé ze zápisu F(x(t)).[2][3][4][5]

Definice

V následujícím textu je y závislá proměnná a x nezávislá proměnná, takže y = y(x) je neznámá funkce proměnné x. Různí autoři používají různé zápisy diferenciálu, podle toho, jaká notace nejlépe vyhovuje zadané úloze. Pro diferenciály se obvykle používá Leibnizova notace (dy/dx,d2y/dx2,...dny/dxn), zatímco pro provádění integrace je vhodnější Newtonova a Lagrangeova notace (y′,y′′, ... y(n)), kterou lze kompaktně zapsat derivace libovolného řádu.

Obecná definice obyčejné diferenciální rovnice

Obyčejnou diferenciální rovnici řádu n lze obecně zapsat ve tvaru:[6]

kde je hledané funkce.

Pokud lze rovnici zapsat ve tvaru

pak říkáme, že je rovnice rozřešena vzhledem k nejvyšší derivaci[7][8].

Existují další klasifikace diferenciálních rovnic:

  • Autonomní – Diferenciální rovnice, které neobsahují nezávislou proměnnou (tj. v předchozím odstavci x) se nazývají autonomní.
  • Lineární – Diferenciální rovnice se nazývá lineární, jestliže F lze zapsat jako lineární kombinaci zavislé proměnné a jejích derivací (v předchozím odstavci funkce y a její derivace):
kde ai(x) a r(x) jsou spojité funkce v x.[7][9][10] Funkce r(x) se nazývá pravá strana a z jejího tvaru vychází další klasifikace:[11][9]
  • Homogenní: Jestliže r(x) = 0 a následně jedno „automatické řešení“ je triviální řešení, y = 0. Řešení lineární homogenní rovnice je komplementární funkce, označovaná yc.
  • Nehomogenní: Jestliže r(x) ≠ 0. Přídavné řešení komplementární funkce je určitý integrál označovaný zde yp. Obecné řešení lineární rovnice lze vyjádřit jako y = yc + yp.
Nelineární rovnice je rovnice, která není lineární.

Soustava obyčejných diferenciálních rovnic

Několik spojených diferenciálních rovnic vytváří soustavu rovnic. Jestliže y je vektor, jehož prvky jsou funkce; y(x) = [y1(x), y2(x),..., ym(x)] a F je vektorová funkce proměnné y a jejích derivací, pak

je explicitní soustava obyčejných diferenciálních rovnic řádu nebo dimenze m. Ve formě sloupcového vektoru:

Funkce obecně nemusí být lineární. V implicitním tvaru lze soustavu obyčejných diferenciálních rovnic zapsat takto:

kde 0 = (0, 0,... 0) je nulový vektor. V maticovém tvaru

Řešení

Je-li dána diferenciální rovnice

Partikulární řešení je libovolná funkce u: I ⊂ ℝ → ℝ, jestliže má na I derivace až do řádu n a

Grafické znázornění partikulárního řešení se nazývá integrální křivka pro F.

Máme-li dvě řešení u: J ⊂ ℝ → ℝ a v: I ⊂ ℝ → ℝ, pak u se nazývá rozšířením v, jestliže IJ a

Řešení, které nemá žádné rozšíření, se nazývá maximální řešení. Řešení definované na celém ℝ se nazývá globální řešení.

Obecné řešení rovnice n-tého řádu obsahuje n nezávislých konstant.

Partikulární řešení lze získat z obecného řešení dosazením hodnot těchto konstant obvykle zvolených tak, aby byly splněny 'počáteční podmínky nebo hraniční podmínky'[12].

Singulární řešení je řešení, které nelze získat dosazením hodnot integračním konstantám v obecném řešení[13].

Teorie obyčejných diferenciálních rovnic

Singulární řešení

Teorie singulárních řešení obyčejných a parciálních diferenciálních rovnic byla předmětem výzkumu již od Leibnitzových dob, ale zvláštní pozornost jí byla věnována až od poloviny 19. století. Cennou, ale málo známou práci na toto téma napsal Houtain (1854). V roce 1873 se vůdčím duchem této teorie stal Jean Gaston Darboux, a svými geometrickými interpretacemi těchto řešení vytvořil obor, do něhož přispěli další matematici, především Felice Casorati a Arthur Cayley. Cayley vytvořil kolem roku 1872 teorii singulárních řešení diferenciálních rovnic první řádu, která byla okolo roku 1900 veřejně přijata.

Redukce na kvadratury

Primitivní pokus o řešení diferenciálních rovnic představovala snaha o redukci na kvadratury. Stejně jako se algebraici 18. století pokoušeli nalézt metodu pro řešení algebraických rovnic n-tého stupně, matematici, kteří prováděli první práce na poli diferenciálních rovnic, věřili, že existuje obecná metoda pro integraci libovolné diferenciální rovnice. Ale Carl Friedrich Gauss v roce 1799 ukázal, že pokud se nepoužijí komplexní čísla, je tato metoda velmi omezená. Matematici se proto začali věnovat studiu funkcí, čímž vznikl nový a velmi plodný obor. Prvním, kdo si uvědomoval důležitost tohoto přístupu, byl Augustin Louis Cauchy. Pak není třeba si klást otázku, zda je možné řešení pomocí známých funkcí nebo jejich integrálů, ale zda daná diferenciální rovnice postačuje pro definici funkce nezávislé proměnné nebo proměnných, a pokud ano, jaké jsou její charakteristické vlastnosti.

Fuchsova teorie

Související informace naleznete také v článku Frobeniova metoda.

Lazarus Fuchs svými dvěma monografiemi (Crelle, 1866, 1868) inspiroval nový přístup, který následně rozpracoval Thomé a Ferdinand Georg Frobenius. K teorii významně přispěl i Collet od roku 1869, i když svou metodu pro integraci nelineárních soustav popsal Bertrandovi již v roce 1868. Alfred Clebsch přispěl v roce 1873 teorií Abelovských integrálů. Protože je lze klasifikovat podle vlastnosti stěžejní křivky, která zůstává nezměněna při racionálních transformacích, Clebsch navrhl klasifikaci transcendentních funkcí definovaných diferenciálními rovnicemi podle invariantních vlastností odpovídajících povrchům při racionálních vzájemně jednoznačných transformacích.

Lieova teorie

Na uspokojivější základ postavil teorii diferenciálních rovnic Sophus Lie. Ve své práci z roku 1870 ukázal, že starší teorie integrace lze sjednotit využitím grup nyní nazývaných jeho jménem Lieovy grupy, a že obyčejné diferenciální rovnice, které připouštějí stejné infinitezimální transformace mají srovnatelně obtížné integrace. Lie také zdůraznil význam kontaktních transformací.

K uznání Lieovy grupová teorie diferenciálních rovnic vedly dva důvody:

  • sjednocuje mnoho známých ad hoc metod pro řešení diferenciálních rovnic
  • poskytuje výkonné nové způsoby pro hledání řešení.

Lieova teorie má aplikace pro obyčejné i parciální diferenciální rovnice.[14]

Obecný přístup k řešení diferenciálních rovnic využívá vlastnosti symetrie diferenciálních rovnic, spojité infinitezimální transformace jednoho řešení na jiné (Lieova teorie). Spojitá teorie grup, Lieovy algebry a diferenciální geometrie se používají pro porozumění struktuře lineární a nelineární (částečné) diferenciální rovnice pro generování integrovatelných rovnic, pro nalezení jejich Laxových dvojic, rekurzivních operátorů, Bäcklundovy transformace a nakonec nalezení přesného analytického řešení diferenciální rovnice.

Metody využívající symetrie se byly uznány pro studium diferenciálních rovnic vznikajících v matematice, fyzice, technice a mnoha jiný disciplínách.

Sturmova–Liouvilleova teorie

Související informace naleznete také v článku Sturmova–Liouvilleova teorie.

Sturmova–Liouvilleova teorie je teorií vlastních hodnot a vlastních funkcí lineárních operátorů definovaných pomocí homogenních lineárních rovnic druhého řádu a je užitečná pro analýzu určitých parciálních diferenciálních rovnic.

Existence a jednoznačnost řešení

Několik matematických vět určuje podmínky existence a jednoznačnosti řešení problémů počáteční hodnoty pro obyčejné diferenciální rovnice, lokálně i globálně. Dvě hlavní věty jsou:

VětaPředpokladZávěr
Peanova existenční větaF spojitápouze lokální existence
Picardova–Lindelöfova větaF Lipschitzovsky spojité zobrazenílokální existence a jednoznačnost

což jsou obojí lokální výsledky.

Zjednodušená věta o lokální existenci a jednoznačnosti

Věta může být jednoduše formulována následovně.[15] Pro rovnici a problém počáteční hodnoty:

jestliže F a ∂F/∂y jsou spojité v uzavřeném pravoúhelníku

v rovině x-y, kde a a b jsou reálná čísla (a, b ∈ ℝ) a × označuje kartézský součin, špičaté závorky znamenají uzavřený interval, pak existuje interval

pro nějaké h ∈ ℝ, na kterém lze nalézt řešení výše uvedené rovnice a problému počáteční hodnoty. Tj. existuje právě jedno řešení. Protože funkce F nemusí být lineární, tento vztah platí pro nelineární rovnice, které mají tvar F(x, y) a může být použito i na soustavy rovnic.

Globální jednoznačnost a maximální definiční obor řešení

Když jsou splněny podmínky Picardovy–Lindelöfovy věty, pak lokální existenci a jednoznačnost lze rozšířit na globální výsledek. Přesněji:[16]

Pro každou počáteční podmínku (x0, y0) existuje jediný maximální (může být i nekonečný) otevřený interval

takový, že libovolné řešení, které vyhovuje této počáteční podmínce je restrikcí řešení, které vyhovuje této počáteční podmínce s definičním oborem Imax.

V případě, že , existují právě dvě možnosti

  • exploze v konečném čase:
  • zachování definičního oboru:

kde Ω je otevřená množina, na níž je F definována, a je její hranice.

Všimněte si, že maximální definiční obor řešení

  • je vždy interval (aby byl jednoznačný)
  • může být menší než ℝ
  • může záviset na zvláštní volbě (x0, y0).
Příklad

To znamená, že F(x, y) = y2, která je C1 a tedy Lipschitzovsky spojitá pro všechna y splňující Picardovu–Lindelöfovu větu.

Ani v tomto jednoduchém případě nemůže být maximální definiční obor řešení celé ℝ, protože řešení je

která má maximální definiční obor:

To jasně ukazuje, že maximální interval může záviset na počátečních podmínkách. Bylo by možné brát definiční obor y jako , ale to by vedlo na definiční obor, která není intervalem, takže opačná hranice než ta, která je daná počáteční podmínkou, by nebyla spojena s počáteční podmínkou a proto by jí nebyla jednoznačně určena.

Maximální definiční obor není ℝ, protože

která je jedním ze dvou možných případů podle výše uvedené věty.

Redukce řádu

Diferenciální rovnici lze obvykle řešit snadněji, jestliže lze snížit její řád.

Redukce na soustavu prvního řádu

Jakákoli diferenciální rovnice řádu n,

může být zapsána jako soustava n diferenciálních rovnic prvního řádu definováním nové rodiny neznámých funkcí

pro i = 1, 2,... n. n-rozměrná soustava diferenciálních rovnic prvního řádu je pak

nebo kompaktněji ve vektorovém vyjádření:

kde

Přehled přesných řešení

Některé diferenciální rovnice mají řešení, které lze napsat v přesném a uzavřeném tvaru. Některé důležité třídy jsou uvedeny v této části.

V následující tabulce jsou P(x), Q(x), P(y), Q(y) a M(x,y), N(x,y) libovolné integrovatelné funkce x, y; b a c jsou dané reálné konstanty, a C1, C2,... jsou libovolné konstanty (obecně komplexní). Diferenciální rovnice jsou ve svém ekvivalentních a alternativních tvarech, které vedou k řešení pomocí integrace.

V integrálním řešení, λ a ε jsou nastrčené integrační proměnné (spojitá analogie indexů při sumaci) a notace ∫xF(λ)dλ znamená integrovat F(λ) vzhledem k λ, a po integraci provést substituci λ = x, bez přidávání konstant (explicitně uvedených).

Diferenciální rovniceMetoda řešeníObecné řešení
Separabilní rovnice
Prvního řádu, separabilní v x a y (obecný případ, speciální případy viz níže)[17]

Separace proměnných (vydělíme P2Q1).
Prvního řádu, separabilní v x[15]

Přímá integrace.
Prvního řádu, autonomní, separabilní v y[15]

Separace proměnných (vydělíme F).
Prvního řádu, separabilní v x a y[15]

Integrovat.
Obecná rovnice prvního řádu
Prvního řádu, homogenní[15]

Substituce y = ux, pak řešíme separací proměnných u a x.
Prvního řádu, separabilní[17]

Separace proměnných (vydělit xy).

Jestliže N = M, pak řešení je xy = C.

Exaktní diferenciální rovnice prvního řádu[15]

kde

Zintegrovat.

kde Y(y) a X(x) jsou funkce, které jsou ve vzorci místo integračních konstant, které vyjádříme tak, aby výsledná funkce F(x, y) vyhovovala počáteční rovnici.

Diferenciální rovnice prvního řádu, která není exaktní[15]

kde

Integrační faktor μ(x, y) vyhovující vztahu

Jestliže μ(x, y) je možné nalézt:

Obecná rovnice druhého řádu
Druhého řádu, autonomní[18]

Znásobíme rovnici výrazem 2dy/dx, provedeme substituci , pak dvakrát integrujeme.
Lineární rovnice (až po řád n)
Prvního řádu, lineární, nehomogenní, funkce koeficienty[15]

Integrační faktor: .
Druhého řádu, lineární, nehomogenní s konstantními koeficienty[19]

Komplementární funkce yc: předpokládáme yc = eαx, provedeme substituci a řešíme polynom v α, nalezneme lineárně nezávislé funkce .

Určitý integrál yp: obecně metodou variace konstant; pro velmi jednoduché r(x) můžeme řešení odhadnout[15].

Jestliže b2 > 4c, pak:

Jestliže b2 = 4c, pak:

Jestliže b2 < 4c, pak:

n-tého řádu, lineární, nehomogenní, s konstantními koeficienty[19]

Komplementární funkce yc: předpokládáme yc = eαx, provedeme substituci a řešíme polynom v α, nalezneme lineárně nezávislé funkce .

Určitý integrál yp: obecně metodou variace konstant, i když pro velmi jednoduché r(x) můžeme řešení odhadnout[15].

Protože αj jsou řešení polynomu stupně n: , pak:

pro αj vesměs různá,

pro každý kj-násobný kořen αj

pro některé αj komplexní, pak položíme α = χj + iγj a pomocí Eulerova vzorce lze některé termy v předchozích výsledcích zapsat ve tvaru

kde ϕj je libovolná konstanta (fázový posuv).

Software pro řešení obyčejných diferenciálních rovnic

  • FuncDesignerBSD licence, používá automatickou derivaci
  • Počítačový algebraický systém Maxima (GNU GPL)
  • COPASI – volně šiřitelný (Artistic License 2.0) softwarový balíček pro řešení a analýzu obyčejných diferenciálních rovnic
  • MATLAB – výpočetní software (MATrix LABoratory)
  • GNU Octave – vyšší programovací jazyk, primárně určený pro numerické výpočty
  • Scilab – software s otevřenými zdrojovými texty pro numerické výpočty
  • Maple
  • Mathematica – programovací jazyk společnosti Wolfram Research
  • Julia
  • SciPy – balíček pro jazyk Python obsahující modul pro řešení obyčejných diferenciálních rovnic.

Související články

Reference

V tomto článku byl použit překlad textu z článku Ordinary differential equation na anglické Wikipedii.

  • HALLIDAY, David; RESNICK, Robert. Physics. 3rd. vyd. New York: Wiley, 1977. ISBN 0-471-71716-9. 
  • HARPER, Charlie. Introduction to Mathematical Physics. New Jersey: Prentice-Hall, 1976. Dostupné online. ISBN 0-13-487538-9. 
  • KREYSZIG, Erwin. Advanced Engineering Mathematics. 3rd. vyd. New York: Wiley, 1972. Dostupné online. ISBN 0-471-50728-8. 
  • POLYANIN, D.; ZAITSEV, V. F. Handbook of Exact Solutions for Ordinary Differential Equations. 2. vyd. Boca Raton: Chapman & Hall/CRC Press, 2003. ISBN 1-58488-297-2. 
  • SIMMONS, George F. Differential Equations with Applications and Historical Notes. New York: McGraw-Hill, 1972. Dostupné online. 
  • TIPLER, Paul. Physics for Scientists and Engineers: Extended version. 3rd. vyd. New York: Worth Publishers, 1991. Dostupné online. ISBN 0-87901-432-6. 
  • BOSCAIN, Ugo; CHITOUR, Yacine. Introduction à l'automatique. [s.l.]: [s.n.], 2011. Dostupné online. (francouzsky) 
  • LAWRENCE, Dresner. Applications of Lie's Theory of Ordinary and Partial Differential Equations. Bristol and Philadelphia: Institute of Physics Publishing, 1999. 
  1. Mathematics for Chemists, D.M. Hirst, Macmillan Press, 1976, (Žádný ISBN) SBN: 333-18172-7
  2. Kreyszig 1972, strana 64
  3. Simmons 1972, strana 1-2
  4. Halliday 1977, strana 78
  5. Tipler 1991, strana 78–83
  6. Simmons 1972, strana 3
  7. a b Harper 1976, strana 127
  8. Kreyszig 1972, strana 2
  9. a b Kreyszig 1972, strana 24
  10. Simmons 1972, strana 47
  11. Harper 1976, strana 128
  12. Kreyszig 1972, strana 78
  13. Kreyszig 1972, strana 4
  14. Lawrence 1999, strana 9
  15. a b c d e f g h i j Elementary Differential Equations and Boundary Value Problems (4th Edition), W.E. Boyce, R.C. Diprima, Wiley International, John Wiley & Sons, 1986, ISBN 0-471-83824-1
  16. Boscain; Chitour 2011, p. 21
  17. a b Mathematical Handbook of Formulas and Tables (3rd edition), S. Lipschutz, M.R. Spiegel, J. Liu, Schuam's Outline Series, 2009, ISC_2N 978-0-07-154855-7
  18. Further Elementary Analysis, R. Porter, G.Bell & Sons (London), 1978, ISBN 0-7135-1594-5
  19. a b Mathematical methods for physics and engineering, K.F. Riley, M.P. Hobson, S.J. Bence, Cambridge University Press, 2010, ISC_2N 978-0-521-86153-3

Literatura

Literatura v češtině:

  • KURZWEIL, Jaroslav. Obyčejné diferenciální rovnice. Praha: SNTL - Nakladatelství technické literatury, 1978. Dostupné online. 

Literatura v angličtině:

  • CODDINGTON, Earl; LEVINSON, Norman. Theory of Ordinary Differential Equations. New York: McGraw-Hill, 1955. Dostupné online. 
  • HARTMAN, Philip. Theory of Ordinary Differential Equations. 2. vyd. [s.l.]: Society for Industrial & Applied Math, 2002. ISBN 0-89871-510-5. 
  • W. Johnson, A Treatise on Ordinary and Partial Differential Equations, John Wiley a Sons, 1913, v University of Michigan Historical Math Collection
  • INCE, E. L. Ordinary Differential Equations. [s.l.]: Dover Publications, 1958. ISBN 0-486-60349-0. 
  • HUREWICZ, Witold. Lectures on Ordinary Differential Equation. [s.l.]: Dover Publications, 2002. Dostupné online. ISBN 0-486-49510-8. 
  • IBRAGIMOV, Nail H. CRC Handbook of Lie Group Analysis of Differential Equations Vol. 1-3. Providence: CRC-Press, 1993. ISBN 0-8493-4488-3. .
  • TESCHL, Gerald. Ordinary Differential Equations and Dynamical Systems. [s.l.]: American Matematický Society, 2012. Dostupné online. ISBN 978-0-8218-8328-0. [nedostupný zdroj]
  • POLYANIN, D.; ZAITSEV, V. F.; MOUSSIAUX. Handbook of First Order Partial Differential Equations. London: Taylor & Francis, 2002. ISBN 0-415-27267-X. 
  • ZWILLINGER, D. Handbook of Differential Equations. 3. vyd. Boston: Academic Press, 1997. 

Externí odkazy

Média použitá na této stránce

Parabolic trajectory.svg
Illustration of a parabolic trajectory.