Методичні вказівки та завдання

МЕТОДИЧНІ ВКАЗІВКИ ТА ЗАВДАННЯ

до лабораторних та самостійних робіт

з курсів „Основи інформатики та обчислювальної техніки”

та „Алгоритмічні мови та програмне забезпечення”

за темою „Чисельні методи в інженерних розрахунках”

 

для студентів технічних спеціальностей денної форми навчання

 

 

Методичні вказівки та завдання до лабораторних та самостійних робіт з курсів „Основи інформатики та обчислювальної техніки” та „Алгоритмічні мови та програмне забезпечення” за темою „Чисельні методи в інженерних розрахунках” для студентів технічних спеціальностей денної форми навчання /Укл.:В.Г.Вишневська, О.В.Корнєєва, Г.В.Романіченко, Л.І.Лозовська, К.В.Пугіна. – Запоріжжя: ЗНТУ,2005. – 65с.

 

Методичні вказівки містять індивідуальні завдання та теоретичні відомості до лабораторних робіт з курсів „Основи інформатики та обчислювальної техніки” та „Алгоритмічні мови та програмне забезпечення”, приклади їх виконання з використанням середовища Excel та програмування мовою VBA для студентів технічних спеціальностей денної форми навчання.

 

Укладачі: В.Г. Вишневська, доцент,

О.В. Корнєєва, асистент,

Г.В. Романіченко, асистент,

Л.І. Лозовська, доцент,

К.В. Пугіна, аспірант.

 

 

Рецензенти: Н.І. Біла, доцент,

О.І. Денисенко, доцент.

 

Експерт: А.О. Шумілов, к.т.н., доцент.

 

Відповідальний

за випуск Г.В. Корніч, професор.

 

Затверджено

на засіданні кафедри

обчислювальної математики

протокол № 8 від 29.04.05 р.

ЗМІСТ

 

1 Лабораторна робота №1: Методи розв’язання алгебраїчних і трансцендентних рівнянь.......................................................................... 5

1.1 Відділення числового проміжку, у якому міститься один корінь рівняння...................................................................................................... 5

1.1.1 Відділення кореня графічно (перший спосіб)..................... 5

1.1.2Другий спосіб відділення кореня......................................... 6

1.2 Уточнення значення кореня рівняння ƒ(x) = 0........................... 6

1.2.1 Метод половинного ділення (метод бісекцій).................... 6

1.2.2 Метод хорд (метод пропорційних чисел)............................ 6

1.2.3 Метод Ньютона (метод дотичних)....................................... 7

1.2.4 Використання пакету аналізу „что - если” Excel................ 8

1.3 Індивідуальні завдання.................................................................. 8

1.4 Приклади виконання лабораторної роботи............................... 12

2 Лабораторна робота №2: Методи розв’язання систем нелінійних рівнянь...................................................................................................... 21

2.1 Відділення першого наближення до рішення............................ 21

2.2 Методи уточнення рішення системи нелінійних рівнянь....... 21

2.2.1 Метод ітерацій...................................................................... 21

2.2.2 Використання пакету аналізу Excel................................... 22

2.3 Індивідуальні завдання................................................................ 23

2.4 Приклади виконання лабораторної роботи............................... 27

3 Лабораторна робота №3: Обчислення інтегралів............................ 33

3.1 Теоретичні відомості................................................................... 33

3.2 Індивідуальні завдання................................................................ 34

3.3 Приклади виконання лабораторної роботи............................... 38

4 Лабораторна робота №4: Наближення (інтерполяція) функцій...... 42

4.1 Теоретичні відомості................................................................... 42

4.2 Індивідуальні завдання............................................................... 43

4.3 Приклади виконання лабораторної роботи............................... 48

5 Лабораторна робота №5: Апроксимація даних (емпіричні формули).................................................................................................. 51

5.1 Теоретичні відомості.................................................................. 51

5.1.1 Визначення параметрів емпіричних формул по способу найменших квадратів у випадку лінійної залежності......................... 51

5.1.2 Визначення параметрів емпіричних формул по способу найменших квадратів у випадку нелінійної залежності...................... 53

5.2 Індивідуальні завдання................................................................ 53

5.3 Приклад виконання лабораторної роботи.................................. 53

