Симплексний метод розв’язування задач ЛП

Приклади розв’язування задач ЛП симллекс методом

Приклад 2.1.

Продукція чотирьох видів А, В, С і Д проходить послідовну обробку на двох верстатах. Тривалість обробки одиниці продукції кожного виду задано таблицею.

Верстат	Тривалість обробки, год, одиниці продукції
А	В	С	Д

Витрати на виробництво одиниці продукції кожного виду визначають як величини, прямо пропорційні до часу використання верстатів (у машино-годинах). Вартість однієї машино-год становить 10 дол. для верстата 1 і 15 дол. — для верстата 2. Можливий час використання верстатів обмежений: для верстата 1 він становить 450 машино-год, а для верстата 2 — 380 машино-год.

Ціна одиниці продукції кожного виду дорівнює відповідно 73, 70, 55 та 45 дол.

Визначити оптимальний план виробництва продукції всіх чотирьох видів, який максимізує загальний чистий прибуток.

 

Побудова математичної моделі. Нехай х_j — план виробництва продукції j -го виду, де j може набувати значень від 1 до 4.

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

для верстата 1 (машино-год);

для верстата 2 (машино-год).

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

.

Отже, математична модель поставленої задачі має такий вигляд:

;

Розв’язування. Розв’яжемо задачу симплекс-методом згідно з розглянутим алгоритмом.

1. Запишемо систему обмежень задачі в канонічному вигляді. Для цього перейдемо від обмежень-нерівностей до строгих рівнянь, увівши до лівої частини обмежень додаткові змінні х ₅ та х ₆:

Ці додаткові змінні за економічним змістом означають можливий, але не використаний для виробництва продукції час роботи верстатів 1 та 2. У цільовій функції Z додаткові змінні мають коефіцієнти, які дорівнюють нулю:

Канонічну систему обмежень задачі запишемо у векторній формі:

де

Оскільки вектори та одиничні та лінійно незалежні, саме з них складається початковий базис у зазначеній системі векторів. Змінні задачі х ₅ та х ₆, що відповідають одиничним базисним векторам, називають базисними, а решту — вільними змінними задачі лінійного програмування. Прирівнюючи вільні змінні до нуля, з кожного обмеження задачі дістаємо значення базисних змінних:

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

.

Оскільки додатні коефіцієнти х ₅ та х ₆ відповідають лінійно незалежним векторам, то за означенням

є опорним планом задачі і для цього початкового плану

.

2. Складемо симплексну таблицю для першого опорного плану задачі.

Базис	С баз	План				–5			θ
х 1	х 2	х 3	х 4	х 5	х 6
← х 5
х 6
Zj – Cj ≥ 0		–8	–10

Елементи останнього рядка симплекс-таблиці є оцінками ∆ _j, за допомогою яких опорний план перевіряють на оптимальність. Їх визначають так:

;

;

;

;

;

.

У стовпчику «План» оцінкового рядка записують значення цільової функції Z, якого вона набуває для визначеного опорного плану: .

3. Після обчислення всіх оцінок опорний план перевіряють на оптимальність. Для цього продивляються елементи оцінкового рядка. Якщо всі ∆ _j ≥ 0 (для задачі на max) або ∆ _j ≤ 0 (для задачі на min), визначений опорний план є оптимальним. Якщо ж в оцінко­вому ряд­ку присутня хоча б одна оцінка, що не задовольняє умову оптималь­ності (від’ємна в задачі на max або додатна в задачі на min), то опорний план є неоптимальним і його можна поліпшити.

У цій задачі в оцінковому рядку дві оцінки та суперечать умові оптимальності, і тому перший визначений опор­ний план є неоптимальним. За алгоритмом симплекс-методу необхідно від нього перейти до іншого опорного плану задачі.

4. Перехід від одного опорного плану до іншого виконують зміною базису, тобто за рахунок виключення з поточного базису якоїсь змінної та включення замість неї нової з числа вільних змінних.

Для введення до нового базису беремо змінну х ₂, оскільки їй відповідає найбільша за абсолютною величиною оцінка серед тих, які не задовольняють умову оптимальності (│–10│>│–8│).

Щоб визначити змінну, яка підлягає виключенню з поточного базису, для всіх додатних елементів стовпчика «х ₂» знаходимо відно­шен­ня і вибираємо найменше значення. Згідно з даними симплексної таблиці бачимо, що , і тому з базису виключаємо змінну х ₅, а число а ₁₂ = 3 називатимемо розв’язувальним елементом. Подальший перехід до нового опорного плану задачі полягає в побудові наступної симплексної таблиці, елементи якої розраховують за методом Жордана—Гаусса.

