- •Міністерство освіти і науки України
- •Математичні моделі економічних задач
- •1.1. Задача планування виробництва
- •1.2. Задача складання раціону (задача про дієту, задача про суміші)
- •1.3. Транспортна задача
- •1.4. Задача про мінімізацію відходів
- •К 2ількість шматків
- •1.5. Задача про призначення
- •Загальна постановка задач лінійного програмування (лп)
- •Перелік питань для самоперевірки
- •Лекція 2
- •Тема 2. Геометрична інтерпретація задач лінійного програмування. Задача лінійного програмування, форми її запису
- •Приведення задачі лп до канонічного виду
- •Приведення задачі лп до симетричного виду
- •Перелік питань для самоперевірки
- •3.1. Визначення вихідного опорного плану
- •3.2. Симплексні таблиці
- •3.3. Поняття про м-метод
- •Перелік питань для самоперевірки
- •Лекція 4
- •Тема 4. Двоїстість у лінійному програмуванні
- •Перелік питань для самоперевірки
- •Лекція 5
- •Тема 5. Методика розв’язування транспортної задачі
- •5.1. Приведення задачі до замкненої форми
- •5.2. Визначення вихідного опорного плану
- •5.3. Метод потенціалів
- •Перелік питань для самоперевірки
- •6.1. Метод відсікань Гоморі
- •Перелік питань для самоперевірки
- •Лекція 6
- •Тема 7. Елементи теорії ігор
- •7.1. Графічний метод
- •7.2. Приведення матричної гри до задачі лінійного програмування
- •Перелік питань для самоперевірки
- •8.2. Задачі нелінійного програмування з нелінійною цільовою функцією та лінійною системою обмежень
- •Перелік питань для самоперевірки
- •Лекція 8
- •Тема 9. Динамічне програмування
- •9.1. Задача про розподіл коштів між підприємствами
- •Рішення
- •9.2. Задача про заміну обладнання
- •Рішення
- •Перелік питань для самоперевірки
- •Список рекомендованої літератури
3.2. Симплексні таблиці
Алгоритм знаходження оптимального плану за допомогою симплекс-таблиць:
Дані задачі, отримані після приведення її системи обмежень до невід’ємного одиничного базису (після визначення опорного плану), заносимо в симплекс-таблицю. У самому верхньому рядку таблиці записуємо імена стовпців Б, cб, b і змінні xj (j=1,..,n) задачі, указуючи під ними відповідні коефіцієнти cj (j=1,..,n) цільової функції. У робочу (виділену) частину таблиці заносимо коефіцієнти aij (i=1,..,m; j=1,..,n) при змінних із системи обмежень. У першому стовпці таблиці записуємо базисні змінні (порядок їхнього запису визначається положенням одиниці відповідного одиничного вектора-стовпця), у другому – коефіцієнти цільової функції при цих базисних змінних, у третьому – вільні члени системи. У першій клітині нижнього рядка записуємо значення цільової функції на даному опорному плані (сума добутків елементів стовпця cб на відповідні елементи стовпця b), в інших n клітинах – оцінки оптимальності. Оцінка оптимальності , що відповідає зміннійxj, дорівнює сумі добутків елементів стовпця cб на відповідні елементи aij стовпця коефіцієнтів при xj мінус cj.
Перевіряємо виконання критерію оптимальності. Якщо всі оцінки оптимальності недодатні (), то даний опорний план оптимальний – рішення закінчене. Якщо серед оцінок є додатні, то будуємо новий опорний план.
Найбільша додатна оцінка оптимальності визначає вільну змінну xs, що вводиться до базису (ключовий стовпець s). Знаходимо мінімальне значення із симплексних відносин: . Якщо мінімуму нема, то задача не має розв’язку і обчислення закінчені. Якщо мінімум є, то обираємо рядок, на якому він досягається (будь-який, якщо їх декілька). Цей рядок є ключовим. Він вказує на базисну змінну, що виводиться з базису. За допомогою елементарних перетворень одержуємо новий опорний план.
Переходимо до п. 1 алгоритму.
Задача 3.2. Розв’язати симплексним методом задачу планування виробництва, математичну модель якої було побудовано в задачі 1.1 (с. 5).
Рішення
1. Приведемо дану задачу до канонічного виду. Для цього до лівої частини кожної нерівності додамо невід’ємні змінні x5, x6, x7 і змінимо знак функції мети на протилежний. Таким чином одержимо задачу
2. Знайдемо вихідний опорний план.
Будуємо розширену матрицю системи обмежень
.
Ця матриця містить одиничний базис – три попарно різних одиничних вектори. Тому базисні змінні – це x5, x6, x7; вільні змінні – x1, x2, x3, x4. Базисний розв’язок: . Він є вихідним опорним планом (усі його компоненти невід’ємні).
Перевіримо даний опорний план на оптимальність.
Заповнюємо першу симплексну таблицю.
Таблиця 3.1
Б |
cб |
b |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
–12 |
–7 |
–18 |
–10 |
0 |
0 |
0 | |||
x5 |
0 |
18 |
1 |
2 |
1 |
0 |
1 |
0 |
0 |
x6 |
0 |
30 |
1 |
1 |
2 |
1 |
0 |
1 |
0 |
x7 |
0 |
40 |
1 |
3 |
3 |
2 |
0 |
0 |
1 |
|
0 |
12 |
7 |
18 |
10 |
0 |
0 |
0 |
В нижньому рядку першої клітини записано значення цільової функції на опорному плані : f() = 0 · 18 + 0 · 30 + 0 · 40 = 0 в інших 7 клітинах – оцінки оптимальності : c1 = 0 · 1 + 0 · 1 + 0 · 1 – (–12) = 12; c1 = 0 · 2 + 0 · 1 + 0 · 3 – (–7) = 7; ; ;
; ;.
Опорний план не є оптимальним, тому що серед оцінок в останньому рядку є додатні оцінки: 12, 7, 18, 10.
3. Переходимо до нового опорного плану. Змінну x3 вводимо до базису, оскільки цій змінній відповідає максимальна додатна оцінка 18. Змінну x7 виводимо з базису, тому що їй відповідає мінімальне значення із симплексних відносин: . За допомогою елементарних перетворень у третьому стовпці розширеної матриці будуємо одиничний вектор, в якому 1 знаходиться в ключовому рядку III.
.
Базисні змінні: x3, x5, x6. Вільні змінні: x1, x2, x4, x7. Базисний розв’язок – новий опорний план.
Заповнюємо другу симплекс-таблицю.
Таблиця 3.2
Б |
cб |
b |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
–12 |
–7 |
–18 |
–10 |
0 |
0 |
0 | |||
x5 |
0 |
14/3 |
2/3 |
1 |
0 |
–2/3 |
1 |
0 |
–1/3 |
x6 |
0 |
10/3 |
1/3 |
–1 |
0 |
–1/3 |
0 |
1 |
–2/3 |
x3 |
–18 |
40/3 |
1/3 |
1 |
1 |
2/3 |
0 |
0 |
1/3 |
|
–240 |
6 |
–11 |
0 |
–2 |
0 |
0 |
–6 |
Опорний план не є оптимальним, тому що в останньому рядку є додатна оцінка 6.
4. Переходимо до нового опорного плану. Змінну x1 вводимо до базису, оскільки цій змінній відповідає додатна оцінка 6. Змінну x5 виводимо з базису, тому що цій змінній відповідає мінімальне значення із симплексних відносин: . За допомогою елементарних перетворень у першому стовпці останньої матриці будуємо одиничний вектор, в якому 1 розташована в ключовому рядку I.
.
Базисні змінні: x1, x3, x6. Вільні змінні: x2, x4, x5, x7. Базисний розв’язок – новий опорний план.
Заповнюємо третю симплекс-таблицю.
Таблиця 3.3
Б |
cб |
b |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
–12 |
–7 |
–18 |
–10 |
0 |
0 |
0 | |||
x1 |
–12 |
7 |
1 |
3/2 |
0 |
–1 |
3/2 |
0 |
–1/2 |
x6 |
0 |
1 |
0 |
–3/2 |
0 |
0 |
–1/2 |
1 |
–1/2 |
x3 |
–18 |
11 |
0 |
1/2 |
1 |
1 |
–1/2 |
0 |
1/2 |
|
–282 |
0 |
–20 |
0 |
4 |
–9 |
0 |
–3 |
Опорний план не є оптимальним, тому що є додатна оцінка 4.
5. Переходимо до нового опорного плану. Змінну x4 вводимо до базису, оскільки цій змінній відповідає додатна оцінка 4. Змінну x3 виводимо з базису, тому що цій змінній відповідає мінімальне значення із симплексних відносин: . За допомогою елементарних перетворень у четвертому стовпці останньої матриці будуємо одиничний вектор, в якому 1 знаходиться в ключовому рядку III.
.
Базисні змінні: x1, x4, x6. Вільні змінні: x2, x3, x5, x7. Базисний розв’язок – новий опорний план.
Заповнюємо четверту симплекс-таблицю.
Таблиця 3.4
Б |
cб |
b |
x1 |
x2 |
x3 |
x4 |
x5 |
x6 |
x7 |
–12 |
–7 |
–18 |
–10 |
0 |
0 |
0 | |||
x1 |
–12 |
18 |
1 |
2 |
1 |
0 |
1 |
0 |
0 |
x6 |
0 |
1 |
0 |
–3/2 |
0 |
0 |
–1/2 |
1 |
–1/2 |
x4 |
–10 |
11 |
0 |
1/2 |
1 |
1 |
–1/2 |
0 |
1/2 |
|
–326 |
0 |
–22 |
–4 |
0 |
–7 |
0 |
–5 |
Опорний план є оптимальним, тому що в останньому рядку всі оцінки недодатні.
Оптимальне значення цільової функції: .
Повернемося до вихідної задачі: ,.
Відповідно до отриманого результату підприємство дістане прибуток у розмірі 326 грн, якщо виготовить 18 одиниць продукції першого виду і 11 одиниць продукції четвертого виду, а продукцію другого і третього видів взагалі не буде випускати.
Зауваження 3.2. Задача має єдиний оптимальний розв’язок, якщо всі оцінки, що відповідають вільним змінним, від’ємні. Якщо серед таких оцінок є нульові, то розв’язок не єдиний. Безліч оптимальних розв’язків описується за допомогою опуклої комбінації двох розв’язків:
,
де інший оптимальний план виконується за допомогою введення до базису вільної змінної, якій відповідає нульова оцінка оптимальності.
Зауваження 3.3. (Геометрична інтерпретація симплексного методу). Опорному плану відповідає кутова точка многокутника розв’язків системи обмежень. Ідея симплексного методу – перебір кутових точок многокутника розв’язків з урахуванням зміни значення функції мети (кожний наступний розв’язок “кращий” за попередній, тобто при ). Такий перебір дозволяє скоротити число кроків при відшукуванні оптимуму.