Метод хорд

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

История

Идея метода восходит к античным математикам. Впервые систематическое описание метода ложного положения (regula falsi) встречается в древнеегипетском папирусе Ахмеса (около 1650 года до н. э.), где он применялся для решения линейных уравнений. В средневековой Европе метод был популяризирован благодаря работам Леонардо Фибоначчи (XIII век) и Луки Пачоли (XV век). В современной форме метод хорд как численный алгоритм для решения нелинейных уравнений был формализован в XVIII—XIX веках в трудах Леонарда Эйлера, Жозефа Луи Лагранжа и Огюстена Луи Коши. С развитием вычислительной техники в XX веке метод получил широкое распространение благодаря своей простоте и наглядности.

Принцип работы

Метод хорд основан на геометрической интерпретации: на заданном отрезке [a, b], где функция f(x) меняет знак (то есть f(a) * f(b) < 0), строится хорда — прямая, соединяющая точки (a, f(a)) и (b, f(b)). Точка пересечения этой хорды с осью абсцисс принимается за новое приближение корня x₁. Затем процесс повторяется на отрезке [a, x₁] или [x₁, b] в зависимости от знака функции в точке x₁.

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

Для отрезка [a, b] с условием f(a) * f(b) < 0 очередное приближение xₙ вычисляется по формуле:

`` xₙ = a - f(a) * (b - a) / (f(b) - f(a)) ``

или в симметричной форме:

`` xₙ = (a f(b) - b f(a)) / (f(b) - f(a)) ``

После вычисления xₙ проверяется знак f(xₙ):

• Если f(a) * f(xₙ) < 0, то корень находится на отрезке [a, xₙ], и b заменяется на xₙ.
• Если f(xₙ) * f(b) < 0, то корень находится на отрезке [xₙ, b], и a заменяется на xₙ.

Процесс продолжается до тех пор, пока длина отрезка не станет меньше заданной точности ε или пока |f(xₙ)| < ε.

Сходимость

Метод хорд обладает линейной сходимостью (скорость сходимости порядка 1), что медленнее, чем у метода Ньютона (квадратичная сходимость), но быстрее, чем у метода деления отрезка пополам (бисекции). Скорость сходимости зависит от формы функции: чем ближе функция к линейной на отрезке, тем быстрее сходимость. Для функций с большой кривизной сходимость может быть замедленной.

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

Для гарантированной сходимости метода хорд достаточно, чтобы:

1. Функция f(x) была непрерывна на отрезке [a, b].
2. На концах отрезка функция принимала значения разных знаков: f(a) * f(b) < 0.
3. Функция была монотонна на отрезке (не обязательно строго).

При нарушении монотонности метод может сходиться к другому корню или расходиться.

Разновидности метода

Метод ложного положения (regula falsi)

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

Модифицированный метод хорд

Для ускорения сходимости применяют модификации, например, метод хорд с поправкой (метод Хордена — Гросса), где на каждом шаге вычисляется не только точка пересечения хорды, но и корректируется положение одной из границ.

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

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

Применение

Метод хорд широко применяется в различных областях науки и техники:

• Инженерные расчёты: решение уравнений теплопроводности, гидравлики, механики.
• Физика: нахождение корней уравнений состояния, моделирование колебаний.
• Экономика: расчёт внутренней нормы доходности, точки безубыточности.
• Компьютерная графика: построение неявных поверхностей, расчёт пересечений.
• Обучение: используется в учебных курсах по численным методам как наглядный пример итерационного алгоритма.

Достоинства и недостатки

Достоинства

• Простота реализации и понимания.
• Не требует вычисления производной функции (в отличие от метода Ньютона).
• Гарантированная сходимость при выполнении условий.
• Устойчивость к вычислительным ошибкам.

Недостатки

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

Пример реализации

Пример реализации метода хорд на языке Python для решения уравнения f(x) = x³ - 2x - 5 = 0 на отрезке [2, 3]:

```python def f(x): return x**3 - 2*x - 5

def chord_method(a, b, eps=1e-6, max_iter=100): if f(a) f(b) >= 0: raise ValueError("Функция должна иметь разные знаки на концах отрезка") for i in range(max_iter): x = a - f(a) (b - a) / (f(b) - f(a)) if abs(f(x)) < eps: return x if f(a) * f(x) < 0: b = x else: a = x return x

root = chord_method(2, 3) print(f"Корень: {root:.6f}") ```

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

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

• Метод хорд является частным случаем метода линейной интерполяции, который используется не только в математике, но и в обработке сигналов, компьютерной графике и статистике.
• В некоторых источниках метод хорд называют «методом ложного положения» (regula falsi), хотя в современной литературе эти термины часто различают: regula falsi — это метод с фиксированной границей, а метод хорд — с переменной.
• Метод хорд был одним из первых алгоритмов, реализованных на механических вычислительных машинах в XIX веке.
• В российских учебниках по численным методам (например, в работах Н. С. Бахвалова, В. И. Лебедева) метод хорд традиционно рассматривается как базовый алгоритм для изучения итерационных процессов.

Источники

• Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. — М.: Бином, 2008.
• Самарский А. А., Гулин А. В. Численные методы. — М.: Наука, 1989.
• Калиткин Н. Н. Численные методы. — М.: Наука, 1978.
• Форсайт Дж., Малькольм М., Моулер К. Машинные методы математических вычислений. — М.: Мир, 1980.
• Демидович Б. П., Марон И. А. Основы вычислительной математики. — М.: Наука, 1970.