6 Лабораторна робота №6: Визначення параметрів емпіричних формул по способу найменших квадратів у випадку нелінійної залежності................................................................................................ 57

6.1 Теоретичні відомості................................................................... 57

6.2 Індивідуальні завдання................................................................ 58

6.3 Приклад виконання лабораторної роботи..................................61

6.3.1 Метод Ейлера........................................................................61

6.3.2 Метод Рунге-Кутта...............................................................62

7 Література............................................................................................ 64

8 Вимоги до оформлення лабораторної роботи................................... 64

8.1 Додаток А...................................................................................... 65

ЛАБОРАТОРНА РОБОТА № 1

ТЕМА: Методи розв’язання алгебраїчних і трансцендентних рівнянь

Розв’язання рівнянь – алгебраїчних і трансцендентних – являє собою одну з істотних задач прикладного аналізу, потреба в якій виникає в найрізноманітніших розділах фізики, техніки і природознавства.

Задача визначення кореня рівняння з одним невідомим

ƒ(x) = 0, (1.1)

де ƒ(x) – безперервна функція, складається з двох етапів:

1) відділення кореня, тобто визначення числового проміжку, у якому міститься один корінь рівняння;

2) уточнення значення кореня шляхом побудови послідовності

x_к = φ (x_к-1), к = 1, 2, 3,...

на основі відповідного методу.

Для уточнення значення кореня існують різні ітераційні методи.

1.1 Відділення числового проміжку, у якому міститься один корінь рівняння

1.1.1 Відділення кореня графічно (перший спосіб)

Якщо рівняння (1.1) зручно представити у вигляді

g (х) – h (х), (1.2)

то абсцису х₀ точки перетинання графіків у = g(х) і у = h (х) можна знайти по кресленню.

Величину х₀ визначити з достатньою точністю графічно не можливо. Тому варто вибрати такий числовий проміжок [a; b] для якого свідомо виконується нерівність a ≤ х₀ ≤ b.

Різні знаки функції при х =а і х = b

ƒ(а) * ƒ(b) ≤ 0 (1.3)

свідчать про наявність кореня в проміжку [a; b].

 

 

Другий спосіб відділення кореня

Цей спосіб містить звичайне табулювання функції у = ƒ (х) на інтервалі існування функції, при цьому ступінь зміни аргументу підбирається значимим. І знову, різні знаки функції при х =а і х = b, тобто ƒ(а) * ƒ(b) ≤ 0 свідчать про наявність кореня в проміжку [a; b].

1.2 Уточнення значення кореня рівняння ƒ(x) = 0

 

Індивідуальні завдання

Для кожної функції визначити числовий проміжок у якому міститься один корінь рівняння, уточнити значення кореня вказаними вище методами.

1.3.1 1.3.2

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) , . 4) .

1.3.3 1.3.4

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) .

 

1.3.5 1.3.6

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) , .

1.3.7 1.3.8

1) ; 1) ;

2) ; 2)

3) ; 3) ;

4) . 4) .

1.3.9 1.3.10

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) .

1.3.11 1.3.12

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) , . 4) .

1.3.13 1.3.14

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) .

1.3.15 1.3.16

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) ,.

1.3.17 1.3.18

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) , . 4) .

1.3.19 1.3.20

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) .

1.3.21 1.3.22

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) , . 4) .

 

1.3.23 1.3.24

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) .

1.3.25 1.3.26

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) , .

1.3.27 1.3.28

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) .

1.3.29 1.3.30

1) ; 1) ;

2) ; 2) ;

3) ; 3) ;

4) . 4) ,

ЛАБОРАТОРНА РОБОТА № 2

Метод ітерацій

Зведемо систему (2.1) до вигляду

, (2.3)

по встановленому початковому наближенню до рішення х₀, у₀, уточнення рішення відбувається за формулами:

 

(2.4)

Умови збіжності методу припускають, що

(2.5)

Рішення системи (2.1) вважається знайденим з точністю ε, якщо виконуються умови:

 

 

Індивідуальні завдання

Для кожної системи нелінійних рівнянь визначити початкове наближення до рішення та уточнити це рішення з точністю ε = 0,001.

 

 

2.3.1

1) 2)

 

2.3.2

1) 2)

 

