- •1. Чисельні методи та обчислювальна математика
- •1.1. Роль математичного моделювання у розв’язанні задач навколишнього світу
- •1.2. Структура похибки розв’язку задачі
- •1.3. Поняття стійкості та коректності
- •1.4. Основні етапи розв’язання прикладних математичних задач
- •1.5. Контрольні питання
- •2. Елементарна теорія похибок
- •2.1. Джерела і класифікація похибок
- •2.2. Абсолютна та відносна похибки
- •2.3. Вірні значущі цифри
- •2.4. Похибки арифметичних операцій
- •2.5. Графи обчислювальних процесів
- •2.6. Контрольні питання
- •2.7. Завдання
- •Т Варіант Варіант Варіант аблиця 2.1
- •Порядок розв’язання задачі
- •Т Варіант Варіант Варіант аблиця 2.2
- •3. Алгебра матриць
- •3.1. Визначник матриці, властивості визначника і правила його обчислення
- •3.2. Обернена матриця
- •3.3. Трикутні матриці
- •3.4. Обернення матриць за допомогою розбиття її на клітини
- •3.6. Поняття рангу матриці
- •3.7. Контрольні питання
- •3.8. Завдання
- •4. Розв’язання систем лінійних алгебраїчних рівнянь
- •4.1. Теорема Кронекера – Капеллі
- •4.2. Розв’язання довільних систем лінійних алгебраїчних рівнянь
- •4.3. Класифікація методів розв’язання систем лінійних алгебраїчних рівнянь та їх порівняльні характеристики
- •4.4. Точні методи розв’язання систем лінійних алгебраїчних рівнянь
- •4.4.1. Метод Крамера
- •4.4.2. Метод Гаусса
- •4.4.3. Метод головних елементів
- •4.4.4. Метод Гаусса–Жордана
- •4.4.5. Схема Халецького
- •4.4.6. Метод квадратних коренів
- •4.5. Ітераційні методи розв’язання систем лінійних алгебраїчних рівнянь
- •4.5.1. Метод ітерацій
- •4.5.2. Умови збіжності ітераційного процесу
- •4.5.3. Оцінювання похибки методу ітерацій
- •4.5.4. Метод Зейделя
- •4.5.5. Умови збіжності процесу Зейделя
- •4.5.6. Приведення системи лінійних рівнянь до вигляду, зручного для ітерацій
- •4.6. Контрольні питання
- •4.7. Завдання
- •Література
4.5.5. Умови збіжності процесу Зейделя
Процес
Зейделя для лінійної системи
так
, як і процес послідовних наближень,
збігається
до єдиного розв’язку при будь–якому
виборі початкового наближення, якщо
хоч би одна із норм матриці
менше одиниці.
Процес Зейделя часто збігається до єдиного розв’язку швидше ніж метод процесу простої ітерації.
Приклад 2. Перевірити чи збігається процес Зейделя для системи, яка розглянута в прикладі 1.
Розв’язання
1) Після приведення системи до нормального виду отримаємо матрицю
2)
Знайдемо
Отже, процес ітерації для даної системи
збігається до єдиного розв’язку, не
дивлячись на те, що
Оцінювання похибки процесу Зейделя.
Нехай
дана лінійна система
.
Якщо
–
точні значення її невідомих, а
–е
наближення, обчислене за методом Зейделя,
то для оцінювання похибки цього методу
застосовується формула
(4.57)
Приклад. Підрахувати, скільки ітерацій за методом Зейделя необхідно виконати, щоб з точністю до 10–4 знайти розв’язок системи
Розв’язання. 1) Приведемо систему до нормального вигляду:
2)
За нульові наближення приймемо стовпчик
вільних членів
і
обчислюємо перші наближення:
3) Матриця
Отже,
Оскільки
і
то
4) за формулою (4.57) визначаємо k:
аналогічно
можна виконувати оцінювання методу
Зейделя за нормою
.
4.5.6. Приведення системи лінійних рівнянь до вигляду, зручного для ітерацій
Процеси послідовних наближень і Зейделя для лінійної системи збігається до єдиного розв’язку незалежно від вибору початкового вектора, якщо
Отже,
для збіжності вищевказаних ітераційних
процесів достатньо, щоб значення
елементів
матриці
при
були невеликими за абсолютною величиною.
Це рівнозначно тому, що якщо для лінійної
системи
модулі діагональних коефіцієнтів
кожного рівняння системи більше суми
модулів всіх наступних коефіцієнтів
(не враховуючи вільних членів), то
ітераційні процеси для цієї системи
збігаються, тобто якщо дана система
,
причому
,
то процеси простої ітерації і Зейделя
для цієї системи збігаються [4].
Застосовуючи елементарні перетворення, лінійну систему можна замінити такою еквівалентною системою , для якої умови збіжності будуть виконані.
Приклад. Привести до виду, зручного для ітерацій, таку систему
Розв’язання
1) із заданої системи виділяємо рівняння з коефіцієнтами, модулі яких більші суми модулів наступних коефіцієнтів системи. Кожне виділене рівняння вписуємо в такий рядок нової системи, щоб найбільший за модулем коефіцієнт виявився діагональним.
В
рівнянні (Б) коефіцієнт при
за модулем більше суми модулів наступних
коефіцієнтів. Приймаємо рівняння (Б)
за друге рівняння нової системи:
2) із залишених невикористаних рівнянь системи утворюємо лінійно незалежні між собою комбінації. Так, за перше рівняння нової системи можна взяти лінійну комбінацію (2В)+(А), тоді матимемо
(I)
За третє рівняння нової системи можна прийняти лінійну комбінацію (2А)–(Б), тобто
(ІІІ)
3)
В результаті отримаємо перетворену
систему лінійних рівнянь (І), (ІІ), (ІІІ),
еквівалентну вихідній і задовольняє
умови збіжності ітераційного процесу:
Приводячи цю систему до нормального вигляду, отримаємо
Тепер слід лише розв’язати систему одним з ітераційних методів.
