Метод Зейделя

В методі Зейделя уточнене значення х₁ зразу ж використовується для обчислення х₂, далі нові значення х₁ і х₂ використовуються для обчислення х₃ і т. д.

Це невелике удосконалення ітераційної процедури дозволяє суттєво збільшити швидкість збіжності.

Будь-яке (k +1)-е наближення в методі Зейделя будується за наступними формулами:

(3.17)

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

Ітерації закінчуються, коли із заданою точністю одержано однакові значення невідомих у двох ітераціях підряд.

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

Отже, умови збіжності можна записати так:

(i = 1, 2, …, n) або (j = 1, 2, …, n).

Як і в методі простої ітерації треба привести СЛАР до виду, який придатний для ітерацій. Для виконання умов збіжності ітераційного процесу достатньо, щоб значення елементів матриці при були невеликими з абсолютної величини. Це рівносильно тому, що якщо для СЛАР модулі діагональних коефіцієнтів кожного рівняння системи більше суми модулів всіх інших коефіцієнтів (без врахування вільних членів), то ітераційні процеси для цієї системи збігаються.

Окремо, на прикладі показується, як виконується еквівалентне перетворення вихідної СЛАР і отримується нормалізована система в загальному випадку.

Вихідна СЛАР:

Виконуються наступні дії:

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

б) інші рівняння нової еквівалентної системи одержуються шляхом складання лінійних незалежних між собою комбінацій. Так, за перше рівняння можна прийняти таку лінійну комбінацію (P)+(R), тоді:

.

За друге рівняння нової системи – таку комбінацію: (2Q)+(R)-(P), тобто

.

В результаті одержано перетворену СЛАР яка еквівалентна вихідній і задовольняє умовам збіжності ітераційного процесу:

Для перевірки цього твердження еквівалентна система приводиться до нормального виду і перевіряється, чи задовольняється хоч одна з умов збіжності.

При цьому використовується такий спосіб: записуються коефіцієнти при невідомих x₁, x₂, x₃ у відповідних рівняннях системи як m ·x, де m – число, що близьке до коефіцієнта при відповідному невідомому і на яке легко розділити коефіцієнти при невідомих і вільні члени.

Наприклад, приймається m = 10. Тоді система, що приводиться до нормального виду, перепишеться так:

або

Матриця і вектор приймають вид

, .

Суми модулів елементів рядків матриці :

0,24 + 0,05 + 0,24 = 0,53;

0,22 + 0,09 + 0,44 = 0,75;

0,18 + 0,25 + 0,54 = 0,97.

Більша із сум 0,97 < 1. Отже, одна з умов збіжності ітераційного процесу виконується. І хоч друга умова не виконується (більша сума модулів елементів стовпців 0,24 + 0,44 + 0,54 = 1,22 > 1) процес ітерації для системи, що розглядається, збігається до єдиного розв‘язку.

Приклад 2. Розрахувати струми в гілках електричного кола (рис. 3.1) методом Зейделя.

Розв'язується СЛАР (3.13) з прикладу 1.

Після еквівалентних перетворювань система має вид (3.14)

Надалі еквівалентна СЛАР приводиться до нормального виду, що придатний до проведення ітерацій, за допомогою коефіцієнта m= 20, аналогічно попередньому прикладу еквівалентних перетворювань, але з урахуванням матричної форми запису системи рівнянь.

(3.18)

З (3.18) видно, що

,

і умова збіжності ітераційного процесу виконується по рядках при .

Нульове наближення, що дорівнює з (3.17):

Задається абсолютна похибка розрахунку

Перше наближення згідно ітераційних формул методу (3.17):

Друге наближення:

Третє наближення:

Аналогічно виконуються кроки наступних ітерацій:

Перевірка умови закінчення ітераційного процесу після 9-го кроку:

За дев'ять кроків ітераційний процес закінчився з заданою точністю.

Струм в гілках схеми (рис. 1.1) за методом Зейделя становить:

Перевірка у вузлі „а” (рис. 1.1) за першим законом Кірхгофа виконується з точністю до |(0,249 + 1,566) – 1,832| = 0,017 А.

Дата добавления: 2015-05-26