Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Matanaliz

.pdf
Скачиваний:
17
Добавлен:
02.06.2015
Размер:
1.56 Mб
Скачать

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

 

 

401

=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}.

402

Ця множина, як замкнена п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

 

403

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вняння. (Зауважимо, що

404

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

 

 

405

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

406

 

 

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

метрикою є стискуючим?

407

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

408

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, . . . називають членами ряду.

409

 

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

410

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]