Материал: 1692

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

у = (x) с прямой у = x. Взяв в качестве начальной произвольную точку x0 [a, b], строим ломаную линию. Абсциссы вершин этой ломаной представляют собой последовательные приближения корня . Из рисунков видно, что если '(x)<0 на отрезке [a, b], то последовательные приближения xn = (xn-1), колеблются около корня , если же производная '(x)>0, то последовательные приближения сходятся к кор-

ню монотонно.

 

a)

б)

Рис. 10

При использовании метода простых итераций основным моментом является выбор функции у = (x), эквивалентной исходной.

На рис. 11 рассмотрен пример, когда условие окончания ите-

рационного процесса y xn 1 выполняется на первом шаге

итерационного процесса, т.е. x1 x0 , из этого следует, что х0 является приближенным значением искомого корня. Однако из рис. 11 видно, что это неверно, т.к. решением задачи является .

Для метода итераций следует подбирать функцию (x) так, чтобы | '(x)| δ <1, в противном случае процесс итерации расходящийся. При этом следует помнить, что скорость сходимости последовательности {xn} к корню тем выше, чем меньше число δ.

Ключевой момент в применении метода простой итерации – эквивалентное преобразование уравнения F(x)=0 к виду (2). Конечно,

31

Рис. 11

такое преобразование имеет смысл только тогда, когда оказывается

выполненным условие

 

при 0 q 1. Если первая (обычно

(х) q

самая простая и напрашивающаяся) попытка представления уравнения в требуемом виде оказалась неудачной, отчаиваться не следует. В ряде случаев можно использовать специальные приемы. Рассмотрим некоторые из них [2,5].

Способ 1. Если f (x) содержит в себе выражение некоторой об-

ратимой на [с; d] функции

y (x), причем такой, что

 

q 1

на

(x)

[c; d], то следует попытаться заменить уравнение на равносильное с использованием обратной для функции : x (y). Этот способ основан на известном соотношении между производными взаимообратных функций (y) 1/ (y) и следствии изнего:

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если

 

(x)

q 1,

то

(y)

 

 

 

 

 

 

 

 

1.

 

 

 

 

 

 

 

 

 

 

 

(x)

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Привести уравнение x3 3x 1 0 к виду,

пригодному для

решения методом простой итерации на интервале [0,8; 2].

Прибавим к правой и левой частям х

и получим:

x x3 3x 1 x.

Проверим условие сходимости:

 

 

 

 

 

 

 

 

 

 

 

 

(x) x

3

 

 

 

 

 

 

 

2

 

2;

 

3x 1 x (x) 3x

 

 

 

 

 

 

х [0,8;2],

 

 

 

 

 

 

(2) 10 1 при

 

 

 

 

 

условие сходимости не выполняется.

Другой вариант уравнения: x x3 1. Проверим условие сходи- 3

мости:

32

(x) 1 x3 1 (x) x2 ;

33

(2) 4 1 при х [0,8; 2],

условие сходимости не выполняется.

Так как ни одно из приведенных нами уравнений не удовлетворяет условию сходимости, то применим описанный способ:

x3 3x 1 x 33x 1;

 

 

 

 

 

 

1

 

 

 

(x)

3

3x 1

 

 

 

 

 

 

 

 

 

 

 

;

 

(x)

 

 

 

 

 

3

(3x 1)2

 

 

 

 

 

 

 

 

(0,8) 0,799 1при х [0,8; 2]

условие сходимости выполняется. Таким образом, для уточнения нужного нам корня методом простой итерации можно использовать уравнение x 33x 1.

Способ 2. В случае, когда способ 1 применить трудно или он не даст нужного результата, можно использовать следующий прием.

Пусть дано уравнение с единственным корнем в [a; b]. Предположим, что на отрезке [с; d] производная f функции F непрерывна, не равна константе и принимает значения одного и того же знака.

Будем считать,

 

 

т.к. в противном случае можно

что f (x) 0,

рассматривать равносильное уравнение: f (x) 0.

 

Введем обозначения:

 

 

1

 

m

 

 

 

 

 

 

 

 

m min f

M max f

k

и q 1-

.

(x),

(x),

[c;d ]

 

[c;d]

 

 

