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í řešení, jež lze sčítat a násobit koeficienty. Řešení těchto rovnic lze obvykle zapsat v analytickém tvaru pomocí elementárních funkcí. Naproti tomu obecné obyčejné diferenciální rovnice nemají aditivní řešení, nejsou lineární, a jejich řešení je mnohem komplikovanější; často je není možné ani vyjádřit v analytickém tvaru pomocí elementárních funkcí. Někdy lze přesná a analytická řešení obyčejných diferenciálních rovnic nalézt ve tvaru řady nebo integrálu. Pro řešení těchto diferenciálních rovnic se při ručních výpočtech nebo na počítači často používají grafické nebo numerické metody, které řešení aproximují a mohou dávat užitečné informace, často dostačující při neexistenci přesných analytických ř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 objevují v nejrůznějších souvislostech v matematice i dalších vědách (přírodních i společenských): Jestliže změny popisujeme matematicky, nejpřesnější popis používá diferenciály a derivace. 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ů. Určité veličiny jsou často definovány jako rychlost změny jiných veličin (časová derivace) nebo jejich gradienty (spády), a jejich vzájemné vztahy lze přirozeně popsat diferenciálními rovnicemi.

Zvláštní matematický obor zahrnuje geometrii a analytickou mechaniku. Vědecké obory zahrnují větší část fyziky a astronomie (nebeská mechanika), geologie (modelování počasí), chemie (dynamika reakcí)[1], biologie (infekční nemoci, genetické změny), ekologie a modelování populace (soutěžení populací), ekonomiku (vývoj kurzů akcií, úrokových sezeb a tržní rovnováha změn ceny).

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

m{\frac {\mathrm {d} ^{2}x(t)}{\mathrm {d} t^{2}}}=F(x(t)),\,

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,d²y/dx²,...dⁿy/dxⁿ), zatímco pro provádění integrace je vhodnější Newtonova a Lagrangeova notace (y′,y′′, ... y⁽ⁿ⁾), 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]

F\left(x,y,y',y'',\ \cdots ,\ y^{(n)}\right)=0

kde $y=y(x)$ je hledané funkce.

Pokud lze rovnici zapsat ve tvaru

f\left(x,y,y',\cdots y^{(n-1)}\right)=y^{(n)}

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 nezávislé na x se nazývají autonomní.
Lineární – Diferenciální rovnice se nazývá lineární, jestliže F lze zapsat jako lineární kombinaci derivací y:

y^{(n)}=\sum _{i=0}^{n-1}a_{i}(x)y^{(i)}+r(x)

kde a_i(x) a r(x) jsou spojité funkce v x.[7][9][10]

Nelineární je taková rovnice, která nemůže být zapsaná v tomto tvaru. Funkce r(x) se nazývá pravá strana a je z jejího tvaru vychází další důležitá 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á y_c.

Nehomogenní: Jestliže r(x) ≠ 0. Přídavné řešení komplementární funkce je určitý integrál označovaný zde y_p. Obecné řešení lineární rovnice lze vyjádřit jako y = y_c + y_p.

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) = [y₁(x), y₂(x),..., y_m(x)] a F je vektorová funkce proměnné y a jejích derivací, pak

\mathbf {y} ^{(n)}=\mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n-1)}\right)

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

{\begin{pmatrix}y_{1}^{(n)}\\y_{2}^{(n)}\\\vdots \\y_{m}^{(n)}\end{pmatrix}}={\begin{pmatrix}F_{1}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n-1)}\right)\\F_{2}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n-1)}\right)\\\vdots \\F_{m}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n-1)}\right)\\\end{pmatrix}}

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

\mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n)}\right)={\boldsymbol {0}}

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

{\begin{pmatrix}F_{1}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n)})\\F_{2}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n)})\\\vdots \\F_{m}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\cdots \mathbf {y} ^{(n)})\\\end{pmatrix}}={\begin{pmatrix}0\\0\\\vdots \\0\\\end{pmatrix}}

Řešení

Je-li dána diferenciální rovnice

