- •Математичні моделі в сапр: методичні вказівки до лабораторних робіт
- •Передмова до методичних вказівок
- •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
- •Література
7.Лабораторна робота NonLinSys на тему „Розв’язання систем нелінійних рівнянь”
Завдання. Графічно відокремити розв’язки системи двох нелінійних рівнянь з двома невідомими і наближено знайти ці розв’язки з точністю ε = 0.001 за методом Ньютона.
Хід виконання. Для заданої системи здійснити наступні дії:
А) Відокремити графічно розв’язки цієї системи, знаходячи точки перетину відношень і на площині, кожне з яких локально можна записати як деякі функції від і від . Знайти початкові (стартові) наближення { , }.
Б) Уточнити один відокремлений розв’язок за методом Ньютона. Для цього:
Знайти якобіан .
Скласти ітераційну схему двовимірного метода Ньютона, тобто
а) Вибрати стартове наближення { , };
б) Записати ітераційний крок: по-перше, розв’язати систему лінійних рівнянь (наприклад, за методом Гауса) відносно невідомих . Ця система лінійних рівнянь має матричний вигляд:
,
де - головна матриця системи, і, по-друге, знайти нове наближення за формулами .
в) Вибрати момент зупинки (за числом ітерацій N або за величиною поточного кроку , де - вибрана раніше векторна норма).
Провести обчислення трьох ітерацій без використання програми на ЕОМ. Обчислення звести в таблицю Таблиця 1 (з наступними графами: Графа 1 – Номера ітерацій k; Графи 2,3 – наближення розв’язку на попередній ітерації; Графи 4,5 – прирости на поточному кроці; Графи 6,7 – нове наближення розв’язку , що одержано на k-й ітерації).
Написати на мові програмування Turbo Pascal власну програму Newton2d знаходження розв’язків системи нелінійних рівнянь за методом Ньютона або скористатись тією, що наведена у допоміжних матеріалах до лабораторної роботи. При цьому праві частини , , , , , потрібно реалізувати у окремому модулі (як тексти відповідних функцій) і процес ітерацій виводити в окремий текстовий файл.
Провести уточнення розв’язку, що відокремлений, за допомогою програми Newton2d для ітерацій.. Порівняти результати таблиці 1 з підрахунками на ЕОМ.
Оформити звіт і захистити лабораторну роботу.
Теоретичні питання
Простір Rn. Норма в Rn. Приклад норм (m-норма, l-норма, k-норма).
Перетворення або . Якобіан перетворення .
Система нелінійних рівнянь у -вимірному просторі. Багатовимірний ітераційний метод Ньютона.
Умови збіжності метода Ньютона.
Системи лінійних рівнянь нормального вигляду у -вимірному просторі, де . Нерухомі точки (фікспункти) відображення .
Нерівність Ліпшиця для відображення . Стискуючі перетворення. Коефіцієнт стиску.
Головна теорема про розв’язок системи лінійних рівнянь нормального виду , де - стискуюче відображення.
Метод простих ітерацій у -вимірному просторі. Відповідна ітераційна схема.
Приклад виконання роботи NonLinSys
Лабораторна робота TrancEq на тему “Розв’язання трансцендентних рівнянь з двома невідомими”.
Варіант 31. Виконав студент групи БІТ-04, Павлік Олексій Васильович. Дата виконання: 05.11.2003
Завдання. Графічно відокремити розв’язки системи нелінійних рівнянь з двома невідомими і уточнити будь-який один з них за методом Ньютона.
Виконання роботи.
Відокремлення коренів системи рівнянь. Перше рівняння записуємо у вигляді , друге – у вигляді . Будуємо (креслення відповідних графіків залишаємо студентові) на одному малюнку ці два відношення як два графіка: перший - як функцію від , другий - як функцію від . Графічно бачимо, що один з наближених розв’язків розташовано біля точки . Тому беремо цю точку за початкове наближення:
Зауваження. На жаль, якщо ви помилились при побудові графіків, то знайдене стартове наближення не буде отримане. Треба бути уважним при побудові графіків тому, що єдиною можливістю перевірки вірності побудови є збіжність ітераційного процесу.
1. Уточнення відокремленого кореня за методом Ньютона без використання ЕОМ. Отже, графічно нами були знайдені стартові значення: для системи : Переходимо до виконання ітераційних кроків.
Крок 1. Підраховуємо значення функцій лівих частин і їх похідних в стартовій (нульовий) точці:
.
Тому маємо якобіан і для приростів і систему лінійних рівнянь .
За формулами Крамера: ,
,
.
Тому знаходимо наступне наближення:
.
.
Так як умови , не виконані, то обчислення продовжуємо.
Крок 2. Підраховуємо значення функцій лівих частин і їх похідних в наступній (першій) точці:
.
Тому маємо якобіан і для приростів і - систему лінійних рівнянь .
За формулами Крамера: ,
Тому ,
.
Так як умови , не виконані, то продовжуємо обчислення.
Крок 3. Підраховуємо значення функцій лівих частин і їх похідних в наступній (другій) точці:
Тому маємо якобіан і для приростів і систему лінійних рівнянь .
За формулами Крамера : ,
,
.
Тому ,
Так як умови , виконані , то обчислення завершуємо .
Зведемо пiдрахунки у таблицю .
Таблиця 1.