Метод Ньютона решения систем нелинейных уравнений

 

Для решения системы (4.3) будем пользоваться методом последовательных приближений.

Предположим, известно k-е приближение

=

одного из изолированных корней векторного уравнения (4.2). Тогда точный корень уравнения (4.2) можно представить в виде

, (4.9)

где - поправка (погрешность корня).

Подставляя выражение (4.9) в (4.2), имеем

= 0. (4.10)

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

, (4.11)

или в развернутом виде,

Из формул (4.11) и (4.12) видно, что под производной следует понимать матрицу Якоби системы функций f1, f2, ..., fn относительно переменных x1, x2, ..., xn, т. е.

= W(x) =,

или в краткой записи

= W(x) = (i, j = 1, 2, …, n),

поэтому формула (4.12) может быть записана в следующем виде:

=0. (4.13)

Если det =, то

. (4.14)

Отсюда видно, что метод Ньютона решения системы (4.1) состоит в построении итерационной последовательности:

, (4.15)

где k = 0, 1, 2, ….

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

Пример 4.1.Найти методом Ньютона приближенное положительное решение системы уравнений

исходя из начального приближения x0 = y0 = z0 =0,5.

Полагая:

х(0) =, f (х) =,

имеем:

f (х) =

Отсюда

f ( х(0) ) =

Составим матрицу Якоби

W(x) =

Имеем

= ,

причем

D = det =

Следовательно, матрица – неособенная. Вычисляем обратную ей матрицу

=

По формуле (4.15) получаем первое приближение

=

= - = + = .

Аналогично находятся дальнейшие приближения. Результаты вычислений приведены в табл. 4.1.

 

Таблица 4.1