Metody řešení diferenčních rovnic
Obsah
1 Rovnice řešitelné sumací 1
1.1 Posun, diference, antidiference, suma . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Nejjednodušší diferenční rovnice . . . . . . . . . . . . . . . . . . . . . . . . . 4
2 Lineární rovnice 7
2.1 Rovnice prvního řádu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2 Rovnice druhého řádu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.2.1 Struktura řešení homogenní rovnice . . . . . . . . . . . . . . . . . . . 9
2.2.2 Řešení nehomogenní rovnice . . . . . . . . . . . . . . . . . . . . . . . . 12
2.2.3 Homogenní rovnice s konstantními koeficienty . . . . . . . . . . . . . . 15
2.2.4 Nehomogenní rovnice, jejíž pravá strana má konstantní koeficienty . . 18
2.2.5 Cauchyho-Eulerova rovnice . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3 Lineární rovnice k-tého řádu . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.3.1 Fundamentální systém řešení homogenní rovnice . . . . . . . . . . . . 23
2.3.2 Nehomogenní rovnice a metoda variace konstant . . . . . . . . . . . . 24
2.3.3 Homogenní rovnice s konstantními koeficienty . . . . . . . . . . . . . . 27
2.3.4 Rovnice s konstantními koeficienty a speciální pravou stranou . . . . . 32
2.4 Cvičení . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3 Další explicitně řešitelné rovnice 39
3.1 Riccatiho a Bernoulliova rovnice . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.2 Homogenní rovnice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.2.1 Implicitní rovnice x(t + 1)2 + a(t)x(t + 1)x(t) + b(t)x(t)2 = 0 . . . . . 46
3.3 Logaritmicky lineární rovnice . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.4 Rovnice řešitelné speciálními substitucemi . . . . . . . . . . . . . . . . . . . . 48
3.4.1 Goniometrické a hyperbolické substituce . . . . . . . . . . . . . . . . . 48
3.4.2 Logistická rovnice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.5 Cvičení . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
ii OBSAH
Kapitola 1
Rovnice řešitelné sumací
1.1 Posun, diference, antidiference, suma
Intervalem celých čísel rozumíme libovolnou z množin
{p, p + 1, p + 2, . . . , q − 1, q} ,
{p, p + 1, p + 2, . . . } , {. . . , q − 2, q − 1, q} , {. . . , −2, −1, 0, 1, 2, . . . } ,
kde p, q ∈ Z, p < q. Pro interval I celých čísel klademe Iκ = I \ max I. To znamená, že
{p, p + 1, p + 2, . . . , q − 1, q}κ
= {p, p + 1, p + 2, . . . , q − 2, q − 1} ,
{. . . , q − 2, q − 1, q}κ
= {. . . , q − 3, q − 2, q − 1}
a v ostatních případech Iκ = I.
Uvažujme posloupnost a, tj. zobrazení intervalu celých čísel do čísel reálných. Členy této
posloupnosti budeme zapisovat standardním symbolem a(t) jako hodnotu posloupnosti a v indexu
(nezávisle proměnné) t.
Pro každou posloupnost a a index t ∈ (Dom a)κ klademe
aσ
(t) = a(t + 1).
Posloupnost aσ se nazývá posun (podrobněji posun vpřed, anglicky shift) posloupnosti a.
Diference posloupnosti a je posloupnost, označená ∆a, definovaná vztahem
∆a(t) = a(t + 1) − a(t).
Diference posloupnosti a je definována na množině (Dom a)κ. Vztah posunu a diference posloupnosti
je zřejmý,
∆a = aσ
− a.
Pro libovolné posloupnosti a, b se stejným definičním oborem, pro libovolné reálné konstanty
α, β a každý index t ∈ (Dom a)κ platí
∆ αa + βb (t) = α∆a(t) + β∆b(t),
1
2 KAPITOLA 1. ROVNICE ŘEŠITELNÉ SUMACÍ
tj. diference je lineární operátor na prostoru posloupností. Dále
∆ ab (t) = ∆a(t) b(t) + a(t + 1) ∆b(t) = ∆a(t) b(t + 1) + a(t) ∆b(t) .
Pokud navíc b(t) = 0 = b(t + 1), pak
∆
a
b
(t) =
∆a(t) b(t) − a(t) ∆b(t)
b(t)b(t + 1)
.
Stručně odvozené rovnice zapíšeme jako
∆(αa + βb) = α∆a + β∆b, ∆ab = b∆a + aσ
∆b = bσ
∆a + a∆b, ∆
a
b
=
b∆a − a∆b
bbσ
.
Antidiference posloupnosti a je posloupnost A se stejným definičním oborem Dom a, pro
kterou platí ∆A = a, podrobněji
∆A(t) = a(t)
pro všechna t ∈ (Dom a)κ. Antidiference není určena jednoznačně: Je-li γ libovolná konstanta
a A antidiference posloupnosti a, pak posloupnost ˜A = A+γ je také antidiferencí posloupnosti
a, neboť
∆ ˜A(t) = ∆ A + γ (t) = A(t + 1) + γ − A(t) + γ = A(t + 1) − A(t) = ∆A(t) = a(t).
Naopak, pokud posloupnosti A a ˜A jsou antidiference posloupnosti a, pak
a(t) = ∆A(t) = ∆ ˜A(t),
neboli
A(t + 1) − A(t) = ˜A(t + 1) − ˜A(t)
pro všechny indexy t ∈ (Dom a)κ. Odtud dostaneme rovnost
(A − ˜A)(t + 1) = A(t + 1) − ˜A(t + 1) = A(t) − ˜A(t) = (A − ˜A)(t),
která zase platí pro všechny indexy t. To znamená, že posloupnost A − ˜A je konstantní.
Dostáváme tak závěr, že dvě posloupnosti A a ˜A jsou antidiferencí téže posloupnosti a právě
tehdy, když se liší o aditivní konstantu.
Nechť posloupnost A, resp. B, je antidiference posloupnosti a, resp. b, a α, β jsou libovolné
konstanty. Pak lineární kombinace αA + βB je antidiferencí posloupnosti αa + βb.
Nechť a je posloupnost. Její (nějakou) antidiferenci označíme Σa (význam této symboliky
se objasní později). Předchozí výsledek lze při tomto označení zapsat ve tvaru
Σ(αa + βb) = αΣa + βΣb. (1.1)
Z druhé formule pro diferenci součinu posloupností plyne další rovnost
Σ(a∆b) = ab − Σ(bσ
∆a). (1.2)
Rovnosti (1.1) a (1.2) chápeme ve smyslu „až na aditivní konstantu , tj. rozdíl posloupností
na levé a pravé straně těchto rovností je konstantní posloupnost.
Diference a antidiference některých posloupností („elementárních posloupností ) jsou shrnuty
v Tabulce 1.1
1.1. POSUN, DIFERENCE, ANTIDIFERENCE, SUMA 3
posloupnost diference antidiference
1. a(t) = 1 0 t
2. a(t) = t 1 1
2 t(t − 1)
3. a(t) = t(r) =
t
j=t−r+1
j, t ≥ 0 rt(r−1) t(r+1)
r + 1
4. a(t) = αt, α = 1 (α − 1)αt αt
α − 1
5. a(t) = cos(αt + β), α ≡ 0, mod 2π −2 sin 1
2α sin(αt + β + 1
2 α)
sin(αt + β − 1
2α)
2 sin 1
2α
6. a(t) = sin(αt + β), α ≡ 0, mod 2π 2 sin 1
2 α cos(αt + β + 1
2α) −
cos(αt + β − 1
2α)
2 sin 1
2α
Tabulka 1.1: Diference a antidiference některých posloupností. Posloupnost na řádku 1. je
konstantní, na ř. 2. aritmetická s diferencí 1, posloupnost na řádku 3. se nazývá faktoriálová,
na řádku 4. je geometrická posloupnost s kvocientem α, na řádcích 5. a 6. jsou posloupnosti
goniometrické.
Sumace (a součiny).
Nechť a je posloupnost, p, q ∈ Dom(a) takové, že p ≤ q. Sumu členů posloupnosti a od p do
q definujeme obvyklým způsobem jako
q
j=p
a(j) = a(q) + a(q − 1) + · · · + a(p + 1) + a(p)
a podobně součin členů posloupnosti a od p do q jako
q
j=p
a(j) = a(q)a(q − 1) · · · a(p + 1)a(p).
Dále klademe
p−1
j=p
a(j) = 0,
p−1
j=p
a(j) = 1.
Pro p, q ∈ Dom a takové, že p > q klademe
q
j=p
a(j) = −
p−1
j=q+1
a(j)
a pokud jsou všechny členy a(q + 1), a(q + 2), . . . , a(p − 2), a(p − 1) nenulové, klademe také
q
j=p
a(j) =


p−1
j=q+1
a(j)


