03.08.2013 Views

Rychlý úvod do Greenových funkcí pro obyčejné diferenciální rovnice

Rychlý úvod do Greenových funkcí pro obyčejné diferenciální rovnice

Rychlý úvod do Greenových funkcí pro obyčejné diferenciální rovnice

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Greenova funkce <strong>pro</strong> dvoubo<strong>do</strong>vé<br />

okrajové úlohy <strong>pro</strong> <strong>obyčejné</strong> <strong>diferenciální</strong> <strong>rovnice</strong><br />

Jan Tomeček<br />

Tento stručný text si klade za cíl co nejrychlejší uvedení <strong>do</strong> teorie <strong>Greenových</strong> <strong>funkcí</strong> <strong>pro</strong> <strong>obyčejné</strong> <strong>diferenciální</strong><br />

<strong>rovnice</strong>. Autor tohoto textu byl inspirován textem M. Tvrdého [2], kde zaujatý čtenář najde více. Další lze<br />

nalézt např. v knize Obyčejné <strong>diferenciální</strong> <strong>rovnice</strong> [1]. Greenovy funkce <strong>pro</strong> obecnější okrajové úlohy lze nalézt<br />

např. v [3].<br />

Označení<br />

Symbolem R m×n budeme rozumět lineární <strong>pro</strong>stor všech matic typu (m, n) nad R. Speciálně, R n×1 (resp.<br />

R 1×n ) bude označovat množinu všech n–dimenzionálních sloupcových (resp. řádkových) vektorů.<br />

1 Soustava <strong>diferenciální</strong>ch rovnic prvního řádu<br />

Uvažujme systém <strong>diferenciální</strong>ch rovnic<br />

x ′ (t) + A(t)x(t) = b(t), (1)<br />

kde A ∈ C (0) ([α, β]; R n×n ), b ∈ C (0) ([α, β]; R n×1 ) spolu s okrajovými podmínkami<br />

kde M, N ∈ R m×n , m ≥ 1, n > 1.<br />

Mx(α) + Nx(β) = 0, (2)<br />

Definice 1 Řekneme, že x ∈ C (1) ([α, β]; R n×1 ) je řešením okrajové úlohy (1), (2), jestliže x splňuje soustavu<br />

<strong>diferenciální</strong>ch rovnic (1) <strong>pro</strong> každé t ∈ (α, β) a splňuje okrajové podmínky (2).<br />

Společně se soustavou (1) budeme pracovat s příslušnou homogenní soustavou, tzn. se soustavou (1) <strong>pro</strong><br />

b ≡ 0, tedy<br />

x ′ (t) + A(t)x(t) = 0. (3)<br />

Uvažujme fundamentální systém řešení soustavy (3), tzn. maticovou funkci X ∈ C (1) ([α, β]; R n×n ) jejíž<br />

sloupce tvoří lineárně nezávislou množinu řešení systému (3) na intervalu [α, β], tzn. platí<br />

Pak <strong>pro</strong> dané τ ∈ [α, β] je funkce<br />

X ′ (t) + A(t)X(t) = 0 <strong>pro</strong> každé t ∈ [α, β].<br />

x(t) = X(t)c + X(t)<br />

t<br />

s parametrem c ∈ R n×1 obecným řešením systému (1).<br />

τ<br />

X −1 (s)b(s) ds, <strong>pro</strong> každé t ∈ [α, β] (4)<br />

Poznámka 2 Z (4) mj. plyne, že každé řešení soustavy (3) (na nějakém intervalu) je ve tvaru X(t)c, kde<br />

c ∈ R n×1 . A z toho také vidíme, že maticová funkce typu n × n, jejíž sloupce jsou řešeními soustavy (3) je ve<br />

tvaru X(t)C, kde C ∈ R n×n .<br />

1


Ve zbytku kapitoly budeme uvažovat jeden fundamentální systém řešení X ∈ C (1) ([α, β]; R n×n ) soustavy (3)<br />

(i když to v našich úvahách nebude příliš důležité).<br />

Lemma 3 Okrajová úloha (1), (2) má řešení právě tehdy, když má řešení soustava algebraických rovnic<br />

β<br />

(MX(α) + NX(β))c = −NX(β) X<br />

α<br />

−1 (s)b(s) ds. (5)<br />

Navíc, řešení x úlohy (1), (2) je dáno ve tvaru (4), kde c ∈ R n×1 je řešením soustavy (5).<br />

