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

Метод секущих — это итерационный численный метод приближённого нахождения корня нелинейного уравнения, относящийся к классу методов решения нелинейных уравнений. В отличие от метода Ньютона, метод секущих не требует вычисления производной функции, а аппроксимирует её с помощью секущей — прямой, проходящей через две последние точки на графике функции. Метод обладает более низкой скоростью сходимости по сравнению с методом Ньютона, но может быть предпочтительнее в случаях, когда аналитическое выражение производной сложно или невозможно получить.

История

Идея использования секущей для приближённого нахождения корня восходит к работам древнегреческих математиков, однако в современном виде метод был формализован в XVIII—XIX веках в рамках развития численных методов анализа. Значительный вклад в его обоснование внесли такие математики, как Исаак Ньютон, Джозеф Рафсон и Карл Фридрих Гаусс, хотя сам метод секущих часто рассматривается как разностный аналог метода Ньютона. В XX веке, с развитием вычислительной техники, метод получил широкое распространение благодаря простоте реализации и отсутствию необходимости в символьном дифференцировании.

Описание метода

Метод секущих основан на замене касательной, используемой в методе Ньютона, секущей — прямой, проходящей через две точки на графике функции, соответствующие двум последовательным приближениям к корню. Пересечение этой секущей с осью абсцисс даёт новое приближение.

Итерационная формула

Пусть требуется найти корень уравнения \( f(x) = 0 \). Задаются два начальных приближения \( x_0 \) и \( x_1 \). Тогда последующие приближения вычисляются по формуле:

\[ x_{n+1} = x_n - f(x_n) \cdot \frac{x_n - x_{n-1}}{f(x_n) - f(x_{n-1})} \]

где \( n = 1, 2, 3, \dots \). Геометрически это означает, что через точки \( (x_{n-1}, f(x_{n-1})) \) и \( (x_n, f(x_n)) \) проводится прямая, и её пересечение с осью \( x \) даёт \( x_{n+1} \).

Геометрическая интерпретация

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

Сходимость и особенности

Скорость сходимости

Метод секущих обладает сверхлинейной сходимостью. Порядок сходимости \( p \) приблизительно равен золотому сечению:

\[ p = \frac{1 + \sqrt{5}}{2} \approx 1.618 \]

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

Условия сходимости

Для гарантии сходимости метода секущих необходимо, чтобы функция \( f(x) \) была непрерывна и имела непрерывную вторую производную в окрестности корня, а также чтобы начальные приближения \( x_0 \) и \( x_1 \) были достаточно близки к истинному корню. Если начальные точки выбраны неудачно, метод может расходиться или зацикливаться. Кроме того, метод может столкнуться с проблемой деления на ноль, если \( f(x_n) = f(x_{n-1}) \), что означает, что секущая параллельна оси абсцисс.

Преимущества и недостатки

Преимущества:

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

Недостатки:

• Необходимость задания двух начальных приближений, что может быть нетривиально.
• Возможна потеря устойчивости при близких значениях функции в двух точках.
• Метод может не сходиться, если функция имеет особенности (разрывы, точки перегиба) в окрестности корня.

Алгоритм

Алгоритм метода секущих для решения уравнения \( f(x) = 0 \) можно описать следующим образом:

1. Задать начальные приближения \( x_0 \) и \( x_1 \), а также требуемую точность \( \varepsilon \).
2. Вычислить \( f(x_0) \) и \( f(x_1) \).
3. Для \( n = 1, 2, \dots \) до достижения сходимости:

• Вычислить \( x_{n+1} \) по итерационной формуле.
• Вычислить \( f(x_{n+1}) \).
• Проверить условие остановки: если \( |x_{n+1} - x_n| < \varepsilon \) или \( |f(x_{n+1})| < \varepsilon \), то остановиться.
• Обновить значения: \( x_{n-1} = x_n \), \( x_n = x_{n+1} \).

1. Вывести \( x_{n+1} \) как приближённое значение корня.

Пример применения

Рассмотрим уравнение \( f(x) = x^3 - 2x - 5 = 0 \). Известно, что корень находится вблизи \( x = 2 \). Выберем начальные приближения \( x_0 = 2 \) и \( x_1 = 3 \).

Итерации:

• \( x_0 = 2 \), \( f(2) = -1 \)
• \( x_1 = 3 \), \( f(3) = 16 \)
• \( x_2 = 3 - 16 \cdot \frac{3 - 2}{16 - (-1)} = 3 - 16 \cdot \frac{1}{17} \approx 2.0588 \), \( f(2.0588) \approx -0.390 \)
• \( x_3 = 2.0588 - (-0.390) \cdot \frac{2.0588 - 3}{-0.390 - 16} \approx 2.0813 \), \( f(2.0813) \approx -0.147 \)
• \( x_4 \approx 2.0948 \), \( f(2.0948) \approx -0.054 \)
• \( x_5 \approx 2.0946 \), \( f(2.0946) \approx 0.000 \)

После нескольких итераций метод сходится к корню \( x \approx 2.0946 \) с точностью до \( 10^{-4} \).

Сравнение с другими методами

Метод Ньютона

Метод Ньютона требует вычисления производной \( f'(x) \) на каждом шаге, что даёт квадратичную сходимость. Однако если производная сложна или её вычисление затратно, метод секущих может быть эффективнее. В методе Ньютона на каждой итерации вычисляются \( f(x_n) \) и \( f'(x_n) \), тогда как в методе секущих — только одно новое значение функции (после первой итерации).

Метод бисекции

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

Метод хорд

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

Применение

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

Интересные факты

• Метод секущих является частным случаем метода Бройдена, который используется для решения систем нелинейных уравнений и относится к классу квазиньютоновских методов.
• В некоторых источниках метод секущих называют «методом хорд», хотя в классической литературе метод хорд (regula falsi) отличается тем, что одна из точек остаётся неподвижной.
• Порядок сходимости метода секущих (\( \approx 1.618 \)) совпадает с золотым сечением, что является редким свойством среди численных методов.

Источники

• Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. — М.: Бином. Лаборатория знаний, 2008.
• Калиткин Н. Н. Численные методы. — М.: Наука, 1978.
• Самарский А. А., Гулин А. В. Численные методы. — М.: Наука, 1989.
• Burden R. L., Faires J. D. Numerical Analysis. — 9th ed. — Brooks/Cole, 2011.