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

smagina-belousova-method

.pdf
Скачиваний:
31
Добавлен:
21.05.2015
Размер:
491.17 Кб
Скачать

 

 

1

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

d)

xn =

 

,...,

 

,0,0,... , α > 0 .

 

 

 

 

nα

nα

 

 

 

 

 

14243

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

3. Сходятся ли в C[0,1],

C1[0,1] последовательности

 

 

 

 

a) xn (t) =

t n+1

 

t n+2

;

b) xn (t) = sin t sin(

t

) .

 

 

 

n +1

n + 2

n

 

 

 

 

 

 

 

 

 

 

4.Сходится ли в L1[0,1] последовательность xn (t) = 1 + nt 2t 2 ?

5.Привести пример последовательности, которая принадлежала бы каждой из рассматриваемой пары пространств и:

a)сходилась бы в m , но не сходилась в l1 ;

b)с сходилась бы в m , но не сходилась в l2 ;

c)ходилась бы в l2 , но не сходилась в l1 .

3. Открытые и замкнутые множества

Пусть

X

- нормированное пространство. Открытым шаром с

центром

в

 

 

точке

x0

и

радиуса

r > 0

называется

множество

B(x0 , r) = {x X :

 

 

 

 

 

x x0

 

 

 

 

 

< r}. Замкнутым шаром с центром в точке x0

и

 

 

 

 

радиуса

r > 0

 

 

 

 

 

называется

множество

B[x0 , r] = {

 

x x0

 

 

 

r}.

Сферой

с

 

 

 

 

 

 

 

 

центром

в

 

 

точке

x0

и

радиуса

r > 0

называется

множество

S(x0 , r) = {x X :

 

 

 

x x0

 

 

 

= r}. Под окрестностью точки x0

понимается любой

 

 

 

 

открытый шар с центром в этой точке.

 

 

 

 

 

 

 

 

 

 

 

Множество

 

 

 

 

 

M X

называется ограниченным,

если

его можно

заключить в некоторый шар (открытый или замкнутый) конечного радиуса. Множество M называется открытым, если любая его точка является внутренней, т.е. для любого x0 M существует такое r > 0 , что

B(x0 , r) M . Точка x0 X называется предельной точкой множества M , если существует последовательность точек xn M , отличных от x0 , такая что xn x0 0 при n → ∞. Множество M называется замкнутым, если оно содержит все свои предельные точки.

Свойства открытых и замкнутых множеств

1. Пространство X открыто и замкнуто одновременно.

11

2.Множество M X замкнуто тогда и только тогда, когда его дополнение CM = X \ M открыто.

3.Пересечение конечного числа и объединение любого числа открытых множеств открыто.

4.Объединение конечного числа и пересечение любого числа замкнутых множеств замкнуто.

Пример 1. Описать B(x0 , r) C[a,b] .

 

 

 

 

 

 

 

 

 

 

 

Решение. Пусть x B(x0 , r) . Это

означает,

что

 

 

 

x x0

 

 

 

C < r

или

 

 

 

 

 

 

 

 

 

 

max

 

x(t) x0 (t)

 

< r . Отсюда r < x(t) x0 (t) < r для

любого t [a,b]

или

 

 

atb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0 (t) r < x(t) < x0 (t) + r . Таким образом,

шар B(x0 , r)

в пространстве

C[a,b]

представляет собой полосу шириной 2r , получающуюся сдвигом графика

функции x0 (t)

вверх и вниз на r .

 

 

 

 

 

Пример 2.

Пусть M = {x = {xk } l2 : xk > 0 k}.

Будет ли оно открытым,

замкнутым в l2 ?

1

 

 

1

 

Решение.

Покажем, что точка x0 = (

,...,

,...) M не является

2

2n

 

 

 

 

 

внутренней, т.е. в любой ее окрестности содержатся точки, не

принадлежащие M . Действительно,

пусть задано произвольное ε > 0 .

Выберем n так, чтобы

 

 

 

1

 

 

 

< ε 2 ,

2k

k =n+1

2

 

 

что возможно в силу стремления к нулю остатка сходящегося ряда. Тогда

для xn = (

1

,

1

,...,

1

,0,0,...) имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

xn x0

 

 

 

2 =

 

< ε 2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2

