
Lektsiyi_Ukhanska_1
.pdf
двоїстих оцінок не змінювалися. Іншими словами, необхідно знайти такі інтервали зміни
bi |
i 1, m , в межах яких оптимальний план двоїстої задачі залишається незмінним. Це |
||||||||||
має місце для тих значень |
bi |
bi , |
при яких стовпчик вектора P0 останньої симплекс- |
||||||||
таблиці прямої задачі не містить від’ємних чисел, тобто компоненти вектора |
|||||||||||
|
|
|
b |
b |
|
|
|
|
|
|
|
|
|
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
B |
1 |
.......... |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|||||
|
|
|
|
b |
b |
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|||
|
|
|
m |
|
|
|
|
|
|||
повинні бути невід’ємними. |
Тут B |
1 |
- матриця, обернена до матриці B , утвореної з |
||||||||
|
компонент векторів базису, який визначає оптимальний план прямої задачі.
Визначимо інтервали стійкості двоїстих оцінок по відношенню до зміни ресурсів кожного типу для останнього прикладу. Знайдемо компоненти вектора
|
|
b |
b |
|
|
|
1 |
1 |
|
|
|
|
1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
B |
1 |
.......... |
|
|
0.5 |
1 |
|||
|
|
||||||||
|
|
|
b |
b |
|
|
|
0 |
1 |
|
|
|
m |
|
|||||
|
|
m |
|
|
|
0 0 1
|
35 |
b |
||
|
|
|
1 |
|
|
|
|
||
|
|
30 b2 |
||
|
|
40 |
b |
|
|
|
|||
|
3 |
|
5 b |
b |
||
|
1 |
|
2 |
|
|
|
|||
|
12.5 0.5 b |
b |
||
|
1 |
2 |
||
|
10 b |
b |
||
|
||||
2 |
|
3 |
.
( компоненти матриці B
початковий базис прямої
1 |
записані у стовпцях |
P5 , P6 , P7 |
, оскільки ці вектори утворюють |
|
задачі). З умови невід’ємності компонент знаходимо:
5 b |
b |
0 |
|||
|
|
1 |
|
2 |
|
|
0.5b |
b |
|||
12.5 |
|||||
|
|
|
|
1 |
2 |
|
b |
|
b 0 |
||
10 |
2 |
||||
|
|
|
|
3 |
0
.
Нехай
b |
b |
1 |
3 |
0
:
тоді
b |
2 |
||
|
|
||
|
|
||
b2 |
|||
|
b |
|
|
|
2 |
||
|
5 |
|
|
|
|
|
12.5 |
|
12.5 b |
2 |
5 |
|
|
|
|
|
|
|
10 |
|
|
|
|
|
17.5 b2
35
.
Якщо кількість ресурсів виду B належить вказаному інтервалові, а кількість решти видів
ресурсів не змінюється, то двоїста задача має такий самий оптимальний план Y |
* |
3; 4; 0 . |
|
Нехай
b2
b3
0
:
тоді
b |
5 |
|
5 b |
25 |
|
|
|
1 |
|
||||
|
|
|
|
|
||
b |
25 |
|
1 |
|
|
|
|
|
|
|
|
||
1 |
|
|
|
|
|
30 b1
60
.
Нехай b2 b1 0 : |
тоді |
b3 10 |
|
30 b3 |
50 . |
|
|
Остання нерівність |
означає: якщо кількість сировини виду C |
буде збільшено чи |
|||||
зменшено в межах 10 одиниць, оптимальний план двоїстої задачі не зміниться. |
|||||||
Знайдемо тепер як зміниться значення цільової функції при оптимальному плані |
|||||||
прямої задачі, |
якщо |
зменшити кількість сировини |
виду A на 3 |
одиниці, а кількість |
|||
сировини видів |
B, C |
збільшити відповідно на 4 і 5 |
одиниць. Спочатку перевіримо, чи |
можна проводити вказані зміни. Для цього знайдемо
30