−1
.
4 KAPITOLA 1. ROVNICE ŘEŠITELNÉ SUMACÍ
Při této konvenci pro libovolné hodnoty p, q, r ∈ Dom a platí
q
j=p
a(j) +
r
j=q+1
a(j) =
r
j=p
a(j),
q
j=p
a(j)
r
j=q+1
a(j) =
r
j=p
a(j).
Suma diferencí posloupnosti a splňuje rovnosti
t−1
j=t0
∆a(j) =
t−1
j=t0
a(j + 1) − a(j) =
t
j=t0+1
a(j) −
t−1
j=t0
a(j) = a(t) − a(t0). (1.3)
Odtud vidíme, že pro posloupnost a a její antidiferenci A platí
t−1
j=t0
a(j) = A(t) − A(t0);
zavedeme-li pro libovolnou posloupnost x označení
x(j)
t
j=t0
= x(t) − x(t0),
můžeme předchozí vásledek přepsat ve tvaru
t−1
j=t0
a(j) = A(j)
t
j=t0
.
Při znalosti antidiference posloupnosti tedy již snadno spočítáme konečný součet členů po-
sloupnosti.
Z rovnosti (1.2) získáme užitečnou formuli
t−1
j=t0
a(j)∆b(j) = a(j)b(j)
t
j=t0
−
t−1
j=t0
b(j + 1)∆a(j); (1.4)
nazýváme ji sumace „per partes .
1.2 Nejjednodušší diferenční rovnice
Počáteční úlohu pro diferenční rovnici prvního druhu ve tvaru
∆x = b(t), x(t0) = x0, (1.5)
kde x je hledaná posloupnost a b je daná posloupnost, můžeme bezprostředně vyřešit sumací
obou jejích stran. Podle (1.3) je totiž
x(t) − x(t0) =
t−1
j=t0
b(j),
takže
x(t) = x0 +
t−1
j=t0
b(j).
Pokud známe antidiferenci B posloupnosti b, řešení úlohy (1.5) zapíšeme ve tvaru
x(t) = x0 + B(j)
t
j=t0
.
1.2. NEJJEDNODUŠŠÍ DIFERENČNÍ ROVNICE 5
Příklady: Najdeme řešení počáteční úlohy:
1. ∆x = t3, x(1) = 1.
Posloupnost na pravé straně rozepíšeme pomocí faktoriálových posloupností,
t3
= t(t − 1)(t − 2) + 3t2
− 2t = t(t − 1)(t − 2) + 3t(t − 1) + t = t(3)
+ 3t(2)
+ t(1)
.
S využitím třetího a druhého řádku tabulky 1.1 dostaneme
t−1
j=1
j3
=
t−1
j=1
j(3)
+ 3
t−1
j=1
j(2)
+
t−1
j=1
j =
j(4)
4
t
j=1
+ 3
j(3)
3
t
j=1
+
j(j − 1)
2
t
j=1
=
=
1
4
t(t − 1)(t − 2) − 0 + 3
1
3
t(t − 1)(t − 2) − 0 +
1
2
t(t − 1) − 0 =
1
4
t2
(t − 1)2
.
Řešení dané úlohy tedy je posloupnost daná předpisem
x(t) = 1 +
t2(t − 1)2
4
.
2. ∆x = 2tt, x(0) = 0.
Pro nalezení sumy posloupnosti na pravé straně rovnice využijeme sumaci „per partes
(1.4) (kde a(t) = t, ∆b(t) = 2t, tj. b(t) = 2t) a čtvrtý řádek tabulky 1.1:
t−1
j=0
2j
j = 2j
j
t
j=0
−
t−1
j=0
2j+1
= 2t
t − 0 −
2j+1
1
t
t=0
= 2t
t − 2t+1
+ 2.
Řešení dané rovnice tedy je x(t) = 2t(t − 2) + 2.
3. x(t + 1) = 1 −
1
t + 1
x(t) +
(−1)t
t + 1
, x(1) = 1
2 .
Danou rovnici upravíme na tvar
(1 + t)x(t + 1) − tx(t) = (−1)t
, tj. ∆ tx(t) = cos tπ.
Nyní zavedeme novou neznámou posloupnost y substitucí y(t) = tx(t). Posloupnost y
vyhovuje diferenční rovnici
∆y(t) = cos tπ
s počáteční podmínkou y(1) = 1 · x(1) = 1
2. Pro řešení této úlohy využijeme pátý řádek
tabulky 1.1:
y(t) = 1
2 +
t−1
j=1
cos tπ = 1
2 +
sin(jπ− π
2 )
2 sin π
2
t
j=1
= 1
2 + 1
2 sin tπ − π
2 − 1 = 1
2 sin(2t − 1)π
2 .
Zpětnou substitucí dostaneme řešení dané úlohy ve tvaru
x(t) =
sin(2t − 1)π
2
2t
=
(−1)t−1
2t
.
6 KAPITOLA 1. ROVNICE ŘEŠITELNÉ SUMACÍ
Kapitola 2
Lineární rovnice
2.1 Rovnice prvního řádu
Lineární diferenční rovnice prvního řádu je tvaru
∆x = a(t)x + b(t), (2.1)
kde a, b jsou nějaké posloupnosti. Pokud je posloupnost b identicky nulová, b ≡ 0, nazýváme
rovnici (2.1) homogenní.
Diferenční rovnici (2.1) můžeme přepsat do tvaru rekurentní formule (diferenční rovnice
druhého typu)
x(t + 1) = a(t) + 1 x(t) + b(t),
nebo
x(t + 1) − q(t)x(t) = b(t), (2.2)
při označení q(t) = a(t) + 1. K rovnici přidáme počáteční podmínku
x(t0) = x0. (2.3)
Z rovností (2.2), (2.3) postupně počítáme:
x(t0 + 1) = q(t0)x0 + b(t0) = x0q(t0) + b(t0),
x(t0 + 2) = q(t0 + 1)x(t0) + b(t0 + 1) = q(t0 + 1) x0q(t0) + b(t0) + b(t0 + 1) =
= x0q(t0)q(t0 + 1) + b(t0)q(t0 + 1) + b(t0 + 1),
x(t0 + 3) = q(t0 + 2) x0q(t0)q(t0 + 1) + b(t0)q(t0 + 1) + b(t0 + 1) + b(t0 + 2) =
= x0q(t0)q(t0 + 1)q(t0 + 2) +
+ b(t0)q(t0 + 1)q(t0 + 2) + b(t0 + 1)q(t0 + 2) + b(t0 + 2) =
= x0
2
i=0
q(t0 + i) +
2
i=0
b(t0 + i)
2
j=i+1
q(t0 + j),
atd. Obecně
x(t0 + l) = x0
l−1
i=0
q(t0 + i) +
l−1
i=0
b(t0 + i)
l−1
j=i+1
q(t0 + j) = x0
t0+l−1
i=t0
q(i) +
t0+l−1
i=t0
b(i)
t0+l−1
j=i+1
q(j).
Tento výsledek lze snadno ověřit úplnou indukcí. Dostáváme tak
7
8 KAPITOLA 2. LINEÁRNÍ ROVNICE
Tvrzení 1. Řešení lineární diferenční rovnice (2.1) s počáteční podmínkou (2.3) je dáno
formulí
x(t) = x0
t−1
i=t0
1 + a(i) +
t−1
i=t0
b(i)
t−1
j=i+1
1 + a(j) ;
řešení lineární rekurentní formule prvního řádu (2.2) s počáteční podmínkou (2.3) je dáno
formulí
x(t) = x0
t−1
i=t0
q(i) +
t−1
i=t0
b(i)
t−1
j=i+1
q(j).
Významné speciální případy lineární rovnice prvního řádu jsou ty, ve kterých je některá
z posloupností a, b konstantní, a ≡ α = −1 nebo b ≡ β.
Důsledky.
1. Řešení lineární diferenční rovnice, resp. rekurentní formule, prvního řádu
∆x = αx + b(t), resp. x(t + 1) = κx(t) + b(t),
s počáteční podmínkou (2.3) je dáno formulí
x(t) = x0(1 + α)t−t0
+
t−1
i=t0
b(i)(1 + α)t−1−i
,
resp.
x(t) = x0κt−t0
+
t−1
i=t0
b(i)κt−1−i
.
2. Řešení lineární diferenční rovnice, resp. rekurentní formule, prvního řádu
∆x = αx + β, resp. x(t + 1) = κx(t) + β,
s počáteční podmínkou (2.3) je dáno formulí
x(t) = x0(1 + α)t−t0
−
β
α
1 − (1 + α)t−t0
= x0 +
β
α
(1 + α)t−t0
−
β
α
,
resp.
x(t) = x0 −
β
1 − κ
κt−t0
+
β
1 − κ
.
(Při odvození druhého důsledku byl využit vzorec pro součet konečné geometrické posloup-
nosti.)
2.2. ROVNICE DRUHÉHO ŘÁDU 9
2.2 Rovnice druhého řádu
Lineární diferenční rovnice druhého typu (rekurentní formule) druhého řádu je tvaru
x(t + 2) + a1(t)x(t + 1) + a0(t)x(t) = b(t), (2.4)
kde a0 je posloupnost taková, že a(t) = 0. Pokud je posloupnost b identicky nulová, b ≡ 0,
nazýváme rovnici (2.4) homogenní.
K rovnice (2.4) přísluší počáteční podmínka
x(t0) = ξ0, x(t0 + 1) = ξ1. (2.5)
Počáteční úloha pro rovnici (2.4) má zřejmě jediné řešení definované na průniku definičních
oborů posloupností a0, a1, b. Na této množině totiž můžeme spočítat člen hledané posloupnosti
x ze dvou „předchozích ,
x(t + 2) = b(t) − a0(t)x(t) − a1(t)x(t + 1),
a díky předpokladu a0(t) můžeme také spočítat člen posloupnosti ze dvou „následujících ,
x(t) =
1
a0(t)
b(t) − a1(t)x(t + 1) − x(t + 2) .
Z počáteční podmínky tedy určíme členy hledané posloupnosti pro všechny přípustné hodnoty
indexu t.
2.2.1 Struktura řešení homogenní rovnice
Pro lineární homogenní rovnici druhého řádu
x(t + 2) + a1(t)x(t + 1) + a0(t)x(t) = 0, (2.6)
platí princip superpozice:
Tvrzení 2. Jsou-li x1 a x2 řešení rovnice (2.6) a c1, c2 konstanty, pak také lineární kombinace
y = c1x1 + c2x2 je řešením této rovnice.
Důkaz:
y(t + 2) + a1(t)y(t + 1) + a0(t)y(t) =
= c1x1(t + 2) + c2x2(t + 2) + a1(t) c1x1(t + 1) + c2x2(t + 1) + a0 c1x1(t) + c2x2(t) =
= c1 x1(t + 2) + a1(t)x1(t + 1) + a0(t)x1(t) +
+ c2 x2(t + 2) + a1(t)x2(t + 1) + a0(t)x2(t) = 0.
Bezprostředním důsledkem tohoto tvrzení je evidentní skutečnost, že nulová posloupnost
je také řešením rovnice. Celkem tak dostáváme, že množina řešení lineární homogenní rovnice
(2.6) je lineární (vektorový) prostor.
Řešení y1 rovnice (2.6) s počáteční podmínkou
x(t0) = 1, x(t0 + 1) = 0
10 KAPITOLA 2. LINEÁRNÍ ROVNICE
a řešení y2 této rovnice s počáteční podmínkou
x(t0) = 0, x(t0 + 1) = 1
jsou evidentně lineárně nezávislá. To znamená, že prostor řešení lineární homogenní rovnice
(2.6) má dimenzi alespoň 2. Naopak, pro jednoznačně určené řešení x rovnice (2.6) s obecnou
počáteční podmínkou (2.5) platí
x = ξ1y1 + ξ2y2;
přitom y1, y2 splnují uvedené počáteční podmínky. To znamená, že libovolné řešení homogenní
rovnice (2.6) lze vyjádřit jako lineární kombinaci dvou posloupností y1 a y2. Prostor řešení
má tedy dimenzi nejvýše 2.
Provedené úvahy shrneme: Množina všech řešení lineární homogenní rovnice druhého řádu
(2.6) tvoří dvourozměrný lineární prostor. Bázi tohoto prostoru nazveme fundamentální systém
řešení homogenní rovnice (2.6).
Tvoří-li posloupnosti y1 a y2 fundamentální systém řešení homogenní rovnice (2.6), pak
posloupnost c1y1 + c2y2 je obecným řešením lineární homogenní rovnice (2.6). Rešení počáteční
úlohy (2.6), (2.5) je totiž tohoto tvaru, přičemž konstanty c1, c2 jsou řešením soustavy
lineárních (algebraických) rovnic
c1y1(t0) + c2y2(t0) = ξ0
c1y1(t0 + 1) + c2y2(t0 + 2) = ξ1.
Známe-li tedy fundamentální systém y1, y2 řešení lineární homogenní rovnice (2.6), můžeme
psát řešení počáteční úlohy (2.6), (2.5) ve tvaru
x(t) =
ξ0y2(t0 + 1) − ξ1y2(t0) y1(t) − ξ0y1(t0 + 1) − ξ1y1(t0) y2(t)
y1(t0)y2(t0 + 1) − y1(t0 + 1)y2(t0)
.
Výpočet druhé složky fundamentálního systému
Předpokládejme, že známe řešení y1 homogenní rovnice (2.6) takové, že y1(t) = 0.
Rovnice (2.6) je diferenční, proto by se ve vyjádření jejího řešení mohla objevovat sumace.
Z tohoto důvodu budeme hledat řešení této rovnice ve tvaru
x(t) = y1(t)
t−1
i=t0
u(i), (2.7)
kde u je zatím neznámá posloupnost. Pak
x(t + 1) = y1(t + 1)
t
i=t0
u(i) = y1(t + 1)
t−1
i=t0
u(i) + u(t) ,
x(t + 2) = y1(t + 2)
t−1
i=t0
u(i) + u(t) + u(t + 1) .
Toto vyjádření dosadíme do rovnice (2.6),
u1(t + 2) + a1(t)y1(t + 1) + a0(t)y1(t)
t−1
i=t0
u(i) + y1(t + 2) + a1y1(t + 1) u(t)+
+ y1(t + 2)u(t + 1) = 0.
2.2. ROVNICE DRUHÉHO ŘÁDU 11
Poněvadž posloupnost y1 je řešením rovnice (2.6), je výraz v první závorce roven nule a výraz
ve druhé závorce je roven −a0(t)y1(t). Podle předpokladu je y1(t + 2) = 0, proto můžeme
tímto členem předchozí rovnost vydělit. Dostaneme
u(t + 1) = a0(t)
y1(t)
y1(t + 2)
u(t).
To je lineární homogenní rovnice prvního řádu pro neznámou posloupnost u, jejíž řešení je
podle Tvrzeni 1 dáno součinem
u(t) = c
t−1
j=t0
a0(j)
y1(j)
y1(j + 1)
,
kde c = u(t0) je nějaká konstanta. Toto vyjádření dosadíme do rovnosti (2.7) a dostaneme
řešení rovnice (2.6) ve tvaru
x(t) = cy1(t)
t−1
i=t0
i−1
j=t0
a0(j)
y1(j)
y1(j + 2)
. (2.8)
Ověříme, že posloupnost daná výrazem (2.8) je pro c = 0 lineárně nezávislá na posloupnosti
y1. Casoratián posloupností y1 a x je
C(t; y1, x) =
y1(t) cy1(t)
t−1
i=t0
i−1
j=t0
a0(j) y1(j)
y1(j+1)
y1(t + 1) cy1(t + 1)
t−1
i=t0
i−1
j=t0
a0(j) y1(j)
y1(j+1) +
t−1
j=t0
a0(j) y1(j)
y1(j+1)
=
= cy1(t)y1(t + 1)
t−1
j=t0
a0(j)
y1(j)
y1(j + 1)
= 0;
nebot podle předpokladů a0(t) = 0, y1(t) = 0. Dostali jsme tedy výsledek:
Tvrzení 3. Nechť posloupnost y1 je první složkou fundamentálního systému řešení rovnice
(2.6) taková, že y1(t) = 0. Pak druhá složka fundamentálního systému řešení je dána výrazem
y2(t) = y1(t)
t−1
i=t0
i−1
j=t0
a0(j)
y1(j)
y1(j + 2)
.
Příklad:
x(t + 2) − (t + 2)x(t + 1) + (t + 2)x(t) = 0, t > 0
V tomto případě je t0 = 1, a0(t) = t+2 a posloupnost y1 daná předpisem y1(t) = t je řešením
dané rovnice, neboť
t + 2 − (t + 2)(t + 1) + (t + 2)t = (t + 2)(1 − t − 1 + t) = 0.
Druhá složka fundamentálního systému řešení tedy je
y2(t) = t
t−1
i=1
i−1
j=1
(j + 2)
j
j + 2
= t
t−1
i=1
(i − 1)! .
12 KAPITOLA 2. LINEÁRNÍ ROVNICE
2.2.2 Řešení nehomogenní rovnice
Pokud jsou „koeficienty a1, a2 v rovnicích (2.4) a (2.6) stejné, řekneme, že rovnice (2.6) je
přidružená homogenní rovnice k lineární rovnici (2.4).
Předpokládejme, že známe fundamentální systém řešení y1, y2 homogenní rovnice (2.6) a
nějaké řešení xN nehomogenní rovnice. Pak posloupnost x, definovaná jako součet obecného
řešení homogenní rovnice (2.6) a posloupnosti xN , tj. x = c1y1 + c2y2 + xN , kde c1, c2 jsou
nějaké konstanty, je řešením nehomogenní rovnice (2.4). Vskutku
x(t + 2) + a1(t)x(t + 1) + a0(t)x(t) =
= c1y1(t + 2) + c2y2(t + 2) + xN (t + 2) +
+ a1(t) c1y1(t + 1) + c2y2(t + 1) + xN (t + 1) + a0(t) c1y1(t) + c2y2(t) + xN (t) =
= c1 y1(t + 2) + a1(t)y1(t + 1) + a0(t)y1(t) +
+ c2 y2(t + 2) + a1(t)y2(t + 1) + a0(t)y2(t) +
+ xN (t + 2) + a1(t)xN (t + 1) + a0(t)xN (t) = 0 + 0 + b(t) = b(t).
Pokud konstanty c1 a c2 zvolíme tak, aby byla splněna počáteční podmínka (2.5), tj. aby
splňovaly soustavu lineárních algebraických rovnic
c1y1(t0) + c2y2(t0) + xN (t0) = ξ0,
c1y1(t0 + 1) + c2y2(t0 + 1) + xN (t0 + 1) = ξ1,
neboli
c1 =
ξ0 − xN (t0) y2(t0 + 1) − ξ1 − xN (t0 + 1) y2(t0)
y1(t0)y2(t0 + 1) − y2(t0)y1(t0 + 1)
,
c2 =
y1(t0) ξ1 − xN (t0 + 1) − y1(t0 + 1) ξ0 − xN (t0)
y1(t0)y2(t0 + 1) − y2(t0)y1(t0 + 1)
,
pak posloupnost x = c1y1 + c2y2 + xN je řešením počáteční úlohy (2.4), (2.5). Touto úvahou
dostáváme
Tvrzení 4. Obecné řešení lineární rovnice (2.4) je součtem obecného řešení přidružené homogenní
rovnice a nějakého řešení nehomogenní rovnice.
Ještě poznamenejme, že výpočet konstant c1, c2 pro konkrétní počáteční úlohu je nejjednodušší,
pokud partikulární řešení xN nehomogenní rovnice volíme tak, aby xN (t0) = 0 a
xN (t0 + 1) = 0.
Nalezení partikulárního řešení nehomogenní rovnice
Budeme hledat řešení rovnice (2.4) s nulovou počáteční podmínkou
x(t0) = 0, x(t0 + 1) = 0. (2.9)
Budeme přitom předpokládat, že známe fundamentální systém řešení y1, y2 přidružené homogenní
rovnice.
Počáteční úlohu (2.4), (2.9) lze interpretovat: Proces, který je popsán nehomogenní rovnicí
(2.4), začíná s nulovými hodnotami a v průběhu času na něho působí nějaké vlivy. Stav procesu
v čase t tedy můžeme popsat jako součet těchto vlivů od počátečního okamžiku t0 do okamžiku
2.2. ROVNICE DRUHÉHO ŘÁDU 13
t − 1 bezprostředně předcházejícímu času t. Přesněji vyjádřeno, řešení počáteční úlohy (2.4),
(2.9) budeme hledat ve tvaru
x(t) =
t−1
i=t0
w(i, t), (2.10)
kde w je nějaká, zatím neurčená, funkce dvou celočíselných proměnných. Myšlenka hledat
řešení nehomogenní rovnice s nulovými počátečními podmínkami ve tvaru součtu (nebo integrálu)
se nazývá Duhamelův princip.
Posloupnost daná rovností (2.10) splňuje první z počátečních podmínek (2.9). Ta druhá
z nich nabývá tvar
0 = x(t0 + 1) = w(t0, t0 + 1).
Tato podmínka bude splněna zejména tehdy, když po funkci w budeme požadovat, aby měla
vlastnost
w(i, i + 1) = 0 (2.11)
pro všechny přípustné hodnoty i.
Z rovnosti (2.10) dále s využitím podmínky (2.11) dostaneme
x(t + 1) =
t
i=t0
w(i, t + 1) =
t−1
i=t0
w(i, t + 1) + w(t, t + 1) =
t−1
i=t0
w(i, t + 1),
x(t + 2) =
t
i=t0
w(i, t + 2) =
t−1
i=t0
w(i, t + 2) + w(t, t + 2).
Toto vyjádření dosadíme do dané rovnice (2.4). Dostaneme
t−1
i=t0
w(i, t + 2) + a1(t)w(i, t + 1) + a0(t)w(i, t) + w(t, t + 2) = b(t).
Tato rovnost bude splněna zejména tehdy, pokud funkce w bude mít vlastnosti
w(i, i + 2) = b(i), pro všechna i, (2.12)
w(i, t + 2) + a1(t)w(i, t + 1) + a0(t)w(i, t) = 0, pro všechna i, t. (2.13)
Nyní budeme proměnnou i považovat za parametr. Získané rovnosti budeme chápat jako
úlohu pro lineární homogenní rovnici druhého řádu (2.13) s počátečními podmínkami (2.11)
a (2.12); hledaná posloupnost je w(i, ·). Rovnice (2.13) je přidruženou homogenní rovnicí
k dané rovnici (2.4). Podle předpokladu tedy známe její fundamentální systém řešení, její
obecné řešení je lineární kombinací bázových posloupností y1, y2. Příslušné koeficienty ovšem
závisí na parametru i, tedy
w(i, t) = c1(i)y1(t) + c2(i)y2(t). (2.14)
Konkrétní vyjádření koeficientů (posloupností) c1, c2 získáme z počátečních podmínek (2.11)
a (2.12):
c1(i)y1(i + 1) + c2(i)y2(i + 1) = 0,
c1(i)y1(i + 2) + c2(i)y2(i + 2) = b(i),
14 KAPITOLA 2. LINEÁRNÍ ROVNICE
tj.
c1(i) =
−b(i)y2(i + 1)
y1(i + 1)y2(i + 2) − y1(i + 2)y2(i + 1)
,
c2(i) =
b(i)y1(i + 1)
y1(i + 1)y2(i + 2) − y1(i + 2)y2(i + 1)
.
Dosazením těchto výrazů do rovnosti (2.14) a pak do (2.10) dostaneme řešení počáteční úlohy
(2.4), (2.9) ve tvaru
x(t) =
t−1
i=t0
b(i)
y2(i + 1)y1(t) − y1(i + 1)y2(t)
y1(i + 2)y2(i + 1) − y1(i + 1)y2(i + 2)
. (2.15)
Při označení Cj(t) =
t−1
i=t0
cj(i), j = 1, 2 můžeme předchozí formuli pro řešení nehomogenní
lineární rovnice druhého řádu přepsat ve tvaru
x(t) = C1(t)y1(t) + C2(t)y2(t),
tedy jako lineární kombinaci bázových posloupností y1, y2; přitom koeficienty nejsou konstantní
ale variabilní, závisí na čase t. Proto se rovnost (2.15) nazývá formule variace kon-
stant.
Příklad: Najdeme řešení počáteční úlohy
x(t + 2) +
1
t
x(t + 1) − 1 +
1
t
x(t) = 4(t + 1), x(1) = 1, x(2) = −1.
Prvním problémem je najít fundamentální systém řešení přidružené homogenní rovnice
x(t + 2) +
1
t
x(t + 1) − 1 +
1
t
x(t) = 0.
Tato rovnice má evidentně konstantní řešení y1 ≡ 1. Druhou složku fundamentálního systému
řešení vypočítáme podle Tvrzení 3:
y2(t) =
t−1
i=1
i−1
j=1
−1 −
1
j
=
t−1
i=1
(−1)i−1
i−1
j=1
j + 1
j
=
t−1
i=1
(−1)i−1 2
1
·
3
2
· · ·
i − 1
i − 2
·
i
i − 1
=
=
t−1
i=1
(−1)i−1
i = 1 − 2 + 3 − 4 + · · · + (−1)t−3
(t − 2) + (−1)t−2
(t − 1) =
=
−1
2(t − 2) + t − 1, t sudé
−1
2(t − 1), t liché
=
1
2t, t sudé
1
2(1 − t), t liché
= 1
4 1 + (−1)t(2t − 1) .
Dostáváme tak obecné řešení xH přidružené homogenní rovnice ve tvaru
xH(t) = A + 1
4B 1 + (−1)t(2t − 1) .
Počáteční hodnoty jsou
1 = xH(1) = A, −1 = xH(2) = A + B.
2.2. ROVNICE DRUHÉHO ŘÁDU 15
Z těchto rovnic snadno vypočítáme A = 1, B = −2. Řešení přidružené homogenní rovnice
které splňuje počáteční podmínky je tedy dáno rovností
xH(t) = 1 − 1
2 1 + (−1)t(2t − 1) = 1
2 − 1
2 (−1)t(2t − 1).
Rešení nehomogenní rovnice, které splnuje nulové počáteční podmínky xN (1) = xN (2) = 0
dostaneme dosazením do vzorce (2.15).
xN (t) =
t−1
i=1
4(i + 1)
1
4 1 − (−1)i(2i + 1) − 1
4 1 + (−1)t(2t − 1)
1
4 (1 − (−1)i(2i + 1)) − 1
4 (1 + (−1)i(2i + 3))
=
=
t−1
i=1
2i + 1 + (−1)t+i
(2t − 1) =
t−1
i=1
(2i + 1) + (−1)t
(2t − 1)
t−1
i=1
(−1)i
=
= 1
2(t − 1)(2t + 2) − (−1)t(2t − 1)1
2 (−1)t + 1) = t(t − 1) − 1
2 − 1
2 (−1)t(2t − 1).
Rešení dané rovnice je součet posloupností xH a xN , tedy
x(t) = t(t − 1) − (−1)t
(2t − 1).
2.2.3 Homogenní rovnice s konstantními koeficienty
Lineární homogenní rovnice prvního řádu
x(t + 1) = κx(t)
má podle důsledku Tvrzení 1 řešení x(t) = κt, tj. geometrickou posloupnost. Toto pozorování
vede k myšlence hledat řešení lineární homogenní rovnice druhého řádu s konstantními
koeficienty
x(t + 2) + α1x(t) + α0x(t) = 0 (2.16)
také ve tvaru geometrické posloupnosti
x(t) = λt
, (2.17)
kde λ je nějaká nenulová konstanta.
Po dosazení výrazu (2.17) do dané rovnice (2.16) dostaneme
λt+2
+ α1λt+1
+ α0λt
= 0 tj. λt
λ2
+ α1λ + α0 = 0.
Podle předpokladu λ = 0 je tato rovnost ekvivalentní s rovnicí
λ2
+ α1λ + α0 = 0. (2.18)
To je kvadratická rovnice, která má kořeny
λ1,2 =
1
2
−α1 ± α2
1 − 4α0 .
16 KAPITOLA 2. LINEÁRNÍ ROVNICE
Z obecného předpokladu kladeného na lineární rovnice druhého řádu plyne α0 = 0, takže oba
kořeny jsou nenulové. Dostáváme tak dvě, ne nutně různá, řešení dané diferenční rovnice ve
tvaru x1(t) = λt
1, x2(t) = λt
2. Casoratián posloupností x1, x2 je
C(t; x1, x2) =
λt
1 λt
2
λt+1
1 λt+1
2
= (λ1λ2)t
(λ2 − λ1).
Pro λ1 = λ2 je C(t; x1, x2) = 0, takže posloupnosti x1, x2 jsou v takovém případě nezávislé.
(Algebraická) rovnice (2.18) se nazývá charakteristická rovnice, její řešení se nazývají
charakteristické kořeny. Rozlišíme tři standardní případy:
Dva reálné různé charakteristické kořeny
Tento případ nastává, pokud α2
1 > 4α0. Bezprostředně můžeme napsat fundamentální systém
řešení rovnice (2.16),
y1(t) = λt
1 = 1
2 α2
1 − 4α0 − α1
t
, y2(t) = λt
2 = −1
2 α2
1 − 4α0 + α1
t
a její obecné řešení ve tvaru
x(t) = Aλt
1 + Bλt
2.
Příklad: Najdeme obecné řešení rovnice
x(t + 2) − 4x(t + 1) + 3x(t) = 0.
Charakteristická rovnice
λ2
− 4λ + 3 = 0
má kořeny λ1 = 3, λ2 = 1, takže obecné řešení dané rovnice je
x(t) = 3t
A + B.
Komplexně sdružené charakteristické kořeny
Tento případ nastává, pokud α2
1 < 4α0. Charakteristické kořeny jsou
λ1,2 = 1
2 −α1 ± i
√
4α0 − α1 ,
v goniometrickém tvaru
λ1,2 = re±iω
= r (cos ω ± i sin ω) ,
kde
r =
√
α0, ω = arccos −
α1
2
√
α0
, ω ∈ (0, π). (2.19)
Dostáváme tak dvě řešení rovnice (2.16) ve tvaru
x1(t) = rt
eiωt
= rt
(cos ωt + i sin ωt), x2(t) = rt
eiωt
= rt
(cos ωt − i sin ωt).
2.2. ROVNICE DRUHÉHO ŘÁDU 17
Lineární kombinace těchto posloupností jsou podle principu superpozice také řešením rovnice
(2.16), zejména tedy posloupnosti
y1(t) =
1
2
x1(t) +
1
2
x2(t) = rt
cos ωt, y2(t) =
1
2i
x1(t) −
1
2i
x2(t) = rt
sin ωt
jsou řešením. Jejich Casoratián
C(t; y1, y2) =
rt cos ωt rt sin ωt
rt+1 cos ω(t + 1) rt+1 sin ω(t + 1)
=
= rt+2
cos ωt sin ω(t + 1) − sin ωt cos ω(t + 1) =
= rt+2
cos ωt(sin ωt cos ω + sin ω cos ωt) − sin ωt(cos ωt cos ω − sin ω sin ωt) =
= rt+2
(cos ωt)2
+ (sin ωt)2
sin ω = rt+2
sin ω = 0,
poněvadž r > 0 a 0 < ω < π. Posloupnosti y1, y2 jsou tedy lineárně nezávislé, tvoří fundamentální
systém řešení. Obecné řešení rovnice (2.16) proto můžeme psát ve tvaru
x(t) = rt
(A cos ωt + B sin ωt), (2.20)
kde r a ω jsou dány rovnostmi (2.19).
Příklad: Najdeme řešení počáteční úlohy
x(t + 2) − 4x(t + 1) + 8x(t) = 0, x(0) = 1, x(1) = −
√
2.
Charakteristická rovnice
λ2
− 4λ + 8 = 0
má komplexně sdružené kořeny
λ1,2 = 2 ± 2i = 2 cos π
4 ± i sin π
4 .
Fundamentální systém řešení dané rovnice je
y1(t) = 2t
cos π
4 t, y2(t) = 2t sin π
4 t
a její obecné řešení je tvaru
x(t) = 2t
A cos π
4 t + B sin π
4 t .
Z počáteční podmínky dostaneme soustavu dvou rovnic
A = 1,√
2A+
√
2B = −
√
2
a z nich hodnoty konstant A = 1, B = −2. Řešení dané úlohy je tedy
x(t) = 2t
cos π
4 t − 2 sin π
4 t = 2t
cos π
4 t − 2t+1 sin π
4 t.
Tvar obecného nenulového řešení rovnice (2.16) s α2
1 < 4α0 můžeme také upravit na tvar
x(t) = rt
A2 + B2
A
√
A2 + B2
cos ωt +
B
√
A2 + B2
sin ωt .
Označíme C =
√
A2 + B2, ϕ = arctg A
B a obecné řešení rovnice (2.16) zapíšeme ve tvaru
x(t) = Crt
(sin ϕ cos ωt + cos ϕ sin ωt) = Crt
sin(ωt + ϕ).
18 KAPITOLA 2. LINEÁRNÍ ROVNICE
Dvojnásobný reálný charakteristický kořen
Tento případ nastává, pokud α2
1 = 4α0. Dvojnásobný kořen charakteristické rovnice (2.18) je
λ = −1
2α1. Jedna složka fundamentálního systému řešení rovnice (2.16) je tedy dána výrazem
y1(t) = −
α1
2
t
.
Druhá složka je podle Tvrzení 3 tvaru
y2(t) = −
α1
2
t
t−1
i=t0
i−1
j=t0
α0
−1
2α1
j
−1
2α1
j+2
= −
α1
2
t
t−1
i=t0
i−1
j=t0
α0 −
α1
2
−2
=
= −
α1
2
t
t−1
i=t0
i−1
j=t0
4
α0
α2
1
= −
α1
2
t
t−1
i=t0
1i−t0
= (t − t0) −
α1
2
t
.
Obecné řešení rovnice (2.16) je tedy v případě α2
1 = 4α0 rovno
x(t) = A −
α1
2
t
+ Bt −
α1
2
t
= (A + Bt) −
α1
2
t
.
Příklad: Najdeme řešení počáteční úlohy
x(t + 2) − 4x(t + 1) + 4x(t) = 0, x(0) = 1, x(1) = 0.
Charakteristická rovnice
λ2
− 4λ + 4 = 0
má dvojnásobný kořen λ = 2. Obecné řešení dané rovnice tedy je
x(t) = 2t
A + 2t
tB.
Z počáteční podmínky dostaneme soustavu rovnic
A = 1,
2A+2B = 0,
která má řešení A = 1, B = −1. Řešení dané úlohy tedy je
x(t) = 2t
(1 − t).
2.2.4 Nehomogenní rovnice, jejíž pravá strana má konstantní koeficienty
Přidružená homogenní rovnice k nehomogenní lineární rovnici druhého řádu
x(t + 2) + α1x(t) + α0x(t) = b(t) (2.21)
je rovnice (2.16), kterou jsme se zabývali v předchozí části. Můžeme tedy napsat její fundamentální
systém řešení a pak řešení nehomogenní rovnice (2.21) najít pomocí formule variace
konstant (2.15).
2.2. ROVNICE DRUHÉHO ŘÁDU 19
Příklad: Najdeme řešení počáteční úlohy
x(t + 2) − 2x(t + 1) + x(t) = t, x(0) = 1, x(1) = −1.
Přidružená homogenní rovnice
x(t + 2) − 2x(t + 1) + x(t) = 0
má charakteristickou rovnici λ2 −2λ+1 = 0, která má dvojnásoný kořen 1. Obecné řešení přidružené
homogenní rovnice je tedy tvaru xH (t) = A+Bt, řešení splňující počáteční podmínku
je
xH (t) = 1 − 2t.
Partikulární řešení nehomogenní rovnice splňující nulové počáteční podmínky je podle formule
variace konstant (2.15)
xN (t) =
t−1
j=0
j
(j + 1) − t
(j + 1) − (j + 2)
= −
t−1
j=1
j(j + 1) − jt = t
t−1
j=1
j −
t−1
j=1
j(j + 1) =
= 1
2 t2(t − 1) − 1
3(t + 1)t(t − 1) = 1
6 t(t − 1)(t − 2).
Celkem tak dostáváme řešení dané úlohy ve tvaru
x(t) = 1 − 2t + 1
6t(t − 1)(t − 2).
Příklad: Najdeme obecné řešení nehomogenní rovnice
x(t + 2) + x(t) = f(t)
s obecnou pravou stranou. K přidružené homogenní rovnici
x(t + 2) + x(t) = 0
Přísluší charakteristická rovnice λ2+1 = 0, která má ryze imaginární kořeny λ1,2 = ±i = e±i π
2 .
Fundamentální systém řešení přidružené homogenní rovnice je tedy
y1(t) = cos π
2 t, y2(t) = sin π
2 t,
obecné řešení této rovnice je
xH(t) = A cos π
2 t + B sin π
2 t.
Řešení nehomogenní rovnice je dáno formulí variace konstant (2.15),
xN (t) =
t−1
j=t0
f(j)
sin π
2 (j + 1) cos π
2 t − cos π
2 (j + 1) sin π
2 t
cos π
2 (j + 2) sin π
2 (j + 1) − cos π
2 (j + 1) sin π
2 (j + 2)
.
Platí
sin π
2 (j + 1) = sin π
2 j cos π
2 + sin π
2 cos π
2 j = cos π
2 j,
20 KAPITOLA 2. LINEÁRNÍ ROVNICE
cos π
2 (j + 1) = cos π
2 j cos π
2 − sin π
2 sin π
2 j = − sin π
2 j,
takže
sin π
2 (j + 2) = cos π
2 (j + 1) = − sin π
2 j, cos π
2 (j + 2) = − sin π
2 (j + 1) = − cos π
2 j,
sin π
2 (j + 1) cos π
2 t − cos π
2 (j + 1) sin π
2 t = cos π
2 j cos π
2 t + sin π
2 j sin π
2 t = cos π
2 (t − j),
cos π
2 (j + 2) sin π
2 (j + 1) − cos π
2 (j + 1) sin π
2 (j + 2) = − cos π
2 j
2
− sin π
2 j
2
= −1.
Odtud dostáváme řešení nehomogenní rovnice ve tvaru
xN = −
t−1
j=t0
f(j) cos π
2 (t − j).
Řešení dané rovnice tedy je
x(t) = A cos π
2 t + B sin π
2 t −
t−1
j=t0
f(j) cos π
2 (t − j).
Fundamentální systém řešení přidružené homogenní rovnice (2.16) k nehomogenní rovnici
(2.21) můžeme pomocí charakteristických kořenů napsat explicitně. Proto také můžeme napsat
pomocí sumace tvar partikulárního řešení nehomogenní rovnice, které splňuje nulové počáteční
podmínky:
Tvrzení 5. Partikulární řešení xN nehomogenní rovnice (2.21) splňující počáteční podmínky
x(t0) = 0, x(t0 + 1) = 0
je jednoho z následujících tvarů:
1. Má-li charakteristická rovnice (2.18) dva reálné různé kořeny λ1, λ2, pak
xN (t) =
λt−1
1
λ1 − λ2
t−1
j=t0
b(j)λ−j
1 −
λt−1
2
λ1 − λ2
t−1
j=t0
b(j)λ−j
2 =
=
λt−2
1
1 − (λ2/λ1)
t−1
j=t0
b(j)λ−j
1 1 −
λ2
λ1
t−j−1
.
2. Má-li charakteristická rovnice (2.18) komplexně sdružené kořeny r(cos ω ± i sin ω), pak
xN (t) =
rt−2
sin ω