k =n+1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. xn B(x0 ,ε) ,

но

xn M , т.е

множество M

не является открытым.

Заметим, что M совпадает с пересечением бесконечного числа открытых

множеств M n

= {x = {xk } l2 : xn > 0}, где n фиксировано.

Покажем,

что M

не

является

 

 

и

замкнутым. Для этого рассмотрим

последовательность точек

xn = (

1

,

1

,

1

 

 

,...) M . Так как

2n

 

 

2n+1

2n+2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

 

 

 

2

=

 

1

=

1

1

=

1

0

 

 

 

 

 

 

 

 

2(n+k 1)

n

k

n1

 

 

 

 

 

 

 

 

l2

k =1 2

 

4

k =0 4

3 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при n → ∞, то x0 = {0} - предельная точка множества M , но x0 M .

Пример 3.

Показать,

что

M = {x(t) C[0,1] : sin t < x(t) <1 +t} - открытое

множество в C[0,1] .

 

 

 

 

 

 

 

 

 

 

 

Решение.

Пусть

x0 M .

 

Так

как

функции x0 (t) sin t и

1 +t x0 (t) непрерывны, то они ограничены на [0,1] . Обозначим

 

α = min(x0 (t) sin t) > 0 ,

β = min(1 +t x0 (t)) > 0

 

 

0t1

 

 

 

 

 

 

 

 

0t1

и положим r = min (α, β) . Тогда для x B(x0 , r) , с одной стороны,

 

x(t) < x0 (t) + r x0 (t) + β x0 (t) +1 +t x0 (t) =1 +t ,

с другой стороны,

 

 

 

 

 

 

 

 

 

 

 

 

x(t) > x0 (t) r x0 (t) α x0 (t) x0 (t) +sin t = sin t .

Таким образом, B(x0 , r) M , т.е.

M - открытое множество.

Пример4. Показать, чтомножество M = {x C[0,1] : x(0) = 0} замкнутов C[0,1] .

Решение. Пусть x0 - предельная точка множества M . Тогда существует

последовательность {xn } M

такая, что xn

x0 и

 

 

 

xn x0

 

 

 

C

0 . Отсюда

 

 

 

 

 

 

 

 

 

следует, что x0 C[0,1]

как

равномерный

предел

 

последовательности

непрерывных функций. Так как из сходимости в C[0,1] вытекает

поточечная сходимость,

то xn (t) x0 (t) для любого

 

t [0,1] .

В частности,

при t = 0 имеем x0 (0) = lim xn (0) = 0 . Это означает, что x0 M .

 

n→∞

 

 

 

 

 

 

 

 

 

 

 

 

Пример 5. Пусть P[a,b] - множество всех алгебраических полиномов, определенных на [a,b] . Является ли оно замкнутым, открытым в C[0,1] ?

 

n

t

k

 

Решение. Рассмотрим многочлен

pn (t) =

 

P[a,b]. Так как он

k!

 

k =0

 

является отрезком ряда Тейлора для функции x0 (t) = et , то pn (t) равномерно на [a,b] сходится к x0 (t) , а значит последовательность {pn } сходится к

13

x0 и по норме пространства C[a,b] . Таким образом, x0 - предельная точка

множества P[a,b], но x0 P[a,b] , т.е. множество P[a,b] не является замкнутым в C[0,1] .

Покажем, что множество P[a,b] не является и открытым в C[0,1] . Для этого рассмотрим многочлен первой степени p0 (t) = t P[a,b]. Пусть задано

произвольное ε > 0 . Покажем,

что в шаре

B( p0 ,ε)

содержатся точки, не

являющиеся полиномами. Положим x(t) = t +sin

t

, где

1

< ε . Тогда имеем

 

 

n

 

 

 

 

 

 

 

 

 

1

 

 

1

 

n

 

 

 

 

 

p0 x

 

 

 

C

= max

sin t

 

< ε ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

atb

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. x B( p0 ,ε) , но x P[a,b] .

Задания для самостоятельного решения

1.Построить на числовой прямой:

а) последовательность замкнутых множеств, объединение которых не является замкнутым;

б) последовательность открытых множеств, пересечение которых не является открытым множеством.

2.Изобразить открытый единичный шар с центром в нуле в пространствах: а) R2 , б) R12 , с) R22 .

