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

Дифференциальные уравнения

..pdf
Скачиваний:
12
Добавлен:
05.02.2023
Размер:
474.71 Кб
Скачать

101

Пусть A Rn Rn — линейный оператор. Тогда Ax + b = 0 — система n линейных уравнений с n неизвестными. Рассмотрим оператор B Rn Rn, действующий по формуле Bx = Ax + b. При b 0 B — оператор, полученный из линейного оператора A сдвигом на вектор b. При b = 0 оператор B совпадает с оператором A. Найдем условия сжимаемости оператора B при различных метриках в пространстве Rn. Для Rn1 имеем:

n

 

 

 

ρi(Bx, By) = ρi

(Ax + b,n

Ayn+ bi) j=

 

 

 

 

n

 

i j

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

n

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

= ∑ (

 

 

 

− ( ) −

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

i 1

 

 

Ax

 

 

 

bi

 

 

 

 

Ay

 

 

 

 

 

 

bi

 

 

i

1 Rj

 

1 ajx

 

 

 

 

 

j

1 ajy R

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

R

 

=

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

a x

 

 

 

y

 

 

 

 

 

 

 

 

 

Ra

 

 

 

x

 

 

 

 

y

 

 

 

 

 

 

 

 

R

 

 

 

 

 

R

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

R

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

n

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

i

 

 

j

 

 

 

 

 

 

j

 

R

 

 

 

 

 

 

 

j

 

 

 

 

 

j

=

 

 

 

 

 

= ∑R∑

 

 

(

 

 

 

 

 

)R

∑∑

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

R

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

=

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

x

 

 

 

y

 

 

 

 

 

R

 

 

max

 

 

 

 

 

 

 

aj

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

i 1

Raj

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rj 1

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

i

1 j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ∑(∑

R

i

)

 

 

j

 

 

j

 

 

 

R

 

 

(∑

 

i

)

 

 

j

 

j

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

1 i 1 R

 

 

 

i

 

 

 

 

 

 

 

 

 

j

 

Rj 1 1 j n

 

i

 

1

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

=

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

=j

 

 

 

 

 

=

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

max

 

 

 

a

j j 1

 

x

y

=

max

 

 

 

 

 

a

j

ρ(

x, y .

 

1 j n i 1

 

 

 

 

 

1 j n i 1

 

 

 

 

 

)

 

 

 

 

 

 

=

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, мы получили условие сжимаемости оператора B, а следова-

тельно, и оператора A.

n

maxaij < 1.

1 j n i=1

Для Rnимеет вид

условие сжимаемости оператора B, а следовательно, и оператора A

n

maxaij < 1.

1 i n j=1

Для Rn2 условие сжимаемости оператора B, а следовательно, и оператора A

nn

∑∑(aij)2 < 1.

i=1 j=1

Соответствующие вычисления предлагается проделать самостоятельно или посмотреть в [9].

Подводя итоги, получаем, что если систему n линейных уравнений с n неизвестными удается записать в форме x = Ax + b с матрицей A, удовлетворяющей одному из полученных условий сжимаемости оператора A, то, по теореме о сжимающем операторе, последовательные приближения xn+1 = Axn + b сходятся к точке x0, являющейся решением данной системы линейных уравнений. Соответствующий процесс называется итерационным.

На этой идее основаны методы простой итерации и его модификации (метод Зейделя).

Пусть теперь функция f (x, y) удовлетворяет условиям теоремы существования и единственности (см. п. 1.4) решения задачи Коши для дифференциального уравнения y= f (x, y) с начальными условиями y(x0) = y0, то есть непрерывна по совокупности переменных в некоторой области D и удовлетворяет в ней условию Липшица по y. Перейдём к эквивалентному интегральному уравнению

102

 

 

 

 

 

 

Приложениe Б. Принцип сжатых отображений

 

=

 

+

x

(

 

)

 

(

 

) =

y

y0

f

x, y

dx. Рассмотрим оператор, действующий по формуле By

x

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

= y0 +

x

 

 

 

 

 

 

 

 

f (t, y(t)) dt. Этот оператор переводит непрерывную функцию в непрерыв-

x0

ную. Получим условия сжимаемости оператора B в метрике пространства C[a, b]. Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rx0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

(

 

) =

[

 

]

 

 

(

 

) −x

 

 

 

( )

=

[

 

 

R

 

 

 

(

 

 

(

 

))−

 

(

 

 

 

 

(

)))

 

R

 

 

By , By

 

y

 

x

 

 

 

 

] R (

 

 

 

 

f

t, y

 

dt

R

 

 

max

 

 

 

y

 

 

 

x

 

 

 

 

max

R∫

 

f

t, y

 

 

t

 

 

 

 

 

 

 

t

 

R

 

 

ρ 1 2

x a, b

 

 

1

 

 

 

 

 

2

 

 

 

 

 

x a, b

 

R

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

max

 

 

f

 

t, y

1

 

t

 

 

 

Rf

t, y

 

t

))

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