b |
|
b |
|
35 3 32, |
1 |
1 |
|
||
b |
2 |
b |
2 |
30 4 34, |
|
|
|
||
b |
|
b |
|
40 5 45, |
3 |
3 |
|
b |
|
b 30; 60 , |
|||
1 |
|
|
1 |
|
17.5;35 , |
b |
2 |
b |
2 |
||
|
|
30;50 , |
|||
b |
|
|
b |
|
|
3 |
|
3 |
|
Отже вказані зміни у використанні сировини можна вносити ( оптимальний план двоїстої
задачі при цьому не зміниться ). Знайдемо приріст функції |
* |
Lmax |
L* max
y |
b y |
b |
y |
b |
|
* |
|
* |
|
* |
|
1 |
1 |
2 |
2 |
3 |
3 |
3 3 4 4 5 0
7
.
Це означає, що значення цільової функції збільшиться на 7 одиниць, тобто можна побудувати такий план виробництва продукції, що прибуток від реалізації буде на 7
одиниць вищим від запланованого при початково заданих кількостях сировини.
Зауважимо, що збільшення кількості сировини виду C не впливає на величину прибутку,
в той час як збільшення кількості сировини виду B на 4 од. веде до збільшення значення
* |
на 16 од., а зменшення сировини виду |
A на 3 од. веде до зменшення |
* |
на 9 од. |
|
Lmax |
Lmax |
||||
|
Якщо значення величин b1 , b2 , b3 |
змінюються одночасно, то дослідження стійкості |
|||
двоїстих оцінок значно ускладнюється. |
|
|
|
|
§ 2.3. Двоїстий симплекс-метод.
Двоїстий сипмлекс-метод використовується для знаходження розв’язку задач лінійного програмування, записаних у формі прямої задачі у випадку, коли вільні члени bi i 1, m в системі обмежень приймають довільні значення ( а не лише невід’ємні ).
Розглянемо задачу лінійного програмування при довільних значеннях bi i 1, m :
n |
|
|
L c j x j max |
||
j 1 |
|
|
n |
|
|
P x |
j |
P |
j |
0 |
|
j 1 |
|
|
(2.7)
(2.8)
x j 0 |
( j 1, n ) |
|
|
|
|
|
|
|
|
(2.9) |
і відповідну їй двоїсту задачу. Нагадаємо, що |
P |
j |
- |
вектори, |
складені з коефіцієнтів при |
|||||
|
|
|||||||||
невідомих у системі обмежень. Припустимо, |
|
що |
серед Pj |
векторів є m одиничних |
||||||
векторів P1 , P2 ,..., Pm , |
|
|
|
|
|
|
|
|
|
|
Pj a1 j , a2 j ,..., amj T , |
j |
|
; |
|
|
P0 b1 , b2 ,..., bm T |
||||
m 1, n |
|
|
Розв’язок системи рівнянь (2.8)
31

X |
0 |
b |
, b |
,..., b |
, 0, 0,...,0 |
|
1 |
2 |
m |
|
не є планом задачі (2.7)-(2.9), тому що серед значень bi i 1, m
Означення. Розв’язок (2.10), що визначається базисом з називається псевдопланом задачі (2.7)-(2.9), якщо
є від’ємні.
m |
одиничних |
(2.10)
векторів,
m |
|
|
|
j ci aij c j |
0 для всіх j 1, n . |
|
|
i 1 |
|
|
|
Теорема 1. Кожному псевдоплану X |
задачі (2.7)-(2.9) з базисом P1 , P2 ,..., Pm |
відповідає |
|
опорний план Y |
відповідної двоїстої задачі з тим же базисом. |
|
|
Якщо деяке значення |
bi 0 , але всі aij 0 , то задача не має розв’язку. Якщо серед |
відповідних значень |
aij |
є від’ємні, то можна перейти до наступного псевдоплану, при |
якому значення цільової функції не зменшується.
Теорема 2. (Ознака оптимальності псевдоплану). Псевдоплан задачі (2.7)-(2.9), всі компоненти якого невід’ємні, є її оптимальним планом.
Отже, пошук розв’язку задач лінійного програмування двоїстим симплекс-методом
полягає у впорядкованому переході від однієї симплекс-таблиці до іншої, доки не буде знайдено оптимальний розв’язок, або не буде доведено, що задача немає розв’язку.
|
Щоб знайти розв’язок задачі лінійного програмування двоїстим симплекс-методом, |
||||
необхідно: |
|
|
|
|
|
1. |
Знайти псевдоплан вихідної задачі. |
|
|
|
|
2. |
Якщо всі bi 0 , то псевдоплан є оптимальним. Якщо серед |
bi є від’ємні, |
то, |
||
|
знайшовши найбільше по модулю від’ємне значення |
bi |
, виводимо вектор з |
||
|
відповідним індексом з базису. |
|
|
|
|
3. |
Визначаємо індекс, при якому модуль відношення j |
до від’ємних aij |
є |
||
|
мінімальним, і відповідний вектор вводимо в базис. |
|
|
|
|
4. |
Процес продовжуємо доти, доки не буде знайдено оптимальний план, тобто доки |
||||
|
всі від’ємні значення bi не будуть виведені з псевдо плану, |
або доки не буде |
|||
|
встановлено, що задача не має розв’язку. |
|
|
|
|
Теорема 3. (Ознака нерозв’язності задачі лінійного програмування на основі псевдо-
плану). Якщо для псевдоплану X з
k k m , що bk 0 і всі коефіцієнти
програмування не має планів.
базисом
a |
kj |
0 |
j |
|
|
|
P1
, P2 ,...Pm
m 1, n
існує такий індекс
, то задача лінійного
32