Důkaz. Z předchozích úvah je zřejmé, že hledané řešení úlohy (1), (2) bude ve tvaru (4) <strong>pro</strong> τ = α. Dosazením<br />

funkce x z (4) <strong>do</strong> okrajové podmínky (2) <strong>do</strong>stáváme<br />

<br />

β<br />

MX(α)c + N X(β)c + X(β) X<br />

α<br />

−1 <br />

(s)b(s) ds = 0.<br />

Jednoduchou úpravou <strong>do</strong>stáváme soustavu (5) a tím i tvrzení lemmatu. <br />

Poznámka 4 Okrajovou úlohu (1), (2) jsme tedy převedli na úlohu řešit soustavu algebraických rovnic (5).<br />

Bude se nám hodit následující tvrzení z lineární algebry:<br />

Soustava algebraických rovnic<br />

Dc = w, (D ∈ R m×n , w ∈ R m×1 )<br />

je řešitelná právě tehdy, když<br />

Z Lemmatu 3 a Poznámky 4 tedy vyplývá následující.<br />

(∀γ T ∈ R 1×m )(γ T D = 0 ⇒ γ T w = 0).<br />

Lemma 5 Okrajová úloha (1), (2) je řešitelná právě tehdy když <strong>pro</strong> každé γ T ∈ R 1×m splňující<br />

platí<br />

γ T (MX(α) + NX(β)) = 0<br />

γ T β<br />

NX(β) X<br />

α<br />

−1 (s)b(s) ds = 0.<br />

Bude se nám také hodit následující jednoduché tvrzení.<br />

Lemma 6 Nechť A ∈ C (0) ([α, β]; R n×n ), M, N ∈ R m×n . Pak jsou následující tvrzení ekvivalentní:<br />

(i) Pro každé b ∈ C (0) ([α, β]; R n×1 ) má úloha (1), (2) jediné řešení.<br />

(ii) Úloha (3), (2) má jen triviální řešení (tzn. nulová funkce x(t) ≡ 0 <strong>pro</strong> každé t ∈ [α, β])<br />

Důkaz. Implikace (i) ⇒ (ii) je zřejmá z faktu, že úloha (3), (2) je speciálním případem úlohy (1), (2) <strong>pro</strong> b ≡ 0<br />

na [α, β]. Tedy úloha (3), (2) je podle předpokladu jednoznačně řešitelná, přitom má vždy alespoň triviální<br />

řešení – je tedy jediné.<br />

Implikaci (ii) ⇒ (i) <strong>do</strong>kážeme násle<strong>do</strong>vně. Nechť <strong>pro</strong> dané b jsou x1, x2 : [α, β] → R n×1 řešení úlohy (1), (2).<br />

Snadno se ukáže, že funkce x1 − x2 je řešením úlohy (3), (2). Z předpokladu (ii) plyne, že x1 − x2 je nulová<br />

funkce na [α, β], tedy x1 = x2 na [α, β]. <br />

Lemma 7 Nechť A ∈ C (0) ([α, β]; R n×n ), M, N ∈ R m×n , h((M, N)) = m. Pak tvrzení (i) a (ii) z Lemmatu 6<br />

jsou ekvivalentní s předpokladem<br />

m = n a det(MX(α) + NX(β)) = 0. (6)<br />

2


Důkaz. Nejprve <strong>do</strong>kážeme implikaci (6) ⇒ (i). Ta plyne z faktu, že za předpokladu platnosti (6) je soustava<br />

(5) jednoznačně řešitelná. Její řešení je ve tvaru<br />

c = −(MX(α) + NX(β)) −1 β<br />

NX(β) X<br />

α<br />

−1 (s)b(s) ds<br />

a tedy řešení úlohy (1), (2) je jednoznačně určeno vztahem (4) <strong>pro</strong> c z předchozí rovnosti.<br />

Dokážeme implikaci (i) ⇒ (6). Z předpokladu (i) a Lemmatu 3 plyne, že soustava<br />

(MX(α) + NX(β))c = 0<br />

má pouze triviální řešení. Pak (z faktu, že dimenze <strong>pro</strong>storu řešení soustavy je rovna rozdílu počtu sloupců<br />

matice soustavy a její hodnosti)<br />

0 = n − h(MX(α) + NX(β)).<br />

Protože matice MX(α) + NX(β) má m řádků, pak zřejmě<br />