sin ω(t − 1)
t−1
j=t0
b(j) cos ωj − cos ω(t − 1)
t−1
j=t0
b(j) sin ωj

 =
=
rt−2
sin ω
t−1
j=t0
b(j) sin ω(t − j − 1).
3. Má-li charakteristická rovnice (2.18) dvojnásobný reálný kořen λ = −1
2α1 = −
√
α0, pak
xN (t) =
1
λ2
t−1
j=t0
b(j)λt−j
(t − j − 1) =
1
α0
t−1
j=t0
b(j) −
α1
2
t−j
(t − j − 1).
Důkaz: Přímým výpočtem na základě výsledků z 2.2.3 a formule variace konstant (2.15).
2.2. ROVNICE DRUHÉHO ŘÁDU 21
2.2.5 Cauchyho-Eulerova rovnice
Jedná se o lineární rovnici druhého řádu ve tvaru
t(t + 1)x(t + 2) + α1tx(t + 1) + α0x(t) = b(t), (2.22)
kde α0 = 0. Tuto rovnici uvažujeme na oboru t > 0. Přidružená homogenní rovnice je
t(t + 1)x(t + 2) + α1tx(t + 1) + α0x(t) = 0 (2.23)
Rešení této rovnice budeme hledat ve tvaru
x(t) =
y(t)
(t − 1)!
kde y je zatím neurčená posloupnost. Při této substituci je
x(t + 1) =
y(t + 1)
t!
=
1
t
y(t + 1)
(t − 1)!
, x(t + 2) =
y(t + 2)
(t + 1)!
=
1
(t + 1)t
y(t + 2)
(t − 1)!
.
Po dosazení tohoto vyjádření do rovnice (2.23) dostaneme lineární diferenční rovnici
y(t + 2) + α1y(t + 1) + α0y(t) = 0
pro neznámou posloupnost y, což je lineární homogenní rovnice s konstantními koeficienty.
S pomocí faktoriálové posloupnost (viz Tab. 1.1) zapíšeme Cauchyho-Eulerovu rovnici a
používanou substituci v kratším tvaru
t(2)
x(t + 2) + α1t(1)
x(t + 1) + α0x(t) = b(t) a x(t) = 0(1−t)
y(t).
Příklad: Najdeme obecné řešení rovnice
t(t + 1)x(t + 2) − x(t) = 1.
Substituce
x(t) =
y(t)
(t − 1)!
ji převádí na rovnici
y(t + 2) − y(t) = (t − 1)!.
Charakteristická rovnice λ2 − 1 = 0 má dva reálné různé kořeny λ1,j = ±1. Podle Tvrzení 5
je partikulární řešení nehomogenní rovnice tvaru
yN (t) =
1
2
t−1
j=1
(j − 1)! 1 − (−1)t−j−1
=
1
2
t−2
j=0
j! 1 − (−1)t−j
,
takže obecné řešení transformované rovnice je
y(t) = A + (−1)t
B +
1
2
t−2
j=0
j! 1 − (−1)t−j
.
Zpětnou substitucí dostaneme řešení dané rovnice ve tvaru
x(t) =
A + (−1)tB
(t − 1)!
+
1
2
t−2
j=0
1 − (−1)t−j
(t − 1)(t − 2) · · · (j + 1)
.
22 KAPITOLA 2. LINEÁRNÍ ROVNICE
2.3 Lineární rovnice k-tého řádu
Jedná se o rovnici
∆k
x + ak−1(t)∆k−1
x + ak−2(t)∆k−2
(t)x + · · · + a1(t)∆x + a0(t) = b(t). (2.24)
O posloupnostech a0, a1, a2, . . . , ak−1, b předpokládáme, že mají stejný definiční obor, označíme
ho D, a pro každé t z tohoto definičního oboru platí
ak−1(t) − ak−2(t) + ak−3(t) − · · · + (−1)k−1
a0(t) = 1. (2.25)
V případě b ≡ 0 se rovnice (2.24) nazývá homogenní, v opačném případě nehomogenní.
Je-li t0 ∈ D, jsou počáteční podmínky pro rovnici (2.24) tvaru
x(t0) = ξ0, x(t0 + 1) = ξ1, . . . , x(t0 + k − 1) = ξk−1. (2.26)
Rovnici (2.24) přepíšeme na rovnici druhého typu:
∆k
x(t) = x(t + k) +
k
j=1
(−1)j k
j
x(t + k − j) = x(t + k) +
k−1
j=0
(−1)k−j k
j
x(t + j),
k−1
j=0
aj(t)∆j
x(t) =
k−1
j=0
aj(t)
j
i=0
(−1)i j
i
x(t+j −i) =
k−1
i=0
k−1
j=i
aj(t)(−1)i j
i
x(t+j −i) =
=
k−1
i=0
k−1−i
j=0
aj+i(t)(−1)i j + i
i
x(t + j) =
k−1
j=0
k−1−i
i=0
ai+j(t)(−1)i j + i
i
x(t + j),
takže levá strana rovnice (2.24) je tvaru
x(t + k) +
k−1
j=0
(−1)k−j k
j
+
k−1−j
i=0
ai+j(t)(−1)i j + i
i
x(t + j).
Označíme
cj(t) = (−1)k−j k
j
+
k−1−j
i=0
ai+j(t)(−1)i j + i
i
pro j = 0, 1, 2, . . . , k − 1
a dostaneme rovnici druhého typu ekvivalentní s rovnicí (2.24) ve tvaru
x(t+k)+ck−1(t)x(t+k−1)+ck−2(t)x(t+k−2)+· · ·+c1(t)x(t+1)+c0(t)x(t) = b(t); (2.27)
podmínka (2.25) zaručí, že c0(t) = 0 pro všechna t ∈ D, takže se skutečně jedná o rovnici ktého
řádu. Z tvaru rovnice (2.27) vidíme, že počáteční úloha (2.27), (2.26), nebo ekvivalentně
úloha (2.24), (2.26), má jediné řešení, které je definováno na množině D.
2.3. LINEÁRNÍ ROVNICE K-TÉHO ŘÁDU 23
2.3.1 Fundamentální systém řešení homogenní rovnice
Lineární homogenní diferenční rovnice k-tého řádu
x(t + k) + ck−1(t)x(t + k − 1) + ck−2(t)x(t + k − 2) + · · · + c1(t)x(t + 1) + c0(t)x(t) = 0 (2.28)
splňuje princip superpozice: jsou-li posloupnosti x1 a x2 řešení rovnice (2.28) a p a q jsou
libovolné reálné konstanty, pak také posloupnost x = px1 + qx2 je řešením rovnice (2.41), tj.
libovolná lineární kombinace řešení této rovnice je jejím řešením. Navíc nulová posloupnost
x ≡ 0 je řešením rovnice (2.41). To znamená, že množina všech řešení lineární homogenní
diferenční rovnice tvoří vektorový prostor.
Pro i ∈ {0, 1, 2, . . . , k − 1} označme yi posloupnost, která je řešením homogenní rovnice
(2.28) s počátečními podmínkami
x(t0 + j) =
1, j = i,
0, j = i,
j = 0, 1, 2, . . . , k − 1.
Pak je zřejmé, že posloupnosti y0, y1, . . . , yk−1 jsou lineárně nezávislé. To znamená, že dimenze
vektorového prostoru řešení je alespoň k.
Nechť y je libovolné řešení homogenní rovnice (2.28). Označme
η0 = y(t0), η1 = y(t0 + 1), . . . , ηk−1 = y(t0 + k − 1).
Lineární kombinace posloupností y0, y1, . . . , yk−1 s koeficienty η0, η1, . . . , ηk−1, tj. posloupnost
η0y0 + η1y1 + · · · + ηk−1yk−1 (2.29)
je podle principu superpozice řešením rovnice (2.28) a splňuje stejné počáteční podmínky,
jako posloupnost y. Z jednoznačnosti řešení počáteční úlohy plyne, že posloupnost y a lineární
kombinace (2.29) jsou shodné. Odtud dále plyne, že prostor řešení lineární homogenní rovnice
(2.28) má dimenzi k a posloupnosti yi, i = 0, 1, 2, . . . , k − 1 tvoří bázi tohoto prostoru.
Z provedených úvah plyne, že platí
Věta 1. Množina všech řešení lineární homogenní difereneční rovnice k-tého řádu (2.28)
tvoří vektorový prostor dimenze k.
Definice 1. Báze vektorového prostoru všech řešení lineární homogenní rovnice (2.28) se
nazývá fundamentální systém řešení.
Posloupnosti z1, z2, . . . , zk tvoří fundamentální systém řešení lineární homogenní diferenční
rovnice (2.28) právě tehdy, když libovolné řešení x této rovnice lze vyjádřit jako jejich lineární
kombinaci, tj. právě tehdy, když existují jednoznačně určené konstanty A1, A2, . . . , Ak takové,
že
x(t) = A1z1(t) + A2z2(t) + · · · + Akzk(t) (2.30)
pro libovolné t z definičního oboru D. Předchozí rovnost je ekvivalentní s rovnostmi
A1z1(t) + A2z2(t) + · · · + Akzk(t) = ξ0
A1z1(t + 1) + A2z2(t + 1) + · · · + Akzk(t + 1) = ξ1
...
...
...
...
A1z1(t + k − 1) + A2z2(t + k − 1) + · · · + Akzk(t + k − 1) = ξk−1
(2.31)
24 KAPITOLA 2. LINEÁRNÍ ROVNICE
a jednoznačná existence konstant A1, A2, . . . , Ak je ekvivalentní s jednoznačnou řešitelností
(2.31) chápané jako systém (algebraických) rovnic pro neznámé A1, A2, . . . , Ak. Determinant
této soustavy je Casoratián posloupností z1, z2, . . . , zk v indexu t,
C(t; z1, z2, . . . , zk) =
z1(t) z2(t) . . . zk(t)
z1(t + 1) z2(t + 1) . . . zk(t + 1)
...
...
...
...
z1(t + k − 1) z2(t + k − 1) . . . zk(t + k − 1)
.
Dostáváme tak závěr:
Věta 2. Posloupnosti z1, z2, . . . , zk tvoří fundamentální systém řešení lineární homogenní
diferenční rovnice (2.28) právě tehdy, když každá z nich je řešením rovnice (2.28) a pro každé
t z definičního oboru D platí
C(t; z1, z2, . . . , zk) = 0,
kde C(t; z1, z2, . . . , zk) je Casoratián posloupností z1, z2, . . . , zk.
2.3.2 Nehomogenní rovnice a metoda variace konstant
Pokud jsou posloupnosti c0, c1, . . . , ck−1 v rovnicích (2.27) a (2.28) stejné, řekneme, že homogenní
lineární diferenční rovnice (2.28) je přidružená k nehomogenní rovnici (2.27).
Je-li posloupnost y řešením nehomogenní rovnice (2.27) a posloupnost z je řešením přidružené
homogenní rovnice (2.28), pak jejich součet x = z + y je opět řešením nehomogenní
rovnice (2.27), neboť
x(t + k) +
k
i=1
ci(t)x(t + i) = z(t + k) + y(t + k) +
k
i=1
ci(t) z(t + i) + y(t + i) =
= z(t + k) +
k
i=1
ci(t)z(t + i) + y(t + k) +
k
i=1
ci(t)y(t + i) = 0 + b(t) = b(t).
Platí tedy
Věta 3. Nechť z1, z2, . . . , zk tvoří fundamentální systém řešení lineární homogenní diferenční
rovnice (2.28) přidružené k nehomogenní rovnici (2.27). Pak každé řešení nehomogenní rovnice
(2.27) je tvaru
x(t) = B1z1(t) + B2z2(t) + · · · + Bkzk(t) + y(t),
kde y je nějaké řešení nehomogenní rovnice a B1, B2, . . . , Bk jsou konstanty.
Nechť posloupnosti z1, z2, . . . , zk tvoří fundamentální systém řešení homogenní rovnice
(2.28) přidružené k nehomogenní rovnici (2.28). Pak je
zi(t + k) = −
k−1
j=0
cj(t)zi(t + j) pro i = 1, 2, . . . , k. (2.32)
2.3. LINEÁRNÍ ROVNICE K-TÉHO ŘÁDU 25
Řešení nehomogenní rovnice (2.28) budeme hledat ve tvaru
x(t) =
k
i=1
ui(t)zi(t), (2.33)
kde u1, u2, . . . , uk jsou zatím neurčené posloupnosti. Hledáme ho tedy jako analogii řešení
homogenní rovnice (2.30); místo konstant A1, A2, . . . , Ak však píšeme posloupnosti — varírujeme
konstanty. Z tohoto důvodu se tato metoda řešení nehomogenní rovnice nazývá metoda
variace konstant.
Nyní můžeme vyjádřit
x(t + 1) =
k
i=1
ui(t + 1)zi(t + 1) =
k
i=1
∆ui(t) zi(t + 1) + ui(t)zi(t + 1) .
Budeme požadovat, aby posloupnosti u1, u2, . . . , uk splňovaly rovnici
k
i=1
∆ui(t) zi(t + 1) = 0.
Pak x(t + 1) =
k
i=1
ui(t)zi(t + 1), takže
x(t + 2) =
k
i=1
ui(t + 1)zi(t + 2) =
k
i=1
∆ui(t) zi(t + 2) + ui(t)zi(t + 2) .
Dále budeme požadovat, aby posloupnosti u1, u2, . . . , uk splňovaly rovnice
k
i=1
∆ui(t) zi(t + 2) = 0,
takže x(t + 2) =
k
i=1
ui(t)zi(t + 2). Takto budeme pokračovat až k požadavku
k
i=1
∆ui(t) zi(t + k − 1) = 0
a vyjádření x(t + k − 1) =
k
i=1
ui(t)zi(t + k − 1).
Celkem tedy požadujeme
k
i=1
∆ui(t) zi(t + j) = 0, j = 1, 2, . . . , k − 1 (2.34)
a dostáváme
x(t + j) =
k
i=1
ui(t)zi(t + j), j = 1, 2, . . . , k − 1. (2.35)
26 KAPITOLA 2. LINEÁRNÍ ROVNICE
V poslední z rovností (2.35), tj. v té, v níž j = k − 1, budeme psát t + 1 místo t a upravíme
ji s použitím (2.32). Dostaneme
x(t + k) =
k
i=1
ui(t + 1)zi(t + k) =
k
i=1
∆ui(t) zi(t + k) +
k
i=1
ui(t)zi(t + k) =
=
k
i=1
∆ui(t) zi(t + k) −
k
i=1
ui(t)
k−1
j=0
cj(t)zi(t + j). (2.36)
Současně posloupnost x má být řešením rovnice (2.27), takže s využitím vztahů (2.35) dosta-
neme
x(t + k) = b(t) −
k−1
j=0
cj(t)x(t + j) = b(t) −
k−1
j=0
cj(t)
k
i=1
ui(t)zi(t + j) =
= b(t) −
k
i=1
ui(t)
k−1
j=0
cj(t)zi(t + j). (2.37)
Porovnáním (2.36) a (2.37) vidíme, že
k
i=1
∆ui(t) zi(t + k) = b(t). (2.38)
Diference posloupností u1, u2, . . . , uk tedy splňují systém rovnic (2.34), (2.38). Přepíšeme ho
do tvaru
z1(t + 1) ∆u1(t) + z2(t + 1) ∆u2(t) + · · · + zk(t + 1) ∆uk(t) = 0
z1(t + 2) ∆u1(t) + z2(t + 2) ∆u2(t) + · · · + zk(t + 2) ∆uk(t) = 0
...
...
...
...
z1(t + k − 1) ∆u1(t) + z2(t + k − 1) ∆u2(t) + · · · + zk(t + k − 1) ∆uk(t) = 0
z1(t + k) ∆u1(t) + z2(t + k) ∆u2(t) + · · · + zk(t + k) ∆uk(t) = b(t).
Determinant této soustavy je Casoratiánem fundamentálního systému řešení homogenní rovnice
(2.28) v indexu t + 1. Je tedy nenulový a soustava je jednoznačně řešitelná. Označíme
w(t) =
z1(t) z2(t) . . . zk(t)
z1(t + 1) z2(t + 1) . . . zk(t + 1)
...
...
...
...
z1(t + k − 1) z2(t + k − 1) . . . zk(t + k − 1)
,
mi(t) =
=
z1(t) z2(t) . . . zi−1(t) zi+1(t) . . . zk(t)
z1(t + 1) z2(t + 1) . . . zi−1(t + 1) zi+1(t + 1) . . . zk(t + 1)
...
...
...
...
...
...
...
z1(t + k − 2) z2(t + k − 2) . . . zi−1(t + k − 2) zi+1(t + k − 2) . . . zk(t + k − 2)
;
2.3. LINEÁRNÍ ROVNICE K-TÉHO ŘÁDU 27
w(t) je Casoratián fundamentálního řešení homogenní rovnice (2.28). Diference posloupností
u1, u2, . . . uk nyní můžeme vyjádřit ve tvaru
∆ui(t) =
(−1)k+ib(t)mi(t + 1)
w(t + 1)
, i = 1, 2, . . . , k.
Odtud a z rovnosti (1.3) dostaneme
ui(t) = ui(t0) + (−1)k+i
t−1
j=t0
b(j)mi(j + 1)
w(j + 1)
, i = 1, 2, . . . , k.
Při označení Bi = ui(t0) můžeme řešení rovnice (2.27) podle vztahu (2.33) psát ve tvaru
x(t) =
k
i=1
Bizi(t) + (−1)k
t−1
j=t0
b(j)
w(j + 1)
k
i=1
(−1)i
mi(j + 1)zi(t).
2.3.3 Homogenní rovnice s konstantními koeficienty
Jedná se o rovnici
∆k
x + αk−1∆k−1
x + αk−2∆k−2
x + · · · + α1∆x + α0 = 0, (2.39)
kde reálné koeficienty α0, α1, . . . , αk−1 splňují rovnost
αk−1 − αk−2 + αk−3 − · · · + (−1)k−1
α0 = 1 (2.40)
analogickou k rovnosti (2.25). Rovnice (2.39) má pro libovolné počáteční podmínky tvaru
(2.26) jediné řešení, které je definováno pro každé t ∈ Z.
Stejně jako v případě obecné lineární rovnice k-tého řádu můžeme rovnici (2.39) přepsat
na rovnici druhého typu
x(t + k) + γk−1x(t + k − 1) + γk−2x(t + k − 2) + · · · + γ1x(t + 1) + γ0x(t) = 0, (2.41)
kde jsme označili
γj = (−1)k−j k
j
+
k−1−j
i=0
αi+j(−1)i j + i
i
pro j = 0, 1, 2, . . . , k − 1.
S pomocí operátoru posunu σ můžeme tuto rovnici přepsat do tvaru
xσk
(t) + γk−1xσk−1
x(t) + γk−2xσk−2
x(t) + · · · + γ1xσ
(t) + γ0x(t) = 0.
Položíme-li γk = 1, můžeme operátorovou rovnici zapsat ještě stručněji
k
i=0
γi ·σi
x ≡ 0. (2.42)
Ze stejných důvodů jako v odstavci 2.2.3 budeme řešení rovnice (2.41) hledat ve tvaru
x(t) = λt, kde λ je zatím neurčená nenulová konstanta. Dosadíme tuto posloupnost do rovnice
(2.41)
λt+k
+ γk−1λt+k−1
+ γk−2λt+k−2
+ · · · + γ1λt+1
+ γ0λt
= 0
28 KAPITOLA 2. LINEÁRNÍ ROVNICE
a po vynásobení výrazem λ−t dostaneme charakteristickou rovnici
λk
+ γk−1λk−1
+ γk−2λk−2
+ · · · + γ1λ + γ0 = 0. (2.43)
Řešení této algebraické rovnice se nazývají charakteristické kořeny. Povšimněme si, že žádný
kořen rovnice (2.43) není nulový, neboť γ0 = 0.
Věta 4. Nechť λp je r-násobný kořen charakteristické rovnice (2.43). Pak každá z posloupností
definovaných vztahem
x(t) = tq
λt
p, q = 0, 1, 2, . . . , r − 1
je řešením lineární homogenní diferenční rovnice (2.41).
Důkaz: Položíme γk = 1 a polynom na levé straně rovnice (2.43) označíme P(λ), tj.
P(λ) =
k
i=1
γiλi
.
Nejprve dokážeme pomocné tvrzení: Ke každému přirozenému číslu s a každému přirozenému
číslu j ∈ {0, 1, 2, . . . , s} existuje polynom ps,j stupně nejvýše s ve dvou proměnných t, λ
takový, že
k
i=0
γi(t + i)s
λi
=
s
j=0
ps,j(t, λ)P(j)
(λ).
Tvrzení dokážeme úplnou indukcí vzhledem k proměnné s. Pro s = 0 je
k
i=0
γi(t + i)0
λi
=
k
i=0
γiλi
= P(λ) = 1P(0)
(λ),
tedy p0,0 ≡ 1.
Indukční krok je obsažen ve výpočtu:
k
i=0
γi(t + i)s+1
λi
=
k
i=0
γi(t + i)s
(t + i)λi
= t
k
i=0
γi(t + i)s
λi
+ λ
k
i=1
γi(t + i)s
iλi−1
=
= t
s
j=0
ps,j(t, λ)P(j)
(λ) + λ
k
i=1
γi(t + i)s d
dλ
λi
=
= t
s
j=0
ps,j(t, λ)P(j)
(λ) + λ
d
dλ
k
i=1
γi(t + i)s
λi
=
= t
s
j=0
ps,j(t, λ)P(j)
(λ) + λ
d
dλ
k
i=0
γi(t + i)s
λi
− γ0ts
=
= t
s
j=0
ps,j(t, λ)P(j)
(λ) + λ
d
dλ