Приклад.
L x |
x |
2 |
2x |
3 |
max |
1 |
|
|
|
|
x x |
|
|
|
x |
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
1 |
|
|
|
2 |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x1 |
x2 |
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
x |
|
2x |
2 |
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
i |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
xi |
1,3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
Перепишемо систему обмежень |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
x1 x2 |
|
x3 8 |
|
|
|
|
|
|
|
x1 x2 x3 8 |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
x1 x2 |
|
|
|
|
|
|
x1 x2 x4 4 |
|
|
|
|||||||||||||||||||||||||
|
|
|
x1 2x2 6 |
|
|
|
|
|
x1 2x2 x5 6 |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
x |
i |
|
0 |
i 1,3 |
|
|
|
|
|
|
|
x |
i |
0 |
i 1,5 |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Побудуємо двоїсту задачу |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
* |
8y |
|
|
4y |
|
|
6y |
|
min |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
L |
|
|
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
y |
y |
2 |
y |
3 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
y y |
|
|
|
2 y |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
2 |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
y |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
y |
i |
0, |
|
i 1,3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Результати розв’язування прямої задачі запишемо у симплекс-таблиці. |
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
№ |
|
|
|
|
|
|
|
|
Б |
|
|
|
|
Сб |
|
Р0 |
1 |
|
|
|
|
1 |
|
|
2 |
|
0 |
0 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Р1 |
|
|
|
|
Р2 |
|
Р3 |
|
Р4 |
Р5 |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
|
|
|
|
|
|
|
|
Р3 |
|
|
|
|
2 |
|
|
8 |
1 |
|
|
|
|
1 |
|
|
1 |
|
0 |
0 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2 |
|
|
|
|
|
|
|
|
|
Р4 |
|
|
|
|
0 |
|
|
-4 |
-1 |
|
|
|
|
1 |
|
|
0 |
|
1 |
0 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
3 |
|
|
|
|
|
|
|
|
|
Р5 |
|
|
|
|
0 |
|
|
-6 |
-1 |
|
|
|
|
|
-2 |
|
|
0 |
|
0 |
1 |
|
|||
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L0=16 |
1 |
|
|
|
|
1 |
|
|
0 |
|
0 |
0 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
|
|
|
|
|
|
|
|
Р3 |
|
|
|
|
2 |
|
|
5 |
1/2 |
|
|
|
|
0 |
|
|
1 |
|
0 |
½ |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
2 |
|
|
|
|
|
|
|
|
|
Р4 |
|
|
|
|
0 |
|
|
-7 |
|
-3/2 |
|
|
|
|
0 |
|
|
0 |
|
1 |
½ |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
3 |
|
|
|
|
|
|
|
|
|
Р2 |
|
|
|
|
1 |
|
|
3 |
1/2 |
|
|
|
|
1 |
|
|
0 |
|
0 |
-1/2 |
|
||||
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
13 |
1/2 |
|
|
|
|
0 |
|
|
0 |
|
0 |
½ |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
|
|
|
|
|
|
|
|
|
Р3 |
|
|
|
|
2 |
|
|
8/3 |
0 |
|
|
|
|
0 |
|
|
1 |
|
1/3 |
2/3 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
2 |
|
|
|
|
|
|
|
|
|
Р1 |
|
|
|
|
1 |
|
|
14/3 |
1 |
|
|
|
|
0 |
|
|
0 |
|
-2/3 |
-1/3 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
3 |
|
|
|
|
|
|
|
|
|
Р2 |
|
|
|
|
1 |
|
|
2/3 |
0 |
|
|
|
|
1 |
|
|
0 |
|
1/3 |
-1/3 |
|
||||
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
32/3 |
0 |
|
|
|
|
0 |
|
|
0 |
|
1/3 |
2/3 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
33