m ≥ h(MX(α) + NX(β)) = n.<br />

Dokážeme, že m = n, sporem. Předpokládejme, že m > n. Pak řádky matice MX(α) + NX(β) jsou lineárně<br />

závislé, tzn. existuje γ T ∈ R 1×m tak, že<br />

γ T (MX(α) + NX(β)) = 0. (7)<br />

Z faktu γ T MX(α) = −γ T NX(β) a regularity matic X(α) a X(β) plyne ekvivalence<br />

Z faktu h((M, N)) = m plyne<br />

Odtud a z (8) plyne γ T N = 0. Položme<br />

Pak<br />

γ T M = 0 ⇐⇒ γ T N = 0. (8)<br />

0 = γ T (M, N) = (γ T M, γ T N).<br />

b(t) = X −1 (t)X −1 (β)(γ T N) T<br />

∀t ∈ [α, β].<br />

γ T β<br />

NX(β) X<br />

α<br />

−1 (s)b(s) ds = γ T N(γ T N) T = 0.<br />

Z (i), (7) a Lemmatu 5 <strong>do</strong>stáváme spor. Dále podle předpokladu (i) má okrajová úloha (1), (2) <strong>pro</strong> každé<br />

b ∈ C (0) ([α, β]; R n×1 ) jediné řešení, což je podle Lemmatu 3 ekvivalentní s faktem<br />

det(MX(α) + NX(β)) = 0.<br />

Definice 8 Platí–li jedna z podmínek (i), (ii) z Lemmatu 6 nebo podmínka (6), řekneme, že úloha (1), (2) není<br />

v rezonanci.<br />

Lemma 9 Nechť úloha (1), (2) není v rezonanci. Pak <strong>pro</strong> každé b ∈ C (0) ([α, β]; Rn×1 ) lze řešení x úlohy (1),<br />

(2) napsat ve tvaru<br />

kde<br />

G(t, s) =<br />

kde D = MX(α) + NX(β).<br />

x(t) =<br />

β<br />

α<br />

G(t, s)b(s) ds, t ∈ [α, β] (9)<br />

X(t)(−D −1 NX(β))X −1 (s) <strong>pro</strong> α ≤ t ≤ s ≤ β,<br />

X(t)(En − D −1 NX(β))X −1 (s) <strong>pro</strong> α ≤ s < t ≤ β,<br />

3<br />

<br />

(10)


Důkaz. Z Lemmat 3 a 7 plyne, že řešení úlohy (1), (2) je ve tvaru<br />

x(t) = −X(t)D −1 β<br />

NX(β)<br />

= −X(t)D −1 NX(β)<br />

=<br />

=<br />

t<br />

α<br />

β<br />

α<br />

α<br />

t<br />

α<br />

X −1 (s)b(s) ds + X(t)<br />

t<br />

α<br />

X −1 (s)b(s) ds<br />

X −1 (s)b(s) ds − X(t)D −1 NX(β)<br />

X(t)(En − D −1 NX(β))X −1 (s)b(s) ds +<br />

G(t, s)b(s) ds, <strong>pro</strong> t ∈ [α, β],<br />

β<br />

t<br />

β<br />

t<br />

X −1 (s)b(s) ds + X(t)<br />

X(t)(−D −1 NX(β))X −1 (s)b(s) ds<br />

t<br />

α<br />

X −1 (s)b(s) ds<br />

kde G je dána vzorcem (10). <br />

V našich dalších úvahách budeme často pracovat s množinou<br />

△ = {(t, t) ∈ R 2 : t ∈ [α, β]}.<br />

Definice 10 Nechť úloha (1), (2) není v rezonanci. Funkci G : [α, β] 2 → R n×n spojitou a omezenou na [α, β] 2 \△<br />

takovou, že <strong>pro</strong> každé b ∈ C (0) ([α, β]; R n×1 ) lze řešení x úlohy (1), (2) psát ve tvaru (9) nazýváme Greenovou<br />

<strong>funkcí</strong> úlohy (1), (2).<br />

Poznámka 11 Z definice je vidět, že Greenova funkce není definovaná jednoznačně, <strong>pro</strong>tože na množině △<br />

může nabývat jakýchkoliv hodnot. Z Lemmatu 9 je vidět, že existuje vždy alespoň jedna Greenova funkce dané<br />

úlohy (<strong>pro</strong> úlohu, která není v rezonanci). V dalším lemmatu <strong>do</strong>kážeme, že každá Greenova funkce je dána<br />

