МЕТОДЫ РЕШЕНИЯ СИСТЕМЫ НЕЛИНЕЙНЫХ УРАВНЕНИЙ

Постановка задачи

Пусть требуется решить систему n нелинейных уравнений:

(1)

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

Систему уравнений (1) можно кратко записать в векторном виде:

. (2)

Уравнение (2) может иметь один или несколько корней в области определения D. Требуется установить существование корней уравнения и найти приближённые значения этих корней. Для нахождения корней обычно применяют итерационные методы, в которых принципиальное значение имеет выбор начального приближения. Начальное приближение иногда известно из физических соображений. В случае двух неизвестных начальное приближение можно найти графически: построить на плоскости (x₁, x₂) кривые f₁(x₁, x₂)=0 и f₂(x₁, x₂)=0 и найти точки их пересечения. Для трех и более переменных (а также для комплексных корней) удовлетворительных способов подбора начального приближения нет.

Рассмотрим два основных итерационных метода решения системы уравнений (1), (2) - метод простой итерации и метод Ньютона.

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

Метод простой итерации

Представим систему (1) в виде

(3)

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

(4)

Алгоритм метода простой итерации состоит в следующем. Выберем некоторое нулевое приближение

Следующее приближение находим по формулам:

или более подробно:

(5)

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

На практике часто вместо последнего условия используют неравенство:

(6)

где - среднеквадратичная норма n-мерного вектора , т.е.

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

Метод Ньютона

В переводной литературе можно встретить название метод Ньютона-Рафсона. Этот метод обладает гораздо более быстрой сходимостью, чем метод простой итерации.

Пусть известно некоторое приближение к корню , так что

Тогда исходную систему (2) можно записать следующим образом:

Разлагая уравнение (7) в ряд Тейлора в окрестности точки и ограничиваясь линейными членами по отклонению , получим:

или в координатной форме:

(8)

Систему (8) можно переписать в виде:

(9)

Полученная система (9) является системой линейных алгебраических уравнений относительно приращений

Значение функций F₁, F₂, …, F_n и их производные в (9) вычисляются при

Определителем системы (9) является якобиан J:

(10)

Для существования единственного решения системы уравнений (9) он должен быть отличен от нуля. Решив систему (9), например, методом Гаусса, найдём новое приближение:

Проверяем условие (6). Если оно не удовлетворяется, находим и якобиан (10) с новым приближением и опять решаем (9), таким образом, находим 2-е приближение и т.д.