|
В |
першій симплекс-таблиці всі значення j 0 і початковий псевдоплан має |
вигляд |
X |
0 0, 0,8, 4, 6 і значення L0 16. Від’ємних елементів тут є два ( -4 і –6 ) і |
більшим по модулю є другий, то розглядаємо елементи 3-го рядка. Серед них є від’ємні.
Це означає, що з базису виводимо вектор
P5
. Знайшовши
|
1 |
, |
1 |
|
|
min |
1 |
2 |
|
||
|
|
|
|
1 2
, робимо
висновок, що в базис вводиться вектор
P2
і ключовий елемент дорівнює -2. Провівши
далі обчислення, з
плані |
X |
* |
14 |
; |
2 |
; |
|
|
|
|
|
||||
|
|
|
|
3 |
|
3 |
|
останньої таблиці маємо, що
8 . 3
L |
|
* |
10 |
2 |
; |
|
max |
L |
min |
3 |
|||
|
|
|
|
|||
|
|
|
|
|
|
при оптимальному
Оптимальний план двоїстої задачі матимемо, розглянувши оцінки опорного плану
вихідної задачі, які відповідають початковому базисові |
P3 , P4 , P5 |
. З останньої симплекс- |
||||||||||||||||||||||||||||
таблиці |
з |
четвертого |
рядка |
і |
стовпців |
векторів |
P3 , P4 |
, P5 |
знаходимо |
компоненти |
||||||||||||||||||||
* |
|
|
* |
|
1 |
|
0 |
|
1 |
|
|
* |
|
2 |
|
0 |
|
2 |
|
|
|
|
* |
|
|
1 |
|
2 |
|
|
2, y |
|
|
, |
|
|
або |
Y |
2; |
|
; |
|
. Оскільки всі |
||||||||||||||||||
y1 0 2 |
2 |
|
3 |
|
3 |
y3 |
|
3 |
|
3 |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
3 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
значення |
yi |
0 |
, |
то всі види сировини використовуються повністю. |
Якщо, |
наприклад, |
кількість першого виду сировини збільшити на 1 одиницю, то можна скласти такий опорний план, що значення цільової функції збільшиться на 2 одиниці.
34

РОЗДІЛ 3. ТРАНСПОРТНА ЗАДАЧА.
§ 3.1. Постановка транспортної задачі.
Загальні методи лінійного програмування ( симплексний, двоїстий симплексний ) в
принципі дають можливість розв’язати будь-яку задачу, проте, як правило цей розв’язок у практично важливих випадках пов’язаний зі значними обчисленнями. Тому виникає питання про виділення окремих класів задач, розв’язок яких можна отримати за допомогою пристосованих для них більш простих спеціальних методів обчислень.
Найбільш широким класом таких задач є так звані транспортні задачі.
Класична транспортна задача – це задача про найбільш економний план перевезення однорідного продукту чи взаємозамінних продуктів із пунктів виробництва в пункти споживання. Математична структура транспортної задачі характерна для великого класу задач лінійного програмування, реальний зміст яких може бути найрізноманітнішим, зовсім не пов’язаним із перевезенням вантажів. Транспортна модель широко використовується для розв’язування задач розміщення виробництва, розподілу
капіталовкладень, задачі оптимального призначення тощо.
Нехай в m пунктах виробництва |
A1 |
, A2 |
,..., Am |
міститься однорідний товар, який |
треба перевезти в n пунктів призначення |
B1 |
, B2 |
,..., Bn |
. Відома кількість одиниць товару в |
кожному пункті постачання і скільки одиниць товару потребує кожний пункт призначення. Крім того, задана собівартість перевезень одиниці товару з кожного пункту постачання у кожний пункт призначення. Необхідно знайти найекономніший план перевезення повного обсягу товарів із пунктів виробництва в пункти призначення ( в
якості критерію оптимальності беруть як правило мінімальні затрати на перевезення всього вантажу або мінімальний час доставки вантажу).
Математична модель транспортної задачі.
Нехай
m n
-кількість пунктів постачання;
-кількість пунктів призначення;
ai
b |
j |
|
cij
- кількість одиниць товару в i -му пункті постачання i 1, m ;
-кількість одиниць товару, що потребує j -тий пункт призначення j
-собівартість перевезення одиниці товару з i -го пункту постачання
пункт призначення B j ;
1, n
Ai в
; j
-тий
35

