§ 2. Методы уточнения корней

1. Метод половинного деления (бисекции, дихотомии)

Считаем, что отделение корней уравнения (2.1) проведено и на отрезке расположен один корень, который необходимо уточнить с погрешностьюe. В качестве начального приближения корня принимаем середину этого отрезка: (рис. 2.5). Затем исследуем значение функции на концах отрезков и . Тот из отрезков, на концах которого принимает значения разных знаков, содержит искомый корень; поэтому его принимаем в качестве нового отрезка (на рис. 2.5 это отрезок ). Вторую половину отрезка , на которой не меняет знак, отбрасываем. В качестве следующего приближения корня принимаем середину нового отрезка и т.д. Таким образом,k-е приближение вычисляется как

. (2.5)

После каждой итерации отрезок, на котором расположен корень, уменьшается вдвое, а после kитераций в раз:

. (2.6)

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

. (2.7)

Поскольку корень принадлежит отрезку , а – середина этого отрезка, то величина всегда будет меньше половины длины отрезка (см. рис. 2.5), т.е.

. (2.8)

Следовательно, условие (2.7) будет выполнено, если

. (2.9)

Таким образом, итерационный процесс нужно продолжать до тех пор, пока не будет выполнено условие (2.9).

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

Однако метод половинного деления довольно медленный. С каждым шагом погрешность приближенного значения уменьшается в два раза, т.е.

, (2.10)

поэтому данный метод является методом с линейной сходимостью.

Вычислим количество итераций N, требуемое для достижения заданной точностиe. Пользуясь выражением (2.6) можно выяснить для каких значенийkбудет выполнено условие (2.9), и взять в качествеNнаименьшее из такихk:

, , (2.11)

где – целая часть числаx. Например, при и получим .

ØЗамечание. При реализации метода следует учитывать, что функция вычисляется с некоторой абсолютной погрешностью . Вблизи корня значения функции малы по абсолютной величине и могут оказаться сравнимы с погрешностью ее вычисления. Другими словами, при подходе к корню мы можем попасть в полосу шумов и дальнейшее уточнение корня окажется невозможным. Поэтому целесообразно задать ширину полосы шумов и прекратить итерационный процесс при попадании в нее. Если принять , то итерационный процесс можно завершать, когда значение функции послеk-й итерации станет меньшим по модулюe., т.е.

. (2.12)

Также необходимо иметь ввиду, что при уменьшении интервала увеличиваются погрешности вычисления его длины за счет вычитания близких чисел.<^3[3]

2. Метод хорд

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

Запишем уравнение прямой, проходящей через точки АиВ:

.

Для точки пересечения прямой с осью абсцисс () получим уравнение

. (2.13)

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

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

(2.14)

или при выполнении условия (2.12).

ØЗамечание. Метод половинного деления и метод хорд очень похожи, в частности, процедурой проверки знаков функции на концах отрезка. При этом второй их них в ряде случаев дает более быструю сходимость итерационного процесса. Однако в некоторых случаях метод хорд может сходится существенно медленнее метода половинного деления. Такая ситуация показана на рис. 2.7. Оба рассмотренных метода не требуют знания дополнительной информации о функции . Например, не требуется, чтобы функция была дифференцируема. Даже для разрывных функций рассмотренные методы обладают гарантированной сходимостью. Более сложные методы уточнения корня используют дополнительную информацию о функции , прежде всего свойство дифференцируемости. Как результат они обычно обладают более быстрой сходимостью, но в то же время, применимы для более узкого класса функций, и их сходимость не всегда гарантирована. Примером такого метода служит метод Ньютона.<