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. +more 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í

Newtonova zákona. +more.

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ů. +more 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í), 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

: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)).

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. +more 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:

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

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í. +more * 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):.

:y^{(n)} = \sum_{i=0}^{n-1} a_i(x) y^{(i)} + r(x)

:kde ai(x) a r(x) jsou spojité funkce v x. Funkce r(x) se nazývá pravá strana a z jejího tvaru vychází další klasifikace: :* 'Homogenní: Jestliže r(x) = 0 a následně jedno „automatické řešení“ je triviální řešení, y = 0. +more Ř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),. +more, ym(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í. +more Ř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'.

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

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í. +more 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. +more 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. +more 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. +more 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.

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ě. 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. +more 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:

Pro každou počáteční podmínku (x0, y0) 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 Imax.

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ě (x0, y0).

Příklad

:y' = y^2

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 :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

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 C1, C2,. jsou libovolné konstanty (obecně komplexní). +more 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).

:

oblast platnosti="col" px" | Diferenciální rovnice	oblast platnosti="col" px" | Metoda řešení	oblast platnosti="col" px" | Obecné řešení
Separabilní rovnice	Separabilní rovnice	Separabilní rovnice
Prvního řádu, separabilní v x a y (obecný případ, speciální případy viz níže) 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 \,\. +more	Separace proměnných (vydělíme P2Q1).	\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 \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 \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 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	Obecná rovnice prvního řádu	Obecná rovnice prvního řádu
Prvního řádu, homogenní \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í 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 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{align} 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{align} \,\. 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í 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{align} 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{align} \, \.
Obecná rovnice druhého řádu	Obecná rovnice druhého řádu	Obecná rovnice druhého řádu
Druhého řádu, autonomní \frac{\mathrm{d}^2y}{\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}^2y}{\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)	Lineární rovnice (až po řád n)	Lineární rovnice (až po řád n)
Prvního řádu, lineární, nehomogenní, funkce koeficienty \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 \frac{\mathrm{d}^2y}{\mathrm{d}x^2} + b\frac{\mathrm{d}y}{\mathrm{d}x} + cy = r(x)\,\.	Komplementární funkce yc: předpokládáme yc = eαx, provedeme substituci a řešíme polynom v α, nalezneme lineárně nezávislé funkce e^{\alpha_j x}. Určitý integrál yp: obecně metodou variace konstant; pro velmi jednoduché r(x) můžeme řešení odhadnout.	y=y_c+y_p Jestliže b2 > 4c, pak: y_c=C_1e^{ \left ( -b+\sqrt{b^2 - 4c} \right )\frac{x}{2}} + C_2e^{-\left ( b+\sqrt{b^2 - 4c} \right )\frac{x}{2}}\,\. Jestliže b2 = 4c, pak: y_c = (C_1x + C_2)e^{-bx/2}\,\. Jestliže b2 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.

\sum_{j=0}^n b_j \frac{\mathrm{d}^j y}{\mathrm{d}x^j} = r(x)\,\! || Komplementární funkce yc: předpokládáme yc = eαx, provedeme substituci a řešíme polynom v α, nalezneme lineárně nezávislé funkce e^{\alpha_j x}.

Určitý integrál yp: obecně metodou variace konstant, i když pro velmi jednoduché r(x) můžeme řešení odhadnout. || 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ý kj-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_je^{\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 * Směrové pole

Reference

Literatura

Literatura v češtině:

*

Literatura v angličtině:

* * * W. Johnson, [url=http://www. +morehti. umich. edu/cgi/b/bib/bibperm. q1=abv5010. 0001. 001]A Treatise on Ordinary and Partial Differential Equations[/url], John Wiley a Sons, 1913, v [url=http://hti. umich. edu/u/umhistmath/]University of Michigan Historical Math Collection[/url] * * * .

Externí odkazy

Hazewinkel, Michiel, ed. (2001), "Differential equation, ordinary", Encyclopedia of Mathematics, Springer, 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). +more * [url=http://eqworld. ipmnet. ru/index. htm]EqWorld: The World of Mathematical Equations[/url], obsahující seznam obyčejných diferenciálních rovnic s řešeními. * [url=http://tutorial. math. lamar. edu/třídy/de/de. aspx]Online Notes / Differential Equations[/url] autor Paul Dawkins, Lamar University. * [url=http://www. sosmath. com/diffeq/diffeq. html]Differential Equations[/url], S. O. S. Mathematics. * [url=http://numericalmethods. eng. usf. edu/mws/gen/08ode/mws_gen_ode_bck_primer. pdf]A primer on analytical solution of differential equations[/url] from the Holistic Numerical Methods Institute, University of South Florida. * [url=http://www. mat. univie. ac. v/~gerald/ftp/book-ode/]Ordinary Differential Equations and Dynamical Systems[/url] lecture notes by Gerald Teschl. * [url=http://www. jirka. org/diffyqs/]Notes on Diffy Qs: Differential Equations for Engineers[/url] Učebnice úvod do diferenciálních rovnic, autor Jiri Lebl UIUC. * [url=http://www. openeering. com/sites/default/files/LHY_Scilab_Tutorial_Part1. pdf]Modeling with ODEs using Scilab[/url] A tutorial on how to model a physical system described by ODE using Scilab standard programming language by Openeering team.

Kategorie:Diferenciální počet Kategorie:Diferenciální rovnice Kategorie:Obyčejné diferenciální rovnice