xij |
- кількість одиниць товару, яку планується перевезти з |
i -го пункту постачання в |
|
j -тий пункт призначення (невідомі величини). |
|
З економічної точки зору задача полягає в тому, що необхідно так планувати
перевезення, щоб їх загальна вартість була мінімальною. Запишемо математичну модель:
m |
n |
|
L cij xij min |
(3.1) |
|
i 1 |
j 1 |
|
n |
|
i 1, m |
xij |
ai |
|
j 1 |
|
|
(3.2)
m |
|
j |
|
|
|
||||
xij |
b j |
1, n |
(3.3) |
||||||
i 1 |
|
|
|
|
|
|
|
|
|
xij 0 i |
|
, j |
|
|
|
||||
1.m |
1, n |
(3.4) |
|||||||
Зауважимо, |
що |
у випадку, коли |
загальна кількість товару, який є у пунктах |
постачання, дорівнює кількості одиниць товару, що потребують пункти призначення, то виконується так звана умова балансу
m |
n |
ai |
b j |
i 1 |
j 1 |
(3.5)
Всі дані транспортної задачі записують у вигляді таблиці
Пункти при- |
|
|
|
|
Пункти значення |
В1 |
... |
Вj |
... |
постачання |
|
|
|
|
|
c |
|
|
|
|
c |
|
|
|
|
А1 |
11 |
|
|
... |
1 j |
|
|
... |
||
|
x |
|
|
|
x |
|||||
|
|
|
11 |
|
|
|
|
|||
|
|
|
|
|
|
|
1 j |
|
||
|
|
|
|
|
|
|
|
|
|
|
... |
|
... |
|
|
... |
|
... |
|
... |
|
|
c |
i1 |
|
|
|
c |
ij |
|
|
|
Аi |
|
|
|
... |
|
|
|
... |
||
|
x |
|
|
|
x |
|
||||
|
|
i1 |
|
|
|
ij |
|
|||
|
|
|
|
|
|
|
|
|
||
... |
|
... |
|
|
... |
|
... |
|
... |
|
Аm |
cm1 |
|
|
... |
cmj |
|
|
... |
||
|
x |
|
|
|
|
xmj |
||||
|
|
m1 |
|
|
|
|
||||
|
|
|
|
|
|
|
||||
Потреби |
|
b |
|
|
... |
|
b |
j |
|
... |
|
1 |
|
|
|
|
|
Вn
c1n
x1n
...
c |
in |
|
|
x |
in |
||
|
|
|
|
|
|
... |
|
||
c |
mn |
|
|
|
|
|
|
|
|
|
x |
mn |
||
|
|
|
||
|
b |
n |
|
|
|
|
|
Запаси
a |
|
|
|
1 |
|
... |
||
a |
i |
|
|
|
|
... |
||
a |
m |
|
|
Означення 1. Планом транспортної задачі називається набір величин xij i 1, m, j 1, n , який задовольняє умови (3.2)-(3.4).
Позначимо план транспортної задачі через X |
xij |
, а величини xij |
називатимемо перевезеннями. |
m,n |
|
|
36