3.В пространстве C[1,1] даны множество M и точка x0 . Проверить, что x0 является для M внутренней точкой:

 

а) M = {x(t) :

 

x(t)

 

1}, x0 (t) = sin t; x0 (t) =

t

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

б) M = {x(t) : x(t) < t}, x0 (t) = t 1.

 

 

 

 

 

4.

Показать, что

 

множество

M = x C[a,b] :

t

< x(t) <1

+ cos t

открыто

 

4

 

вC[a,b] .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.

Показать,

 

что

множества

 

 

M1 = {x C[a,b] : x(t) 0},

 

M 2 = {x C[a,b] :1 x(t) 2} замкнуты в C[a,b] .

 

 

 

 

6.Доказать, что множество M = {x C[a,b] : x(a) =1} замкнуто в C[a,b] .

7.Будет ли замкнутым в C[a,b] множество многочленов степени точно

k ?

14

4. Принцип сжимающих отображений

Теорема. Пусть выполняются следующие условия:

1)X - банахово пространство, Q - замкнутое множество в нем;

2)Отображение f переводит Q в себя;

3) Отображение f

является сжимающим,

 

 

 

то есть существует 0 < q <1

такое, что для любых x, y Q выполняется неравенство

 

 

 

 

 

 

 

f (x) f ( y)

 

 

 

q

 

 

 

x y

 

 

 

 

 

 

.

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда отображение f

имеет в Q единственную неподвижную точку x* , то

есть такую точку, что

f (x* ) = x* . Эта точка может быть получена методом

последовательных приближений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim xn = x* ,

 

 

где xn = f (xn1 ),

n =1,2,...

 

n→∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

начиная с произвольного

 

 

 

x0 Q .

 

 

 

При

 

 

 

этом

справедлива оценка

погрешности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn x*

 

 

 

 

qn

 

 

 

 

 

f (x0 ) x0

 

 

 

.

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание. При проверке условия сжатия (1) часто бывает удобно пользоваться формулой Лагранжа конечных приращений

f (x) f ( y) = f (ξ)(x y), x ξ y .

(3)

 

Пример 1. Показать, что для вычисления a можно пользоваться формулой

xn+1 =

 

1

(xn +

a

) , начиная с любого x0

a .

 

 

2

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Выберем замкнутое в R множество Q =[

a, ) и рассмотрим

на нем отображение

 

 

 

 

1

 

 

 

a

 

1

 

 

a

 

 

 

f (x) =

 

 

(x

+

 

) . Так

как

f (x) =

 

 

(1

 

) > 0 , то

f

 

 

2

x

2

x2

монотонно возрастает

и, следовательно,

f (x) f ( a ) =

a ,

то есть

f

переводит Q в себя. Кроме того, f

- сжимающее отображение, так как

 

15

 

 

 

 

 

 

f (x) f ( y)

 

=

1

(1

 

a

)

 

x y

 

1

 

 

 

x y

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому существует единственное решение x* Q

 

 

 

уравнения f (x) = x , то

есть уравнения

1

(x +

a

) = x или

x2

= a . Таким образом, последовательность

 

 

2

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn действительно сходится к

a .

При этом справедлива оценка скорости

сходимости

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn a

 

 

 

 

1

 

 

 

x1

x0

 

1

 

 

 

 

x02 a

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

1

 

2n

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

