- •Поліном Лагранжа
- •Завдання для самостійної роботи
- •Тема: Апроксимація функцій Теоретичні відомості Емпіричні формули
- •Визначення параметрів емпіричної залежності
- •Метод найменших квадратів
- •Розв’язок
- •Завдання для самостійної роботи
- •Розв’язок
- •Метод Сімпсона
- •Розв’язок
- •Завдання для самостійної роботи
- •Лабораторна рОбота №3
- •Теоретичні відомості Чисельне розв’язання систем лінійних алгебраїчних рівнянь
- •Метод Гаусса
- •Метод Гаусса-Зейделя
- •Розв’язок
- •Завдання для самостійної роботи
- •Лабораторна рОбота №4 Тема: Нелінійні рівняння Теоретичні відомості
- •Метод хорд
- •Метод Ньютона (метод дотичних)
- •Розв’язок
- •Комбінований метод хорд і дотичних
- •Метод ітерацій або метод послідовних наближень
- •Розв’язок
- •Завдання для самостійної роботи
- •Лабораторна рОбота №5 Тема: Чисельне розв’язання звичайних диференційних рівнянь Теоретичні відомості
- •Метод Ейлера
- •Модифікації методу Ейлера
- •Метод Рунге-Кутта
- •Завдання для самостійної роботи
- •Список літератури
ЗМІСТ
Лабораторна РОБОТА №1 4
Тема: Інтерполяція та екстраполювання функцій 4
Теоретичні відомості 4
Лінійна і квадратична інтерполяції 4
Поліном Лагранжа 5
Завдання для самостійної роботи 5
Тема: Апроксимація функцій 7
Теоретичні відомості 7
Емпіричні формули 7
Визначення параметрів емпіричної залежності 7
Метод найменших квадратів 8
Завдання для самостійної роботи 11
Лабораторна рОбота №2 13
Тема: Чисельне інтегрування 13
Теоретичні відомості 13
Методи прямокутників і трапецій 13
Метод Сімпсона 16
Завдання для самостійної роботи 18
Лабораторна рОбота №3 21
Тема: Системи лінійних рівнянь 21
Теоретичні відомості 21
Чисельне розв’язання систем лінійних алгебраїчних рівнянь 21
Метод Гаусса 21
Метод Гаусса-Зейделя 24
Завдання для самостійної роботи 26
Лабораторна рОбота №4 31
Тема: Нелінійні рівняння 31
Теоретичні відомості 31
Метод хорд 32
Метод Ньютона (метод дотичних) 35
Комбінований метод хорд і дотичних 37
Метод ітерацій або метод послідовних наближень 38
Завдання для самостійної роботи 39
Лабораторна рОбота №5 41
Тема: Чисельне розв’язання звичайних диференційних рівнянь 41
Теоретичні відомості 41
Метод Ейлера 42
Модифікації методу Ейлера 44
Метод Рунге-Кутта 46
Завдання для самостійної роботи 48
СПИСОК ЛІТЕРАТУРИ 50
Лабораторна РОБОТА №1
Тема: Інтерполяція та екстраполювання функцій
Теоретичні відомості
Лінійна і квадратична інтерполяції
Простим типом локальної інтерполяції, яку часто використовують, є лінійна (або кусочно-лінійна) інтерполяція. Вона полягає у тому, що задані точки з'єднуються прямолінійними відрізками і функція наближається ламаною з вершинами в даних точках.
Рівняння кожного відрізку ламаної в загальному випадку різне. Оскільки є інтервалів , то для кожного з них у якості рівняння інтерполяційного багаточлена використовується рівняння прямої, яка проходить через дві точки. Зокрема, для i-го інтервалу можна написати рівняння прямої, що проходить через точки і наступному вигляді:
. |
Звідси маємо:
, , , . |
(1) |
Отже, при використанні лінійної інтерполяції спочатку необхідно визначити інтервал, в який потрапляє значення аргументу , а потім підставити його у формулу (1) і знайти наближене значення функції в цій точці.
Розглянемо тепер випадок квадратичної інтерполяції, коли інтерполяційна функція на відрізку приймається як поліном другого ступеню. Таку інтерполяцію називають також параболічною:
,
|
(2) |
Рівняння поліному другого ступеню (2) містить три невідомі коефіцієнти , для визначення яких необхідно мати три рівняння. Ними служать умови проходження параболи (2) через три точки , , . Ці умови можна записати у наступному виді:
|
(3) |
Поліном Лагранжа
Перейдемо до випадку глобальної інтерполяції, тобто до побудови інтерполяційного багаточлена єдиного для всього відрізку . Шукатимемо інтерполяційний багаточлен у вигляді лінійної комбінації багаточленів ступеня :
. |
(4) |
При цьому необхідно, щоб кожен багаточлен дорівнював нулю у всіх вузлах інтерполяції, за винятком одного ( -го), де він повинен дорівнювати одиниці. Легко перевірити, що цим умовам при = 0 відповідає багаточлен виду:
. |
(5) |
Дійсно . При чисельник виразу (5) дорівнює нулю. По аналогії з (5) отримаємо:
, , ……………………………………………… , , .......................................................................... .
|
(6) |
Підставляючи в (4) вирази (5), (6), знаходимо:
. |
(7) |
Формула (7) визначає інтерполяційний багаточлен Лагранжа.
Завдання для самостійної роботи
Знайти наближене значення функції при заданому значенні аргументу за допомогою інтерполяційного багаточлена Лагранжа, якщо функція задана:
1. у не рівновіддалених вузлах;
2. у рівновіддалених вузлах таблиці.
№ варіанту |
№ таблиці |
|
|
№ варіанту |
№ таблиці |
|
|
1 |
1 |
0,702 |
1,3832 |
16 |
4 |
0,665 |
0,1944 |
2 |
2 |
0,102 |
0,1264 |
17 |
5 |
0,774 |
0,2232 |
3 |
3 |
0,526 |
0,1521 |
18 |
6 |
0,332 |
1,4396 |
4 |
4 |
0,616 |
0,1838 |
19 |
1 |
0,736 |
1,3934 |
5 |
5 |
0,896 |
0,2121 |
20 |
2 |
0,203 |
0,1334 |
6 |
6 |
0,314 |
1,4179 |
21 |
3 |
0,552 |
0,1543 |
7 |
1 |
0,512 |
1,3926 |
22 |
4 |
0,537 |
0,1676 |
8 |
2 |
0,114 |
0,1315 |
23 |
5 |
0,955 |
0,2263 |
9 |
3 |
0,453 |
0,1611 |
24 |
6 |
0,275 |
1,4236 |
10 |
4 |
0,478 |
0,1875 |
25 |
1 |
0,608 |
0,3866 |
11 |
5 |
0,812 |
0,2165 |
26 |
2 |
0,154 |
0,1285 |
12 |
6 |
0,235 |
1,4258 |
27 |
3 |
0,436 |
0,1625 |
13 |
1 |
0,645 |
1,3862 |
28 |
4 |
0,673 |
0,2038 |
14 |
2 |
0,125 |
0,1232 |
29 |
5 |
0,715 |
0,2244 |
15 |
3 |
0,482 |
0,1662 |
30 |
6 |
0,186 |
1,4315 |
Таблиці до завдання 1
Таблиця 1 |
Таблиця 2 |
Таблиця 3 |
||||||||||||||||||||||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||||||||||||||
Таблиця 4 |
Таблиця 5 |
Таблиця 6 |
||||||||||||||||||||||||||||||||||||||||||
|
|
|
Таблиці до завдання 2
Таблиця 1 |
Таблиця 2 |
Таблиця 3 |
||||||||||||||||||||||||||||||||||||||||||
|
|
|
||||||||||||||||||||||||||||||||||||||||||
Таблиця 4 |
Таблиця 5 |
Таблиця 6 |
||||||||||||||||||||||||||||||||||||||||||
|
|
|