x x [a, b]x0

 

(

 

 

(

 

)) −R (

 

 

2(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

x [a, b]x0

 

 

 

(

 

 

 

 

(

)

 

 

 

 

 

 

0 x [a, b]

 

 

 

(

 

) −

 

2(

 

 

)

 

 

 

 

 

 

 

1

) −

 

 

2

dt

 

x

y

1

x

y

x

.

 

 

 

 

 

 

max

L y

 

t

 

y

 

t

 

 

 

L x

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, если L x

 

x0

 

 

 

 

1, то оператор B — сжимающий. Тогда, по теоре-

 

 

 

 

 

 

 

 

 

 

решение интегрального уравнения y

=

y

 

x

f

(

x, y

dx,

 

 

 

 

 

 

 

 

 

 

 

 

ме о сжимающем операторе, − <

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 + x0

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а следовательно, и задачи Коши y= f (x, y), y(x0) = y0, существует и единственно

11

на интервале (x0 L, x0 + L).

Приложение В

ТАБЛИЦА ИНТЕГРАЛОВ

 

 

 

0αdx

=

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 dx

 

 

 

+

C.

 

 

 

 

 

1.

 

 

 

 

 

xα

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

 

 

dx

= x

 

 

 

 

 

 

 

 

 

3.

x dx

 

 

 

 

 

 

 

+

 

 

 

 

 

 

C, α

 

 

 

 

 

1.

 

 

 

 

 

 

 

 

 

 

4.

 

 

 

 

 

 

ln x

 

 

 

C.

 

 

 

 

= α

 

 

1 +

 

≠ −

 

 

 

 

 

 

 

 

 

 

 

x =

 

 

 

 

 

 

 

 

 

 

dx

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ̃ x

 

 

 

 

 

+

 

 

 

 

 

 

5.

 

 

 

 

 

=

 

 

 

 

 

 

 

+

 

x = −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

x2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5a.

 

 

+

 

 

 

 

 

 

 

= a

 

 

 

 

 

 

 

 

 

 

a +

 

 

 

 

= −a

 

 

 

 

a +

̃

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a dx x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg

 

 

 

C

 

 

 

 

 

arcctg

 

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

+

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x̃

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.

 

 

 

 

 

 

 

 

 

 

 

 

 

arcsin x

 

 

C

 

 

 

 

 

arccos x

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

+

+

= −

 

 

 

 

 

+

+

̃

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xa

 

 

 

 

 

ax =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6a.

 

 

 

 

 

 

 

 

 

 

 

 

arcsin

 

 

 

 

C

 

 

 

arccos

 

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

 

a dx

 

 

ln a

 

 

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7а.

 

 

e dx

 

 

ex

 

 

 

C.

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

+

 

 

 

 

 

 

8.

 

 

cos x dx

=

 

sin x

 

+

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9.

 

sin x dx

= −

cos x

+

C.

10.

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11.

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg x

+

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

ctg x

+

C.

 

 

cos2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin2x

12.

 

sh x dx=

 

 

 

 

 

 

 

 

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13.

 

ch x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

ch x

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

=

sh x

+

C.

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

14.

 

 

 

 

 

 

 

= −

cth x

 

+

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15.

 

 

 

 

 

x

=

th x

+

C.

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

ax

 

 

 

 

 

 

 

 

 

 

 

 

eax

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16.

 

eax cos bx dx =

 

(b sin bx + a cos bx) + C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2eax b2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

17.

 

e

 

sin bx dx

 

 

 

 

 

 

 

 

+

 

 

 

 

a sin bx

 

b cos bx

 

C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= a

2

+ b

2

(

) +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учебное издание

Ельцов Александр Александрович Ельцова Тамара Александровна

ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

Учебное пособие

Корректор Осипова Е. А. Компьютерная верстка Перминова М. Ю.

Подписано 09.10.13. Формат 60х84/8. Усл. печ. л. 12,09. Тираж 500 экз. Заказ

Издано в ООО «Эль Контент» 634029, г. Томск, ул. Кузнецова д. 11 оф. 17

Отпечатано в Томском государственном университете систем управления и радиоэлектроники.

634050, г. Томск, пр. Ленина, 40 Тел. (3822) 533018.