Материал: 5544

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

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

x 1.

 

 

 

 

 

 

 

 

f 0 02

 

 

 

 

 

 

Попробуем

взять

 

x

0

 

и

 

При

этом

 

0 3 3

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f 1 12

 

3 3. Выбор неудачен, поскольку при поиске x2

произойдёт де-

1

ление на 0. Секущая горизонтальна и не пересекает ось OX.

 

 

 

 

 

 

 

 

 

2 , тогда

 

 

f 2 22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возьмём x

 

 

 

 

 

2 3 0,414 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак,

 

x0 0, f 0 3, x1

 

2, f 2 0,414 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2-й шаг. Находим x2

 

2 0,414 2

 

 

2 0

 

 

 

2 0,414 2

 

2

 

 

 

2,320 4 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,414 2 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,585 8

 

 

при этом f 2,320 4 2,320 42

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,320 4 3 0,8611. Далее,

 

 

 

 

 

 

 

x3

2,320 4 0,8611

 

 

 

2,320 4 2

 

 

2,104 0 , f 2,104 0 0,023 7 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,8611 0,414 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

2,104 0 0,023 7

 

2,104 0 2,320 4

2,109 8 ,

f 2,109 8 0,0013;

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,023 7

0,8611

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x5

2,109 8 0,001 3

 

 

 

2,109 8 2,104 0

 

 

2,1011.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,001 3 0,023 7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3-й шаг. Поскольку

 

x5 x4

 

 

 

 

 

2,1101 2,109 8

 

0,000 3 0,001, точность

до-

 

 

 

 

стигнута. Более того,

f 2,1101 0,000 096 0,0001 . Дальнейшие действия с 4

знаками не изменят результат:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x6 2,1101 0,000 096

 

 

 

2,1101 2,109 8

 

 

2,1101.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,000 096 0,0013

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ответ: решение уравнения с точностью 0,001 – число x 2,110 .

 

 

 

 

 

При программировании метода надо учесть,

что знаменатель f xn f xn 1

может быть близок к 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисления можно оформить в виде таблицы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

xn 1

xn

 

 

 

 

f xn 1

 

 

 

 

 

 

 

 

f xn

 

 

xn xn 1

 

 

f xn f xn 1

 

 

 

xn 1

 

 

1

 

0

 

 

2

 

 

 

 

 

 

 

–3

 

 

 

 

 

–0,414 2

 

 

2

 

 

 

 

2,585 8

 

 

2,320 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

2,320 4

 

–0,414 2

 

 

 

 

0,861 1

 

 

0,320 4

 

 

1,275 3

 

 

2,104 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

2,320 4

2,104 0

 

 

 

 

0,861 1

 

 

 

 

–0,023 7

 

 

–0,216 4

 

 

–0,884 8

 

 

2,109 8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

2,104 0

2,109 8

 

–0,023 7

 

 

 

–0,001 3

 

 

0,005 8

 

 

0,022 4

 

 

2,110 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

2,109 8

2,110 1

 

–0,001 3

 

 

 

–0,000 096

 

 

0,000 3

 

 

–0,001 204

 

2,110 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В каждой строке, начиная со 2-й, три значения берутся из предыдущей. Поясним сущность описанного метода хорд, т.е. применяемую формулу (2.4).

Предположим, что функция f(x) уравнения (2.1) на исходном отрезке [a,b] имеет непрерывные производные и постоянного знака ( <0 или

>0, f (x) <0 или f (x) >0), причём её значения f(a) и f(b) на концах этого отрезка

– противоположных знаков.

Тогда на отрезке [a,b] будет строгая монотонность, сохраняться выпуклость (вогнутость), а кривая пересечёт ось ОХ только один раз, т.е. x c будет единственным корнем уравнения (2.1) на [a,b]. Рекомендуем читателю сделать четыре рисунка возможных ситуаций.

Дугу кривой y f x на отрезке [a,b] заменим хордой, соединяющей точки М1(a, f(a)) и М2(b, f(b)). За первое приближение х1 к корню с примем абсциссу точки пересечения этой хорды с осью ОХ. Очевидно, что при этом будет выполняться неравенство c x1 b a . Теперь укажем формулу для х1. Для этого со-

