R´
ownania R´
o ˙zniczkowe Zwyczajne
wyk lad dla student´
ow na kierunku automatyka i robotyka - wersja robocza (6 stycze´
n 2006
)
Bogus law Bo˙zek
Wydzia l Matematyki Stosowanej AGH
1
Spis tre´
sci
Rozdzia l 1. Wprowadzenie
. . . . . . . . . . . . . . . . . . . . . . . . . . .
5
Rozdzia l 2. Elementy analizy funkcjonalnej
. . . . . . . . . . . . . . .
9
Rozdzia l 3. Twierdzenia o istnieniu i jednoznaczno´
sci
. . . . . . . .
11
Rozdzia l 4. Proste typy r´
owna´
n r´
o ˙zniczkowych skalarnych
. . . . .
15
4.1.
R´
ownanie r´
o˙zniczkowe o zmiennych rozdzielonych . . . . . . . . . . . . .
15
4.2.
R´
ownanie jednorodne . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
16
4.3.
R´
ownanie r´
o˙zniczkowe zupe lne . . . . . . . . . . . . . . . . . . . . . . . .
17
4.3.1. Czynnik ca lkuj
֒
acy . . . . . . . . . . . . . . . . . . . . . . . . . .
18
4.4.
R´
ownanie Clairauta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
19
Rozdzia l 5. Liniowe r´
ownania r´
o ˙zniczkowe
. . . . . . . . . . . . . . . .
21
5.1.
R´
ownania i uk lady r´
owna´
n r´
o˙zniczkowych liniowych . . . . . . . . . . . .
21
5.2.
Skalarne r´
ownanie liniowe rz
֒
edu pierwszego . . . . . . . . . . . . . . . .
23
5.3.
R´
ownanie Bernoulliego . . . . . . . . . . . . . . . . . . . . . . . . . . . .
24
5.4.
R´
ownanie Riccatiego . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
24
5.5.
R´
ownanie Lagrange’a . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
25
5.6.
Skalarne r´
ownanie r´
o˙zniczkowe liniowe
n
-tego rz
֒
edu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
26
5.7.
Obni˙zanie rz
֒
edu r´
ownania liniowego . . . . . . . . . . . . . . . . . . . . .
29
5.7.1. Wz´
or Liouville’a . . . . . . . . . . . . . . . . . . . . . . . . . . .
29
5.7.2. R´
ownania wy˙zszych rz
֒
ed´
ow . . . . . . . . . . . . . . . . . . . . .
30
5.8.
Niejednorodne r´
ownanie r´
o˙zniczkowe liniowe
n
-tego rz
֒
edu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
30
5.9.
R´
ownanie liniowe n-tego rz
֒
edu o sta lych wsp´
o lczynnikach . . . . . . . . .
31
5.10. Metoda przewidywa´
n . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
32
5.11. Uk lad skalarnych r´
owna´
n r´
o˙zniczkowych liniowych rz
֒
edu pierwszego . . .
33
5.12. Uk lady r´owna´
n liniowych o sta lych
wsp´
o lczynnikach . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
37
5.12.1. Metoda warto´sci i wektor´
ow w lasnych . . . . . . . . . . . . . . .
37
5.12.2. Sprowadzanie macierzy uk ladu do postaci Jordana . . . . . . . .
38
5.13. R´
ownanie ruchu harmonicznego . . . . . . . . . . . . . . . . . . . . . . .
42
3
Spis tre´sci
Rozdzia l 6. Rozwi
֒
azania w postaci szereg´
ow funkcyjnych
. . . . . .
43
6.1.
Rozwi
֒
azania w postaci szereg´
ow pot
֒
egowych . . . . . . . . . . . . . . . .
43
6.1.1. Uk lad r´
owna´
n liniowych rz
֒
edu pierwszego o sta lych
wsp´
o lczynnikach . . . . . . . . . . . . . . . . . . . . . . . . . . .
44
6.1.2. Skalarne r´
ownania r´
o˙zniczkowe rz
֒
edu pierwszego i drugiego . . .
44
6.2.
R´
ownania r´
o˙zniczkowe liniowe rz
֒
edu drugiego – szeregi Frobeniusa
. . .
46
Rozdzia l 7. Stabilno´
s´
c rozwi
֒
aza´
n r´
owna´
n r´
o ˙zniczkowych
. . . . . . .
49
7.1.
Podstawowe definicje . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
49
7.2.
Twierdzenie Lapunowa . . . . . . . . . . . . . . . . . . . . . . . . . . . .
51
7.3.
Problem Routha–Hurwitza . . . . . . . . . . . . . . . . . . . . . . . . . .
53
7.4.
Punkty osobliwe r´
ownania r´
o˙zniczkowego zupe lnego . . . . . . . . . . . .
54
Rozdzia l 8. Transformata Laplace’a
. . . . . . . . . . . . . . . . . . . . .
57
8.1.
Podstawowe definicje i twierdzenia
. . . . . . . . . . . . . . . . . . . . .
57
8.2.
Wyznaczanie transformaty r´
ownania r´
o˙zniczkowego . . . . . . . . . . . .
58
8.3.
Wyznaczanie funkcji na podstawie jej transformaty . . . . . . . . . . . .
59
Rozdzia l 9. Dodatek
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
63
9.1.
Tablice transformat Laplace’a . . . . . . . . . . . . . . . . . . . . . . . .
63
9.2.
Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne . . . . . .
66
9.3.
Przyk ladowe tematy zada´
n egzaminacyjnych – studia zaoczne . . . . . .
83
Bibliografia
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
93
Rozdzia l 1
Wprowadzenie
R´ownaniem r´o˙zniczkowym nazywamy zwi
֒
azek mi
֒
edzy pewn
֒
a nieznan
֒
a funkcj
֒
a,
a jej pochodnymi; gdy funkcja niewiadoma jest funkcj
֒
a jednej zmiennej, to m´owimy
o r´ownaniu r´o˙zniczkowym zwyczajnym, w przeciwnym wypadku o r´ownaniu
r´o˙zniczkowym cz
֒
astkowym. Zwi
֒
azek postaci
F (t, x(t), x
′
(t), . . . , x
(n)
(t)) = 0
nazywamy r´ownaniem r´o˙zniczkowym zwyczajnym n-tego rz
֒
edu, je´sli lewa strona
istotnie zale˙zy od x
(n)
. Nie musi oba zale˙ze´c od x i t. Przyk ladowo r´ownanie
x
′′′
+ t(x
′
)
30
− e
x sin t
= 0
jest r´ownaniem r´o˙zniczkowym rz
֒
edu trzeciego. Funkcja x mo˙ze by´c funkcj
֒
a ska-
larn
֒
a, albo wektorow
֒
a.
R´ownania r´o˙zniczkowe w zagadnieniach technicznych powstaj
֒
a na og´o l w wy-
niku stosowania nast
֒
epuj
֒
acych metod post
֒
epowania:
a) Przedstawiania praw fizyki w postaci matematyczno-analitycznej.
b) Przedstawiania zwi
֒
azk´ow geometrycznych w postaci analitycznej.
c) Rugowania parametr´ow z n-parametrowej rodziny funkcji i n r´owno´sci.
Ad a) Niech v : R
× R
3
⊃ [t
0
, T ]
× R
3
∋ (t, x) → v(t, x) ∈ R
3
b
֒
edzie zadanym
polem pr
֒
edko´sci. R´ownanie x
′
= v(t, x) opisuje ruchy cz
֒
astek unoszonych w polu
v. Je´sli dodatkowo przyj
֒
a´c warunek x(t
0
) = x
0
, to x(t) jest po lo˙zeniem w chwili
t tej cz
֒
astki, kt´ora w chwili t
0
znajdowa la si
֒
e w punkcie x
0
.
Ad b) Niech y = f (x). Wielko´s´c
ρ(A) =
(1 + (y
′
)
2
)
3
2
|y
′′
|
(A)
5
Rozdzia l 1. Wprowadzenie
nazywamy promienie krzywizny, a jej odwrotno´s´c
1
ρ
(A) krzywizn
֒
a w punkcie A.
R´ownanie r´o˙zniczkowe
|y
′′
|
(1 + (y
′
)
2
)
3
2
= a
R
∋ a ≥ 0
jest zadem r´ownaniem r´o˙zniczkowym, kt´orego rozwi
֒
azaniem s
֒
a krzywe o sta lej
krzywi´znie r´ownej a.
Ad c) Rozwa˙zmy rodzin
֒
e okr
֒
eg´ow
(x
− a)
2
+ (y
− b)
2
= R
2
,
(1.1)
gdzie a, b, R parametry. Za´o˙zmy, ˙ze y = y(x). R´o˙zniczkuj
֒
ac trzykrotnie zwi
֒
azek
(1.1) dostajemy
x
− a + (y − b)y
′
= 0
1 + (y
′
)
2
+ (y
− b)y
′′
= 0
3y
′
y
′′
+ (y
− b)y
′′′
= 0.
Ruguj
֒
ac z tych r´owna´
n wszystkie trzy parametry dostajemy r´ownanie r´o˙zniczkowe
rodziny okr
֒
eg´ow:
3y
′
(y
′′
)
2
− 1 + (y
′
)
2
y
′′′
= 0.
Bez nale˙zytej precyzji mo˙zemy przyj
֒
a´c w tej cwili, ˙ze r´ownaniem r´o˙zniczkowym
nazywamy r´ownanie postaci
F (t, x, x
′
, . . . , x
n
) = 0.
(1.2)
Je´sli funkcja ϕ : [a, b]
→ R klasy C
n
spe lnia to˙zsamo´sciowo r´owno´s´c
F (t, ϕ(t), ϕ
′
(t), . . . , ϕ
n
(t)) = 0 w [a, b],
to ϕ nazywamy ca lk
֒
a szczeg´oln
֒
a r´ownania r´o˙zniczkowego. Gdy ϕ jest funkcj
֒
a ele-
mentarn
֒
a, to m´owimy ˙ze (1.2) ma rozwi
֒
azanie efektywne. Na przyk lad r´ownanie
x
′′
(t) + ω
2
x(t) = 0
ma rozwi
֒
azania efektywne
ϕ
1
(t) = sin ωt i ϕ
2
(t) = cos ωt.
Z kolei r´ownanie Riccatiego
dx
dt
= a x
2
+ b t
n
a, b sta le,
n
∈ N
ma rozwi
֒
azanie efektywne (niestety) tylko dla pewnych n.
Je´sli rozwi
֒
azanie mo˙zna wyznaczy´c przez sko´
nczon
֒
a liczb
֒
e ca lkowaN, to m´owimy,
6
˙ze tak przedstawione rozwi
֒
azania s
֒
a rozwi
֒
azaniami przez kwadratur
֒
e. Na przyk lad
r´ownanie
x
′
=
sin t
t
ma rozwi
֒
azanie
x(t) =
Z
sin t
t
dt + C.
Niestety s
֒
a r´ownania, kt´ore nie s
֒
a rozwi
֒
azywalne przez kwadratur
֒
e. Przyk ladem
takiego r´ownania jest r´ownanie Bessela
t
2
x
′′
+ tx
′
+ t
2
− n
2
x = 0.
Mo˙zna dla niego poda´c rozwi
֒
azanie w postaci szereg´ow funkcyjnych. W szczeg´olno´sci
funkcje
I
0
(t) =
∞
X
k=0
(
−1)
k
(k!)
2
t
2
2k
,
Y
0
(t) = 2
∞
X
k=0
(
−1)
k
(k!)
2
t
2
2k
ln
t
2
+ C
−
k
X
ν=1
1
ν
!
,
gdzie C = 0.5772157 . . . jest sta l
֒
a Eulera, s
֒
a rozwi
֒
azaniami r´ownania Bessela dla
n = 0. Funkcje I
0
i Y
0
nosz
֒
a nazw
֒
e funkcji Bessela 1-go i 2-go rodzaju rz
֒
edu 0.
Nie ka˙zde r´ownanie r´o˙zniczkowe ma rozwi
֒
azanie. R´ownanie
1 +
dx
dt
2
= 0
nie ma rozwi
֒
aza´
n rzeczywistych, ma jednak rozwi
֒
azanie zespolone
x(t) = it.
R´ownanie
exp
dx
dt
= 0
w og´ole nie ma rozwi
֒
aza´
n, bo funkcja C
∋ z → e
z
∈ C nie ma zer. Z kolei
r´ownanie
x
′
= f (t, x),
gdzie prawa strona jest ci
֒
ag
֒
a ma niesko´
nczenie wiele rozwi
֒
aza´
n
Rozdzia l 2
Elementy analizy funkcjonalnej
Za l´o˙zmy, ˙ze X
6= ∅.
Definicja 1.
Funcj
֒
e ρ : X
×X → [0, ∞) nazywamy metryk
֒
a, wtedy i tylko wtedy,
gdy
1.
∀
x,y∈X
ρ(x, y) = 0
⇐⇒ x = y,
2.
∀
x,y∈X
ρ(x, y) = ρ(y, x),
3.
∀
x,y,z∈X
ρ(x, z)
≤ ρ(x, y) + ρ(y, z).
Definicja 2.
Je´sli X
6= ∅ i ρ : X × X → R metryka, to par
֒
e (X, ρ) nazywamy
przestrzeni
֒
a metryczn
֒
a.
Niech X b
֒
edzie przestrzeni
֒
a wektorow
֒
a nad cia lem K (K = R, lub K = C).
Definicja 3.
Funkcj
֒
e
k · k : X → [0, ∞) nazywamay norm
֒
a, wtedy i tylko wtedy,
gdy
1.
kxk = 0 ⇐⇒ x = 0,
2. ∀
α∈K
∀
x∈X
kαxk = |α|kxk,
3.
∀
x,y∈X
kx + yk ≤ kxk + kyk.
Definicja 4.
Par
֒
e (X,
k · k) nazywamy przestrzeni
֒
a unormowan
֒
a.
Uwaga 1.
Ka˙zda norma indukuje metryk
֒
e wed lug wzoru
ρ(x, y) :=
kx − yk,
tote˙z ka˙zda przestrze´
n unormowana jest przestrzeni
֒
a metryczn
֒
a.
9
Rozdzia l 2. Elementy analizy funkcjonalnej
Definicja 5.
Niech (X, ρ) - przestrze´
n metryczna. Ci
֒
ag
{x
n
}
n∈N
⊂ X nazywamy
ci
֒
agiem Cauchy’ego (ci
֒
agiem fundamentalnym) wtedy i tylko wtedy, gdy
∀
ε>0
∃
k∈N
∀
m>k
∀
n>k
ρ (x
m
, x
n
) < ε.
Definicja 6.
Niech (X, ρ) - przestrze´
n metryczna. M´owimy, ˙ze ci
֒
ag
{x
n
}
n∈N
⊂
X jest zbie˙zny do granicy g
∈ X wtedy i tylko wtedy, gdy ci
֒
ag liczbowy ρ (x
n
, g)
ma granic
֒
e r´own
֒
a 0, tj.
lim
n→∞
x
n
= g
⇐⇒ lim
n→∞
ρ (x
n
, g) = 0
⇐⇒ ∀
ε>0
∃
k∈N
∀
N
∋n>k
ρ (x
n
, g) < ε
Definicja 7.
M´owimy, ˙ze ci
֒
ag
{x
n
}
n∈N
⊂ X jest zbie˙zny w przestrzeni metrycz-
nej (X, ρ) wtedy i tylko wtedy, gdy istnieje g
∈ X, takie ˙ze lim
n→∞
x
n
= g.
Twierdzenie 1.
Ka˙zdy ci
֒
ag zbie˙zny w przestrzeni metrycznej (X, ρ) jest ci
֒
agiem
Cauchy’ego.
Twierdzenie odwrotne nie jest prawdziwe.
Przyk lad 1.
Ci
֒
ag
1
n
n∈N
jest zbie˙zny do zera w przestrzeni metrycznej (R, ρ
E
),
gdzie ρ
E
jest metryk
֒
a euklidesow
֒
a. Jest on zatem w my´sl poprzedniego twierdzenia
ci
֒
agiem Cauchy’ego. Niech X := (0, 1) i niech d b
֒
edzie restrykcj
֒
a metryki ρ
E
do
X
× X. Przestrze´n (X, d) jest przestrzeni
֒
a metryczn
֒
a, a rozwa˙zany ci
֒
ag w tej
przestrzeni nie jest zbie˙zny, gdy˙z 0
6∈ X.
Definicja 8.
Przestrze´
n metryczn
֒
a (X, ρ) nazywamy zupe ln
֒
a, wtedy i tylko wtedy,
gdy ka˙zdy ci
֒
ag Cauchy’ego
{x
n
}
n∈N
⊂ X jest zbie˙zny (do elementu przestrzeni
X).
Definicja 9.
Przestrze´
n unormowan
֒
a zupe ln
֒
a nazywamy przestrzeni
֒
a Banacha.
Twierdzenie 2.
(Banacha o odwzorowaniach zw
֒
e˙zaj
֒
acych)
Je´sli
- (X,
k · k) przestrze´n Banacha,
- T : X
→ X q-zw
֒
e˙zaj
֒
ace tzn.
∃
q∈[0,1)
∀
x,y∈X
kT (x) − T (y)k ≤ qkx − yk,
to
— T ma jedyny punkt sta ly tzn.
∃! x
⋆
∈ X :
T (x
⋆
) = x
⋆
.
— Ponadto, je´sli x
0
∈ X, x
n+1
:= T (x
n
), to
ρ (x
⋆
, x
p
)
≤
q
p
1
− q
ρ (x
1
, x
p
)
dla
p
∈ N.
Rozdzia l 3
Twierdzenia o istnieniu i
jednoznaczno´
sci
Twierdzenie 3.
Je´sli
1. t
0
∈ I = [a, b] ⊂ R,
x
0
∈ B = B (x
0
, R)
⊂ U ∈ topX,
f
∈ C(I × U, X),
2. funkcja f : I
× U ∋ (t, x) → f(t, x) ∈ X spe lnia warunek Lipschitza
wzgl
֒
edem drugiej zmiennej na zbiorze I
× B tzn.:
∃
L>0
∀
t∈I
∀
y,z∈B
kf(t, y) − f(t, z)k ≤ ky − zk,
3. rozwa˙zamy r´ownanie r´o˙zniczkowe postaci:
(RR)
x
′
(t) = f (t, x(t)) t
∈ I,
(WPC) x (t
0
) = x
0
,
to r´ownanie (RR) z zadanym warunkiem pocz
֒
atkowym Cauchy’ego (WPC) ma
dok ladnie jedno rozwi
֒
azanie x = x(t) na przedziale J = I
∩ [t
0
− r, t
0
+ r], gdzie
r :=
+
∞ gdy R = +∞ czyli B = X
R
M
gdy R < +
∞
i M := sup
{kf(t, y)k : t ∈ T, y ∈ B}.
Definicja 10.
Niech
(X, d), (Y, ρ) przestrzenie metryczne,
U
⊂ X,
11
Rozdzia l 3. Twierdzenia o istnieniu i jednoznaczno´sci
f : I
× U ∋ (t, x) → f(t, x) ∈ Y .
M´owimy, ˙ze f spe lnia lokalnie warunek Lipschitza wzgl
֒
edem zmiennej x, je˙zeli
∀
t
0
∈I
∀
x
0
∈U
∃
J∈top(t
0
)
∃
B=B(x
0
,R)
∃
L=L(J,B)
∀
t∈J
∀
y,z∈B∩U
ρ (f (t, y), f (t, z))
≤ L · d(y, z).
Twierdzenie 4.
Je˙zeli U
∈ topX, f = f(t, x) ∈ C(I × U, X), f spe lnie lokalnie
warunek Lipschitza wzgl
֒
edem zmiennej x, to dla ka˙zdego (t
0
, x
0
)
∈ I ×U r´ownanie
x
′
= f (t, x) z warunkiem pocz
֒
atkowym Cauchy’ego x (t
0
) = x
0
ma dok ladnie jedno
rozwi
֒
azanie okre´slone w pewnym otoczeniu punktu t
0
.
Twierdzenie 5.
(Zasada identyczno´sci) Przyjmijmy za lo˙zenia poprzedniego twier-
dzenia. Niech P przedzia l, P
⊂ I. Niech x = x(t), y = y(t) b
֒
ed
֒
a dwoma
rozwi
֒
azaniami tego samego r´ownania r´o˙zniczkowego x
′
= f (t, x) okre´slonymi na
P i spe lniaj
֒
acymi warunki pocz
֒
atkowe Cauchy’ego x (t
1
) = x
0
, y (t
2
) = y
0
. Je´sli
istnieje taki punkt p
∈ P , w kt´orym x(p) = y(p), to x(t) = y(t) dla t ∈ P .
Twierdzenie 6.
Niech Y = X
n
, U
∈ topY , f ∈ C(I × U, X) i niech f = f(t, y)
spe lnia lokalnie warunek Lipschitza wzgl
֒
edem zmiennej y. Wtedy dla ka˙zdego
t
0
∈ I, dla ka˙zdego x
0
= (x
01
, . . . , x
0n
)
∈ U r´ownanie r´o˙zniczkowe
x
(n)
= f t, x, x
′
, . . . , x
(n−1)
z warunkiem pocz
֒
atkowym Cauchy’ego
x
(j)
(t
0
) = x
0j
j = 0, 1, . . . , n
− 1
ma dok ladnie jedno rozwi
֒
azanie x = x(t) w pewnym otoczeniu punktu t
0
.
Dow´
od
. R´ownanie sprowadzamy do uk ladu r´owna´
n. Niech y
1
:= x oraz
y
′
1
= y
2
=: f
1
(t, y
1
, . . . , y
n
)
y
′
2
= y
3
=: f
2
(t, y
1
, . . . , y
n
)
. . .
. . .
y
′
n−1
= y
n
=: f
n−1
(t, y
1
, . . . , y
n
)
y
′
n
= f (t, y
1
, . . . , y
n
) =: f
n
(t, y
1
, . . . , y
n
)
Uk lad ten mo˙zna zapisa´c w postaci
Y
′
=
F(t, Y),
gdzie
Y = (y
1
, . . . , y
n
)
T
,
F = (f
1
, . . . , f
n
)
T
.
c.k.d
Definicja 11.
Rozwi
֒
azanie okre´slone na ca lym przedziale I okre´slono´sci r´ownania
r´o˙zniczkowego nazywamy rozwi
֒
azaniem globalnym tego r´ownania.
Twierdzenie 7.
(o rozwi
֒
azaniu globalnym) Niech t
0
∈ I = |a, b| ⊂ R i niech
f
∈ C(I × X, X) i niech dane b
֒
edzie r´ownanie
x
′
= f (t, x),
t
∈ I
z warunkiem pocz
֒
atkowym Cauchy’ego
12
x (t
0
) = x
0
.
Je´sli
∀
J=[a
′
,b
′
]⊂I
∃
L=L(J)>0
∀
t∈J
∀
x,y∈X
kf(t, x) − f(t, y)k ≤ L kx − yk ,
to powy˙zsze r´ownanie r´o˙zniczkowe z dowolnie zadanym warunkiem pocz
֒
atkowym
Cauchy’ego ma dok ladnie jedno rozwi
֒
azanie globalne tj. okre´slone na przedziale
I.
Podobne twierdzenie ma miejsce dla uk lad´ow r´owna´
n.
Przyk lad 2.
R´ownanie x
′
= x
2
nie spe lnia za lo˙ze´
n powy˙zszego twierdzenie. Ca lka
og´olna tego r´ownania jest okre´slona wzorem x(t) =
−
1
t+C
(C
∈ R) i nie jest
okre´slona na X = R.
Rozdzia l 4
Proste typy r´
owna´
n r´
o ˙zniczkowych
skalarnych
4.1. R´
ownanie r´
o ˙zniczkowe o zmiennych rozdzielonych
R´ownanie r´o˙zniczkowe postaci
x
′
(t) =
f (t)
g(x)
,
(4.1)
gdzie f
∈ C(I, R), g ∈ C(J, R), x ∈ C
1
(I, J), I, J przedzia ly, t
∈ I, g(x) 6= 0 dla
x
∈ J nazywamy r´ownaniem o zmiennych rozdzielonych. R´ownanie to mo˙zemy
zapisa´c w postaci
g(x)x
′
(t) = f (t).
Niech G = G(x) oraz F = F (t) b
֒
ed
֒
a dowolnymi funkcjami pierwotnymi odpo-
wiednio funkcji g = g(x) i f = f (t). W´owczas r´ownanie (4.1) mo˙zna przepisa´c w
postaci
d
dt
(G
◦ x)(t) =
d
dt
F (t),
czyli
d
dt
[G(x)
− F (t)] = 0, x = x(t), t ∈ I.
Poniewa˙z I przedzia l, to r´ownanie to na podstawie twierdzenia Lagrange’a jest
r´ownowa˙zne r´ownaniu
G(x)
− F (t) = C, x = x(t), t ∈ I, C ∈ R,
15
Rozdzia l 4. Proste typy r´
owna´
n r´
o˙zniczkowych skalarnych
kt´ore mo˙zemy zapisa´c w postaci
Z
g(x)dx =
Z
f (t)dt,
x = x(t).
(4.2)
4.2. R´
ownanie jednorodne
R´ownaniem r´o˙zniczkowym jednorodnym nazywamy r´ownanie postaci
x
′
= f
x
t
,
(4.3)
gdzie t
∈ I, x = x(t), f ∈ C(J, R), I, J - przedzia ly. Podstawienie
x(t) = ty(t)
sprowadza r´ownanie (4.3) do r´ownania r´o˙zniczkowego
y
′
=
f (y)
− y
t
o zmiennych rozdzielonych.
Dodatkowo nale˙zy sprawdzi´c, czy rozwi
֒
azaniem
r´ownania (4.3) jest funkcja x(t) := y
0
t, gdzie y
0
, jest rozwi
֒
azaniem r´ownania
f (y
0
)
− y
0
= 0.
R´ownanie
dx
dt
= f
a
1
t + b
1
x + c
1
a
2
t + b
2
x + c
2
,
(4.4)
gdzie f jest funkcj
֒
a ci
֒
ag l
֒
a oraz a
1
b
2
− a
2
b
1
6= 0 mo˙zna przez stosown
֒
a zmian
֒
e
zmiennych sprowadzi´c do r´ownania jednorodnego. Je´sli bowiem wektor (¯
t, ¯
x) jest
rozwi
֒
azaniem uk ladu r´owna´
n
a
1
b
1
a
2
b
2
t
x
=
−c
1
−c
2
to zmiana zmiennych
t = ¯
t + ξ,
x = ¯
x + η
przy kt´orej
dη
dξ
=
d(x−¯
x)
dt
dt
dξ
=
dx
dt
sprowadza r´ownanie (4.4) do r´ownania jednorod-
nego
dη
dξ
= f
a
1
ξ + b
1
η
a
2
ξ + b
2
η
= f
a
1
+ b
1
η
ξ
a
2
+ b
2
η
ξ
!
=: g
η
ξ
.
Gdy a
1
b
2
− a
2
b
1
= 0, to istnieje takie λ
∈ R, ˙ze a
2
t + b
2
x = λ (a
1
t + b
1
x) lub
a
1
t + b
1
x = λ (a
2
t + b
2
x). R´ownanie (4.4) przekszta lca si
֒
e w r´ownanie postaci
x
′
= ¯
f (a
1
t + b
1
x)
lub x
′
= ¯
f (a
2
t + b
2
x) .
Podstawienie odpowiednio
u(t) = a
1
t + b
1
x(t) lub u(t) = a
2
t + b
2
x(t)
sprowadza je do r´ownania o zmiennych rozdzielonych.
16
4.3. R´
ownanie r´
o˙zniczkowe zupe lne
4.3. R´
ownanie r´
o ˙zniczkowe zupe lne
Niech D
⊂ R
2
b
֒
edzie obszarem tj. zbiorem otwartym i sp´ojnym. Niech
P, Q
∈ C(D, R) oaz Q(t, x) 6= 0 dla (t, x) ∈ D.
Definicja 12.
R´ownanie r´o˙zniczkowe
x
′
=
−
P (t, x)
Q(t, x)
(4.5)
czyli
P (t, x)dt + Q(t, x)dx = 0
(4.6)
nazywamy zupe lnym wtedy i tylko wtedy, gdy istnieje taka funkcja U
∈ C
1
(D, R),
˙ze
d
(t,x)
U = P (t, x)dt + Q(t, x)dx dla (t, x)
∈ D.
(4.7)
Poniewa˙z zbi´or D jest obszarem, zatem je´sli (4.6) jest r´ownaniem r´o˙zniczkowym
zupe lnym, to ca lka og´olna tego r´ownania ma posta´c
U (t, x) = C,
C
∈ R.
Z twierdzenia Poincar`e’go wynika nast
֒
epuj
֒
ace
Twierdzenie 8.
Je´sli D jest obszarem ´sci
֒
agalnym w R
2
, P, Q
∈ C(D, R) oraz
∂P
∂x
=
∂Q
∂t
w D, to (4.6) jest r´ownaniem r´o˙zniczkowym zupe lnym,
przy czym:
Definicja 13.
Obszar D nazywamy ´sci
֒
agalnym w R
2
wtedy i tylko wtedy, gdy
istniej
֒
a obszar obszar gwia´zdzisty G
⊂ R
2
oraz dyffeomorfizm h : G
→ D (tzn. h
bijekcja, H, h
−1
klasy C
1
).
Definicja 14.
Zbi´or G
⊂ R
2
nazywamy zbiorem gwia´zdzistym wtedy i tylko
wtedy, gdy
∃
x
0
∈G
∀
x∈G
[x
0
, x]
⊂ G,
Wiedz
֒
ac, ˙ze (4.6) zupe lne z warunku (4.7) mamy:
∂U
∂t
= P,
∂U
∂x
= Q.
Ca lkuj
֒
ac pierwszy z tych zwi
֒
azk´ow wzgl
֒
edem zmiennej t dostajemy:
U (t, x) =
Z
P (t, x)dt + C(x) (t, x)
∈ D.
Z kolei
Q(t, x) =
∂U
∂x
(t, x) =
Z
∂P (t, x)
∂x
dt + C
′
(x),
17
Rozdzia l 4. Proste typy r´
owna´
n r´
o˙zniczkowych skalarnych
sk
֒
ad
C
′
(x) = Q(t, x)
−
Z
∂P
∂x
(t, x)dt.
i w konsekwencji
C(x) =
Z
Q(t, x)dx
−
Z Z
∂P
∂x
(t, x)dt
dx.
Ostatecznie
U (t, x) =
Z
P (t, x)dt +
Z
Q(t, x)dx
−
Z Z
∂P
∂x
(t, x)dt
dx,
tak wi
֒
ec rozwi
֒
azanie og´olne r´ownania r´o˙zniczkowego zupeLnego (4.6) wyra˙za si
֒
e
wzorem:
Z
P (t, x)dt +
Z
Q(t, x)dx
−
Z Z
∂P
∂x
(t, x)dt
dx = C,
C
∈ R.
(4.8)
4.3.1. Czynnik ca lkuj
֒
acy
Je˙zeli r´ownanie (4.6) nie spe lnia warunku
∂P
∂x
=
∂Q
∂t
w zadanym obszarze
´sci
֒
agalnym D, to szukamy takiej funkcji µ = µ(t, x)
∈ C
1
(D, R), aby
∂(µP )
∂x
=
∂(µQ)
∂t
(t, x)
∈ D
(4.9)
Definicja 15.
Funkcj
֒
e µ
∈ C
1
(D, R), dla kt´orej zachodzi warunek (4.9) nazy-
wamy czynnikiem ca lkuj
֒
acym r´ownania (4.6).
Twierdzenie 9.
Je´sli funkcje P, Q
∈ C
1
(D, R) i D obszar ´sci
֒
agalny, to istnieje
µ
∈ C
1
(D, R) czynnik ca lkuj
֒
acy r´ownania (4.6).
Efektywne wyznaczenie czynnika ca lkuj
֒
acego jest mo˙zliwe zawsze, gdy zale˙zy on
od jednej zmiennej oraz w sytuacji, gdzy µ = µ(ω(t, x)), gdzie ω(t, x) jest znan
֒
a
funkcj
֒
a klasy C
1
(D, R). W pozosta lych przypadkach jest to zagadnienie trudne
cz
֒
esto niemo˙zliwe do zrealizowania.
Za´o˙zmy zatem, ˙ze istnieje czynnik ca lkuj
֒
acy r´ownania (4.6) postaci µ =
µ(ω(t, x)). Warunek
∂(µP )
∂x
=
∂(µQ)
∂t
(t, x)
∈ D
jest r´ownowa˙zny warunkowi
µ
′
∂ω
∂x
P + µ
∂P
∂x
= µ
′
∂ω
∂t
Q + µ
∂Q
∂t
,
kt´ory mo˙zna zapisa´c w postaci
µ
′
µ
=
∂Q
∂t
−
∂P
∂x
∂ω
∂x
P
−
∂ω
∂t
Q
.
(4.10)
18
4.4. R´
ownanie Clairauta
Poniewa˙z lewa strona, z za lo˙zenia, zale˙zy od ω(t, x), zatem warunkiem istnienia
czynnika ca lkuj
֒
acego postaci µ = µ(ω(t, x)) jest aby prawa strona r´ownania (4.10)
by la zale˙zna od ω(t, x). Wtedy te˙z dostajemy wz´or:
ln
|µ(ω)| =
Z
∂Q
∂t
−
∂P
∂x
∂ω
∂x
P
−
∂ω
∂t
Q
(ω)
!
dω =: χ(ω)
z kt´orego wynika, ˙ze ka˙zda z funkcji
e
µ(t, x) := µ(ω(t, x)) = Ce
χ(ω(t,x))
(C
∈ R \ {0})
(4.11)
jest szukanym czynnikiem ca lkuj
֒
acym.
Poszukuj
֒
ac czynnika ca lkuj
֒
acego nale˙zy rozpocz
֒
ac od najprostszych przy-
padk´ow tj. ω(t, x) = t lub ω(t, x) = x, potem rozwa˙zy´c kolejno ω(t, x) = t + x,
ω(t, x) = t
− x, ω(t, x) = tx, ω(t, x) =
t
x
. Gdy nie przyniesie to rezultatu szanse
na znalezienie czynnika ca lkuj
֒
acego s
֒
a znikome.
Przyk lad 3.
Istnieje czynnik ca lkuj
֒
acy µ = µ(t) r´ownania (t + t
2
+ x
2
) dt +
xdx = 0, gdy˙z
µ
′
(t)
µ(t)
= 2. Rozwi
֒
azuj
֒
ac ostatnie r´ownanie dostajemy
d
dt
ln
|µ(t)| = 2
i w konsekwencji µ(t) = Ce
2t
(C
∈ R\{0}) jest szukanym czynnikiem ca lkuj
֒
acym.
4.4. R´
ownanie Clairauta
Definicja 16.
R´ownaniem Clairauta nazywamy r´ownanie r´o˙zniczkowe
x
− tx
′
− f (x
′
) = 0,
(4.12)
gdzie t
∈ I, I - przedzia l, x ∈ C
2
(I, J), J - przedzia l, f
∈ C
1
(J, R) i funkcja f
nie jest postaci f (τ ) = Aτ + B.
R´o˙zniczkuj
֒
ac (4.12) stronami dostajemy:
x
′
− x
′
− tx
′′
− f
′
(x
′
) x
′′
= 0
czyli
x
′′
(t + f
′
(x
′
)) = 0.
Je´sli istnieje x = x(t) rozwi
֒
azanie r´ownania (4.12) klasy C
2
(I, R), to
x
′′
= 0 lub t + f
′
(x
′
) = 0.
Je´sli x
′′
(t) = 0, to x
′
(t) = C, x(t) = Ct + b. Wstawiaj
֒
ac funkcj
֒
e x(t) = Ct + b
do r´ownania (4.12) dostajemy b = f (C). Tak wi
֒
ec ka˙zda prosta
x(t) = Ct + f (C),
C
∈ J
(4.13)
jest rozwi
֒
azaniem (4.12).
Rozdzia l 4. Proste typy r´
owna´
n r´
o˙zniczkowych skalarnych
W sytuacji t + f
′
(x
′
) = 0, traktujemy pochodn
֒
a x
′
jak parametr i oznaczamy
go symbolem p. Tak wi
֒
ec t =
−f
′
(p). R´ownanie (4.12) mo˙zemy przepisa´c w
postaci x = tp + f (p) =
−f
′
(p)p + f (p). R´ownanie parametryczne
t
=
−f
′
(p)
x = f (p)
− pf
′
(p)
(4.14)
jest r´ownaniem obwiedni rodziny prostych (4.13).
20
Rozdzia l 5
Liniowe r´
ownania r´
o ˙zniczkowe
5.1. R´
ownania i uk lady r´
owna´
n r´
o ˙zniczkowych liniowych
Niech (X,
k · k) przestrze´n Banacha, I = |a, b| ⊂ R - dowolny przedzia l,
L(X, X) :=
{T : X → X :
T
operator liniowy i ci
֒
ag ly
}. Niech
A : I
∋ t → A(t) ∈ L(X, X) ci
֒
ag le,
g
∈ C(I, X),
x = x(
·) ∈ C
1
(I, X).
Definicja 17.
R´ownaniem r´o˙zniczkowym liniowym jednorodnym rz
֒
edu pierw-
szego (RRLJ) nazywamy r´ownanie postaci
x
′
(t) = A(t) (x(t)) ,
t
∈ I,
(5.1)
kr´otko x
′
= A(t)x, x = x(t), t
∈ I.
Definicja 18.
R´ownaniem r´o˙zniczkowym liniowym niejednorodnym rz
֒
edu pierw-
szego (RRLN) nazywamy r´ownanie postaci
x
′
(t) = A(t) (x(t)) + g(t),
t
∈ I,
(5.2)
kr´otko x
′
= A(t)x + g(t), x = x(t), t
∈ I.
Definicja 19.
W sytuacji X = R
n
(RRLJ), (RRLN) nazywamy uk ladem r´owna´
n
r´o˙zniczkowych liniowych.
Definicja 20.
R´ownanie r´ozniczkowe
x
(n)
= A(t) x, x
′
, . . . , x
(n−1)
+ g(t),
(5.3)
21
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
gdzie x = x(t), t
∈ I, I - przedzia l, A ∈ C (I, L (X
n
, X)), g
∈ C(I, X), X
- przestrze´
n Banacha, nazywamy r´ownaniem r´o˙zniczkowym liniowym rz
֒
edu n -
tego. Je´sli g = 0, to r´ownanie (5.3) nazywamy r´ownaniem jednorodnym, w prze-
ciwnym wypadku niejednorodnym.
Jak wiadomo z wcze´sniejszych rozwa˙za´
n, r´ownanie to mo˙zna sprowadzi´c do
r´ownania rz
֒
edu pierwszego w przestrzeni Banacha X
n
.
Twierdzenie 10.
(Twierdzenie o istnieniu rozwi
֒
azania globalnego) Standardowe
r´ownanie r´o˙zniczkowe liniowe (5.2) ma zawsze rozwi
֒
azanie globalne przy dowol-
nym warunku pocz
֒
atkowym Cauchy’ego.
Twierdzenie 11.
Zbi´or rozwi
֒
aza´
n r´ownania r´o˙zniczkowego liniowego jednorod-
nego (5.1) (ca lka og´olna) jest przestrzeni
֒
a liniow
֒
a.
Dow´
od
. Wystarczy pokaza´c, ˙ze je´sli funkcje x i y s
֒
a rozwi
֒
azaniami (5.1) to ich
dowolna kombinacja liniowa tak˙ze. Niech α, β
∈ R. Mamy
(αx + βy)
′
= αx
′
+ βy
′
= αA(t)x + βA(t)y =
= A(t)(αx) + A(t)(βy) = A(t)(αx + βy)
c.k.d
Twierdzenie 12.
Rozwi
֒
azanie og´olne r´ownania r´o˙zniczkowego liniowego niejed-
norodnego (5.2) jest sum
֒
a rozwi
֒
azania szczeg´olnego (5.2) i rozwi
֒
azania og´olnego
r´ownania r´o˙zniczkowego liniowego jednorodnego(5.1), a dok ladniej:
Ka˙zde rozwi
֒
azanie (5.2) jest sum
֒
a pewnego ustalonego rozwi
֒
azania (5.2) i
pewnego rozwi
֒
azania (5.1).
Dow´
od
. Niech
R :=
x
∈ C
1
(I, X) :
x
′
= A(t)x + g(t)
,
Y
:=
y
∈ C
1
(I, X) :
y
′
= A(t)y
.
Ustalmy e
x
∈ R i zdefiniujmy
Z :=
z
∈ C
1
(I, X) :
z = e
x + y, y
∈ Y
= e
x + Y.
Mamy pokaza´c, ˙ze R = Z.
Udowodnimy najpierw, ˙ze Z
⊂ R.
We´zmy z
∈ Z. Z definicji zbioru Z wynika, ˙ze istnieje y ∈ Y , ˙ze z = ex + y.
Poniewa˙z
z
′
= e
x
′
+ y
′
= (A(t)e
x + g(t)) + A(t)y = A(t) (e
x + y) + g(t) = A(t)z + g(t)
zatem z
∈ R.
22
5.2. Skalarne r´
ownanie liniowe rz
֒
edu pierwszego
Teraz udowodnimy, ˙ze R
⊂ Z.
We´zmy x
∈ R. Wektor x mo˙zemy zapisa´c w postaci x = ex+(x − ex). Zdefiniujmy
y := x
− ex. Zauwa˙zmy, ˙ze
y
′
= (x
− ex)
′
= x
′
= e
x
′
= (A(t)x + g(t))
− (A(t)ex + g(t)) =
= A(t)x
− A(t)ex = A(t) (x − ex) = A(t)y,
co oznacza, ˙ze y
∈ Y . W takim razie x ∈ Z.
c.k.d
5.2. Skalarne r´
ownanie liniowe rz
֒
edu pierwszego
Skalarne r´ownanie liniowe rz
֒
edu pierwszego
x
′
+ f (t)x = 0,
(5.4)
gdzie x = x(t), t
∈ I, I - przedzia l, f ∈ C(I, R), jest r´ownaniem o zmiennych
rozdzielonych. Ca lk
֒
a og´oln
֒
a tego r´ownania jest rodzina funkcji
x(t) = Ce
−
R
f (t)dt
C = const
∈ R, t ∈ I.
Ca lk
֒
e szczeg´oln
֒
a r´ownania niejednorodnego
x
′
+ f (t)x = g(t) t
∈ I,
(5.5)
mo˙zemy znale´z´c metod
֒
a uzmienniania sta lej. Przypu´s´cmy bowiem, ˙ze istnieje
rozwi
֒
azanie r´ownania (5.5) postaci
x(t) = C(t)e
−
R
f (t)dt
= C(t)e
−F (t)
,
gdzie F (t) :=
R
f (t)dt. Je´sli funkcja ta jest rozwi
֒
azaniem r´ownania (5.5), to
g(t) = x
′
+f (t)x = C
′
(t)e
−F (t)
+C(t) (
−F
′
(t)) e
−F (t)
+f (t)C(t)e
−F (t)
= C
′
(t)e
−F (t)
,
sk
֒
ad
C
′
(t) =
g(t)
e
−F (t)
= g(t)e
F (t)
.
Rozwi
֒
azaniem tego r´ownania jest funkcja
C(t) =
Z
g(t)e
F (t)
dt,
t
∈ I.
Tak wi
֒
ec ca lk
֒
a szczeg´oln
֒
a r´ownania (5.5) jest funkcja
x(t) =
Z
g(t)e
R
f (t)dt
dt
e
−
R
f (t)dt
dt
23
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
5.3. R´
ownanie Bernoulliego
R´ownaniem r´o˙zniczkowym Bernoulliego nazywamy r´ownanie postaci
x
′
+ f (t)x = g(t)x
p
,
p = const
∈ R \ {1},
(5.6)
przy czym w stosunku do funkcji f i g przyjmujemy takie same za lo˙zenia jak w
przypadku r´ownania liniowego. Przez zmian
֒
e zmiennych
y(t) := x
1−p
(t)
r´ownanie to mo˙zna sprowadzi´c do r´ownanie r´o˙zniczkowego liniowego. Zauwa˙zmy
bowiem, ˙ze skoro y
′
= (1
− p)x
−p
x
′
, to obustronnie mno˙z
֒
ac r´ownanie (5.6) przez
(1
− p)x
−p
dostajemy
(1
− p)x
−p
x
′
+ (1
− p)f(t)x
1−p
= (1
− p)g(t),
czyli r´ownanie r´o˙zniczkowe liniowe niejednorodne
y
′
+ (1
− p)f(t)y = (1 − p)g(t).
5.4. R´
ownanie Riccatiego
R´ownaniem r´o˙zniczkowym Riccatiego nazywamy r´ownanie postaci
x
′
= a(t)x
2
+ b(t)x + c(t),
(5.7)
gdzie a, b, c : I
→ R ci
֒
ag le, I - przedzia l otwarty.
Z poprzednich twierdze´
n latwo pokaza´c, ˙ze ka˙zdy punkt zbioru I
×R jest punk-
tem globalnej jednoznaczno´sci. Gdy a(t) = 0, to r´ownanie (5.7) jest r´ownaniem
r´o˙zniczkowym liniowym, a gdy c(t) = 0 r´ownaniem Bernoulliego.
Specjalnym r´ownaniem Riccatiego nazywamy szczeg´olny przypadek r´ownanoa
(5.7) a mianowicie
x
′
= c
1
x
2
+ c
2
t
n
c
1
, c
2
∈ R.
Nawet dla tego ostatniego r´ownania mo˙zna poda´c efektywne metody dla pewnych
warto´sci wyk ladnika n. W og´olnym przypadku zachodzi natomiast nast
֒
epuj
֒
ace:
Twierdzenie 13.
Niech I = (α, β)
⊂ R. Je´sli ϕ jest ca lk
֒
a szczeg´oln
֒
a r´ownania
(5.7) okre´slon
֒
a na I, to dla ka˙zdego rozwi
֒
azania x tego r´ownania okre´slonego w
przedziale
△ ⊂ I funkcja okre´slona wzorem:
y(t) := x(t)
− ϕ(t) (t ∈ △),
jest rozwi
֒
azaniem r´ownania Bernoulliego
y
′
= [b(t) + 2a(t)ϕ(t)] y + a(t)y
2
(5.8)
24
5.5. R´
ownanie Lagrange’a
i na odwr´ot, dla ka˙zdego rozwi
֒
azania y r´ownania (5.8) okre´slonego w
△ funkcja
x zdefiniowana wzorem:
x(t) = ϕ(t) + y(t) (t
∈ △)
jest rozwi
֒
azaniem r´ownania (5.7).
Dow´
od
. Niech ϕ i x b
֒
ed
֒
a dwoma rozwi
֒
azaniami r´ownania (5.7), czyli
ϕ
′
= a(t)ϕ
2
+ b(t)ϕ + c(t),
x
′
= a(t)x
2
+ b(t)x + c(t).
W´owczas
y
′
= x
′
− ϕ
′
= a(t)x
2
+ b(t)x + c(t)
− a(t)ϕ
2
+ b(t)ϕ + c(t)
=
= a(t) x
2
− ϕ
2
+ b(t) (x
− ϕ) = a(t) (x + ϕ) (x − ϕ) + b(t) (x − ϕ) =
= (a(t) (x + ϕ) + b(t)) (x
− ϕ) = (b(t) + a(t)x + a(t)ϕ) (x − ϕ) =
= (b(t) + 2a(t)ϕ + a(t)x
− a(t)ϕ) (x − ϕ) =
= (b(t) + 2a(t)ϕ + a(t) (x
− ϕ)) (x − ϕ) = (b(t) + 2a(t)ϕ) y + a(t)y
2
.
Tak wi
֒
ec
y
′
= (b(t) + 2a(t)ϕ) y + a(t)y
2
.
c.k.d
Przyk lad 4.
Rozwa˙zmy r´ownanie Riccatiego
x
′
− 2tx + x
2
= 5
− t
2
,
kt´orego ca lk
֒
a szczeg´oln
֒
a jest funkcja ϕ(t) = t + 2. Przepisuj
֒
ac to r´ownanie w
postaci x
′
= (
−1)x
2
+ (2t)x + (5
− t
2
), widzimy, ˙ze a(t) =
−1, b(t) = 2t, c(t) =
5
− t
2
. Skojarzone r´ownanie Bernoulliego przybiera wi
֒
ec posta´c
y
′
= [2t + 2(
−1)(t + 2)] y + (−1)y
2
=
−4y − y
2
.
Jego rozwi
֒
azaniem og´olnym jest rodzina funkcji
y(t) = Ce
4t
(C
∈ R), tak wi
֒
ec
rozwi
֒
azaniem r´ownania wyj´sciowego jest rodzina funkcji
x(t) = Ce
4t
+t+2(C
∈
R
).
5.5. R´
ownanie Lagrange’a
R´ownaniem Lagrange’a nazywamy r´ownanie postaci:
x = a (x
′
) t + f (x
′
) .
(5.9)
Zak ladamy, ˙ze funkcje a, f
∈ C
1
(J, R), x
∈ C
2
(I, J), I, J przedzia ly. Je´sli funkcja
a jest funkcj
֒
a identyczno´sciow
֒
a, to r´ownanie Lagrange’a jest r´ownaniem Cla-
irauta. Przyjmijmy zatem dalej, ˙ze a(p)
6= p dla wszystkich p ∈ J. R´o˙zniczkuj
֒
ac
25
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
r´ownanie (5.9) stronami i podstawiaj
֒
ac za pocchodn
֒
a x
′
now
֒
a funkcj
֒
e p = p(t)
mo˙zemy to r´ownanie przekszta lci´c do postaci:
x
′
= a
′
(x
′
) x
′′
t + a (x
′
) + f
′
(x
′
) x
′′
p = a
′
(p) p
′
t + a (p) + f
′
(p) p
′
p = (a
′
(p)t + f
′
(p))
dp
dt
+ a(p),
dp
dt
=
p
− a(p)
a
′
(p)t + f
′
(p)
.
Zamieniaj
֒
ac role zmiennych p i t mamy
dt
dp
=
a
′
(p)t + f
′
(p)
p
− a(p)
=
a
′
(p)
p
− a(p)
t +
f
′
(p)
p
− a(p)
,
czyli r´ownanie r´o˙zniczkowe niejednorodne
dt
dp
+
a
′
(p)
a(p)
− p
t =
f
′
(p)
p
− a(p)
,
z niewiadom
֒
a funkcj
֒
a t = t(p). Po wyznaczeniu tego rozwi
֒
azania wstawiamy
je do wyj´sciowego r´ownania (5.9), w kt´orym w miejsce pochodnej x
′
wstawiamy
parametr p. Ostatecznie
t
= t(p)
x = a (p) t(p) + f (p) .
(5.10)
jest rozwi
֒
azaniem r´ownania (5.9) w postaci parametrycznej.
5.6. Skalarne r´
ownanie r´
o ˙zniczkowe liniowe
n-tego rz
֒
edu
Definicja 21.
Skalarnym r´ownaniem r´o˙zniczkowym jednorodnym n-tego rz
֒
edu
(SRRLJ) nazywamy r´ownanie
x
(n)
+ a
n−1
(t)x
(n−1)
+ . . . + a
1
(t)x
′
+ a
0
(t)x = 0,
(5.11)
w kt´orym a
j
(t)
∈ C(I, R), (j = 0, 1, . . . , n − 1), I - przedzia l.
Niech
L(t) :=
d
n
dt
n
+ a
n−1
(t)
d
n−1
dt
n−1
+ . . . + a
1
(t)
d
dt
+ a
0
(t),
t
∈ I,
w´owczas r´ownanie (5.11) mo˙zna zapisa´c w zwi
֒
ez lej postaci
L(t)x = 0,
t
∈ I.
(5.12)
26
5.6. Skalarne r´
ownanie r´
o˙zniczkowe liniowe
n
-tego rz
֒
edu
Definicja 22.
Wro´
nskianem funkcji x
1
, . . . , x
n
∈ C
n−1
(I, R) nazywamy funkcj
֒
e
W (x
1
, . . . , x
n
) (t) := det
x
(k−1)
j
(t)
k = 1, . . . , n
j = 1, . . . , n
!
(5.13)
Twierdzenie 14.
a) Je´sli wro´
nskian W (x
1
, . . . , x
n
) (t
0
)
6= 0 dla pewnego t
0
∈ I,
to funkcke x
1
, . . . , x
n
s
֒
a liniowo niezale˙zne.
b) Niech x
1
, . . . , x
n
b
֒
ed
֒
a rozwi
֒
azaniami r´ownania (5.11). Je´sli x
1
, . . . , x
n
s
֒
a
liniowo niezale˙zne, to ich wro´
nskian W (x
1
, . . . , x
n
) (t))
6= 0 dla ka˙zdego t ∈ I.
Dow´
od
Kolejno udowodnimy obie cz
֒
e´sci twierdzenia.
ad a) (nie wprost)
Przyjmijmy, ˙ze x
1
, . . . , x
n
∈ C
n−1
(I, R) s
֒
a liniowo zale˙zne. Zatem istniej
֒
a takie
sta le C
1
, . . . , C
n
∈ R, ˙ze
P
n
j=1
C
2
j
6= 0 oraz
n
X
j=1
C
j
x
j
(t) = 0 dla t
∈ I.
(5.14)
R´o˙zniczkuj
֒
ac t
֒
e r´owno´s´c sukcesywnie wzgl
֒
edem zmiennej t dostajemy zwi
֒
azek
n
X
j=1
C
j
x
(k−1)
j
= 0 k = 1, . . . , n, t
∈ I.
Poniewa˙z W (x
1
, . . . , x
n
) (t
0
)
6= 0, zatem uk lad (5.14) ma tylko rozwi
֒
azanie zerowe
C
1
= C
2
= . . . = C
n
= 0 wbrew za lo˙zeniu.
ad b) (nie wprost) Przypu´s´cmy, ˙ze istnieje taki punkt t
0
∈ I : W (x
1
, . . . , x
n
) (t
0
) =
0.
Po l´o˙zmy
a
j
k
:= x
(k−1)
j
(t
0
) ,
j, k
∈ {1, . . . , n} i zdefiniujmy macierz A := a
j
k
.
Niech wektor C = (C
1
, . . . , C
n
)
T
b
֒
edzie niezerowym rozwi
֒
azaniem uk ladu
AC = 0.
Takie rozwi
֒
azanie istnieje, gdy˙z
det A = det a
j
k
= W (x
1
, . . . , x
n
) (t
0
) = 0.
We´zmy
x = x(t) :=
n
X
j=1
C
j
x
j
(t).
Funkcja ta jest rozwi
֒
azaniem r´ownania (5.11) bo jest kombinacj
֒
a liniow
֒
a rozwi
֒
aza´
n
x
j
(j = 1, . . . , n). Zauwa˙zmy, ˙ze x(t) spe lnia warunek pocz
֒
atkowy Cauchy’ego:
x
(k−1)
(t
0
) =
n
X
j=1
C
j
x
k−1
j
(t
0
) = 0.
27
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
Z drugiej strony funkcja sta la r´owna zero te˙z spe lnoa powy˙zszy warunek pocz
֒
atkowy
i jest rozwi
֒
azaniem r´ownania (5.11). Wobec jedyno´sci rozwi
֒
azania problemu
pocz
֒
atkowego dla r´ownania (5.11) i wobec liniowej niezale˙zno´sci x
1
, . . . , x
n
mamy
C
1
= C
2
= . . . = C
n
= 0 co przeczy za lo˙zeniu.
c.k.d
Wniosek 1.
Je˙zeli x
1
, . . . , x
n
s
֒
a rozwi
֒
azaniami r´ownania (5.11), to
∀
t∈I
W (x
1
, . . . , x
n
) (t) = 0,
lub
∀
t∈I
W (x
1
, . . . , x
n
) (t)
6= 0.
Definicja 23.
Zbi´or
{x
1
, . . . , x
n
} liniowo niezale˙znych rozwi
֒
aza´
n szczeg´olnych
r´ownania (5.11) nazywamy fundamentalnym uk ladem rozwi
֒
aza´
n (SRRLJ) rz
֒
edu
n.
Twierdzenie 15.
Ka˙zde r´ownanie r´o˙zniczkowe liniowe jednorodne rz
֒
edu n-tego
(5.11) ma fundamentalny uk lad rozwi
֒
aza´
n.
Dow´
od
Niech A = a
j
k
∈ R
n
2
b
֒
edzie dowoln
֒
a macierz
֒
a nioeosobliw
֒
a i niech
t
0
∈ I. Wiadomo, ˙ze r´ownanie (5.11) ma rozwi
֒
azania globalne przy zadanych
warunkach pocz
֒
atkowych Cauchy
֒
ego
x
(k−1)
j
(t
0
) = a
j
k
,
k = 1, . . . , n.
Oznaczmy je symbolami x
j
, (j = 1, . . . , n). Z konstrukcji tych rozwi
֒
aza´
n wynika,
˙ze
W (x
1
, . . . , x
n
) (t
0
) = det A
6= 0
i wobec poprzedniego twierdzenia rozwi
֒
azania x
1
, . . . , x
n
tworz
֒
a fundamentalny
uk lad rozwi
֒
aza´
n.
c.k.d
Twierdzenie 16.
Je˙zeli rozwi
֒
azania x
1
, . . . , x
n
tworz
֒
a fundamentalny uk lad rozwi
֒
aza´
n
jednorodnego r´ownania r´o˙zniczkowego liniowego rz
֒
edu n (5.11), to rodzina funkcji
x =
n
X
j=1
C
j
x
j
,
gdzie C
j
, (j = 1, . . . , n) jest rozwi
֒
azaniem og´olnym tego r´ownania.
Dow´
od
Nale˙zy pokaza´c, ˙ze dla dowolnego rozwi
֒
azania szczeg´olnego x spe lniaj
֒
acego
warunek pocz
֒
atkowy Cauchy’ego
x
(k−1)
(t
0
) = x
0k
(k = 1, . . . , n)
28
5.7. Obni˙zanie rz
֒
edu r´
ownania liniowego
istniej
֒
a sta le
C
j
(j = 1, . . . , n)
takie, ˙ze
x =
P
n
j=1
C
j
x
j
.
Rozwa˙zmy uk lad r´owna´
n
n
X
j=1
C
j
x
(k−1)
j
(t
0
) = x
0k
(k = 1, . . . , n).
Macierz tego uk ladu jest nieosobliwa, bo jej wyznacznik jest r´owny W (x
1
, . . . , x
n
) (t
0
)
6=
0. Niech rozwi
֒
azaniem tego uk ladu b
֒
edzie wektor e
C =
e
C
1
, . . . , e
C
n
T
. Latwo
zauwa˙zy´c, ˙ze sk ladowe e
C
j
tego wektora s
֒
a poszukiwanymi sta lymi.
c.k.d
5.7. Obni ˙zanie rz
֒
edu r´
ownania liniowego
5.7.1. Wz´
or Liouville’a
Rozwa˙zmy teraz jednorodne r´ownanie r´o˙zniczkowe liniowe (5.11) rz
֒
edu dru-
giego. Mo˙zna pokaza´c nast
֒
epuj
֒
ace twierdzenie Liouville’a:
Twierdzenie 17.
Je´sli x
1
, x
2
stanowi
֒
a uk lad fundamentalny rozwi
֒
aza´
n jedno-
rodnego r´ownania r´o˙zniczkowego liniowego (5.11) rz
֒
edu drugiego, to
∃
C∈R
W (x
1
, x
2
) (t) = C exp
−
Z
a
1
(t) dt
.
Je´sli x
1
jest znanym rozwi
֒
azaniem r”wnania (5.11), to drugie rozwi
֒
azanie
niezale˙zne mo˙zna znale˙z´c nast
֒
epuj
֒
acym sposobem:
∀
t∈R
x
1
(t) x (t)
x
′
1
(t) x
′
(t)
6= 0,
(5.15)
x
1
x
′
− x
′
1
x = C exp
−
Z
a
1
(t) dt
,
x
1
x
′
− x
′
1
x
x
2
1
=
1
x
2
1
C exp
−
Z
a
1
(t) dt
,
d
dt
x
x
1
=
1
x
2
1
C exp
−
Z
a
1
(t) dt
,
x
x
1
=
Z
1
x
2
1
C exp
−
Z
a
1
(t) dt
dt,
x (t) = x
1
(t)
Z
1
x
2
1
(t)
C exp
−
Z
a
1
(t) dt
dt + C
1
. (5.16)
29
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
5.7.2. R´
ownania wy ˙zszych rz
֒
ed´
ow
Druga metoda, bardziej uniwersalna metoda, to zastosowanie podstawienia:
x (t) = x
1
(t) y (t) .
(5.17)
Ma bowiem miejsce nast
֒
epuj
֒
ace
Twierdzenie 18.
Je˙zeli
e
x(t)
6= 0
jest rozwi
֒
azanie jednorodnego liniowego
r´ownania r´o˙zniczkowego (5.11) rz
֒
edu n, to po podstawieniu x(t) = e
x(t)y(t) otrzy-
mujemy r´ownanie, kt´orego rz
֒
ad mo˙zna obni˙zy´c do rz
֒
edu n
− 1.
5.8. Niejednorodne r´
ownanie r´
o ˙zniczkowe liniowe
n-tego rz
֒
edu
Definicja 24.
Niejednorodnym r´ownaniem r´o˙zniczkowym liniowym rz
֒
edu n na-
zywamy r´ownanie postaci
L(t)x = g(t),
(5.18)
gdzie g : R
⊃ I → R jest funkcj
֒
a ci
֒
ag l
֒
a.
Za l´o˙zmy, ˙ze znamy uk lad fundamentalny
{x
1
, . . . , x
n
} skojarzonego jednorod-
nego r´ownania r´o˙zniczkowego (5.12). Ca lk
֒
e szczeg´oln
֒
a r´ownania niejednorodnego
(5.18) znajdziemy metod
֒
a uzmienniania sta lych (metod
֒
a Lagrange’a).
Zak ladamy, ˙ze poszukiwane rozwi
֒
azanie jest postaci
x(t) =
n
X
j=1
C
j
(t)x
j
(t).
Funkcje C
j
(t) wyznaczamy rozwi
֒
azuj
֒
ac uk lad r´owna´
n r´o˙zniczkowych
x
1
(t)
. . .
x
n
(t)
x
′
1
. . .
x
′
n
(t)
...
...
x
(n−2)
1
. . . x
(n−2)
n
(t)
x
(n−1)
1
. . . x
(n−1)
n
(t)
C
′
1
(t)
C
′
2
(t)
...
C
′
n−1
(t)
C
′
n
(t)
=
0
0
...
0
g(t)
.
Rozwi
֒
azuj
֒
ac powy˙zszy uk lad dostajemy n r´owna´
n o zmiennych rozdzielonych
C
′
j
(t) = F
j
(t) (j = 1, . . . , n),
gdzie funkcje F
j
s
֒
a okre´slone wzorami Cramera.
Twierdzenie 19.
(Zasada superpozycji) Je´sli funkcja x
1
(t) jest rozwi
֒
azaniem
r´ownania L(t)x = g
1
(t), a x
2
(t) rozwi
֒
azaniem L(t)x = g
2
(t), to x
1
(t) + x
2
(t) jest
rozwi
֒
azaniem r´ownania L(t)x = g
1
+ g
2
(t).
Uzasadnienie tego faktu zostanie przedstawiony przy omawianiu metody uzmien-
niania sta lych dla uk ladu r´owna´
n r´o˙zniczkowych liniowych.
30
5.9. R´
ownanie liniowe n-tego rz
֒
edu o sta lych wsp´
o lczynnikach
5.9. R´
ownanie liniowe
n-tego rz
֒
edu o sta lych
wsp´
o lczynnikach
Rozwa˙zamy r´ownanie postaci
x
(n)
+ a
n−1
x
(n−1)
+ . . . + a
1
x
′
+ a
0
x = 0,
(5.19)
w kt´orym a
j
∈ R, (j = 0, 1, . . . , n − 1). Niech
L :=
d
n
dt
n
+ a
n−1
d
n−1
dt
n−1
+ . . . + a
1
d
dt
+ a
0
,
w´owczas r´ownanie (5.19) mo˙zna zapisa´c kr´otko
Lx = 0.
(5.20)
Przewidujemy rozwi
֒
azanie r´ownania (5.19) w postaci x(t) = e
λt
, gdzie λ
∈ C.
Po wsrawieniu pochodnych x
(j)
(t) = λ
j
e
λt
do (5.19) i wydzieleniu przez e
λt
do-
stajemy:
λ
n
+ a
n−1
λ
n−1
+ . . . + a
0
= 0.
(5.21)
Wniosek 2.
Funkcja x(t) = e
λt
jest rozwi
֒
azaniem r´ownania r´o˙zniczkowego (5.19)
wtedy i tylko wtedy, gdy λ jest pierwiastkiem r´ownania (5.21) zwanego r´ownaniem
charakterystycznym
.
Uwaga 2.
Funkcja zespolona x(t) jest rozwi
֒
azaniem r´ownania r´o˙zniczkowego
(5.19) wtedy i tylko wtedy, gdy
ℜe x(t) oraz ℑm x(t) s
֒
a rozwi
֒
azaniami tego r´ownania.
Niech λ
1
, . . . , λ
n
∈ C b
֒
ed
֒
a wszystkimi pierwiastkami r´ownania charaktery-
stycznego (5.21), przy czym pierwiastek k-krotny wyst
֒
epuje w tym ci
֒
agu k razy.
Funkcje x
j
(t) = e
λ
j
t
maj
֒
a wro´
nskian
W (x
1
, . . . , x
n
) (t) = e
(λ
1
+...+λ
n
)t
1
1
. . . 1
λ
1
λ
2
. . . λ
n
λ
2
1
λ
2
2
. . . λ
2
n
...
...
...
λ
n−1
1
λ
n−1
2
. . . λ
n−1
n
=
= e
(λ
1
+...+λ
n
)t
n
Y
k=1
n
Y
j=k+1
(λ
j
− λ
k
) .
Macierz wyznacznika wyst
֒
epuj
֒
acego w ostatnim wzorze nasi nazw
֒
e macierzy Van-
dermonde’a.
Mog
֒
a zaistnie´c cztery przypadki:
1. Wielomian charakterystyczny ma n r´o˙znych pierwiastk´ow rzeczywistych tj.:
∀
i∈{1,...,n}
λ
i
∈ R oraz
∀
i,j∈{1,...,n}
i
6= j ⇒ λ
i
6= λ
j
.
31
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
Wtedy ∀
t∈R
W (x
1
, . . . , x
n
) (t)
6= 0, zatem rodzina funkcji
x(t) =
n
X
j=1
C
j
e
λ
j
t
jest ca lk
֒
a og´oln
֒
a r´ownania (5.19).
2. Wielomian charakterystyczny ma n r´o˙znych pierwiastk´ow, ale nie wszystkie
pierwiastki s
֒
a rzeczywiste tj.:
∀
i∈{1,...,n}
λ
i
∈ C oraz
∀
i,j∈{1,...,n}
i
6= j ⇒ λ
i
6= λ
j
.
Niech np. λ
m
= a + ib b
֒
edzie jednym z pierwiastk´ow zespolonych. Po-
niewa˙z wielomian charakterystyczny (5.21) ma wsp´o lczynniki rzeczywiste, za-
tem r´ownie˙z λ
m
= a
− ib musi by´c pierwiastkiem tego wielomianu. Mo˙zna
bez szkody dla og´olno´sci przyj
֒
a´c, ˙ze jest to kolejny pierwiastek na li´scie pier-
wiastk´ow tj. λ
m+1
= λ
m
. Par
֒
e liniowo niezale˙znych rozwi
֒
aza´
n zespolonych
y
1
(t) = e
λ
m
t
,
y
2
(t) = e
λ
m+1
t
= e
λ
m
t
zast
֒
epujemy par
֒
a liniowo niezale˙znych rozwi
֒
aza´
n rzeczywistych
x
m
(t) =
ℜe e
λ
m
t
= e
at
cos(bt),
x
m+1
(t) =
ℑm e
λ
m
t
= e
at
sin(bt).
3. Wielomian charakterystyczny (5.21) ma wszystkie pierwiastki rzeczywiste, ale
s
֒
a w´sr´od nich pierwiastki wielokrotne. W tej sytuacji W e
λ
1
t
, . . . , e
λ
n
t
= 0.
Niech λ
m
∈ R b
֒
edzie pierwiastkiem krotno´sci k > 1. W´owczas funkcje
t
0
e
λ
m
t
= e
λ
m
t
, t
1
e
λ
m
t
, . . . , t
k−1
e
λ
m
t
s
֒
a liniowo niezale˙zne, ponadto ka˙zda z nich jest rozwi
֒
azaniem (5.19). Jak
latwo bowiem sprawdzi´c bezpo´srednim rachunkiem
d
dt
− λ
t
s
e
λt
= st
s−1
e
λt
,
sk
֒
ad wniosek, ˙ze je´sli λ jest pierwiastkiem k krotnym i s
≤ k − 1, to
d
dt
− λ
k
t
s
e
λt
= 0.
5.10. Metoda przewidywa´
n
W przypadku niejednorodnego r´ownania r´o˙zniczkowego rz
֒
edu n o sta lych
wsp´o lczynnikach
x
(n)
+ a
n−1
x
(n−1)
+ . . . + a
1
x
′
+ a
0
x = g(t),
(5.22)
32
5.11. Uk lad skalarnych r´
owna´
n r´
o˙zniczkowych liniowych rz
֒
edu pierwszego
mo˙zliwe jest skonstruowanie ca lki szczeg´olnej tego r´ownania, je´sli
g(t) = e
at
(p
k
(t) cos bt + q
m
(t) sin bt) ,
gdzie p
k
i q
m
s
֒
a wielomianami odpowiednio stopnia k i m. Rozwi
֒
azanie szczeg´olne
przewidujemy w postaci
x(t) = e
at
t
p
(r
l
(t) cos bt + s
l
(t) sin bt) ,
gdzie:
— p jest krotno´sci
֒
a pierwiastka a + ib wielomianu charakterystycznego r´ownania
jednorodnego skojarzonego z (5.22); gdy a+ib nie jest pierwiastkiem, to p = 0,
— l = max
{k, l},
— r
l
, s
l
wielomiany stopnia l.
Wsp´o lczynniki wielomian´ow r
l
, s
l
dobieramy metod
֒
a wsp´o lczynnik´ow nieozna-
czonych.
5.11. Uk lad skalarnych r´
owna´
n r´
o ˙zniczkowych liniowych
rz
֒
edu pierwszego
Rozwa˙zamy uk lad r´owna´
n r´o˙zniczkowych rz
֒
edu pierwszego postaci
x
′
j
(t) =
n
X
k=1
a
k
j
(t)x
k
(t) + g
j
(t) (j = 1, . . . , n),
(5.23)
czyli
x
′
(t) = A(t)x(t) + g(t),
(5.24)
gdzie
x(t) =
x
1
(t)
...
x
n
(t)
,
A(t) =
a
1
1
(t) . . . a
n
1
(t)
...
...
a
1
n
(t) . . . a
n
n
(t)
,
g(t) =
g
1
(t)
...
g
n
(t)
.
Przyjmujemy za lo˙zenia regularno´sciowe takie jak w teorii dotycz
֒
acej zagadnie´
n
liniowych. W tym przypadku oznacza to, ˙ze
∀
j,k∈{1,...,n}
I
∋ t → g
j
(t),
I
∋ t → a
k
j
(t)
∈ C (I, R)
∀
j∈{1,...,n}
I
∋ t → x
j
(t)
∈ C
1
(I, R) ,
gdzie I
⊂ R jest przedzia lem.
Niech M
∈ R
n×n
. Definiujemy
M
0
:= I
M
1
:= M
M
j
:= M
· M
j−1
33
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
oraz
e
M
:=
∞
X
k=0
M
k
k!
.
Szereg ten jest zbie˙zny w R
n
2
dla ka˙zdej macierzy M . Wynika to st
֒
ad, ˙ze wobec
oszacowania
M
k
=
M · M
k−1
≤ kMk
M
k−1
≤ kMk
k
mamy nier´owno´s´c
∞
X
k=0
M
k
k!
≤
∞
X
k=0
kMk
k
k!
.
Szereg
P
∞
k=0
kMk
k
k!
jest zbie˙zny, a zatem szereg
P
∞
k=0
M
k
k!
jest zbie˙zny, gdzy˙z
w przestrzeniach Banacha zachodzi twierdzenie, ˙ze szereg kt´ory jest zbie˙zny
wzgl
֒
edem normy (czyli jezt zbie˙zny bezwzgl
֒
ednie) jest zbie˙zny.
Twierdzenie 20.
Je´sli macierze M, N
∈ R
n×n
s
֒
a przemienne, to znaczy gdy
M N = N M , to e
M +N
= e
M
e
N
.
Wniosek 3.
Dla dowolnej macierzy M
∈ R
n×n
:
e
M
−1
= e
−M
.
Dow´
od.
Macierze M i
−M s
֒
a przemienne, a zatem e
M
e
−M
= e
M −M
= e
0
= I.
Mno˙z
֒
ac ten zwi
֒
azek lewostronnie przez e
M
−1
dostajemy tez
֒
e.
c.k.d.
Niech A(t) = a
k
j
(t)
j,k=1,...,n
. Wprowadzamy oznaczenie
Z
A(t) dt :=
Z
a
k
j
(t) dt
j,k=1,...,n
.
Twierdzenie 21.
Je´sli macierze A(t) i
R
A(t) dt s
֒
a przemienne, to funkcja
x(t) := e
R
A(t) dt
C,
(5.25)
gdzie C = (C
1
, . . . , C
n
)
T
∈ R
n
jest rozwi
֒
azaniem jednorodnego uk ladu r´owna´
n
r´o˙zniczkowych liniowych rz
֒
edu pierwszego
x
′
(t) = A(t)x(t).
(5.26)
Dow´
od.
Policzmy:
x
′
(t) =
e
R
A(t) dt
C
′
=
e
R
A(t) dt
′
C =
∞
X
k=0
R
A(t) dt
k
k!
!
′
C =
=
∞
X
k=1
k
R
A(t) dt
k−1
R
A(t)dt
′
k!
!
C =
∞
X
k=1
R
A(t) dt
k−1
A(t)
(k
− 1)!
!
C =
= A(t)
∞
X
k=0
R
A(t) dt
k
k!
!
C = A(t)e
A(t)
C = A(t)x(t).
34
5.11. Uk lad skalarnych r´
owna´
n r´
o˙zniczkowych liniowych rz
֒
edu pierwszego
c.k.d.
Wz´or (5.25) ma niewielkie znaczenie praktyczne, gdy macierz uk ladu A(t)
zale˙zy istotnie od zmiennej t.
Uwaga 3.
Je´sli macierz uk ladu (5.26) jest sta la tj. A(t) = A, to
R
A(t) dt =
R
A dt = tA, a zatem macierze A i
R
A dt = tA s
֒
a przemienne. W konsekwencji
rozwi
֒
azaniem uk ladu
x
′
= Ax
jest funkcja
x(t) = e
tA
C
Jak si
֒
e dalej oka˙ze efektywne obliczenie macierzy e
tA
b
֒
edzie mo˙zliwe.
W podobny spos´ob jak przedstawiony powy˙zej, mo˙zna pokaza´c, ˙ze funkcja
x(t) = e
R
A(t) dt
Z
e
−
R
A(t) dt
g(t) dt + e
R
A(t) dt
C
(C
∈ R
n
),
jest rozwi
֒
azaniem og´olnym niejednorodnego uk ladu (5.24). Wektor C dla rozwi
֒
azania
spe lniaj
֒
acego warunek pocz
֒
atkowy Cauchy’ego x (t
0
) =
0
x ma posta´c:
C = e
−
R
A(t) dt 0
x −
Z
e
−
R
A(t) dt
g(t) dt.
Twierdzenie 22.
Niech funkcje x
k
j
∈ C
1
(I, R) (j, k = 1, . . . , n), niech x
k
ozna-
cza wektor x
k
:= x
k
1
, . . . , x
k
n
T
i niech
D x
1
, . . . , x
n
(t) := det
x
k
j
(t)
j,k=1,...,n
.
a) Je´sli D (x
1
, . . . , x
n
) (t)
6= 0 dla pewnego t
0
∈ I, to x
1
, . . . , x
n
s
֒
a liniowo
niezale˙zne.
b) Je´sli x
1
, . . . , x
n
s
֒
a liniowo niezale˙znymi rozwi
֒
azaniami jednorodnego uk ladu
(5.26), to
∀
t∈I
D (x
1
, . . . , x
n
) (t)
6= 0.
Dow´
od.
Ad a). (Nie wprost) Przyjmijmy, ˙ze x
1
, . . . , x
n
liniowo zale˙zne tzn.
istniej
֒
a takie sta le C
1
, . . . , C
n
, ˙ze
P
n
k=1
C
2
k
6= 0 oraz ∀
t∈I
P
n
k=1
C
k
x
k
(t) = 0.
To jednak oznacza, ˙ze det
x
k
j
(t)
j,k=1,...,n
= D (x
1
, . . . , x
n
) (t
0
) = 0, wbrew
za lo˙zeniu.
Ad b). (Nie wprost) Dla dowodu nie wprost przyjmijemy, ˙ze D (x
1
, . . . , x
n
) (t
0
) =
0 dla pewnego t
0
∈ I. Niech wektor (C
1
, . . . , C
n
)
T
b
֒
edzie niezerowym rozwi
֒
azaniem
uk ladu
x
1
1
(t
0
) , . . . x
n
1
(t
0
)
...
...
x
1
n
(t
0
) , . . . x
n
n
(t
0
)
·
C
1
...
C
n
=
0
...
0
.
35
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
Zdefiniujmy funkcj
֒
e x(t) jako
x(t) :=
n
X
k=1
C
k
x
k
(t),
t
∈ I.
Jako kombinacja liniowa rozwi
֒
aza´
n x
k
funkcja x jest rozwi
֒
azaniem uk ladu (5.24).
Ponadto spe lnia ona warunek pocz
֒
atkowy Cauchy’ego
x (t
0
) = 0.
Funkcja y(t)
≡ 0 jest r´ownie˙z rozwi
֒
azaniem uk ladu (5.24) spe lniaj
֒
acym ten sam
warunek pocz
֒
atkowy. Wobec jednoznaczno´sci rozwi
֒
azania funkcje te musz
֒
a by´c
r´owne, czyli x = 0. Oznacza to jednak wbrew za lo˙zeniu, ˙ze funkcje x
1
, . . . , x
n
s
֒
a
liniowo zale˙zne.
c.k.d.
Uwaga 4.
Je˙zeli x
1
, . . . , x
n
s
֒
a rozwi
֒
azaniami uk ladu (5.26), to
∀
t∈I
D x
1
, . . . , x
n
(t) = 0,
lub
∀
t∈I
D x
1
, . . . , x
n
(t)
6= 0.
Definicja 25.
Zbi´or
{x
1
, . . . , x
n
} liniowo niezale˙znych rozwi
֒
aza´
n uk ladu (5.26)
nazywamy fundamentalnym uk ladem rozwi
֒
aza´
n.
Twierdzenie 23.
Ka˙zdy jednorodny uk lad r´owna´
n r´o˙zniczkowych liniowych ma
fundamentalny uk lad rozwi
֒
aza´
n i je´sli funkcje x
1
, . . . , x
n
tworz
֒
a fundamentalny
uk lad rozwi
֒
aza´
n, to rodzina odwzorowa´
n x(t) =
P
n
k=1
C
k
x
k
(t), gdzie C
k
∈ R jest
rozwi
֒
azaniem og´olnym tego uk ladu.
Je´sli
{x
1
, . . . , x
n
} tworz
֒
a fundamentalny uk lad rozwi
֒
aza´
n (5.26), to ca lk
֒
e
szczeg´oln
֒
a niejednorodnego uk ladu (5.24) znajdujemy metod
֒
a uzmienniania sta lych.
Przewidujemy j
֒
a w postaci
x(t) :=
n
X
k=1
C
k
(t)x
k
(t).
Dalej mamy x
′
(t) :=
P
n
k=1
C
′
k
(t)x
k
(t) +
P
n
k=1
C
k
(t) x
k
′
(t) i po wstawieniu do
r´ownania otrzymujemy:
n
X
k=1
C
′
k
(t)x
k
(t) +
n
X
k=1
C
k
(t) x
k
′
(t) = A(t)
n
X
k=1
C
k
(t)x
k
(t) + g(t),
czyli
n
X
k=1
C
′
k
(t)x
k
(t) = g(t),
36
5.12. Uk lady r´
owna´
n liniowych o sta lych
wsp´
o lczynnikach
to jest
x
1
1
(t) x
2
1
(t)
· · · x
n
1
(t)
x
1
2
(t) x
2
2
(t)
· · · x
n
2
(t)
...
...
...
x
1
n
(t) x
2
n
(t)
· · · x
n
n
(t)
C
′
1
(t)
C
′
2
(t)
...
C
′
n
(t)
=
g
1
(t)
g
2
(t)
...
g
n
(t)
.
Poniewa˙z dla wszystkich t
∈ I :
D (x
1
, . . . , x
n
) (t)
6= 0, st
֒
ad powy˙zszy uk lad
ma dok ladnie jedno rozwi
֒
azanie okre´slone wzorami Cramera
C
′
k
(t) = p
k
(t) (k = 1, . . . , n).
Ka˙zde z tych r´owna´
n jest r´ownaniem o zmiennych rozdzielonych zatem
C
k
(t) =
Z
p
k
(t) dt + M
k
,
gdzie M
k
∈ R, (k = 1, . . . , n).
Ostatecznie
x(t) =
n
X
k=1
M
k
x
k
(t) +
n
X
k=1
Z
p
k
(t) dt
x
k
(t).
5.12. Uk lady r´
owna´
n liniowych o sta lych
wsp´
o lczynnikach
Zak ladamy teraz, ˙ze macierz uk ladu (5.26) jest macierz
֒
a sta l
֒
a tj. a
k
j
(t)
≡ a
k
j
∈
R
. Jak wiadomo z wcze´sniejszych rozwa˙za´
n, rozwi
֒
azanie tego uk ladu jest postaci
x(t) = e
tA
C,
gdzie C
∈ R
n
.
5.12.1. Metoda warto´
sci i wektor´
ow w lasnych
Je´sli w
6= 0 jest wektorem w lasnym macierzy A tj. istnieje λ ∈ C :
Aw = λw
i we´zmiemy x(t) = y(t)
· w, gdzie y(t) ∈ C
1
(R, R), to po podstawieniu x do
r´ownania (5.26) dostajemy y
′
(t)w = λy(t)w co daje (y
′
(t)
− λy(t)) w = 0. Wobec
w
6= 0 mamy y
′
(t) = λy(t) r´ownanie o zmiennych rozdzielonych z rozwi
֒
azaniem
y(t) = Ce
λt
,
t
∈ R.
Jak wiadomo zbi´or rozwi
֒
aza´
n uk ladu (5.26) jest przestrzeni
֒
a wektorow
֒
a n-wymiarow
֒
a.
Poszukujemy zatem fundamentalnego uk ladu rozwi
֒
aza´
n. Mo˙zemy rozwa˙zy´c przy-
padki:
1. Ka˙zdej warto´sci w lasnej λ
j
o krotno´sci k
j
odpowiada k
j
liniowo niezale˙znych
wektor´ow w lasnych w
j,1
, . . . , w
j,k
j
macierzy A
(j = 1, . . . , p,
k
1
+ k
2
+
. . . + k
p
= n). Poniewa˙z wektory w lasne odpowiadaj
֒
ace r´o˙znym warto´sciom
w lasnym s
֒
a liniowo niezale˙zne, wi
֒
ec dla
x
j,s
(t) := e
λ
j
t
w
j,s
(s = 1, . . . , k
j
, j = 1, . . . , p)
37
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
wyznacznik
D x
1,1
, . . . , x
p,k
p
(t) = e
(k
1
λ
1
+...+k
p
λ
p
)t
det w
j,s
i
6= 0,
gdzie i = 1, . . . , n, s = 1, . . . , k
j
, j = 1, . . . , p. W konsekwencji funkcje x
j,s
tworz
֒
a fundamentalny uk lad rozwi
֒
aza´
n.
Je´sli warto´s´c w lasna i wektor w lasny s
֒
a zespolone tj. np. λ
1
= λ
2
= λ i w =
u+iv = (u
1
+ iv
1
, . . . , u
n
+ iv
n
)
T
jest wektorem w lasnym odpowiadaj
֒
acym λ
1
i w = u
− iv = (u
1
− iv
1
, . . . , u
n
− iv
n
)
T
wektorem w lasnym odpowiadaj
֒
acym
λ
2
, to poniewa˙z r´owno´s´c Aw = λw poci
֒
aga r´owno´s´c Aw = Aw = Aw = λw =
λw, zatem zamiast zespolonych rozwi
֒
aza´
n
y
1
= e
λ
1
t
w,
y
2
= e
λ
2
t
w
bierzemy
x
1
=
ℜe y
1
= e
tℜe λ
(u cos (t
ℑm λ) − v sin (tℑm λ)) ,
x
2
=
ℑm y
1
= e
tℜe λ
(u sin (t
ℑm λ) + v cos (tℑm λ)) .
2. Niech warto´sci w lasnej np. λ
1
= λ o krotno´sci k odpowiada tylko r liniowo
niezale˙znych wektor´ow w lasnych, gdzie r < k. Tak jest wtedy, gdy
rz
֒
ad (A
− λI) = n − r > n − k.
Poszukujemy rozwi
֒
azania og´olnego odpowiadaj
֒
acego warto´sci w lasnej λ po-
staci
x(t) = e
λt
P (t),
gdzie P (t) = (P
1
(t), . . . , P
n
(t))
T
i P
j
jest wielomianem stopnia k
− 1, j =
1, . . . , n przy czym w rozwi
֒
azaniu og´olnym powinno wyst
֒
api´c k sta lych do-
wolnych.
3. Rozwi
֒
azanie og´olne uk ladu (5.26) jest sum
֒
a rozwi
֒
aza´
n szczeg´olnych odpowia-
daj
֒
acych poszczeg´olnym warto´sciom w lasnym.
5.12.2. Sprowadzanie macierzy uk ladu do postaci Jordana
Przypadek szczeg´
olny
Je´sli A jest diagonalizowaln
֒
a rzeczywist
֒
a macierz
֒
a
wymiaru n, tj. istnieje macierz podobie´
nstwa P taka, ˙ze P
−1
AP = D, gdzie
D jest macierz
֒
a diagonaln
֒
a, to podstawiaj
֒
ac x = P y sprowadzamy uk lad x
′
(t) =
Ax(t),
t
∈ I
(
∈ topR) do postaci y
′
(t) = Dy(t), kt´orego rozwi
֒
azaniem jest
y(t) = e
Dt
C = e
d
ii
t
C =
C
1
e
d
11
t
...
C
n
e
d
nn
t
,
38
5.12. Uk lady r´
owna´
n liniowych o sta lych
wsp´
o lczynnikach
a zatem
x (t) = P
C
1
e
d
11
t
...
C
n
e
d
nn
t
.
W sczeg´olno´sci, je´sli macierz A ma n r´o˙znych warto´sci w lasnych λ
i
(i = 1, . . . , n),
to jest diagonalizowalna i D = diag
{λ
1
, . . . , λ
n
}. Wtedy te˙z macierz podo-
bie´
nstwa P jest r´owna
P = (v
1
, . . . , v
n
) ,
gdzie v
i
jest wektorem w lasnym odpowiadaj
֒
acym warto´sci w lasnej λ
i
(i = 1, . . . , n).
Przypadek og´
olny
Niech A b
֒
edzie dan
֒
a rzeczywist
֒
a macierz
֒
a kwadratow
֒
a
wymiaru n. Niech λ
r
(r = 1, . . . , q) b
֒
ed
֒
a warto´sciami w lasnymi tej macie-
rzy, przy czym przyjmujemy, ˙ze warto´s´c w lasna λ
r
ma krotno´s´c k
r
. Oczywi´scie
P
q
r=1
k
r
= n. Niech P b
֒
edzie tak
֒
a macierz
֒
a nieosobliw
֒
a, ˙ze macierz
J =
P
−1
AP
jest macierz
֒
a Jordana, tzn.
J =
J
11
0
· · · 0
· · · 0
· · · 0
0
J
12
· · · 0
· · · 0
· · · 0
...
...
. ..
...
0
0
J
1i(1)
0
0
...
...
. ..
...
0
0
0
J
q1
0
...
...
. .. ...
0
0
· · · 0
· · · 0
· · · J
q,i(q)
,
gdzie
J
rj
=
λ
r
0
0
· · · 0
0
1
λ
r
0
0
0
0
1
λ
r
0
0
...
. .. ... ...
0
0
0
· · · λ
r
0
0
0
0
· · · 1
λ
r
lub
J
rj
= (λ
r
)
(macierze J
rj
nazywamy klatkami Jordana). Oznaczmy przez k
rj
liczb
֒
e wierszy i
kolumn macierzy J
rj
. Obliczaj
֒
ac wielomian charakterystyczny macierzy J, r´owny
wielomianowi charakterystycznemu macierzy A latwo mo˙zna si
֒
e przekona´c, ˙ze
maj
֒
a miejsce nast
֒
epuj
֒
ace r´owno´sci:
k
r
=
i(r)
X
j=1
k
rj
(r = 1, . . . , q) .
Liczby k
rj
mo˙zna wyznaczy´c np. metod
֒
a przedstawion
֒
a w [6].
39
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
Niech D
m
(λ) oznacza najwi
֒
ekszy wsp´olny dzielnik wszystkich minor´ow stop-
nia m macierzy A
− λI. Mo˙zna pokaza´c, ˙ze D
m
(λ) dzieli si
֒
e przez D
m−1
(λ).
Zatem z dok ladno´sci
֒
a do czynnika a, takiego ˙ze
|a| = 1:
D
n
(λ) = (λ
− λ
1
)
u
11
(λ
− λ
2
)
u
21
. . . (λ
− λ
q
)
u
q1
,
D
n−1
(λ) = (λ
− λ
1
)
u
12
(λ
− λ
2
)
u
22
. . . (λ
− λ
q
)
u
q2
,
. . .
........................................................
D
1
(λ) = (λ
− λ
1
)
u
1n
(λ
− λ
2
)
u
2n
. . . (λ
− λ
q
)
u
qn
,
przy czym u
i1
≥ u
i2
≥ u
i3
≥ . . . ≥ u
in
co mo˙zna zapisa´c kr´otko u
ik
≥ u
ij
dla k
≤ j . Nie wykluczamy przypadku, gdy pewne u
ik
= 0. W tym przypadku
jednak u
ij
= 0 dla wszystkich j
≥ k. Przy tych oznaczeniach:
k
11
= u
11
− u
12
,
k
12
= u
12
− u
13
, . . . ,
k
ij
= u
ij
− u
i,j+1
, . . .
Mamy w´owczas i (r) = max
{j :
k
rj
6= 0}.
Je˙zeli k
m
= 1 dla pewnego m, to i (m) = 1 oraz J
m1
= (λ
m
) jest macierz
֒
a
wymiaru 1
× 1. Bez straty og´olno´sci mo˙zemy przyj
֒
a´c, ˙ze je´sli istniej
֒
a pierwiastki
jednokrotne, to maj
֒
a one kolejne numery rozpoczynaj
֒
ace si
֒
e od 1. Macierz Jor-
dana J jest wi
֒
ec postaci
λ
1
0
· · · 0
0
· · · 0
0
λ
2
· · · 0
0
· · · 0
... ... ...
...
0
0
λ
p
0
0
0
0
0
J
p+1,1
0
... ...
. .. ...
0
0
· · · 0
0
· · · J
q,i(q)
.
Dowodzi si
֒
e, ˙ze je´sli
A = P JP
−1
,
to
e
At
= P e
Jt
P
−1
.
Z kolei
e
Jt
=
e
λ
1
t
· · · 0
0
· · · 0
...
. ..
...
0
e
λ
p
t
0
0
0
0
e
J
p+1,1
t
0
...
. .. ...
0
0
· · ·
0
· · · e
J
q,i(q)
t
,
gdzie λ
1
, . . . , λ
p
s
֒
a jednokrotnymi pierwiastkami wielomianu charakterystycznego.
40
5.12. Uk lady r´
owna´
n liniowych o sta lych
wsp´
o lczynnikach
Niech J
rj
wymiaru k
rj
b
֒
edzie jedn
֒
a z klatek Jordana odpowiadaj
֒
acych warto´sci
w lasnej λ
r
o krotno´sci k
r
. Wprost z definicji mo˙zna pokaza´c, ˙ze
e
J
rj
t
= e
λ
r
t
1
0
· · · 0
t
1!
1
· · · 0
t
2
2!
t
1!
0
...
. .. ...
t
krj −1
(k
rj
−1)!
t
krj −2
(k
rj
−2)!
· · · 1
.
Niech P b
֒
edzie macierz
֒
a sprowadzaj
֒
ac
֒
a macierz A do postaci Jordana tj. J =
P
−1
AP . Ostatni zwi
֒
azek jest r´ownowa˙zny r´owno´sci P J = AP . Wprowadzaj
֒
ac
now
֒
a funkcj
֒
e niewiadom
֒
a y (t) okre´slon
֒
a r´owno´sci
֒
a
x(t) = P y(t),
sprowadzamy ostatni URRLJ do r´ownowa˙znego uk ladu
y
′
(t) = Jy(t),
kt´orego rozwi
֒
azaniem og´olnym jest funkcja
y(t) = e
Jt
C, C
∈ R
n
.
Tak wi
֒
ec rozwi
֒
azaniem og´olnym wyj´sciowego URRLJ jest funkcja
x(t) = P e
Jt
C, C
∈ R
n
.
Przyk lad 5.
Rozwa˙zmy uk lad r´owna´
n:
x
′
(t) =
1 1 2
0 1 1
0 0 2
x (t)
Jak latwo sprawdzi´c λ
1
= 1,
k
1
= 2
λ
2
= 2,
k
2
= 1
P =
1 1 3
1 0 1
0 0 1
J =
1 0 0
1 1 0
0 0 2
Stosuj
֒
ac standardowe podstawienie x (t) = P y (t) rozwi
֒
azujemy uk lad r´owna´
n
y
′
(t) = Jy(t). Jego rozwi
֒
azaniem jest
y(t) = e
Jt
C =
e
t
1 0
t 1
0
0
e
2t
C
1
C
2
C
3
=
e
t
0
0
te
t
e
t
0
0
0 e
2t
C
1
C
2
C
3
,
41
Rozdzia l 5. Liniowe r´
ownania r´
o˙zniczkowe
zatem
x(t) =
1 1 3
1 0 1
0 0 1
e
t
0
0
te
t
e
t
0
0
0 e
2t
C
1
C
2
C
3
=
=
(1 + t) e
t
e
t
3e
2t
e
t
0
e
2t
0
0
e
2t
C
1
C
2
C
3
.
5.13. R´
ownanie ruchu harmonicznego
R´ownanie ruchu pod dzia laniem si ly elastycznej, tj. r´ownanie ruchu harmo-
nicznego jest opisane r´ownaniem r´o˙zniczkowym wektorowym:
m
..
r
=
−k
2
r
,
gdzie r = r (t) = (x (t) , y (t) , z (t)). Podstawiaj
֒
ac ω =
k
√
m
dostajemy uk lad
separowanych r´owna´
n skalarnych:
..
x +ω
2
x = 0
..
y +ω
2
y = 0
..
z +ω
2
z = 0
.
Ca lka og´olna pierwszego z nich ma posta´c:
x (t) = C
1
sin ωt + C
2
cos ωt = A sin (ωt + γ) ,
gdzie A =
p
C
2
1
+ C
2
2
, γ = arctan (C
1
/C
2
). Latwo zauwa˙zy´c, ˙ze rozwi
֒
azanie x (t)
jest okresowe o okresie T =
2π
ω
. Sta l
֒
a ω nazywamy cz
֒
esto´sci
֒
a ko low
֒
a lub pulsacj
֒
a,
ν =
1
T
cz
֒
esto´sci
֒
a, ωt + γ faz
֒
a, za´s γ sta l
֒
a fazow
֒
a.
Je˙zeli na punkt materialny opr´ocz si ly elastycznej
−k
2
x dzia la dodatkowa si la
−ρ
.
x
(ρ > 0), to otrzymujemy drgania t lumione.
Rozdzia l 6
Rozwi
֒
azania w postaci szereg´
ow
funkcyjnych
Jak wiadomo nie zawsze mo˙zna efektywnie rozwi
֒
aza´c r´ownanie r´o˙zniczkowe,
nie zawsze mo˙zna otrzyma´c rozwi
֒
azanie przez sko´
nczon
֒
a liczb
֒
e kwadratur. Cza-
sami trzeba si
֒
egn
֒
a´c do sposob´ow bardziej wyrafinowanych - jednym z nich jest
wyra˙zenie rozwi
֒
azania w postaci szeregu funkcyjnego. Poni˙zej om´owione s
֒
a dwa
przypadki takiego post
֒
epowania.
6.1. Rozwi
֒
azania w postaci szereg´
ow pot
֒
egowych
Niech b
֒
edzie dane zagadnienie pocz
֒
atkowe Cauchy’ego
x
′
= f (t, x)
(t
∈ I) ,
x (t
0
) = x
0
,
gdzie I
⊂ R przedzia l, taki ˙ze t
0
∈
◦
I, x : I
∋ t → x(t) ∈ U ⊂ R, U zbi´or
otwarty w R, x
0
∈ U. Twierdzenie Cauchy-Kowalewskiej orzeka, ˙ze je´sli funkcja
f : I
× U → R jest analityczna w otoczeniu punktu (t
0
, x
0
), to istnieje dok ladnie
jedno analityczne rozwi
֒
azanie tego r´ownania w pewnym otoczeniu punktu t
0
.
43
Rozdzia l 6. Rozwi
֒
azania w postaci szereg´
ow funkcyjnych
6.1.1. Uk lad r´
owna´
n liniowych rz
֒
edu pierwszego o sta lych
wsp´
o lczynnikach
Rozwa˙zamy uk lad r´owna´
n r´o˙zniczkowych rz
֒
edu pierwszego (5.23) o sta lych
wsp´o lczynnikach czyli uk lad postaci
x
′
j
(t) =
n
X
k=1
a
k
j
x
k
(t) + g
j
(t) (j = 1, . . . , n).
(6.1)
Przyjmijmy, ˙ze
g
j
(t) =
∞
X
ν=0
c
jν
(t
− t
0
)
ν
(j = 1, . . . , n).
Szukamy rozwi
֒
azania x : R
⊃ I → R
n
, kt´orego wszystkie sk ladowe s
֒
a szeregami
pot
֒
egowymi o ´srodku w punkcie t
0
:
x
j
(t) =
∞
X
ν=0
b
jν
(t
− t
0
)
ν
(j = 1, . . . , n).
Podstawiaj
֒
ac szeregi g
j
(t), x
j
(t) i
x
′
j
(t) =
∞
X
ν=0
(ν + 1)b
j,ν+1
(t
− t
0
)
ν
(j=1,. . . ,n) do r´ownania (6.1), przegrupowuj
֒
ac wyrazy i korzystaj
֒
ac z definicji
r´owno´sci szereg´ow pot
֒
egowych dostajemy zwi
֒
azki rekurencyjne na wsp´o lczynniki
szereg´ow x
j
(t):
b
j,ν+1
=
1
ν + 1
n
X
k=1
a
k
j
b
kν
+ c
jν
!
.
(6.2)
Wsp´o lczynniki b
j0
s
֒
a wyznaczone przez warunek pocz
֒
atkowy Cauchy’ego.
Je´sli g
j
(t)
≡ 0, czyli uk lad (6.1) jest jednorodny, to wyznaczaj
֒
ac kolejno
n rozwi
֒
aza´
n uk ladu (6.1) z warunkiem pocz
֒
atkowym x(t
0
) = e
i
, gdzie e
i
jest
wersorem i-tej osi, otrzymujemy fundamentalny uk lad rozwi
֒
aza´
n.
6.1.2. Skalarne r´
ownania r´
o ˙zniczkowe rz
֒
edu pierwszego i drugiego
Nieliniowe r´
ownanie r´
o ˙zniczkowe rz
֒
edu drugiego
Rozwa˙zmy przypadek szczeg´olny, r´ownanie skalarne postaci:
x
′′
= w(g(t), x, x
′
)
t
∈ (t
0
, T )
x (t
0
) = x
0
,
x
′
(t
0
) = x
1
,
gdzie w(p
1
, p
2
, p
3
) jest wielomianem stopnia co najwy˙zej drugiego
w(p
1
, p
2
, p
3
) = a
0
+
3
X
i=1
a
i
p
i
+
3
X
i,j=1
i<=j
a
ij
p
i
p
j
,
44
6.1. Rozwi
֒
azania w postaci szereg´
ow pot
֒
egowych
o wsp´o lczynnikach rzeczywistych, a g jest funkcj
֒
a analityczn
֒
a w otoczeniu punktu
t
0
. Przyjmijmy, ˙ze funkcja g ma rozwini
֒
ecie w szereg pot
֒
egowy
g(t) =
∞
X
k=0
g
k
(t
− t
0
)
k
,
a szukana funkcja x(t) rozwini
֒
ecie
x(t) =
∞
X
k=0
c
k
(t
− t
0
)
k
.
Pierwsza i druga pochodna funkcji szukanej maj
֒
a zatem rozwini
֒
ecia
x
′
(t) =
∞
X
k=0
(k + 1)c
k+1
(t
− t
0
)
k
,
x
′′
(t) =
∞
X
k=0
(k + 2)(k + 1)c
k+2
(t
− t
0
)
k
.
Iloczyny
g(t)g(t),
g(t)x(t),
g(t)x
′
(t),
x(t)x(t),
x(t)x
′
(t),
x
′
(t)x
′
(t)
po prawej stronie r´ownania r´o˙zniczkowego s
֒
a iloczynami Cauchy’ego:
∞
X
k=0
u
k
(t
− t
0
)
k
!
∞
X
k=0
v
k
(t
− t
0
)
k
!
=
∞
X
k=0
k
X
j=0
u
j
v
k−j
!
(t
− t
0
)
k
.
Ostatecznie dostajemy r´owno´s´c dw´och szereg´ow pot
֒
egowych:
P
∞
k=0
(k + 2)(k + 1)c
k+2
(t
− t
0
)
k
=
P
∞
k=0
((δ
0k
a
0
+ a
1
g
k
+ a
2
c
k
+ a
3
(k + 1)c
k+1
) +
+
P
k
j=0
(a
11
g
j
g
k−j
+ a
12
g
j
c
k−j
+ a
13
(k + 1
− j)g
j
c
k+1−j
+ a
22
c
j
c
k−j
+
+ a
23
(k + 1
− j)c
j
c
k+1−j
+ a
33
(j + 1)(k + 1
− j)c
j+1
c
k+1−j
)) (t
− t
0
)
k
,
kt´ora przez por´ownanie wsp´o lczynnik´ow przy tych samych pot
֒
egach (t
− t
0
) pro-
wadzi do niesko´
nczonego uk ladu r´owna´
n algebraicznych o niewiadomych
c
k
(k
∈
N
).
Uwzgl
֒
edniaj
֒
ac warunki pocz
֒
atkowe mamy
c
0
= x
0
,
c
1
= x
1
.
Kolejne wsp´olczynniki c
k
mo˙zna wyznaczy´c rekurencyjnie:
c
k+2
=
1
(k + 1)(k + 2)
δ
0k
a
0
+ a
1
g
k
+ a
2
c
k
+ (k + 1)a
3
c
k+1
+
k
X
j=0
S
kj
!
(k
∈ N),
gdzie
S
kj
:= a
11
g
j
g
k−j
+ c
k−j
(a
12
g
j
+ a
22
c
j
) +
+(k + 1
− j)c
k+1−j
(a
13
g
j
+ a
23
c
j
+ (j + 1)a
33
c
j+1
) ,
a δ
ij
jest delt
֒
a Kroneckera.
45
Rozdzia l 6. Rozwi
֒
azania w postaci szereg´
ow funkcyjnych
Jednorodne liniowe r´
ownanie r´
o ˙zniczkowe rz
֒
edu drugiego
Wyznaczenie rozwi
֒
azania r´ownania liniowego jednorodnego rz
֒
edu drugiego
x
′′
+ p(t)x
′
+ q(t)x = 0
t
∈ (t
0
, T )
x (t
0
) = x
0
,
x
′
(t
0
) = x
1
,
ze wsp´o lczynnikami p(t), q(t) analitycznymi w otoczeniu punktu t
0
wygl
֒
ada po-
dobnie do przedstawionego powy˙zej. Je´sli
p(t) =
∞
X
k=0
a
k
(t
− t
0
)
k
,
q(t) =
∞
X
k=0
b
k
(t
− t
0
)
k
,
to r´ownanie rekurencyjne na wsp´o lczynniki c
k
ma posta´c:
c
k+2
=
−
1
(k + 1)(k + 2)
k
X
j=0
((j + 1)a
k−j
c
j+1
+ b
k−j
c
j
)
(k
∈ N),
(Punkt t
0
, w otoczeniu kt´orego wsp´o lczynniki r´ownania liniowego jednorodnego
s
֒
a funkcjami analitycznymi, nazywamy punktem nieosobliwym tego r´ownania.)
Bior
֒
ac kolejno dwa warunki pocz
֒
atkowe Cauchy’ego x (t
0
) = x
0
, x
′
(t
0
) = x
1
oraz x (t
0
) = x
0
, x
′
(t
0
) = x
1
takie, ˙ze det
x
0
x
1
x
0
x
1
6= 0, mo˙zna wygenerowa´c
dwa liniowo niezale˙zne rozwi
֒
azania tego r´ownania i jego rozwi
֒
azanie og´olne. Naj-
pro´sciej przyj
֒
a´c x
0
= 1, x
1
= 0 oraz x
0
= 0 i x
1
= 1.
Specjalne r´
ownanie Riccatiego
Jeszcze jednym przyk ladem niech b
֒
edzie spos´ob wyznaczenia rozwi
֒
azania spe-
cjalnego r´ownania Riccatiego x
′
(t) = ax
2
(t) + bt
n
z warunkiem pocz
֒
atkowym
Caychy’ego x(0) = x
0
, gdzie a, b
∈ R, n ∈ N. Dla prostoty przyjmijmy
n = 2. Post
֒
epowanie takie jak wy˙zej prowadzi do wzor´ow rekurencyjnych na
wsp´olczynniki rozwi
֒
azania x(t) =
P
∞
k=0
c
k
t
k
:
c
0
= x
0
,
c
1
= ac
2
0
,
c
2
= ac
0
c
1
,
c
3
=
1
3
a 2c
0
c
2
+ c
2
1
+
1
3
c
0
,
c
ν+1
=
1
ν + 1
a
ν
X
k=0
c
k
c
ν−k
ν = 3, 4, 5, . . . .
6.2. R´
ownania r´
o ˙zniczkowe liniowe rz
֒
edu drugiego –
szeregi Frobeniusa
Niech b
֒
edzie dane liniowe r´ownanie r´o˙zniczkowe rz
֒
edu drugiego
x
′′
+ p(t)x
′
+ q(t)x = 0
(t
∈ I) ,
46
6.2. R´
ownania r´
o˙zniczkowe liniowe rz
֒
edu drugiego – szeregi Frobeniusa
gdzie I
⊂ R przedzia l, taki ˙ze t
0
∈
◦
I, x : I
∋ t → x(t) ∈ U ⊂ R, U ∈ topR.
Wiadomo, ˙ze zbi´or rozwi
֒
aza´
n r´ownania jednorodnego jest przestrzeni
֒
a wek-
torow
֒
a dwuwymiarow
֒
a. W przypadku, gdy p(t) = const, q(t) = const, w pro-
sty i znany spos´ob mo˙zna wypisa´c wzory dw´och liniowo niezale˙znych rozwi
֒
aza´
n
tego r´ownania i w konsekwencji dla zadanego warunku pocz
֒
atkowego Cauchy’ego
wyznaczy´c rozwi
֒
azanie problemu pocz
֒
atkowego. Gdy t
0
jest punktem nieosobli-
wym r´ownania tj. p(t), q(t) s
֒
a funkcjami analitycznymi w otoczeniu punktu t
0
,
to mo˙zna wyznaczy´c rozwi
֒
azanie tego problemu w postaci szeregu pot
֒
egowego
o ´srodku w punkcie t
0
, jak to zosta lo pokr´otce opisane powy˙zej, a tak˙ze wy-
znaczy´c dwa szeregi pot
֒
egowe, kt´orych sumy s
֒
a dwoma liniowo niezale˙znymi
rozwi
֒
azaniami r´ownania jednorodnego. Gdy funkcje p(t), q(t) nie s
֒
a analityczne
w otoczeniu punktu t
0
, to punkt ten nazywamy punktem osobliwym r´ownania, a
nazywamy go punktem osobliwym regularnym, je´sli funkcje (t
− t
0
) p(t), (t
− t
0
)
2
q(t)
s
֒
a analityczne w otoczeniu t
0
.
Niech t
0
b
֒
edzie regularnym punktem osobliwym rozwa˙zanego r´ownania i niech
funkcje (t
− t
0
) p(t), (t
− t
0
)
2
q(t) analityczne w otoczeniu
|t − t
0
| < R maj
֒
a roz-
wini
֒
ecia w szeregi pot
֒
egowe:
(t
− t
0
) p(t) =
∞
X
k=0
p
k
(t
− t
0
)
k
,
(t
− t
0
)
2
q(t) =
∞
X
k=0
q
k
(t
− t
0
)
k
.
Niech λ
1
, λ
2
b
֒
ed
֒
a pierwiastkami r´ownania
λ(λ
− 1) + p
0
λ + q
0
= 0,
zwanego r´ownaniem indeksowym (wyznaczaj
֒
acym), gdzie p
0
= lim
t→t
0
(t
− t
0
) p(t),
q
0
= lim
t→t
0
(t
− t
0
)
2
q(t). W jednym z mo˙zliwych przypadk´ow, w sytuacji gdy
λ
1
, λ
2
∈ R, λ
1
> λ
2
, λ
1
− λ
2
6∈ N rozwa˙zane r´ownanie ma dwa liniowo niezale˙zne
rozwi
֒
azania w przedziale (t
0
, t
0
+ R) postaci:
x
1
(t) = (t
− t
0
)
λ
1
∞
X
k=0
a
k
(t
− t
0
)
k
,
x
2
(t) = (t
− t
0
)
λ
2
∞
X
k=0
b
k
(t
− t
0
)
k
.
Bior
֒
ac dowolne a
0
6= 0, kolejne wsp´o lczynniki a
k
(k = 1, 2, . . .) wyznaczamy z
zale˙zno´sci:
a
k
=
−
P
k
j=1
(p
j
(k
− j + λ
1
) + q
j
) a
k−j
(k + λ
1
) (k + λ
1
− 1) + p
0
(k + λ
1
) + q
0
.
Podobnie, bior
֒
ac dowolne b
0
6= 0, kolejne wsp´o lczynniki b
k
(k = 1, 2, . . .) wyzna-
czamy z zale˙zno´sci:
b
k
=
−
P
k
j=1
(p
j
(k
− j + λ
2
) + q
j
) b
k−j
(k + λ
2
) (k + λ
2
− 1) + p
0
(k + λ
2
) + q
0
.
47
Rozdzia l 6. Rozwi
֒
azania w postaci szereg´
ow funkcyjnych
Gdy λ
1
, λ
2
∈ R, λ
1
= λ
2
, rozwi
֒
azanie szczeg´olne x
2
(t) ma posta´c:
x
2
(t) = x
1
(t) ln (t
− t
0
) + (t
− t
0
)
λ
1
∞
X
k=0
b
k
(t
− t
0
)
k
,
natomiast, gdy λ
1
, λ
2
∈ R, λ
1
≥ λ
2
, λ
1
− λ
2
∈ N jest postaci:
x
2
(t) = Cx
1
(t) ln (t
− t
0
) + (t
− t
0
)
λ
2
∞
X
k=0
b
k
(t
− t
0
)
k
,
gdzie sta la C mo˙ze by´c r´owna zeru.
Podobne wzory mo˙zna wyprowadzi´c dla zespolonych pierwiastl´ow r´ownania
indeksowego.
Rozdzia l 7
Stabilno´
s´
c rozwi
֒
aza´
n r´
owna´
n
r´
o ˙zniczkowych
7.1. Podstawowe definicje
Definicja 26.
Niech X b
֒
edzie przestrzeni
֒
a Banacha. Niech dane b
֒
edzie (RR):
x
′
= f (t, x) z (WPC): x (t
0
) = x
0
, gdzie t
∈ I, I przedzia l , x
0
∈ U ∈ top X
Za l´o˙zmy, ˙ze
∀
y
0
∈U
(RR)
z
(W P C) :
x (t
0
) = y
0
ma rozwi
֒
azanie x (t, y
0
) okre´slone na maksymalnym przedziale istnienia J (y
0
) =
[t
0
, R (t
0
, y
0
)).
1. Rozwi
֒
azanie x (
·, x
0
) nazywamy stabilnym, lub stabilnym w sensie Lapu-
nowa, je˙zeli
∀
ε>0
∃
δ>0
:
ky
0
− x
0
k < δ ⇒ kx (t, y
0
)
− x (t, x
0
)
k < ε
dla t
∈ J (x
0
)
∩ J (y
0
).
2. M´owimy, ˙ze rozwi
֒
azanie x (t, x
0
) jest lokalnie asymptotycznie stabilne, je˙zeli
I = [0, +
∞), rozwi
֒
azanie jest stabilne i ponadto ma w lasno´s´c lokalnego przyci
֒
agania,
tzn.
∃
δ>0
:
ky
0
− x
0
k < δ ⇒
⇒
J (y
0
) = [t
0
, +
∞) , lim
t→∞
kx (t, y
0
)
− x (t, x
0
)
k = 0
.
W skr´ocie piszemy: x (t, x
0
) jest LAS.
49
Rozdzia l 7. Stabilno´s´c rozwi
֒
aza´
n r´
owna´
n r´
o˙zniczkowych
3. M´owimy, ˙ze rozwi
֒
azanie x (t, x
0
) jest globalnie asymptotycznie stabilne, je˙zeli
jest stabilne i ponadto ma w lasno´s´c globalnego przyci
֒
agania, tzn.
∀
y
0
∈U
: J (y
0
) = [t
0
, +
∞) , lim
t→∞
kx (t, y
0
)
− x (t, x
0
)
k = 0.
W skr´ocie piszemy: x (t, x
0
) jest GAS.
Przyk lad 6.
R´ownianie x
′
−x = 0 z warunkiem pocz
֒
atkowym Cauchy’ego x(0) =
x
0
ma rozwi
֒
azanie postaci x(t, x
0
) = x
0
e
t
. Rozwi
֒
azanie x(t, 0) = 0 nie jest
stabilne, bo dla r > 0 mamy
sup
{|x(t, x
0
)
− 0| :
t
≥ 0, |x
0
− 0| < r} = +∞.
Przyk lad 7.
R´ownanie mx
′′
+ 2px
′
+ kx = 0 z warunkiem pocz
֒
atkowym Cau-
chy’ego x(0) = A, x
′
(0) = υ gdy p
2
< km, p > 0, m > 0 ma rozwi
֒
azanie postaci
x(t, A, υ) = Ce
−qt
sin(ωt + ϕ), gdzie q =
p
m
, ω =
q
k
m
− q
2
, C =
q
A
2
+
qA+υ
ω
2
,
ϕ = arccos
A
C
. Rozwi
֒
azanie zerowe jest, co oczywiste, lokalnie asymptotycznie
stabilne.
Twierdzenie 24.
Rozwi
֒
azanie x (t, x
0
) = p (t) r´ownania x
′
= f (t, x) jest sta-
bilne (asymptotycznie stabilne) wtedy i tylko wtedy, gdy rozwi
֒
azanie y (t) = 0
r´ownania y
′
= g (t, y) := f (t, y + p (t))
− f (t, p (t)) jest stabilne (asymptotycznie
stabilne).
Dow´od. Niech p (t) := x (t, x
0
) stabilne rozwi
֒
azanie r´ownania x
′
= f (t, x).
Funkcja y (t) := x (t)
− p (t) spe lnia r´ownanie y
′
(t) = f (t, x (t))
− f (t, p (t)) =
f (t, y (t) + p (t))
− f (t, p (t)) =: g (t, y (t)). Funkcje f i g s
֒
a tej samej klasy
regularno´sci.
Inaczej
. Niech p (t) = x (t, x
0
), x (t) = x (t, y
0
), p
′
= f (t, p), x
′
= f (t, x).
Tak wi
֒
ec x
′
− p
′
= f (t, x)
− f(t, p). Zdefiniujmy z := x − p. Mamy z
′
=
f (t, z + p)
− f(t, p) =: g(t, z) czyli z
′
= g(t, z). Skoro
ky
0
− x
0
k < δ =⇒ kx (t, y
0
)
− x (t, x
0
)
k < ε
zatem
kz
0
k < δ =⇒ kx (t) − p (t)k = kz(t)k < ε
co jest r´ownowa˙zne
kz
0
− 0k < δ =⇒ kz(t) − 0k < ε
Przyk lad 8.
Rozwa˙zmy uk lad r´owna´
n r´o˙zniczkowych
x
′
1
= (x
1
− 1) (x
2
− 1)
x
′
2
= x
1
x
2
− 2
50
7.2. Twierdzenie Lapunowa
kt´ore mo˙zna zapisa´c jako jedno r´ownanie w postaci wektorowej
x
′
= f (x) :=
(x
1
− 1) (x
2
− 1)
x
1
x
2
− 2
gdzie x(t) =
x
1
(t)
x
2
(t)
. Latwo zauwa˙zy´c, ˙ze funkcje x(t) =
1
2
, x(t) =
2
1
s
֒
a rozwi
֒
azaniami przyk ladowego r´ownania (jego po lo˙zeniami r´ownowagi).
Stabilno´s´c pierwszego z tych rozwi
֒
aza´
n jest r´ownowa˙zna stabilno´sci rozwi
֒
azania
zerowego r´ownania
y
′
= f
y +
1
2
− f
1
2
=
y
1
(y
2
+ 1)
y
1
y
2
+ 2y
1
+ y
2
a stabilno´s´c drugiego z nich stabilno´sci rozwi
֒
azania zerowego r´ownania
y
′
= f
y +
2
1
− f
2
1
=
(y
1
+ 1) y
2
y
1
y
2
+ y
1
+ 2y
2
Uwaga 5.
Stabilno´s´c nie implikuje przyci
֒
agania i odwrotnie.
Przyk lad 9.
Rozwa˙zmy r´ownanie x
′′
+ x = 0 z warunkiem pocz
֒
atkowym Cau-
chy’ego x (0) = x
0
, x
′
(0) = 0. Jego rozwi
֒
azaniem jest funkcja x (t) = x
0
cos t.
Rozwi
֒
azanie zerowe jest wi
֒
ec stabilne, ale nie ma w lasno´sci przyci
֒
agania.
Przyk lad 10.
Rozwi
֒
azaniem uk ladu
x
1
x
2
′
=
x
2
−x
1
z warunkiem pocz
֒
atkowym
x
1
x
2
(0) =
x
0
0
jest funkcja
x
1
x
2
(t) =
x
0
cos t
−x
0
sin t
. Tak wi
֒
ec zerowe
rozwi
֒
azanie jest stabilne, ale nie ma w lasno´sci przyci
֒
agania.
7.2. Twierdzenie Lapunowa
Twierdzenie 25.
(Lapunowa) Niech dany b
֒
edzie skalarny uk lad r´owna´
n r´o˙zniczkowych
x
′
j
= f
j
(t, x)
(j = 1, . . . , n)
gdzie t
∈ [t
0
, +
∞), x = x (t) = (x
1
(t) , . . . , x
n
(t))
∈ U ∈ topR
n
, f = (f
1
, . . . , f
n
)
∈
C
1
([t
0
, +
∞) × U, R
n
). Za l´o˙zmy,˙ze
— f (t, 0, . . . , 0) = 0
— a
jk
:=
∂f
j
∂x
k
(t, 0, . . . , 0)
∈ R
j, k = 1, . . . , n
— det
(a
jk
− λδ
jk
)
j,k=1,...,n
= 0 =
⇒ Re λ < 0
51
Rozdzia l 7. Stabilno´s´c rozwi
֒
aza´
n r´
owna´
n r´
o˙zniczkowych
—
∃ M : U −→ R, ˙ze lim
x−→0
M (x) = 0
oraz
f
j
(t, x)
−
n
P
k=1
a
jk
x
k
≤
M (x)
kxk dla t ≥ t
0
, x
∈ U, j = 1, . . . , n.
Wtedy rozwi
֒
azanie zerowe x (
·, 0) = 0 powy˙zszego uk ladu jest lokalnie asympto-
tycznie stabilne tzn. ∃
r>0
ky
0
k < r =⇒ { maksymalny przedzia l J (y
0
) istnienia
rozwi
֒
azania x (
·, y
0
) jest r´owny [t
0
, +
∞) } oraz ∀
ε>0
∃
δ>0
:
ky
0
k < δ =⇒ lim
t−→+∞
x (t, y
0
) = 0 i
kx (t, y
0
)
k < ε dla t ≥ t
0
.
Wniosek 4.
Niech dany b
֒
edzie skalarny uk lad r´owna´
n r´o˙zniczkowych
x
′
j
= f
j
(x)
(j = 1, . . . , n)
gdzie x = x (t) = (x
1
(t) , . . . , x
n
(t))
∈ U ∈ topR
n
, f = (f
1
, . . . , f
n
)
∈ C
1
(U, R
n
).
Za l´o˙zmy, ˙ze
— f (0, . . . , 0) = 0
— a
jk
:=
∂f
j
∂x
k
(0, . . . , 0)
∈ R
j, k = 1, . . . , n
— det
(a
jk
− λδ
jk
)
j,k=1,...,n
= 0 =
⇒ Re λ < 0
Wtedy rozwi
֒
azanie zerowe x (
·, 0) = 0 powy˙zszego uk ladu jest lokalnie asymp-
totycznie stabilne.
Przyk lad 11.
Rozwa˙zmy pierwsze r´ownanie z przyk ladu 8 tj.
y
′
= g (y) :=
y
1
(y
2
+ 1)
y
1
y
2
+ 2y
1
+ y
2
Latwo zauwa˙zy´c, ˙ze g(0) = 0 natomiast
∂g
i
∂x
j
(0)
=
y
2
+ 1
y
1
y
2
+ 2 y
1
+ 1
|y
1
=0, y
2
=0
=
1 0
2 1
. Macierz ta ma warto´s´c w lasn
֒
a λ = 1 o krotno´sci k = 2, a zatem
rozwi
֒
azanie zerowe powy˙zszego r´ownania nie jest lokalnie asymptotycznie stabilne
i nie jest stabilne.
Wniosek 5.
Rozwa˙zmy r´ownanie skalarne x
′
= f (x). Je´sli f (0) = 0 oraz
f
′
(0) < 0, to rozwi
֒
azanie zerowe tego r´ownania jest lokalnie asymptotycznie
stabilne.
Wniosek 6.
Rozwa˙zamy uk lad r´owna´
n r´o˙zniczkowych liniowych x
′
= Ax. Je´sli
wszystkie warto´sci w lasne macierzy A maj
֒
a ujemne cz
֒
e´sci rzeczywiste, to rozwi
֒
azanie
zerowe rozwa˙zanego uk ladu jest lokalnie asymptotycznie stabilne.
52
7.3. Problem Routha–Hurwitza
Twierdzenie 26.
Rozwi
֒
azanie zerowe uk ladu r´owna´
n r´o˙zniczkowych liniowych
x
′
= Ax jest stabilne, gdy Re λ
≤ 0 dla ka˙zdej warto´sci w lasnej λ macierzy A, a
w przypadku Re λ = 0, krotno´s´c tej warto´sci w lasnej jest r´owna 1.
7.3. Problem Routha–Hurwitza
Niech b
֒
edzie dany wielomian W o wsp´o lczynnikach rzeczywistych. Poda´c
takie warunki na jego wsp´o lczynniki, aby pierwiastki wielomianu W le˙za ly w
lewej p´o lp laszczy´znie p laszczyzny zespolonej.
Niech W (λ) = λ
n
+a
1
λ
n−1
+
· · ·+a
n−1
λ+a
n
= 0, gdzie a
j
∈ R (j = 1, . . . , n).
Twierdzenie 27.
(Warunek konieczny)
∀
i∈{1,...,n}
a
i
> 0 . Je˙zeli n
≤ 2, to ten
warunek jest warunkiem wystarczaj
֒
acym.
Twierdzenie 28.
(Warunek Routha–Huwitza) Warunkiem koniecznym i wy-
starczaj
֒
acym na to, aby wszystkie pierwiastki wielomianu W mia ly ujemne cz
֒
e´sci
rzeczywiste jest, aby wszystkie minory g l´owne macierzy Hurwitza
a
1
1
0
0
0
0
· · ·
0
0
0
a
3
a
2
a
1
1
0
0
0
0
0
a
5
a
4
a
3
a
2
a
1
1
0
0
0
...
. ..
...
0
0
0
0
0
0
a
n
a
n−1
a
n−2
0
0
0
0
0
0
· · ·
0
0
a
n
by ly dodatnie.
Twierdzenie 29.
Je´sli W (λ) = a
0
λ
n
+ a
1
λ
n−1
+
· · · + a
n−1
λ + a
n
= 0 jest wielo-
mianem Hurwitza tzn. wszystkie jego pierwiastki maj
֒
a ujemne cz
֒
e´sci rzeczywiste,
to V (λ) := λ
n
W
1
λ
= a
n
λ
n
+ a
n−1
λ
n−1
+ . . . + a
1
λ + a
0
jest tak˙ze wielomianem
Hurwitza.
Przyk lad 12.
Wyznaczy´c obszar asymptotycznej stabilno´sci dla uk ladu
dx
dt
=
−x + αy
dy
dt
= βx
− y + αz
dz
dt
= βy
− z,
gdzie α, β s
֒
a parametrami rzeczywistymi.
53
Rozdzia l 7. Stabilno´s´c rozwi
֒
aza´
n r´
owna´
n r´
o˙zniczkowych
W rozwa˙zanym przypadku wielomian charakterystyczny jest r´owny
W (λ) =
−1 − λ
α
0
β
−1 − λ
α
0
β
−1 − λ
= λ
3
+ 3λ
2
+ (3
− 2αβ)λ + (1 − 2αβ).
Macierz Hurwitza dla tego wielomianu ma posta´c
3
1
0
1
− 2αβ 3 − 2αβ
3
0
0
1
− 2αβ
.
Jej minory g l´owne s
֒
a r´owne:
△
1
= 3,
△
2
= 8
− 4αβ, △
3
= (8
− 4αβ)(1 − 2αβ).
Jak latwo zauwa˙zy´c s
֒
a one wszystkie dodatnie dla αβ <
1
2
.
7.4. Punkty osobliwe r´
ownania r´
o ˙zniczkowego zupe lnego
Rozwa˙zmy r´ownanie
P (t, x)dt + Q(t, x)dx = 0
(7.1)
okre´slone w obszarze ´sci
֒
agalnym D
⊂ R
2
, gdzie P, Q
∈ C
1
(D, R). Poprzednio
zak ladali´smy, ˙ze
|P (t, x)| + |Q(t, x)| > 0. Przy tych za lo˙zeniach mo˙zna by lo
powy˙zsze r´ownanie sprowadzi´c do postaci
x
′
=
−
P (t, x)
Q(t, x)
,
x = x(t),
lub
t
′
=
−
Q(t, x)
P (t, x)
,
t = t(x)
r´owna´
n maj
֒
acych jednoznaczne rozwi
֒
azanie przy zadanych WPC.
Definicja 27.
Je´sli istnieje taki punkt (t
0
, x
0
)
∈ D w kt´orym
P (t
0
, x
0
) = Q (t
0
, x
0
) = 0
to taki punkt nazywamy punktem osobliwym r´ownania r´o˙zniczkowego (7.1).
Przez punkt osobliwy mo˙ze przechodzi´c wiele krzywych ca lkowych, lub ˙zadna
krzywa ca lkowa.
Przyk lad 13.
Rozwi
֒
azaniem og´olnym r´ownania
2t dx
− x dt = 0, (x, t) ∈ R
2
jest rodzina krzywych
t = Cx
2
,
C
∈ R.
Punkt (0, 0) jest punktem osobliwym przez kt´ory przechodzi niesko´
nczenie wiele
ca lek — jest to tzw. punkt w
֒
ez lowy.
54
7.4. Punkty osobliwe r´
ownania r´
o˙zniczkowego zupe lnego
Przyk lad 14.
Rozwi
֒
azaniem og´olnym r´ownania
2at dt + 2bx dx = 0,
(x, t)
∈ R
2
,
a, b > 0
jest rodzina krzywych
at
2
+ bx
2
+ C,
C
∈ R
+
.
Punkt (0, 0) jest punktem osobliwym tzw. punktem wirowym.
Przyk lad 15.
Rozwi
֒
azaniem og´olnym r´ownania
2at dt
− 2bx dx = 0, (x, t) ∈ R
2
,
a, b > 0
jest rodzina krzywych
at
2
− bx
2
+ C,
C
∈ R.
Przez punkt osobliwy (0, 0) przechodz
֒
a dwie krzywe ca lkowe:
x = x(t) =
r
a
b
t,
x = x(t) =
−
r
a
b
t,
t
∈ R.
Punkt (0, 0) jest to tzw. punktem siod lowy.
Przyk lad 16.
R´ownanie
(2t + x) dt + (2x
− t) dx = 0, (x, t) ∈ R
2
0
ma rozwi
֒
azanie og´olne, kt´ore we wsp´o lrz
֒
ednych biegunowych ma posta´c
r = Ce
ϕ/2
C
∈ R
+
.
Punkt osobliwy (0, 0) jest w tym wypadku tzw. punktem asymptotycznym.
Rozdzia l 8
Transformata Laplace’a
8.1. Podstawowe definicje i twierdzenia
Niech ϕ(t) b
֒
edzie funkcj
֒
a zmiennej niezale˙znej t
∈ R. za´s s := σ + iω liczb
֒
a
zespolon
֒
a.
Definicja 28.
Transformat
֒
a Laplace’a (transformat
֒
a) funkcji ϕ(t) nazywamy funkcj
֒
e
ϕ(s) (zmiennej niezale˙znej s) okre´slon
֒
a wzorem
ϕ(s) :=
Z
∞
0
e
−st
ϕ(t)dt
(8.1)
Aby transformata funkcji ϕ(t) by la okre´slona, wystarczy aby ca lka (8.1) istnia la
dla pewnego zbioru warto´sci s, przy czym dla pozosta lych s ca lka ta mo˙ze nie
istnie´c. Mo˙ze si
֒
e zdarzy´c, ˙ze ca lka (8.1) nie istnieje dla ˙zadnej warto´sci s. W tym
przypadku przekszta lcenie Laplace’a nie jest mo˙zliwe.
Przyk lad 17.
Niech ϕ(t)
≡ 1. Latwo policzy´c
ϕ(s) =
R
∞
0
e
−st
dt = lim
A→+∞
R
A
0
e
−st
dt =
=
lim
A→+∞
−
1
s
R
e
−σA
(cos ωA−i sin ωA)
1
z
dz
z
=
= lim
A→+∞
−
1
s
e
−σA
(cos ωA
− i sin ωA) − 1
=
=
1
s
, gdy σ > 0
nie istnieje , gdy σ
≤ 0.
(8.2)
Uwaga 6.
Mo˙zna pokaza´c, ˙ze je´sli ϕ(t) jest w przedziale 0
≤ t ≤ ∞ ograniczona,
albo ro´snie ze wzreostem t jak t
α
lub e
αt
, gdzie α > 0, to jej transformata istnieje.
57
Rozdzia l 8. Transformata Laplace’a
Twierdzenie 30.
(Transformata pochodnej) Niech ψ(t) =
dϕ
dt
. W´owczas ψ(s) =
sϕ
(s)
− ϕ(0), gdzie symbol ϕ(0) oznacza granic
֒
e prawostronn
֒
a w zerze funkcji ϕ.
Dow´
od.
Niech ψ(t) =
dϕ
dt
. Z definicji
ψ(s) =
Z
∞
0
e
−st
ψ(t) dt =
Z
∞
0
e
−st
dϕ
dt
(t) dt =
=
e
−st
ϕ(t)
|
∞
0
+ s
Z
∞
0
e
−st
ϕ(t) dt.
Je´sli
ℜe(s) na tyle du˙ze, ˙ze lim
t→∞
e
−st
ϕ(t) = 0, to ψ(s) = sϕ(s)
− ϕ(0). Je´sli
ϕ(t) ograniczona, lub wzrost ϕ(t) jest wielomianowy (funkcja ro´snie jak t
α
), to
wystarczy przyj
֒
a´c σ > 0, je´sli ϕ(t) ro´snie jak funkcja e
αt
, to wystarczy przyj
֒
a´c
σ > α.
c.k.d.
Ostatni wz´or jest prawdziwy, gdy funkcja ϕ jest ci
֒
ag la; je´sli nie, a konkretnie,
je´sli ma nieci
֒
ag lo´sci skokowe, to we wzorze tym pojawi
֒
a si
֒
e dodatkowe sk ladniki.
W szczeg´olnym przypadku, gdy ϕ(0) = 0 dostajemy ψ(s) = sϕ(s). Otrzymany
rezultat latwo uog´olni´c.
Twierdzenie 31.
Je´sli ψ(t) =
d
n
ϕ
dt
n
, to ψ(s) = s
n
ϕ(s)
− s
n−1
ϕ(0)
− s
n−2
ϕ
′
(0)
−
. . .
−ϕ
(n−1)
(0), gdzie symbol ϕ
(k)
(0) oznacza granic
֒
e prawostronn
֒
a w zerze funkcji
ϕ
(k)
.
Twierdzenie 32.
Je´si ψ(t) :=
R
t
0
ϕ(τ )dτ , to ψ(s) =
ϕ(s)
s
.
Dow´
od.
Zauwa˙zmy, ˙ze ϕ(t) =
dψ
dt
(t),
ψ(0) = 0. Tak wi
֒
ec na podstawie
wzoru na transformat
֒
e pochodnej ϕ(s) = sψ(s), sk
֒
ad bezpo´srednio wynika teza
twierdzenia.
c.k.d.
Twierdzenie 33.
Transformata Laplace’a jest operatorem liniowym.
8.2. Wyznaczanie transformaty r´
ownania r´
o ˙zniczkowego
Niech dane b
֒
edzie r´ownanie r´o˙zniczkowe
x
′
+ ax = f (t),
a
∈ R
z warunkiem pocz
֒
atkowym Cauchy’ego x(0) = x
0
. Mno˙z
֒
ac obie strony r´ownania
przez e
−st
i ca lkuj
֒
ac w granicach od 0 do +
∞ mo˙zemy napisa´c
Z
∞
0
e
−st
x
′
(t) dt + a x(s) = f (s).
Korzystaj
֒
ac ze wzoru na transformat
֒
e pochodnej i warunku pocz
֒
atkowego dosta-
jemy r´ownanie
(s + a)x(s)
− x
0
= f (s),
58
8.3. Wyznaczanie funkcji na podstawie jej transformaty
sk
֒
ad
x(s) =
f (s) + x
0
s + a
,
gdzie f (s) =
R
∞
0
e
−st
f (t) dt.
Analogicznie dla r´ownania
x
′′
+ ax
′
+ bx = f (t),
a, b
∈ R
mno˙z
֒
ac je obustronnie przez e
−st
i ca lkuj
֒
ac w granicach od 0 do
∞ dostajemy
x s
2
+ as + b
= f (s) + s x(0) + x
′
(0) + a x(0),
sk
֒
ad
x =
f (s) + (s + a) x(0) + x
′
(0)
s
2
+ as + b
.
T
֒
e sam
֒
a metod
֒
e mo˙zemy zastosowa´c do uk ladu r´owna´
n o wsp´o lczynnikach sta lych.
Przyk ladowo rozwa˙zmy
x
′
+ a
1
x + b
1
y
′
+ c
1
y = f
1
(t)
x
′
+ a
2
x + b
2
y
′
+ c
2
y = f
2
(t).
Mno˙zymy ka˙zde z tych r´owna´
n przez e
−st
i ca lkujemy w przedziale od 0 do +
∞.
W konsekwencji po przekszta lceniach otrzymujemy uk lad r´owna´
n algebraicznych
s + a
1
b
1
s + c
1
s + a
2
b
2
s + c
2
x(s)
y(s)
=
f
1
(s) + x(0) + b
1
y(0)
f
2
(s) + x(0) + b
2
y(0)
.
Je´sli macierz tego uk ladu jest nieosobliwa, to rozwi
֒
azanie tego uk ladu jest okre´slone
wzorami Cramera.
Transformat
֒
e Laplace’a mo˙zna r´ownie˙z z powodzeniem stosowa´c do pewnych
r´owna´
n r´o˙zniczkowo–ca lkowych np. do r´ownania
x
′
(t) + ax(t) + b
Z
t
0
x(τ ) dτ = f (t).
Og´olnie mo˙zna bez k lopotu poda´c wzory na transformat
֒
e dowolnego r´ownania
liniowego rz
֒
edu n-tego o sta lych wsp´o lczynnikach i uk ladu r´owna´
n r´o˙zniczkowych
o macierzy liczbowej.
8.3. Wyznaczanie funkcji na podstawie jej transformaty
Za l´o˙zmy, ˙ze dana jest funkcja ϕ(s). Zajmiemy si
֒
e problemem wyznaczenia
ϕ(t):
Z
∞
0
e
−st
ϕ(t)dt = ϕ(s)
(8.3)
R´ownanie ca lkowe powy˙zej nazywa si
֒
e r´ownaniem Laplace’a.
59
Rozdzia l 8. Transformata Laplace’a
Twierdzenie 34.
Dla danych a
i
∈ R, ϕ
i
(s), (i = 1, . . . , n) mamy
Z
∞
0
e
−st
n
X
i=1
a
i
ϕ
i
(t)
!
dt =
n
X
i=1
a
i
ϕ
i
(s).
Stosunkowo latwo jest rozwi
֒
aza´c r´ownanie Laplace
֒
a w przypadku, gdy prawa
strona tego r´ownania jest funkcj
֒
a wymiern
֒
a.
Twierdzenie 35.
Je´sli
ϕ(s) =
U (s)
V (s)
gdzie U (s) i V (s) s
֒
a wielomianami, przy czym st.U (s) = m < st.V (s) = n oraz
V (s) = (s
− s
1
) . . . (s
− s
n
), przy czym s
i
6= s
j
je´sli i
6= j, to
ϕ(t) =
n
X
k=1
U (s
k
)
V
′
(s
k
)
e
s
k
t
.
(8.4)
Dow´
od.
Iloraz
U (s)
V (s)
mo˙zna przedstawi´c w postaci sumy u lamk´ow prostych
U (s)
V (s)
=
n
X
k=1
c
k
s
− s
k
.
Mno˙z
֒
ac obie strony przez s
− s
1
mamy
(s
− s
1
)
U (s)
V (s)
= c
1
+ (s
− s
1
)
n
X
k=2
c
k
s
− s
k
.
Przechodz
֒
ac obustronnie z s do granicy w s
1
i stosuj
֒
ac regu l
֒
e de l’Hospitala
dostajemy
U (s
1
)
V
′
(s
1
)
= c
1
.
Podobnie obliczamy warto´sci pozosta lych wsp´o lczynnik´ow c
i
. Tak wi
֒
ec rozk lad
na u lamki proste ma posta´c:
U (s)
V (s)
=
n
X
k=1
U (s
k
)
V
′
(s
k
)
·
1
s
− s
k
.
Latwo si
֒
e przekona´c, ˙ze r´ownanie
Z
∞
0
e
−st
ψ(t)dt =
1
s
− s
k
ma rozwi
֒
azanie
ψ(t) = e
s
k
t
.
Wobec tych fakt´ow i liniowo´sci transformaty otrzymujemy tez
֒
e twierdzenia.
c.k.d.
60
8.3. Wyznaczanie funkcji na podstawie jej transformaty
Twierdzenie 36.
(Twierdzenie o rozk ladzie) Niech
ϕ(s) =
U (s)
s W (s)
,
gdzie U (s) i W (s) s
֒
a wielomianami odpowiednio stopni m i n, przy czym m
≤ n.
Zak ladamy, ˙ze W (0)
6= 0 i wielomian W nie ma pierwiastk´ow wielokrotnych tj.
W (s) = (s
− s
1
) . . . (s
− s
n
), przy czym s
i
6= s
j
dla i
6= j. Wtedy
ϕ(t) =
U (0)
W (0)
+
n
X
i=1
U (s
i
)
s
i
W
′
(s
i
)
e
s
i
t
.
Dow´
od.
Przyjmuj
֒
ac V (s) = s W (s) mo˙zemy na podstawie poprzedniego twier-
dzenia napisa´c
ϕ(t) =
n
X
i=0
U (s
i
)
d
dt
[s W (s)]
|s=s
i
e
s
i
t
,
przy czym s
0
:= 0. Z kolei
d
dt
[s W (s)]
|s=s
i
= W (s
i
) + s
i
W
′
(s
i
) .
Dla i = 0 drugi sk ladnik jest r´owny zeru, a dla i
6= 0 zeruje si
֒
e pierwszy sk ladnik,
tak wi
֒
ec
d
dt
[s W (s)]
|s=s
0
= W (s
0
) = W (0),
d
dt
[s W (s)]
|s=s
i
= s
i
W
′
(s
i
)
dla i
6= 0.
Podstawienie tych wzor´ow do (8.4) ko´
nczy dow´od.
c.k.d.
Twierdzenie 37.
(Twierdzenie o przesuni
֒
eciu rzeczywistym) Niech ϕ(s) b
֒
edzie
transformat
֒
a funkcji ϕ(t), a ψ niech b
֒
edzie funkcj
֒
a zdefiniowan
֒
a wzorem:
ψ(t) :=
0
dla t < t
0
,
ϕ (t
− t
0
) dla t > t
0
.
W´owczas
ψ(s) = e
−st
0
ϕ(s).
Dow´
od.
Z definicji transformaty:
ψ(s) =
Z
∞
0
e
−st
ψ(t)dt =
Z
∞
t
0
e
−st
ϕ (t
− t
0
) dt =
=
Z
∞
0
e
−s(ξ+t
0
)
ϕ(ξ)dξ = e
−st
0
Z
∞
0
e
−sξ
ϕ(ξ)dξ = e
−st
0
ϕ(s).
c.k.d
61
Rozdzia l 8. Transformata Laplace’a
Twierdzenie 38.
(Twierdzenie o przesuni
֒
eciu zespolonym) Niech ψ(t) := e
−λt
ϕ(t),
gdzie λ
∈ R, lub λ ∈ C. Wowczas ψ(s) = ϕ(s + λ).
Dow´
od.
Wprost z definicji:
ψ(s) =
Z
∞
0
e
−st
e
−λt
ϕ(t) dt =
Z
∞
0
e
−(s+λ)t
ϕ(t) dt = ϕ(s + λ).
c.k.d
Twierdzenie 39.
(Twierdzenie o splocie) Niech ψ(t) :=
R
t
0
ϕ
1
(τ )ϕ
2
(t
− τ) dτ.
W´owczas ψ(s) = ϕ
1
(s)ϕ
2
(s).
Obserwacja 1.
Je´sli ψ(t) :=
d
dt
R
t
0
ϕ
1
(τ )ϕ
2
(t
− τ) dτ, to ψ(s) = sϕ
1
(s)ϕ
2
(s).
Rozdzia l 9
Dodatek
9.1. Tablice transformat Laplace’a
Transformat
֒
a Laplace’a (transformat
֒
a) funkcji ϕ(t) nazywamy funkcj
֒
e ϕ(s)
(zmiennej niezale˙znej s
∈ C) okre´slon
֒
a wzorem
ϕ (s) =
Z
∞
0
e
−st
ϕ(t)dt.
Pot
֒
egi
ϕ (t)
ϕ (s)
1
1
s
t
1
s
2
t
n
n!
s
n+1
, n ∈ N
t
−1/2
r
π
s
t
1/2
√
π
2s
3/2
t
α
Γ (α + 1)
s
α+1
, α > −1
63
Rozdzia l 9. Dodatek
Funkcje trygonometryczne
ϕ (t)
ϕ (s)
sin kt
k
s
2
+ k
2
cos kt
s
s
2
+ k
2
sin
2
kt
2k
2
s (s
2
+ 4k
2
)
cos
2
kt
s
2
+ 2k
2
s (s
2
+ 4k
2
)
t sin kt
2ks
(s
2
+ k
2
)
2
t cos kt
s
2
− k
2
(s
2
+ k
2
)
2
2 (1 − cos kt)
t
ln
s
2
+ k
2
s
2
sin at
t
arctan
a
s
ϕ (t)
ϕ (s)
sin kt + kt cos kt
2ks
2
(s
2
+ k
2
)
2
sin kt − kt cos kt
2k
3
(s
2
− k
2
)
2
1 − cos kt
k
2
s (s
2
+ k
2
)
kt − sin kt
k
3
s
2
(s
2
+ k
2
)
a sin bt − b sin at
ab (a
2
− b
2
)
1
(s
2
+ a
2
) (s
2
+ b
2
)
cos bt − cos at
a
2
− b
2
s
(s
2
+ a
2
) (s
2
+ b
2
)
sin at cos bt
t
1
2
arctan
a + b
s
+
1
2
arctan
a − b
s
Funkcje hiperboliczne
ϕ (t)
ϕ (s)
sinh kt
k
s
2
− k
2
cosh kt
s
s
2
− k
2
sinh
2
kt
2k
2
s (s
2
− 4k
2
)
cosh
2
kt
s
2
− 2k
2
s (s
2
− 4k
2
)
ϕ (t)
ϕ (s)
t sinh kt
2ks
(s
2
− k
2
)
2
t cosh kt
s
2
+ k
2
(s
2
− k
2
)
2
2 (1 − cosh kt)
t
ln
s
2
− k
2
s
2
Funkcje wyk ladnicze
ϕ (t)
ϕ (s)
e
at
1
s − a
te
at
1
(s − a)
2
t
n
e
at
n!
(s − a)
n+1
, n ∈ N
e
bt
− e
at
t
ln
s−a
s−b
ϕ (t)
ϕ (s)
1
√
πt
e
−a
2
/4t
e
−a
√
s
√
s
a
2
√
πt
3
e
−a
2
/4t
e
−a
√
s
e
at
− e
bt
a − b
1
(s − a) (s − b)
ae
at
− be
bt
a − b
s
(s − a) (s − b)
64
9.1. Tablice transformat Laplace’a
Funkcje wyk ladnicze i trygonometryczne
ϕ (t)
ϕ (s)
e
at
sin kt
k
(s − a)
2
+ k
2
e
at
cos kt
s − a
(s − a)
2
+ k
2
Funkcje wyk ladnicze i hiperboliczne
ϕ (t)
ϕ (s)
e
at
sinh kt
k
(s − a)
2
− k
2
e
at
cosh kt
s − a
(s − a)
2
− k
2
Funkcje trygonometryczne i hiperboliczne
ϕ (t)
ϕ (s)
sin kt sinh kt
2k
2
s
s
2
+ 4k
4
sin kt cosh kt
k s
2
+ 2k
2
s
4
+ 4k
4
cos kt sinh kt
k s
2
− 2k
2
s
4
+ 4k
4
cos kt cosh kt
s
3
s
4
+ 4k
4
Funkcja Bessela
ϕ (t)
ϕ (s)
J
0
(kt)
1
√
s
2
+ k
2
Uog´
olniona funkcja b l
֒
edu
ϕ (t)
ϕ (s)
erfc
a
2
√
t
= 1 − erf
a
2
√
t
e
−a
√
s
s
2
q
t
π
e
−a
2
/4t
− a erfc
a
2
√
t
e
−a
√
s
s
√
s
e
ab
e
b
2
t
erfc
b
√
t +
a
2
√
t
e
−a
√
s
√
s
√
s + b
−e
ab
e
b
2
t
erfc
b
√
t +
a
2
√
t
+ erfc
a
2
√
t
be
−
a
√
s
s
(
√
s+b
)
Delta Diraca
ϕ (t)
ϕ (s)
δ (t)
1
δ (t − t
0
)
e
−st
0
65
Rozdzia l 9. Dodatek
Funkcja Heaviside’a
ϕ (t)
ϕ (s)
ϕ (t − a) H (t − a)
e
−as
ϕ (s)
H (t − a)
e
−
as
s
przy czym H(t) :=
0
dla
t < 0
1
dla
t ≥ 0
.
Og´
olne prawa
ϕ (t)
ϕ (s)
e
at
ϕ (t)
ϕ (s − a)
ϕ (t − a) H (t − a)
e
−as
ϕ (s)
ϕ
(n)
(t)
s
n
ϕ (s) − s
(n−1)
ϕ (0) − . . . − ϕ
(n−1)
(0)
t
n
ϕ (t)
(−1)
n d
n
ds
n
ϕ (s)
R
t
0
ϕ (τ ) ψ (t − τ) dτ
ϕ (s) ψ (s)
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia
dzienne
Pisemny egzamin z r´owna´
n r´o˙zniczkowych jest dwucz
֒
e´sciowy. Cz
֒
e´s´c pierwsza
ma na celu sprawdzenie bieg lo´sci rachunkowej, a cz
֒
e´s´c druga, umownie zwana
jest cz
֒
e´sci
֒
a ,,teoretyczn
֒
a” i nie ma ona charakteru wy l
֒
acznie rachunkowego. Czas
trwania egzaminu z cz
֒
e´sci zadaniowej: 110 minut. Czas trwania egzaminu z cz
֒
e´sci
teoretycznej: 50 minut. Ka˙zde zadanie jest punktowane w skali 0
− 10 punkt´ow.
Poni˙zej zaprezentowane s
֒
a zestawy zada´
n egzaminacyjnych z kilku sesji. S
֒
a one
reprezentatywne, je´sli chodzi o poziom trudno´sci temat´ow. W poszczeg´olnych
latach zmienia si
֒
e jednak cz
֒
esciowo zakres wyk ladanego materia lu materia lu, a
wi
֒
ec i tematyczny zakres zada´
n.
9 czerwiec 2001
Cz
֒
e´s´c zadaniowa:
1. Rozwi
֒
a˙z r´ownanie Ricattiego
x
′
= 2t
2
+
1
t
x
− 2x
2
wiedz
֒
ac, ˙ze jedn
֒
a z jego ca lek jest wielomian stopnia pierwszego.
2. Wyznacz rozwi
֒
azanie og´olne r´ownania
t
2
(t + 1)x
′′
− 2x = 0
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja x
1
(t) = 1 +
1
t
.
3. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
+ 3x
′
+ 2x = e
−t
cos
2
t.
Wskaz´owka. Tak przekszta l´c praw
֒
a stron
֒
e, aby mo˙zliwe by lo zastosowanie
metody przewidywa´
n.
66
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
4. Metod
֒
a Frobeniusa znajd´z fundamentalny uk lad rozwi
֒
aza´
n r´ownania
2tx
′′
+ (1 + t)x
′
+ x = 0.
5. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
−3
1
2
−4
x +
3t
e
−t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
1
−1
.
6. Zbadaj stabilno´s´c po lo˙ze´
n r´ownowagi uk ladu r´owna´
n:
dx
dt
= y
− x
2
− x
dy
dt
= 3x
− x
2
− y.
Cz
֒
e´s´c teoretyczna:
1. Dla jakich warto´sci paramatr´ow a i b rozwi
֒
azanie zerowe r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne.
2. Znajd´z krzyw
֒
a o tej w lasno´sci, ˙ze trapez utworzony przez osie wsp´o lrz
֒
ednych
Ox i Oy, styczn
֒
a do krzywej i prost
֒
a prostopad l
֒
a do osi Ox w punkcie
styczno´sci, ma sta le pole r´owne 3a
2
.
3. Rozstrzygnij dla jakich a i b rozwi
֒
azania r´ownania x
′′
+ ax
′
+ bx = 0 s
֒
a
ograniczone na ca lej prostej?
20 czerwiec 2001
Cz
֒
e´s´c zadaniowa:
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
3t
2
(1 + ln x) dt =
2x
−
t
3
x
dx
2. Wyznacz rozwi
֒
azanie og´olne r´ownania
tx
′′
− (2t + 1)x
′
+ (t + 1)x = 0
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja postaci e
αt
.
3. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
+ tx
′
− 2t
2
+ 1
x = 0.
67
Rozdzia l 9. Dodatek
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
0
−1 1
0
0 1
−1
0 1
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
1
1
2
1
2
.
5. Wyznacz obszar asymptotycznej stabilno´sci dla uk ladu
dx
dt
=
−x + αy
dy
dt
= βx
− y + αz
dz
dt
= βy
− z,
gdzie α, β s
֒
a parametrami rzeczywistymi.
Cz
֒
e´s´c teoretyczna:
1. Dla jakich warto´sci paramatr´ow a i b rozwi
֒
azanie zerowe r´ownania
x
IV
+ 2x
′′′
+ ax
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne. Zaprezentuj rozwi
֒
azanie w postaci gra-
ficznej.
2. Wyznacz krzywe, dla kt´orych odcinek stycznej zawarty mi
֒
edzy osiami wsp´o lrz
֒
ednych
ma sta l
֒
a d lugo´s´c d.
3. Oblicz e
A
, gdzie A jest macierz
֒
a uk ladu z zadania (4) w cz
֒
e´sci zadaniowej tj.
A =
0
−1 1
0
0 1
−1
0 1
.
13 wrzesie´
n 2001
Cz
֒
e´s´c zadaniowa:
1. Rozwi
֒
a˙z problem pocz
֒
atkowy Cauchy’ego
t
2
+ x
2
dt
− 2tx dx = 0, x(4) = 0.
2. Rozwi
֒
a˙z r´ownanie
t
x
+ 1
dt +
t
x
− 1
dx = 0.
3. Znajd´z ca lk
֒
e og´oln
֒
a r´ownania
x
(6)
+ 2x
(4)
+ x
(2)
= 0.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
5
−1 −4
−12
5
12
10
−3 −9
x
68
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
1
1
1
.
5. Znajd´z uk lad fundamentalny rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych unor-
mowanych w punkcie t
0
= 0 r´ownania:
x
′′
+
1
1
− t
x = 0
i okre´sl rozwi
֒
azanie og´olne.
Cz
֒
e´s´c teoretyczna:
1. Dla jakich warto´sci paramatr´ow a i b rozwi
֒
azanie zerowe r´ownania
x
IV
+ 2x
′′′
+ ax
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne. Zaprezentuj rozwi
֒
azanie w postaci gra-
ficznej.
2. Wyznacz krzywe, dla kt´orych odcinek stycznej zawarty mi
֒
edzy osiami wsp´o lrz
֒
ednych
ma sta l
֒
a d lugo´s´c d.
3. Oblicz e
A
, gdzie A jest macierz
֒
a uk ladu z zadania (4) w cz
֒
e´sci zadaniowej tj.
A =
0
−1 1
0
0 1
−1
0 1
.
27 wrzesie´
n 2001
Cz
֒
e´s´c zadaniowa:
1. Rozwi
֒
a˙z r´ownanie
x
′
= 2
x + 2
t + x
− 1
2
2. Odgadnij rozwi
֒
azanie szczeg´olne, a nast
֒
epnie rozwi
֒
a˙z r´ownanie Riccatiego
x
′
− 2tx + x
2
= 5
− t
2
3. Wiedz
֒
ac, ˙ze funkcja x (t) =
1
t
jest rozwi
֒
azaniem szczeg´olnym r´ownania
2t
2
x
′′
+ 3tx
′
− x = 0 rozwi
֒
a˙z r´ownanie
2t
2
x
′′
+ 3tx
′
− x =
1
t
(obni˙zaj
֒
ac jego rz
֒
ad jednym z dw´och poznanych sposob´ow) a nast
֒
epnie wska˙z
jego ca lk
֒
e spe lniaj
֒
ac
֒
a warunki pocz
֒
atkowe x (1) = 1, x
′
(1) =
−
4
3
.
4. Znajd´z ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n:
x
′
=
−1 2
1
1
x +
2e
t
0
69
Rozdzia l 9. Dodatek
5. Rowi
֒
a˙z r´ownanie
x
′′
+ 3x
′
+ 2x =
1
e
t
+ 1
Cz
֒
e´s´c teoretyczna:
1. Znajd´z krzyw
֒
a o tej w lasno´sci, ˙ze trapez utworzony przez osie uk ladu wsp´o lrz
֒
ednych
Ox, Oy, styczn
֒
a do krzywej i prost
֒
a prostopad l
֒
a do osi Ox w punkcie styczno´sci,
ma sta le pole r´owne 3a
2
.
2. Dla jakich a i b r´ownanie x
′′
+ax
′
+bx = 0 ma przynajmniej jedno rozwi
֒
azanie
x (t)
6= 0 takie, ˙ze lim
t→+∞
x (t) = 0.
3. Zbadaj stabilno´s´c wszystkich po lo˙ze´
n r´ownowagi uk ladu
x
′
= ln (y
2
− x)
y
′
= x
− y − 1
Definicja. Niech X przestrze´
n Banacha, f : X
⊃ U → X, u : R ⊃ I → X,
U
∈ topX, I ∈ topR. Po lo˙zeniem r´ownowagi uk ladu u
′
= f (u) nazywamy
ω
∗
∈ U takie, ˙ze f (ω
∗
) = 0.
10 czerwiec 2002
Cz
֒
e´s´c zadaniowa:
1. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania r´o˙zniczkowego
t(x
′
+ x
2
) = x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(1) = 1.
2. Wyznacz rozwi
֒
azanie og´olne r´ownania
tx
′′
− x
′
− 4t
3
x = 0,
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja e
t
2
.
3. Znajd´z dwa liniowo niezale˙zne rozwi
֒
azania szczeg´olne r´ownania
x
′′
+
2
t
x
′
+ x = 0.
w postaci szereg´ow Frobeniusa w otoczeniu punktu osobliwego regularnego
t
0
= 0.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
−1 −6
3
5
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
2
2
.
5. Wyznacz wszystkie po lo˙zenia r´ownowagi uk ladu
x
′
= xy
y
′
= x
2
+ y
2
− 4
i zbadaj ich stabilno´s´c.
70
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
6. Przy pomocy transformaty Laplace’a rozwi
֒
a˙z r´ownanie
x
′′
− 2x
′
+ x = 1 + t,
x(0) = 0,
x
′
(0) = 0.
Cz
֒
e´s´c teoretyczna:
1. Rozwa˙zamy dwuwymiarowy uk lad r´owna´
n:
x
′
= ax + by
y
′
= cx + dy,
gdzie a, b, c, d
∈ R. Wyka˙z, ˙ze je´sli jedno z jego rozwi
֒
aza´
n jest funkcj
֒
a okre-
sow
֒
a, to wszystkie rozwi
֒
azania, opr´ocz rozwi
֒
azania zerowego, s
֒
a funkcjami
okresowymi.
2. Wyznacz r´ownanie krzywej przechodz
֒
acej przez punkt (1, 1), dla kt´orej pole
tr´ojk
֒
ata utworzonego przez o´s Ot, styczn
֒
a i wektor wodz
֒
acy punktu styczno´sci
jest sta le i r´owna si
֒
e 1.
3. Dla jakich warto´sci paramatr´ow a i b rozwi
֒
azanie zerowe r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne.
17 czerwiec 2002
Cz
֒
e´s´c zadaniowa:
1. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
xdx = (tdx + xdt)
√
1 + x
2
.
2. Rozwi
֒
a˙z r´ownanie
x
′
− 2tx + x
2
= 5
− t
2
.
3. Znajd´z dwa liniowo niezale˙zne rozwi
֒
azania szczeg´olne r´ownania
t(t
− 1)x
′′
+ (1 + t)x
′
− x = 0.
w postaci szereg´ow Frobeniusa w otoczeniu punktu osobliwego regularnego
t
0
= 0, lub t
0
= 1.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
5
3
−3 −1
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
1
−1
.
5. Wyznacz wszystkie po lo˙zenia r´ownowagi uk ladu
x
′
=
−x + y
y
′
= x + y
− 2xy
i zbadaj ich stabilno´s´c.
71
Rozdzia l 9. Dodatek
6. Przy pomocy transformaty Laplace’a wyznacz ca lk
֒
e sczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
−2y + 3t
y
′
= 2x + 4
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 2, y(0) = 3.
Cz
֒
e´s´c teoretyczna:
1. Jakie warunki musz
֒
a spe lnia´c warto´sci i wektory w lasne macierzy uk ladu:
x
′
= ax + by
y
′
= cx + dy,
(a, b, c, d
∈ R), aby jego rozwi
֒
azanie u(t) = (x(t), y(t)) spe lniaj
֒
ace warunek
pocz
֒
atkowy x(0) = y(0) = 1 mia lo w lasno´s´c:
a) lim
t→∞
u(t) = (0, 0),
b) lim
t→∞
ku(t)k = ∞,
c) u jest funkcj
֒
a ograniczon
֒
a.
2. Wyznacz r´ownanie r´o˙zniczkowe rodziny krzywych x = e
Ct
i r´ownanie r´o˙zniczkowe
rodziny krzywych ortogonalnych do danych.
3. Oblicz e
A
dla macierzy:
A =
−2 −4
1
2
.
16 wrzesie´
n 2002
Cz
֒
e´s´c zadaniowa:
1. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
(1 + t + x + tx) x
′
= 1.
2. Rozwi
֒
a˙z r´ownanie
dx = x
2
e
t
− x
dt.
3. Znajd´z dwa liniowo niezale˙zne rozwi
֒
azania szczeg´olne r´ownania
t(t
− 1)x
′′
+ (
−1 + 3t)x
′
+ x = 0.
w postaci szereg´ow Frobeniusa w otoczeniu punktu osobliwego regularnego
t
0
= 0, lub t
0
= 1.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
3 2
−5 1
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
−1
1
.
72
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
5. Wyznacz wszystkie po lo˙zenia r´ownowagi uk ladu
x
′
= 3
−
p
4 + x
2
+ y
y
′
= ln (x
2
− 3)
i zbadaj ich stabilno´s´c.
6. Przy pomocy transformaty Laplace’a wyznacz ca lk
֒
e sczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
−x + y + e
t
y
′
= x
− y + e
t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 1, y(0) = 1.
Cz
֒
e´s´c teoretyczna:
1. Jakie warunki musz
֒
a spe lnia´c warto´sci i wektory w lasne macierzy uk ladu:
x
′
= ax + by
y
′
= cx + dy,
(a, b, c, d
∈ R), aby jego rozwi
֒
azanie u(t) = (x(t), y(t)) spe lniaj
֒
ace warunek
pocz
֒
atkowy x(0) = y(0) = 1 mia lo w lasno´s´c:
a) lim
t→∞
u(t) = (0, 0),
b) lim
t→∞
ku(t)k = ∞,
c) u jest funkcj
֒
a ograniczon
֒
a.
2. Wyznacz r´ownanie r´o˙zniczkowe rodziny hiperbol x =
C
t
i r´ownanie r´o˙zniczkowe
rodziny krzywych ortogonalnych do danych.
3. Oblicz e
A
dla macierzy:
A =
3
−1
2
0
.
9 czerwiec 2003
Cz
֒
e´s´c zadaniowa:
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
6txdt + (4x + 9t
2
)dx = 0.
2. Wyznacz rozwi
֒
azanie og´olne r´ownania
dx
dt
= e
2t
+ (1 + 2e
t
)x + x
2
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja postaci x
1
(t) =
−e
t
.
3. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
+ e
t
x
′
− x = 0.
73
Rozdzia l 9. Dodatek
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
0 8
0
0 0
−2
2 8
−2
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
1
0
0
.
5. Korzystaj
֒
ac z transformaty Laplace’a znajd´z ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
r´o˙zniczkowych
(
d
2
x
dt
2
+
d
2
y
dt
2
= t
2
d
2
x
dt
2
−
d
2
y
dt
2
= 4t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) = 8, x
′
(0) = y(0) = y
′
(0) =
0.
Cz
֒
e´s´c teoretyczna:
1. Dla jakich warto´sci paramatr´ow a i b rozwi
֒
azanie zerowe r´ownania
x
′
= x + ay + y
2
y
′
= bx
− 3y − x
2
.
jest lokalnie asymptotycznie stabilne. Zaprezentuj rozwi
֒
azanie w postaci gra-
ficznej.
2. Znajd´z krzyw
֒
a x = x(t) o tej w lasno´sci, ˙ze tr´ojk
֒
at utworzony przez o´s Ot,
styczn
֒
a do krzywej oraz promie´
n wodz
֒
acy w punkcie styczno´sci jest tr´ojk
֒
atem
r´ownoramiennym.
3. Oblicz e
A
, gdzie A jest macierz
֒
a uk ladu z zadania (4) w cz
֒
e´sci zadaniowej tj.
A =
0 8
0
0 0
−2
2 8
−2
x
16 czerwiec 2003
Cz
֒
e´s´c zadaniowa:
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
(t
2
+ 2tx
− x
2
)dt + (x
2
+ 2tx
− t
2
)dx = 0,
wiedz
֒
ac, ˙ze ma ono czynnik ca lkuj
֒
acy postaci µ = µ(t + x).
2. Rozwi
֒
a˙z r´ownanie
2xx
′
= t(x
′2
+ 4).
3. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
− t
3
x
′
+ (t + 1)x = 0.
74
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
4. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′′
− x
′′
+ 4x
′
− 4x = 3e
2t
− 4 sin 2t.
5. Korzystaj
֒
ac z transformaty Laplace’a rozwi
֒
a˙z r´ownanie
x(t) = 3t
2
− e
−t
−
Z
t
0
x(τ )e
t−τ
dτ .
Cz
֒
e´s´c teoretyczna:
1. Dla jakich paramatr´ow a i b zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne. Zaznacz wyznaczony zbi´or na p laszczy´znie
Oab.
2. Znajd´z rodzin
֒
e krzywych ortogonalnych do krzywych rodziny
x
2
= Ce
t
+ t + 1,
gdzie C
∈ R.
3. Przeprowad´z dyskusj
֒
e dla jakich rzeczywistych parametr´ow p i q wszystkie
rozwi
֒
azania r´ownania x
′′
+ px
′
+ qx = 0 s
֒
a ograniczone na ca lej prostej?
Zaznacz wyznaczony zbi´or na p laszczy´znie Opq.
22 wrzesie´
n 2003
Cz
֒
e´s´c zadaniowa:
1. Znajd´z ca lk
֒
e og´oln
֒
a r´ownania
x
′′
+
1
2
t
−
1
t
x
′
− x = 0,
je´sli x
1
(t) = t
2
jest jego ca lk
֒
a szczeg´oln
֒
a.
2. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
− tx
′
− 2t
2
+ 1
x = 0.
3. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n
x
′
= 2y
− x,
y
′
= 4y
− 3x +
e
3t
e
2t
+1
.
4. Korzystaj
֒
ac z transformaty Laplace’a rozwi
֒
a˙z r´ownanie
x(t) = 3t
2
− e
−t
−
Z
t
0
x(τ )e
t−τ
dτ .
75
Rozdzia l 9. Dodatek
5. Sprowad´z r´ownanie
x
′′
+ 2rx
′
+ ω
2
x = 0
(r > 0, ω > 0)
z warunkami pocz
֒
atkowymi x(0) = x
0
, x
′
(0) = x
1
, do r´ownowa˙znego mu
uk ladu r´owna´
n rz
֒
edu pierwszego oraz zbadaj stabilno´s´c po lo˙zenia r´ownowagi
tego uk ladu.
Cz
֒
e´s´c teoretyczna:
1. Znajd´z krzywe, dla kt´orych tr´ojk
֒
at utworzony przez o´s Oy, styczn
֒
a i wektor
wodz
֒
acy punktu styczno´sci jest r´ownoramienny (o podstawie na osi Oy).
2. Przeprowad´z dyskusj
֒
e dla jakich rzeczywistych parametr´ow p i q wszystkie
rozwi
֒
azania r´ownania x
′′
+ px
′
+ qx = 0 s
֒
a ograniczone na ca lej prostej?
Zaznacz wyznaczony zbi´or na p laszczy´znie Opq.
3. Wyznacz obszar asymptotycznej stabilno´sci dla uk ladu
dx
dt
=
−x + αy
dy
dt
= βx
− y + αz
dz
dt
= βy
− z,
gdzie α, β s
֒
a parametrami rzeczywistymi.
7 czerwiec 2004
Cz
֒
e´s´c pierwsza:
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
(t
2
+ tx + 3x
2
)dt
− (t
2
+ 2tx)dx = 0.
2. Rozwi
֒
a˙z r´ownanie
x
′
=
−
4
t
2
−
1
t
x + x
2
,
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja x(t) =
2
t
.
3. Metod
֒
a Frobeniusa znajd´z fundamentalny uk lad rozwi
֒
aza´
n r´ownania
(1
− t
2
)x
′′
− 2tx
′
+ 30x = 0
w otoczeniu punktu t
0
=
−1 (grupa A), t
0
= 1 (grupa B).
4. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′′
− x
′′
+ 4x
′
− 4x = 3e
2t
− 4 cos 2t.
5. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n
x
′
=
0 1 1
1 0 1
2 2 1
x.
Cz
֒
e´s´c druga:
76
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
1. Zdefiniuj transformat
֒
e Laplace’a funkcji. Oblicz transformat
֒
e Laplace’a ca lki
szczeg´olnej uk ladu r´owna´
n
2x
′
+ y
′
− 2x
= 1
x
′
+ y
′
− 3x − 3y = 2
z warunkiem pocz
֒
atkowym x(0) = 0, y(0) = 0.
2. Obni˙z rz
֒
ad r´ownania
x
′′
− x
′
tan t + 2x = 0
wiedz
֒
ac, ˙ze funkcja x
1
(t) = sin t jest jego ca lk
֒
a szczeg´oln
֒
a.
3. Stosuj
֒
ac twierdzenie Lapunowa zbadaj stabilno´s´c rozwi
֒
azania zerowego uk ladu
r´owna´
n
x
′
= tan(y
− x)
y
′
= 2
y
− 2 cos
π
3
− x
.
15 czerwiec 2004
Cz
֒
e´s´c pierwsza:
1. Rozwi
֒
a˙z zagadnienie pocz
֒
atkowe
tx
2
x
′
+ x
3
= 1,
x(1) = 2.
2. Rozwi
֒
a˙z r´ownanie Lagrange’a
tx
′
(x
′
+ 2) = x.
3. Metod
֒
a Frobeniusa znajd´z fundamentalny uk lad rozwi
֒
aza´
n r´ownania
tx
′′
− (2t − 1)x
′
+ (t
− 1)x = 0
w otoczeniu punktu t
0
= 0.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′′′
− 2x
′′
+ x
′
= te
t
+ 5,
x(0) = 2,
x
′
(0) = 2,
x
′′
(0) =
−1.
5. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n
x
′
=
5
−4 0
1
0 2
0
2 5
x.
Cz
֒
e´s´c druga:
1. Stosuj
֒
ac transformat
֒
e Laplace’a rozwi
֒
a˙z problem pocz
֒
atkowy
x
′′
− 6x
′
+ 9x = t
2
e
3t
,
x(0) = 2,
x
′
(0) = 6.
2. Wiedz
֒
ac, ˙ze jedno z rozwi
֒
aza´
n r´ownania Riccatiego
x
′
− 2tx + x
2
= 5
− t
2
jest wielomianem, sprowad´z to r´ownanie do r´ownania Bernoulliego.
77
Rozdzia l 9. Dodatek
3. Wyznacz warto´sci parametr´ow a i b, dla kt´orych zerowe rozwi
֒
azanie uk ladu
x
′
= x + ay + y
2
y
′
= bx
− 3y − x
2
jest lokalnie asymptotycznie stabilne.
13 wrzesie´
n 2004
Cz
֒
e´s´c pierwsza:
1. Rozwi
֒
a˙z zagadnienie pocz
֒
atkowe
x
′
− 9t
2
x = (t
5
+ t
2
)x
2
3
,
x(0) = 0.
2. Rozwi
֒
a˙z r´ownanie
(t sin x + x cos x) dt + (t cos x
− x sin x) dx = 0.
3. Znajd´z krzyw
֒
a, kt´orej styczne tworz
֒
a z osiami wsp´o lrz
֒
ednych tr´ojk
֒
at o po-
wierzchni 2a
2
.
4. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych o
´srodku w punkcie t
0
= 0 r´ownania:
x
′′
+ tx
′
− (2t
2
+ 1)x = 0.
5. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n
x
′
=
−3
2 2
−3 −1 1
−1
2 0
x.
Cz
֒
e´s´c druga:
1. Stosuj
֒
ac transformat
֒
e Laplace’a rozwi
֒
a˙z problem pocz
֒
atkowy
x
′′
+ 4x
′
+ 13x = te
−t
,
x(0) = 0,
x
′
(0) = 2.
2. Wiedz
֒
ac, ˙ze jedno z rozwi
֒
aza´
n r´ownania Riccatiego
x
′
=
−x
2
+ 1 + t
2
jest wielomianem stopnia pierwszego, sprowad´z to r´ownanie do r´ownania Ber-
noulliego.
3. Stosuj
֒
ac twierdzenie Lapunowa zbadaj stabilno´s´c rozwi
֒
azania zerowego uk ladu
r´owna´
n
x
′
= ln (3e
y
− 2 cos x)
y
′
= 2e
x
−
3
√
8 + 12y.
24 wrzesie´
n 2004
Cz
֒
e´s´c pierwsza:
78
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
1. Rozwi
֒
a˙z r´ownanie
x
′
=
−
4
t
2
−
1
t
x + x
2
,
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja ϕ(t) =
2
t
.
2. Znajd´z krzywe, dla kt´orych odcinek odci
֒
ety na osi rz
֒
ednych Ox (w uk ladzie
wsp´o lrz
֒
ednych Otx) przez styczn
֒
a, jest r´owny kwadratowi rz
֒
ednej punktu
styczno´sci.
3. Znajd´z ca lk
֒
e szczeg´oln
֒
a r´ownania
(t + 1)x
′′
− (2 − t)x
′
+ x = 0
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 2, x
′
(0) =
−1 w postaci szeregu
pot
֒
egowego o ´srodku w punkcie t
0
= 0.
4. Metod
֒
a warto´sci i wektor´ow w lasnych, lub przez sprowadzenie macierzy uk ladu
do postaci Jordana, wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
−1 −2
3
4
x +
3
3
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy
x(0) =
−4
5
.
5. Metod
֒
a transformaty Laplace’a wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
2x
′
+ y
′
− y = t
x
′
+ y
′
= t
2
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 1, y(0) = 0.
Cz
֒
e´s´c druga:
1. Obni˙z rz
֒
ad r´ownania r´o˙zniczkowego
(1 + 2t)x
′′
+ 4tx
′
− 4x = 0
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja ϕ(t) = e
−2t
.
2. Zbadaj dla jakich parametr´ow a i b, zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne.
3. Wyznacz r´ownanie r´o˙zniczkowe rodziny parabol x(t) = at
2
+ bt.
9 czerwiec 2005
Cz
֒
e´s´c pierwsza:
1. Wyznacz krzyw
֒
a le˙z
֒
ac
֒
a w pierwszej ´cwiartce uk ladu wsp´o lrz
֒
ednych, dla kt´orej
styczne tworz
֒
a z osiami uk ladu tr´ojk
֒
aty o sta lym polu r´ownym 2.
79
Rozdzia l 9. Dodatek
2. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania r´o˙zniczkowego
x
′
=
2t
2
x + x
3
2t
3
− tx
2
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(1) = 1.
3. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania
(2t
− t
2
)x
′′
+ (t
2
− 2)x
′
+ 2(1
− t)x = 0
spe lniaj
֒
ac
֒
a warunki pocz
֒
atkowe x(1) = 0, x
′
(1) = 1, je´sli dana jest ca lka
szczeg´olna tego r´ownania x
1
(t) = e
t
.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu
x
′
=
2
−1
1
0
x +
0
2e
t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
1
2
.
5. Wyznacz wszystkie po lo˙zenia r´ownowagi uk ladu
x
′
1
= x
1
x
2
x
′
2
= x
2
1
+ x
2
2
− 4
i zbadaj ich stabilno´s´c.
(Po lo˙zeniem r´ownowagi uk ladu autonomicznego x
′
= f (x), f : R
n
→ R
n
,
x : R
⊃ I → R
n
nazywamy wektor x
0
∈ R
n
taki, ˙ze f (x
0
) = 0.)
Cz
֒
e´s´c druga:
1. Podaj definicj
֒
e transformaty Laplace’a i stosuj
֒
ac j
֒
a, rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
x
′
− x = sin t, x(0) = 1.
2. Metod
֒
a szereg´ow pot
֒
egowych rozwi
֒
a˙z r´ownanie x
′
= x
2
+ t
2
, x(0) = 1. Wy-
znacz cztery pierwsze wyrazy rozwini
֒
ecia.
3. Rozwa˙zamy dwuwymiarowy uk lad r´owna´
n:
x
′
= ax + by
y
′
= cx + dy,
gdzie a, b, c, d
∈ R. Wyka˙z, ˙ze je´sli jedno z jego rozwi
֒
aza´
n jest funkcj
֒
a okre-
sow
֒
a, to wszystkie rozwi
֒
azania, opr´ocz rozwi
֒
azania zerowego, s
֒
a funkcjami
okresowymi.
13 czerwiec 2005
Cz
֒
e´s´c pierwsza:
1. Wyznacz krzyw
֒
a, dla kt´orej rz
֒
edna punktu przeci
֒
ecia dowolnej stycznej z
osi
֒
a Ox (w uk ladzie Otx) jest o dwie jednostki mniejsza od odci
֒
etej punktu
styczno´sci.
80
9.2. Przyk ladowe tematy zada´
n egzaminacyjnych – studia dzienne
2. Rozwi
֒
a˙z r´ownanie
(tx
− t
2
)x
′
+ x
2
− 3tx − 2t
2
= 0
znajduj
֒
ac czynnik ca lkuj
֒
acy zale˙zny od jednej zmiennej.
3. Znajd´z ca lk
֒
e og´oln
֒
a niejednorodnego liniowego r´ownania r´o˙zniczkowego
tx
′′
+ 2x
′
− tx = e
t
wiedz
֒
ac, ˙ze ca lk
֒
a szczeg´oln
֒
a skojarzonego r´ownania jednorodnego jest funkcja
x
1
(t) =
e
t
t
.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu
dx
dt
+
dy
dt
−4y = 1, x+
dy
dt
−3y = t
2
spe lniaj
֒
ac
֒
a
warunek pocz
֒
atkowy Cauchy’ego x(0) = 2, y(0) =
−2.
5. Wyznacz wszystkie po lo˙zenia r´ownowagi uk ladu
x
′
1
= x
2
− x
2
1
− x
1
x
′
2
= 3x
1
− x
2
1
− x
2
i zbadaj ich stabilno´s´c.
(Po lo˙zeniem r´ownowagi uk ladu autonomicznego x
′
= f (x), f : R
n
→ R
n
,
x : R
⊃ I → R
n
nazywamy wektor x
0
∈ R
n
taki, ˙ze f (x
0
) = 0.)
Cz
֒
e´s´c druga:
1. Metod
֒
a szereg´ow pot
֒
egowych wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′′
+ t
2
x
′
+ x = t,
x(0) = 1,
x
′
(0) = 0.
Wyznacz pi
֒
e´c pierwszych wyraz´ow rozwini
֒
ecia.
2. Dla jakich a
∈ R r´ownanie x
′′
+x
′
+ax = 0 ma przynajmniej jedno rozwi
֒
azanie
x(t)
6= 0 takie, ˙ze lim
t→+∞
x(t) = 0.
3. Podaj definicj
֒
e transformaty Laplace’a i stosuj
֒
ac j
֒
a, rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
x
′
+ 4x = 1, x(0) = 1.
12 wrzesie´
n 2005
Cz
֒
e´s´c pierwsza:
1. Wyznacz krzyw
֒
a, kt´orej styczne odcinaj
֒
a na osiach wsp´o lrz
֒
ednych odcinki,
kt´orych suma d lugo´sci jest r´owna 2a.
2. Rozwi
֒
a˙z r´ownanie
dx
dt
= 1
− t − x + tx
2
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja x
1
(t) = 1.
3. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych o
´srodku w punkcie t
0
= 0 r´ownania x
′′
+
1
1−t
x = 0.
4. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu
x
′
= 2y
− x
y
′
= 4y
− 3x +
e
3t
e
2t
+1
.
81
Rozdzia l 9. Dodatek
5. Wyznacz wszystkie po lo˙zenia r´ownowagi uk ladu
x
′
1
= e
x
2
− e
x
1
x
′
2
=
p
3x
1
+ x
2
2
− 2
i zbadaj ich stabilno´s´c.
(Po lo˙zeniem r´ownowagi uk ladu autonomicznego x
′
= f (x), f : R
n
→ R
n
,
x : R
⊃ I → R
n
nazywamy wektor x
0
∈ R
n
taki, ˙ze f (x
0
) = 0.)
Cz
֒
e´s´c druga:
1. Ze wzoru Liouville’a wyznacz drug
֒
a liniowo niezale˙zn
֒
a ca lk
֒
e r´ownania
t
2
(t + 1)x
′′
− 2x = 0
wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja x
1
(t) = 1 +
1
t
.
2. Dla jakich a
∈ R r´ownanie x
′′
+x
′
+ax = 0 ma przynajmniej jedno rozwi
֒
azanie
x(t)
6= 0 takie, ˙ze lim
t→+∞
x(t) = 0.
3. Podaj definicj
֒
e transformaty Laplace’a i stosuj
֒
ac j
֒
a, rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
x
′
+ 2x = t, x(0) =
−1.
22 wrzesie´
n 2005
Cz
֒
e´s´c pierwsza:
1. W pierwszej ´cwiartce uk ladu wsp´o lrz
֒
ednych Otx znajd´z r´ownanie krzywej
przechodz
֒
acej przez punkt (1, 1), dla kt´orej pole tr´ojk
֒
ata utworzonego przez
o´s Ot, styczn
֒
a i wektor wodz
֒
acy punktu styczno´sci jest sta le i r´owna si
֒
e 1.
2. Stosuj
֒
ac podstawienie x(t) = z
2
(t) sprowad´z r´ownanie
t
3
(x
′
− t) = x
2
do r´ownania jednorodnego i rozwi
֒
a˙z go.
3. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych o
´srodku w punkcie t
0
= 0 r´ownania x
′′
+ (1
− t)x = 0.
4. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu
x
′
= x
− y − z
y
′
= x + y
z
′
= 3x + z.
5. Z pomoc
֒
a twierdzenia Lapunowa zbadaj stabilno´s´c zerowego rozwi
֒
azania uk ladu
r´owna´
n
x
′
1
= ln(3e
y
− 2 cos x)
x
′
2
= 2e
x
−
3
√
8 + 12y.
Cz
֒
e´s´c druga:
1. Zredukuj r´ownanie Ricattiego
x
′
=
−
4
t
2
−
1
t
x + x
2
do r´ownania Bernoulliego, wiedz
֒
ac, ˙ze jego ca lk
֒
a szczeg´oln
֒
a jest funkcja x
1
(t) =
2
t
.
82
9.3. Przyk ladowe tematy zada´
n egzaminacyjnych – studia zaoczne
2. Wiedz
֒
ac ˙ze rozwi
֒
azaniem og´olnym uk ladu
x
′
=
1 3
5 3
x
jest
x(t) = c
1
1
−1
e
−2t
+ c
2
3
5
e
6t
,
gdzie c
1
, c
2
∈ R, wyznacz rozwi
֒
azanie og´olne uk ladu
x
′
=
1 3
5 3
x +
1
0
.
3. Opisz metod
֒
e redukcji rz
֒
edu jednorodnego r´ownania r´o˙zniczkowego liniowego
o sta lych wsp´o lczynnikach, je´sli znana jest jego ca lka szczeg´olna.
9.3. Przyk ladowe tematy zada´
n egzaminacyjnych – studia
zaoczne
Pisemny egzamin z r´owna´
n r´o˙zniczkowych jest jednocz
֒
e´sciowy. Czas trwania
egzaminu: zazwyczaj 120 minut. Ka˙zde zadanie jest punktowane w skali 0
− 10
punkt´ow. Poni˙zej zaprezentowane s
֒
a zestawy zada´
n egzaminacyjnych z kilku
sesji.
16 czerwiec 2002
1. Rozwi
֒
a˙z r´ownanie jednorodne
x
′
=
x
t + x
2. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania Bernoulliego
t(x
′
+ x
2
) = x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(1) = 1.
3. Wyznacz 4 pierwsze wyrazy rozwini
֒
ecia rozwi
֒
azania problemu pocz
֒
atkowego
x(x
′
+ 1) = t,
x(0) = 1
w szereg pot
֒
egowy w otoczeniu punktu t
0
= 0.
4. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
+ 3x
′
+ 2x = t.
5. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
−1 −6
3
5
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
2
2
.
83
Rozdzia l 9. Dodatek
6. Dla jakich warto´sci paramatr´ow a i b rozwi
֒
azanie zerowe r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne. Zaznacz wyznaczony zbi´or na p laszczy´znie
Oab.
20 wrzesie´
n 2002
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
dx
dt
+
x
t
=
−tx
2
.
2. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
+ x
′
+ tx = 0.
3. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
0
−1
3
4
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
−1
1
.
4. Korzystaj
֒
ac z transformaty Laplace’a znajd´z ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′′
− x
′
= sin t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) = 1, x
′
(0) = 0.
5. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
− 2x
′
+ x = t + e
t
.
6. Dla jakich paramatr´ow a i b zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ 2x
′
+ bx = 0
jest lokalnie asymptotycznie stabilne. Zaznacz wyznaczony zbi´or na p laszczy´znie
Oab.
28 wrzesie´
n 2002
1. Rozwi
֒
a˙z r´ownanie
1 + e
t
xx
′
= e
t
.
2. Rozwi
֒
a˙z r´ownanie
x
′
=
tx
− x
2
t
2
.
84
9.3. Przyk ladowe tematy zada´
n egzaminacyjnych – studia zaoczne
3. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
− x = t
2
− t + 1.
4. Znajd´z ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
1 1
4 1
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy postaci
x(1) =
1
0
.
5. Dla jakich paramatr´ow a i b zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne. Zaznacz wyznaczony zbi´or na p laszczy´znie
Oab.
14 czerwiec 2003
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
6txdt + (4x + 9t
2
)dx = 0.
2. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
+ tx
′
− x = 0.
3. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
0
−3
1
2
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
1
0
.
4. Korzystaj
֒
ac z transformaty Laplace’a znajd´z ca lk
֒
e szczeg´oln
֒
a r´ownania
dx
dt
− 3x = e
2t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) = 1.
5. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
+ 2x
′
+ x = te
t
.
85
Rozdzia l 9. Dodatek
6. Dla jakich paramatr´ow a i b zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne. Zaznacz wyznaczony zbi´or na p laszczy´znie
Oab.
20 wrzesie´
n 2003
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe
dx
dt
+
x
t
=
−tx
2
.
2. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
+ x
′
+ tx = 0.
3. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
0
−1
3
4
x
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) =
−1
1
.
4. Korzystaj
֒
ac z transformaty Laplace’a znajd´z ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′′
− x
′
= sin t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy Cauchy’ego x(0) = 1, x
′
(0) = 0.
5. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
− 2x
′
+ x = t + e
t
.
6. Dla jakich paramatr´ow a i b zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ 2x
′
+ bx = 0
jest lokalnie asymptotycznie stabilne. Zaznacz wyznaczony zbi´or na p laszczy´znie
Oab.
3 pa´
zdziernik 2003
1. Znajd´z ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′
=
−e
x+t+1
,
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) =
−1.
86
9.3. Przyk ladowe tematy zada´
n egzaminacyjnych – studia zaoczne
2. Znajd´z fundamentalny uk lad rozwi
֒
aza´
n w postaci szereg´ow pot
֒
egowych, unor-
mowany w punkcie t
0
= 0, r´ownania
x
′′
− x
′
+ tx = 0.
3. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n
x
′
=
0
−1
2
2
x.
4. Korzystaj
֒
ac z transformaty Laplace’a rozwi
֒
a˙z zagadnienie pocz
֒
atkowe
x
′′
+ x = t,
x(0) = 0,
x
′
(0) = 1.
5. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
− 2x
′
+ x = e
−t
.
6. Dla jakich paramatr´ow a i b zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ bx
′
+ x = 0
jest lokalnie asymptotycznie stabilne. Zaznacz wyznaczony zbi´or na p laszczy´znie
Oab.
12 czerwiec 2004
1. Rozwi
֒
a˙z jednorodne r´ownanie r´o˙zniczkowe
(t
2
+ tx + 3x
2
)dt
− (t
2
+ 2tx)dx = 0.
2. Znajd´z ca lk
֒
e szczeg´oln
֒
a w postaci szeregu pot
֒
egowego, unormowanego w
punkcie t
0
= 0, r´ownania
x
′′
+ tx
′
− x = 0, x(0) = 1, x
′
(0) = 0.
3. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′′
− x
′′
+ 4x
′
− 4x = 3e
2t
.
4. Wyznacz ca lk
֒
e szczeg´oln
֒
a uk ladu r´owna´
n
x
′
=
3
−3
2
−2
x +
4
−1
z warunkiem pocz
֒
atkowym x(0) =
0
0
.
87
Rozdzia l 9. Dodatek
5. Stosuj
֒
ac transformat
֒
e Laplace’a oblicz x-sow
֒
a sk ladow
֒
a ca lki szczeg´olnej uk ladu
r´owna´
n
2x
′
+ y
′
− 2x
= 1
x
′
+ y
′
− 3x − 3y = 2
z warunkiem pocz
֒
atkowym x(0) = 0, y(0) = 0.
6. Dla jakiej warto´sci parametru a zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ x
′′
+ 2x
′
+ x = 0
jest lokalnie asymptotycznie stabilne.
17 wrzesie´
n 2004
1. Rozwi
֒
a˙z problem pocz
֒
atkowy Cauchy’ego
x
′
sin t = x ln x,
x(π/2) = 1.
2. Znajd´z ca lk
֒
e szczeg´oln
֒
a w postaci szeregu pot
֒
egowego, unormowanego w
punkcie t
0
= 0, r´ownania
x
′′
+ x
′
− t
2
x = 0,
x(0) =
−1, x
′
(0) = 1.
3. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′′
− 2x
′
+ 2x = te
−t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = x
′
(0) = 0.
4. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n
x
′
= 2x
− y,
y
′
= x + 2e
t
.
5. Stosuj
֒
ac transformat
֒
e Laplace’a oblicz y-kow
֒
a sk ladow
֒
a ca lki szczeg´olnej uk ladu
r´owna´
n
2x
′
+ y
′
− 2x
= 1
x
′
+ y
′
− 3x − 3y = 2
z warunkiem pocz
֒
atkowym x(0) = 0, y(0) = 0.
6. Dla jakich warto´sci parametr´ow a i b, zerowe rozwi
֒
azanie r´ownania
x
IV
+ 2x
′′′
+ 4x
′′
+ ax
′
+ bx = 0
jest lokalnie asymptotycznie stabilne.
27 wrzesie´
n 2004
1. Rozwi
֒
a˙z problem pocz
֒
atkowy Cauchy’ego
x
′
= x ln x,
x(0) = e.
88
9.3. Przyk ladowe tematy zada´
n egzaminacyjnych – studia zaoczne
2. Znajd´z ca lk
֒
e szczeg´oln
֒
a w postaci szeregu pot
֒
egowego, unormowanego w
punkcie t
0
= 0, r´ownania
x
′′
+ (t + 1)x = 0,
x(0) = 1,
x
′
(0) = 1.
3. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′′
− 6x
′
+ 9x = t
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 0, x
′
(0) = 1.
4. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu r´owna´
n
x
′
=
x + 2y,
y
′
=
−
1
2
x + y.
5. Dla jakich warto´sci parametr´ow a i b, zerowe rozwi
֒
azanie r´ownania
x
IV
+ 2x
′′′
+ 4x
′′
+ ax
′
+ bx = 0
jest lokalnie asymptotycznie stabilne.
6. Stosuj
֒
ac transformat
֒
e Laplace’a rozwi
֒
a˙z problem pocz
֒
atkowy Cauchy’ego
x
′′
+ x = sin t,
x(0) = 1,
x
′
(0) =
−1.
11 czerwiec 2005
1. Rozwi
֒
a˙z r´ownanie
tx
x
′
+ x
2
= 1.
2. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania
x
′
−
x
1
− t
2
= t + 1
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 0.
3. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
(2t
− t
2
)x
′′
+ (t
2
− 2)x
′
+ 2(1
− t)x = 0,
je´sli dana jest ca lka szczeg´olna tego r´ownania x
1
(t) = e
t
.
4. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu
x
′
=
2
−1
1
0
x +
0
2e
t
.
5. Metod
֒
a szereg´ow pot
֒
egowych wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania r´o˙zniczkowego
x
′′
+ tx
′
− x = 0
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 1, x
′
(0) = 0. Wyznacz cztery pierw-
sze niezerowe wyrazy szeregu.
89
Rozdzia l 9. Dodatek
6. Zbadaj dla jakich parametr´ow a i b zerowe rozwi
֒
azanie r´ownania
x
IV
+ ax
′′′
+ 4x
′′
+ 2x
′
+ bx = 0
jest lokalnie asymptotycznie stabilne.
21 wrzesie´
n 2005
1. Rozwi
֒
a˙z jednorodne r´ownanie r´o˙zniczkowe
2t
3
x
′
= x(2t
2
− x
2
).
2. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe zupe lne
x
t
dt + (x
2
+ ln t) dx = 0.
3. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
− 5x
′
= 3t
2
.
4. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu
x
′
=
4 1
−2 1
x +
−e
2t
0
.
5. Metod
֒
a szereg´ow pot
֒
egowych wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania r´o˙zniczkowego
x
′′
− (t + 1)x = 0
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 1, x
′
(0) = 1. Wyznacz cztery pierw-
sze niezerowe wyrazy szeregu.
6. Zbadaj dla jakiego parametru a zerowe rozwi
֒
azanie r´ownania
x
IV
+ 2x
′′′
+ 3x
′′
+ 2x
′
+ ax = 0
jest lokalnie asymptotycznie stabilne.
24 wrzesie´
n 2005
1. Rozwi
֒
a˙z r´ownanie r´o˙zniczkowe o zmiennych rozdzielonych
t(1 + x
2
)dt + x(1 + t
2
)dx = 0.
2. Wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania Bernoulliego
x
′
− 2tx = 2t
3
x
2
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 1.
90
9.3. Przyk ladowe tematy zada´
n egzaminacyjnych – studia zaoczne
3. Wyznacz ca lk
֒
e og´oln
֒
a r´ownania
x
′′
− 5x
′
= e
t
.
4. Wyznacz ca lk
֒
e og´oln
֒
a uk ladu
x
′
=
3
−2
1
1
x.
5. Metod
֒
a szereg´ow pot
֒
egowych wyznacz ca lk
֒
e szczeg´oln
֒
a r´ownania r´o˙zniczkowego
x
′′
+ tx = 0
spe lniaj
֒
ac
֒
a warunek pocz
֒
atkowy x(0) = 1, x
′
(0) = 2. Wyznacz cztery pierw-
sze niezerowe wyrazy szeregu.
6. Zbadaj dla jakiego parametru a zerowe rozwi
֒
azanie r´ownania
x
IV
+ 2x
′′′
+ 3x
′′
+ ax
′
+ x = 0
jest lokalnie asymptotycznie stabilne.
Bibliografia
[1] F.Bierski, Funkcje zespolone, Szeregi i przekszta lcenia Fouriera, Przekszta lcenia
ca lkowe Laplace’a, Przekszta lcenia Laurenta (Z), wyd. pi
֒
ate poprawione, Uczel-
niane Wydawnictwa Naukowo-Dydaktyczne AGH, Krak´
ow 1999.
[2] B.P.Conrad, Differential Equations, A Systems Approach, Pearson Education,
Inc., Upper Saddle River, New Jersey, 2003.
[3] B.P.Demidowicz, Matematyczna teoria stabilno´sci, Wyd. Naukowo-Techniczne,
Warszawa 1972.
[4] L.Dru˙zkowski, Analiza Matematyczna dla fizykow, Cz
֒
e´s´c II, Wybrane zagadnie-
nia, Wyd. UJ, Krak´
ow 1997.
[5] A.F.Filippow, Zbi´
or zada´
n z r´
owna´
n r´
o˙zniczkowych, Izd. Nauka, Moskwa 1973.
[6] I.M. Gelfand, Wyk lady z algebry liniowej, wyd. 3, PWN, Warszawa 1977.
[7] R.Gutowski, Rownania r´
o˙zniczkowe zwyczajne, Wyd. Naukowo-Techniczne, War-
szawa 1971.
[8] M.I.Kontorowicz, Rachunek operatorowy i procesy w uk ladach elektrycznych,
Wyd. Naukowo-Techniczne, Warszawa 1968.
[9] N.M.Matwiejew, Metody ca lkowania r´
owna´
n r´
o˙zniczkowych zwyczajnych, PWN,
Warszawa 1972.
[10] J.Niedoba, W.Niedoba, R´
ownania r´
o˙zniczkowe zwyczajne i cz
֒
astkowe, Zadania z
matematyki, Wydanie trzecie, Uczelniane Wydawnictwa Naukowo-Dydaktyczne
AGH, Krak´
ow 2001.
[11] J.Ombach, Wyk lady z r´
owna´
n r´
o˙zniczkowych, Wyd. UJ, Krak´
ow 1996.
[12] A.Palczewski, R´
ownania r´
o˙zniczkowe zwyczajne (teoria i metody numeryczne
z wykorzystaniem komputerowego systemu oblicze´
n symbolicznych), Wyd.
Naukowo-Techniczne, Warszawa 1999.
[13] A.Pelczar, J.Szarski, Wst
֒
ep do teorii r´
owna´
n r´
o˙zniczkowych, Cz
֒
e´s´c I, PWN, War-
szawa 1987.
[14] A.Pelczar, Wst
֒
ep do teorii r´
owna´
n r´
o˙zniczkowych, Cz
֒
e´s´c II, PWN, Warszawa
1989.
[15] K.K.Ponomariew, Uk ladanie i rozwi
֒
azywanie r´
owna´
n r´
o˙zniczkowych w zagadnie-
niach technicznych, Wyd. Naukowo-Techniczne, Warszawa 1969.
[16] W.Stankiewicz, J.Wojtowicz, Zadania z matematyki dla wy˙zszych uczelni tech-
nicznych, Cz
֒
e´s´c II, PWN, Warszawa 1983.
93
Bibliografia
[17] F.G.Tricomi,Differential Equations, Blackie&Son Limited, 1961.
[18] D.G.Zill, Differential Equations with Boundary-Value Problems, PWS-KENT Pu-
blishing Company, Boston, 1986.
94