- •Економічна інформатика
- •1) Лінійні і нелінійні задачі:
- •2) Дискретні та неперервні задачі:
- •3) Детерміновані та стохастичні (ймовірностні) задачі:
- •4) Статичні (однокрокові) та динамічні (багатокрокові) задачі:
- •Зміст виконання завдання
- •Зміст виконання завдання
- •Критерій оптимальності – мінімум затрат праці - запишемо як
- •Зміст виконання завдання
- •Скласти план вантажних перевезень з мінімальним вантажообігом.
- •Втрати живої ваги при перевезенні худоби, кг на 1 т
- •Площі попередників озимої пшениці, га
- •Площа сортів озимої пшениці, га
- •Середня урожайність озимої пшениці за попередниками, ц з 1 га
- •Зміст виконання завдання
- •Зм 3. Двоїстість у лінійному програмуванні
- •Зм 4. Цілочислове програмування
- •6.1. Метод відтинання Гоморі
- •6.2. Метод гілок і меж
- •Зм 5. Оптимізаційні економіко-математичні моделі
- •7.1. Моделювання виробничих систем в рослинництві
- •7.2. Моделювання виробничих систем в тваринництві
- •7.3. Моделювання виробництва і реалізації продукції
- •Зміст виконання завдання
- •Зм 6. Оптимізаційні задачі управління запасами Завдання 8. Детерміновані та стохастичні моделі управління запасами
- •Детермінована статична однономенклатурна модель управління запасами без дефіциту
- •Стохастична модель управління запасами за умови, що попит характеризується нормальним законом розподілу
- •Стохастична модель управління запасами за умови штрафу за дефіцит
- •Зм 7. Моделі задач масового обслуговування
- •Зміст виконання завдання
- •Зм 8. Задачі упорядкування та координації
- •Зміст виконання завдання
- •Зм 9. Задачі та моделі заміни обладнання Завдання 11. Моделювання заміни обладнання
- •Отже, рекурентне співвідношення для періоду т буде мати вигляд:
- •Якщо обладнання після списання реалізується, то рекурентне свіввідношення має вигляд
- •Зміст виконання завдання
- •Зм 10. Задачі з умовами невизначеності та конфлікту
- •Задачі для самостійного розв’язування
- •Зміст виконання завдання
- •Зміст виконання завдання
- •Зм 11. Багатокритеріальні задачі
- •Зміст виконання завдання
- •Додаток а Приклад використання надбудови SimplexWin для розв’язування задач лінійного програмування в симплексних таблицях
- •Додаток б Приклад використання Excel для розв" язання симплексних задач лінійного програмування за допомогою надбудови "Поиск решения"
- •Додаток в Приклад використання Excel для розв’язання транспортних задач лінійного програмування (тзлп) за програмою "Поиск решения"
- •Список рекомендованої літератури Підручники та навчальні посібники
- •Електронні ресурси
- •Марченко Володимир Петрович Економічна інформатика
Втрати живої ваги при перевезенні худоби, кг на 1 т
М'ясокомбінати |
Господарства | ||||
S1 |
S2 |
S3 |
S4 |
S5 | |
D1 |
20 |
15 |
25 |
10 |
50 |
D2 |
15 |
10 |
30 |
20 |
30 |
D3 |
25 |
30 |
15 |
35 |
55 |
Задача 4.2.Потрібно скласти такий план розміщення сортів озимої пшениці за попередниками, щоб очікуваний валовий збір зерна був максимальним. Площі попередників та посівні площі сортів озимої пшениці наведені в таблицях 4.4 та 4.5, а середня урожайність сортів озимої пшениці за попередниками в таблиці 4.6.
Таблиця 4.4
Площі попередників озимої пшениці, га
Варіант (за передостанньою цифрою шифру Р) |
Попередники | ||||
Ярі зернобобові |
Кукурудза на зелений корм |
Багаторічні трави |
Насінники багаторічних трав |
Ранні овочі | |
S1 |
S2 |
S3 |
S4 |
S5 | |
0 |
300 |
300 |
300 |
50 |
50 |
1 |
200 |
250 |
400 |
50 |
100 |
2 |
200 |
150 |
500 |
100 |
50 |
3 |
300 |
100 |
400 |
100 |
100 |
4 |
250 |
350 |
300 |
60 |
40 |
5 |
400 |
150 |
250 |
80 |
120 |
6 |
350 |
200 |
350 |
50 |
50 |
7 |
300 |
150 |
450 |
40 |
60 |
8 |
300 |
200 |
350 |
70 |
80 |
9 |
250 |
250 |
350 |
60 |
90 |
Таблиця 4.5
Площа сортів озимої пшениці, га
Варіант (за останньою цифрою шифру К) |
Сорти озимої пшениці | |||
Тіра |
Одеська-267 |
Донецька -48 |
Коломак-5 | |
D1 |
D2 |
D3 |
D4 | |
0 |
250 |
450 |
200 |
100 |
1 |
300 |
300 |
250 |
150 |
2 |
400 |
300 |
100 |
200 |
3 |
350 |
300 |
250 |
100 |
4 |
300 |
250 |
300 |
150 |
5 |
100 |
350 |
300 |
250 |
6 |
150 |
400 |
150 |
300 |
7 |
200 |
400 |
200 |
200 |
8 |
350 |
200 |
350 |
100 |
9 |
450 |
100 |
200 |
250 |
Таблиця 4.6
Середня урожайність озимої пшениці за попередниками, ц з 1 га
Попередники |
Сорти озимої пшениці | ||||
Тіра |
Одеська-267 |
Донецька 48 |
Коломак 5 | ||
D1 |
D2 |
D3 |
D4 | ||
Ярі зернобобові |
S1 |
47 |
48 |
49 |
52 |
Кукурудза на зелений корм |
S2 |
44 |
47 |
48 |
45 |
Багаторічні трави |
S3 |
52 |
56 |
51 |
54 |
Насінники багаторічних трав |
S4 |
50 |
56 |
48 |
54 |
Ранні овочі |
S5 |
46 |
49 |
47 |
48 |
Зміст виконання завдання
1. Запис умов задачі за індивідуальним варіантом.
2. Побудова початкового опорного плану задачі.
3. Розв'язок задачі в транспортних таблицях та на ПЕОМ за допомогою надбудови MS Excel ”Поиск решения” (додаток В) .
4. Висновки за результатами розв'язку задачі.
Зм 3. Двоїстість у лінійному програмуванні
Завдання 5. Двоїсті задачі лінійного програмування
Теоретична частина. Кожну пряму задачу лінійного програмування можна перетворити у двоїсту за такими правилами:
Кожному обмеженню прямої задачі відповідає одна змінна двоїстої задачі.
Кількість обмежень двоїстої задачі дорівнює кількості змінних прямої задачі.
Коефіцієнти при змінних двоїстої задачі отримуються шляхом транспонування матриці коефіцієнтів при змінних в обмеженнях прямої задачі.
Якщо в прямій задачі обмеження мають вигляд менше або дорівнює () і на змінні накладена умова невід’ємності, то обмеження двоїстої задачі мають вигляд більше або дорівнює (), а на її змінні також накладено умову невід’ємності і навпаки.
Обсяги обмежень прямої задачі є оцінками змінних двоїстої задачі.
Оцінки змінних прямої задачі є обсягами обмежень двоїстої задачі.
7. Пошук максимуму цільової функції прямої задачі замінюється на пошук мінімуму цільової функції двоїстої задачі і навпаки.
Порівнюючи оптимальні розв'язки прямої і двоїстої задач можна помітити, що розв'язок прямої задачі одночасно дає і розв'язок двоїстої задачі, а саме: значення цільових функцій прямої і двоїстої задач однакові Zmax=Wmin, а коефіцієнти при змінних в рядку цільової функції прямої задачі дорівнюють значенням базисних змінних двоїстої задачі. При цьому коефіцієнти при додаткових змінних у рядку цільової функції прямої задачі дорівнюють значенням основних змінних двоїстої задачі, а коефіцієнти цільової функції при основних змінних прямої задачі дорівнюють значенням додаткових змінних двоїстої задачі.
Для аналізу, як правило, використовується кінцева симплексна таблиця (оптимальний розв'язок) прямої задачі, у якій коефіцієнти при невідомих в рядку цільової функції є двоїстими оцінками, а коефіцієнти при невідомих в стовпчиках – коефіцієнтами структурних зрушень. При цьому двоїсті оцінки базисних змінних завжди дорівнюють нулю, а коефіцієнтами структурних зрушень на перетині однойменних стовпчика і рядка дорівнюють одиниці, а інші нулю, і в економіко-математичному аналізі оптимальних розв'язків не використовуються
Доведено, що при введенні до оптимального розв'язку одиниці небазисної змінної двоїсті оцінки показують величину зміни значення цільової функції в одиницях виміру цільової функції, а коефіцієнти структурних зрушень в колонках небазисних змінних показують величину зміни значень відповідних базисних невідомих в одиницях виміру цих змінних. При цьому для задач максимізації:
1) при введенні в базис одиниці основної небазисної змінної додатні коефіцієнти структурних зрушень показують величину зменшення базисних змінних, а від’ємні - збільшення. Значення цільової функції зменшується на величину двоїстої оцінки;
2) при введенні в базис одиниці додаткової небазисної змінної при обмеженнях виду менше або дорівнює () додатні коефіцієнти структурних зрушень показують величину збільшення базисних змінних, а від’ємні - зменшення. Значення цільової функції збільшується на величину двоїстої оцінки;
3) при обмеженнях виду більше або дорівнює () додатні коефіцієнти структурних зрушень показують величину зменшення базисних змінних, а від’ємні - збільшення. Значення цільової функції зменшується на величину двоїстої оцінки.
Для задач мінімізації цільової функції знаки змінюються на зворотні.
Слід відзначити, що наведені показники дійсні тільки для умов конкретної задачі і мають певні межі застосування. Для визначення максимальних меж введення в базис небазисних змінних потрібно значення базисних змінних поділити на відповідні коефіцієнти структурних зрушень і вибрати найменше відношення. При цьому в задачах максимізації цільової функції для основних небазисних змінних значення базисних змінних потрібно поділити на додатні коефіцієнти структурних зрушень, для додаткових небазисних змінних при обмеженнях виду менше або дорівнює () – на від’ємні, а для додаткових небазисних змінних при обмеженнях виду більше або дорівнює (≥)– знов на додатні.. Для задач мінімізації цільової функції знаки змінюються на зворотні.
Властивості коефіцієнтів структурних зрушень і двоїстих оцінок дозволяють без повторного розв’язання симплексної задачі отримати нові варіанти оптимального розв'язку шляхом введення в базис (в розрахованих межах) небазисних невідомих.
Розрахунки виконуються за формулою
Xj* Xj´ Ij ∆Xj
де Xj* - значення базисної змінної в новому варіанті оптимального розв"язку;
Xj´ - значення базисної змінної у початковому оптимальному розв"язку;
ij - коефіцієнти структурних зрушень та двоїсті оцінки небазисної змінної, яка вводиться в оптимальний розв’язок;
∆Xj - величина небазисної змінної, яка вводиться в оптимальний розв’язок.
Використовуючи властивості двоїстих оцінок можна також визначити доцільність виробництва нового виду продукції. Для цього потрібно затрати ресурсів необхідні для виробництва одиниці продукції перемножити на двоїсті оцінки цих ресурсів. Якщо загальна оцінка всіх ресурсів менша ціни одиниці цієї продукції, то виробляти таку продукцію доцільно.
Контрольні питання
1. За якими правилами пряму задачу лінійного програмування можна перетворити у двоїсту ?
2. Де знаходяться двоїсті оцінки та коефіцієнтів структурних зрушень в кінцевій симплексній таблиці прямої задачі ?
3. Як впливають двоїсті оцінки на значення цільової функції ?
4. Як впливають коефіцієнти структурних зрушень на значення відповідних базисних змінних ?
5. Як визначаються максимальні межі введення в базисний розв'язок небазисних змінних ?
6. Як отримати новий варіант оптимального розв'язку задачі за допомогою двоїстих оцінок та коефіцієнтів структурних зрушень ?
7. Як визначити доцільність введення в оптимальний розв'язок нової небазисної змінної ?
Приклад 5.1. За умовами прикладу 2.1 записати двоїсту до неї задачу і розв'язати її за допомогою симплексного методу. Використовуючи кінцеву симплексну таблицю прямої задачі провести економіко-математичний аналіз, а саме: пояснити економічну суть двоїстих оцінок і коефіцієнтів структурних зрушень небазисних змінних.
Розв'язання. Пряма задача має такий математичний запис:
Zmax = 90x1 + 50x2
при обмеженнях: 1) x1 + x2 < 400; 2) 50x1+30x2 < 18000; 3) 6x1+ 3x2 < 1800
Тоді, двоїста задача буде мати вигляд:
Wmin=400y1+18000y2+1800y3
при обмеженнях: 1) у1 + 50у2 + 6у3 90;2) у1 + 30у2 + 3у3 50
Розв'язок двоїстої задачі: Wmin = 28000, у1 =10; у2 =0; у3= 13,3; r1 = 0; r2 = 0,
де r1 і r2 – додаткові змінні двоїстої задачі.
Приведемо оптимальний розв'язок (кінцеву симплексну таблицю) прямої задачі (див. приклад 2.1).
Базисні змінні |
Значення базисних змінних |
x1 |
x2 |
s1 |
s2 |
s3 | ||||||
x2 |
200 |
0 |
1 |
2 |
0 |
-0,333 | ||||||
s2 |
2000 |
0 |
0 |
-10 |
1 |
-6,667 | ||||||
x1 |
200 |
1 |
0 |
-1 |
0 |
0,333 | ||||||
Zmax |
28000 |
0 |
0 |
10 |
0 |
13,333 |
Нагадаємо, що:
х1 – площа кукурудзи на зерно, га;
х2 – площа ячменю, га;
s1 – площа невикористаної ріллі, га;
s2 – кількість невикористаних трудових ресурсів, люд.-год.;
s3 – кількість невикористаних мінеральних добрив, ц;
Zmax - обсяг виробництва кормів, ц к. од.
Порівняння оптимальних розв'язків прямої (рядок Zmax симплексної таблиці) і двоїстої задач показує, що значення цільових функцій прямої і двоїстої задач однакові Zmax = Wmin=28000; коефіцієнти при змінних в Zmax-рядку прямої задачі дорівнюють значенням базисних змінних двоїстої задачі. При цьому: y1 = 10; y2 = 0;
у3 = 13.333; r1 = 0; r2 = 0. Аналіз кінцевої симплексної таблиці прямої задачі показує, що двоїсті оцінки базисних змінних х1, х2 та s2 дорівнюють нулю. Двоїста оцінка небазисної змінної s1 дорівнює 10. Це означає, що кожний додатковий гектар ріллі дасть додатково 10 ц к. од. (відзначимо, що за результатами розв'язку прямої задачі s1 = 0, тобто всі 400 га ріллі були використані повністю).
Двоїста оцінка змінної s3 дорівнює 13,333. Це означає, що кожний додатковий центнер мінеральних добрив дасть змогу збільшити виробництво кормів на 13,3 ц к. од. (за результатами розв'язку прямої задачі мінеральні добрива також використовуються повністю, s3 = 0). Отже двоїсті оцінки додаткових змінних при обмеженнях виду менше або дорівнює () показуютьміру дефіцитності виробничих ресурсів. Дійсно, трудові ресурси повністю не використані (s2=2000), тому двоїста оцінка цієї змінної дорівнює нулю. Двоїсті оцінки основних небазисних змінних показують неефективність введення цих змінних в базисний розв'язок.
Коефіцієнти структурних зрушень для змінної s1 (2; -10; -1) означають, що коли у виробництво додатково залучити 1 га ріллі, то площа ячменю (х2) збільшиться на 2 га, а площа кукурудзи (x1) зменшиться на 1га, залишок трудових ресурсів в (s2) зменшиться на 10 люд.-год.. Відповідно при залученні у виробництво 1 ц мінеральних добрив (s3) площа кукурудзи (х1) збільшиться на 0,333 га, а площа ячменю (х2) зменшиться на 0,333 га. Залишок трудових ресурсів (s2) зменшиться на 6,667 люд.-год..
Максимальна межа введення в базисний розв'язок небазисної змінної S1 дорівнює 200 (min2000/10=200; 200/1=200), а для небазисної змінної S3 - 300 min200/0,333=600; 2000/6,667=300).
Приклад 5.2. За умовами прикладу 2.1 визначити новий варіант оптимального розв'язку задачі при додатковому залученні у виробництво 10 га ріллі або S1= +10.
Розв’язання. Для розв’язання задачі використаємо двоїсту оцінку і коефіцієнти структурних зрушень небазисної змінної S1 кінцевої симплексної таблиці прямої задачі. Розрахунки виконаємо в табл. 5.1.
Таблиця 5.1
Базисні змінні |
Значення базисних змінних |
Коефіцієнти структурних зрушень і двоїста оцінка s1 |
Добутки коефіцієнтів структурних зрушень і двоїстої оцінки при s1 на 10 |
Новий варіант оптимального розв'язку задачі*) |
x2 |
200 |
2 |
210 = 20 |
200+20=220 |
s2 |
2000 |
-10 |
-1010 =-100 |
2000+(-100)=1900 |
x1 |
200 |
-1 |
-110 = -10 |
200+(-10)=190 |
Zmax |
28000 |
10 |
1010 = 100 |
28000+100=28100 |
*)Примітка. Для основних та додаткових змінних при обмеженнях виду більше або дорівнює () добутки віднімаються.
Висновки. В новому варіанті оптимального розв'язку задачі площа кукурудзи (х1) дорівнює 190 га, ячменю (х2) – 220 га, невикористані трудові ресурси (S2) складають 1900 люд.-год.. Виробництво кормів становить 28100ц к. од. (Zmax =28100).
Приклад 5.3. За умовами прикладу 2.1 визначити доцільність вирощування гороху при затратах праці 35 люд.-год. і добрив 2 ц на 1 га та при виході кормів 40 ц к. од. з 1 га.
Розв’язання. Для визначення доцільності вирощування гороху потрібно обсяги виробничих ресурсів, які затрачені на 1 га гороху, перемножити на двоїсті оцінки цих ресурсів. Якщо сумарна оцінка цих ресурсів менша за вихід кормових одиниць з 1 га гороху, то вирощування гороху рентабельне: 1х10 + 0х35 + 2х13,3 = 36,6 < 40.
Висновки. Вирощування гороху рентабельне.
Задачі для самостійного розв'язання
Задача 5. За умовами задачі 2.1 для самостійного розв"язування:
а) записати до неї двоїсту задачу і розв'язати її на ПЕОМ за допомогою надбудови MS Excel ”Поиск решения” (додаток А) ;
б) використовуючи кінцеву симплексну таблицю задачі 2.1 провести економіко-математичний аналіз, а саме: *)
- пояснити економічну суть двоїстих оцінок і коефіцієнтів структурних зрушень небазисних змінних;
- визначити максимальні межі введення в базисний розв'язок небазисних змінних;
- отримати нові варіанти оптимальних розв'язків задачі шляхом введення в базисний розв'язок небазисних змінних (величина небазисної змінної не повинна перевищувати визначених меж);
в) визначити доцільність вирощування гороху при затратах праці 30 + К люд.- год. і мінеральних добрив 2 ц на 1 га та урожайності 30 + Р ц зерна з 1 га .
*)Примітка. Розрахунки виконати для всіх небазисних змінних кінцевої симплексної таблиці прямої задачі.