Методы приближенного решения уравнений. Метод касательных

Классический метод Ньютона или касательных заключается в том, что если xn – некоторое приближение к корню x* уравнения , то следующее приближение определяется как координата пересечения касательной с осью ОХ, проведенной в точке xn.

Уравнение касательной к функции f(x) в точке xn имеет вид:

В уравнении касательной положим y=0 и x=xn+1.

Тогда алгоритм последовательных вычислений в методе Ньютона состоит в следующем:

(1)

Сходимость метода очень быстрая.

Здесь x0, x1, x2 … последовательные приближения к корню x*, полученные в результате применения метода Ньютона. Ясно, что сходимость последовательности {xk} к корню зависит от свойств функции f(x) и не всегда имеет место. Так, легко представить, что уже приближение x1 не попадает на исходный интервал и процесс итераций останавливается. Приведём полезную теорему, гарантирующую в некоторых случаях сходимость метода Ньютона.

Теорема

Если , причём отличны от нуля (и, следовательно, сохраняют определённые знаки при x ∈ [a, b]), то, исходя из начального приближения x0 ∈ [a, b], удовлетворяющего условию , можно вычислить методом Ньютона по формуле (1) единственный корень x* уравнения (1) с любой степенью точности.

Замечание. Практическим критерием окончания вычислений является выполнение условия |xn+1 − xn| < ε, где ε – требуемая точность вычисления корня.

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

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

 

Обоснование

Чтобы численно решить уравнение f(x)=0, его необходимо привести к эквивалентному уравнению: , где — сжимающее отображение.

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

В предположении, что точка приближения «достаточно близка» к корню, и что заданная функция непрерывна, окончательная формула такова:

С учётом этого функция определяется:

При некоторых условиях эта функция в окрестности корня осуществляет сжимающее отображение[1], и алгоритм нахождения численного решения уравнения f(x)=0 сводится к итерационной процедуре вычисления:

По теореме Банаха последовательность приближений стремится к корню уравнения f(x)=0

 

Метод секущих

Родственный метод секущих является «приближённым» методом Ньютона и позволяет не вычислять производную. Значение производной в итерационной формуле заменяется её оценкой по двум предыдущим точкам итераций:

Таким образом, основная формула имеет вид

Этот метод по сравнению с методом Ньютона имеет меньшую скорость сходимости.

 

Лабораторная работа №1.

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