s
j=0
ps,j(t, λ)P(j)
(λ) − γ0ts

 =
2.3. LINEÁRNÍ ROVNICE K-TÉHO ŘÁDU 29
= t
s
j=0
ps,j(t, λ)P(j)
(λ) + λ
s
j=0
∂ps,j(t, λ)
∂λ
P(j)
(λ) + ps,j(t, λ)P(j+1)
(λ) =
=
s
j=0
tps,j(t, λ) + λ
∂ps,j(t, λ)
∂λ
P(j)
(λ) + λ
s+1
j=1
ps,j−1(t, λ)P(j)
(λ) =
= tps,0(t, λ) + λ
∂ps,0(t, λ)
∂λ
P(λ)+
+
s
j=1
tps,j(t, λ) + λ
∂ps,j(t, λ)
∂λ
+ λps,j−1(t, λ) P(j)
(λ)+
+ λps,s(t, λ)P(s+1)
(λ).
Stačí tedy položit
ps+1,0(t, λ) = tps,0(t, λ) + λ
∂ps,0(t, λ)
∂λ
,
ps+1,j(t, λ) = tps,j(t, λ) + λ
∂ps,j(t, λ)
∂λ
+ λps,j−1(t, λ) pro j = 1, 2, . . . , s,
ps+1,s+1(t, λ) = λps,s(t, λ)
a pomocné tvrzení je dokázáno.
Nechť nyní λp je r-násobný kořen charakteristické rovnice. Pak je také kořenem derivací
polynomu P až do řádu r − 1, tj.
P(j)
(λp) = 0 pro j = 0, 1, 2, . . . , r − 1.
Nyní pro x(t) = tqλt
p, q ∈ {0, 1, 2, . . . , r − 1}, podle pomocného tvrzení platí
k
i=0
γix(t + i) =
k
i=0
γi(t + i)q
λt+i
p = λt
k
i=0
γi(t + i)q
λi
p = λt
q
j=0
pq,j(t, λ)P(j)
(λ) = 0.
Důsledek 1. Nechť λc = a(cos ϕ + i sin ϕ) je r-násobný komplexní kořen charakteristické
rovnice (2.43). Pak každá z posloupností definovaných některým ze vztahů
x1(t) = tq
at
cos tϕ, x2(t) = tq
at
sin tϕ, q = 0, 1, 2, . . . , r − 1
je řešením lineární homogenní diferenční rovnice (2.41).
Důkaz: Poněvadž polynom na levé straně rovnice (2.43) má reálné koeficienty, je také komplexně
sdružené číslo ¯λc = a(cos ϕ − i sin ϕ) kořenem charakteristické rovnice (2.43) a má
stejnou násobnost r. Podle Věty 4 (v níž jsme nepředpokládali, že by kořen charakteristické
rovnice byl reálný), je každá z posloupností definovaných vztahem
˜x1(t) = tq
at
(cos tϕ + i sin tϕ), ˜x2(t) = tq
at
(cos tϕ − i sin tϕ), q = 0, 1, 2, . . . , r − 1
30 KAPITOLA 2. LINEÁRNÍ ROVNICE
řešením rovnice (2.41). Podle principu superpozice jsou také posloupnosti
x1(t) =
1
2
˜x1(t) + ˜x2(t) , x2(t) =
1
2i
˜x1(t) − ˜x2(t)
řešením této rovnice.
Důsledek 2. Každému reálnému r-násobnému charakteristickému kořenu λ odpovídá r řešení
lineární homogenní rovnice (2.41)
λt
, tλt
, t2
λt
, . . . , tr−1
λt
a každému komplexnímu r-násobnému charakteristickému kořenu a(cos ϕ + i sin ϕ) odpovídá
2r řešení lineární homogenní rovnice (2.41)
at
cos tϕ, tat
cos tϕ, t2
at
cos tϕ, . . . , tr−1
at
cos tϕ,
at
sin tϕ, tat
sin tϕ, t2
at
sin tϕ, . . . , tr−1
at
sin tϕ.
Důsledek 3. Nechť
λ1, λ2, . . . , λk1
jsou všechny jednoduché reálné různé charakteristické kořeny,
λk1+1, λk1+2, . . . , λk2
jsou všechny reálné různé charakteristické kořeny, které mají násobnosti rk1+1, rk1+2, . . . , rk2
(v tomto pořadí) a
ak2+1(cos ϕk2+1 + i sin ϕk2+1), ak2+2(cos ϕk2+2 + i sin ϕk2+2), . . . , ak3 (cos ϕk3 + i sin ϕk3 )
jsou všechny komplexní charakteristické kořeny takové, že žádné dva z nich nejsou komplexně
sdružené a mají násobnosti rk2+1, rk2+2, . . . , rk3 (v tomto pořadí). Přitom samozřejmě platí
k1 +
k2
i=k1+1
ri + 2
k3
i=k2+1
ri = k.
Pak posloupnost definovaná vztahem
x(t) =
k1
i=1
Aiλt
i +
k2
i=k1+1
ri−1
j=0
Bijtj
λt
i +
k3
i=k2+1
ri−1
j=0
Cijtj
at
i cos tϕi +
k3
i=k2+1
ri−1
j=0
Dijtj
at
i sin tϕi,
(2.44)
kde Ai, Bij, Cij, Dij jsou konstanty, je řešením lineární homogenní rovnice (2.41).
Nechť existuje charakteristický kořen, jehož modul (absolutní hodnota) je větší, než moduly
všech ostatních charakteristických kořenů. Takový charakteristický kořen musí být reálný
a jednoduchý, můžeme ho tedy označit λ1. Platí
|λ1| > |λi| pro i = 2, 3, . . . , k2, |λ1| > ai pro i = k2 + 1, k2 + 2, . . . , k3.
2.3. LINEÁRNÍ ROVNICE K-TÉHO ŘÁDU 31
Charakteristický kořen λ1 s těmito vlastnostmi nazveme ryze dominantní. Nyní pro řešení
x(t) rovnice (2.41) definované vztahem (2.44) za předpokladu A1 = 0 platí
lim
t→∞
x(t)
A1λt
1
= lim
t→∞

1 +
k1
i=2
Ai
A1
λi
λ1
t
+
k2
i=k1+1
ri−1
j=0
Bij
A1
tj λi
λ1
t
+
+
k3
i=k2+1
ri−1
j=0
Cij
A1
tj ai
λ1
t
cos tϕi +
k3
i=k2+1
ri−1
j=0
Dij
A1
tj ai
λ1
t
sin tϕi

 = 1.
Dostáváme tak
Důsledek 4. Pokud existuje ryze dominantní charakteristický kořen λ1 a konstanta A1 v řešení
(2.44) rovnice (2.41) je nenulová, pak toto řešení je asymptoticky ekvivalentní s geometrickou
posloupností s kvocientem λ1.
Řekneme, že charakteristický kořen je dominantní, pokud jeho modul není menší než
modul jakéhokoliv charakteristického kořene, tj. dominantní charakteristický kořen má maximální
modul. Označme tento maximální modul symbolem Λ.
Nechť jsou charakteristické kořeny označeny jako v Důsledku 3 a navíc platí
|λ1| ≥ |λ2| > |λ3| ≥ |λ4| ≥ · · · ≥ |λk1 |,
|λk1+1| ≥ |λk1+2| ≥ · · · ≥ |λk2 |, ak2+1 ≥ ak2+2 ≥ · · · ≥ ak3 .
Položme
l1 =



