ПРИБЛИЖЕННОЕ РЕШЕНИЕ АЛГЕБРАИЧЕСКИХ И ТРАНСЦЕНДЕНТНЫХ УРАВНЕНИЙ
Отделение корней
Пусть дано уравнение
В дальнейшем в некоторых случаях нам понадобится существование и непрерывность первой производной f ' (х) или даже второй производной f " (х), что будет оговорено в соответствующих местах.
Всякое значение ξ, обращающее функцию f (x) в нуль, т. е. такое, что
Мы будем предполагать, что уравнение (1) имеет лишь изолированные корни, т. е. для каждого корня уравнения (1) существует окрестность, не содержащая других корней этого уравнения.
Приближенное нахождение изолированных действительных корней уравнения (1) обычно складывается из двух этапов:
- 1) отделение корней, т.е. установление возможно тесных промежутков [α, β], в которых содержится один и только один корень уравнения (1);
- 2) уточнение приближенных корней, т.е. доведение их до заданной степени точности.
Теорема 1. Если непрерывная функция f (x) принимает значения разных знаков на концах отрезка [α, β], т. е. f(α)·f (β) < 0, то внутри этого отрезка содержится по меньшей мере один корень уравнения f (x) = 0, т. е. найдется хотя бы одно число ξ ∈ (α, β) такое, что f (ξ) = 0 (смотри рисунок.).
Корень ξ заведомо будет единственным, если производная f ' (x) существует и сохраняет постоянный знак внутри интервала (α, β), т. е. если f '(x) > 0 (или f '(x) < 0) при α < x < β (смотри рисунок.).
Процесс отделения корней начинается с установления знаков функции f (x) в граничных точках х = α и х = β области её существования.
Затем определяются знаки функции f (x) в ряде промежуточных точек x = α1, α2, ..., выбор которых учитывает особенности функции f (x). Если окажется, что f (αk)·f (αk+1) < 0, то в силу теоремы 1 в интервале (αk, αk+1) имеется корень уравнения f (x) = 0. Нужно тем или иным способом убедиться, является ли этот корень единствешшм. Для отделения корней практически часто бывает достаточно провести процесс половинного делення, приближенно деля данный интервал (α, β) на две, четыре, восемь и т. д. равных частей (до некоторого шага) и определяя знаки функции f (x) в точках делений. Полезно помнить, что алгебраическое уравнение n-й степени
Пример 1. Отделить корни уравнения
| x | f(x) | x | f(x) |
| - ∞ | - | 1 | - |
| -3 | - | 3 | + |
| -1 | + | + ∞ | + |
| 0 | + |
Пример 2. Отделить корни уравнения
Имеем f (- ∞) > 0; f (1) < 0 (-); f (+ ∞) > 0 (+). Следовательно, уравнение (3) имеет только два действительных корня, из которых один лежит в интервале (- ∞, 1), а другой - в интервале (1, + ∞).
Пример 3. Определить число действительных корней уравнения
Дадим теперь оценку погрешности приближенного корня.
Теорема 2. Пусть ξ - точный, а х* - приближенный корни уравнения f (x) = 0, находящиеся на одном и том же отрезке [α, β], причем | f ' (x) | ≥ m1 > 0 при α ≤ x ≤ β. В частности, за m1 можно взять наименьшее значение | f ' (x) | при α ≤ x ≤ β. В таком случае справедлива оценка
. (5)
Отсюда, так как f (ξ) = 0 и | f ' (c) | ≥ m1, получим:
.
Решение. Имеем f (x*) = 2,2153 - 1,22 - 1 = 0,0047. Так как при x = 1,23 получаем
Графическое решение уравнений
Действительные корни уравнения
Пример 5. Графически решить уравнение
.
Нахождение корней уравнения (6) упрощается, если одна из функций φ (х) или ψ(x) линейная, т. е., например, φ (х) = а х + b. В этом случае корни уравнения (6) находятся как абсциссы точек пересечения кривой у = φ(x) и прямой у = а х + b. Особенно выгодным оказывается этот прием при решении ряда однотипных уравнений, отличающихся только коэффициентами а и b линейной функции. Здесь графическое построение сводится к нахождению точек пересечения фиксированного графика y = ψ (x) различными прямыми. К указанному типу, очевидно, относятся трехчленные уравнения
Отметим, что хотя графические методы решения уравнений весьма удобны и сравнительно просты, но они, как правило, применимы лишь для грубого определения корней. Особенно неблагоприятным в смысле потери точности является случай, когда линии пересекаются под очень острым углом и практически сливаются по некоторой дуге.
Метод половинного деления
Пусть дано уравнение
= 0, то
вляется корнем уравнения. Если
≠ 0, то выбираем ту из половин
или
, на концах которой функция f (x) имеет противоположные знаки. Новый суженный отрезок [al , b1] снова делим пополам и проводим то же рассмотрение и т. д. В результате получаем на каком-то этапе или точный корень уравнения (1), или же бесконечную последовательность вложенных друг в друга отрезков [a1, b2], [а2, b2], ..., [аn, bn], ... таких, что
. (9)
.
. (10)
Метод половинного деления практически удобно применять для грубого нахождения корня данного уравнения, так как при увеличении точности значительно возрастает объем вычислительной работы.
Заметим, что метод половинного деления легко реализуется на электронных счетных машинах. Программа вычисления составляется так, чтобы машина находила значение правой части уравнения (1) в середине каждого из отрезков [аn, bn] ( n = 1, 2, .. .) и выбирала соответствующую половину его.
Пример 7. Методом половинного деления уточнить корень уравнения
Решение. Последовательно имеем:
f (0,5) = 0,06 + 0,25 - 0,5 - 1= - 1,19;
f (0,75) = 0,32 + 0,84 - 0,75 - 1 = - 0,59;
f (0,875) = 0,59 + 1,34 - 0,88 - 1 = + 0,05;
f (0,8125) = 0,436 + 1,072 - 0,812 - 1 = - 0,304;
f (0,8438) = 0,507 + 1,202 - 0,844 - 1 = - 0,135;
f (0,8594) = 0,546 + 1,270 - 0,859 - 1 = - 0,043 и т.д.
Способ пропорциональных частей (метод хорд)
Укажем более быстрый способ нахождения корня ξ уравнения f (x) = 0, лежащего на заданном отрезке [а, b] таком, что f (a)·f (b) < 0.Пусть для определенности f(a) < 0 и f (b) > 0. Тогда, вместо того чтобы делить отрезок [а, b] пополам, более естественно разделить его в отношении - f (a) : f (b). Это дает нам приближенное значение корня
. (12)
Геометрически способ пропорциональных частей эквивалентен замене кривой у = f (х) хордой, проходящей через точки A(a, f (a)] и B [b, f (b)] (смотри рисунок). В самом деле, уравнение хорды АВ есть
(1')
Для доказательства сходимости процесса предположим, что корень отделен и вторая производная f "(x) сохраняет постоянный знак на отрезке [а, b].
Пусть для определенности f " (х) > 0 при а ≤ х ≤ b (случай f " (x) < 0 сводится к нашему, если записать уравнение в виде - f (x) = 0). Тогда кривая у = f (х) будет выпукла вниз и, следовательно, расположена ниже своей хорды АВ. Возможны два случая: 1) f (а) > 0 (смотри рисунок) и 2) f (а) < 0 (смотри рисунок).
В первом случае конец а неподвижен и последовательные приближения: х0 = b;
(13)
(14)
Совершенно так же переходом к пределу в равенстве (14) доказывается, что ξ* = ξ для второго случая.
Для оценки точности приближения можно воспользоваться формулой (5)
,Приведем еще формулу, позволяющую оценивать абсолютную погрешность приближенного значения хn, если известны два последовательных приближения хn-1 и хn.
Будем предполагать, что производная f ' (х) непрерывна на отрезке [а, b], содержащем все приближения, и сохраняет постоянный знак, причем
.
. (16)
, (17)
Решение. Прежде всего отделяем корень. Так как