jednoznačně ve tvaru (10) až na funkční hodnoty bodů z množiny △, a nezávisí na volbě fundamentální matice<br />

X.<br />

Lemma 12 Nechť úloha (1), (2) není v rezonanci. Greenova funkce je definována jednoznačně až na množinu<br />

△.<br />

Důkaz. Mějme Greenovy funkce G1, G2 : [α, β] 2 → Rn×n úlohy (1), (2). Pak <strong>pro</strong> libovolné b ∈ C (0) ([α, β]; Rn×1 )<br />

platí<br />

x(t) =<br />

β<br />

α<br />

G1(t, s)b(s) ds =<br />

je jediným řešením úlohy (1), (2). Z toho plyne<br />

β<br />

a <strong>pro</strong> každé t ∈ [α, β]. Dokážeme, že<br />

α<br />

β<br />

α<br />

G2(t, s)b(s) ds, t ∈ [α, β]<br />

(G1(t, s) − G2(t, s))b(s) ds = 0 <strong>pro</strong> každé b ∈ C (0) ([α, β]; R n×1 ) (11)<br />

G1 = G2 na [α, β] 2 \ △. (12)<br />

Zvolme t ∈ (α, β) (důkaz <strong>pro</strong> t ∈ {α, β} se <strong>pro</strong>vede analogicky). Dokážeme, že z (11) plyne G1(t, s) = G2(t, s)<br />

<strong>pro</strong> všechna s ∈ [α, t) ∪ (t, β]. Podle definice jsou funkce G1(t, ·) a G2(t, ·) spojité na [α, t) a (t, β], tedy to samé<br />

platí i <strong>pro</strong> jejich rozdíl G1(t, ·) − G2(t, ·). Nechť i ∈ {1, . . . , n}. Označme i–tý sloupec matice G1(t, s) − G2(t, s)<br />

jako gi(t, s). Dokážeme, že gi(t, s) = 0 <strong>pro</strong> každé s ∈ [α, t) ∪ (t, β].<br />

Zřejmě existuje p0 ∈ N tak, že<br />

t − 1<br />

1<br />

> α a t + < β <strong>pro</strong> každé p ≥ p0.<br />

p p<br />

Pro každé p ∈ N, p ≥ p0 definujeme vektorovou funkci bp ∈ C (0) ([α, β]; R n×1 ) takovou, že<br />

bp(s) = gi(t, s) T , s ∈ [α, t − 1 1<br />

] ∪ [t + , β],<br />

p p<br />

4


jejíž složky jsou na množině (t − 1 1<br />

p , t + p ) lineární (snadno lze určit funkční předpis, je ovšem složitější na zápis).<br />

Z omezenosti <strong>funkcí</strong> G1, G2 plyne omezenost funkce gi na [α, β] 2 \ △. Z toho plyne, že existuje K > 0 tak, že<br />

Z (11) tedy <strong>pro</strong> každé p ≥ p0 vyplývá<br />

Protože 1<br />

t+ p<br />

t− 1<br />

p<br />

|gi(t, s)| < K a |bp(s)| ≤ K, s ∈ [α, t) ∪ (t, β], p ≥ p0.<br />

0 =<br />

+<br />

β<br />

α<br />

1 t+ p<br />

t− 1<br />

p<br />

gi(t, s)bp(s) ds =<br />

<br />

<br />

gi(t, s)bp(s) ds<br />

≤<br />

t− 1<br />

p<br />

α<br />

gi(t, s) · bp(s) ds +<br />

t+ 1<br />

p<br />

t− 1<br />

p<br />

β<br />

gi(t, s) · g T i (t, s) ds<br />

t+ 1<br />

p<br />

gi(t, s) · g T i (t, s) ds.<br />

K 2 ds = K 2<br />

<br />

t + 1<br />

<br />

1<br />

− t + =<br />

p p<br />

2K2<br />

p<br />

<strong>pro</strong> každé p ≥ p0, pak limitním přechodem v (13) <strong>pro</strong> p → ∞ <strong>do</strong>staneme<br />

0 =<br />

β<br />

α<br />

gi(t, s) · g T i (t, s) ds.<br />

Z toho a spojitosti gi(t, s) · gT i (t, s) na intervalech [α, t) a (t, β] plyne<br />

a tedy<br />