2, Λ = |λ2|,
1, Λ = |λ1| > |λ2|,
0, Λ > |λ1|,
l2 =
max i ∈ {k2 + 1, k2 + 2, . . . , k3} : ai = Λ , Λ = ak2+1,
k2, Λ > ak2+1.
Nechť dominantní charakteristické kořeny jsou jednoduché, tj. |λk1+1| < Λ a pokud l2 > k2
tak max ri : i ∈ {k2 + 1, k2 + 2, . . . , l2} = 1. Označme
y(t) =
l1
i=1
Ai(sgn λi)t
+
l2
i=k2+1
(Ci0 cos tϕi + Di0 sin tϕi).
Pak
x(t)
Λt
− y(t) =
k1
i=l1+1
Ai
λi
Λ
t
+
k2
i=k1+1
ri−1
j=0
Bijtj λi
Λ
t
+
+
k3
i=l2+1
ri−1
j=0
Cijtj ai
Λ
t
cos tϕi +
k3
i=l2+1
ri−1
j=0
Dijtj ai
Λ
t
sin tϕi.
Limita pro t → ∞ posloupnosti na pravé straně této rovnosti je rovna 0. To — zhruba řečeno
— znamená, že „pro dostatečně velké t se řešení rovnice (2.41) chová jako posloupnost y .
32 KAPITOLA 2. LINEÁRNÍ ROVNICE
Poněvadž pro libovolné t platí nerovnosti
−∞ < −
l1
i=1
|Ai|−
l2
i=k2+1
|Ci0|+|Di0| ≤ y(t) ≤ ∞ < −
l1
i=1
|Ai|+
l2
i=k2+1
|Ci0|+|Di0| < ∞,
je
−∞ < m = lim inf
t→∞
y(t) ≤ lim sup
t→∞
= M < ∞,
pro řešení x(t) rovnice (2.41) definované rovností (2.44) platí
mΛt
≤ x(t) ≤ MΛt
.
2.3.4 Rovnice s konstantními koeficienty a speciální pravou stranou
Uvažujme nehomogenní lineární diferenční rovnici k-tého řádu druhého typu
x(t + k) + γk−1x(t + k − 1) + · · · + γ1x(t + 1) + γ0x(t) = b(t) (2.45)
a k ní přidruženou lineární homogenní rovnici (2.41). Označme polynomiální operátor posunu
z levé strany operátorové rovnice (2.42) symbolem PΣ; homogenní rovnici (2.41) tedy můžeme
zapsat ve tvaru
PΣ
x(t) ≡ 0 nebo PΣ
x = 0,
a nehomogenní rovnici ve tvaru
PΣ
x(t) = b(t) nebo PΣ
x = b.
Definice 2. Nechť p ∈ P je posloupnost, β0, β1, . . . , βl ∈ R jsou konstanty takové, že β0 =
0 = βl, a nechť RΣ je polynomiální operátor posunu, RΣ =
l
i=0
βi ·σi
. Řekneme, že operátor
RΣ je anihilátor posloupnosti p, pokud
RΣ
p ≡ 0.
Podle této terminologie je PΣ anihilátorem každého řešení homogenní rovnice (2.41).
Nechť existuje anihilátor QΣ =
l
i=0
βi ·σi
posloupnosti b, která je na pravé straně nehomogenní
rovnice (2.45). To znamená, že b je řešením nějaké lineární homogenní rovnice s konstatntními
koeficienty, takže podle Důsledku 2 je posloupnost b lineární kombinací výrazů
κt, tmκt, cos tψ, sin tψ, tn cos tψ, tn sin tψ, tmκt cos tψ, tmκt sin tψ. Nechť dále y je řešením
nehomogenní rovnice (2.45). Pak platí
QΣ
PΣ
y ≡ 0. (2.46)
To znamená, že řešení nehomogenní lineární rovnice k-tého řádu je současně řešením lineární
rovnice (k + l)-tého řádu.
Nechť λ1, λ2, . . . , λp, p ≤ k, jsou charakteristické kořeny homogenní rovnice
PΣ
x ≡ 0
2.3. LINEÁRNÍ ROVNICE K-TÉHO ŘÁDU 33
b(t) tvar řešení
at C1at
tm C0 + C1t + C2t2 + · · · + Cmtm
tmat at C0 + C1t + C2t2 + · · · + Cmtm
sin ψt, cos ψt C1 sin ψt + C2 cos ψt
at sin ψt, at cos ψt at (C1 sin ψt + C2 cos ψt)
attm sin ψt, attm cos ψt at [(C0 + C1t + · · · + Cmtm) sin ψt+
+ (D0 + D1t + · · · + Dmtm) cos ψt]
Tabulka 2.1: Tvary řešení nehomogenní rovnice (2.45) pro různé pravé strany b.
a µ1, µ2, . . . , µq, q ≤ l, jsou charakteristické kořeny homogenní rovnice
QΣ
x ≡ 0.
Nyní rozlišíme dva případy.
Případ 1: {λ1, λ2, . . . , λp} ∩ {µ1, µ2, . . . , µq} = ∅. V tomto případě můžeme psát řešení
nehomogenní rovnice podle tabulky 2.1. Takové obecně zapsané řešení dosadíme do rovnice
(2.45) a vypočítáme konstanty Cj, Dj.
Případ 2: {λ1, λ2, . . . , λp}∩{µ1, µ2, . . . , µq} = ∅. V tomto případě nejprve najdeme obecné
řešení homogenní rovnice (2.46) a vynecháme v něm všechny členy, které se vyskytují v obecném
řešení přidružené homogenní rovnice (2.41). Tím dostaneme řešení nehomogenní rovnice
(2.45) s neurčitými koeficienty, které určíme dosazením do původní rovnice (2.45).
Příklad:
Najdeme obecné řešení nehomogenní rovnice druhého řádu
x(t + 2) − 4x(t + 1) + 3x(t) = t2
. (2.47)
Charakteristická rovnice přidružené homogenní rovnice
λ2
− 4λ + 3 = 0,
má množina charakteristických kořenů Λ = {1, 3}, její obecné řešení je proto posloupnost
daná předpisem
x(t) = A · 3t
+ B.
Anihilátor řešení přidružené homogenní rovnice je polynomiální operátor posunu
PΣ
= ·σ2
− 4 ·σ
+3 idP
(koeficienty v operátoru PΣ jsou stejné, jako koeficienty v charakteristickém polynomu).
Pravá strana dané rovnice je řešením lineární homogenní rovnice, která má trojnásobný
charakteristický kořen µ = 1. Nejjednodušší taková algebraická rovnice je (µ − 1)3 = 0.
Posloupnost na pravé straně rovnice (2.47) je tedy řešením lineární homogenní diferenční
rovnice
x(t + 3) − 3x(t + 2) + 3x(t + 1) − x(t) = 0, (2.48)
její anihilátor je
QΣ
= (·σ
− idP)3
= ·σ3
− 3 ·σ2
+3 ·σ
− idP .
34 KAPITOLA 2. LINEÁRNÍ ROVNICE
Množina charakteristických kořenů rovnice (2.48) je jednoprvková M = {1}. To znamená, že
Λ ∩ M = {1} = ∅, tj. nastává Případ 2.
Posloupnost na levé straně relace (2.46) je v tomto konkrétním případě dána výrazem
QΣ
PΣ
y (t) = QΣ
y(t + 2) − 4y(t + 1) + 3y(t) =
= y(t + 5) − 4y(t + 4) + 3y(t + 3) − 3 y(t + 4) − 4y(t + 3) + 3y(t + 2) +
+3 y(t + 3) − 4y(t + 2) + 3y(t + 1) − y(t + 2) − 4y(t + 1) + 3y(t) =
= y(t + 5) − 7y(t + 4) + 18y(t + 3) − 22y(t + 2) + 13y(t + 1) − 3y(t).
Homogenní rovnice
y(t + 5) − 7y(t + 4) + 18y(t + 3) − 22y(t + 2) + 13y(t + 1) − 3y(t) = 0 (2.49)
má charakteristickou rovnmici
λ5
− 7λ4
+ 18λ3
− 22λ2
+ 13λ − 3 = 0,
po úpravě
(λ − 3)(λ − 1)4
= 0.
Tato rovnice má jednoduchý kořen λ = 3 a čtyřnásobný kořen λ = 1, obecné řešení homogenní
rovnice (2.49) tedy je posloupnost daná předpisem
y(t) = A · 3t
+ B + Ct + Dt2
+ Et3
. (2.50)
Vynecháním členů, které se vyskytují v obecném řešení přidružené homogenní rovnice k rovnici
zadané dostaneme řešení nehomogenní rovnice (2.47) ve tvaru s neurčitými koeficienty
Ct + Dt2
+ Et3
,
které určíme dosazením do dané rovnice:
C(t + 2) + D(t + 2)2
+ E(t + 2)3
− 4 C(t + 1) + D(t + 1)2
+ E(t + 1)3
+ 3 Ct + Dt2
+ Et3
=
= −6Et2
− 4Dt − 2C + 4E = 2t2
.
Porovnáním koeficientů dostaneme soustavu lineárních algebraických rovnic
−6E = 2
D = 0
−C +2E = 0,
která má řešení E = −1
3, D = 0, C = −2
3. Celkem dostáváme obecné řešení rovnice (2.47) ve
tvaru
x(t) = A · 3t
+ B − 2
3t − 1
3 t3.
Ještě poznamenejme, že uvedený postup měl především ilustrovat obecnou teorii. Pro
praktické počítání je příliš zdlouhavý; zejména vyjadřování příslušných anihilátorů není pro
nalezení výsledků nutné. Bezprostředně lze totiž psát, že řešení nehomogenní rovnice (2.47)
je řešením lineární homogenní rovnice pátého řádu, jejíž charakteristická rovnice je tvaru
(λ − 1)3
λ2
− 4λ + 3 = 0.
2.3. LINEÁRNÍ ROVNICE K-TÉHO ŘÁDU 35
Poté napíšeme obecné řešení takové rovnice – což je posloupnost (2.50), dosadíme do rovnice
zadané a tak určíme tři z pěti konstant obecného řešení.
Stručně lze říci, že „speciální pravá strana rovnice (2.45) je taková posloupnost b, která
je řešením nějaké lineární homogenní rovnice. V takovém případě vezmeme všechny charakteristické
kořeny rovnice, jejímž řešením posloupnost b je (ty jsou vidět již z tvaru posloupnosti
b podle Důsledků 1–3 Věty 4), a všechny charakteristické kořeny lineární homogenní rovnice
přidružené k rovnici (2.45); kořeny bereme včetně násobností. Napíšeme obecné řešení homogenní
rovnice odpovídající všem těmto charakteristickým kořenům, dosadíme ho do dané
rovnice a určíme všechny konstanty, které určit lze.
Příklad:
Najdeme obecné řešení rovnice
x(t + 2) + x(t) = t cos 1
2πt.
Charakteristická rovnice přidružené homogenní rovnice λ2 + 1 = 0 má komplexně sdružené
kořeny λ1,2 = ±i = cos 1
2π ± i sin 1
2 π. Posloupnost na pravé straně dané rovnice je řešením
lineární homogenní rovnice, která má dvojnásobný kořen λ = i; abychom zůstali v reálném
oboru, vezmeme dvojici komplexně sdružených dvojnásobných kořenů λ = ±i. Máme tedy
celkem trojnásobné komplexně sdružené kořeny λ = ±i = cos 1
2π ± i sin 1
2 π a obecný tvar
řešení
x(t) = A cos 1
2πt + B sin 1
2πt + Ct cos 1
2 πt + Dt sin 1
2πt + Et2 cos 1
2πt + Ft2 sin 1
2 πt.
Tuto posloupnost dosadíme do dané rovnice. Poněvadž
cos 1
2 π(t + 2) = cos 1
2πt cos π − sin 1
2πt sin π = − cos 1
2πt,
sin 1
2 π(t + 2) = sin 1
2πt cos π + cos 1
2πt sin π = − sin 1
2πt,
dostaneme
x(t + 2) + x(t) = −A cos 1
2πt − B sin 1
2 πt − C(t + 2) cos 1
2 πt − D(t + 2) sin 1
2 πt−
− E(t2
+ 4t + 4) cos 1
2πt − F(t2 + 4t + 4) sin 1
2πt+
+ A cos 1
2πt + B sin 1
2πt + Ct cos 1
2πt + Dt sin 1
2πt + Et2 cos 1
2πt + Ft2 sin 1
2 πt =
= −(4Et + 2C + 4E) cos 1
2 πt − (4Ft + 2D + 4F) sin 1
2πt.
Porovnáním koeficientů dostáváme −4E = 1, 2C +4E = 0, 4F = 0, 2D+4F = 0, tedy C = 1
2 ,
E = −1
4, D = F = 0 a koeficienty A, B zůstávají neurčené. Obecné řešení dané rovnice je
dáno výrazem
x(t) = A cos 1
2πt + B sin 1
2πt + 1
2t cos 1
2π − 1
4 t2 cos 1
2π = A + 1
2 t − 1
4t2 cos 1
2πt + B sin 1
2 πt,
kde A, B jsou konstanty.
Pro nehomogenní lineární diferenční rovnice tvaru (2.45) s pravou stranou ve tvaru polynomu
můžeme zformulovat ještě jednodušší pravidlo pro hledání řešení: Nechť m je stupeň
polynomu na pravé straně rovnice (2.45) a n je násobnost charakteristického kořene λ = 1 přidružené
homogenní rovnice (samozřejmě, že může být n = 0, pokud charakteristická rovnice
nemá kořen 1). Pak řešení nehomogenní rovnice je tvaru polynomu stupně n + m.
36 KAPITOLA 2. LINEÁRNÍ ROVNICE
2.4 Cvičení
V úlohách 1–6 najděte řešení dané rovnice s počáteční podmínkou x(0) = 1.
1. x(t + 1) = (t + 1)x(t)
2. x(t + 1) = 3tx(t)
3. x(t + 1) = e2tx(t)
4. x(t + 1) = 1
2x(t) + 2
5. x(t + 1) = x(t) + et
6. x(t + 1) =
t
t + 1
x(t) + 4
7. Najděte obecné řešení rovnice (t + 1)x(t + 1) = tx(t) v prostorech posloupností P−∞ a
P1.
8. Nechť je na konci každého období do banky ukládána částka $ 200 a banka platí každé
období úrok 0,8%. Jaká je uložená částka po n obdobích.
9. Dluh $ 12 000 má být amortizována splátkami $ 380 na konci každého měsíce, plus jedna
závěrečná splátka menší. Úrok 12% p.a. je připisován každý měsíc. Určete čas splácení
(v měsících) a závěrečnou splátku.
10. Úvěr $ 80 000 má být splacen pravidelnými měsíčními splátkami. Úroková sazba je
10% p.a. Jaká je měsíční splátka, aby byl dluh splacen do 30 měsíců.
11. Hypotéka na 30let má úrokovou sazbu 8% p.a. Jste schopni splácet $ 1 000 měsíčně.
Kolik si můžete půjčit?
12. Pokud organismus žije, je v jeho tkáních zastoupení radioaktivního uhlíku 14C stejné
jako v atmosféře. Po uhynutí organismu se radioaktivní uhlík rozkládá s poločasem
rozpadu 5 700 let. Ve vzorku je 14C zastoupen ze 70% ve srovnání s atmosférou. Jak je
vzorek starý?
13. Slon se dožije průměrně 65 let. Jedna samice má mládě jednou za 4 roky, poměr samců
a samic mezi novorozenými slůňaty je 1 : 1. Kolik zvířat je nutné každý rok odstřelit,
aby na území, jehož úživnost je několik tisíc slonů žilo trvale 250 slonic a 50 slonů?
V úlohách 14–25 najděte obecné řešení dané rovnice.
14. x(t + 2) − 16x(t) = 0
15. x(t + 2) + 16x(t) = 0
16. ∆3x(t) = 0
17. σ2
+ 2
2
x(t) = 0
18. (σ − 3)2 σ2
+ 4 x(t) = 0
19. x(t + 2) + 8x(t + 1) + 12x(t) = et
20. x(t + 2) − x(t + 1) − 6x(t) = 5 · 3t
21. x(t + 2) + x(t + 1) − 6x(t) = 5 · 3t
22. x(t + 2) + x(t + 1) − 12x(t) = t · 2t
23. x(t + 2) + 4x(t) = 8 · 2t cos 1
2πt
24. x(t + 2) − 5x(t + 1) + 6x(t) = t + 1
25. x(t + 2) − 5x(t + 1) + 4x(t) = 4t − t2
26. Najděte řešení počáteční úlohy x(t + 3) − 7x(t + 2) + 16x(t + 1) − 12x(t) = 0, x(0) = 0,
x(1) = 1, x(2) = 1.
2.4. CVIČENÍ 37
27. Napište lineární diferenční rovnici druhého řádu, jejímž řešením je posloupnost 1
8, 1
2, 1,
3, 7, 19, 47, . . . .
V úlohách 28–31 najděte řešení daného počátečního problému.
28.
x(t + 1) = −x(t) + y(t),
y(t + 1) = 2y(t),
x(0) = 1, y(0) = 2.
29. x(t + 1) =


1 −2 −2
0 0 −1
0 2 3

 x(t), x(0) =


1
1
0

 .
30. x(t + 1) =
2 1
0 2
x(t) +
t
1
, x(0) =
1
0
.
31. x(t + 1) =


4 1 2
0 2 −4
0 1 6

 x(t), x(0) =


ξ0
η0
ζ0

 .
Výsledky:
1. t!
2. 3
t(t−1)
2
3. et(t−1)
4. 4 −
3
2t
5. 1 +
et
− 1
e − 1
6. x(t) =
c, t = 0
2(t + 1), t ≥ 1
7. x ≡ 0 v P−∞, x(t) =
c
t
v P1
8. Uložená částka x(t) v t-tém období je řešením počáteční úlohy x(t + 1) = 1,008x(t) + 200,
x(0) = 0, tedy x(n) = 25 000 1,008n
− 1 .
9. Dluh x(t) v t-ém měsíci je řešením počáteční úlohy x(t + 1) = 12
√
1,12 x(t) − 380, x(0) = 12 000.
Závěrečná splátka x(37) = 268,99.
10. Nesplacený úvěr x(t) v t-ém měsíci je řešením počáteční úlohy x(t + 1) = 12
√
1,1 x(t) − b, x(0) =
80 000. Měsíční splátka b = 3 008,9
.
= 3 010.
11. Nesplacená hypotéka x(t) v t-ém měsíci je řešením počáteční úlohy x(t+1) = 12
√
1,08 x(t)−1 000,
x(360) = 0. Hypotéka x(0) = 136 283,5.
12. 5 700
ln 2 ln 100
70
.
= 2933
13. 1 425
52
.
= 27,4 samic a 1 585
52
.
= 30, 5 samců, tj. něco mezi 27 a 28 slonicemi a 30 a 31 slony každý
rok.
14. 4t
A + (−1)t
B
15. 4t
A cos 1
2 πt + B sin 1
2 πt
16. A + Bt + Ct2
17.
√
2t (A + Bt) cos 1
2 πt + (C + Dt) sin 1
2 πt
18. 3t
(A + Bt) + 2t
C cos 1
2 πt + D sin 1
2 πt
19. (−2)t
(A + 3t
B) +
et
(e + 2)(e + 6)
20. 3t−1
(A + t) + (−2)t
B
21. 2t
A + (−3)t
B + 5
2 3t−1
22. 3t
A + (−4)t
B − 1
6 t + 5
3 2t
23. 2t
A sin 1
2 πt + (B − t) cos 1
2 πt
24. 3t
A + 2t
B + 1
2 t + 5
4
25. A + B + 1
12 t 4t
+ 1
9 t3
− 1
18 t2
+ 7
54 t
38 KAPITOLA 2. LINEÁRNÍ ROVNICE
26. (3 + 2t)2t
− 3t+1
27. x(t + 2) − x(t + 1) − 4x(t) = 0
28. x(t) = 1
3 2t+1
− (−1)t
, y(t) = 2t+1
29. x(t) =