F\left(x,y,y',\cdots ,y^{(n)}\right)=0

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

F(x,u,u',\ \cdots ,\ u^{(n)})=0\quad x\in I.

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 I ⊂ J a

u(x)=v(x)\quad x\in I.\,

Ř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

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 $f=0$ 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

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ěta	Předpoklad	Závěr
Peanova existenční věta	F spojitá	pouze lokální existence
Picardova–Lindelöfova věta	F 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:

y'=F(x,y)\,,\quad y_{0}=y(x_{0})

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

R=\langle x_{0}-a,x_{0}+a\rangle \times \langle y_{0}-b,y_{0}+b\rangle

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

I=\langle x_{0}-h,x_{0}+h\rangle \subset \langle x_{0}-a,x_{0}+a\rangle

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 (x₀, y₀) existuje jediný maximální (může být i nekonečný) otevřený interval

I_{max}=(x_{-},x_{+}),x_{\pm }\in \mathbb {R} ,x_{0}\in I_{max}

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 I_max.

V případě, že $x_{\pm }\nrightarrow \pm \infty$ , existují právě dvě možnosti

exploze v konečném čase: $\lim _{x\to x_{\pm }}\|y(x)\|\rightarrow \infty$
zachování definičního oboru: $\lim _{x\to x_{\pm }}\in \partial {\bar {\Omega }}$

kde Ω je otevřená množina, na níž je F definována, a $\partial {\bar {\Omega }}$ 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ě (x₀, y₀).

Příklad

y'=y^{2}

To znamená, že F(x, y) = y², která je C¹ 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

y(x)={\frac {y_{0}}{(x_{0}-x)y_{0}+1}}

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

{\begin{cases}\mathbb {R} &y_{0}=0\\(-\infty ,x_{0}+{\frac {1}{y_{0}}})&y_{0}>0\\(x_{0}+{\frac {1}{y_{0}}},+\infty )&y_{0}<0\end{cases}}

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 $\mathbb {R} \smallsetminus (x_{0}+1/y_{0})$ , 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

\lim _{x\to x_{\pm }}\|y(x)\|\rightarrow \infty \,,

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,

F\left(x,y,y',y'',\ \cdots ,\ y^{(n-1)}\right)=y^{(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í

y_{i}=y^{(i-1)}.\!

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

{\begin{array}{rcl}y_{1}'&=&y_{2}\\y_{2}'&=&y_{3}\\&\vdots &\\y_{n-1}'&=&y_{n}\\y_{n}'&=&F(x,y_{1},\cdots ,y_{n}).\end{array}}

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

\mathbf {y} '=\mathbf {F} (x,\mathbf {y} )

kde

\mathbf {y} =(y_{1},\cdots ,y_{n}),\quad \mathbf {F} (x,y_{1},\cdots ,y_{n})=(y_{2},\cdots ,y_{n},F(x,y_{1},\cdots ,y_{n})).

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 C₁, C₂,... 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í rovnice	Metoda ř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] $P_{1}(x)Q_{1}(y)+P_{2}(x)Q_{2}(y)\,{\frac {\mathrm {d} y}{\mathrm {d} x}}=0\,\!$ $P_{1}(x)Q_{1}(y)\,\mathrm {d} x+P_{2}(x)Q_{2}(y)\,\mathrm {d} y=0\,\!$	Separace proměnných (vydělíme P₂Q₁).	$\int ^{x}{\frac {P_{1}(\lambda )}{P_{2}(\lambda )}}\,\mathrm {d} \lambda +\int ^{y}{\frac {Q_{2}(\lambda )}{Q_{1}(\lambda )}}\,\mathrm {d} \lambda =C\,\!$
Prvního řádu, separabilní v x[15] ${\frac {\mathrm {d} y}{\mathrm {d} x}}=F(x)\,\!$ $\mathrm {d} y=F(x)\,\mathrm {d} x\,\!$	Přímá integrace.	$y=\int ^{x}F(\lambda )\,\mathrm {d} \lambda +C\,\!$
Prvního řádu, autonomní, separabilní v y[15] ${\frac {\mathrm {d} y}{\mathrm {d} x}}=F(y)\,\!$ $\mathrm {d} y=F(y)\,\mathrm {d} x\,\!$	Separace proměnných (vydělíme F).	$x=\int ^{y}{\frac {\mathrm {d} \lambda }{F(\lambda )}}+C\,\!$
Prvního řádu, separabilní v x a y[15] $P(y){\frac {\mathrm {d} y}{\mathrm {d} x}}+Q(x)=0\,\!$ $P(y)\,\mathrm {d} y+Q(x)\,\mathrm {d} x=0\,\!$	Integrovat.	$\int ^{y}P(\lambda )\,{\mathrm {d} \lambda }+\int ^{x}Q(\lambda )\,\mathrm {d} \lambda =C\,\!$
Obecná rovnice prvního řádu
Prvního řádu, homogenní[15] ${\frac {\mathrm {d} y}{\mathrm {d} x}}=F\left({\frac {y}{x}}\right)\,\!$	Substituce y = ux, pak řešíme separací proměnných u a x.	$\ln(Cx)=\int ^{y/x}{\frac {\mathrm {d} \lambda }{F(\lambda )-\lambda }}\,\!$
Prvního řádu, separabilní[17] $yM(xy)+xN(xy)\,{\frac {\mathrm {d} y}{\mathrm {d} x}}=0\,\!$ $yM(xy)\,\mathrm {d} x+xN(xy)\,\mathrm {d} y=0\,\!$	Separace proměnných (vydělit xy).	$\ln(Cx)=\int ^{xy}{\frac {N(\lambda )\,\mathrm {d} \lambda }{\lambda [N(\lambda )-M(\lambda )]}}\,\!$ Jestliže N = M, pak řešení je xy = C.
Exaktní diferenciální rovnice prvního řádu[15] $M(x,y){\frac {\mathrm {d} y}{\mathrm {d} x}}+N(x,y)=0\,\!$ $M(x,y)\,\mathrm {d} y+N(x,y)\,\mathrm {d} x=0\,\!$ kde ${\frac {\partial M}{\partial x}}={\frac {\partial N}{\partial y}}\,\!$	Zintegrovat.	${\begin{aligned}F(x,y)&=\int ^{y}M(x,\lambda )\,\mathrm {d} \lambda +\int ^{x}N(\lambda ,y)\,\mathrm {d} \lambda \\&+Y(y)+X(x)=C\end{aligned}}\,\!$ 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] $M(x,y){\frac {\mathrm {d} y}{\mathrm {d} x}}+N(x,y)=0\,\!$ $M(x,y)\,\mathrm {d} y+N(x,y)\,\mathrm {d} x=0\,\!$ kde ${\frac {\partial M}{\partial x}}\neq {\frac {\partial N}{\partial y}}\,\!$	Integrační faktor μ(x, y) vyhovující vztahu ${\frac {\partial (\mu M)}{\partial x}}={\frac {\partial (\mu N)}{\partial y}}\,\!$	Jestliže μ(x, y) je možné nalézt: ${\begin{aligned}F(x,y)&=\int ^{y}\mu (x,\lambda )M(x,\lambda )\,\mathrm {d} \lambda +\int ^{x}\mu (\lambda ,y)N(\lambda ,y)\,\mathrm {d} \lambda \\&+Y(y)+X(x)=C\\\end{aligned}}\,\!$
Obecná rovnice druhého řádu
Druhého řádu, autonomní[18] ${\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}=F(y)\,\!$	Znásobíme rovnici výrazem 2dy/dx, provedeme substituci $2{\frac {\mathrm {d} y}{\mathrm {d} x}}{\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}={\frac {\mathrm {d} }{\mathrm {d} x}}\left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}\,\!$ , pak dvakrát integrujeme.	$x=\pm \int ^{y}{\frac {\mathrm {d} \lambda }{\sqrt {2\int ^{\lambda }F(\epsilon )\,\mathrm {d} \epsilon +C_{1}}}}+C_{2}\,\!$
Lineární rovnice (až po řád n)
Prvního řádu, lineární, nehomogenní, funkce koeficienty[15] ${\frac {\mathrm {d} y}{\mathrm {d} x}}+P(x)y=Q(x)\,\!$	Integrační faktor: $e^{\int ^{x}P(\lambda )\,\mathrm {d} \lambda }$ .	$y=e^{-\int ^{x}P(\lambda )\,\mathrm {d} \lambda }\left[\int ^{x}e^{\int ^{\lambda }P(\epsilon )\,\mathrm {d} \epsilon }Q(\lambda )\,{\mathrm {d} \lambda }+C\right]$
Druhého řádu, lineární, nehomogenní s konstantními koeficienty[19] ${\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+b{\frac {\mathrm {d} y}{\mathrm {d} x}}+cy=r(x)\,\!$	Komplementární funkce y_c: předpokládáme y_c = e^αx, provedeme substituci a řešíme polynom v α, nalezneme lineárně nezávislé funkce $e^{\alpha _{j}x}$ . Určitý integrál y_p: obecně metodou variace konstant; pro velmi jednoduché r(x) můžeme řešení odhadnout[15].	$y=y_{c}+y_{p}$ Jestliže b² > 4c, pak: $y_{c}=C_{1}e^{\left(-b+{\sqrt {b^{2}-4c}}\right){\frac {x}{2}}}+C_{2}e^{-\left(b+{\sqrt {b^{2}-4c}}\right){\frac {x}{2}}}\,\!$ Jestliže b² = 4c, pak: $y_{c}=(C_{1}x+C_{2})e^{-bx/2}\,\!$ Jestliže b² < 4c, pak: $y_{c}=e^{-b{\frac {x}{2}}}\left[C_{1}\sin {\left({\sqrt {\left\|b^{2}-4c\right\|}}{\frac {x}{2}}\right)}+C_{2}\cos {\left({\sqrt {\left\|b^{2}-4c\right\|}}{\frac {x}{2}}\right)}\right]\,\!$
n-tého řádu, lineární, nehomogenní, s konstantními koeficienty[19] $\sum _{j=0}^{n}b_{j}{\frac {\mathrm {d} ^{j}y}{\mathrm {d} x^{j}}}=r(x)\,\!$	Komplementární funkce y_c: předpokládáme y_c = e^αx, provedeme substituci a řešíme polynom v α, nalezneme lineárně nezávislé funkce $e^{\alpha _{j}x}$ . Určitý integrál y_p: obecně metodou variace konstant, i když pro velmi jednoduché r(x) můžeme řešení odhadnout[15].	$y=y_{c}+y_{p}$ Protože α_j jsou řešení polynomu stupně n: $\prod _{j=1}^{n}\left(\alpha -\alpha _{j}\right)=0\,\!$ , pak: pro α_j vesměs různá, $y_{c}=\sum _{j=1}^{n}C_{j}e^{\alpha _{j}x}\,\!$ pro každý k_j-násobný kořen α_j $y_{c}=\sum _{j=1}^{n}\left(\sum _{\ell =1}^{k_{j}}C_{\ell }x^{\ell -1}\right)e^{\alpha _{j}x}\,\!$ 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 $C_{j}e^{\alpha _{j}x}=C_{j}e^{\chi _{j}x}\cos(\gamma _{j}x+\phi _{j})\,\!$ kde ϕ_j je libovolná konstanta (fázový posuv).

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

FuncDesigner – BSD 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

Okrajová úloha
Diferenční rovnice
Laplaceova transformace aplikovaná na diferenciální rovnice
Seznam článků o dynamických systémech a diferenciálních rovnicích
Maticové diferenciální rovnice
Metoda neurčitých koeficientů
Numerické řešení obyčejných diferenciálních rovnic
Separace proměnných
Rungeova–Kuttova metoda

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.

Mathematics for Chemists, D.M. Hirst, Macmillan Press, 1976, (Žádný ISBN) SBN: 333-18172-7
Kreyszig 1972, strana 64
Simmons 1972, strana 1-2
Halliday 1977, strana 78
Tipler 1991, strana 78–83
Simmons 1972, strana 3
Harper 1976, strana 127
Kreyszig 1972, strana 2
Kreyszig 1972, strana 24
Simmons 1972, strana 47
Harper 1976, strana 128
Kreyszig 1972, strana 78
Kreyszig 1972, strana 4
Lawrence 1999, strana 9
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
Boscain; Chitour 2011, p. 21
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
Further Elementary Analysis, R. Porter, G.Bell & Sons (London), 1978, ISBN 0-7135-1594-5
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

Obrázky, zvuky či videa k tématu obyčejná diferenciální rovnice na Wikimedia Commons
Hazewinkel, Michiel, ed. (2001), "Differential equation, ordinary", Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4, id=p/d031910
Differential Equations on the Open Directory Project - Science/Math/Differential_Equations/|Differential Equations (obsahuje seznam programového vybavení pro řešení diferenciálních rovnic).
EqWorld: The World of Mathematical Equations, obsahující seznam obyčejných diferenciálních rovnic s řešeními.
Online Notes / Differential Equations autor Paul Dawkins, Lamar University.
Differential Equations, S.O.S. Mathematics.
A primer on analytical solution of differential equations from the Holistic Numerical Methods Institute, University of South Florida.
Ordinary Differential Equations and Dynamical Systems^{[nedostupný zdroj]} lecture notes by Gerald Teschl.
Notes on Diffy Qs: Differential Equations for Engineers Učebnice úvod do diferenciálních rovnic, autor Jiri Lebl UIUC.
Modeling with ODEs using Scilab A tutorial on how to model a physical system described by ODE using Scilab standard programming language by Openeering team.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[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

[#1-7] Harper 1976, strana 127

[8] Kreyszig 1972, strana 2

[#2-9] 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

[EDEBVP-15] 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

[MHFT-17] 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

[MMPE-19] 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