f (x1) = - 0,173;

f (x2) = - 0,036;

f (x2) = - 0,0072;
.
Метод Ньютона (метод касательных)
.
. (19)
Выберем, например, х0 = b, для которого f (x0)·f "(x0) > 0. Проведем касательную к кривой y = f (x) в точке В0[х0, f(x0)]. В качестве первого приближения х1 корня ξ возьмем абсциссу точки пересечения этой касательной с осью Ох. Через точку B1[x1, f (x1)) снова проведем касательную, абсцисса точки пересечения которой даст нам второе приближение х2 корня ξ и т. д.. Очевидно, что уравнение касательной в точке Вn [хn, f(xn)] (n = 0, 1, 2, ...) есть
.
Теорема. Если f (a)·f (b) < 0, причем f '(x) u f "(x) отличны от нуля и сохраняют определенные знаки при а ≤ х ≤ b, то, исходя из начального приближения х0 ∈ [а, b], удовлетворяющего неравенству (20), можно вычислить методом Ньютона (формула (19)) единственный корень ξ уравнения (1) с любой степенью точности.
Доказательство. Пусть, например, f (а) < 0, f (b) > 0, f '(x) > 0, f "(х) > 0 при a ≤ x ≤ b (остальные случаи рассматриваются аналогично). Согласно неравенству (20) имеем f (x0) > 0 (например, можно принять х0 = b).
Методом математической индукции докажем, что все приближения xn > ξ (n = 0, 1, 2, ...) и, следовательно, f (xn) > 0. В самом деле, прежде всего, х0 > ξ.
Пусть теперь хn > ξ. Положим
, (21)
Так как f "(х) > 0, то имеем:
,
Из формулы (19), учитывая знаки f (xn) и f '(хn), имеем хn+1 < хn (n = 0, 1, ...), т. е. последовательные приближения x0, xl, ..., хn, ... образуют ограниченную монотонно убывающую последовательность. Следовательно, существует
.Переходя к пределу в равенстве (19), будем иметь:
.
Поэтому, применяя метод Ньютона, следует руководствоваться следующим правилом: в качестве исходной точки х0 выбирается тот конец интервала (а, b), которому отвечает ордината того же знака, что и знак f " (х).
Замечание 1. Если: 1) функция f (x) определена и непрерывна при - ∞ < x < ∞; 2) f (a)·f (b) < 0; 3) f ' (х) ≠ 0 при а ≤ х ≤ b; 4) f "(х) существует всюду и сохраняет постоянный знак, то при применении метода Ньютона для нахождения корня уравнения f (х) = 0, лежащего в интервале (а, b), за начальное приближение х0 можно принять любое значение с ∈ [а, b]. В частности, можно взять х0 = а или х0 = b.
Действительно, пусть, например, f '(x) > 0 при a ≤ x ≤ b, f "(х) > 0 и х0 = с, где a ≤ c ≤ b.
Если f (с) = 0, то корень ξ = с и задача, таким образом, решена.
Если f (с) > 0, то справедливо приведенное выше рассуждение и процесс Ньютона с начальным значением с сходится к корню ξ ∈ (a, b).
Наконец, если f (с) < 0, то находим значение
.
,
Таким образом,
Аналогичное рассмотрение можно провести для других комбинаций знаков производных f ' (х) и f " (х).
Замечание 2. Из формулы (19) видно, что чем больше численное значение производной f ' (х) в окрестности данного корня, тем меньше поправка, которую нужно прибавить к n - му приближению, чтобы получить (n + 1) - е приближение. Поэтому метод Ньютона особенно удобно применять тогда, когда в окрестности данного корня график функции имеет большую крутизну. Но если численное значение производной f '(x) близ корня мало, то поправки будут велики, и вычисление корня по этому методу может оказаться очень долгим, а иногда и вовсе невозможным. Следовательно, если кривая y = f (x) вблизи точки пересечения с осью Ох почти горизонтальна, то применять метод Ньютона для решения уравнения f {x) = 0 не рекомендуется.
Для оценки погрешности n-го приближения хn можно воспользоваться общей формулой (5)
,
Выведем еще одну формулу для оценки точности приближения хn. Применяя формулу Тейлора, имеем:
,(23)
(24)
Заметим, что в общем случае совпадение с точностью доедвух последовательных приближений xn-1 и xn вовсе не гарантирует, что с той же точностью совпадает значение xn и точный корень ξ (смотри рисунок).
Установим также формулу, связывающую абсолютные погрешности двух последовательных приближений хn и xn+1. Из формулы (21) получаем:
. (25)
.
и | ξ - xn | < 10-m,
Пример 1. Вычислить методом Ньютона отрицательный корень уравнения f (x) = x4 - 3 x2 + 75 x - 10 000 = 0 с пятью верными знаками.
Решение. Полагая в левой части уравнения х = 0, -10, - 100, , получим f (0) = - 10 000, f (- 10) = - 1050, f (- 100) ≈ + 108.
Следовательно, искомый корень ξ находится в интервале - 100 < ξ < - 10. Сузим найденный интервал. Так как f (- 11) = 3453, то - 11 < ξ < - 10. В этом последнем интервале f '(х) < 0 и f "(х) > 0. Так как f ( - 11) > 0 и f "( - 11) > 0, то можем принять за начальное приближение х0 = - 11. Последовательные приближения хn (n = 1, 2, ...) вычисляем по следующей схеме:
| n | xn | f (xn) | f ' (xn) | ![]() |
| 0 | - 11 | 3453 | - 5183 | 0,7 |
| 1 | - 10,3 | 134,3 | - 4234 | 0,03 |
| 2 | - 10,27 | 37,8 | - 4196 | 0,009 |
| 3 | - 10,261 | 0,2 | - | - |
Комбинированный метод
Отсюда, в частности, вытекает, что цифры, общие для хn и хn*, обязательно принадлежат точному корню ξ. Теоретически здесь возможны четыре случая:
- 1) f '(х) > 0; f "(x) > 0;
- 2) f '(х) > 0; f "(x) < 0;
- 3) f '(х) < 0; f "(x) > 0;
- 4) f '(х) < 0; f "(x) < 0.