3 − 2t+1
2 − 2t
−2 + 2t+1


30. x(t) =
(t + 2)2t−1
− t
2t
− 1
31. x(t) =


1 1 1
0 −2 −1
0 1 1




4 0 0
0 4 1
0 0 4


t 

1 0 −1
0 −1 −1
0 1 2




ξ0
η0
ζ0

 = 4t


ξ0
η0
ζ0

+(η0 +2ζ0) t 4t−1


1
−2
1


Kapitola 3
Další explicitně řešitelné rovnice
V prvních třech částech této kapitoly uvedeme některé typy nelineárních diferenčních rovnic,
u nichž je známa substituce převádějící je na rovnice lineární. V poslední části ukážeme
speciální rovnice, které byly získány volbou goniometrických nebo hyperbolických funkcí na
místě transformující funkce. Řešení těchto rovnic bude užitečné pro nalezení explicitního řešení
logistické rovnice
x(t + 1) = rx(t) 1 − x(t) . (3.1)
pro několik speciálních hodnot parametru r.
3.1 Riccatiho a Bernoulliova rovnice
Riccatiho diferenční rovnice je tvaru
p(t)x(t + 1)x(t) + x(t + 1) − 1 + q(t) x(t) = r(t), (3.2)
kde p je nenulová regresivní posloupnost. Rovnici můžeme přepsat ve tvaru rekurentní formule
x(t + 1) =
1 + q(t) x(t) + r(t)
1 + p(t)x(t)
(3.3)
nebo explicitní diferenční rovnice prvního typu
∆x =
−p(t)x2 + q(t)x + r(t)
1 + p(t)x
.
S využitím operátorů posunu a diference můžeme rovnici (3.2) přepsat do tvaru
pxxσ
+ x + ∆x − (1 + q)x = r
a z něho vyjádřit diferenci hledané posloupnosti
∆x = −pxxσ
+ qx + r.
Tato rovnice je diskrétní analogií Riccatiho diferenciální rovnice x′ = −px2 + qx + r.
Riccatiho diferenční rovnici řešíme pomocí substituce
x(t) =
1
p(t)
∆y(t)
y(t)
=
y(t + 1) − y(t)
p(t)y(t)
. (3.4)
39
40 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Dosadíme do rekurentní formule (3.3) a postupně upravujeme:
y(t + 2) − y(t + 1)
p(t + 1)y(t + 1)
=
1 + q(t)
y(t + 1) − y(t)
p(t)y(t)
+ r(t)
1 +
y(t + 1) − y(t)
y(t)
y(t + 2) − y(t + 1)
p(t + 1)y(t + 1)
=
1 + q(t) y(t + 1) − y(t) + r(t)p(t)y(t)
p(t)y(t + 1)
y(t + 2) − y(t + 1) =
p(t + 1)
p(t)
1 + q(t) y(t + 1) −
p(t + 1)
p(t)
(1 + q(t) − r(t)p(t)) y(t).
Odtud vidíme, že posloupnost y je řešením lineární homogenní diferenční rovnice druhého
řádu
y(t + 2) − 1 +
p(t + 1)
p(t)
1 + q(t) y(t + 1) +
p(t + 1)
p(t)
(1 + q(t)) − r(t)p(t + 1) y(t) = 0,
kterou můžeme také zapsat stručněji pomocí operátorů posunu a diference
∆2
y + 1 −
pσ
p
(1 + q) ∆y − pσ
ry = 0.
Tvrzení 6. Riccatiho diferenční rovnice (3.2) pro neznámou posloupnost x se substitucí (3.4)
transformuje na lineární homogenní rovnici druhého řádu pro neznámou posloupnost y.
Příklad:
x(t + 1) =
2x(t) + 3
3x(t) + 2
, x(0) = x0
Zavedeme substituci
x(t) =
∆y(t)
3
2 y(t)
=
2
3
y(t + 1)
y(t)
−
2
3
,
dosadíme do dané rovnice a postupně ji upravíme
2
3
y(t + 2)
y(t + 1)
−
2
3
=
4
3
y(t + 1)
y(t)
− 4
3 + 3
2
y(t + 1)
y(t)
− 2 + 2
,
4
y(t + 2) − y(t + 1)
y(t + 1)
=
4y(t + 1) + 5y(t)
y(t + 1)
.
Daná rovnice se tedy transformuje na lineární homogenní rovnici druhého řádu
4y(t + 2) − 8y(t + 1) − 5y(t) = 0.
Její charakteristická rovnice 4λ2 − 8λ − 5 = 0 má dva reálné různé kořeny
λ1,2 =
8 ±
√
64 + 80
8
=
5
2
−1
2
.
3.1. RICCATIHO A BERNOULLIOVA ROVNICE 41
Obecné řešení lineární diferenční rovnice tedy je y(t) = A 5
2
t
+ B −1
2
t
.
Označme y0 = y(0). Pro počáteční hodnoty dále platí x0 = 2
3
y(1)
y(0)
− 2
3 , a z toho vypočítáme
y(1) = 1
2(3x0 + 2)y0.
Z těchto podmínek dostaneme systém (algebraických) rovnic pro konstanty A, B,
y0 = y(0) = A + B,
3x0 + 2
2
y0 = y(1) =
5
2
A −
1
2
B,
tj.
A + B = y0
5A − B = (3x0 + 2)y0.
Z něho vypočítáme A = 1
2 y0(1 + x0), B = 1
2 y0(1 − x0). Řešení úlohy pro lineární rovnici je
y(t) =
y0
2t+1
(1 + x0)5t
+ (1 − x0)(−1)t
.
Zpětnou substitucí tedy dostaneme řešení zadané úlohy ve tvaru
x(t) =
2
3
y(t + 1)
y(t)
− 1 =
2
3
1
2
(1 + x0)5t+1 + (1 − x0)(−1)t+1
(1 + x0)5t + (1 − x0)(−1)t
− 1 =
=
(1 + x0)5t − (1 − x0)(−1)t
(1 + x0)5t + (1 − x0)(−1)t
=
1 + x0
1 + x0 + (1 − x0) −1
5
t −
1 − x0
1 − x0 + (1 + x0)(−5)t
.
Pokud r ≡ 0, tj. na pravé straně rovnice (3.2) je nula, můžeme použít jednodušší substituci.
V tomto případě položíme
x(t) =
1
z(t)
, (3.5)
dosadíme do rovnice (3.2) a vynásobíme výrazem z(t)z(t + 1). Dostaneme
p(t) + z(t) − 1 + q(t) z(t + 1) = 0.
Je-li přitom posloupnost q regresivní, upravíme tuto rovnici na tvar lineární diferenční rovnice
prvního řádu
z(t + 1) =
1
1 + q(t)
z(t) +
p(t)
1 + q(t)
nebo ∆z = −
q(t)
1 + q(t)
z +
p(t)
1 + q(t)
.
Tato rovnice má podle Tvrzení 1 řešení
z(t) = z0
t−1
i=t0
1
1 + q(i)
+
t−1
i=t0
p(i)
1 + q(i)
t−1
j=i+1
1
1 + q(j)
=
=

z0 +
t−1
i=t0
p(i)
i−1
j=t0
1 + q(j)


