
- •Основи чисельних методів математики (з використанням Excel) Передмова
- •Розділ 1. Методи обчислень: предмет, основні поняття та застосування
- •§ 1. Предмет і застосування
- •§ 2. Основні поняття
- •1. Похибки наближень.
- •2. Граничні похибки. Похибки функції.
- •3. Похибки розв'язку.
- •4. Стійкість і коректність.
- •Питання, тести
- •Розділ 2. Інтерполяція функцій
- •§1. Задача інтерполювання
- •§2. Інтерполяційна формула Лагранжа
- •§3. Поділені різниці. Формула Ньютона з поділеними різницями
- •§4. Інтерполяційна формула за допомогою Excel
- •§5. Інтерполювання за схемою Ейткіна
- •§6. Скінчені різниці. Інтерполяційні формули Ньютона для рівновіддалених вузлів
- •§7. Інтерполювання із скінченими різницями за допомогою Excel
- •§8. Інші методи інтерполювання
- •Питання, тести
- •Завдання
- •Розділ 3. Чисельне диференціювання та інтегрування.
- •§ 1. Однобічні формули чисельного диференціювання
- •§ 2. Оцінки похибки чисельного диференціювання
- •§ 3. Чисельне інтегрування. Квадратурні формули
- •§ 4. Квадратурні формули Ньютона – Котеса
- •§ 5. Узагальнені квадратурні формули.
- •§ 6. Метод подвійного перерахунку.
- •1. R2n ( f ) ≈ (правило Рунге) (14)
- •§ 7. Метод кратного перерахунку за допомогою Excel
- •Питання, тести
- •Завдання
- •Розділ 4. Чисельні методи розв‘язування рівнянь з однією змінною
- •§ 1. Відокремлення коренів
- •§ 2. Метод дихотомії (поділу відрізка пополам)
- •§ 3. Ітераційні методи та оператор стиску.
- •§ 4. Похибки ітераційного процесу
- •§ 5. Реалізація методу простої ітерації за допомогою електронних таблиць
- •§ 6. Метод Ньютона. Порядок збіжності ітераційного процесу.
- •§ 7. Метод лінійного інтерполювання.
- •§ 8. Інші приклади ітераційних методів.
- •Питання, тести
- •Завдання
- •Розділ 5. Методи розв’язування систем лінійних рівнянь
- •§ 1. Метод Гаусса
- •Метод Гаусса в матричній формі
- •Елементарні операції над матрицею:
- •§ 2. Метод Гаусса за допомогою Excel
- •§ 3. Матричні операції в Excel
- •3. Множення матриць.
- •§ 4. Метод простої ітерації для слр
- •§ 5. Метод Зейделя
- •Питання, тести
- •Завдання
- •Розділ 6. Методи лінійного програмування
- •§ 1. Оптимізаційні задачі. Математичне програмування
- •§ 2. Геометричний зміст задач лінійного програмування. Графічний метод
- •§3. Канонічна форма задачі лінійного програмування. Опорні розв’язки
- •§4. Симплекс – таблиця
- •§5. Симплекс – метод.
- •§6. Розв’язування задач лінійного програмування за допомогою excel
- •§7 Приклади
- •§8. Пошук початкового опорного розв’язку. Метод штучного базису
- •Властивості допоміжної задачі.
- •Питання, тести
- •Завдання
- •Розділ 7. Чисельні методи розв’язування звичайних диференціальних рівнянь
- •§ 1. Метод Ейлера
- •§ 2. Метод Ейлера за допомогою Excel
- •§ 3. Методи Рунге – Кутта
- •§ 4. Подвійний перерахунок для методів Рунге – Кутта
- •§ 5. Кратний перерахунок для методів Рунге – Кутта за допомогою Excel
- •§ 6. Методи Рунге – Кутта з вищими порядками похибки
- •Питання, тести
- •Завдання
- •Іменний покажчик
- •Предметний покажчик
- •Література
Завдання
1.
Розв’язати СЛР
методом Гаусса.
2. Знайти
матрицю X =
– 5
,
використавши
вбудовані оператори
Excel.
3. Знайти
матрицю, обернену до матриці А =
методом Гаусса.
4. Розв’язати СЛР матричним методом, використавши вбудовані оператори Excel.
5. Розв’язати CЛР методом простої ітерації.
Розділ 6. Методи лінійного програмування
§ 1. Оптимізаційні задачі. Математичне програмування
Нехай f(x) – функція n змінних, визначена на множині V Rn , підмножина Ω V. Оптимізаційна задача задається трійкою (f, V, Ω) , де f(x) називається цільовою функцією, а Ω – множиною допустимих розв’язків оптимізаційної задачі; вона полягає у пошуку мінімуму (найменшого) або максимуму (найбільшого) значення цільової функції f(x) на допустимій множині Ω.
Означення 1. Оптимальним розв’язком задачі (f, V, Ω) називається таке х0 є Ω, що f(x0) ≤ f(x) (або f(x0) ≥ f(x) ) для всіх х є Ω. Оптимізаційну задачу називають нерозв’язною, якщо вона не має оптимальних розв’язків. Розв’язати оптимізаційну задачу означає або знайти всі її оптимальні розв’язки, або ж довести її нерозв’язність.
Зокрема, задача мінімізації (максимізації) (f, V, Ω) буде нерозв’язною, якщо f(x) не обмежена знизу (зверху) на допустимій множині Ω.
Методи розв’язування оптимізаційної задачі (f, V, Ω) часто називають методами математичного програмування. Сам термін програмування в цьому випадку зумовлений тим, що саме математичне програмування є підґрунтям мікроекономіки – одної з двох найважливіших складових економічної науки і шуканою у таких задачах насправді є деяка програма дій. Необхідною передумовою будь – якого використання методів обчислень є створення математичної моделі досліджуваного об’єкта чи явища. Першим кроком при розв’язуванні задач економіки є побудова їх економіко – математичної моделі, тобто такої математичної задачі, до розв’язання якої зводиться поставлена економічна задача. Розглянемо деякі приклади.
1. Для виготовлення продукції двох типів (А і Б) на заводі використовують сировину трьох видів (1, 2 и 3). Для виробництва кожного виробу типу А витрачається 10 одиниць сировини виду 1, 20 одиниць виду 2 і 15 виду 3; для кожного виробу типу Б – 20, 10 и 15 відповідно. Запаси сировини складають: 1 виду 100 одиниць, 2 виду – 100, 3 виду – 90. Скласти щоденний план виробництва, що забезпечує найбільший прибуток, якщо кожний виріб типу А дає 200 гривен прибутку, Б – 300 гривен.
Зручно дані поставленої задачі сформувати у наступну таблицю.
Вироби |
Прибуток |
Сировина |
||
1 |
2 |
3 |
||
А |
200 |
10 |
20 |
15 |
Б |
300 |
20 |
10 |
15 |
Запаси сировини |
100 |
100 |
90 |
Побудуємо математичну модель цієї задачі. Позначимо через хА та хБ кількість одиниць продукції А і Б відповідно, яку передбачається виготовити по плану. За умовою задачі цільова функція f – прибуток, який треба максимізувати: f = 200хА + 300хБ. При цьому сировини 1 буде витрачено 10хА + 20хБ. Кількість витраченої сировини не може перевищувати її запасу, тому 10хА + 20хБ ≤ 100. Аналогічно для сировини 2 20хА + 10хБ ≤ 100, для сировини 3 15хА + 15хБ ≤ 90. Змінні хА і хБ повинні бути невід’ємними, оскільки це кількості продукції. Отже економіко – математична модель формулюється так: треба знайти такі значення змінних хА та хБ , за яких задовольняється система нерівностей
і цільова функція f = 200 хА + 300 хБ набуває найбільшого можливого значення.
Знайти такі хА та хБ і означає знайти план виробництва. Це простіший приклад задачі планування виробництва, яка найбільш традиційно формулюється так у загальному вигляді. Підприємство виробляє деяку однорідну продукцію і використовує для цього n різних технологічних засобів і m різних ресурсів. За одиницю часу використання j – го технологічного засобу виготовляється сj одиниць продукції, при цьому витрати і – го ресурсу aij одиниць (і = 1, 2, …, m, j = 1, 2, …, n). Треба так спланувати роботу підприємства, щоби кількість продукції була найбільшою при умові, що витрати і – го ресурсу не можуть перебільшувати bi одиниць (і = 1, 2, …, m).
Нехай
хj
≥ 0 (j = 1, 2,
…, n) –
час використання
підприємством j
– го технологічного засобу. Замість
прибутку у попередньому прикладі тут
цільова функція f
– це кількість виробленої
продукції: f
=
;
все інше так само, лише у загальному
вигляді. За планом буде витрачено
одиниць і–го ресурсу. Звідси
≤
bi
(і = 1, 2, …, m).
Отже ми отримуємо таку економіко –
математичну модель: треба знайти такі
значення змінних хj
(j = 1, 2, …,
n), за яких
задовольняється система нерівностей
і цільова функція f = набуває найбільшого можливого значення.
2. Задача про раціон. У господарстві є два вида кормів вартістю 20 та 30 гривен за одиницю корму відповідно. У першому кормі міститься 2 одиниці вітаміну А та 3 одиниці вітаміну В, у другому 5 одиниць А та 2 одиниці В. Раціон повинен містити не менш як 9 одиниць А та 8 одиниць В. Скласти найдешевший раціон, який задовольняє цим вимогам.
Через хА і хВ позначимо відповідно кількість першого і другого кормів у раціоні. Тут цільова функція f – це вартість раціону: f = 20хА + 30хВ. За таким планом раціон містить 2хА + 5хВ одиниць вітаміну А, за вимогою задачі 2хА + 5хВ ≥ 9. Аналогічно для вітаміну В 3хА + 2хВ ≥ 8. Змінні хА і хВ повинні бути невід’ємними, оскільки це кількості продукції. Отже математична модель формулюється так: треба знайти такі значення змінних хА та хВ , за яких задовольняється система нерівностей
і цільова функція f = 20хА + 30хВ набуває найменшого можливого значення.
У загальному вигляді задача про раціон традиційно формулюється так. У господарстві є n різних кормів, кожен з яких містить m різних поживних речовин. Відомо, що одиниця j – го корму коштує сj грошових одиниць (j = 1, 2, …, n) і містить aij одиниць і – ої поживної речовини (і = 1, 2, …, m). Треба скласти раціон так, щоб він задовольняв мінімальну потребу bi (і = 1, 2, …, m) у кожній речовині при найменшій вартості раціону.
Позначимо кількість j – го корму в раціоні через хj ≥ 0 (j = 1, 2, …, n), цільова функція f – це вартість раціону: f = . За таким планом раціон містить одиниць і – ої поживної речовини, звідки ≥ bi (і = 1, 2, …, m). Отже ми отримуємо таку математичну модель: треба знайти такі значення змінних хj (j = 1, 2, …, n), за яких задовольняється система нерівностей
і цільова функція f = набуває найменшого можливого значення.
3. Транспортна задача. Нехай з трьох пунктів відправлення Р1, Р2 , Р3 треба перевезти однорідний вантаж до трьох пунктів призначення М1, М2 , М3, в тому числі з пункту Р1 – 12 т, з пункту Р2 – 8 т, з пункту Р3 – 10 т. Вантаж повинен надійти за призначенням у пункт М1 – 6 т, пункт М2 – 9 т, пункт М3 – 15 т. Вартість перевезення тони вантажу з пункту Рi до пункту Мj (i, j = 1,2,3) наведена у наступній таблиці.
Пункт відправлення |
Пункт призначення |
||
М1 |
М2 |
М3 |
|
Р1 |
1 |
3 |
4 |
Р2 |
2 |
5 |
3 |
Р3 |
6 |
7 |
4 |
Треба, щоби весь вантаж був вивезений з пунктів відправлення, повністю був задоволений попит пунктів призначення і загальна вартість перевезень була мінімальною.
Побудуємо математичну модель цієї задачі. Позначимо через хij вагу вантажу (в тонах), запланованого для перевезення з пункту Рi до пункту Мj (i, j = 1,2,3). Вартість перевезення багажу дорівнює f = х11 + 3х12 + 4х13 + 2х21 + 5х22 + 3х23 + 6х31 + 7х32 + 4х33 – цільова функція задачі. Весь вантаж буде вивезений з пунктів відправлення, якщо х11 + х12 + х13 = 12, х21 + х22 + х23 = 8, х31 + х32 + х33 = 10. Попит пунктів призначення буде повністю задоволений, якщо х11 + х21 + х31 = 6, х12 + х22 + х32 = 9, х13 + х23 + х33 = 15. Отже, серед невід’ємних розв’язків системи
треба знайти такий, за якого функція f набуває найменшого значення.
У загальному випадку транспортну задачу звичайно формулюють так. Нехай однорідний вантаж, що знаходиться у m постачальників у кількості аi (і = 1, 2, …, m) одиниць відповідно, треба перевезти n споживачам у кількості bj (j = 1, 2, …, n) одиниць. Вартість перевезення одиниці вантажу від і – го постачальника до j – го споживача становить сij одиниць. Спланувати перевезення вантажу так, щоби весь вантаж був вивезений від постачальників, повністю був задоволений попит споживачів і вартість перевезень була мінімальною.
Позначимо
через хij
кількість одиниць вантажу, запланованого
для перевезення від постачальника і
до споживача j.
Математична модель транспортної задачі
має такий вигляд. Знайти найменше
значення функції f
=
за обмежень
У всіх наведених прикладах цільова функція f є лінійною, а множина допустимих розв’язків оптимізаційної задачі визначається системою лінійних рівнянь і нерівностей.
Означення 2. Оптимізаційна задача (f, V, Ω), у якій цільова функція f є лінійною на V, а Ω є множиною розв’язків деякої системи лінійних рівнянь і нерівностей називають задачею лінійного програмування. Систему лінійних рівнянь і нерівностей, яка визначає Ω, називають системою обмежень задачі лінійного програмування.