6.3.5. Система нелінійних рівнянь

На відміну від СЛАР для систем нелінійних рівнянь (СНР) не існує прямих методів розв’язку, а тому їх розв’язують лише ітераційними способами. В загальній формі СНР записується так:

, (6.3)

або в векторній формі , деХ – вектор невідомих, – вектор-функція.

Для одержання ітераційної формули створення процесу прямої ітерації приведемо систему (6.3) до вигляду:

що в векторній формі записують .

Задавши початкове значення вектора невідомих , одержуємо ітераційний процес:і т.д., допоки різниця норм векторівХ на сусідніх ітераціях не стануть менше наперед заданого малого числа :

.

Недоліком таких процесів є те, що початкове значення потрібно вибирати лише взоні збіжності поблизу точки розв’язку. Цю зону визначають із фізичних властивостей процесів чи об’єктів, режим роботи яких описується даною системою нелінійних рівнянь. В двовимірному просторі зону збіжності можна визначати графічно.

Наступним обмеженням застосування метода прямої ітерації є те, що для створення збіжного ітераційного процесу перетворення потрібно здійснити таким чином, щоб

Наведемо приклад розв’язку СНР для двох невідомих:

Найбільшого поширення в інженерній практиці при розв’язуванні СНР набув метод Ньютона, який має ряд переваг перед іншими ітераційними методами. В першу чергу – це його значна швидкість збіжності. Для побудови ітераційного процесу за цим методом вектор-функцію розкладемо вn-вимірному просторі в ряд Тейлора:

Згідно (6.3) . Далі в цьому ряду

–вектор-функція, розміщена в зоні збіжності;

–матриця Якобі, яка складається із елементів, що являють собою частинні похідні від усіх рівнянь системи по усім невідомим;

 – вектор-нев’язка, яка наближає вектор до точки розв’язку системи;– матриця Гессе, що складається із частинних похідних другого порядку. В двовимірному випадку маємо ці вектори на рисунку 26.

Рисунок 26 – Геометрична інтерпретація метода Ньютона

Враховуючи ітераційний процес, залишимо в ряду лише два перших елемента і одержимо значення вектора :

Звідси ітераційна формула буде мати вигляд:

або:

.

Таким чином ітераційний процес по методу Ньютона реалізуються схемою: і т.д.

Метод збігається до точки розв’язку дуже швидко (2-3 ітерації), але і для нього існує проблема вибору початкового вектора .

Застосуємо метод для попередньої задачі.

Перевіримо, розв’язок інструментарієм системи MathCad:

Варіанти завдань

Таблиця 13

№ вар	Нелінійне рівняння f(x)	№ вар	Нелінійне рівняння f(x)
1		14
2		15
3		16
4		17
5		18
6		19
7		20
8		21
9		22
10		23
11		24
12		25
13		26

Таблиця 14

№ вар	Система нелінійних рівнянь	№ вар	Система нелінійних рівнянь
1		14
2		15
3		16
4		17
5		18
6		19
7		20
8		21
9		22
10		23
11		24
12		25
13		26