ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.07.2025
Просмотров: 227
Скачиваний: 0
14
1.2.5 Метод Ньютона
Пусть на отрезке [a; b] отделен корень уравнения (1), функция f(x) дважды дифференцируема, а f ′(x) и f ′′(x) сохраняют постоянные знаки на указанном интервале.
В методе Ньютона точное значение ξ корня уравнения (1) заменяется приближенным значением x [a; b], где x – абсцисса пересечения касательной, проведенной к кривой y = f(x) в точке С [a; b]. Уравнение этой касательной
y – f(С) = f ′(С)(x – С). |
(17) |
|||||
Так как при x = x y = 0, то из (17) получаем |
||||||
f ( C ) |
||||||
= С – |
. |
(18) |
||||
x |
||||||
f ' ( C ) |
||||||
Остается решить вопрос о выборе точки С таким образом, чтобы x [a; b]. Выбор точки С определяется четырьмя случаями, изображенными на рис. 6.
y |
|||||||
f(x) |
|||||||
a |
C=b |
||||||
0 |
ξ |
x |
|||||
y |
|||||||
y |
2) f ′(x)<0, f ′′(x)<0, |
x |
|||||
C=b, f(C)<0 |
|||||||
f(x) |
|||||||
f(x) |
|||||||
a |
|||||||
0 |
a |
ξ |
0 |
b |
С=b |
x |
|
y |
x |
||||||
1) f ′(x)>0, f ′′(x)>0, C=b, f(C)>0 |
|||||||
3) f ′(x)>0, f ′′(x)<0, |
|||||||
x |
|||||||
C=a, f(C)<0 |
|||||||
ξ |
b |
x x |
|||||
0 |
a |
||||||
4) f |
Рис. 6 |
||||||
′(x)<0, f ′′(x)>0, |
|||||||
C=a, f(C)>0
x
15
Обычно принимают C = a или C = b, смотря по тому, в какой из этих точек знак функции совпадает со знаком второй производной, т.е. выбирают так, чтобы произведение sgn f(C) × sgn f ¢¢ (C) было положительно. Можно показать, что в этом случае a < x < b. Полученное значение x можно использовать для дальнейшего уточнения корня, беря интервал [a x ] (случаи 1 – 2) или интервал [ x ; b] (случаи 2 – 4).
Геометрический смысл метода Ньютона состоит в замене дуги y = f(x) касательной, проведенной к одной из крайних точек так, что имеет место итерационная формула
xk = xk–1 – |
f ( xk − 1 ) |
(k = 1, 2, …), |
f ' ( xk − 1 ) |
основанная на формуле (18).
Для оценки погрешности в методе Ньютона можно воспользоваться форму-
лой
|xk – z| £ |
M |
|xk – xk–1|2 , |
|
2m |
|||
min |
max |
||
где m = x [α ; |
β ] |f ¢(x)|, M = x [α ;β ] |f ¢¢(x)|. |
||
1.2.5.1Алгоритм метода Ньютона
1.Установить значения a, b, e – границы отрезка отделения корня и принятую точность приближения.
2. |
Определить точку (h, M = max | f '' ( x ) | ). |
|
x [ a ; b ] |
||
3. |
Вычислить m = |
min |
x [ a ; b ] |f '(x)|. |
||
4. Установить начальное приближение x0 Î [a; b]:
если sgn f ''(h) = sgn f(a), то |
x0 = a, |
иначе |
x0 = b. |
5.Начать цикл уточнения корня.
5.1.Вычислить очередное приближение x1 = x0 –
5.2.Вычислить d = (x1 – x0)2.
5.3.Принять x0 = x1.
6.Конец цикла, если d < 2Mm e .
7.Вывод результата x1.
8.Конец алгоритма.
f ( x0 ) . f ' ( x0 )
16
1.2.6 Метод хорд
Пусть на отрезке [a; b] отделен корень уравнения (1), функция f(x) дважды дифференцируема, а f ′(x) и f ′′(x) сохраняют постоянные знаки на указанном интервале.
Метод заключается в том, что на интервале [a;b] отделения корня ξ уравнения (1) дуга кривой y = f(x) заменяется стягивающей ее хордой и в качестве приближенного значения корня x принимается точка пересечения хорды с осью абсцисс (рис. 7).
Уравнение хорды определяем как уравнение прямой, проходящей через точки A(a, f(a)) и B(b, f(b)), имеющее вид
y
B(b, f(b)) f(b)
α |
|||
0 |
ξ |
β |
x |
f(a) |
A(a, f(a)) |
||
x − |
a |
= |
y − f ( a ) |
. |
b − |
a |
f ( b ) − f ( a ) |
||
(19) |
||||
Если x = x , то y = 0, тогда из (19) следует
Рис. 7. Геометрическая интерпретация метода хорд
x
x − a |
= − |
f ( a ) |
, |
b − a |
f ( b ) − f ( a ) |
откуда находим приближенное значение корня
x = a – |
f ( a ) |
(b – a). |
(20) |
f ( b ) − f ( a ) |
Для нахождения последующих приближений определяется отрезок, на концах которого функция имеет разные знаки. Если sgn f(a) = sgn f( x ), то полагается a = x , иначе b = x и повторяются вычисления по формуле (20).
Замечание. Вычисление приближенного значения x корня ξ уравнения (1)
выполняется с недостатком, т. е. x – ξ < 0, если на отрезке [a; b] имеют места неравенства:
f ′(x) > 0 и f ′′(x) > 0 – функция f(x) возрастающая, график f(x) вогнутый;
f ′(x) < 0 и f ′′(x) < 0 – функция f(x) убывающая, график f(x) выпуклый.
Вычисление приближенного значения x корня ξ уравнения (1) выполняется с избытком, т. е. x – ξ > 0, если на отрезке [a; b] имеют место неравенства: