- •Математичні моделі в сапр: методичні вказівки до лабораторних робіт
- •Передмова до методичних вказівок
- •1.Лабораторна робота EvalExpr на тему "Знаходження похибок обчислень"
- •Теоретичні питання.
- •Приклад виконання роботи EvalExpr
- •Варіанти індивідуальних завдань до роботи EvalExpr
- •2.Лабораторна робота DetKram на тему “Розв’язання слр за формулами Крамера”
- •Теоретичні питання
- •Приклад виконання роботи DetKram
- •Варіанти індивідуальних завдань до роботи DetKram
- •Допоміжні матеріали до роботи DetKram Інструкція користувачу до додатку InteractiveDet
- •Приклад змісту вхідного файла Det.Inp:
- •Вихідні файли програми.
- •Натискання клавіші “Введення з файлу”.
- •Натискання клавіші “Зменшити у рядку”.
- •Натискання клавіші “Зменшити у стовпці”.
- •Натискання клавіші “Зменшити розміри”.
- •Завершення роботи додатку.
- •3.Лабораторна робота InvMatr на тему “Розв’язання слр за методом оберненої матриці”
- •Теоретичні питання
- •Приклад виконання роботи InvMatr
- •Варіанти індивідуальних завдань до роботи InvMatr
- •Допоміжні матеріали до роботи InvMatr Програма AdjMatr
- •Інструкція користувачу до додатку InteractiveAdjMatr
- •4.Лабораторна робота Gauss на тему „Наближене розв`язання слр методами виключення невідомих”
- •Теоретичні питання
- •Приклад виконання роботи Gauss
- •2В. Обчислення компонент вектора нев’язки.
- •3A. Знаходження наближеної оберненої матриці за методом Гауса.
- •3B. Знаходження компонент наближеного розв’язку .
- •3C. Знаходження компонент вектору нев’язки
- •5. Знаходження міри обумовленості головної матриці системи.
- •Варіанти індивідуальних завдань до роботи Gauss
- •Додаткові матеріали до роботи Gauss Програма NormVect знаходження різних норм вектора
- •Програма NormMatr знаходження різних норм матриці
- •Процедура g3v33l4
- •Модуль GaussRow
- •Iнструкцiя користувачу до програми InteractiveGauss
- •5.Лабораторна робота LinSysIt на тему “Ітераційні методи розв’язання слр”
- •Теоретичні питання
- •Приклад виконання роботи LinSysIt
- •Варіанти індивідуальних завдань до роботи LinSysIt
- •Допоміжні матеріали до роботи LinSysIt Програма NormForm
- •Програма SysIter
- •Програма Zeidel
- •6.Лабораторна робота TransEq на тему "Розв'язання трансцендентних рівнянь з одним невідомим"
- •Теоретичні питання
- •Приклад виконання роботи TransEq
- •Розрахункова таблиця знаходження кореня
- •Розрахункова таблиця знаходження кореня
- •Розрахункова таблиця знаходження кореня
- •Уточнення відокремленого кореня за методом простих ітерацій.
- •Розрахункова таблиця знаходження кореня рівняння за методом простих ітерацій
- •Розрахункова таблиця знаходження кореня рівняння за методом простих ітерацій
- •Варіанти індивідуальних завдань до роботи TransEq
- •Додаткові матеріали до роботи TransEq Інструкція користувачу для роботи з додатком TranscEq
- •Робота з додатком TranscEq
- •7.Лабораторна робота NonLinSys на тему „Розв’язання систем нелінійних рівнянь”
- •Теоретичні питання
- •Приклад виконання роботи NonLinSys
- •Розрахункова таблиця знаходження роз'вязку
- •Варіанти індивідуальних завдань до роботи NonLinSys
- •Допоміжні матеріали до роботи NonLinSys Інструкція користувачу для роботи з програмою Newton2
- •Текст програми Newton2
- •Текст модуля Funct
- •Текст модуля NewtonS
- •8.Лабораторна робота SelfEq на тему ”Знаходження характеристичного рівняння”
- •Теоретичні питання
- •Приклад виконання роботи SelfEq
- •Варіанти індивідуальних завдань до роботи SelfEq
- •9.Лабораторна робота SelfVect на тему “Знаходження власних векторів”
- •Теоретичні питання.
- •Приклад виконання роботи SelfVect
- •Варіанти індивідуальних завдань до роботи SelfVect
- •Література
3.Лабораторна робота InvMatr на тему “Розв’язання слр за методом оберненої матриці”
Завдання. Для заданої СЛР знайти головну матрицю системи, а також її обернену матрицю. Здійснити перевірку того, що знайдена матриця дійсно є оберненою і знайти розв’язок СЛР за методом оберненої матриці. Провести перевірку знайденого розв’язку.
Хід виконання. Для заданої системи лінійних рівнянь (СЛР):
Записати СЛР у матричному виді. Знайти її головну та розширену матриці.
Знайти визначник головної матриці СЛР.
Знайти доповняльні алгебраїчні мінори головної матриці СЛР і обернену матрицю до головної матриці СЛР. Здійснити перевірку того, що знайдена матриця дійсно є оберненою матрицею до головної матриці СЛР.
Знайти розв’язок системи за допомогою оберненої матриці. Здійснити перевірку знайденого розв’язку.
Знайти розв’язок СЛР з цілими коефіцієнтами за допомогою стандартної програми на ЕОМ або написати власну програму.
Оформити звіт та захистити лабораторну роботу.
Теоретичні питання
Квадратна матриця.
Визначник матриці, його властивості.
Ранг матриці, його властивості.
Обернена матриця, приєднана матриця та зв'язок між ними
Транспонована матриця, її властивості.
Доповняльні алгебраїчні мінори.
Алгоритм знаходження елементів оберненої матриці.
Приклад виконання роботи InvMatr
Лабораторна робота InvMatr на тему “Розв’язання СЛР за методом оберненої матриці”.
Варіант 33. Виконав студент групи БІТ-1-05, Іванов Олексій Вікторович. Дата виконання: 5.10.2005.
Завдaння. Розвязати систему лінійних рівнянь.
(1.3.3)
за допомогою оберненої матриці. Провести перевірку знайденого розвязку.
Виконання роботи. Систему (1.3.3) можна записати у матричному вигляді:
(1.3.4)
або , де - вектор-рядок невідомих , - вектор-рядок правої частини. - головна матриця системи. 1.Знаходження головного визначника. Використовуюч додаток Windows, InteractiveDet, одержуємо .
2.Знаходження доповняльних (алгебраїчних) мінорів. Будь-який доповняльний мінор одержуємо з головного визначника викреслюванням -го рядка і -го стовпця. Тому маємо:
, |
, |
, |
, |
, |
, |
, |
, |
, |
, |
, |
, |
, |
, |
, |
. |
3.Знаходження матрицю алгебраїчних доповнень. Алгебраїчне доповнення до елементу матриці системи (1.3.3) знаходиться за формулою ( ), то матриця алгебраїчних доповнень дорівнює
5. Знаходження оберненої матриці. Обернена матриця знаходиться за формулою
(1.3.5)
Тому
6. Перевірка знайденої оберненої матриці. Треба здійснити перевірку двох співвідношень і (за означенням оберненої матриці). Нехай , .
1) Перевірка першого співвідношення . Знайдемо матрицю Для її елементів маємо
,
,
,
,
,
,
,
,
,
,
,
,
,
,
.
Тому і перевірка першого з співвідношень означення оберненої матриці здійснена.
2) Перевірка другого співвідношення. Знайдемо матрицю . Для її елементів маємо
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
,
.
Тому і перевірка другого з співвідношень означення оберненої матриці здійснена.
Як наслідок, .
Тут
,
,
,
.
Отже, маємо вектор розв'язку
3.Перевірка знайденого вектора розв’язку. Перевірку знайденого вектора розв’язку здійснюємо підстановкою знайдених значень невідомих у ліві частини рівнянь системи і одержання значень правих частин. При маємо:
1) ;
2) ;
3)-5 ;
4)
5) .
Праві частини рівнянь системи дорівнюють лівим частинам. Тому вектор дійсно є розв’язком системи (1.3.3 ). Перевірка здійснена.
Висновки до лабораторної роботи. Точні розв’язки СЛР з цілими коефіцієнтами можуть бути числами не цілими, але, як це бачимо з формул для знаходження оберненої матриці ,є завжди числа дробові. Тому, при знаходженні розв’язку потрібно вести точні (без заокруглювань) підрахунки з дробами (у вигляді чисельник / знаменник), щоб одержати точні значення виразів, які ми підраховуємо.