2.3.3

1) 2)

2.3.4

1) 2)

 

2.3.5

1) 2)

 

2.3.6

1) 2)

2.3.7

1) 2)

 

2.3.8

1) 2)

2.3.9

1) 2)

2.3.10

1) 2)

2.3.11

1) 2)

2.3.12

1) 2)

2.3.13

1) 2)

2.3.14

1) 2)

2.3.15

1) 2)

2.3.16

1) 2)

2.3.17

1) 2)

2.3.18

1) 2)

2.3.19

1) 2)

2.3.20

1) 2)

2.3.21

1) 2)

2.3.22

1) 2)

2.3.23

1) 2)

2.3.24

1) 2)

2.3.25

1) 2)

2.3.26

1) 2)

2.3.27

1) 2)

2.3.28

1) 2)

2.3.29

1) 2)

2.3.30

1) 2)

ЛАБОРАТОРНА РОБОТА № 3

ТЕМА: Обчислення інтегралів

3.1 Теоретичні відомості

 

Точне обчислення визначеного інтеграла за формулою Ньютона-Лейбниця не завжди представляється можливим чи доцільним. У цих випадках, а також у випадках, коли підінтегральна функція задана табличним способом, доцільно проводити обчислювання приблизно. До найпростіших методів відносяться формула прямокутників, формула парабол (Сімпсона).

При обчисленнях за формулами прямокутників результати будуть тим точніші, чим більше число відрізків розбивки проміжку інтегрування. Найпростіша - формула лівих прямокутників:

, (3.1)

де

Для забезпечення наданої степені точності ε значення n визначається за умовою:

(3.2)

За формулою парабол маємо:

(3.3)

де

Згідно екстраполяції по Річардсону більш точне значення інтеграла можна визначити за формулою:

(3.4)

де , - два наближених значення ,знайдених по формулі параболи для n₁ та n₂ (n₂ > n₁).

 

Індивідуальні завдання

Для кожного варіанту:

1) обчислити інтеграл по формулі прямокутників з точністю ε = 0,01;

2) по формулі парабол обчислити інтеграл, уточнювання провести згідно формулі Річардсона.

 

3.2.1

1) ; 2) .

 

3.2.2

1) ; 2) .

 

3.2.3

1) ; 2) .

3.2.4

1) ; 2) .

 

3.2.5

1) ; 2) .

3.2.6

1) ; 2) .

3.2.7

1) ; 2) .

3.2.8

1) ; 2) .

3.2.9

1) ; 2) .

3.2.10

1) ; 2) .

3.2.11

1) ; 2) .

3.2.12

1) ; 2) .

3.2.13

1) ; 2) .

3.2.14

1) ; 2) .

3.2.15

1) ; 2) .

3.2.16

1) ; 2) .

3.2.17

1) ; 2) .

3.2.18

1) ; 2) .

3.2.19

1) ; 2) .

3.2.20

1) 2) .

3.2.21

1) ; 2) .

3.2.22

1) ; 2) .

3.2.23

1) ; 2) .

3.2.24

1) ; 2) .

3.2.25

1) ; 2) .

3.2.26

1) ; 2) .

3.2.27

1) ; 2) .

3.2.28

1) ; 2) .

3.2.29

1) ; 2) .

3.2.30

1) ; 2) .

 

3.3 Приклади виконання лабораторної роботи

 

Приклад 3.3.1 Обчислити інтеграл за формулою прямокутників з точністю ε=0,01.

 

Вирішимо задачу, використовуючи середовище Excel. Треба обчислити функцію на проміжку [0,7;1,3].

Кількість відрізків розбивки проміжку інтегрування визначимо з умови (3.2):

,

 

,

 

М=2, тоді ; .

Розмістимо початкові значення та таблицю значень x, у у такому вигляді, як показано на рис. 3.1.

 

Рисунок 3.1 – Таблиця початкових значень та значень х, у

 

Згідно протоколу рішення, який показано на рис. 3.2, отримаємо відповідь: 0,406865.

 

 

Рисунок 3.2 – Протокол рішення інтегралу за формулою прямокутників

 

Приклад 3.3.2 Обчислити інтеграл за формулою

парабол, уточнювання провести згідно формули Річардсона.