gi(t, s) · g T i (t, s) = 0 <strong>pro</strong> každé s ∈ [α, t) ∪ (t, β]<br />

gi(t, s) = 0 <strong>pro</strong> každé s ∈ [α, t) ∪ (t, β].<br />

Dokázali jsme tedy rovnost <strong>funkcí</strong> G1 a G2 na množině [α, β] 2 \ △. Z tohoto faktu a spojitosti G1 a G2 na<br />

[α, β] 2 \ △ již plyne (12). <br />

Následující věta udává vlastnosti Greenovy funkce, kterými je <strong>do</strong>konce charakterizována.<br />

Věta 13 Nechť úloha (1), (2) není v rezonanci. Pak maticová funkce G : [α, β] 2 → R n×n je Greenovou <strong>funkcí</strong><br />

úlohy (1), (2) právě tehdy, když<br />

(I) Funkce G je spojitá a omezená na [α, β] 2 \ △.<br />

(II) Pro každé s ∈ (α, β) platí<br />

lim G(t, s) − lim G(t, s) = En.<br />

t→s+ t→s−<br />

(III) Pro každé s ∈ (α, β) jsou sloupce maticové funkce G(·, s) řešením soustavy (3) na intervalech (α, s) a<br />

(s, β).<br />

(IV) Pro každé s ∈ (α, β) platí<br />

MG(α, s) + NG(β, s) = 0.<br />

Důkaz. Nejprve <strong>do</strong>kážeme nutnost podmínek. Nechť G je Greenova funkce úlohy. Dokážeme, že pak splňuje<br />

(I)–(IV). Vlastnost (I) je splněna okamžitě z definice. Vzhledem k Poznámce 11 a Lemmatu 12 můžeme bez<br />

újmy na obecnosti psát (10), kde X je fundamentální matice soustavy (3) na [α, β]. Pak <strong>pro</strong> každé s ∈ (α, β)<br />

platí (díky spojitosti X −1 na [α, β])<br />

lim G(t, s) − lim<br />

t→s+ t→s− G(t, s) = X(t)(En − D −1 NX(β))X −1 (t)<br />

− X(t)(−D −1 NX(β))X −1 (t) = En,<br />

5<br />

(13)


tedy (II) platí. Vlastnost (III) plyne z (10). Zbývá ověřit vlastnost (IV). Pro s ∈ (α, β) platí<br />

MG(α, s) + NG(β, s)<br />

= MX(α)(−D −1 NX(β))X −1 (s) + NX(β)(I − D −1 NX(β))X −1 (s)<br />

= −MX(α)D −1 NX(β)X −1 (s) + NX(β)X −1 (s) − NX(β)D −1 NX(β)X −1 (s)<br />

= −(MX(α) + NX(β))D −1 NX(β)X −1 (s) + NX(β)X −1 (s)<br />

= −NX(β)X −1 (s) + NX(β)X −1 (s)<br />

= 0.<br />

Nyní ověříme postačitelnost podmínek (I)–(IV). Nechť G : [α, β] 2 → R n×n je maticová funkce splňující (I)–<br />

(IV). Dokážeme, že jde o Greenovu funkci úlohy (1), (2). Z podmínky (III) a Poznámky 2 plyne, že existují<br />

maticové funkce C1, C2 : (a, b) → R n×n tak, že<br />

G(t, s) =<br />

Z podmínky (II) <strong>do</strong>stáváme <strong>pro</strong> s ∈ (α, β)<br />

a z (IV) plyne<br />

Z těchto dvou maticových rovnic vyjádříme<br />

X(t)C1(s) <strong>pro</strong> α < s < t ≤ β,<br />

X(t)C2(s) <strong>pro</strong> α ≤ t < s < β.<br />

En = lim G(t, s) − lim<br />

t→s+ t→s− G(t, s) = X(s)C2(s) − X(s)C1(s).<br />

MX(α)C1(s) + NX(β)C2(s) = 0.<br />

C1(s) = −D −1 NX(β)X −1 (s) a C2(s) = (En − D −1 NX(β))X −1 (s)<br />

<strong>pro</strong> s ∈ (α, β). Z (I) a spojitosti X −1 na [α, β] také plyne<br />

C1(β) = −D −1 NX(β)X −1 (β) a C2(α) = (En − D −1 NX(β))X −1 (α).<br />

Tímto je tedy funkce G určena jednoznačně na množině [α, β] 2 \ △ na níž je rovna funkci z (10). Vzhledem k<br />

