VI. Решения нелинейных уравнений.

§1. Постановка задачи. Основные определения. Рассмотрим нелинейное уравнение f(x)=0. Корнем этого уравнения называется значение , при котором. Решение уравнения заключается в нахождение его корней. Кореньназывается простым, если . Корень называется кратным, если . Целое числоm называется кратностью корня , еслидляk=1,2,3...(m-1), а . Случайk=1 соответствует простому корню.

Рассмотрим график некоторой функции y=f(x) (x[a,b]), который представлен на рис.1. Из определения следует, что корень является простым, если график функции y=f(x) пересекает ось 0x в точке подyглом 0 , и кратным, если он касается оси 0x в точке , т.к. имеем=0.

Существуют следующие методы решения нелинейных уравнений: аналитические, графические и численные. Аналитическими методами удается воспользоваться только для уравнений определенного вида, в общем случае они не применимы. Графические методы естественно обладают большой погрешностью. Поэтому основными методами являются численные методы.

Численное решение задачи нахождения корней нелинейного уравнения осуществляется в два этапа: этапа локализации корней и этапа итерационного уточнения корней.

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

Для локализации используются различные методы: аналитические, графические, таблицы. Аналитические и графические методы применяются для простых уравнений, например, для уравнения: . Для более сложных уравнений строятся таблицы и определяются значенияx_i и x_l+I, при которых функция y=f(x) меняет знак (поиск простых корней), или производные меняют знак (поиск кратных корней). Отсюда сразу можно сделать вывод, что задача нахождения простых корней существенно проще, чем задача отыскания кратных корней.

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

В основе алгоритма лежит итерационная формула (название происходит от латинского слова “iteracio”-повторение). Для нахождения корня с точностью  используется та или иная итерационная формула. Итерационный метод называется одношаговым, если для вычисления очередного приближения x_n используется только x_n-1 приближение, и k–шаговым, если используются k предыдущих приближений: x_n-1,x_n-2,...,x_n-k.

Критерий сходимости. Для сходимости итерационного процесса необходимо и достаточно выполнение следующего условия:

(6.1)

где c и p - некоторые константы, число p называется порядком сходимости метода.

При p=1 и с<1, метод обладает линейной сходимостью. Если p>1, то имеем сверхлинейную сходимость, если p=2, то имеем квадратичную сходимость. Если для всех n выполняется условие: , гдеq<1, то говорят о сходимости со скоростью геометрической прогрессии.

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

Для каждого итерационного процесса существует свой критерий окончания в виде неравенства: где()-функция от , при выполнение которого всегда имеет место неравенство: .

§2. Метод Ньютона (метод касательных). Пусть корень локализован на отрезке [a,b]. Выбирается точка x₀ на этом отрезке и в этой точке строится касательная к графику функции y=f(x). За новое приближение x₁ принимается точка, в которой касательная пересекает ось 0x. Далее процесс повторяется. Таким образом, на каждом этапе решение исходного уравнения заменяется на решение линейного уравнения (уравнения касательной), т.е. проводится линеаризация исходного уравнения.

Запишем уравнение касательной в точке x_n:

(6.2)

Полагая в, что при x=x_n+1 имеем y=0 получаем формулу метода Ньютона:

(6.3)

Критерий окончания метода Ньютона:

(6.4)