ставим уравнение прямой, проходящей через две точки М1

и М2:

 

 

 

 

y f (a)

 

 

x a

.

(2.5)

 

 

 

 

f (b) f (a)

 

 

 

 

 

 

 

b a

 

Положив в этом уравнении

 

y 0 , найдём точку пересечения прямой (следо-

вательно, и хорды) с осью ОХ:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (a)

 

 

x a

.

 

 

 

f (b) f (a)

 

 

 

 

 

 

b a

 

Отсюда очевидным образом имеем, что

 

 

 

 

 

 

 

 

x1

a

b a

f (a).

 

 

 

 

 

 

 

f (b) f (a)

 

 

 

 

 

 

 

 

 

 

Последнее равенство и приводит к расчётной формуле (2.4), применяемой на отрезке [xn 1, xn ].

Уравнение (2.5) приводится к виду

y f (a)

f (b) f (a)

(x a)

b a

 

 

и, следовательно, представляет собой линейную функцию y kx d . Это означает, что на отрезке [a,b] этой функцией заменена функция f(x) рассматриваемого уравнения (2.1). Поэтому метод хорд (метод секущих) называют ещё методом

линейной интерполяции.

2.3. Метод касательных

Метод касательных применяют для функций, дифференцируемых на интервале около предполагаемого решения. Метод особенно удобен, если функция быстро меняется вблизи корня. Распространённое название «метод Ньютона»

32

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

Схема метода касательных

Надо найти корень с уравнения (2.1) с точностью .

1)Находим общую формулу для вычисления производной f x ;

2)выбираем x0 – начальное приближение;

3)каждое следующее приближение находим по формуле

x

x

 

 

f xn

 

.

(2.6)

 

f xn

 

n 1

 

n

 

 

 

 

 

 

 

 

 

Формула (2.6) является расчётной формулой метода касательных;

4)вычисления прекращаем, когда xn 1 xn ;

5)в ответе указываем xn 1 .

Точку x0 следует выбирать так, чтобы в ней значение функции и 2-й произ-

водной были одного знака: f x0 f x0 0 . В этом случае касательная пересекает ось абсцисс (ось OX) около корня функции, т.е. примерно там же, где ось пересекается с графиком функции. Таким образом, метод заключается в замене на каждом шаге функции f(x) на касательную и в поиске корня этой касательной.

Действительно, уравнение касательной к кривой y f x в точке xn ; f xn

имеет вид

y f (xn ) f (xn )(x xn ).

Положив в этом равенстве y 0 , найдём координату х пересечения этой кривой с ОХ. Это х и берём за xn 1 , т.е. за следующее приближение к корню уравне-

ния (2.1). Это и приводит к формуле (2.6).

Если же условие f x0 f x0 0 нарушено, т.е. знаки f x0 и f x0 различны, то касательная и график пересекают ось ОХ в точках, расположенных далеко одна от другой, и метод расходится.

Метод касательных устойчив к ошибкам: поиск корня можно начинать с любой точки, а случайная (не систематическая) ошибка лишь приводит к дополнительному шагу приближения и не опасна.

Пример. Решим уравнение x3 6 x с точностью 0,001.

Решение. Здесь f x x3 x 6 , поскольку сначала надо свести уравнение к виду f x 0 , в нашем случае x3 x 6 0 .

33

1-й шаг. Дифференцируем f x , получаем

 

f x 3x2 1.

И функция, и производная определены при всех x.

2-й шаг. Выбираем любое

x

0

.

Пусть

x

0

0 ,

тогда f 0 03 0 6 6 и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f 0 3 02 1 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3-й шаг. Находим x

0 6

 

6 ,

f 6 63

6 6 216 и f 6 3 62 1 109 .

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На следующем шаге x

 

6

216

4 , вычисления удобно проводить в таблице:

2

109

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

xn

 

 

f xn

 

f xn

 

xn 1

 

 

0

 

 

0

 

 

 

–6

 

 

1

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

6

 

 

 

216

 

 

109

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

4

 

 

 

62

 

 

49

 

2,734 7

 

 

 

 

 

 

 

 

 

 

 

3

2,734 7

 

17,186 4

23,435 8

 