t−1
i=t0
1
1 + q(i)
,
kde z0 = z(t0) = x(t0)−1. Platí tedy:
42 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Tvrzení 7. Je-li r ≡ 0, pak má Riccatiho rovnice řešení
x(t) =
x0
t−1
i=t0
1 + q(i)
1 + x0
t−1
i=t0
p(i)
i−1
j=t0
1 + q(j)
.
Rovnice (3.2) s r ≡ 0 a s regresivní posloupností q se v literatuře objevuje v rozmanitých
tvarem. Ukážeme některé z nich. Rovnici v takovém případě můžeme přepsat na tvar
p(t)
1 + q(t)
+
1
1 + q(t)
1
x(t)
−
1
x(t + 1)
= 0
a při označení
a(t) =
1
1 + q(t)
, b(t) =
p(t)
1 + q(t)
dostaneme
1
x(t + 1)
−
1
x(t)
= a(t) − 1
1
x(t)
+ b(t),
neboli
∆
1
x
= a(t) − 1
1
x
+ b(t) (3.6)
případně
∆x1−2
= a(t) − 1 x1−2
+ b(t). (3.7)
S pomocí operátoru posunu můžeme rovnici (3.6) přepsat ve tvaru
1
xσ
−
1
x
= (a − 1)
1
x
+ b.
Vynásobením výrazem xxσ dostaneme rovnici ve tvaru
x − xσ
= (a − 1)xσ
+ bxxσ
.
Z ní můžeme vyjádřit
∆x = (1 − a − bx)xσ
= (1 − a) 1 −
b
1 − a
x xσ
nebo
xσ
=
x
a + bx
. (3.8)
Poslední rovnici vynásobíme jmenovatelem pravé strany a upravíme na tvar
xσ
=
x
a
(1 − bxσ
) ,
ze kterého dostaneme jiné vyjádření diference hledané posloupnosti
∆x = x
1
a
− 1 −
b
a
xσ
=
1 − a
a
x 1 −
b
1 − a
xσ
.
3.1. RICCATIHO A BERNOULLIOVA ROVNICE 43
Bernoulliova diferenční rovnice je tvaru
∆x1−α
= a(t) − 1 x1−α
+ b(t), (3.9)
kde α = 1 je nějaké reálné číslo. Bernoulliovu diferenční rovnici můžeme také vyjádřit ve
tvaru rekurentní formule
x(t + 1) = a(t)x(t)1−α
+ b(t)
1/(1−α)
.
Porovnáním s rovnicí (3.7) vidíme, že Riccatiho rovnice (3.2) s r ≡ 0 je speciálním případem
Bernoulliovy rovnice (3.9) s parametrem α = 2.
Tvar Bernoulliovy rovnice bezprostředně ukazuje, že substituce
x(t)1−α
= z(t), tj. x(t) = z(t)1/(1−α)
(3.10)
transformuje Bernoulliovu diferenční rovnici (3.9) na lineární nehomogenní rovnici prvního
řádu
∆z = a(t) − 1 z + b(t), tj. z(t + 1) = a(t)z(t) + b(t). (3.11)
Tvrzení 8. Bernoulliova diferenční rovnice (3.9) pro neznámou posloupnost x se substitucí
(3.10) transformuje na lineární nehomogenní rovnici prvního řádu (3.11) pro neznámou posloupnost
z.
Příklad: Bevertonovu-Holtovu rovnici
x(t + 1) = x(t)
rK
K + (r − 1)x(t)
(3.12)
modelující vývoj velikosti populace v prostředí s omezenými zdroji můžeme přepsat ve tvaru
x(t + 1) = x(t)
r
1 + r−1
K x(t)
.
Jedná se o rovnici (3.8), tj. rovnici, která je současně Riccatiho i Bernoulliova. Můžeme ji tedy
vyřešit substitucí (3.5). Tuto substituci nyní odvodíme intuitivně, z úvahy o modelovaném
ději.
Budeme hledat řešení nenulové, tj. chceme modelovat nevyhynulou populaci. V tom případě
můžeme napsat rovnost převrácených hodnot obou stran rovnice (3.12)
1
x(t + 1)
=
1
r
1
x(t)
+
r − 1
rK
.
Nyní pro zjednodušení označíme z posloupnost převrácených hodnot posloupnosti x. Podle
předchozí rovnosti vidíme, že posloupnost z splňuje rekurentní formuli
z(t + 1) =
1
r
z(t) +
r − 1
rK
.
To je lineární rekurentní formule prvního řádu s konstantními koeficienty. Proto podle 2. důsledku
Tvrzení 1 můžeme obecný člen posloupnosti z vyjádřit ve tvaru
z(t) = z(t0)
1
r
t−t0
+
r − 1
rK
r−(t−t0) − 1
1 − r
r =
Kz(t0) + rt−t0 − 1
Krt−t0
,
44 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
přitom z(t0) = 1/x(t0). Můžeme tedy napsat obecný člen řešení Bevertonovy-Holtovy rovnice
s počáteční hodnotou x(t0) = x0 jako převrácenou hodnotu posloupnosti z, tj.
x(t) =
Krt−t0 x0
K + (rt−t0 − 1) x0
=
Kx0
x0 + (K − x0)r−(t−t0)
. (3.13)
Snadno ověříme, že touto formulí je skutečně zadáno řešení počáteční úlohy pro BevertonovuHoltovu
rovnici s počáteční hodnotou x(t0) = x0. Navíc takto zadaná posloupnost je v případě
x0 = 0 konstantní nulová, x ≡ 0; vyjadřuje tedy také řešení úlohy s počáteční hodnotou
x(t0) = 0.
Nyní můžeme snadno vyšetřit kvalitativní vlastnosti řešení Bevertonovy-Holtovy rovnice:
• Pro r = 1 nebo x0 = 0 je řešení x ≡ x0.
• Pokud r > 1, pak lim
t→∞
r−(t−t0) = 0, takže pro každou počáteční podmínku x0 = 0 řešení
x splňuje
lim
t→∞
x(t) = lim
t→∞
Kx0
x0 + (K − x0)r−(t−t0)
= K.
• Pokud r ∈ (0, 1), pak lim
t→∞
r−(t−t0) = ∞, takže pro každou počáteční podmínku x0 ∈ R
platí
lim
t→∞
x(t) = lim
t→∞
Kx0
x0 + (K − x0)r−(t−t0)
= 0.
Tyto výsledky dobře odpovídají ekologické intuici: pokud je vnitřní koeficient růstu r větší
než 1, tj. pokud v neomezeném prostředí má populace porodnost větší než úmrtnost, pak se
její velikost ustálí na kapacitě prostředí; pokud je úmrtnost větší než porodnost, populace
vymře.
3.2 Homogenní rovnice
Homogenní diferenční rovnice prvního řádu je rovnice tvaru
f t,
x(t + 1)
x(t)
= 0, (3.14)
kde f je funkce, která není konstantní ve druhé proměnné. Povšimněme si, že lineární homogenní
rovnici x(t + 1) = q(t)x(t) můžeme přepsat jako
x(t + 1)
x(t)
− q(t) = 0,
takže je skutečně speciálním případem rovnice (3.14); slovo „homogenní je použito opráv-
něně.
Substituce
y(t) =
x(t + 1)
x(t)
(3.15)
převede rovnici (3.14) na rovnici
f t, y(t) = 0,
3.2. HOMOGENNÍ ROVNICE 45
ze které vyjádříme y(t) = g(t) a řešení dané rovnice (3.14) hledáme jako řešení lineární
homogenní rovnice x(t + 1) = g(t)x(t).
Pokud hledáme kladná řešení rovnice (3.14), můžeme použít substituce
z(t) = ln x(t),
která převádí danou rovnici na implicitní diferenční rovnici
f t, ∆z(t) = 0.
Homogenní diferenční rovnice k-tého řádu je rovnice tvaru
F t,
x(t + k)
x(t + k − 1)
,
x(t + k − 1)
x(t + k − 2)
, . . . ,
x(t + 1)
x(t)
= 0,
kde F je funkce, která není konstantní ve druhé a poslední proměnné. Tuto rovnici převede
substituce (3.15) na diferenční rovnici (k − 1)-ního řádu druhého typu
F t, y(t + k − 1), y(t + k − 2), . . . , y(t) = 0.
Příklad: Najdeme řešení homogenní rovnice druhého řádu
x(t + 1) =
x(t − 1)x(t)
x(t − 1) + x(t)
.
Rovnici vynásobíme jmenovatelem zlomku na pravé straně a vydělíme výrazem x(t)x(t − 1).
Dostaneme
1 +
x(t)
x(t − 1)
x(t + 1)
x(t)
= 1.
Substituce (3.15) převede tuto rovnici na tvar
1 + y(t − 1) y(t) = 1,
který je ekvivalentní s 1 + y(t) y(t + 1) = 1, neboli
y(t + 1)y(t) + y(t + 1) = 1.
To je Riccatiho rovnice. Proto zavedeme novou posloupnost z substitucí
y(t) =
z(t + 1) − z(t)
z(t)
.
Po dosazení a úpravě dostaneme
z(t + 2) − z(t + 1)
z(t + 1)
z(t + 1) − z(t)
z(t)
+ 1 = 1,
z(t + 2) − z(t + 1) − z(t) = 0,
což je lineární homogenní rovnice druhého řádu.
46 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
3.2.1 Implicitní rovnice x(t + 1)2
+ a(t)x(t + 1)x(t) + b(t)x(t)2
= 0
Tato rovnice má očividně řešení x ≡ 0. Budeme hledat také řešení nenulová. Rovnici vydělíme
výrazem x(t)2 a tím ji převedeme na tvar rovnice homogenní
x(t + 1)
x(t)
2
+ a(t)
x(t + 1)
x(t)
+ b(t) = 0.
Pokud posloupnosti a, b splňují relaci a(t)2 ≥ 4b(t) pro všechny indexy, položíme
p(t) =
1
2
−a(t) + a(t)2 − 4b(t) a q(t) =
1
2
−a(t) − a(t)2 − 4b(t) .
a pravou stranu rovnice přepíšeme jako součin dvou výrazů
x(t + 1)
x(t)
− p(t)
x(t + 1)
x(t)
− q(t) = 0.
Odtud vidíme, že řešení každé z lineárních homogenních diferenčních rovnic prvního řádu
x1(t + 1) = p(t)x1(t) a x2(t + 1) = q(t)x2(t)
je také řešením původní rovnice. Tato řešení jsou
x(t) = x0
t−1
i=t0
p(i) a x(t) = x0
t−1
i=t0
q(i).
Povšimněme si, že nulové řešení je v tomto vyjádření zahrnuto pro x0 = 0. Pokud a2 = 4b,
pak p = q = −1
2a a obě řešení splývají, x(t) = x0(−1
2)t−t0
t−1
i=t0
a(i).
Tvrzení 9. Počáteční úloha pro implicitní rovnici tvaru
x(t + 1)2
+ a(t)x(t + 1)x(t) + b(t)x(t)2
= 0, x(t0) = x0
má pro x0 = 0 a a2 = 4b dvě řešení. V opačném případě je jednoznačně řešitelná.
Příklad: Najdeme všechna řešení rovnice v implicitním tvaru
x(t + 1)2
− 3x(t)x(t + 1) + 2x(t)2
= 0.
Rovnici postupně upravíme:
x(t + 1)
x(t)
2
− 3
x(t + 1)
x(t)
+ 2 = 0,
x(t + 1)
x(t)
− 2
x(t + 1)
x(t)
− 1 = 0.
Dostáváme tak dvě homogenní lineární rovnice x(t + 1) = 2x(t) a x(t + 1) = x(t). Daná
rovnice má tedy dvě řešení, konkrétně
x(t) = x02t−t0
a x(t) = x0,
kde x0 = x(t0). Tato řešení splývají pro x0 = 0.
3.3. LOGARITMICKY LINEÁRNÍ ROVNICE 47
3.3 Logaritmicky lineární rovnice
Jedná se o rovnici
x(t + k)rk(t)
x(t + k − 1)rk−1(t)
· · · x(t + 1)r1(t)
x(t)r0(t)
= b(t);
přitom r0, r1, . . . , rk jsou posloupnosti takové, že r0(t) = 0 = rk(t) pro všechna t z definičního
oboru. Substitucí
x(t) = ez(t)
(3.16)
tj. z(t) = ln x(t) převedeme uvažovanou rovnici na tvar
erk(t)z(t+k)+rk−1(t)z(t+k−1)+···+r1(t)z(t+1)+r0(t)z(t)
= b(t),
a dále zlogaritmováním na lineární rovnici k-tého řádu
z(t + k) +
rk−1(t)
rk(t)
z(t + k − 1) + · · · +
r1(t)
rk(t)
z(t + 1) +
r0(t)
rk(t)
z(t) =
ln b(t)
rk(t)
.
Povšimněme si, že z transformačního vztahu (3.16) plyne, že řešení původní rovnice musí být
kladné. Uvedený postup tedy můžeme použít pouze v případě, že počáteční hodnoty hledané
posloupnosti splňují podmínky
x(t0) = x0 > 0, x(t0 + 1) = x1 > 0. . . . , x(t0 + k − 1) = xk−1 > 0.
Příklad:
x(t + 2) =
x(t + 1)
x(t)
2
.
Rovnici přepíšeme ve tvaru
x(t + 2)x(t + 1)−2
x(t)2
= 1
a zavedeme substituci z(t) = ln x(t). Dostaneme lineární homogenní rovnici druhého řádu
z(t + 2) − 2z(t + 1) + 2z(t) = 0.
Její charakteristická rovnice λ2 − 2λ + 2 = má komplexně sdružené kořeny
λ1,2 =
2 ±
√
4 − 8
2
= 1 ± i.
Modul a argument charakteristických kořenů jsou
|λ| =
√
1 + 1 =
√
2, arg λ = arctg 1 =
1
4
π.
To znamená, že obecné řešení lineární rovnice je
z(t) =
√
2 t
A cos
πt
4
+ B sin
πt
4
a obecné řešení dané rovnice je
x(t) = exp
√
2 t
A cos
πt
4
+ B sin
πt
4
.
48 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
3.4 Rovnice řešitelné speciálními substitucemi
3.4.1 Goniometrické a hyperbolické substituce
Ukážeme několik speciálních rovnic, u kterých lze najít explicitní řešení pomocí goniometrické
nebo hyperbolické substituce. U všech těchto rovnic budeme uvažovat také počáteční
podmínku
x(t0) = x0. (3.17)
Uvedené rovnice byly získány pomocí známých vztahů pro goniometrické nebo hyperbolické
funkce násobného argumentu. Je z nich zřejmé, jak lze odvozovat další explicitně
řešitelné rovnice. Navíc téměř libovolnou transformací hledané posloupnosti lze z uvedených
rovnic získat další rovnice, které jsou opět explicitně řešitelné. Tuto skutečnost ukážeme na
příkladech
Rovnice x(t + 1) = 2x(t)2 − 1
Řešení uvažované úlohy je pro libovolné t ≥ t0 určeno jednoznačně, neboť se jedná o rekurentní
formuli prvního řádu s počáteční podmínkou. Přitom je na pravé straně rovnosti výraz
definovaný pro jakoukoliv hodnotu x(t).
Z rovnice a z počáteční podmínky (3.17) plyne, že hodnota řešení x(t0 − 1) musí splňovat
rovnici x0 = 2 [x(t0 − 1)]2
− 1. Tato algebraická rovnice pro neznámou x(t0 − 1) nemá reálné
řešení, pokud x0 < −1, a má dvě různá reálná řešení, pokud x0 > −1. Obecně tedy úloha
není jednoznačně řešitelná pro t < t0. Proto budeme řešení hledat pouze pro t ≥ t0.
Pokud počáteční hodnota splňuje nerovnost |x0| ≤ 1, položíme x(t) = cos y(t). S využitím
známých vztahů pro goniometrické funkce
(cos α)2
+ (sin α)2
= 1 a cos 2α = (cos α)2
− (sin α)2
dostaneme
cos y(t + 1) = x(t + 1) = 2x(t)2
− 1 = 2 cos y(t)
2
− 1 = cos y(t)
2
− sin y(t)
2
= cos 2y(t).
To znamená, že y(t + 1) je řešením goniometrické rovnice
cos y(t + 1) = cos 2y(t),
a tedy
y(t + 1) = ±2y(t) + 2kπ, k ∈ Z.
Každá z tohoto spočetného systému lineárních nehomogenních diferenčních rovnic prvního
řádu má podle 2. důsledku Tvrzení 1 řešení tvaru
y(t) = y0(±2)t−t0
+ 2kπ
(±2)t−t0 − 1
±2 − 1
, k ∈ Z,
kde y0 = y(t0), tj. cos y0 = x0, y0 = arccos x0. Druhý sčítanec na pravé straně rovnosti
můžeme upravit na tvar
2kπ
(±2)t−t0 − 1
±2 − 1
= 2kπ
t−t0−1
i=0
(±2)i
.
3.4. ROVNICE ŘEŠITELNÉ SPECIÁLNÍMI SUBSTITUCEMI 49
Součet celých čísel je celé číslo a to znamená, že druhý sčítanec je sudým násobkem π, tj.
y(t) = y0(±2)t−t0
+ 2lπ
pro nějaké l ∈ Z. Zpětnou substitucí a úpravou s využitím součtového vzorce
cos(α + β) = cos α cos β − sin α sin β
dostaneme
x(t) = cos y0(±2)t−t0
+ 2lπ = cos y0(±2)t−t0
cos 2lπ − sin y0(±2)t−t0
sin 2lπ =
= cos (±1)t−t0
2t−t0
y0 = cos 2t−t0
y0 ,
neboť cosinus je sudá funkce. Řešení úlohy
x(t + 1) = 2x(t)2
− 1, x(t0) = x0 ∈ [−1, 1] (3.18)
je tedy tvaru
x(t) = cos 2t−t0
arccos x0 .
Pokud je |x0| > 1, položíme x(t) = cosh y(t). S využitím známého vztahu pro hyperbolický
cosinus1
cosh 2α = 2(cosh α)2
− 1
dostaneme
cosh y(t + 1) = x(t + 1) = 2x(t)2
− 1 = 2 cosh y(t)
2
− 1 = cosh 2y(t).
Hodnota y(t + 1) je tedy řešením rovnice cosh y(t + 1) = cosh 2y(t). Poněvadž hyperbolický
cosinus je sudá funkce, která je ryze monotonní na každém z intervalů (−∞, 0] a [0, ∞), platí
y(t + 1) = ±2y(t),
což jsou dvě rekurentní formule pro geometrickou posloupnost, jedna má kvocient 2, druhá
−2. Posloupnost tedy můžeme vyjádřit ve tvaru
y(t) = y0(±2)t−t0
= (±1)t−t0
2t−t0
y0,
kde y0 = y(t0), tj. cosh y0 = x0,
y0 = argcosh x0 = ln |x0| + x2
0 − 1 .
Poněvadž hyperbolický cosinus je sudá funkce, dostaneme řešení úlohy s počáteční hodnotou
x0 ∈ (−∞, −1) ∪ (1, ∞) ve tvaru
x(t) =
x0, t = t0,
cosh 2t−t0 ln |x0| + x2
0 − 1 , t > t0.
1
cosh 2α = 1
2
e2α
+ e−2α
= 1
2
eα
+ e−α 2
− 2 = 2 1
2
(eα
+ e−α
)
2
− 1 = 2(cosh α)2
− 1
50 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Příklad:
x(t + 1) = 2x(t) 2x(t) − 1 , x(0) = 1
8
Nejprve upravíme pravou stranu rovnice tak, aby byla tvaru f(2X2 − 1) pro nějakou funkci
f a nějaký výraz X závisející na x(t); použijeme doplnění na úplný čtverec:
4x(t)2
− 2x(t) = 2x(t) − 1
2
2
− 1
4 = 1
2 2 2x(t) − 1
2
2
− 1 + 1
4 .
Odtud vidíme, že daná diferenční rovnice je ekvivalentní s rovnicí
2x(t + 1) − 1
2 = 2 2x(t) − 1
2
2
− 1.
Můžeme tedy použít substituci y(t) = 2x(t) − 1
2, která převádí danou úlohu na počáteční
úlohu ve tvaru
y(t + 1) = 2y(t)2
− 1, y(0) = −1
4,
která má řešení
y(t) = cos 2t
arccos −1
4 = cos 2t
π − arccos 1
4 =
−1
4, t = 0
cos 2t arccos 1
4 , t > 0.
Řešení dané úlohy je tedy pro t > 0 dáno výrazem
x(t) = 1
2 cos 2t
arccos 1
4 + 1
2 = 1
4 + 1
2 cos 2t
arccos 1
4 .
Rovnice x(t + 1) = 2x(t) 1 ± x(t)2
Řešení uvažované úlohy je pro t ≥ t0 určeno jednoznačně, neboť se jedná o rekurentní formuli
prvního řádu s počáteční podmínkou a odmocninu považujeme v reálném oboru za jednoznačnou
funkci. Řešení ovšem v případě znaménka „− pod odmocninou nemusí být definováno
pro každé t ≥ t0; je-li totiž v takovém případě |x(t)| > 1, pak není x(t + 1) definováno.
Z rovnice a z počáteční podmínky plyne, že pro hodnotu x(t0 − 1) řešení by mělo platit
x0 = 2x(t0 − 1) 1 ± x(t0 − 1)2,
nebo po snadné úpravě
±4 [x(t0 − 1)]4
+ 4 [x(t0 − 1)]2
− x2
0 = 0,
takže by mělo platit
x(t0 − 1) =
−4 ± 16 ∓ 16x2
0
±8
=
1
2
1 ± x2
0 ∓ 1 ;
hodnota x(t0 − 1) tedy není určena jednoznačně. Z tohoto důvodu má smysl uvažovat řešení
pouze pro t ≥ t0.
Pokud je pod odmocninou na pravé straně rovnice znaménko „+ , tedy pokud je rovnice
tvaru
x(t + 1) = 2x(t) 1 + x(t)2, (3.19)
3.4. ROVNICE ŘEŠITELNÉ SPECIÁLNÍMI SUBSTITUCEMI 51
zavedeme substituci x(t) = sinh y(t) a využijeme známých vlastností hyperbolických funkcí
sinh 2α = 2 sinh α cosh α, (cosh α)2
− (sinh α)2
= 1, cosh α > 0.
Pak je
sinh y(t + 1) = x(t + 1) = 2x(t) 1 + x(t)2 = 2 sinh y(t) 1 + sinh y(t)
2
=
= 2 sinh y(t) cosh y(t) = sinh 2y(t).
Poněvadž hyperbolický sinus je prostá funkce, implicitní diferenční rovnice
sinh y(t + 1) = sinh y(t)
je ekvivalentní s explicitní rovnicí y(t + 1) = 2y(t) a její řešení je tvaru
y(t) = 2t−t0
y0,
kde y0 = y(t0), tj. x0 = sinh y0, y0 = argsinh x0 = ln x0 + x2
0 + 1 .
Zpětnou substitucí dostaneme řešení úlohy (3.19), (3.17) pro t ≥ t0 ve tvaru
x(t) = sinh 2t−t0
argsinh x0 =
1
2
x0 + 1 + x2
0
2t−t0
− x0 + 1 + x2
0
−2t−t0
.
(3.20)
Rovnice se znaménkem „− pod odmocninou na pravé straně, tedy rovnice
x(t + 1) = 2x(t) 1 − x(t)2 (3.21)
může mít řešení pouze pro počáteční podmínku x0 ∈ [−1, 1], pro |x| > 0 není pravá strana
rovnice definována. Navíc pro všechny hodnoty řešení musí platit |x(t)| ≤ 1. Pokud je x0 = 0,
bude řešením úlohy (3.21), (3.17) konstantní posloupnost x ≡ 0.
Dále si můžeme všimnout, že pro řešení úlohy platí
x(t + 1)x(t) = 2x(t)2
1 − x(t)2 ≥ 0,
neboť odmocninu v reálném oboru chápeme jako nezápornou funkci. To znamená, že řešení
rovnice nemění znaménko, tj. sgn x(t) = sgn x0 pro všechna t ≥ t0.
Toto pozorování umožňuje zavést substituci
x(t) = | sin y(t)| sgn x0. (3.22)
Využijeme známých vlastností goniometrických funkcí
sin 2α = 2 sin α cos α, (cos α)2
+ (sin α)2
= 1
a pro x0 = 0 dostaneme
| sin y(t + 1)| = 2| sin y(t)| 1 − sin y(t)
2
= 2| sin y(t)| · | cos y(t)| =
= |2 sin y(t) cos y(t)| = | sin 2y(t)|.
52 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Řešíme tedy goniometrickou rovnici s absolutní hodnotou | sin y(t + 1)| = | sin y(t)| pro neznámou
y(t + 1). Řešení této rovnice může být řešením některé ze dvou goniometrických
rovnic
sin y(t + 1) = sin 2y(t), sin y(t + 1) = − sin 2y(t)
pro neznámou y(t + 1). První z těchto rovnic má dvě řešení
y(t + 1) = 2y(t) + 2kπ, y(t + 1) = −2y(t) + (2k + 1)π, k ∈ Z,
druhá má také dvě řešení
y(t + 1) = −2y(t) + 2kπ, y(t + 1) = 2y(t) + (2k + 1)π, k ∈ Z.
To znamená, že posloupnost y splňuje některou z nekonečného systému lineárních rekurentních
formulí prvního řádu
y(t + 1) = ±2y(t) + kπ, k ∈ Z,
nebo ekvivalentně lineárních diferenčních rovnic ∆y = (±2 − 1)y + kπ. Podle 2. důsledku
Tvrzení 1 je řešení těchto rovnic tvaru
y(t) = y0(±2)t−t0
+ kπ
(±2)t−t0 − 1
±2 − 1
= y0(±2)t−t0
+ kπ
t−t0−1
i=0
(±2)i
,
kde y0 = y(t0) = arcsin x0. Sumu na pravé straně rovnosti lze vyjádřit jako
t−t0−1
i=0
(±2)i
=
0, t = t0,
1 ± 2 + 4 ± · · · + (±2)t−t0−1, t > t0,
což znamená, že pro t > t0 je rovna lichému celému číslu. Proto pro t > t0 platí
y(t) = (±2)t−t0
y0 + (2l + 1)π, l ∈ Z
a tedy také
x(t) = | sin y(t)| sgn x0 =
= sin (±2)t−t0
y0 cos (2l + 1)π + cos (±2)t−t0
y0 sin (2l + 1)π sgn x0 =
= −(∓1)t−t0
sin 2t−t0
y0 + 0 sgn x0 = (±1)t−t0
sin 2t−t0
y0 sgn x0 =
= sin 2t−t0
arcsin x0 sgn x0.
Dostáváme tak výsledek, že řešení počáteční úlohy (3.21), (3.17) s x0 ∈ [−1, 1] je pro t ≥ t0
dáno formulí
x(t) = sin 2t−t0
arcsin x0 sgn x0; (3.23)
toto řešení nemění znaménko a pro libovolné t ≥ t0 splňuje nerovnost |x(t)| ≤ 1.
Příklad:
x(t + 1) =
x(t)2
4 x(t) + 1
, x(0) = −3
3.4. ROVNICE ŘEŠITELNÉ SPECIÁLNÍMI SUBSTITUCEMI 53
Zavedeme substituci x(t) = −
1
y(t)2
. Pak
1
y(t + 1)2
= −x(t + 1) = −
x(t)2
4 x(t) + 1
=
−1
4y(t)4 −
1
y(t)2
+ 1
=
1
4y(t)2 1 − y(t)2
.
Tedy y(t + 1)2 = 4y(t)2 1 − y(t)2 , neboli
y(t + 1) = 2y(t) 1 − y(t)2.
Řešení této rovnice s počáteční podmínkou y(0) =
√
3
3 je podle předchozího výsledku dáno
výrazem y(t) = sin 2t arcsin
√
3
3 . Řešení dané úlohy tedy je
x(t) = −

 1
sin 2t arcsin
√
3
3