(1

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2. Доказать, что уравнение 2tet =1 имеет единственное решение в интервале (0,1) . Оценить количество приближений для его вычисления с

точностью до 102 .

 

 

 

 

 

 

 

 

 

 

 

Решение. Запишем уравнение в виде t =

et

и рассмотрим отображение

 

 

et

 

 

 

 

 

 

 

 

 

 

 

2

 

 

f (t) =

. Очевидно, что f

монотонно убывает и для t [0,1] выполняются

 

2

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

неравенства

0 <

 

 

= f (1) f (t) f (0) =

<1, т.е.

f отображает отрезок [0,1]

в

2e

2

 

 

 

 

 

 

et

 

 

 

 

 

 

 

себя. Так

как

 

 

 

1

, то

согласно

формуле Лагранжа f

-

 

f (t)

=

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сжимающее отображение и, следовательно, имеет вид на [0,1] единственную неподвижную точку, т.е. данное уравнение имеет единственное решение t* [0,1]. Так как t = 0, t =1 не удовлетворяет уравнению, то t* (0,1) . Это решение может быть получено методом последовательных приближений, начиная с любого t0 (0,1) . Справедлива оценка погрешности

 

 

 

tn t*

 

1

 

t1

t0

 

=

 

1

 

et0

t0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2n1

2n1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частности, для

t0 = 0 имеем

 

 

tn

t*

 

2n <102 . Отсюда 2n >102 , или

 

 

n lg 2 2 , т.е. n >[

lg 2

] = 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

Задания для самостоятельного решения

 

1.Доказать,

 

что

последовательность

цепных

дробей

2, 2 +

1

, 2 +

1

 

 

,... имеет предел и найти его.

 

 

 

 

1

 

 

 

2

 

2 +

 

 

 

 

 

 

 

 

 

 

2

 

 

любых a [0,1] последовательные приближения

2. Доказать,

что

для

xn+1 = xn (xn2

a) / 2,

x0 = 0

сходятся к

a .

 

 

3. При каком a > 0 последовательность

 

 

 

 

 

 

 

 

 

 

 

a,

a + a , a +

a + a ,...

 

 

имеет предел? Найдите его.

4. Доказать, что последовательные приближения:

a) tn+1

=

 

tn

, t0 =1 ; b)

tn+1

=

(5 +tn2 )

, t0 = 5

2

+tn

 

 

 

 

 

 

2tn

имеют предел и найти его.

5. Является ли сжимающим отображение

x

 

 

0.7 x

 

+ 0.8x

 

 

 

 

1

 

1

 

2

 

x

2

 

 

0.2x

+ 0.05x

2

 

 

 

 

1

 

 

 

 

плоскости в себя в R2 , R12 , R22 .

6.Показать, что следующие уравнения имеют единственное решение:

a) x =

1

1 + x2 +

1 sin x ; b) t(3 +sin t) =1 ; c)

x =

1

arctg x +

1

cos x .

8

2

4

 

 

4

 

 

 

4. Решение интегральных уравнений Фредгольма и Вольтерры второго рода методом последовательных приближений

Уравнение

x(t) = λb K (t, s)x(s)ds + f (t)

(1)

a

 

17

называется интегральным уравнением Фредгольма второго рода. Здесь

λ - числовой параметр, f (t), K(t, s) - заданные функции переменных t, s [a,b], x(t) - искомая функция. Рассмотрим два случая:

a) f C[a,b] и K C[a,b]×[a,b]; b) f L2 [a,b] и K L2 [a,b]×[a,b] .

Метод последовательных приближений решения уравнений (1) состоит в построении последовательности

 

xn+1 (t) = λb K (t, s)xn (s)ds + f (t), n = 0,1, 2,...,

 

(2)

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

начиная с

произвольно выбранной функции x0 C [a,b] в

 

случае а)

и

x0 L2 [a,b]

в случае b). Если выполнено

условие сжатия

 

q <1, где

q

вычисляется в случаях а) и b) по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

b b

1

 

 

 

a) q =

 

λ

 

maxatb

 

K (t, s)

 

ds,

b)

 

λ

 

(∫∫

 

K (t, s)

 

2 dtds)

 

,

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

a a

 

 

 

то xn (t) сходятся к решению x* (t)

в соответствующих пространствах. При

этом справедлива оценка погрешности

 

 

 

 

 

 

 

 

 

 

 

xn

x*

 

qn

 

 

 

x1 x0

 

 

 

.

 

 

 

 

 

 

 

 

 

 

1 q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнением Вольтерры называется частный случай уравнения Фредгольма вида

x(t) = λt

K (t, s)x(s)ds + f (t) .

(3)

a

 

 

Последовательные приближения, построенные для уравнения (3), в отличие от уравнения Фредгольма, сходятся при любом значении λ к единственному решению уравнения Вольтерры.

Пример 1. Решить интегральное уравнение

x(t) =

1

1 tsx(s)ds +

5

t .

2

6

 

0

 

18

Решение. Проверим, что выполнено условие разрешимости этого уравнения в пространстве С[0,1] . Действительно,

 

 

 

q =

1

max 1 tsds =

1

 

1

=