Властивості транспортної задачі.
1.Транспортна задача (3.1)-(3.4) завжди має розв’язок, якщо виконується балансова умова (3.5).
2. Кількість лінійно-незалежних рівнянь системи (3.2)-(3.3) дорівнює
m
n
1
.
Позначимо
|
|
|
m j |
|
|
T |
|
|
|
|
|
||
P |
|
|
0,0,...0,1,0,...0,1, 0,...,0 |
|
|
|
ij |
|
|
|
|
|
|
|
|
|
||||
|
|
|
i |
n j |
|
|
- вектори розміру
m
n
, компонентами
яких є коефіцієнти при невідомих |
xij |
||||||
вигляді |
|
|
|
|
|
|
|
m |
n |
|
|
|
|
|
|
|
x |
ij |
P |
P |
, |
||
|
ij |
|
0 |
|
|||
i 1 |
j 1 |
|
|
|
|
|
|
системи
P0
(3.2)-(3.3).
a |
, a |
2 |
,..., a |
1 |
|
|
Перепишемо систему (3.2)-(3.3) у
m ,b1 ,...bn T .
Означення 2. План X xij |
транспортної задачі називається опорним, якщо вектори |
|||||
|
m,n |
|
|
|
|
|
Pij , які відповідають додатним компонентам |
xij |
плану |
X , утворюють |
|||
лінійно незалежну систему. |
|
|
|
|||
Означення 3. Опорний план |
X xij |
m,n |
транспортної задачі називається невиродженим, |
|||
|
|
|
|
|
|
|
якщо він містить m n 1 |
додатних величин |
xij . |
Якщо цей план містить |
менше, ніж m n 1 додатних величин xij , то він називається виродженим.
Означення 4. Транспортна задача називається невиродженою, якщо всі її опорні плани невироджені. Якщо хоча б один опорний план транспортної задачі є виродженим, то така задача називається виродженою.
Означення 5. Якщо виконується балансова умова (3.5), то така задача називається
збалансованою ( або транспортною задачею з правильним балансом ), а
відповідна математична модель називається закритою. Якщо умова (3.5) не виконується, то транспортна задача називається незбалансованою ( за-
дачею з неправильним балансом ). Останню можна звести до задачі з правильним балансом.
§ 3.2. Методи пошуку початкових опорних планів транспортної задачі.
Як і для інших задач лінійного програмування, ітераційний процес по пошуку оптимального плану транспортної задачі починається з деякого опорного плану. Найбільш поширеними методами пошуку початкового опорного плану транспортної задачі є метод північно-західного кута, метод мінімального елемента і метод апроксимації Фогеля. Суть всіх цих методів полягає у послідовному заповненні на протязі m n 1 кроків клітинок у
37

таблиці умов транспортної задачі. Кожна заповнена клітинка називається зайнятою або
базовою, а незаповнена – вільною. Якщо в даному розв’язку змінна xij дорівнює деякому
числу |
a , то це число записуємо у відповідній клітинці |
(i, j) ; якщо |
xij |
0 |
, то таку |
клітинку залишаємо вільною. Таким чином, кожному опорному розв’язкові буде відповідати m n 1 базових клітинок і mn m n 1 вільних клітинок.
1. Метод північно-західного кута ( діагональний метод ).
Заповнення таблиці починають із клітинки (1;1). На кожному кроці додатними
перевезеннями заповнюють тільки верхню ліву клітинку незаповненої частини таблиці,
поміщаючи в цю клітинку максимально можливу кількість продукту. На першому кроці
|
min a |
|
. Якщо a |
b |
|
|
a , |
|
0 j |
|
|
( матимемо нульові |
|
визначаємо x |
,b |
, то |
x |
x |
2, n |
||||||||
11 |
1 |
1 |
1 |
1 |
|
11 |
1 |
1 j |
|
|
|
|
|
значення в першому рядку). Якщо a1 |
b1 , то |
|
b1 , |
|
0 i |
|
( матимемо нульові |
||||||
x11 |
xi1 |
2, m |
|||||||||||
значення в першому стовпчику). Якщо a1 |
b1 |
, то нульовими перевезеннями заповнюють і |
перший рядок і перший стовпчик. На першому кроці додатними перевезеннями заповнюється лише одна клітинка - (1;1).
|
Якщо на першому кроці заповнено 1-й рядок, то на другому кроці визначаємо |
|||
x21 |
min a2 |
,b1 x11 |
. |
Якщо на першому кроці заповнено 1-й стовпчик, то знаходимо |
x12 |
min a1 |
x11 ,b2 |
. |
Якщо заповнено і 1-й стовпчик і 1-й рядок, то знаходимо |
x22 |
min a2 ,b2 . Залежно від знайденої величини нульовими елементами заповнюють або |
другий рядок, або другий стовпчик, або кожний з них. Аналогічно заповнюємо всю таблицю.
Нехай |
N - кількість додатних перевезень плану, побудованого методом північно- |
західного кута. Якщо N m n 1 |
, то опорний план є невиродженим; якщо |
N m n 1 |
, |
то опорний план є виродженим. |
|
|
|
Таким чином, якщо ми побудуємо невироджений опорний план транспортної задачі, то
в матриці X xij |
додатними будуть тільки |
m n 1 |
елементів, а решта будуть |
|
m,n |
|
|
дорівнювати нулю.
Побудований таким чином опорний план буде далеким від оптимального, оскільки при його побудові не враховувалась вартість перевезень cij .
Приклад. Знайти початковий опорний план транспортної задачі, що характеризується матрицею перевезень, заданою таблицею.
38