Lemmatu 9 a Lemmatu 12 je důkaz hotov. <br />

Poznámka 14 Označme G(t, s) = (gij(t, s)) n i,j=1 Greenovu funkci úlohy (1), (2). Z vlastnosti (II) z Věty 13 je<br />

vidět, že funkce gij <strong>pro</strong> i = j lze spojitě <strong>do</strong>definovat na celou množinu [α, β] 2 . Funkce gii vzhledem k (II) nelze<br />

<strong>do</strong>definovat spojitě, ale spojitě lze <strong>do</strong>definovat na uzávěr její restrikce vzhledem k množinám {(t, s) ∈ R2 : α ≤<br />

s < t ≤ β} a {(t, s) ∈ R2 : α ≤ t < s ≤ β}.<br />

2 Diferenciální <strong>rovnice</strong> n–tého řádu<br />

Uvažujme lineární <strong>diferenciální</strong> rovnici n–tého řádu<br />

a0(t)u (n) + a1(t)u (n−1) + a2(t)u (n−2) + . . . + an(t)u = f(t), (14)<br />

kde ai ∈ C (0) ([α, β]; R), i = 0, . . . , n, a0(t) = 0 <strong>pro</strong> všechna t ∈ [α, β], f ∈ C (0) ([α, β]; R) spolu s okrajovou<br />

podmínkou<br />

Mξ(α) + Nξ(β) = 0 (15)<br />

kde M, N ∈ R m×n , ξ(t) = (u(t), u ′ (t), . . . , u (n−1) ) T <strong>pro</strong> t ∈ {α, β}. Ze stejného důvodu jako v předchozí kapitole<br />

si můžeme <strong>do</strong>volit předpokládat<br />

h((M, N)) = m.<br />

V celé kapitole budeme mít tedy tyto předpoklady<br />

ai ∈ C (0) ([α, β]; R), i = 0, . . . , n, a0(t) = 0 <strong>pro</strong> ∀t ∈ [α, β],<br />

f ∈ C (0) ([α, β]; R), M, N ∈ R m×n , h((M, N)) = m.<br />

6


Definice 15 Řekneme, že u ∈ C (n) ([α, β]; R) je řešením okrajové úlohy (14), (15), jestliže u splňuje <strong>diferenciální</strong><br />

rovnici (14) <strong>pro</strong> t ∈ (α, β) a splňuje okrajové podmínky (15).<br />

Společně s rovnicí (14) budeme pracovat s příslušnou homogenní rovnicí tzn. s rovnicí (14) <strong>pro</strong> f ≡ 0, tedy<br />

a0(t)u (n) + a1(t)u (n−1) + a2(t)u (n−2) + . . . + an(t)u = 0. (16)<br />

Lemma 16 Nechť u je řešení <strong>rovnice</strong> (14). Pak vektorová funkce<br />

ψu(t) = (u(t), u ′ (t), . . . , u (n−1) (t)) T , t ∈ [α, β]<br />

je řešením soustavy (1) na [α, β] <strong>pro</strong><br />

⎛<br />

0 1 0 . . . 0<br />

⎜<br />

0 0 1 . . . 0<br />

A(t) = ⎜<br />

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .<br />

⎝ 0 0 0 . . . 1<br />

− an an−1<br />

− − a0 a0<br />

an−2<br />

. . . − a0<br />

a1<br />

⎞ ⎛ ⎞<br />

0<br />

⎟ ,<br />

⎠<br />

a0<br />

⎜<br />

0 ⎟<br />

b(t) = ⎜<br />

. . ⎟ .<br />

⎝ 0 ⎠<br />

f<br />

(17)<br />

Zobrazení ψ množiny všech řešení <strong>rovnice</strong> (14) definovaných na intervalu J ⊂ [α, β] na množinu všech řešení<br />

soustavy (1) definovaných na intervalu J ⊂ [α, β] <strong>pro</strong> (17) je izomorfismus.<br />

Řešení u <strong>rovnice</strong> (14) splňuje okrajové podmínky (15) právě tehdy, když ψu splňuje okrajové podmínky (2).<br />

Definice 17 Nechť úloha (14), (15) není v rezonanci. Funkci g ∈ C (0) ([α, β] 2 ; R) nazveme Greenovou <strong>funkcí</strong><br />

