- •Список прийнятих скорочень
- •Тема 1. Методи розв'язання систем лінійних рівнянь
- •Лекція 1. Метод Гауса
- •Концепція методів
- •Метод Гауса
- •Верхня трикутна система лінійних рівнянь
- •Метод виключення Гауса й вибір головного елемента
- •Схема єдиного ділення
- •Лекція 2. Ітераційні методи
- •Метод ітерацій
- •Зауваження про точність розрахунку
- •Достатня умова
- •Приведення лінійної системи до виду зручному для ітерації.
- •Метод Зейделя
- •Тема 2. Методи рішення нелінійних рівнянь
- •Лекція 3. Метод половинного ділення
- •Наближене рішення нелінійних рівнянь
- •Відділення корінь
- •Метод половинного ділення
- •Лекція 4. Метод Ньютона
- •Методика рішення задачі
- •Помилка ділення на нуль.
- •Швидкість збіжності.
- •Модифікації методу Ньютона.
- •Спрощений метод Ньютона
- •Метод Ньютона-Бройдена
- •Метод січних
- •Тема 3. Чисельне інтегрування
- •Лекція 5. Метод трапецій
- •Постановка задачі
- •Формула трапецій
- •Погрішність формули трапецій
- •Загальна формула трапецій
- •Лекція 6. Метод Сімпсона
- •Формула Сімпсона
- •Залишковий член формули Сімпсона
- •Загальна (узагальнена) формула Сімпсона
- •Тема 4. Обробка експериментальних даних
- •Лекція 7. Інтерполяція
- •Постановка задачі
- •Линейная інтерполяція
- •Квадратична інтерполяція
- •Інтерполяційна формула Лагранжа.
- •Обчислення Лагранжевых коефіцієнтів
- •Інтерполяція сплайном
- •Лекція 8. Метод найменших квадратів
- •Постановка задачі
- •Метод найменших квадратів
- •Лінійна апроксимація (інтерполяція)
- •Коефіцієнт лінійної кореляції
- •Квадратична апроксимація
- •Додатка
- •Транспонування
- •Обчислення визначника матриці
- •Знаходження зворотної матриці
- •Додавання й вирахування матриць
- •Множення матриці на число
- •Множення матриць
- •Ітераційні методи рішення рівнянь
- •Стандартні форми рівнянь
- •Пошук корінь графічним методом
- •Простий ітераційний метод здогаду й перевірки
- •Подання рівняння у формі 2
- •Пряма підстановка
- •Ітерації в осередку
- •Введення в надбудову Пошук рішення
- •Активування надбудови Пошук рішення
- •Установка надбудови Пошук рішення
- •Застосування надбудови Пошук рішення
- •Додаток 3. Контрольні питання
- •Додаток 4. Список лабораторних робіт
- •Частина 1. Обчислювальна техніка
- •Частина 2. Чисельні методи
- •Список літератури.
- •Основна література
- •Додаткова література
- •Інтернет-ресурси
x1 = −0,4 + 0 x1 − 0,2 x2 + 0,1 x3 − 0,2 x4 |
|
|
= 0,2 + 0,2 x1 + 0 x2 − 0,2 x3 + 0 x4 |
x2 |
|
|
= −0,4 + 0,2 x1 − 0,4 x2 + 0 x3 + 0,2 x4 |
x3 |
|
|
= −1,111 + 0,333 x1 + 0 x2 + 0 x3 + 0 x4 |
x4 |
до якої можна застосувати метод ітерацій.
Метод Зейделя
Метод Зейделя являє собою деяку модифікацію методу ітерацій. Основна ідея методу полягає в тому, що при обчисленні (k + 1)-го наближення
невідомої xi враховувати вже обчислені (k + 1)-е наближення невідомих
x1 , x2 ,...,xi−1 .
Нехай дана наведена система лінійних рівнянь:
n |
|
|
xi = βi + ∑αij xj |
(i = 1,2,...,n) |
(2.6) |
j=1
Вибираємо початкове наближення (розумне): x1(0) , x2(0) ,..., xn(0)
Далі припускаючи, що k -е наближення відомо, будемо будувати (k + 1)-е наближення корінь по наступних формулах:
x1(k+1)
x2(k+1)
|
n |
|
= β1 |
+ ∑αij x(jk ) |
|
|
j=1 |
|
|
+α21 x1(k+1) |
n |
= β2 |
+ ∑α2 j x(jk ) |
|
|
|
j=2 |
........................................ |
|
(2.7) |
|
i−1 |
n |
||
|
|||
xi(k+1) = βi + ∑αij x(jk+1) |
+ ∑αij x(jk ) |
|
|
j=1 |
j=i |
|
|
........................................ |
|
|
|
n−1 |
|
|
|
xn(k+1) = βn + ∑αnj x(jk+1) +αnn xn(k ) |
(k = 0,1,2,...) |
j=1
Або в розгорнутому виді:
x1(k+1) =α11 x1(k ) +α12 x2(k ) +α13 x3(k ) + ...+α1n xn(k ) + β1
x2(k+1) =α21 x1(k+1) +α22 x2(k ) +α23 x3(k ) + ...+α2n xn(k ) + β2
(2.8)
x3(k+1) =α31 x1(k+1) +α32 x2(k+1) +α33 x3(k ) + ...+α3n xn(k ) + β3
M
xn(k+1) =αn1 x1(k+1) +αn2 x2(k+1) +αn3 x3(k+1) + ...+αnn−1 xn(k−+11) +αnn xn(k ) + βn
20
Відзначимо, що зазначена вище теорема збіжності для методу ітерацій залишається вірної й для методу Зейделя.
Звичайно, але не завжди, метод Зейделя дає кращу збіжність Приклад 3. Методом Зейделя вирішити систему рівнянь:
10 x1 + x2 + x3 = 122 x1 + 10 x2 + x3 = 13
2 x1 + 2 x2 + 10 x3 = 14
Приведемо систему до виду, зручному для ітерацій:
x1 = 1,2 − 0,1 x2 − 0,1 x3x2 = 1,3 − 0,2 x1 − 0,1 x3x3 = 1,4 − 0,2 x1 −0,2 x2
Як нульові наближення корінь візьмемо:
x1(0) = 1,2 |
x2(0) = 0 |
x3(0) = 0 |
Застосовуючи метод Зейделя послідовно одержимо: Перший крок:
x1 = 1,2 − 0,1 0 − 0,1 0 = 1,2
x2 = 1,3 − 0,2 1,2 − 0,1 0 = 1,06
x3 = 1,4 − 0,2 1,2 − 0,2 1,06 = 0,948
Другий крок:
x1 = 1,2 − 0,1 1,06 − 0,1 0,948 = 0,9992
x2 = 1,3 − 0,2 0,9992 − 0,1 0,948 = 1,00536
x3 = 1,4 − 0,2 0,9992 − 0,2 1,00536 = 0,999098
і т.д.
Результати обчислень із точністю до чотирьох знаків наведені в таблиці 2.2
Таблиця 2.2 Обчислення рішення системи лінійних рівнянь методом Зейделя
k |
x1(k ) |
|
x2(k ) |
x3(k ) |
0 |
1,2000 |
|
0,0000 |
0,0000 |
1 |
1,2000 |
|
1,0600 |
0,9480 |
2 |
0,9992 |
|
1,0054 |
0,9991 |
3 |
0,9996 |
|
1,0001 |
1,0001 |
4 |
1,0000 |
|
1,0000 |
1,0000 |
5 |
1,0000 |
|
1,0000 |
1,0000 |
Точні значення корінь: x1 = 1 , |
x2 = 1, x3 = 1 |
|
Приклад розрахунку представлений на рис. 2.1
21
22
Рис. 2.1
Приклад
розрахунку по методу Зейделя в
Microsoft Excel