Для вирішення задачі з’ясуємо початкові дані: a = 1,2; b = 1,6, , оберемо n₁ = 10, n₂ = 15, які розмістимо на листі Excel з ім’ям „Сімпсона”. На цьому ж листі розмістимо кнопку Command Button1, як показано на рис. 3.3 та створимо процедуру рішення.

 

Рисунок 3.3 – Зразок листа Excel з ім’ям „Сімпсона”

 

Текст процедури:

Зразок протоколу рішення наведено на рис. 3.4.

 

 

Рисунок 3.4 – Протокол рішення інтегралу за формулою парабол

 

Згідно протоколу рішення маємо:

= 0б08389.

ЛАБОРАТОРНА РОБОТА №4

Теоретичні відомості

Задача наближення функції виникає, коли для функції, даної при дискретних значеннях аргументу у вигляді таблиці (ці значення називаються вузлами інтерполяції) необхідно знайти значення функції в проміжних крапках. Накладаючи вимогу, щоб наближена функція у вузлах співпадала з табличними значеннями (рис. 4.1), одержуємо задачу інтерполяції.

 

 

Рисунок 4.1 – Графік наближеної функції

 

Нехай в результаті спостережень за ходом деякого процесу побудована таблиця:

 

x	x0	x1	x2	…	xn
f(x)	f(x0)	f(x1)	f(x2)	…	f(xn)

 

Тобто, функція f(x) задана таблицею значень для кінцевої безлічі значень х.

Якщо необхідно знайти значення f(x) для проміжного значення аргументу, то будують функцію φ(x), просту для обчислень і таку, що для заданих x₀, x₁, x₂,..., x_n приймає значення f(x₀), f(x₁), f(x₂),..., f(x_n).

В інших точках відрізка [x_0, x_n] вважаємо, що φ(x) приблизно визначає функцію f(x) з тим чи іншим ступенем точності.

Найчастіше, функцію φ(x) представляють у вигляді алгебраїчного багаточлена деякого ступеня.

Найпростіша інтерполяція – це лінійна, тобто, коли невідому аналітичну залежність f(x) замінюють відрізками прямих, які проходять через відповідні вузли інтерполяції. В цьому випадку потрібно визначити якому відрізку належить надане х^* і за формулою лінійної інтерполяції знаходять f(x^*). Якщо x_i <= x^* <= x_i+1, то відповідна пряма проходить через вузли (x_i, f(x_і)), (x_i+1, f(x_і+1)):

(4.1)

 

Точність підрахунків в цьому випадку незначна, тому що враховується вплив тільки 2-ох вузлів інтерполяції. Частіше будують багаточлен P_n(x) ступеня n, що в (n+1) даних точках x₀, x₁, x₂,..., x_n. приймає дані значення y₀ = f(x₀), y₁ = f(x₁), …, y_n = f(x_n), тобто

f(x_і) = P_n(x_і), (і = 0, 1, 2,..., n).

Відзначимо, що двох різних інтерполяційних багаточленів одного і того же ступеня n існувати не може. Цим умовам задовольняє інтерполяційний багаточлен Лагранжа:

(4.2)

 

Тоді

 

Індивідуальні завдання

 

Для кожного варіанту обчислити наближене значення функції, яка задана таблицею, для наданого х^*, використовуючи лінійну інтерполяцію по Лагранжу.

 

4.2.1 х^* = 1,50

 

x	1,21	1,29	1,45	1,61	1,94	2,22
y	3,54	4,11	4,78	4,33	4,01	3,66

4.2.2 х^* = 0,45

 

x	0,24	0,39	0,65	0,88	0,94	1,12
y	2,18	1,98	1,73	1,13	1,01	1,34

 

4.2.3 х^* = 0,66

 

x	0,28	0,42	0,56	0,63	0,84	0,98
y	5,26	6,21	6,66	6,93	7,25	5,99

 

4.2.4 х^* = 3,2

 

x	2,10	2,19	2,35	2,51	2,82	3,33
y	8,54	8,11	7,78	7,03	6,77	7,86

 

4.2.5 х^* = 3,50

 

x	3,12	3,28	3,54	3,72	3,88	4,12
y	7,04	6,50	6,11	5,23	4,81	6,12

 