Друга симплексна таблиця має такий вигляд:

Базис	С баз	План				–5			θ
х 1	х 2	х 3	х 4	х 5	х 6
х 5			2/3		4/5	2/3	1/3
←х 6			5/3		–5/3	2/3	–2/3
Zj – Cj ≥ 0		–4/3		40/3	35/3	10/3

У цій таблиці спочатку заповнюють два перших стовпчики «Базис» і «С _баз», а решту елементів нової таблиці розраховують за розглянутими далі правилами:

1. Розв’язувальний (напрямний) рядок необхідно поділити на розв’язувальний елемент і здобуті числа записати у відповідний рядок нової симплексної таблиці.

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

3. Якщо в напрямному рядку є нульовий елемент, то відповід­ний стовпчик переписують у нову симплексну таблицю без змін.

4. Якщо в напрямному стовпчику є нульовий елемент, то відповідний рядок переписують у нову таблицю без змін.

Усі інші елементи наступної симплексної таблиці розраховують за правилом прямокутника.

Щоб визначити будь-який елемент нової таблиці за цим правилом, необхідно в попередній симплексній таблиці скласти умов­ний прямокутник, вершини якого утворюються такими числами:

1 — розв’язувальний елемент;

2 — число, що стоїть на місці елемента нової симплексної таблиці, який ми маємо розрахувати;

3 та 4 — елементи, що розміщуються в двох інших протилежних вершинах умовного прямокутника.

Необхідний елемент нової симплекс-таблиці визначають так:

Наприклад, визначимо елемент , який розміщується в новій таблиці в другому рядку стовпчика «х ₄». Складемо умовний прямокутник:

Тоді . Це значення записуємо в стовпчик «х ₄» другого рядка другої симплексної таблиці.

Аналогічно розраховують усі елементи нової симплексної таблиці, у тому числі елементи стовпчика «План» та оцінкового рядка. Наявність двох способів визначення оцінок опорного плану (за правилом прямокутника та за відповідною формулою) дає змогу контролювати правильність арифметичних обчислень на кожному кроці симплекс-методу.

Після заповнення нового оцінкового рядка перевіряємо виконан­ня умови оптимальності Z_j – C_j ≥ 0 для другого опорного плану. Цей план також неоптимальний, оскільки . Використовуючи процедуру симплекс-методу, визначаємо третій опорний план задачі, який наведено у вигляді таблиці:

Базис	С баз	План				–5
х 1	х 2	х 3	х 4	х 5	х 6
х 2						2/5	3/5	–2/5
х 1					–1	2/5	–2/5	3/5
Zj – Cj ≥ 0					61/5	14/5	4/5

В оцінковому рядку третьої симплексної таблиці немає від’ємних чисел, тобто всі ∆ _j ≥ 0 і задовольняють умову оптимальності. Це означає, що знайдено оптимальний план задачі:

або

Х ^* = (48; 118; 0; 0; 0; 0);

.

Отже, план виробництва продукції, що передбачає випуск
48 одиниць продукції А та 118 од. продукції В, оптимальний
і дає найбільший прибуток 1564 дол. При цьому час роботи верстатів використовується повністю (х ₅ = х ₆ = 0).

Задачу можна розв’язати симплекс-методом, узявши не три симплексні таблиці, а одну, в якій послідовно записувати всі ітерації.

 

Приклад 2.2.

Розв’язати попередню задачу із додатковою умовою: продукція С має виготовлятися в кількості не менш як 9 одиниць.

Розв’язування. Математичну модель сформульованої задачі запишемо так:

Застосовуючи для розв’язування поставленої задачі симплекс-метод, спочатку записуємо систему обмежень у канонічній формі, а далі — у векторній:

Зауважимо, що нерівність типу «≥» у рівняння перетворюємо введенням у ліву частину обмеження додаткової змінної зі знаком «–».

Векторна форма запису:

Серед записаних векторів є лише два одиничні — та , а базис у тривимірному просторі має складатися з трьох одиничних векторів. Ще один одиничний вектор можна дістати, увівши в третє обмеження з коефіцієнтом +1 штучну змінну х ₈, якій відповідатиме одиничний вектор .

Тепер можемо розглянути розширену задачу лінійного програмування:

На відміну від додаткових змінних штучна змінна х ₈ має в цільовій функції Z коефіцієнт + М (для задачі на min) або – М (для задачі на max), де М — досить велике додатне число.

У розширеній задачі базисними змінними є х ₅, х ₆, х ₈, а решта змінних вільні. Початковий опорний план задачі:

Складемо першу симплексну таблицю задачі:

