Файл: Итенберг С.И. Дифференциальное и интегральное исчисление функций одной переменной учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 14.07.2024

Просмотров: 136

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

равенства (3.3). Применив формулу (3.4) еще раз и считая, что теперь известно приближенное значение х-,, найдем новое прибли­ жение к корню

Л«Э

Л 1

 

 

 

Поступая аналогично еще раз, найдем

 

 

Л д

Л о

 

*

 

 

и вообще

 

 

 

 

 

- т £ т < »

= 0,

1, 2,

. . . ) •

(3.5)

f

(*«)

 

 

 

хп+\

Теперь естественно спросить,

при

каких

условиях точка

(« = 0 , 1 , 2 , . . .) действительно будет ближе к точке х, чем точка

хп.

Ответ на этот вопрос дает следующая теорема, которую мы приве­

дем

без доказательства.

 

 

 

 

 

 

 

 

 

 

 

 

Теорема. Если

функция

f (х)

удовлетворяет

на

интервале,

[а,

Ъ] перечисленным в

начале

настоящего

параграфа

условиям

и

величины f (хп)

и f"

П)

(п =

0,

1,

2, . . .)

одного знака,

то

точка

хп+\

 

расположена ближе

к точке х,

чем точка

хп.

 

 

 

 

 

Обычно вычисления следует вести до тех пор, пока в пределах

требуемой точности

не получится: хп+\ ^

хп.

 

 

 

 

 

 

 

Метод Ньютона имеет простое геометрическое истолкование.

Построим на промежутке [а,

Ь] график функции f (х). По условию

/'

(х) и /" (х)

непрерывны,

знакопостоянны' и

отличны от

нуля.

Геометрически

это

означает,

что в

любой

точке интервала

[а,

Ы

к

кривой у — f (х)

можно провести касательную,

что

кривая

эта

на интервале

[а, Ъ ] не имеет экстремумов

и либо только выпукла,

либо только вогнута. Так как на интервале [с,

Ъ ] имеется единст­

венный корень уравнения, то кривая в единственной точке х пе­

ресекает ось абсцисс. Выберем на кривой произвольную

точку

MQ

с

абсциссой

х0, но такую,

что

/ 0) f"

( х 0 ) > 0

(это значит,

что,

/

0) и /" 0)

одного знака),

и

проведем

к этой

кривой

касатель­

ную:

У f (*о) — F (*i>) (* — Хо).

Полагая в последнем уравнении у = 0, найдем уравнение для

определения точки хх пересечения

касательной

с осью Ох:

— f(x0):=f'(Xo)(x1

— x0),

 

откуда получим формулу (3.4).

Нам известно, что точка хг будет лежать ближе к х, чем точка х0.

Теперь проведем касательную к кривой у = / (х) в точке Мх

с абс­

циссой хх и найдем

 

 

Луо -У Л 1

/ f a )

 

— точку пересечения новой касательной с осью Ох. Точка хг

будет

лежать ближе к х, чем точка хг, и т. д.

 

67


Пользуясь геометрической терминологией, метод Ньютона на­ зывают м е т о д о м к а с а т е л ь н ы х . На рис. 43 наглядно показано построение касательных и вычисление приближенных значений искомых корней уравнения. На этом рисунке представ­ лен случай, когда на интервале [a, b ] f" (х) > 0, а так как / О, то за точку х0 мы взяли точку Ь.

Пример. С помощью метода Ньютона найти с точностью до 0 , 0 0 0 1

корень

уравнения

f (х) = х3 — 5 =

0.

 

 

 

В § 3

. 1 мы нашли интервал

[ 2 ; 2 , 5 ] , где изолирован

корень

заданного

уравнения

(рис. 4 2 ) . На этом интервале /' (х) = ЗА:3 2

и /" (х)

=

6л: не

Рис. 43

изменяют знак и не обращаются в нуль, причем на этом интервале /" (л-) > 0 . За нулевое приближение к искомому корню примем число х0 = 2 , 1 , что можно сделать, так как f, ( 2 , 1 ) > 0 . Вычисляем:

-

/ Ы

- 2 , 1 - Д . " - 2 - 2 , 1 - 5

 

ГЫ

3 - 2 . 1 3 2

= 2 | Q 9 4 6 >

, 2 = X l _ Ж .

= 2 , 0 9 4 6 - * 1 « » - 4 , 1 8 Я 2 - 5 = 2,0945 .

Г(хх)

1 3 , 1 6 1 9 — 2

Очевидно, что число х = 2 , 0 9 4 5 можно принять за приближенное зна­

чение корня.

3.3. М Е Т О Д Х О Р Д

Предположим, что левая часть уравнения f (х) = 0 на интер­ вале [a, b ] удовлетворяет условиям, указанным в начале предыду­ щего параграфа. Тогда численное значение действительного корня заданного уравнения можно найти с помощью так называемого ме­ тода хорд.

Проведем через две точки (a, f (а)) и (b, f (b)) прямую

y~f ^

= х ~ а и

найдем точку х±

f(b)—f(a)

Ь — а

 

 

 

положив,

у = 0 . Получим

 

 

 

0 — f (а)

Хл—а

' v

~ а

 