1

<1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 0t1 0

 

2

 

2

 

4

 

Положим

x0

(t) 0 . Тогда x1

(t) =

5

t, x2 (t) =

5

t(1 +

 

1

) .

 

По индукции докажем,

 

6

6

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что справедливо представление

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn (t) =

 

5

t(1 +

1

+... +

 

 

 

1

 

) .

(4)

 

 

 

6

6

6(n1)

 

 

 

 

 

 

 

 

 

 

 

 

 

Для n = 2 эта формула верна. Предположим, что для всех номеров до n -го порядка она справедлива и проверим ее для n +1. Имеем

 

 

xn+1 (t) =

1

t

1

s

5

s(1 +

 

1

+... +

 

1

 

)ds +

 

5

t =

 

 

 

 

 

2

0

6

6

6(n1)

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

1

 

 

 

 

 

1

 

 

 

1

1

2

 

 

5

 

 

5

 

 

1

 

 

 

 

1

 

 

t(1

+

 

+... +

 

 

 

 

 

)

 

s

 

 

ds +

 

t

=

 

 

t(1 +

 

 

+... +

 

 

).

6

6

6

(n1)

 

2

 

 

6

6

6

 

6

n

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, представление (4) справедливо. Поэтому

 

5

1

 

5

 

1

 

 

x* (t) = lim xn (t) =

t

=

t

 

= t .

6

k

6

1 1

n→∞

k =0

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

Проверка показывает, что решение найдено верно. Выпишем оценку погрешности n -го приближения:

xn x*

 

 

 

 

1

 

 

 

 

5

t

 

=

10

 

.

 

 

 

 

 

 

 

C[0,1]

 

 

 

 

 

 

 

n

 

 

4

 

(1

4)

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

1

 

 

 

6

C[0,1]

 

9 4

 

 

Пример 2. Решить уравнение

x(t) =1 +31 ts2 x(s)ds .

(5)

0

Решение. Условие разрешимости в C[0,1] не выполняется, так как

19

q = 3max t1

s2 ds = 3

1

=1.

3

0t1 0

 

 

Но условие разрешимости в L2 [0,1] выполнено, поскольку

1

1

1

 

 

 

1

 

q = 3(∫∫t 2 s4 dtds)

 

= 3(

1

 

1

)

 

<1,

2

2

3

5

0

0

 

 

 

 

 

 

 

что гарантирует существование

единственного решения x* L2 [0,1]

уравнения (5). Найдем его методом последовательных приближений.

Положим

x0 (t) =1,

тогда

x1 (t) =1 +t ,

x2 (t) =1 +t(1 +

3

) .

С

помощью

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

математической

 

 

индукции

 

нетрудно

 

 

 

убедиться,

что

n1

 

3

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn (t) =1 +t(

)k

=1 +t

 

 

=1 + 4t . Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

k =0

4

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

3

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x* (t) =1 +t

(

)k

=1 +t

 

=1 + 4t .

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

k =0

4

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

Оценим погрешность n -го приближения. Получим

 

 

 

 

 

 

 

 

 

 

 

 

xn x*

 

(

3 )n

15

t

 

= (

3)

n1

 

1 .

 

 

 

 

 

 

L

2

 

 

 

 

 

 

 

 

L

15

15

3

2[ 0 ,1]

 

 

5

 

 

 

15 3

 

 

 

 

 

 

 

 

 

2[ 0 ,1]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 3. Решить уравнение Вольтерры

t

x(t) = t (t s)x(s)ds .

0

Решение. Это уравнение имеет единственное решение x* (t) . Вычислим последовательные приближения. Положим x0 (t) = 0 , тогда x1 (t) = t ,

t

 

 

 

 

 

t

2

 

t

3

 

 

 

1

 

 

1

 

t

3

.

x2 (t) = t (t s)sds = t t

 

+

 

= t

t 3 (

 

) = t

 

 

 

 

 

 

3

3!

0

 

 

 

 

2

3

 

2

 

 

 

 

Покажем, используя метод индукции, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

(t) = t

t 3

+

t 5

... + (1)n1

 

t 2n1

 

 

.

 

(6)

 

 

 

(2n 1)!

 

 

3!

5!

 

 

 

 

 

 

 

 

 

 

 

 

 

20

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