Базис	С баз	План				–5				– М	θ
х 1	х 2	х 3	х 4	х 5	х 6	х 7	х 8
х 5											112,5
х 6
←х 8	– М								–1
Zj – Cj ≥ 0		–8	–10
–9 М			– М				М

Розраховуючи оцінки першого опорного плану, дістаємо

Z ₀ = 0 – 9 M, Z ₁ – C ₁ = –8, Z ₂ – C ₂ = –10, Z ₃ – C ₃ = 0 – М і т. д.
Як бачимо, значення оцінок складаються з двох частин, одна з яких містить М, а інша — просто число. Тому для зручності роз­биваємо оцінковий рядок на два. У перший оцінковий рядок за­писуємо просто число, а в другий — число з коефіцієнтом М.

Оцінки першого плану не задовольняють умову оптимальності, і тому він є неоптимальним. Згідно з алгоритмом, розглянутим у за­да­чі 2.41, виконуємо перехід до наступного опорного плану задачі.

Подальше розв’язування задачі наведене у вигляді таблиці:

Базис	С баз	План				–5				– М	θ
х 1	х 2	х 3	х 4	х 5	х 6	х 7	х 8
← х 5										–4
х 6										–1	185,5
х 3									–1		—
Zj – Cj ≥ 0		–8	–10
								М
х 2			2/3			2/3	1/3		4/3	–4/3
← х 6			5/3			2/3	–2/3		–5/3	5/3
х 3									–1		—
Zj – Cj ≥ 0		–4/3			35/3	10/3		40/3	–40/3
								М
х 2						2/5	3/5	–2/5		–2
х 1						2/5	–2/5	3/5	–1
х 3									–1
Zj – Cj ≥ 0					61/5	14/5	4/5		–12
								М

 

Оптимальним планом задачі є вектор

Х ^* = (57; 100; 9; 0; 0; 0; 0),

Отже, оптимальним є виробництво 57 одиниць продукції А, 100 одиниць продукції В і 9 одиниць продукції С. Тоді прибуток буде найбільшим і становитиме 1456 дол.

 

Приклад 2.3.

Продукція фабрики випускається у вигляді паперових рулонів стандартної ширини — 2 м. За спеціальним замовленням споживачів фабрика постачає також рулони інших розмірів, розрізуючи стандартні рулони.

Типові замовлення на рулони нестандартних розмірів наведено в таблиці:

Замовлення	Потрібна ширина рулонів, м	Потрібна кількість рулонів
	0,8
	1,0
	1,2

Визначити оптимальний варіант розкрою стандартних рулонів, за якого спеціальні замовлення, що надходять, задовольняють повністю з мінімальними відходами паперу.

Розв’язування. Аби виконати спеціальні замовлення, які надійшли, розглянемо п’ять можливих варіантів розкрою стандартних рулонів, що можуть використовуватися в різних комбінаціях. Варіанти розкрою наведено в таблиці:

Потрібна ширина рулона, м	Кількість нестандартних рулонів за варіантами
0,8
1,0
1,2
Кількість відходів, м	0,4		0,2		0,8

Нехай x_j — кількість стандартних рулонів паперу, які буде розрізано j -способом, .

Обмеження задачі безпосередньо пов’язані з вимогою забезпечити повністю необхідну кількість нестандартних рулонів за спеціальними замовленнями. Якщо використовуватимуться всі подані в таблиці способи розкрою, дістанемо такі результати.

Кількість рулонів шириною 0,8 м

2 х ₁ + х ₂ + х ₃ = 150.

Кількість рулонів шириною 1 м

х ₃ + 2 х ₄ = 200.

Кількість рулонів шириною 1,2 м

х ₂ + х ₅ = 300.

Цільова функція задачі — це мінімальні загальні втрати паперу під час розрізування стандартних рулонів на рулони нестандар­тної ширини. Математично вона має такий вигляд:

Математична модель задачі в загальному вигляді запису­ється так:

Для розв’язування цієї задачі застосуємо процедуру симплекс-методу. Оскільки задачу сформульовано в канонічній формі, запишемо її відразу у векторній формі:

де

У системі векторів маємо лише один одиничний вектор . Тому в перше та друге обмеження введемо штучні змінні х ₆ та х ₇. Розширена задача матиме вигляд:

 

Процес розв’язування задачі симплекс-методом подано у виг­ляді таблиці:

Базис	С баз	План	0,4		0,2		0,8	М	M	θ
х 1	х 2	х 3	х 4	х 5	х 6	х 7
х 6	М									—
← х 7	M
х 5	0,8									—
Zj – Cj ≥ 0		–0,4	0,8	–0,2
350 M	2 M	M	2 M	2 M
← х 6	M
х 4					1/2					—
х 5	0,8									—
Zj – Cj ≥ 0		–0,4	0,8	–0,2
150 M	2 M	M	M
← х 1	0,4			1/2	1/2
х 4					1/2					—
х 5	0,8
Zj – Cj ≥ 0
х 2
х 4					1/2
х 5	0,8		–2		–1
Zj – Cj ≥ 0		–2		–1