úlohy (14), (15), jestliže <strong>pro</strong> každé f ∈ C (0) ([α, β]; R) je řešení úlohy (14), (15) ve tvaru<br />

u(t) =<br />

β<br />

α<br />

g(t, s)f(s) ds ∀t ∈ [α, β]. (18)<br />

Lemma 18 Nechť úloha (14), (15) není v rezonanci. Pak existuje jediná Greenova funkce g úlohy (14), (15).<br />

Je charakterizována vlastnostmi<br />

(I’) Funkce g, ∂g<br />

∂t , . . ., ∂n−2 g<br />

∂t n−2 jsou spojité a omezené na [α, β] 2 , funkce ∂n−1 g<br />

∂t n−1 je spojitá a omezená na [α, β] 2 \△.<br />

(II’) Pro každé s ∈ (α, β) platí<br />

lim<br />

t→s+<br />

∂n−1g (t, s) − lim<br />

∂tn−1 t→s−<br />

∂n−1g 1<br />

(t, s) =<br />

∂tn−1 a0(s) .<br />

(III’) Pro každé s ∈ (α, β) je funkce g(·, s) řešením soustavy (16) na intervalech (α, s) a (s, β).<br />

(IV’) Pro každé s ∈ (α, β) platí<br />

M(g(α, s), ∂g<br />

∂t (α, s), . . . , ∂n−1g (α, s))T<br />

∂tn−1 + N(g(β, s), ∂g<br />

∂t (β, s), . . . , ∂n−1g ∂tn−1 (β, s))T = 0.<br />

Důkaz. Uvažujme úlohu (14), (15). Podle Lemmatu 16 je tato úloha ekvivalentní s úlohou (1), (2) <strong>pro</strong> (17).<br />

Protože podle předpokladu úloha (14), (15) není v rezonanci plyne opět z Lemmatu 16, že příslušná úloha (1),<br />

(2) (<strong>pro</strong> (17)) také není v rezonanci. Pro dané f a z ní odvozené b (z (17)) existuje jediné řešení x úlohy (1),<br />

(2), které lze zapsat pomocí Greenovy funkce G úlohy (1), (2) ve tvaru (9). Z Lemmatu 16 plyne existence<br />

jednoznačně určeného řešení u úlohy (1), (2) tak, že<br />

x(t) = ψu(t) = (u(t), u ′ (t), . . . , u (n−1) (t)) T , t ∈ [α, β].<br />

Označíme–li G(t, s) = (gij(t, s)) n,n<br />

i,j=1 , pak <strong>pro</strong> první řádek vektorové rovnosti (9) platí<br />

u(t) = x1(t) =<br />

β<br />

α<br />

n<br />

g1i(t, s)bi(s) ds =<br />

i=1<br />

7<br />

β<br />

α<br />

a0<br />

g1n(t, s) f(s)<br />

a0(s) ds


(poslední rovnost plyne z (17)). Shrneme–li naše úvahy, k funkci f jsme našli jednoznačně určené řešení u úlohy<br />

(14), (15) ve tvaru (18). Je vidět, že Greenovou <strong>funkcí</strong> úlohy (14) a (15) je funkce<br />

g(t, s) = g1,n(t, s)<br />

, (t, s) ∈ [α, β]<br />

a0(s)<br />

2 ,<br />

kterou můžeme s ohledem k Poznámce 14 chápat jako spojitou na [α, β] 2 . Dokážeme, že je charakterizována<br />

vlastnostmi (I)–(IV). Nejprve nutnost. Vzhledem k Poznámce 14 lze považovat funkce g1,n, g2,n, . . . , gn−1,n<br />

za spojité funkce na [α, β] 2 a funkci gn,n jako spojitou na [α, β] 2 \ △, přitom <strong>pro</strong> každé s ∈ (α, β) je vektorová<br />

funkce<br />

(g1,n(·, s), g2,n(·, s), . . . , gn,n(·, s)) T<br />

řešením systému (3) na intervalech (α, s) a (s, β) <strong>pro</strong> (17), z jehož speciálního tvaru plyne<br />

Tedy<br />

gi+1,n(t, s) = ∂i g1,n<br />

∂t i (t, s), (t, s) ∈ (α, β) 2 \ △, i = 1, . . . , n − 1.<br />

∂ig ∂ti (t, s) = gi+1,n(t, s)<br />

, (t, s) ∈ [α, β]<br />

a0(s)<br />

2 , i = 1, . . . , n − 2<br />

