2. Розв’язування систем лiнiйних рiвнянь.
Нехай маємо систему лiнiйних рiвнянь
|
|
|
Axt = bt, |
|
де |
a21 |
|
|
|
|
|
|
a22 |
. . . a2n |
|
x = (x1, x2, . . . , xn), |
|
a11 |
a12 |
. . . a1n |
|
|
|
A = |
... |
... |
... ... |
, |
b = (b1, b2, . . . , bn). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
an1 |
an2 |
. . . ann |
|
|
|
|
|
|
|
|
|
|
Перепишемо його у виглядi Axt + xt −bt = xt або A xt −bt = xt, де
|
|
a21 |
a22 + 1 . . . |
a2n |
|
|
|
a11 + 1 a12 |
. . . |
a1n |
. |
A= |
... |
... ... ... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
an1 |
an2 |
. . . ann + 1 |
|
|
|
|
|
|
|
|
Будемо вважати, що для будь-яких x = (x1, x2, . . . , xn), y = (y1, y2, . . . , yn) Rn
d(x, y) = max |xk − yk|.
k=1,n
Розглянемо вiдображення f : Rn → Rn, яке дiє за таким прави-
лом: |
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
t t |
. |
|
|
x = (x1, x2, . . . , xn) f(x) = ( A x |
− b ) |
|
Тодi для будь-яких x = (x1, x2, . . . , xn), y = (y1, y2, . . . , yn) |
|
d(f(x), f(y)) = d |
( A x − b ) |
, ( A y |
− b ) |
= |
|
|
|
|
|
t |
t t |
t |
t |
t |
|
|
= d |
=1 a1j xj |
− b1 |
, j=1 |
a2j xj − b2, . . . , j=1 anj xj − bn! |
, |
|
n |
|
n |
|
|
|
n |
|
|
|
Xj |
|
X |
|
|
X |
|
|
=1 a1j yj − b1, j=1 |
a2j yj − b2, . . . , j=1 anj yj − bn!! = |
|
n |
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
n |
|
|
|
|
Xj |
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
X |
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
X |
|
|
|
|
|
j)| 6 i=1,n |
Xj |
|
|
|
|
|
|
|
j| 6 |
= i=1,n |
a |
ij ( |
j − |
|
a |
ij |
|
|| |
j − |
|
| j=1 |
y |
=1 | |
|
y |
max |
|
|
|
|
x |
|
|
|
|
|
max |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Xj |
|
|
|
X |
|
|
|
|
|
max x |
|
|
|
y |
|
|
a |
|
= max |
|
a |
|
d(x, y). |
|
− |
max |
=1 | |
ij | |
|
|
|
| |
ij | |
6 i=1,n |
| |
j |
|
|
j| i=1,n |
|
|
i=1,n j=1 |
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отже, якщо max |
jP |
|
|
|
|
, то вiдображення f є стискуючим |
a |
|
|
|
< |
|
=1 | |
|
| |
|
|
i=1,n |
|
|
ij |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
i рiвняння f(x) = x має єдиний розв’язок, тобто задана система рiвнянь має єдиний розв’язок, який можна шукати методом послiдовних наближень.
3. Розв’язування диференцiальних рiвнянь першого поряд-
ку.
Нехай у диференцiального рiвняння
y0 = f(x, y)
функцiя f(x, y) визначена i неперервна на прямокутнику D = {(x, y) | |x − x0| 6 a, |y − y0| 6 b} i задовольняє на ньому умови Лiпшиця по змiннiй y, тобто iснує L > 0 таке, що для будь-яких точок (x, y1), (x, y2) D виконується нерiвнiсть
|
|f(x, y1) − f(x, y2)| 6 L|y1 − y2| |
|
i нехай h |
= min{ |
a, |
b |
, |
1 |
} |
, де M |
= |
max |
| |
f |
x, y |
)| |
. |
|
|
|
M |
|
L |
|
(x,y) D |
( |
|
|
Позначимо через X множину всiх неперервних на вiдрiзку [x0 −h; x0 + h] функцiй, кожна з яких набуває значень з вiдрiзка
[y0 − b; y0 + b], тобто X = {f | f C[x0−h;x0+h], |f(x) − y0| 6 b}.
Ця множина, як замкнена пiдмножина повного метричного простору C[x0−h;x0+h] з рiвномiрною метрикою, є повним метричним простором з цiєю ж метрикою.
|
|
|
Побудуємо вiдображення F множини X |
у множину |
C[x0−h;x0+h] у такийx спосiб: кожнiй функцiї y(x) |
X вiднесе- |
мо функцiю y0 + Z |
f(t, y(t))dt, яка є неперервною на вiдрiзку |
x0
[x0 − h; x0 + h] як сума сталої та iнтеграла зi змiнною верхнею межею вiд неперервної функцiї. Крiм того, оскiльки
|
x |
|
x |
|
|y0 + xZ0 |
f(t, y(t))dt − y0| = |xZ0 |
f(t, y(t))dt| 6 |
x |
|
x |
|
|
6 |xZ0 |
|f(t, y(t))|dt| 6 M|xZ0 |
dt| 6 Mh 6 b, |
то образ кожної функцiї y(x) з X належить множинi X, тобто вiдображення F є вiдображенням X в себе. А оскiльки для будь-яких y1(x), y2(x) X
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
d(F (y1(x)), F (y2(x))) = d(y0 + xZ0 |
f(t, y1(t))dt, y0 + |
|
|
x |
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
)) − |
|
+xZ0 |
f |
( |
t, y |
2( )) ) = x [x0−h;x0+h] |xZ0 |
| |
f |
( |
t, y |
1( |
|
|
|
|
|
t dt |
max |
|
|
|
|
|
|
|
t |
|
|
|
− |
|
|
|
|
|
|
))| | 6 |
x |
| |
|
|
|
|
) − |
|
2( )| |
|
6 |
f |
( |
t, y |
2 |
( |
x [x0−h;x0+h]xZ0 |
y |
1( |
y |
dt |
|
|
|
t |
|
dt |
L max |
|
|
|
t |
|
|
|
|
t |
|
6 |
|
x0+h |
−h;x0+h] | |
|
|
( ) − |
|
|
( |
)| |
|
= |
L |
xZ0 |
x [x0 |
y |
1 |
y |
2 |
dt |
|
|
max |
|
t |
|
t |
|
|
x0+h
Z
= L d(y1(x), y2(x))dt = L · hd(y1(x), yy(x)),
x0
то, врахувавши, що L · h < 1, маємо, що вiдображення F є стискуючим. Отже, рiвняння
x |
|
y(x) = y0 + xZ0 |
f(t, y(t))dt |
має єдиний розв’язок Y (x), який задовольняє умову Y (x0) = y0. А оскiльки задача вiдшукання розв’язку цього iнтегрального рiвняння еквiвалентна задачi вiдшукання розв’язку диференцiального рiвняння y0 = f(x, y), причому цей розв’язок має задовольняти початкову умову Y (x0) = y0, то задачу Кошi для останнього рiвняння можна розв’язувати методом послiдовних наближень.
Як iллюстрацiю розглянемо диференцiальне рiвняння
y0 = y + ln x
i будемо шукати розв’язок, який задовольняє початкову умову Y (1) = 1, тобто будемо шукати iнтегральну криву, що проходить через точку M(1, 1). Легко переконатись, що
x
y = ex(C + Z1 |
e−t ln t dt) |
загальний розв’язок цього |
рiвняння. (Зауважимо, що |
Z |
e |
−t |
ln t dt |
у |
скiнченному виглядi |
через |
елементарнi |
фун- |
|
|
1 |
|
|
|
кцiї не подається). При C = e− |
|
маємо iнтегральну криву |
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
y = ex−1 + ex Z1 |
e−t ln t dt, |
|
|
яка проходить через точку M0(1, 1).
У прямокутнику D = {(x, y) | |x − 1| 6 12 , |y − 1| 6 12 } функцiя f(x, y) = y + ln x неперервна i для всiх (x, y) D |f(x, y)| 6 2.
Для будь-яких (x, y1), (x, y2) D |f(x, y1) − f(x, y2)| = |y1 − y2|, тобто ця функцiя задовольняє умови Лiпшиця по змiннiй y з константою L = 1. Оскiльки
min 12, 14, 1 = 14 < 1,
то iснує єдина функцiя y = Y (x), визначена i неперервна разом з своєю похiдною на вiдрiзку [34 ; 54 ], яка є розв’язком заданого рiвняння i задовольняє умову Y (1) = 1. Застосувавши метод послiдовних наближень, маємо:
y0(x) = 1,
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1(x) = 1 + Z1 |
(1 + ln t)dt = 1 + x ln x, |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y2(x) = 1 + Z1 |
|
|
|
|
|
5 |
|
x2 |
x2 |
(1 + ln t + t ln t)dt = |
|
|
|
|
+ x ln x − |
|
ln x − |
|
, |
4 |
2 |
4 |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y3(x) = 1 + Z1 |
5 |
|
t2 |
t2 |
|
|
|
|
( |
|
+ t ln t + |
|
ln t − |
|
|
|
+ ln t)dt = |
|
|
4 |
2 |
|
|
4 |
|
|
41 |
|
x |
x2 |
5 |
x3 + x ln x + |
|
x2 |
|
|
|
|
x3 |
|
|
|
|
|
|
|
|
|
|
= |
|
+ |
|
− |
|
− |
|
|
|
|
ln x + |
|
|
|
|
ln x. |
|
|
|
|
|
|
|
36 |
4 |
4 |
36 |
2 |
6 |
|
|
|
|
|
|
|
|
Врахувавши, що для всiх x [x0 − h; x0 + h] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
M |
|
|
∞ |
(Lh)k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|Y (x) − yn(x)| 6 |
|
|
kX |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k! , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L |
=n+1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
у нашому випадку маємо: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|Y (x) − y3(x)| 6 2 k=4 |
|
= 2 |
√4 e − 1 |
− |
|
|
− |
|
|
|
|
− |
|
|
|
|
|
< |
4kk! |
4 |
4 |
· |
2! |
4 |
· |
3! |
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< |
3 · 10−4. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Завдання для самоконтролю.
1.Нехай на множинi X задано двi метрики d1 i d2. Будемо називати їх еквiвалентними, якщо для будь-якої послi-
довностi (xn) lim d1(xn, x0) = 0 тодi i тiльки тодi, коли
n→∞
lim d2(xn, x0) = 0. Довести, що якщо iснують α1 > 0 i
n→∞
α2 > 0 такi, що
α1d1(x, y) 6 d2(x, y) 6 α2d1(x, y)
для будь-яких x, y X, то метрики d1 i d2 еквiвалентнi, пручому з повноти одного з просторiв (X, d1), (X, d2) випливає повнота другого.
2. Довести, що метричнi простори (Rn, di), де
d1 |
(x, y) = v |
n |
(xk |
− |
yk)2 |
, |
|
uk=1 |
|
|
|
|
uX |
|
|
|
|
t
|
|
n |
|
|
|
|
|
|
Xk |
|
− yk|, |
d2(x, y) = |
|xk |
|
|
=1 |
|
|
|
|
d |
, |
max x |
k − |
y |
k| |
|
3(x y) = k=1,n |
| |
|
повнi, а їх метрики еквiвалентнi.
3.Довести, що у повному метричному просторi кожна послiдовнiсть вкладених одна в одну замкнених куль, радiуси яких прямують до нуля, має перетином одноелементну множину.
4.Знайдiть нерухомi точки вiдображення f
а) R2 → R2 за правилом y1 = 2x2 + 3x22, y2 = 2x1 + 3x21; б) R3 → R3 за правилом y1 = x2x3, y2 = 2x1x3,
y3 = 3x1x2;
в) R4 → R4 за правилом y1 = 1 − x2 − x3 − x4,
y2 = 4x1 + 3x3 + x4 − 4, y3 = 2 − 2x1 + 3x2 + x4, y4 = x4.
5. Якi з вiдображень
а) f(x) = x2 на [0; 13];
б) |
u = 0, 7x + 0, 8y, v = 0, 2x + 0, 05y на R2; |
|
|
|
x |
|
в) |
F (y) = |
1 |
Z0 |
y(x)dx на C[0;2] з рiвномiрною |
|
3 |
метрикою є стискуючим?
6. Переконайтесь, що система
|
x + 5y |
z = |
8, |
|
10x − 2y + |
z = |
9, |
− |
|
|
|
4x + 2y + 8z = 32 |
|
|
|
|
має єдиний розв’язок, i знайдiть його з точнiстю до 0,01 методом послiдовних наближень, взявши за нулеве наближення точку (0,0,0).
7.Методом послiдовних наближень розв’язати функцiональнi (iнтегральнi) рiвняння:
|
|
|
x |
|
а) |
f(x) = 1 + Z0 |
f(t)dt, f0(x) ≡ 0; |
|
|
|
x |
|
б) |
f(x) = x − Z0 |
(x − t)f(t)dt, f0(x) ≡ 0; |
|
|
x |
f(t)dt + ex − 2(e − 1), f0(x) ≡ 0. |
в) f(x) = 2 Z0 |
|
1 |
|
|
1 |
27 ЛЕКЦIЯ: Числовi ряди з дiйсними та компле-
ксними членами
Числовий ряд, його збiжнiсть. Геометрична прогресiя та гармонiйний ряд. Властивостi збiжних рядiв. Ознаки збiжностi рядiв з невiд’ємними членами. Абсолютно ти умовно збiжнi ряди та їх властивостi.
Лiтература. [1] ч. 1, с. 272–275, 283–292; [3] т. 1, с. 477–514; [9], ч. 2, с. 3–25.
Нехай маємо послiдовнiсть (zn) комплексних чисел. Для кожного n iснує число
n
X z1 + z2 + · · · + zn = zk,
k=1
яке є сумою n перших членiв заданої послiдовностi, тобто послiдовнiсть (zn) породжує нову послiдовнiсть (Sn), де
n
X
Sn = zk.
k=1
Застосування до послiдовностi (Sn) основної операцiї аналiзу дає можливiсть поняття суми узагальнити на деякi випадки „сум“ нескiнченного числа доданкiв. Таким узагальненням є поняття ряду i його суми.
Означення 27.1. Символ виду
|
∞ |
|
z1 + z2 + · · · + zn + · · · або |
X |
(27.1) |
zn |
n=1
називається числовим рядом (або просто рядом), а числа z1, z2, . . . , zn, . . . називають членами ряду.
|
n |
|
|
Послiдовнiсть (Sn), де Sn |
Xk |
zk, називають послiдовнiстю |
= |
|
=1 |
|
|
часткових сум ряду (27.1), а ряд |
|
|
|
|
|
∞ |
zn+1 + zn+2 + · · · |
+ zn+k + · · · або |
kX |
zn |
=n+1
називають залишком ряду (27.1) пiсля n-го члена.
Означення 27.2. Якщо iснує
lim Sn = S,
n→∞
то ряд (27.1) називається збiжним, а число S називається його сумою i записують
∞
X
zn = S.
n=1
У противному разi ряд (27.1) називається розбiжним.
|
∞ |
|
X |
Приклад 1. Дослiдити на збiжнiсть ряд |
aqn−1, де a, q |
|
n=1 |
C (геометрична прогресiя з першим членом a i знаменником q). Розв’язання. Оскiльки для заданого ряду
Sn = a + aq + · · · + aqn−1,
а qSn = aq + aq2 + · · · + aqn, то, вiднявши почлено цi рiвностi, |
маємо |
|
|
|
|
a − aqn |
|
|
|
|
|
|
Sn = |
. |
|
|
|
Якщо |q| < 1, то nlim qn = 0 i |
|
1 − q |
|
|
|
|
|
|
|
|
|
|
|
→∞ |
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
aqn |
a |
lim S |
|
= lim |
|
|
|
|
|
|
|
|
. |
|
|
− q |
− 1 − q = |
|
n→∞ |
n |
n→∞ 1 |
1 − q |