Метод секущих (метод хорд)

2.6 Метод секущих (метод хорд)

В этом и следующем разделе рассмотрим модификации метода Ньютона.

Как видно из формулы (2.13), метод Ньютона требует для своей реализации вычисления производной, что ограничивает его применение. Метод секущих лишен этого недостатка. Если производную заменить ее приближением:

f '(x_n) » ,

то вместо формулы (2.13) получим

x_{n +1} = x_n –. . (2.20)

Это означает, что касательные заменены секущими. Метод секущих является двухшаговым методом, для вычисления приближения x_{n +1}  необходимо вычислить два предыдущих приближения x_n и x_{n – 1} , и, в частности, на первой итерации надо знать два начальных значения x₀ и x₁.

Формула (2.20) является расчетной формулой метода секущих. На рис. 2.9 приведена геометрическая иллюстрация метода секущих.

Рис. 2.9

Очередное приближение x_{n +1} получается как точка пересечения с осью OX секущей, соединяющей точки графика функции f(x) с координатами (x_{n -1, f}(x_{n - 1})) и (x_n , f(x_n)).

Сходимость метода. Сходимость метода секущих устанавливает следующая теорема.

Теорема 2.4 Пусть x^*  – простой корень уравнения f(x) = 0, и в некоторой окрестности этого корня функция f дважды непрерывно дифференцируема, причем f"(x) ¹ 0. Тогда найдется такая малая s-окрестность корня x^*, что при произвольном выборе начальных приближений x₀ иx₁ из этой окрестности итерационная последовательность, определенная по формуле (2.20) сходится и справедлива оценка:

|x_{n + 1} – x^*| £ C |x_n – x^*| ^p, n ³ 0, p =  » 1.618. (2.21)

Сравнение оценок (2.15) и (2.21) показывает, что p < 2, и метод секущих сходится медленнее, чем метод Ньютона. Но в методе Ньютона на каждой итерации надо вычислять и функцию, и производную, а в методе секущих – только функцию. Поэтому при одинаковом объеме вычислений в методе секущих можно сделать примерно вдвое больше итераций и получить более высокую точность.

Так же, как и метод Ньютона, при неудачном выборе начальных приближений (вдали от корня) метод секущих может расходиться. Кроме того применение метода секущих осложняется из-за того, что в знаменатель расчетной формулы метода (2.20) входит разность значений функции. Вблизи корня эта разность мала, и метод теряет устойчивость.

Критерий окончания. Критерий окончания итераций метода секущих такой же, как и для метода Ньютона. При заданной точности e > 0 вычисления нужно вести до тех пор, пока не будет выполнено неравенство

|x_n – x_{n – 1}| < e. (2.22)

Пример 2.4.

Применим метод секущих для вычисления положительного корня уравнения 4(1 – x²) – e^x = 0 с точностью e = 10^-3.

Корень этого уравнения находится на отрезке [0, 1], так как f (0) = 3 > 0, а f (1) = –e < 0. Подсчитаем вторую производную функции: f "(x) = –8 – e^x. Условие f(x)f " (x) ³ 0 выполняется для точки b = 1. В качестве начального приближения возьмем x₀ = b = 1. В качестве второго начального значения возьмем x₁ = 0.5. Проведем вычисления по расчетной формуле (2.20). Результаты приведены в табл. 2.4.

Таблица 2.4

2.7 Метод ложного положения

Рассмотрим еще одну модификацию метода Ньютона.

Пусть известно, что простой корень x^*  уравнения f(x) = 0 находится на отрезке [a, b] и на одном из концов отрезка выполняется условие f(x)f"(x) ³ 0. Возьмем эту точку в качестве начального приближения. Пусть для определенности это будет b. Положим x₀ = a. Будем проводить из точки B = (b, f(b)) прямые через расположенные на графике функции точки B_n с координатами (x_n, f(x_n), n = 0, 1, … . Абсцисса точки пересечения такой прямой с осью OX есть очередное приближение x_n+1.

Геометрическая иллюстрация метода приведена на рис. 2.10.

Рис. 2.10

Прямые на этом рисунке заменяют касательные в методе Ньютона (рис. 2.8). Эта замена основана на приближенном равенстве

f '(x_n) » . (2.23)

Заменим в расчетной формуле Ньютона (2.13) производную f '(x_n) правой частью приближенного равенства (2.23). В результате получим расчетную формулу метода ложного положения:

x_{n +1} = x_n –.. (2.24)

Метод ложного положения обладает только линейной сходимостью. Сходимость тем выше, чем меньше отрезок [a, b].

Критерий окончания. Критерий окончания итераций метода ложного положения такой же, как и для метода Ньютона. При заданной точности e > 0 вычисления нужно вести до тех пор, пока не будет выполнено неравенство

|x_n – x_{n – 1}| < e. (2.25)

Пример 2.5.

Применим метод ложного положения для вычисления корня уравнения x³ + 2x – 11 = 0 с точностью e = 10^-3.

Корень этого уравнения находится на отрезке [1, 2], так как f (1) = –8 < 0, а f (2) = 1 > 0. Для ускорения сходимости возьмем более узкий отрезок [1.9, 2], поскольку f (1.9) < 0, а f (2) > 0. Вторая производная функции f (x) = x³ + 2x – 11 равна 6x. Условие f(x)f"(x) ³ 0 выполняется для точки b = 2. В качестве начального приближения возьмем x₀ = a = 1.9. По формуле (2.24) имеем

x₁ = x₀ –. = 1.9 +  » 1.9254.

Продолжая итерационный процесс, получим результаты, приведенные в табл. 2.5.

Таблица 2.5

Тема 3. Решение систем линейных алгебраических уравнений

Похожие работы

Лекции по вычислительной математике

Скачать

10356

0

0

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

Основные задачи вычислительной математики

Скачать

7721

1

1

... если  - предельная абсолютная погрешность приближённого числа , то  (1.2) отсюда следует, что  (1.3) Значение предельной абсолютной погрешности, обычно, подбирается интуитивно по смыслу задачи. Пример 2: Определить предельную абсолютную погрешность числа , заменяющего число , точное значение которого нам неизвестно. Так как мы знаем, что , ...

Выдающиеся личности в истории вычислительной техники. Августа Ада Лавлейс

Скачать

37333

0

0

... удивили меня…, хоть речь идёт обо мне самой. Они действительно написаны прекрасным стилем, который превосходит стиль самого очерка" /2/. 2.3. Рождение первенца и критическое перенапряжение Августа Ада Лавлейс работает с большим напряжением. В письмах к Бэббиджу она неоднократно жалуется на утомление, болезни, плохое самочувствие. Наконец, 6 августа Бэббидж отсылает Аде свои последние замечания ...

Роль женщин в развитии вычислительной техники

Скачать

54819

0

0

... в Украине, бывшем Советском Союзе и за рубежом научная школа теоретического программирования. В 2001-м году ее не стало... Но не только в научном плане велика роль женщин в развитии вычислительной техники. Со временем образуется огромное количество различных фирм по разработке и продаже программного и аппаратного обеспечения. Следовательно, разыгрываются человеческие трагедии капиталистического ...