- •1. Чисельні методи та обчислювальна математика
- •1.1. Роль математичного моделювання у розв’язанні задач навколишнього світу
- •1.2. Структура похибки розв’язку задачі
- •1.3. Поняття стійкості та коректності
- •1.4. Основні етапи розв’язання прикладних математичних задач
- •1.5. Контрольні питання
- •2. Елементарна теорія похибок
- •2.1. Джерела і класифікація похибок
- •2.3. Вірні значущі цифри
- •2.2. Абсолютна та відносна похибки
- •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.4. Точні методи розв’язання систем лінійних алгебраїчних рівнянь
4.4.1. Метод Крамера
Нехай дана система лінійних рівнянь, у якій число рівнянь дорівнює числу невідомих:
,
(4.6)
де
–
відповідно матриця системи, стовпець
вільних членів і стовпець невідомих.
Припустимо,
що визначник системи
Якщо тепер замінити послідовно у визначнику d стовпці коефіцієнтів при невідомих xj (j = 1, 2, ..., n) стовпцем вільних членів bi то вийдуть відповідно визначники
Теорема Крамера. Система п лінійних рівнянь з п невідомими, визначник якої відмінний від нуля, завжди сумісна і має єдиний розв’язок, що обчислюється за формулами:
(4.7)
Формули (4.7) називаються формулами Крамера.
Приклад.
Розв’язати за формулами Крамера систему
лінійних рівнянь
1)
Обчислюємо визначник системи:
2) Обчислюємо визначники, складені з коефіцієнтів при невідомих xl, x2, x3:
3) Використовуючи формули Крамера (4.7), знаходимо розв’язок системи: x1 = d1/d = 6/6 = 1; x2 = d2/d = 12/6 = 2; x3 =–12/6 =–2.
Зазначимо, що розв’язання СЛАР за формулами Крамера дуже громіздке, тому на практиці системи, зазвичай, розв’язують іншими методами, які будуть розглянуті нижче.
4.4.2. Метод Гаусса
Метод послідовного виключення невідомих або метод Гаусса розв’язання системи (4.1.) полягає в такому [1 – 4].
Припускаючи,
що
(це завжди можна зробити за рахунок
нумерації рівнянь), множимо перше
рівняння на
і додаємо до другого. Одержуємо
рівняння, у якому коефіцієнт при x1
дорівнює нулю. Помноживши перше рівняння
на
і додавши до третього, одержуємо
рівняння,
що
також не містить
x1.
Аналогічним
шляхом можна перетворити всі інші
рівняння, у результаті чого виходить
еквівалентна вихідній
система рівнянь
.
(4.8)
де
– деякі нові коефіцієнти.
Припускаючи,
що
,
і залишаючи незмінними перші
два
рівняння системи (4.8), перетворимо її
так, щоб у кожному
з інших рівнянь коефіцієнт при x2
дорівнював нулю. Після кінцевого числа
перетворень систему (4.8) можна привести
до одного з таких виглядів:
,
де
(4.9)
,
де
(4.10)
,
де
(4.11)
Система (4.9) має єдиний розв’язок; значення хп знаходять з останнього рівняння цієї системи, значення хп–1 – з передостаннього, значення х1 – з першого.
Система
(4.10) має нескінченну
множину рішень.
З останнього рівняння можна виразити
одне з невідомих (наприклад, хk)
через
інші n
– k
невідомих
,
що
входять
у це рівняння. З передостаннього рівняння
виразиться хk–1,
з
першого рівняння визначиться
х1
через ті ж невідомі. В отриманих формулах
невідомі
можуть
приймати будь–які
значення.
Система (4.11) несумісна, оскільки ніякі значення невідомих не можуть задовольнити її останньому рівнянню.
Зауваження 1. Розв’язуючи СЛАР методом Гаусса, перетворення роблять не над рівняннями, а над матрицями, складеними з коефіцієнтів при невідомих і вільних членів. Так, системі (4.1) відповідає матриця
,
(4.12)
системі (4.9) – «трикутна» матриця
,
(4.13)
системі (4.10) – «трапецієподібна» матриця
,
(4.14)
системі (4.11) – матриця
,
(4.15)
яка містить рядок, в якому дорівнюють нулю всі елементи, за винятком одного, що відповідає вільному члену. Вертикальною рискою відділений стовпець, складений з вільних членів.
Приклад
1.
Розв’язати СЛАР
.
Складемо
матрицю з коефіцієнтів при невідомих
і вільних членів:
.
Помноживши
перший рядок на (–3) і додавши до другого
рядка, одержуємо
матрицю
.
Помноживши
перший рядок на (–5) і додавши до третього
рядка, одержуємо
матрицю
.
Помноживши
другий рядок матриці А2
на
i додавши до третього рядка, одержуємо
матрицю
,
якій
відповідає
система рівнянь
.
Розв’язавши систему, знаходимо: х3=–3, x2=1, x1=2. Отже, вихідна система має той же розв’язок: x1=2, x2=1, х3= –3.