f(b)-f{a)

Ь-а

 

 

 

пересечения

ее с осью

Ь — а

, , .

— Т77Т 7 7 - Г - / ( А ) f(b)-f(a)

Ох,

с\

( 3 - 6

)

68


или

хх = b-

f{b)-f{a)

•f(b).

(3.7)

 

 

 

Число xx принимаем за приближенное значение корня уравне­ ния. Далее следует вычислить значение / х) и определить, в ка­ ком из .интервалов [а, хх], [хх, Ь] находится корень уравнения. Если f (а) и / х) разных знаков, то для дальнейших вычислений следует применять формулу типа (3.6), если же / х) и / (Ь) разных знаков, то надо воспользоваться формулой типа (3.7). На рис. 44 показан геометрический смысл метода хорд.

Пример. Решить

уравнение

х3 — 5 =

0 с помощью метода хорд.

Корень данного уравнения отделен в интервале

[2; 2,5] (рис. 42). Вычислим

первое приближенное

значение

корня

хх

 

h

п

 

2 5 2

 

 

- / (а) = 2

^

— • ( — ! ) = 2,0755.

 

 

 

5,625

+ 1

а)

Рис. 44

Вычисления показывают / (2,0755) <0, а так как / (2,5) >0, то корень лежит внутри интервала [2,0755; 2,5]. Находим второе приближение

 

x„ = 2,0755

^

A

U

/

0 D

(-0,2104) = 2,0908.

 

 

5,625 +

0,2104

 

Так

как f, (2,0908) < 0, то

 

 

 

 

 

 

 

х3 = 2,0908 -

9 с .

 

.9

попа

 

 

А

и

а и

в

(— 0,0418) = 2,0938;

 

 

5,625 +

0,0418

 

четвертое

приближение

/(2,0938)<0;

 

 

 

 

 

 

 

 

д-4 = 2,0938 — — — 2

' 0

9

3 8

 

(— 0,0084) = 2,0944;

 

 

5,625 +

0,0084

 

пятое приближение

/ (2,0944) <

0;

9 Я

 

9 (1Q44

 

 

хъ = 2,0944

 

(0,0017) = 2,0945.

 

 

 

0,0017

 

 

5,625 +

 

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

В предыдущих вычислениях мы считали, что корень уравнения найден с нужной степенью точности, если в пределах допустимой

69



погрешности хп ^ xn.v\. Однако, вообще говоря, это условие еще не гарантирует, что корень найден с достаточной точностью. Мо­ жет случиться, что два значения корня близки друг другу, но да­ леки от истинного значения корня уравнения.

3.4. МЕТОД ИТЕРАЦИЙ

 

 

Рассмотрим уравнение

 

 

 

.

/(х) = 0.

(3.8)

Если записать его в форме

 

 

 

х = у{х)

 

(3.9)

(что всегда можно "сделать, прибавив, например, к обеим частям исходного уравнения по х) и если взято любое число х0 из интер­ вала, в котором корень изолирован, то численное значение дейст­ вительного корня может быть найдено с помощью следующего приема, носящего название способа итерации.

 

Подставив в правую

часть

уравнения (3.9)

число

х0,

 

найдем

Х\

= Ф 0);

далее х2

=

ср г),

ха

= ср (х2 ), . . . , хп

=

ср (*„_,)

до тех пор, пока в пределах принятой точности не будет хп

^

xn+i.

 

Выясним

условия,

при которых

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

п]

= хх,

хг,

. . . , хп,

. . . стремится к корню х уравнения

(3.9), или при ко­

торых, как говорят, итерационный процесс сходится. Предполо­

жим, что функция

ф (х) дифференцируема в интервале изоляции

корня. Если х0 и хг

— соответственно нулевое и первое приближе­

ния корня х, то х =

ф (х), х± ф 0). Вычитая из первого равен-

. ства второе и применив к правой части полученного равенства фор­

мулу (2.5) конечных приращений Лагранжа, найдем

 

 

х—х1 = ф {х) — ц> (х„) =

(х—х0)

ср' (с,,),

(3.10)

где.с0 — число, лежащее

между х и

х0.

 

 

Аналогично,

х—х2=

(х—хх)

ф'

г),

 

 

 

 

 

 

 

 

х—х3

= (х—х2)

ф' (с2 ),

 

(3.11)

 

 

 

 

 

 

 

 

 

 

* - * n

= ( * - * ^ - i ) 4 > ' ( c » - i ) '

j

 

где xz, х3,

. . . , х

— приближенные значения корня х, а сг,

с2 , . . .

C_i числа, лежащие

соответственно

внутри

интервалов

[х, хг\,

[х, х 2 ] ,

. . . , [х,

хп_1].

Перемножая

почленно равенства (3.10)

и (3.11),

получим

 

 

 

 

 

 

 

 

х—хп

= (х-х0)-<р'

(с„) .ф'(С 1 ) • • - Ф ' ^ ) : ,

 

Если в рассматриваемом интервале наибольшее по модулю зна­ чение производной функции ф (х) не превышает единицы: max ф' |(х)| = т < 1 , то из предыдущего равенства вытекает оценка

| х—хп | < | х—х01т.".

70