Згідно з останньою симплексною таблицею запишемо оптимальний план задачі:

Х ^* = (0; 150; 0; 100; 150),
min Z = 120.

Визначений оптимальний план передбачає таке: щоб у повному обсязі виконати спеціальні замовлення, які надходять на паперову фабрику, необхідно розрізати 150 стандартних рулонів другим способом, 100 рулонів — четвертим і 150 — п’ятим. За такого оптимального варіанта розкрою величина відходів паперу буде найменшою і становитиме 120 м.

2.2. Приклади та завдання для самостійної роботи

2.1 Фірма має можливість рекламувати свою продукцію, використовуючи для цього телебачення, радіо та газети. Витрати на рекламу в бюджеті фірми обмежені сумою 8000 дол. на місяць. Досвід минулих років показав, що 1 дол., витрачений на телерекламу, дає фірмі прибуток у розмірі 10 дол., а витрачений на рекламу по радіо та в газетах — відповідно 4 та 8 дол.

Фірма має намір витрачати на теле- та радіорекламу не більш як 70 % рекламного бюджету, а витрати на газетну рекламу не повинні більш як удвічі перевищувати витрати на радіорекламу.

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

 

2.2 Розв’язати попередню задачу, якщо вимоги до розподілу рекламного бюджету фірми такі: витрати на радіо­рекламу мають становити не менш як 25 % рекламного бюджету,
а на газетну рекламу — не менш як 50 % витрат на телебачення.

2.3 Промислове підприємство виготовляє продукцію трьох видів А, В і С, для чого використовує два види ресурсів 1 і 2, запаси яких становлять відповідно 4000 та 6000 од. Витрати ресурсів на одиницю продукції кожного виду наведено в таблиці:

Ресурс	Витрати ресурсів на одиницю продукції, ум. од., за видами
А	В	С

Аналіз умов збуту продукції показав, що мінімальний попит на продукцію підприємства для продукції А, В і С відповідно становить 200, 200 та 150 од. Але співвідношення випуску продукції А, В і С має бути 3: 2: 5. Прибуток від реалізації одиниці продукції виду А становить 30 дол., продукції В — 20 дол., а продукції С — 50 дол.

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

 

2.4 Господарство планує вирощувати три сільськогосподарські культури (пшеницю, картоплю та гречку) і може виділити для цього 300 га земельних угідь. Для успішного вирощування сільськогосподарські культури потребують внесення комплексного мінерального добрива, запас якого в господарстві обмежений — 120 т.

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

Показник	Сільськогосподарська культура
Пшениця	Картопля	Гречка
Урожайність, ц/га
Закупівельна ціна, ум. од./ц
Норма внесення добрива, кг/га

Площа земельних угідь, що відводяться під вирощування гречки, має не перевищувати 40 га.

Визначити такий план розподілу посівної площі господарства, який дає найбільший дохід від вирощування сільськогосподарських культур.

 

2.5 Фірма планує організувати виробництво двох видів продукції А та В, але має для цього обмежений інвестиційний фонд у розмірі 5000 дол. У разі потреби цю суму можна збільшити на 10 000 дол. за рахунок банківського кредиту, процентна ставка за використання якого становить 20 %. Втрати, пов’язані з виробництвом одиниці продукції А, дорівнюють 50 дол., а одиниці продукції В — 100 дол.

Очікуваний прибуток фірми від реалізації одиниці продукції А становить 100 дол., а одиниці продукції В — 150 дол.

Фірма має попереднє замовлення на виробництво не менш як 100 одиниць продукції А та 50 одиниць продукції В.

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

 

2.6 Фірма має капітал 300 000 дол., який може використовуватися для фінансування проектів 1 та 2. Реалізація проек­ту 2 гарантує отримання щороку прибутку в розмірі 1 дол. на кожний вкладений долар. Проект 1 гарантує прибуток у розмірі 3 дол. за кожний інвестований долар, але через два роки. У разі фінансування проекту 1 період інвестицій має бути кратним двом рокам.

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

Розв’язати наведені нижче задачі лінійного програмування симплекс-методом.

 

2.7	2.8
2.9	2.10
2.11	2.12
2.13	2.14

 

2.15	2.16
2.17	2.18
2.19	2.20
2.21	2.22
2.23	2.24
2.25	2.26
2.27	2.28
2.29	2.30
2.31	2.32
2.33	2.34