2
.
Rovnice x(t + 1) =
2x(t)
1 − x(t)2
, x(t + 1) =
x(t)2 − 1
2x(t)
Opět má smysl řešit počáteční úlohu pouze pro t ≥ t0.
V případě první rovnice zavedeme substituci x(t) = tg y(t) a využijeme vzorec pro tangens
dvojnásobného argumentu
tg 2α =
2 tg α
1 − (tg α)2
.
Pak je
tg y(t + 1) = x(t + 1) =
2x(t)
1 − x(t)2
=
2 tg y(t)
1 − tg y(t)
2 = tg 2y(t).
Řešíme tedy goniometrickou rovnici tg y(t + 1) = tg 2y(t) pro neznámou y(t + 1). Dostaneme
y(t + 1) = 2y(t) + kπ, k ∈ Z.
Tato lineární nehomogenní rovnice prvního řádu má řešení
y(t) = y02t−t0
+ kπ 2t−t0
− 1 ,
kde y0 = y(0) = arctg x0. Zpětnou substitucí dostaneme
x(t) = tg 2t−t0
arctg x0 + (2t−t0
− 1)kπ =
=
tg 2t−t0 arctg x0 + tg (2t−t0 − 1)kπ
1 − tg 2t−t0 arctg x0 tg (2t−t0 − 1)kπ
= tg 2t−t0
arctg x0 .
Řešení počáteční úlohy
x(t + 1) =
2x(t)
1 − x(t)2
, x(t0) = x0
54 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
je dáno výrazem
x(t) = tg 2t−t0
arctg x0 .
Druhou rovnici řešíme analogicky, použijeme substituci x(t) = cotg y(t).
Příklad:
x(t + 1) = 2
x(t) − 1
x(t) 2 − x(t)
+ 1, x(0) = 1
2.
Rovnici postupně upravujeme
1 − x(t + 1) = 2
1 − x(t)
x(t) 2 − x(t)
,
1 − x(t + 1) = 2
1 − x(t)
1 − 1 − x(t) 1 + 1 − x(t)
,
1 − x(t + 1) = 2
1 − x(t)
1 − 1 − x(t)
2 .
Tento zápis rovnice ukazuje, že substituce y(t) = 1 − x(t) rovnici transformuje na uvažovaný
tvar. Řešení úlohy je tedy dáno relací
1 − x(t) = tg 2t
arctg(1 − x0) = tg 2t−1 π
3
,
neboť 1 − x0 = 1
2, arctg 1
2 = 1
6π. Řešení dané úlohy tedy je
x(t) = 1 − tg
2t−1
3
π .
3.4.2 Logistická rovnice
Logistická diferenční rovnice je rovnice tvaru
x(t + 1) = ax(t) 1 − x(t) , nebo ∆x = x(a − 1 − ax),
kde a ∈ R, a = 0.
Rekurentní vztah pro hledanou posloupnost x ukazuje, že logistická rovnice s počáteční
podmínkou x(t0) = x0 má jednoznačné řešení pro t ≥ t0. Z rekurentního vztahu však obecně
nelze jednoznačně vyjádřit hodnotu x(t) v závislosti na x(t + 1). Z rovnice
ax(t)2
− ax(t) + x(t + 1) = 0
totiž vychází
x(t) =
1
2
1 ± 1 −
4x(t + 1)
a
.
Odtud plyne, že logistická rovnice s počáteční podmínkou x(t0) = x0 má reálné řešení pro
nějaké indexy t < t0 pouze v případě, že x0 ≤ 1
4a. Toto řešení však obecně není vyjádřeno
jednoznačně. Jedinou výjimkou je případ a = 2, x0 = 1
2 ; pak je konstantní posloupnost x ≡ 1
2
3.4. ROVNICE ŘEŠITELNÉ SPECIÁLNÍMI SUBSTITUCEMI 55
jednoznačným řešením počáteční úlohy definovaným na celé množině Z. V případě x0 = 1
4 a,
a = 2 totiž dostaneme x(t0 − 1) = 1
2 = x0 a hodnota x(t0 − 2) již není určena jednoznačně.
Z vyjádření diference hledané posloupnosti x bezprostředně plyne, že konstantní posloup-
nosti
x ≡ 0 a x ≡ 1 −
1
a
(3.24)
jsou řešeními logistické rovnice definovanými na celé množině Z. Tyto posloupnosti ovšem
obecně nelze považovat za řešení logistické rovnice s počáteční hodnotou x0 = 0 nebo x0 =
1 − 1/a.
Počáteční úlohu pro logistickou rovnici vyřešíme pouze ve třech speciálních případech.
Případ a = 2: Budeme řešit počáteční úlohu
x(t + 1) = 2x(t) 1 − x(t) , x(t0) = x0, (3.25)
která je ekvivalentní s úlohou
∆x = x(1 − 2x), x(t0) = x0.
Nejprve zavedeme substituci
y(t) = 1 − 2x(t), tj. x(t) =
1
2
1 − y(t) . (3.26)
Po dosazení a úpravách postupně dostaneme
1
2
1 − y(t + 1) = 1 − y(t) 1 −
1
2
1 − y(t)
1
2
1 − y(t + 1) =
1
2
1 − y(t)2
y(t + 1) = y(t)2
y(t + 1)y(t)−2
= 1.
To je logaritmicky lineární rovnice. Jejím logaritmováním dostaneme
ln y(t + 1) − 2 ln y(t) = 0,
což je lineární homogenní rovnice z(t + 1) = 2z(t) pro posloupnost
z(t) = ln y(t). (3.27)
To znamená, že z(t) = z02t−t0 , kde z0 = ln y(t0) = ln 1 − 2x(t0) . Odtud dostaneme
y(t) = ez(t)
= exp 2t−t0
ln(1 − 2x0) = (1 − 2x0)2t−t0
.
Zpětnou substitucí dostaneme řešení úlohy ve tvaru
x(t) =
1
2
1 − (1 − 2x0)2t−t0
. (3.28)
56 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Při řešení úlohy jsme mlčky předpokládali, že výraz ln(1 − 2x0) je definován, tedy že x0 < 1
2 .
Výraz na pravé straně rovnosti (3.28) je však definován pro každé x0 ∈ R. Platí totiž
(1 − 2x0)2t−t0
= |1 − 2x0| sgn(1 − 2x0)
2t−t0
=
1 − 2x0, t = t0,
exp 2t−t0 ln |1 − 2x0| , t = t0.
Přímým výpočtem se můžeme přesvědčit, že rovností (3.28) je skutečně definováno řešení
úlohy (3.25) pro libovolnou počáteční hodnotu.
Postup hledání tvaru (3.28) řešení počáteční úlohy (3.25) pomocí substitucí (3.26) a (3.27)
můžeme popsat ve zhuštěné formě: substitucí
1 − 2x(t) = exp z(t)
najdeme řešení úlohy (3.25) ve tvaru
1 − 2x(t) = exp 2t−t0
ln(1 − 2x0) .
Kvalitativní vlastnosti řešení úlohy (3.25). Pro x0 ∈ (0, 1) je |1 − 2x0| < 1, takže řešení
úlohy s takovými počátečními hodnotami splňuje
lim
t→∞
x(t) = 1
2 − 1
2 lim
t→∞
(1 − 2x0)2t−t0
= 1
2.
Pro x0 ∈ {0, 1} a t > t0 je (1 − 2x0)2t−t0
= 1, neboť číslo 2t−t0 je sudé. Proto řešení úlohy
(3.25) s počáteční podmínkou x0 ∈ {0, 1} splňuje rovnost x(t) = 0 pro každé t > t0.
Pro x0 > 1 nebo x0 < 0 je |1 − 2x0| > 1 a proto lim
t→∞
x(t) = −∞.
Pro x0 = 1
2 je řešení úlohy (3.25) rovno x ≡ 1
2 v souladu s (3.24). Toto řešení je definováno
na celé množině Z, jak již bylo předesláno.
Případ a = 4: Nejprve budeme hledat nezáporné řešení počáteční úlohy
x(t + 1) = 4x(t) 1 − x(t) , x(t0) = x0. (3.29)
To znamená, že budeme předpokládat, že x0 ∈ [0, 1] a zavedeme substituci
x(t) = y(t)2
, tj. y(t) = x(t) . (3.30)
Po dosazení do rekurentní formule v (3.29) dostaneme
y(t + 1)2
= 4y(t)2
1 − y(t)2
= 2y(t) 1 − y(t)2
2
,
tedy
y(t + 1) = 2y(t) 1 − y(t)2 .
Počáteční podmínka bude y(t0) =
√
x0 ≥ 0. Jedná se tedy o rovnici tvaru (3.21) s nezápornou
počáteční hodnotou, kterou řešíme substitucí
y(t) = | sin z(t)|. (3.31)
3.4. ROVNICE ŘEŠITELNÉ SPECIÁLNÍMI SUBSTITUCEMI 57
Podle (3.23) dostaneme řešení ve tvaru
y(t) = sin 2t−t0
arcsin
√
x0 .
Zpětnou substitucí dostaneme řešení x(t) = y(t)2 úlohy (3.29) vyjádřené formulí
x(t) = sin 2t−t0
arcsin
√
x0
2
. (3.32)
Přímým výpočtem můžeme ověřit, že tato posloupnost je skutečně řešením úlohy (3.29).
Postup hledání řešení tvaru (3.32) počáteční úlohy (3.29) s x0 ∈ [0, 1] pomocí substitucí
(3.30) a (3.31) můžeme opět zformulovat ve zhuštěné podobě: Substitucí
x(t) = | sin z(t)|
najdeme řešení úlohy (3.29) s x0 ∈ [0, 1] ve tvaru
x(t) = sin 2t−t0
arcsin
√
x0 .
Kvalitativní vlastnosti řešení úlohy (3.29) s x0 ∈ [0, 1].
Nejprve ukážeme, že logistická rovnice s parametrem a = 4 má periodická řešení libovolné
periody. Buď tedy n libovolné přirozené číslo a položíme
x0 = sin
2n−1
2n + 1
π
2
.
Dosazením do (3.32) dostaneme
x(t0 + n) =

sin

2n
arcsin sin
2n−1
2n + 1
π
2




2
= sin 2n 2n−1
2n + 1
π
2
=
= sin (2n
+ 1 − 1)
2n−1
2n + 1
π
2
= sin 2n−1
π −
2n−1
2n + 1
π
2
= − sin
2n−1
2n + 1
π
2
= x0.
Pro n = 1 je
sin
2n−1
2n + 1
π
2
= sin 1
3π
2
=
√
3
2
2
= 3
4
v souladu s (3.24).
V případě rovnice (3.25) libovolné řešení s počáteční hodnotou x0 ∈ (0, 1) konvergovalo
ke konstantnímu nenulovému řešení. Rovnice (3.29) tuto vlastnost nemá, periodická řešení
samozřejmě nekonvergují. Ovšem existují taková řešení, která ke 3
4 konvergují. Uvažme řešení
s počáteční hodnotou
x0 = sin
1
3 · 2n
π
2
a buď k ∈ N0 libovolné. Pak platí
x(t0 + k) = sin 2k 1
3 · 2n
π
2
=



3
4, k ≥ n,
sin 1
3·2n−k π
2
= 3
4 , k < n.
58 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Toto řešení tedy konverguje ke 3
4 a to tak, že po n krocích této hodnoty dosáhne a zůstane
konstantní.
Počáteční úloha pro logistickou rovnici s parametrem a = 4 má také řešení, které je
nenulové pouze pro konečně mnoho indexů, tj. řešení, které „po konečně mnoha krocích
vymizí . Buď opět n ∈ N libovolné číslo a položme
x0 = sin
π
2n
2
.
Pro libovolné k ∈ N0 platí
x(t0 + k) = sin 2k π
2n
2
= sin 2k−n
π
2
.
To znamená, že x(t0 + k) = 0 pro k ≥ 0 a x(t0 + k) > 0 pro k ∈ {0, 1, . . . , n − 1}.
Nyní hledejme řešení úlohy (3.29) s počáteční hodnotou x0 splňující nerovnost x0 > 1
nebo x0 < 0, tj. x0 − 1
2 > 1
2 .
Nejprve si všimněme, že pro x0 > 1 je x(t0 + 1) = 4x0(1 − x0) < 0. Dále pro x(t) < 0 je
také
x(t + 1) = 4x(t) 1 − x(t) = −4|x(t)| 1 + |x(t)| < 0,
což znamená, že pokud v nějakém t1 je řešení úlohy (3.29) záporné, pak je záporné pro každé
t ≥ t1. Celkem tak dostáváme, že pro počáteční hodnotu x0 splňující nerovnost x0 − 1
2 > 1
2 ,
řešení úlohy (3.29) splňuje nerovnost x(t) < 0 pro všechna t ≥ t0 +1. Budeme tedy řešit úlohu
x(t + 1) = 4x(t) 1 − x(t) , x(t0 + 1) = x1 < 0. (3.33)
Poněvadž řešení této úlohy je záporné, můžeme použít substituci
x(t) = −y(t)2
, tj. |y(t)| = −x(t). (3.34)
Dosazením do rovnice v (3.33) dostaneme
y(t + 1)2
= −x(t + 1) = −4x(t) 1 − x(t) = 4y(t)2
1 + y(t)2
,
neboli
|y(t + 1)|2
= 2|y(t)| 1 + |y(t)|2.
To je diferenční rovnice tvaru (3.19) pro posloupnost |y|. Příslušná počáteční podmínka je
y(t0 + 1)| =
√
−x1. Tuto úlohu řešíme substitucí |y(t)| = sinh z(t) a podle (3.20) dostaneme
její řešení ve tvaru
|y(t)| = sinh 2t−t0−1
argsinh |y(t0 + 1)| = sinh 2t−t0−1
argsinh
√
−x1 .
Zpětnou substitucí (3.34) napíšeme řešení úlohy (3.33) ve tvaru
x(t) = − sinh 2t−t0−1
argsinh
√
−x1
2
. (3.35)
Tento výsledek můžeme ještě upravit. Nejprve využijeme skutečnosti, že −x1 = 4x0(x0 − 1)
a proto
argsinh
√
−x1 = argsinh 4x2
0 − 4x0 = ln 4x2
0 − 4x0 + 4x2
0 − 4x0 + 1 =
= ln (2x0 − 1)2 − 1 + (2x0 − 1)2 = ln |1 − 2x0| + (2x0 − 1)2 − 1 =
= argcosh |1 − 2x0|;
3.4. ROVNICE ŘEŠITELNÉ SPECIÁLNÍMI SUBSTITUCEMI 59
dále využijeme vzorec pro hyperbolický sinus polovičního argumentu
sinh
α
2
2
=
cosh α − 1
2
a po dosazení dostaneme
x(t) = −1
2 cosh 2t−t0
argcosh |1 − 2x0| − 1 .
Odtud vyjádříme
1 − 2x(t) = cosh 2t−t0
argcosh |1 − 2x0| . (3.36)
Řešení úlohy (3.29) s počáteční hodnotou x0, splňující nerovnost x0 − 1
2 > 1
2 jsme pro t > t0
dostali v implicitním tvaru (3.36). Již snadno ověříme, že rovností
|1 − 2x(t)| = cosh 2t−t0
argcosh |1 − 2x0| (3.37)
je dáno řešení úlohy (3.29) s počáteční hodnotou splňující x0 − 1
2 > 1
2 pro každé t ≥ t0.
Přímým výpočtem můžeme také ukázat, že substitucí
|1 − 2x(t)| = cosh z(t)
dostaneme řešení úlohy (3.29) s počáteční hodnotou x0 splňující nerovnost x0 − 1
2 > 1
2 ve
tvaru (3.37).
Kvalitativní vlastnosti řešení úlohy (3.29) s počáteční hodnotou x0 splňující nerovnost
x0 − 1
2 > 1
2 . Z vyjádření řešení ve tvaru (3.35) vidíme, že pro libovolné řešení x úlohy platí
lim
t→∞
= −∞
a posloupnost x je ryze klesající.
Případ a = −2: Budeme řešit počáteční úlohu
x(t + 1) = −2x(t) 1 − x(t) , x(t0) = x0. (3.38)
Nejprve zavedeme substituci
y(t) = x(t) − 1
2, tj. x(t) = y(t) + 1
2. (3.39)
Po dosazení dostaneme
y(t + 1) = x(t + 1) − 1
2 = −2x(t) 1 − x(t) − 1
2 = −2 y(t) + 1
2 1 − y(t) − 1
2 − 1
2 =
= −2 y(t) + 1
2 −y(t) + 1
2 − 1
2 = 2 y(t)2
− 1
4 − 1
2 = 2y(t)2
− 1.
Posloupnost y je tedy řešením počáteční úlohy
y(t + 1) = 2y(t)2
− 1, y(t0) = x0 − 1
2.
To je první z rovnic řešitelná goniometrickou nebo hyperbolickou substitucí, viz 3.4.1. Řešení
rovnice hledáme pro různé počáteční hodnoty různými substitucemi.
60 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Nechť nejprve y(t0) = x0 − 1
2 ∈ [−1, 1], tj. x0 ∈ −1
2, 3
2 . Substitucí
y(t) = cos z(t) (3.40)
dostaneme řešení ve tvaru y(t) = cos 2t−t0 arccos y0 . Zpětnou substitucí dostaneme řešení
původní úlohy
x(t) = 1
2 + cos 2t−t0
arccos x0 − 1
2 .
Pokud |y(t0)| = x0 − 1
2 > 1, tj. x0 < −1
2 nebo x0 > 3
2, použijeme substituci
y(t) = cosh z(t) (3.41)
a dostaneme řešení ve tvaru
x(t) = 1
2 + cosh 2t−t0
ln |x0 − 1
2| + x0 − 1
2
2
− 1 .
Řešení logistické rovnice s parametrem a = −2 pomocí substitucí (3.39) a (3.40) nebo
(3.41) můžeme opět shrnout:
Řešení úlohy (3.38) s počáteční podmínkou x(t0) = x0 ∈ [−1
2, 3
2 ] najdeme substitucí
x(t) − 1
2 = cos z(t)
ve tvaru
x(t) − 1
2 = cos 2t−t0
arccos x0 − 1
2 ;
řešení úlohy (3.38) s počáteční podmínkou |x0 − 1
2 | > 1 najdeme substitucí
x(t) − 1
2 = cosh z(t)
ve tvaru
x(t) − 1
2 = cosh 2t−t0
argcosh x0 − 1
2 .
„Zobecňující poznámka: Logistickou rovnici
x(t + 1) = ax(t) 1 − x(t)
lze pro některé hodnoty parametru a a některé počáteční hodnoty x0 řešit substitucí
f x(t) = ϕ z(t) ,
kterou dostaneme řešení logistické rovnice v implicitním tvaru
f x(t) = ϕ 2t−t0
ϕ−1
f(x0) .
Hodnoty parametru a a počáteční hodnoty x0, pro které jsme tímto postupem našli řešení
logistické rovnice jsou shrnuty v tabulce 3.1.
3.5. CVIČENÍ 61
parametr počáteční hodnota f(ξ) ϕ(ζ)
a = 2 x0 ∈ R 1 − 2ξ eζ
a = 4 x0 ∈ [0, 1]
√
ξ | sin ζ|
a = 4 x0 > 1 nebo x0 < 0 |1 − 2ξ| cosh ζ
a = −2 x0 ∈ [−1
2, 3
2 ] ξ − 1
2 cos ζ
a = −2 x0 < −1
2 nebo x0 > 3
2 ξ − 1
2 cosh ζ
Tabulka 3.1: Hodnoty parametru a počáteční podmínky, pro které je řešení diskrétní logistické
rovnice x(t+1) = ax(t) 1−x(t) s počáteční podmínkou x(t0) = x0 implicitně dáno rovností
f x(t) = ϕ 2t−t0 ϕ−1 f(x0) .
3.5 Cvičení
V úlohách 1–6 najděte obecné řešení rovnice.
1. x(t + 1)2 − (2 + t)x(t + 1)x(t) − 2tx(t)2 = 0
2. x(t + 1)x(t) − x(t + 1) + x(t) = 0
3. x(t + 1)x(t) − 2
3x(t + 1) + 1
6 x(t) = 5
18
4. x(t + 1) = x(t)2
5. x(t + 1) = 2x(t) 1 − x(t)2
6. x(t + 1) = 1
2 x(t) −
a
x(t)
, a > 0
V úlohách 7–10 najděte řešení rovnice s počáteční podmínkou x(0) = x0.
7. x(t + 1)2 − 2x(t + 1)x(t) − 3x(t)2 = 0
8. x(t + 1) = 5 −
6
x(t)
9. x(t + 1) =
x(t) + a
x(t) + 1
, a > 0
10. x(t + 1) =
2 − x(t)2
2 1 − x(t)
11. Řešte počáteční úlohu x(t + 2) =
x(t + 1)3
x(t)2
, x(0) = x0, x(1) = x1.
12. Ukažte, že k libovolnému kladnému přirozenému číslu n existuje hodnota x0 = x0(n)
taková, že řešení úlohy
x(t + 1) = 2x(t) x(t) − 1 , x(0) = x0
má periodu n.
62 KAPITOLA 3. DALŠÍ EXPLICITNĚ ŘEŠITELNÉ ROVNICE
Výsledky:
1. 2t
c, x(t0)(−1)t−t0
t−1
i=t0
i,
2.
1
c − t
,
3.
5 − 2c(−6)t
6 1 + c(−6)t
, −1
3
4. c2t
5. sin 2t
c
6.
√
a cotg 2t
c
7. x03t
, x0(−1)t
8.
3x0 − 6
x0 − 2 + (x0 + 3) 2
3
t +
2x0 + 6
x0 + 3 + (x0 − 2) 3
2
t
9.
√
a
(x0 +
√
a) (1 +
√
a)
t
+ (x0 −
√
a) (1 −
√
a)
t
(x0 +
√
a) (1 +
√
a)
t
− (x0 −
√
a) (1 −
√
a)
t
10. 1 − cotg (2t
arccotg(1 − x0))
11. x0
x1
x0
2t
−1
12. Například x0 = 1
2 + cos
2π
2n − 1