|
|
Пункти призна- |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
Пункти |
чення |
|
|
|
В1 |
|
|
В2 |
|
|
В3 |
|
|
|
В4 |
|
В5 |
|
|
Запаси |
|
|
|
|||
|
|
постачання |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А1 |
|
|
|
3 |
|
|
1 |
|
|
3 |
|
|
|
4 |
|
3 |
|
|
|
|
65 |
|
|
|
|
|
|
|
|
|
|
|
27 |
|
|
38 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А2 |
|
|
|
2 |
|
|
3 |
|
|
1 |
|
|
|
2 |
|
3 |
|
|
|
|
68 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
15 |
|
21 |
|
|
32 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
А3 |
|
|
|
3 |
|
|
5 |
|
|
2 |
|
|
|
2 |
|
4 |
|
|
|
|
40 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
30 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
Потреби |
|
|
|
|
27 |
|
|
53 |
|
|
21 |
|
|
|
42 |
|
30 |
|
|
173 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
Запаси пункту A1 |
становлять 65 одиниць, а потреби пункту |
|
B1 |
- 27 одиниць, тому |
||||||||||||||||||||||
можна завезти з пункту |
A1 в пункт |
B1 27 одиниць продукції |
x11 |
min 27;65 27 . |
||||||||||||||||||||||||
Потреби |
B1 |
задовольняться повністю, а запаси |
A1 |
зменшаться до 65-27=38 од., тобто |
||||||||||||||||||||||||
тепер запаси пункту A1 |
|
становлять 38 од. Потреби пункту призначення B2 - 53 од. Після |
||||||||||||||||||||||||||
того як у пункт B2 |
завеземо 38 од. товару x12 |
min 65 27;53 min 38;53 38 , |
ми |
|||||||||||||||||||||||||
вичерпаємо запаси пункту |
A1 і одночасно зменшимо потреби пункту B2 до 15 од., |
які |
||||||||||||||||||||||||||
можна задовольнити запасами пункту |
A2 |
x22 |
min 68;53 38 min 68;15 15 . Тоді |
|||||||||||||||||||||||||
запаси товару в пункті |
A2 зменшаться до 68-15=53 од. Потреби пункту |
B3 - 21 од. Завізши |
||||||||||||||||||||||||||
21 |
|
од. |
продукції |
з |
пункту |
|
A2 |
в |
пункт |
B3 |
, |
його |
потреби |
|
будуть |
задоволені |
||||||||||||
x23 |
min 68 15; 21 min 43; 21 21, |
а в пункті |
A2 залишиться 53-21=32 |
од. товару. |
||||||||||||||||||||||||
Пункт призначення |
B4 |
|
потребує 42 од. продукції, |
а задовольнити їх можемо 32 од. з |
||||||||||||||||||||||||
пункту постачання |
A2 |
x24 |
min 68 21 15; 42 min 32; 42 32 |
і 10 од. з пункту |
A3 |
|||||||||||||||||||||||
x34 |
min 40; 42 32 min 40;10 10 . Тоді |
в |
пункті |
A3 |
залишиться 40-10=30 |
од. |
||||||||||||||||||||||
продукції, яку завеземо в пункт B5 |
x35 |
min 40 10;30 min 30;30 30 .. |
|
|
|
|||||||||||||||||||||||
|
|
Число |
базових |
клітинок дорівнює |
m n 1 3 5 1 7 |
, |
що узгоджується з |
|||||||||||||||||||||
умовою оцінки опорності плану. Отже, опорний план має вигляд |
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
27 |
38 |
0 |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X |
0 |
15 |
21 |
32 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
0 |
0 |
0 |
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
30 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
і значення лінійної функції
L 3 27 1 38 3 15 1 21 2 32 4 30 2 10 389.
39