- •Математичні моделі в сапр: методичні вказівки до лабораторних робіт
- •Передмова до методичних вказівок
- •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
- •Література
Теоретичні питання.
Власні вектори квадратної матриці.
Знаходження власних векторів різними методами (Данилевського,. Крилова, за означенням).
Система лінійних рівнянь для знаходження компонент власного вектора.
Використання методу обвідних мінорів для знаходження компонент власних векторів (знаходження лінійно незалежних рівнянь із відповідної системи лінійних рівнянь і тих компонент власного вектора, яки будуть вважатись незалежними і яки породжують власний підпростір).
Приклад виконання роботи SelfVect
Лабораторна робота SelfVal на тему “Знаходження характеристичного полінома, власних чисел і власних векторів квадратної матриці”.
Варіант 31. Виконав студент групи БІТ1-01, Довгий Олексій Вікторович. Дата виконання: 15.2.2005.
Завдання. Знайти власні значення і власні вектори квадратної матриці
за будь-яким методом (Данилевського, Крилова Левер'є,) або з використанням методу обвідних визначників.
Виконання роботи
Розвинення характеристичного визначника. Знайдемо характеристичний визначник матриці (1.8.33) за методом Левер'є з використанням формул Ньютона і властивостей слідів матриць.
З матриці (1.8.33) знаходимо , .
Знаходимо квадрат матриці :
.
За допомогою цього квадрату матриці знаходимо
, .
Знаходимо куб матриці А:
, ,
Тому
Отже, знайдені значення коефіцієнтів характеристичного полінома , тобто
(1.8.34)
Перевірка:
що відрізняється від лише знаком.
Знаходження всіх дійсних власних значень матриці А. Знайдемо всі власні значення матриці А, тобто корені характеристичного полінома .
Так як поліном є непарного степеня, то він завжди має хоча б один дійсний корінь. Спочатку шукаємо цілі розв'язки цього рівняння. За теорією всі цілі розв'язки алгебраїчного рівняння з цілими коефіцієнтами і старшим коефіцієнтом рівним одиниці є дільниками вільного коефіцієнта. Тому єдиними можливими цілими коренями рівняння можуть бути значення . Безпосередня перевірка дає:
,
,
, тобто- корінь рівняння.
Тому характеристичний визначник ділиться на . Здійснивши це ділення, отримаємо представлення: . Так як квадратне рівняння має корені і , то характеристичне рівняння має три кореня, розташовані за зростанням: .
Знаходження власного вектора матриці А, що відповідає власному числу . При отримуємо вироджену характеристичну систему лінійних рівнянь у вигляді системи рівнянь . Для розв'язання цієї системи треба за методом обвідних визначників знайти ранг системи і тій визначние, який визначає цій ранг. У випадку, що розглядається, ранг дорівнює 2 і визначається визначником (те, що випливає з означення власного вектора). Для розв'язання системи можна видалити рівняння, коефіцієнти якого не належать до визначника, з допомогою якого визначено ранг системи, і знаходити невідомі з системи через невідомі, коефіцієнти яких не містяться в цьому визначнику. Отже, можна відкинути останнє рівняння і розв'язувати систему відносно невідомих і , тобто . Тому, , . Якщо тут вибрати , то отримаємо власний вектор .
Перевірка: .
Отже, знайдені власне значення і власний вектор .
Висновки до лабораторної роботи. Нами знайдені дійсний корінь характеристичного рівняння та його власний вектор .