M

 

M

 

 

 

 

 

 

 

 

 

Ясно,

что

0 q 1.

Заменим

равносильное уравнение

уравнением эквивалентным ему

 

 

 

 

xx k f (x)

ипокажем, что для функции g(x) x k f (x) на [c; d] имеет место ус-

ловиесходимости.

Для x [c, d] справедливы неравенства: 0 m f (x) M . Разделим их почленно на число М и для разностей между единицей и полученными дробями получимнеравенство:

0 1 f (x) 1 m q,

M M

откуда ивытекает, что

0 g (x) 1 k f (x) q

при всех x [c, d].

33

Пример № 1. Привести уравнение lnx

1

0 к виду, пригодному

x2

 

 

для решения методом простой итерации на интервале [1,4; 1,7].

Так как условие сходимости не выполняется, то применим второй способ приведения уравнения:

 

 

 

 

 

 

 

 

 

f

 

 

1

 

2

;

 

 

 

 

 

 

 

 

 

 

 

x3

 

 

 

 

 

 

 

 

 

 

(x) x

f

 

 

1

 

 

 

2

 

0,71 0,73 1,44;

1,4

1,43

 

(1,4)

 

f

 

 

1

 

 

 

2

 

0,59 0,40 0,99;

1,7

 

1,73

 

(1,7)

 

 

 

m min f

 

 

 

 

 

 

 

 

 

m 0,99;

 

 

(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M 1,44;

 

M max f

(x)

 

 

k

1

 

 

 

 

1

 

 

 

 

0,69;

 

 

M

 

 

 

 

 

 

 

 

 

 

1,44

 

 

 

 

 

 

 

q 1 m 1 0,99 0,31 1. M 1,44

Таким образом, для уточнения нужного нам корня методом простой итерации можно использовать уравнение

x x 0,31(ln x

1

) .

 

x2

 

 

 

 

 

 

Пример № 2. Привести уравнение

lnx

 

1

0 к виду, пригодному

 

 

x2

 

ln10

 

 

для решения методом простой итерации на интервале [1,7; 2,1].

Так как условие сходимости не выполняется, то применим второй способ приведения уравнения:

 

 

 

 

1

 

 

2

 

 

 

 

x ln10 x3 ;

 

 

f (x)

f

 

1

 

 

2

0,26 0,40 0,66;

1,7 ln10

1,73

(1,7)

f

 

1

 

 

2

0,20 0,22 0,42;

 

3

(2,1)

2,1 ln10

 

 

 

2,1

 

 

 

 

34

 

 

 

 

 

 

m 0,42;

m min f (x)

 

 

 

 

 

 

M 0,66;

M max f (x)

k

1

 

1

1,52;

 

 

0,66

 

 

M

 

 

 

q 1 m 1 0,42 0,36 1. M 0,66

Таким образом, для уточнения нужного нам корня методом простой итерации можно использовать уравнение

x x 1,52(

lnx

 

1

).

ln10

x2

 

 

 

1

 

Пример № 3. Привести уравнение

1

x

e

0 к виду, пригодно-

x

 

 

1

x

 

 

му для решения методом простой итерации на интервале [0,3;0,7]. Так как условие сходимости не выполняется, то применим вто-

рой способ приведения уравнения:

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x) (1 x)2 x2 e ;

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

f

 

 

 

 

 

 

 

 

e

0,3

 

 

 

(0,3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

315,52;

(1 0,3)2

0,32

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

f

 

 

 

 

 

 

 

 

 

e

0,7

 

 

(0,7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30,73;

 

(1 0,7)2

 

 

0,72

 

 

 

 

 

 

m min f (x)

 

 

 

m 30,73;

 

 

 

 

 

 

 

 

 

 

M 315,52;

M max f (x)

 

 

k

1

 

 

 

1

0,003;

 

 

 

 

 

 

 

 

 

 

 

315,52

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q 1 m 1 30,73 0,9 1. M 315,52

Таким образом, для уточнения нужного нам корня методом простой итерации можно использовать уравнение

 

1 x

 

1

 

 

 

 

x

 

x x 0,003

 

e

 

 

.

1 x

 

 

 

 

 

 

 

 

 

 

 

35