4.2.6 х^* = 0,85

 

x	0,25	0,54	0,65	0,80	0,94	1,12
y	0,54	0,91	1,78	2,33	2,81	3,66

 

4.2.7 х^* = 1,74

 

x	1,26	1,39	1,55	1,71	2,08	2,20
y	1,58	4,06	3,98	3,46	3,12	2,42

 

4.2.8 х^* = 4,51

 

x	4,02	4,38	4,44	4,58	4,74	4,98
y	10,54	4,82	4,93	5,14	5,92	7,12

 

4.2.9 х^* = 1,90

 

x	1,22	1,69	1,85	2,22	2,46	2,72
y	0,52	0,98	1,78	2,33	2,01	1,66

 

4.2.10 х^* = 3,82

 

x	3,41	3,59	3,75	3,89	4,00	4,18
y	4,50	4,19	3,81	4,33	3,66	2,84

 

4.2.11 х^* = 2,22

 

x	2,12	2,42	2,65	2,71	2,99	3,12
y	1,50	1,92	2,78	2,03	3,16	2,88

 

4.2.12 х^* = 5,32

 

x	5,13	5,29	5,45	5,61	5,84	6,22
y	1,54	1,11	1,78	1,33	1,01	2,66

 

4.2.13 х^* = 1,35

 

x	1,21	1,29	1,45	1,61	1,92	1,22
y	3,54	4,11	4,78	4,33	4,01	3,66

 

4.2.14 х^* = 0,80

 

x	0,22	0,42	0,63	0,78	0,82	0,99
y	6,52	6,88	7,55	6,15	7,92	8,02

 

4.2.15 х^* = 2,85

 

x	2,22	2,39	2,55	2,72	2,84	3,02
y	2,54	2,11	2,78	3,33	2,01	4,66

 

4.2.16 х^* = 0,60

 

x	0,28	0,49	0,65	0,71	0,96	1,22
y	3,54	4,11	4,78	4,33	4,01	3,66

 

4.2.17 х^* = 1,52

 

x	1,12	1,29	1,45	1,61	1,94	2,22
y	3,54	4,11	3,08	4,33	4,01	6,06

 

4.2.18 х^* = 3,50

 

x	3,12	3,29	3,45	3,66	3,84	4,02
y	3,54	2,11	1,78	1,33	2,01	0,54

 

4.2.19 х^* = 1,44

 

x	1,18	1,39	1,55	1,71	2,94	3,32
y	6,54	4,11	3,78	6,33	3,51	3,04

 

4.2.20 х^* = 0,56

 

x	0,12	0,31	0,45	0,61	0,98	1,32
y	2,54	3,22	3,58	3,13	4,01	4,26

 

4.2.21 х^* = 1,82

 

x	1,22	1,29	1,55	1,71	2,04	2,32
y	3,04	2,77	2,44	2,08	2,84	1,86

 

4.2.22 х^* = 2,50

 

x	2,32	2,49	2,65	2,81	3,04	3,18
y	3,54	4,11	4,78	4,33	4,01	5,66

 

4.2.23 х^* = 3,35

 

x	3,12	3,29	3,45	3,61	3,94	3,22
y	4,54	5,11	5,78	5,33	6,02	6,66

 

4.2.24 х^* = 1,50

 

x	1,32	1,29	1,44	1,66	1,98	2,26
y	3,54	3,11	2,78	2,33	4,01	2,16

 

4.2.25 х^* = 1,56

 

x	1,28	1,39	1,49	1,61	1,94	2,22
y	0,56	0,88	1,78	2,33	2,01	2,66

 

4.2.26 х^* = 1,08

 

x	0,07	1,29	1,45	1,61	1,94	2,22
y	3,54	4,11	4,78	4,33	4,22	3,66

 

4.2.27 х^* = 3,33

 

x	3,12	3,39	3,48	3,61	3,94	3,22
y	3,54	3,11	2,78	2,33	4,01	2,16

 

4.2.28 х^* = 1,44

 

x	0,88	1,26	1,35	1,58	1,74	1,98
y	2,46	2,71	2,98	3,33	3,00	3,37

 

4.2.29 х^* = 1,50

 

x	1,42	1,59	1,65	1,71	1,94	2,22
y	5,12	4,