
книги из ГПНТБ / Косенко Б.Ф. Многоэтапная транспортная задача монография
.pdfТ а б л и ц а 3
Вершины
La
La
Ь2
ТА
1 А
l a
*1 *2
4 |
0 |
8 |
4 |
2 |
14 |
2 |
3 |
13 |
0 |
3 |
3 |
0 |
0 |
0 |
Вершины Хь
L B |
3 |
/? |
3 |
Ч |
2 |
•"4 |
|
00 |
|
l a |
0 |
0 |
0 |
|
Вершины |
х» |
X. |
LB |
г
LB |
3 |
0 |
6 |
LB |
3 |
4 |
18 |
Ь2 |
|
|
|
L B |
0 |
4 |
12 |
LB |
0 |
0 |
0 |
•
|
Т а б л и ц а 4 |
ib |
l ab |
|
|
0 |
6 |
2,5 |
13,5 |
3 |
13 |
3 |
9 |
0 |
0 |
|
• |
Получили две области определения функций LA и LB (рис. 5)г соответственно многоугольники Мj и М2.
Рис. 5. Сопряжение в жесткой цепи.
20
|
|
|
|
|
|
|
|
в в |
|
Области определения соприкасаются по прямой L4L5 на |
|||||||||
ютрезке L^Lz. Интуитивно можно утверждать, |
что |
оптимальное |
|||||||
решение лежит в вершине пятиугольника (области Мх), |
так |
как |
|||||||
L2 принадлежит М2 и M3f a Lz—M2 и Mv |
в то время |
как |
ни |
||||||
одна из вершин |
пятиугольника |
(области |
М2) |
не |
принадлежит |
||||
пятиугольнику |
области Мх. |
Посмотрим, так ли это (табл. 5). |
|||||||
|
|
|
|
|
|
|
Т а б л и ц а |
5 |
|
Вершины |
*1 |
*2 |
1-Л |
, Вершины |
Хз |
*4 |
^в |
||
La |
4 |
0 |
8 |
LB |
|
10 |
0 |
20 |
|
4 |
4 |
2 |
14 |
LB |
|
6,5 |
2,8 |
|
1.4 |
Ь2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
1 А |
2 |
3 |
13 |
LB |
|
4 |
4 |
20 |
|
|
0 |
3 |
3 |
LB |
|
0 |
4 |
12 |
|
|
0 |
0 |
0 |
Т В |
|
8 |
0 |
16 |
|
|
Lb ' |
|
Совместное оптимальное решение в данном случае можно
чюлучить только в точке Таким образом, если области определения имеют общий от
резок прямой, то совмещенное оптимальное решение может быть на одном или на обоих концах отрезка. В этом случае
.для получения совмещенного оптимального решения нужно
отойти |
от оптимального плана в одном |
из звеньев, |
именно |
в том, |
вершина которого (определяющая |
оптимальное |
решение) |
.не принадлежит общему отрезку, или в обоих звеньях, если ни
одна из вершин, определяющих оптимальное решение, |
не при |
надлежит общему отрезку областей определения |
линейных |
форм. |
|
Пример 4. В условиях предыдущих примеров изменим вторую
группу неравенств таким образом, |
чтобы области |
определения |
•функций La и LB имели бы только |
одну общую точку (рис. 6). |
|
12<2Аг3 + 2Аг4< 18; |
(39) |
|
2б4< 4 . |
|
|
Две области определения функций LA и LB имеют лишь одну общую точку, в которой можно ожидать сопряженное опти мальное решение (табл. 6).
Анализ данных (табл. 6) позволяет сделать вывод о том, что ;хотя оптимальные планы ЬА и LB находятся соответственно
АВ
ввершинах 12 и Z,2, они не удовлетворяют условию (32)—(34).
Единственная точка L2, удовлетворяющая этому условию принадлежит одновременно пятиугольнику М х и четырехуголь нику М2 (рис. 6).
Таким образом, при нахождении совмещенного оптимального*
решения двух звеньев необходимо, по крайней |
мере,* |
в |
одном* |
|
из них отойти от оптимального решения, |
если |
вершина |
много |
|
угольника ограничений, определяющая |
оптимальное |
решение,. |
Рис. 6. Сопряжение в жесткой цепи.
принадлежит обоим многоугольникам. Если же общая точка* обеих областей существования функции не определяет оптималь ности одного из звеньев, то в этом случае для нахождения! совмещенного оптимального плана необходимо отойти от опти мального решения в обоих звеньях.
Т а б л и ц а б
|
|
т |
11 |
|
|
|
|
Вершины |
X, |
! |
Вершины |
|
X* |
^ В |
|
-л |
j |
* 3 |
'
тА |
4 |
0 |
8 |
L B |
10 |
0 |
20 |
тА |
Л |
г» |
14 |
i f |
6 |
4 |
24 |
Ь 2 |
т |
Z |
|||||
ТА |
2 |
3 |
13 |
|
2 |
4 |
16 |
L4 |
1 |
3 |
i f |
6 |
0 |
12 |
|
0 |
9 |
||||||
La |
|
|
|
|
|
|
|
La |
0 |
0 |
0 |
|
|
|
|
|
|
|
|
|
• |
|
Пример 5. Сохранив систему неравенств (34), видоизменим неравенства (35) с таким расчетом, чтобы области определения^ функций La и Lb не имели общих точек (рис. 7).
22
|
|
|
|
|
8 > X i > 4; |
|
|
||
|
|
|
|
|
0 -'С Х 3 |
2. |
|
|
|
Определим оптимальное решение (табл. 7). |
|||||||||
Данные табл. 7 свидетель |
7С2 ! ОС(f |
|
|||||||
ствуют |
о том, |
что в этом слу |
|
||||||
чае совместное решение во |
|
|
|
||||||
обще невозможно, так как не |
|
|
|
||||||
существует таких целочислен |
|
|
|
||||||
ных значений Х ъ и |
Х4, кото |
|
|
|
|||||
рые удовлетворяли бы выпол |
|
|
|
||||||
нению требования (32) и (33). |
|
|
|
||||||
В |
общем |
случае |
имеется |
|
|
|
|||
три основных варианта: |
|
|
|
||||||
— Системы линейных урав |
|
|
|
||||||
нений, |
определяющих область |
|
|
|
|||||
существования |
решения по |
|
|
|
|||||
каждому звену, тождественны. |
|
|
|
||||||
В этом случае |
имеем |
гибкий |
|
|
|
||||
двузвенник (рис. 8,а), в кото |
|
|
|
||||||
ром |
оптимальное |
решение |
|
|
|
||||
в одном звене точно соответ |
** 1 ~~П " |
I |
■■т |
||||||
ствует |
оптимальному |
реше |
|||||||
LAS 1 |
г |
з |
|||||||
нию во |
втором и оба |
сопря |
|||||||
жены. |
|
Тогда |
М х= М2= Мг, |
|
|
|
(4 0 )
-------- г-
4 5
Рис. 7. Сопряжение в предельно жесткой цепи.
— Системы линейных уравнений высекают для каждого звена области определения (Мг и М2) линейной формы, которые не тождественны. В этом случае оптимальное решение по цепи
|
|
|
|
|
|
Т а б л и ц а 7 |
|
Вершины |
|
*2 |
La |
Вершины |
*8 |
* 4 |
L B |
1 А |
4 |
0 |
8 |
LB |
2 |
4 |
16 |
La |
4 |
2 |
14 |
LB |
2 |
8 |
28 |
l 2 |
|
|
|
L2 |
|
|
|
l a |
2 |
3 |
13 |
J B |
0 |
8 |
24 |
La |
0 |
3 |
3 |
LB |
0 |
4 |
12 |
l a |
0 |
0 |
0 |
|
|
|
|
Lb |
|
|
|
|
|
|
|
* |
23 |
находится в вершине (рис. 8,6) многоугольника (ЛТ), образован-
„л- |
° |
/И 1 |
или |
ного пересечением многоугольников М• х и М2, a |
at = |
||
М3 . |
|
|
усло |
а2==ЛГ a i _ П° отношению к условиям входа и а2 — по |
виям выхода). Так как Ж3< М 1 и МН< М 2, а М3 ф 0, то 0<а<Т,
что соответствует полу гибкому двузвеннику.
# Х2 -Х* |
'в; |
• |
Рис. 8. Геометрический смысл коэффициента гибкости:
а—гибкий двузвенник; б—полугибкий двузвенник; в—жесткий двузвенник.
—Системы линейных уравнений высекают в каждом звене
области определения (рис. 8,в), |
которые |
не пересекаются, |
а только соприкасаются. Тогда М 3 = 0; откуда |
а —0, а значе |
|
ние а = 0 соответствует жесткому |
двузвеннику. |
24
Таким образом, геометрический анализ позволяет сделать следующие выводы:
— максимум (минимум) линейной формы совмещенного опти мального плана, как правило, меньше (больше) максимума (ми нимума) хотя бы одной из линейных форм звеньев, составляю щих цепь; это объясняется тем, что для получения сопряжен ного оптимального распределения необходимо в плане, по крайней мере одного из звеньев, отойти от экстремального решения;
— сопряженное оптимальное |
решение |
лежит |
в |
общей для |
|||||
обоих звеньев части многоугольника (на |
прямой или в вершине); |
||||||||
— не все точки общего многоугольника |
определения сопря |
||||||||
женного решения дают целочисленное решение. |
|
|
|
||||||
В общей задаче линейного многоэтапного |
программирования |
||||||||
имеем дело не с одним, |
а с г -|- 1 —выпуклым |
многогранником |
|||||||
при г + 1 звене, которые, |
являясь звеньями |
одной цепи, после |
|||||||
довательно попарно сопряжены. |
Каждый многогранник A, F, 'F, |
||||||||
. .., |
В соответственно высекается в т, |
ku |
k2, ... |
, |
п- мерных |
||||
пространствах линейными ограничениями (29) и (30). |
|
||||||||
Рассмотрим |
множество А (звено AQF) и допустим, что оно |
||||||||
ограничено. Тогда А является выпуклым |
многогранником т- |
||||||||
мерного пространства, поскольку имеется т линейно |
независи |
||||||||
мых |
векторов. |
Размерность |
многогранника |
не |
превосходит |
||||
т—fc, так как |
он принадлежит |
общей части гиперплоскостей, |
|||||||
отвечающих линейно независимым уравнениям (29) и (30). |
|||||||||
Линейная форма задачи |
определяет в /72-мерном пространстве |
||||||||
семейство параллельных гиперплоскостей |
|
|
|
|
|
С1Х\ + С2Х2+. . . + скхк = L ( —оо < L < оо),
а |
коэффициенты линейной формы задают |
вектор |
C = (clt c2t |
||
. . |
ск), указывающий направление возрастания L. |
При |
этом |
||
вектор С перпендикулярен семейству гиперплоскостей. |
в |
сто |
|||
|
Передвигая гиперплоскость параллельно |
самой себе |
|||
рону возрастания (уменьшения) Z., можно прийти к такому |
|||||
положению, когда при дальнейшем перемещении она |
не |
будет |
иметь общих точек с многогранником Л. Следовательно, экстре
мальное |
решение лежит |
в |
общих точках многогранника А |
и гиперплоскости, в случае |
расположения многогранника по одну |
||
сторону от последней. Если |
задача имеет единственное реше |
||
ние, то |
многогранник А |
и |
предельная гиперплоскость имеют |
лишь одну общую точку, совпадающую с вершиной много гранника.
При наличии г -f 1 звена в цепи оптимальное решение двух сопряженных звеньев будет находиться в области, образованной пересечением соответствующих многогранников. Размеры этой области будут определяться гибкостью цепи.
25
Если звенья гибкие, то многогранник звена А , образованный ограничениями (29), и многогранник звена В, образованный ограничениями (30), в силу того, что
|
|
|
|
/=1 |
|
|
|
|
|
|
|
|
|
|
|
|
т. е. |
а > 1 , совпадут и в этом случае |
область |
их |
пересечения |
||||||||||||
F совпадет с каждым из этих многогранников, |
откуда |
следует, |
||||||||||||||
что F принадлежит А и В. |
|
|
|
|
|
|
|
|
|
соответствую |
||||||
Для полугибких звеньев область пересечения |
||||||||||||||||
щих им многогранников А и В, |
вследствие |
того, |
что |
|
|
|
||||||||||
|
|
|
т |
|
|
к |
|
|
|
|
|
|
|
|
|
|
|
f e < 2 а , < У f e , |
|
|
|
|
|
|
|
|
|||||||
|
|
|
i= 1 |
|
е= \ |
|
|
|
|
|
|
|
|
|
||
будет представлять собой ^-мерный |
многогранник |
Z7, |
меньший |
|||||||||||||
каждого из них, т. е. F |
принадлежит А и В, а общая |
вершина |
||||||||||||||
/?, определяющая сопряженное оптимальное решение, |
принад |
|||||||||||||||
лежит Л, В и F. |
|
|
|
|
|
|
|
|
|
|
|
решение |
по |
|||
В жестком звене совмещенное оптимальное |
|
|||||||||||||||
цепи лежит в вершинах |
многоугольника, |
образованного каса |
||||||||||||||
нием многогранников А и В, |
при этом по условию |
|
|
|
|
|||||||||||
|
|
к |
|
|
т |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
f e |
= |
I I |
a i |
|
|
|
|
|
|
|
|
|
|
|
|
e= l |
|
2 = 1 |
|
|
|
|
|
|
|
|
|
|
|
|
R принадлежит А и В, а так |
как |
а = |
0, |
то |
решение |
по обоим |
||||||||||
звеньям (А и В) однозначно и взаимное |
влияние не имеет места, |
|||||||||||||||
т. е. сумма линейных форм |
звеньев |
равна |
линейной |
форме |
||||||||||||
цепи. ' |
двузвеннике |
|
совмещенное |
оптималь |
||||||||||||
В предельно жестком |
|
|||||||||||||||
ное |
решение может быть, |
в |
зависимости |
от |
величины а, |
или |
||||||||||
в вершинах линии касания, |
или в точке |
касания, |
так как |
|
||||||||||||
|
|
|
к |
|
т |
|
|
|
|
|
|
|
|
|
|
|
|
о < 2 Л < 2 « / - |
|
|
|
|
|
|
|
|
|||||||
|
|
|
е—\ |
|
/= 1 |
|
|
|
|
|
|
|
|
|
||
В абсолютно жесткой цепи |
|
(а = |
0) |
совмещенного |
оптималь |
ного решения нет, вследствие того, что у множеств А и В нет общих точек, так как
2 Л = о.
е=1
26 |
♦ |
ОСНОВНЫЕ ЗАВИСИМОСТИ МНОГОЭТАПНОЙ ТРАНСПОРТНОЙ
ЗАДАЧИ
Из постановки и геометрической интерпретации задачи много этапного линейного программирования нетрудно заметить, что оно во многом опирается на линейное программирование, его теоремы и утверждения, но, имея более широкую постановку, требует доказательства ряда положений.
Теорема 1. Основное свойство многоэтапного программирова ния заключается в том, что совпадение в транзитном пункте коммуникаций с экстремальными показателями линейной формы, принадлежащих двум смежным звеньям, возможно, но маловеро ятно.
Действительно, в звене AOF выбор входных коммуникаций определяется условиями экстремальности линейной формы, ко торая, в свою очередь, является функцией показателя линейной формы Cie, но число возможных коммуникаций равно тк.
этого числа избираются только те, которые удовлетворяют усло виям экстремальности линейной формы. Практически из тк ком муникаций выбирается не более т -Т k — 1 (условие опорности),. так как даже при наличии равных Cie случайным образом выби
рается |
одна. |
о выходных |
коммуникациях |
То же самое можно сказать и |
|||
(с |
а если это так, то можно |
утверждать, |
что совпадение |
входных и выходных коммуникаций в избранном транзитном пункте определяется вероятностью
1
P(mkn) = mkrti ’
а вероятность совпадения коммуникаций в одном из транзитных пунктов будет равна
Р { т к п ) = - ш - |
<42> |
Из этого основного свойства следует:
- При достаточно большом числе транзитных пунктов мало вероятно совпадение оптимального распределения по цепи с рас пределениями, полученными независимо по звеньям. Следователь но, для сопряжения звеньев требуется балансировка, а это неиз бежно приводит к необходимости отхода от оптимального рас пределения, хотя бы в одном из звеньев цепи транспортировки.
— Сумма максимумов (минимумов) линейных форм, определен ных для каждого последующего звена цепи с учетом получен ного распределения в предшествующем звене, меньше (больше} или равна сумме максимумов (минимумов) линейных форм, опре деленных для этих же звеньев независимо, т. е.
2J
max L ( A Q F Q B ) < max L (A © F) + max L( BQ F);
min L ( A Q F Q B ) > min L (A © F) + min A (BQF).
— Сумма максимумов (минимумов) линейных форм,определеншых последовательно для каждого звена цепи, меньше (больше) линейной формы для цепи в целом
min L ( A Q F , F Q В) > min L (A Q F Q В).
Теорема 2. Если в гибкой цепи (а > |
1) из г + |
1 звена в ком |
|
муникациях x ie.. ./ избраны |
экстремальные значения показате |
||
лей линейной формы Cie...j |
и задача решена относительно их, |
||
то план X = (л'1 ... 1, • • *, Xie... j) является оптимальным. |
|||
Доказательство. Для |
определения |
набора |
экстремальных |
-Cie... у можно воспользоваться методом потенциалов [7] или ал горитмом теории графов [9], с помощью которых можно уста
новить кратчайшие расстояния в транспортной |
сети, имеющей |
||
k. klf • • |
• kr вариантов коммуникаций |
для одного входа и вы |
|
хода. |
|
|
|
После определения экстремальных значений находим распре |
|||
деление, |
при котором линейная форма |
получит |
экстремальное |
значение: |
|
|
|
L (А © F Q В) = СтХт -\- • • • + CiejXiejA- * * *+ cmknxmkn'i
^ ==(Xlll' * * *» Xiej>* ’ *> Xmkn)’
a план будет опорным.
Выбор X с набором Ck будет оптимальным, если каждый из элементов c-iej, не входящих в Ck, имеет неотрицательную оценку
Д О О. Так как |
для |
любого элемента, |
не входящего в набор, |
|||
справедливо неравенство | c-iej — |
| > |
0, |
переход |
к любому |
||
другому набору C'k не может уменьшить |
значение скалярного |
|||||
произведения по сравнению со значением |
L, |
соответствующего |
||||
выбору Ck. А если это так, то план является оптимальным. |
||||||
Таким образом: |
(а > 1) экстремальные значения |
показателя |
||||
— в гибкой цепи |
||||||
линейной формы |
с^- |
однозначно |
определяют выбор |
сквозной |
•коммуникации; это следствие совершенно очевидно;
— отсутствие ограничений по пропускной способности в тран зитных пунктах ( а >1) сводит задачу определения оптимального плана в гибкой цепи многоэтапной транспортной задачи к оты сканию экстремальных значений показателя линейной формы с последующим решением матрицы ||А£|| одним из методов ли
нейного программирования. |
|
|
*Теорема 3. |
В жесткой цепи (а = 0), состоящей из r-f-1 звена, |
|
совокупность |
оптимальных планов звеньев: |
' |
28
составляет оптимальное решение для цепи, если Д*>0, Д е > 0 ,----
• • . , Ду > 0. |
|
цепи условия транспортировки |
Доказательство. В жесткой |
||
определяются равенством |
|
|
т |
к |
п |
/=1 |
^=i |
y=t |
а если это так, то достаточно доказать справедливость его для' одного звена, показав, что оптимальный план, определенный независимо для одного звена, есть сопряженная часть оптималь ного плана цепи.
Рассмотрим любой план первого звена
X = (х1у • • • , X/, • • • , хт) ф X .
Как и компоненты всякого плана, составляющие х* должны, удовлетворять условиям задачи:
А = |
TYI / - ч - / |
|
|
|
|
|
|
^ |
XiAe, х*>0; i — 1, 2, • • .• , m. |
|
|||||
|
i= |
1 |
|
|
|
|
|
Разложим векторы согласно формуле |
|
||||||
1 |
|
|
|
|
к |
|
(44) > |
< |
|
A t ^ X i e A , |
|||||
|
|
|
|
|
е=\ |
|
|
по векторам базиса, определяемого |
опорным планом |
||||||
|
|
|
ТП'—' |
/*/С |
|
|
|
|
|
F = Y , |
Xi |
( Ц х и Ае |
|
||
|
|
|
i—\ |
|
\е=1 |
|
|
Векторы А ь Л2, |
• |
• •, |
Ak |
составляют линейно |
независимую |
||
систему. Следовательно, |
вектор |
ограничений F |
может быть- |
единственным образом представлен как линейная комбинация. этих векторов. Сравнивая уравнения (44) и (45), получим
А\Хг А 2х 2-\-• • • -\-Amxm — В,
откуда заключаем, что
ТП ‘— '
£ xAei = Хе.
/=1
29;