Literatura Modely založené na změnách ODR 1. řádu
Drsná matematika III – 7. týden
Diferenciální rovnice
Jan Slovák
Masarykova univerzita
Fakulta informatiky
31.10. – 4. 11. 2016
Literatura Modely založené na změnách ODR 1. řádu
Obsah přednášky
1 Literatura
2 Modely založené na změnách
Lineární a nelineární modely
Tlumený oscilátor
3 ODR 1. řádu
Rovnice se separovanými proměnnými
Systémy ODR prvního řádu
Obecné věty
Literatura Modely založené na změnách ODR 1. řádu
Kde je dobré číst?
Zuzana Došlá, Roman Plch, Petr Sojka, Diferenciální počet
funkcí více proměnných s programem Maple, MU Brno, 1999,
273 s.
J. Slovák, M. Panák, M. Bulant, Matematika drsně a svižně,
Muni Press, Brno 2013, v+773 s., elektronická edice
www.math.muni.cz/Matematika_drsne_svizne
Riley, K.F., Hobson, M.P., Bence, S.J. Mathematical Methods
for Physics and Engineering, second edition, Cambridge
University Press, Cambridge 2004, ISBN 0 521 89067 5, xxiii
+ 1232 pp.
Literatura Modely založené na změnách ODR 1. řádu
Lineární model prvního řádu
Derivace pracuje s okamžitými změnami studovaných veličin. Ze
stejných důvodů jsme kdysi zaváděli diference a právě vztahy mezi
hodnotami veličin a změnami těch samých nebo jiných veličin vedly
k rovnicím.
Nejjednodušším modelem bylo úročení vkladů nebo půjček a totéž
pro tzv. Malthusiánský model populace. Přírůstek byl úměrný
hodnotě.
V rámci spojitého modelování by stejný požadavek vedl na rovnici
vztahující derivaci funkce y (x) s její hodnotou
y (x) = r · y(x)
s konstantou úměrnosti r. Je snadné uhodnout řešení této rovnosti
y(x) = C erx
s libovolnou konstantou C.
Tuto konstantu určíme jednoznačně volbou tzv. počáteční
hodnoty y0 = y(x0) v nějakém bodě x0.
Literatura Modely založené na změnách ODR 1. řádu
Pokud je část růstu v našem modelu dána konstatním působením
nezávislém na hodnotě y nebo x (např. paušální poplatky za vedení
účtu nebo přirozený úbytek populace třeba v důsledku porážek na
jatkách), přidáme na pravé straně konstantu s:
y (x) = r · y(x) + s.
Zjevně bude řešením této rovnice funkce
y(x) = C erx
−
s
r
.
K tomuto závěru je velice lehké dojít, pokud si uvědomíme, že
množinou všech řešení první (homogenní) rovnice je jednorozměrný
vektorový prostor, zatímco řešení druhé (nehomogenní) rovnice se
obdrží přičtením kteréhokoliv jednoho jejího řešení ke všem řešením
předchozí rovnice. Lze pak snadno najít konstantní řešení y(x) = k
pro k = −s
r .
Literatura Modely založené na změnách ODR 1. řádu
Nelineární model prvního řádu
U diferenčních rovnic jsme diskutovali tzv. logistický model
populačního růstu založený na předpokladu, že poměr změny
velikosti populace p(n + 1) − p(n) a její velikosti p(n) je v afinní
závislosti na samotné velikosti populace.
Nyní můžeme stejně zavést spojitý model pro populaci p(t) závislou
na čase t jako
p (t) = p(t) −
r
K
p(t) + r ,
tj. při hodnotě p(t) = K pro velkou konstantu K je přírůstek
nulový, zatímco pro p(t) blízké nule je poměr rychlosti růstu
populace k její velikosti blízký r (což je malé číslo v řádu setin
vyjadřující rychlost růstu populace za dobrých podmínek.)
Derivováním ověříme, že následující funkce je řešením pro každou
konstantu C ∈ R
p(t) =
K
1 + CK e−rt
.
Literatura Modely založené na změnách ODR 1. řádu
Srovnání diskrétního a spojitého modelu
y
100
80
60
40
20
t
0
200150100500
100
80
60
40
20
x
200150100500
Srovnáním červeného grafu řešení s K = 100, r = 0, 05 a C = 1
(volba C odpovídá p(0) = 1) s řešením diferenční rovnice (napravo)
vidíme, že skutečně oba přístupy k modelování populací dávají
docela podobné výsledky. Pro srovnání je do levého obrázku zeleně
vkresleno řešení Malthusiánského modelu s odpovídajícími daty.
Literatura Modely založené na změnách ODR 1. řádu
Zkusme si popsat jednoduchý model pro pohyb nějakého tělesa
upnutého k jednomu bodu silnou pružinou. Je-li y(t) výchylka
našeho tělesa od bodu y0 = y(0) = 0, pak lze uvažovat, že
zrychlení y (t) v čase t bude úměrné velikosti výchylky, avšak s
opačným znaménkem. Dostáváme tedy tzv. rovnici oscilátoru
y (t) = −y(t).
Tato rovnice odpovídá systému rovnic
x (t) = −y(t), y (t) = x(t).
Literatura Modely založené na změnách ODR 1. řádu
Řešením takového systému je
x(t) = R cos(t − τ), y(t) = R sin(t − τ)
s libovolnou nezápornou konstantou R, která určuje maximální
amplitudu, a konstantou τ, která určuje fázový posun.
Pro určení jednoznačného řešení potřebujeme proto znát nejen
počáteční polohu y0, nýbrž také rychlost pohybu v tomto okamžiku.
Těmito dvěma údaji bude určena jak amplituda tak fázový posun
jednoznačně.
Předpokládejme, že vlivem vlastností materiálu pružiny bude také
působit síla úměrná okamžité rychlosti pohybu našeho objektu, s
opačným znaménkem než má amplituda.
Literatura Modely založené na změnách ODR 1. řádu
To vyjádříme dodatečným členem s první derivací a naše rovnice je
y (t) = −y(t) − αy (t),
kde α je konstanta, která vyjadřuje velikost tlumení. Na obrázku
jsou fázové diagramy pro řešení s dvěmi různými počátečními
podmínkami. Nalevo je nulové tlumení, napravo je α = 0.3
0
5
-3-3
10-2
-2
t-1
-1
0 15
0
1x(t)
y(t)
1
2
203
2
3
Tlumeneoscilace
0
5
-3-3
10-2
-2
t-1
-1
0 15
0
1x(t)
y(t)
1
2
203
2
3
Tlumeneoscilace
Samotné oscilace jsou vyjádřeny hodnotami na ose y, hodnoty x
zobrazují rychlost pohybu.
Literatura Modely založené na změnách ODR 1. řádu
Definition
Diferenciální rovnice prvního řádu je
F(y , y, x) = 0
s nějakou pevnou funkcí F, která každé trojici reálných čísel přiřadí
jedno reálné číslo.
Zápis připomíná implicitně zadané funkce y(x), nicméně navíc je tu
závislost na derivaci hledané funkce y.
Obvykle pracujeme s rovnicemi, které jsou vyřešeny vzhledem k
derivaci, tj.
y = f (x, y),
dále se omezíme jen na tento případ.
Všimněme si také konvence, že z kontextu víme, která je nezávislá
proměnná pro funkci y.
Literatura Modely založené na změnách ODR 1. řádu
Taková rovnice zadává pro každou hodnotu (x, y) v rovině vektor
(1, f (x, y)), tj. rychlost se kterou nám rovnice grafu řešení přikazuje
pohybovat se rovinou. Např. pro logistický model
p = p −
r
K
p + r
0
x
200150100500
y(x)
100
80
60
40
20
Intuitivně lze na základě takových obrázků očekávat, že pro každou
počáteční podmínku bude existovat právě jedno řešení naší rovnice.
Literatura Modely založené na změnách ODR 1. řádu
Užitečným typem rovnic, pro který známe elementární postup k
řešení jsou tzv. rovnice se separovanými proměnnými:
y (x) = f (x) · g(y(x))
pro dvě dostatečně hladké funkce jedné reálné proměnné f a g.
Obecné řešení tu lze získat integrací, tj. nalezením primitivních
funkcí
G(y) =
dy
g(y)
, F(x) = f (x)dx.
Pak totiž spočtením funkce y(x) z implicitně zadaného vztahu
F(x) + C = G(y) s libovolnou konstantou C vede k řešení.
Skutečně, derivováním této rovnosti (s použitím pravidla pro
derivování složené funkce G(y(x)) dostaneme skutečně
1
g(y) · y (x) = f (x).
Literatura Modely založené na změnách ODR 1. řádu
Příklad rovnice
Najdeme řešení rovnice
y (x) = x · y(x).
Přímým výpočtem dostaneme ln |y(x)| = 1
2 x2 + C. Odtud to
vypadá (alespoň pro kladná y) na
y(x) = e
1
2
x2+C
= D · e
1
2
x2
,
kde D je nyní libovolná kladná konstanta.
Ve skutečnosti dostaneme všechna řešení, když uvážíme všechna
reálná D
Literatura Modely založené na změnách ODR 1. řádu
y(x)
x
3
3
2
1
2
0
-1
1
-2
-3
0-1-2-3
y(x)
x
3
3
2
1
2
0
-1
1
-2
-3
0-1-2-3
Na obrázku jsou vynesena dvě řešení, která ukazují na nestabilitu
rovnice vůči počátečním podmínkám: Jestliže pro libovolné x0
volíme y0 blízké nule, pak se nám dramaticky mění chování
výsledného řešení. Navíc si povšimněme konstatního řešení
y(x) = 0, které odpovídá počáteční podmínce y(x0) = 0.
Literatura Modely založené na změnách ODR 1. řádu
Jestliže lehce pozměníme rovnici na
y (x) = 1 − x · y(x),
narazíme naopak na stabilní chování viditelné na následujícím
obrázku. Tuto rovnici už ale neumíme řešit pomocí separace
proměnných.
y(x)
x
3
6
2,5
2
5
1,5
1
4
0,5
0
3210
y(x)
x
3
6
2,5
2
5
1,5
1
4
0,5
0
3210
Zato umíme stejným postupem řešit logistický model populace.
Literatura Modely založené na změnách ODR 1. řádu
Obecná lineární rovnice první řádu
Jde o rovnici tvaru y = a(t)y + b(t) se spojitými koeficienty a(t) a
b(t).
V případě b(t) = 0 snadno najdeme řešení y(t) = y0F(t, t0), kde
F(t, s) = e
t
s a(x)dx
Obecné řešení s počáteční podmínkou y(t0) = y0 je pak
y(t) = y0F(t, t0) +
t
t0
F(t, s)b(s)ds.
Vzpomeňte podobnost s obecnou lineární diferenční rovnicí prvního
řádu.
Literatura Modely založené na změnách ODR 1. řádu
Tranformace souřadnic
homogenní rovnice y = f (y
t ) převedeme transformací z = y
t pro
t = 0 na
z =
1
t
(f (z) − z).
rovnice Bernoulliho typu y = f (t)y + g(t)yr , kde r = 0, 1, jsou
transformací z = y1−r převedeny na obecnou lineární rovnici.
Literatura Modely založené na změnách ODR 1. řádu
Existence i jednoznačnost řešení skutečně platí pro všechny
rozumné funkce f , my si výsledek sformulujeme pro dosti velkou
třídu rovnic takto:
Theorem (O existenci a jednoznačnosti řešení ODR)
Nechť funkce f (x, y) : R2 → R má spojité parciální derivace. Pak
pro každý bod (x0, y0) ∈ R2 existuje maximální interval
[x0 − a, x0 + b], kde a, b ∈ R jsou kladná čísla, a právě jedna funkce
y(x) : R → R, která je řešením rovnice
y = f (x, y)
a splňuje y(x0) = y0.
Literatura Modely založené na změnách ODR 1. řádu
Náznak důkazu
Funkce y(x) je řešením naší rovnice tehdy a jen tehdy, když
y(x) = y0 +
x
x0
y (x) dx = y0 +
x
x0
f (x, y(x)) dx.
Pravá strana tohoto výrazu je ovšem, až na konstantu, integrální
operátor
L(y)(x) = y0 +
x
x0
f (x, y(x)) dx
a při řešení diferenciální rovnice hledáme pevný bod pro tento
operátor L, tj. chceme najít funkci y = y(x) s L(y) = y.
Důkaz spočívá v odhadu, že pro dostatečně malý interval kolem x0
bude takový operátor zmenšovat vzdálenost funkcí. Z obecné věty o
kontrakci pak vyplývá hledané tvrzení.
Literatura Modely založené na změnách ODR 1. řádu
Na řešení rovnice y (x) = f (x, y) lze také pohlížet jako na hledání
(parametrizované) křivky (x(t), y(t)) v rovině, kde jsme již předem
pevně zvolili parametrizaci proměnné x(t) = t. Takto můžeme
nejen zapomenout na tuto pevnou volbu pro jednu proměnnou x,
nýbrž hlavně přibrat libovolný počet proměnných.
Například v rovině můžeme psát takový systém ve tvaru
x (t) = f (t, x(t), y(t)), y (t) = g(t, x(t), y(t))
se dvěmi funkcemi f , g : R3 → R se spojitými derivacemi. Obdobně
pro více proměnných.
Literatura Modely založené na změnách ODR 1. řádu
Jednoduchým příkladem v rovině může sloužit systém rovnic
x (t) = −y(t), y (t) = x(t).
Snadno lze uhádnout (nebo aspoň ověřit), že řešením takového
systému je
x(t) = R cos t, y(t) = R sin t
s libovolnou nezápornou konstantou R a křivky řešení budou právě
parametrizované kružnice o poloměru R.
Na takové systémy umíme přímo rozšířit platnost věty o
jednoznačnosti a řešení:
Literatura Modely založené na změnách ODR 1. řádu
Theorem (O existenci a jednoznačnosti řešení systémů ODR)
Nechť funkce fi (t, x1, . . . , xn) : Rn+1 → R, i = 1, . . . , n všechny
mají spojité parciální derivace. Pak pro každý bod
(t0, z1, . . . , zn) ∈ R2 existuje interval [t0 − a, t0 + a], s a ∈ R
kladným, a právě jedna funkce y(t) : R → Rn, která je řešením
systému rovnic
x1(x) = f1(t, x1(t), . . . , xn(t)), . . . , xn(x) = fn(t, x1(t), . . . , xn(t))
s počáteční podmínkou
x1(t0) = z1, . . . , xn(t0) = zn.
Navíc je i závislost řešení na počátečních podmínkách a případných
dalších parametrech diferencovatelná.
Literatura Modely založené na změnách ODR 1. řádu
Ve skutečnosti je možné se omezit na tzv. autonomní systémy
rovnic, tj. ty kde pravá strana není explicitně závislá na čase t.
Obecný případ ve větě se z nich dostane stejným způsobem, jako
jsme pracovali v rovině s jednou rovnicí y (x) = f (x, y).
Stejně tak je možné všechny dodatečné parametry λ zahrnout mezi
proměnné, jestliže požadujeme (vektorovou) rovnost λ = 0.
Dříve diskutovaný systém dvou autnomních rovnic, jehož řešení
jsou parametrizované kružnice je dobrým příkladem.
Samotný systém rovnic (a každý obecný také) si pak můžeme
představit jako „pole vektorů“ v rovině, prostoru atd., které udává
„tok prostředí“. Řešením je pak křivka, která odpovídá skutečnému
toku jedné „částečky“ tohoto prostředí.
Literatura Modely založené na změnách ODR 1. řádu
Model „Lotka – Voltera“ pro dravce a kořist
Jako o něco složitější příklad si uveďme klasický populační model
„dravec – kořist“, který zavedli ve dvacátých létech minulého století
pánové Lotka a Volterra.
Označíme x(t) vývoj počtu jedinců v populaci kořisti a y(t) totéž
pro dravce. Přepokládáme, že přírůstek kořisti by se řídil
Malthusiánským modelem (tj. exponenciální růst), kdyby nebyli
loveni. U dravce naopak očekáváme, že by bez kořisti pouze
přirozeně vymíral (tj. exponenciální pokles stavů). Přitom ale ještě
musíme uvážit interakci dravce s kořistí, kterou očekáváme přímo
úměrnou počtu obou. Dostáváme tak tzv. Lotka–Volterra model
x (t) = αx(t) − βy(t)x(t)
y (t) = −γy(t) + δx(t)y(t)
Literatura Modely založené na změnách ODR 1. řádu
Stabilita řešení
Theorem
Nechť x(t) a y(t) jsou dvě pevně zvolená řešení
x (t) = f (t, x(t)), y (t) = g(t, y(t))
systémů rovnic, zadaná počátečními podmínkami x(t0) = x0 a
y(t0) = y0. Potom
|x(t) − y(t)| ≤ |x0 − y0| eC|t−t0|
+B
C eC|t−t0|
−1 ,
kde
C = sup
x=y; (t,x), (t,y)∈U
|f (t, x) − f (t, y)|
|x − y|
B = sup
(t,x)∈U
|f (t, x) − g(t, x)|