Читайте также:
|
|
В этом и следующем разделе рассмотрим модификации метода Ньютона.Как видно из формулы (2.13), метод Ньютона требует для своей реализации вычисления производной, что ограничивает его применение. Метод секущих лишен этого недостатка. Если производную заменить ее приближением:
f '(xn)» ,
то вместо формулы (2.13) получим
xn +1 = xn –. . (2.20)
Это означает, что касательные заменены секущими. Метод секущих является двухшаговым методом, для вычисления приближения xn +1 необходимо вычислить два предыдущих приближения xn и xn – 1, и, в частности, на первой итерации надо знать два начальных значения x0 и x1.
Очередное приближение xn +1 получается как точка пересечения с осью OX секущей, соединяющей точки графика функции f(x) с координатами (xn -1, f(xn - 1)) и (xn, f(xn)).
Сходимость метода. Сходимость метода секущих устанавливает следующая теорема.
Теорема 2.4 Пусть x* – простой корень уравнения f(x) = 0, и в некоторой окрестности этого корня функция f дважды непрерывно дифференцируема, причем f"(x) № 0. Тогда найдется такая малая s-окрестность корня x*, что при произвольном выборе начальных приближений x0 и x1 из этой окрестности итерационная последовательность, определенная по формуле (2.20) сходится и справедлива оценка:
|xn + 1 – x*| Ј C |xn – x*| p, n і 0, p = » 1.618. (2.21)
Сравнение оценок (2.15) и (2.21) показывает, что p < 2, и метод секущих сходится медленнее, чем метод Ньютона. Но в методе Ньютона на каждой итерации надо вычислять и функцию, и производную, а в методе секущих – только функцию. Поэтому при одинаковом объеме вычислений в методе секущих можно сделать примерно вдвое больше итераций и получить более высокую точность.Так же, как и метод Ньютона, при неудачном выборе начальных приближений (вдали от корня) метод секущих может расходиться. Кроме того применение метода секущих осложняется из-за того, что в знаменатель расчетной формулы метода (2.20) входит разность значений функции. Вблизи корня эта разность мала, и метод теряет устойчивость.
Критерий окончания. Критерий окончания итераций метода секущих такой же, как и для метода Ньютона. При заданной точности e > 0 вычисления нужно вести до тех пор, пока не будет выполнено неравенство
|xn – xn – 1| < e. (2.22)
Дата добавления: 2015-05-05; просмотров: 13 | Поможем написать вашу работу | Нарушение авторских прав |