2,001 4

 

 

 

 

 

 

 

 

 

 

 

 

4

2,001 4

 

4,018

2

13,016 8

 

1,692 7

 

 

 

 

 

 

 

 

 

 

 

 

5

1,692 7

 

0,542

7

9,595 7

 

1,636 2

 

 

 

 

 

 

 

 

 

 

 

 

6

1,636 2

 

0,016

6

9,031 5

 

1,634 4

 

 

 

 

 

 

 

 

 

 

 

 

7

1,634 4

 

0,000

3

9,013 8

 

1,634 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Последний столбец необязателен, новое приближение можно считать под предыдущим. На седьмом шаге x7 x8 1,634 4 . Очевидно, это значение полу-

чим и на всех последующих шагах. Вычисления закончены. Для уточнения корня необходимо вести вычисления с большим числом цифр после запятой.

Правильность решения проверим подстановкой в исходное уравнение:

1,634 43 4,365 9 , 6 1,634 4 4,365 6 ,

разность 4,365 9 4,365 6 0,000 3 0,001 и не превосходит точности вычислений. Ответ: x 1,634 – решение уравнения x3 6 x с точностью 0,001.

2.4. Комбинированный метод

Комбинированный метод получается при объединении метода секущих и метода касательных. Пусть дано уравнение (2.1). Возьмём две точки, обозначим их x0 и t0 . Пусть для определённости x0 t0 и функция возрастает выпуклостью вниз. Проведём касательную к графику f x в точке с абсциссой t0 . Получим точку t1 – новое приближение, найденное методом касательных.

34

Затем соединим точки графика с абсциссами x0 и t1 . Получим x1 – новое приближение, найденное методом секущих, причём не на отрезке x0 ; t0 , а на от-

резке x0 ; t1 . Этот отрезок короче, поэтому x1 быстрее перемещается к решению уравнения, чем при обычном методе хорд. В этом и есть смысл и цель метода. Соответствующие формулы пересчёта:

 

tn 1 tn

 

 

f tn

 

,

 

 

 

 

 

 

xn 1 xn f

xn

 

x

n

t

n 1

 

,

 

 

f

tn

 

 

 

 

 

 

f

xn

f tn 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где x0 , t0

– произвольное начальное приближение.

 

 

 

 

 

 

 

 

 

 

Метод заведомо сходится, если f x0 f t0 0

и f t0 f t0 0 . Формулы пере-

счёта не меняются, если x0 t0

и (или) функция

 

f x ведёт себя иначе, чем ука-

зано в предположении.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 1. Решим уравнение x3 x 5 0 с точностью 0,001.

 

 

 

 

 

 

Решение.

Здесь

f x x3

x 5 .

 

 

Заметив,

 

что

 

f 1 13

1 5 3 0 и

f 2 23

2 5 5 0 ,

видим,

что в качестве начального приближения можно

взять x0

1 и t0 2 , а корень находится на отрезке 1;2 . Также нам понадобятся

f x 3x2 1 и f x 6x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

f 2 6 2 12 0 , то

 

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

сходимости метода

касательных (функция и 2-я производная одного знака).

 

 

 

 

 

 

 

Сначала уточняем правый конец отрезка.

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим

f 2 3 22 1 13 , тогда

t

1

2 5 /13 1,615 4 ,

где 5 f 2 . По-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

скольку

f 1,615 4 1,615 43 1,615 4 5 0,830 8 , получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 1 3

 

 

1 1,615 4

1,481 9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3 0,830 8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и новый отрезок 1,4819;1,615 4 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Повторяя цикл вычислений для нового отрезка, находим

 

 

 

 

 

 

 

f 1,615 4 3 1,615 42 1 8,828 6 и t

 

 

1,615 4

0,830 8

1,5213 .

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8,828 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим значения функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f 1,481 9 1,481 93 1,481 9 5 0,263 8 ,

f 1,521 3 1,521 33 1,521 3 5 0,042 1 ,

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

1,4819 0,2638

 

1,4819 1,5213

 

1,515 9 ;

 

 

 

 

 

 

 

2

0,2638 0,0421

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

новый отрезок 1,515 9;1,5213 .

35

Источник: https://studfile.net/preview/16711204/