a ∂n−1g ∂tn−1 (t, s) = gn,n(t,s)<br />

a0(s) <strong>pro</strong> (t, s) ∈ [α, β] 2 \ △ (v krajních bodech se rozumějí jednostranné derivace). Z toho<br />

plyne, že g, ∂g<br />

∂t , . . ., ∂n−2g ∂tn−2 jsou spojité na [α, β] 2 , ∂n−1g ∂tn−1 je omezená a spojitá na [α, β] 2 \ △ (mimo jiné také díky<br />

předpokladu nenulovosti a spojitosti a0 na kompaktním intervalu). Tedy (I’) platí. Dále z Poznámky 14 plyne<br />

lim<br />

t→s+<br />

∂n−1g (t, s) − lim<br />

∂tn−1 t→s−<br />

∂n−1g 1<br />

(t, s) =<br />

∂tn−1 a0(s)<br />

<br />

lim<br />

t→s+ gn,n(t, s) − lim<br />

t→s− gn,n(t,<br />

<br />

s)<br />

což je (II’). Vlastnost (III’) plyne z homogenity <strong>rovnice</strong> (16). Stačí již ověřit jen (IV’). Platí<br />

= 1<br />

a0(s) ,<br />

M(g(α, s), ∂g<br />

∂t (α, s), . . . , ∂n−1g ∂tn−1 (α, s))T + N(g(β, s), ∂g<br />

∂t (β, s), . . . , ∂n−1g (β, s))T<br />

∂tn−1 = M 1<br />

a0(s) (g1,n(α, s), g2,n(α, s), . . . , gn,n(α, s)) T + N 1<br />

a0(s) (g1,n(β, s), g2,n(β, s), . . . , gn,n(β, s)) T<br />

= 0.<br />

Poslední rovnost plyne z faktu, že vektory (g1,n(α, s), g2,n(α, s), . . . , gn,n(α, s)) T a (g1,n(β, s), g2,n(β, s), . . . , gn,n(β, s)) T<br />

jsou sloupce matic G(α, s) a G(β, s) a z vlastnosti (IV) z Věty 13. Nyní <strong>do</strong>kážeme, že každá funkce g splňující<br />

(I’)–(IV’) je Greenova funkce úlohy (14), (15). Mějme takovou funkci g dánu. Definujme<br />

˜g1,n(t, s) = a0(s)g(t, s), (t, s) ∈ [α, β] 2 .<br />

Vzhledem k podmínce (I’) existují spojité parciální derivace<br />

∂ i ˜g1,n<br />

∂ti (t, s) = a0(s) ∂ig (t, s)<br />

∂ti <strong>pro</strong> i = 1, . . . , n − 1 na příslušných množinách (viz (I’)), označíme je ˜gi+1,n(t, s). Nechť G označuje Greenovu<br />

funkci ze začátku tohoto důkazu (jde o Greenovu funkci úlohy (1), (2)). Označme symbolem ˜ G maticovou<br />

funkci, která vznikne z maticové funkce G nahrazením n–tého sloupce vektorovou <strong>funkcí</strong><br />

(˜g1,n, ˜g2,n, . . . , ˜gn,n) T .<br />

Z vlastností (I’)–(IV’) <strong>pro</strong> funkci g jednoduše plyne, že funkce ˜ G splňuje podmínky (I)–(IV), tedy jde o Greenovu<br />

funkci úlohy (1), (2). Protože je Greenova funkce určena jednoznačně (až na množinu △), je jasné, že G = ˜ G (až<br />

na množinu △) a tedy ˜g1,n = g1,n (vzhledem k Poznámce 14 můžeme brát obě funkce jako spojité na [α, β] 2 ).<br />

Z první části důkazu víme, že právě g1,n(t,s)<br />

a0(s) je Greenova funkce úlohy (14), (15) což není nic jiného než funkce<br />

g. <br />

8


References<br />

[1] Kurzweil, J., Obyčejné <strong>diferenciální</strong> <strong>rovnice</strong>, Praha, 1978.<br />

[2] Lineární okrajové úlohy [online], <strong>do</strong>stupné z<br />

http://www.math.cas.cz/ tvrdy/okrul.ps [citováno 21. 9. 2009]<br />

[3] Kiguradze, I., Okrajové úlohy <strong>pro</strong> systémy lineárních obyčejných <strong>diferenciální</strong>ch rovnic, Brno, 1997.<br />

9

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!