книги из ГПНТБ / Большие системы и управление
..pdf
|
|
|
|
|
|
|
|
|
н о |
|
|
|
|
|
|
жением (9 ). Если |
I + I |
= |
к |
, то члены матрица а! |
|
содержат |
|||||||||
все |
элементарные |
пути длины |
к |
меаду любыми двумя заданными |
|||||||||||
пунктами |
i |
и |
j |
(i ,j = |
1,2, |
|
т + п + п 1 ) . |
Следует заметить, |
|||||||
что число ненулевых элементов в |
А1К остается постоянным |
для |
|||||||||||||
определенного |
значения |
к |
я |
уменьшается для больших значений к • |
|||||||||||
Так, |
А |
состоит |
полностью из нулевых элементов для всех |
|
|||||||||||
к & т + п + п 1 (рис Л З ) . |
|
|
|
|
|
|
|
|
|
|
|||||
|
Очевидно, что минимальный путь по числу участков, ведущий |
||||||||||||||
из пункта |
I |
в |
пункт j |
, должен быть элементарным. |
Для опре |
||||||||||
деления минимального пути, |
ведущего из |
пункта |
с |
в пункт j |
, |
||||||||||
надо произвести следующие операции, |
|
|
|
|
|
|
|||||||||
|
I# Полагаем к= I . |
|
|
|
|
|
|
|
|
|
|
||||
|
2. Строим А,(к). |
|
|
|
|
|
|
|
|
|
|
||||
|
3» а) Если элемент |
|
|
равен нулю и /с< т + п + п'- 1 , |
то |
||||||||||
увеличиваем |
к |
на |
единицу |
и возвращаемся к пункту |
2; |
б) |
Если |
||||||||
|
|
|
|
|
|
|
|
|
элемент ( i , j ) равен нулю и |
||||||
|
|
|
|
|
|
|
|
|
к - т + |
7, то |
пути из |
пунк |
|||
|
|
|
|
|
|
|
|
|
та |
i в пункт |
J. |
не |
сущест |
||
|
|
|
|
|
|
|
|
|
вует; в) если |
элемент ( L , j ) |
|||||
|
|
|
|
|
|
|
|
|
не равен нулю, то он пред |
||||||
|
|
|
|
|
|
|
|
|
ставляет искомый путь (или |
||||||
|
|
|
|
|
|
|
|
|
пути). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, |
принимая |
||||
Рис. 14. Пример транспортной |
|
во внимание изложенное, |
мож |
||||||||||||
|
но предложить следующий спо |
||||||||||||||
|
|
|
сети |
|
|
|
|
||||||||
|
|
|
|
|
|
|
соб |
отыскания путей между |
|||||||
|
|
|
|
|
|
|
|
|
заданными пунктами, В соответствии со структурой транспортной сети строится ее
математический эквивалент - матрица инциденций А1 , отражающая
связи между заданными пунктами. |
|
, |
л 12 |
J т+п+п |
||
_ |
|
|
|
л |
||
По описанным выше правилам строятся матрицы |
А 0 |
А |
|
|||
d А) |
|
|
|
|
с и |
|
Тогда множество М |
элементарных путей между пунктами |
|||||
у определяется как |
сумма |
элементов а ' ц 1(к= 1,2,.. . 7/n+W -/)r*e. |
||||
|
/1( Ч ) |
_ m+n+n’-i |
аТ - |
|
|
(Ю) |
|
|
£ |
|
|
|
где суш а понимается в теоретико-множественном |
смысле. Рассмот |
рим иллюстративный пример. Определим множества |
элементарных пу- |
|
|
|
|
|
|
|
|
I l l |
|
|
|
|
|
|
|
|
тей для сети, изображенной на рис* 14* |
|
Исходная матрица |
А |
в |
||||||||||||
этом случае имеет |
вид: |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
0 |
|
и п |
|
и 13 |
|
0 |
0 |
|
\ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
и 21 |
|
0 |
|
и г ь |
|
и г ч |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
А = |
|
|
и 31 |
|
У 32 |
|
0 |
|
и з ч |
и з 5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
V |
0 |
|
U 42 |
<Чз |
0 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
0 |
|
0 |
|
US3 |
|
U S k |
0 |
/ |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||||||
Элементы |
(г) |
|
|
определяются как элементы произведе |
||||||||||||
а ц матрицы А |
||||||||||||||||
ния |
А( А1 и равны: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2) |
:= 0 ; |
|
|
|
|
(г; |
|
|
|
|
(2) |
6/;2 1 |
|
|
||
<*п |
|
|
|
а 72 |
= и 13 и32 ’ |
|
|
а гз |
= |
|
|
|||||
|
|
|
|
|
U23 ’ |
|
||||||||||
(2) |
|
|
|
|
|
(2; |
|
|
|
|
(2; |
= ^23 U311 |
|
|||
0 /4 =--ип и24^ и13изь'-> ^/5 |
~ а д 35 ’ |
|
|
^27 |
|
|||||||||||
(2) |
|
|
|
|
|
(2) |
|
|
|
|
(2) |
|
|
|
||
® 22 = 0 ; |
|
|
|
^23 ~ *J21U13V и2^ичз • |
а 24 “ ^23 изч-7 |
|||||||||||
(2) |
|
|
|
|
• |
о (2)- |
|
21 7 |
|
|
(2) |
|
|
|
||
^ 25 = игз U35 V У2-, |
|
? |
а з/ ■- У32 У |
|
|
^32 |
= U31U12^ U34- |
|||||||||
O'(2) |
= о-7 |
|
|
|
|
„(2) |
|
|
|
|
„(2) |
изч |
|
|
||
ц33 |
|
|
|
|
|
а 34 = U32 U2b ? |
|
|
°35 “ |
|
|
|||||
(2) |
|
|
|
|
|
(2) |
|
|
|
|
(2) |
|
|
|
||
<4/ = и« |
U27 V |
|
“ з/ i |
|
~ иьз ^32 ’ |
|
|
а 43 ~ и^г “ 23 V |
||||||||
V U i+5 и5Ъ 7 |
|
(2) |
|
|
'«(2)-/ |
|
|
|
VС2)-/_ |
|
|
|||||
а Ч.ц. = 0 » |
°Ь5 |
~ ЦЬзи35 7 |
^57 |
= U53 U31 7 |
||||||||||||
a 5 Z ~ U 5 3 U 3 2 ^ U 5 k U ^ Z :> а 5 3 ~ U 5 4 - U k 3 ’ |
^ 5 ^ ~ U 5 3 U 3 ^ ° 5 5 |
|
||||||||||||||
|
Элементы матрицы А |
,3 |
| 4 |
определяются аналогично, |
а |
эле- |
||||||||||
|
|
, А |
||||||||||||||
|
|
л15 |
|
/ (5) |
|
/л |
|
|
|
|
|
|
|
|
|
|
менты штршщ A |
|
q - |
= и . |
,4 |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
^.2 |
|3 |
можно определить множество |
|||||||||
Получив матрицы А |
• |
А , |
А , |
|||||||||||||
путей между любой парой пунктов: |
|
|
|
|
|
|
|
|
||||||||
М |
~ |
{ и ) 3 и 3 5 |
|
7 и 1 2 и 2 3 и 3 5 ' |
U 1 2 U 2 4 U h 5 7 |
U 1 3 u 3 4 - ^ 4 5 ’ |
и п и 2 ь и i+3 U 3 s 7 и 1з а з г и г ь и ± 5 •> и п и г ъ и зч- u * + s \ 9
112
м- { ^ 2 k ’ и г г аЗЫ U27 U13 U3*+? U23U35 U5iff UZ1U13U35U5bl
Если транспортную сеть представляется возможным изобразить в виде ориентированного графа (перевозки на каждом участке се ти совершаются в заданном направлении) без контуров, то отыска ние путей между заданными пунктами сводится к следующему* Для перечисления путей между двумя заданными пунктами 10 и ] 0 не обходимо перенумеровать узлы транспортной сети таким образом,
что если |
участок |
соединяет пункт L с пунктом j , то всегда |
. Алгоритм такого упорядочения пунктов на сети приводится |
||
в работе |
[4 ]. |
|
Рис.15. Пример ориентированной сети
После проведения нумерации определяется совокупность путей от одного пункта до другого по рекуррентной формуле, сущность которой состоит в том, что множество различных путей, проходя
щих через пункт J 0 и заданный пункт |
с0, равно объединению всех |
||||
путей, |
проходящих через |
пункты, |
для которых пункт |
является |
|
входом, |
т .е . это можно записать |
следующим образом: |
|
||
|
(i0 |
,Jo> = и м (I07i) |
|
||
|
М |
1-U;: 6 Ujn |
U Lio |
( И ) |
|
|
|
Чо |
Ао |
|
|
где и£0 - множество дуг, входящих в вершину J0 . Используя рекуррентную формулу ( I I ) , можно определить все элементарные пути между заданными пунктами [5 J.
и з
Р а с с м о т р и м п р и м е р о п р е д е л е н и я в с е й п у т е й м е ж д у п у н к т а м и I
и 8 д л я с е т и , и з о б р а ж е н н о й н а р и с . 1 5 .
Ми’г)= ц п -, Л?(7,3= и 13 ; М( Ш = М (иг)и 2, У М (ЬЗ)и3, V t / „ =
= u n u zlf\Jut3u3^\lu1lt’, M (h5>= M (,’‘f,u ltSW (,’3)u35 =
^>Jnu2ku^ slu,3u3fu^s ^иши‘и Уи13и35 ь H ° ’S>= М °’г)игеУ
W MAhb)a^6 = и}2 и26Уи12и^и^6 Уи13и3^и^6 \/u }¥u^6 ;
t (b7) |
„ М ) |
(he) |
|
Uz4-Ui+6 ^67 ^ |
|
M |
— M |
U27 V Л7 U67 - U12U 27V UJ2 u26 U67 |
|||
^ U13 У34 Ul+6 U67 V UJb Ui+6 U67 ’ M |
U58 ^ ^ |
Ut78~ |
~ Ujz u2if UifS u58 ^ U13 U3b Ub5 U58 ^ ^ lb Ub5 US8^ U13 U35 U58^U12U27Ur?8V
|
V u 12u 26 u 67 u 7JB. |
|
|
Таким образом, между пунктами I |
и 8 (рис*15) |
существует |
|
9 различных путей. |
|
|
|
Этот алгоритм может быть использован дня перечисления мно- |
|||
(I1in |
различных путей, проходящих через заданный уча |
||
жества ъ |
|||
сток сети (i\ j |
') при начальном пункте |
i 0 и конечном пункте j 0: |
|
|
г а',у> = м ИоЛ ‘) и М (У ’* ° ) , |
(12) |
Описанные методы перечисления путей для задач больших раз мерностей трудно реализуемы на ЦВМ» требуют выполнения такого объема вычислений, который часто не может быть выполнен в при емлемые сроки.
Поэтому для отыскания оптимального маршрута между заданны ми пунктами может быть применен метод случайного поиска, кото рый обладает следующими преимуществами:
-устойчивость к размерности задачи;
-статистический характер метода снижает влияние сбоев ЦВМ на результаты решения;
-програшы для реализации метода на ЦВМ сравнительно про
сты;
|
|
114 |
|
|
- решение |
может |
быть получено о |
наперед заданной точ |
|
ности)* |
|
случае отыскания |
множества |
путей меж |
Танин образом» в |
||||
ду заданными |
пунктами методом случайного поиска |
оптималь |
ный маршрут получается с некоторой вероятностно при за данном числе испытаний.
115
Кандидат технических наук Л.И. ГРИГОРЬЕВА
МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ТРАНСПОРТНОЙЗАДАЧИ ДО КРИТЕРИЮ "МИНИМУМПОТЕРЬ ГРУЗА В ПУТИПРИ ДОСТАВКЕ В ЗАДАННЫЕ СРОКИ"
Заданные пункты отправления и назначения соединены транс портной сетью, которая позволяет из любого пункта отправления попасть в любой пункт назначения (транспортная сеть связная).
Если транспортная сеть не является связней, то поставщики и по требители разбиваются тем самым на изолированные группы, для каждой из которых составляются индивидуальные планы перевозок. На сети, как правило, имеются промежуточные пункты. Предпола гается, что пропускная способность транспортной сети не огра ничена. С каждым участком л пунктом сети связывается величина, характеризующая надежность (вероятность потери груза при про хождении данного элемента транспортной сети). Перевозимый груз предполагается однородным и измеряемым одними и теми же едини цами (вагон, автомашина). Для пунктов назначения задается вре мя, в течение которого (начиная с заданного начального момен та) груз должен быть доставлен потребителю. Требуется так по строить план перевозок, чтобы груз доставлялся к потребителям в заданные сроки при минимуме потерь в пути.
Математически это можно записать в следующем виде. Имеет ся т пунктов отправления и п пунктов назначения. В каждом
пункте отправления i ( I = 1,2, ♦♦.,/77) |
находится запас a-Lеди |
ниц груза. В пункте назначения J ( j |
= 1,2, .. . , п) требуется |
bj единиц этого же груза. Для промежуточных пунктов j'Q = 1,2у.-,п) количество ввозимого груза равно количеству вывозимого. Неко торые пункты отправления и назначения также могут быть пункта ми, через которые груз частично перевозится.
lie
Пункта отправления, назначения и промежуточные пункта сое динены транспортной сетью, состоящей из N участков* Транс портную сеть можно представить либо в виде неориентированного
графа |
G (W, U) , где W - |
множество вершин графа (соответст |
вует |
множеству пунктов на |
сети ), a U - множество дуг графа |
(соответствует участкам сети ). На рис.16 изображена сеть, где W = AVBVC , А - множество пунктов отправления, В - множество пунктов назначения, С - множество промежуточных пунктов* Транс
портная сеть может быть также представлена в |
виде матрицы ин- |
||||||||
циденций |
\\ais || ( l ,s |
= t,2,...,m+n + n'). Наличие |
связи между пунк |
||||||
том |
I и |
пункте»! s |
фиксирует |
элемент |
als |
= |
u ls |
, |
элемент |
als - |
0 |
соответствует случаю, |
когда |
ь не |
связан |
с |
пунктом s . |
Рис *16* Изображение транспортной сети в виде графа
Для каждого |
участка г ( г= 1 |
,2 , |
.. * , /V) задано среднее |
время |
|
t r прохождения участка г , |
а |
для каждого пункта V (^ |
= 1 ,2 ,..* , |
||
р + т +■П 1 ) |
задано среднее время пребывания груза* Для пункта |
||||
назначения |
задается время |
|
tj |
доп - время, в течение |
которо |
го груз должен быть доставлен |
в пункт назначения. (За начало |
||||
отсчета можно считать, например, время подачи заявки на под |
|||||
воз) . |
|
|
|
|
- |
117
Обозначим через |
некоторый путь |
(вообще говоря, один И8 |
|||||||||||||
множества возможных маршрутов от пункта |
i |
до пункта J ), |
т .е . |
||||||||||||
последовательность |
участков |
транспортной сети |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
L1L г |
|
|
|
к- 1& |
|
(индексы |
|
i 1 lj, |
|
|
• ,lk_1}'jBce различны и пробегают номера тех |
||||||||||
пунктов, через которые проходит путь), |
соединяющую пункт |
L с |
|||||||||||||
пунктом |
j . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Требуется |
определить маршруты (i,j) |
от пунктов отправления |
|||||||||||||
до пунктов назначения и объемы перевозок на них х - так, чтобы |
|||||||||||||||
выполнялись следующие ограничения: |
|
|
|
|
|||||||||||
1. |
|
|
0 (условие |
неотрицательности объемов перевозок). |
|||||||||||
т |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 . £ |
|
X ’ • = Ь: |
л ~ 1, 2, . |
п |
(условие |
полного обеспечения не- |
|||||||||
i=i |
l* |
J |
|
|
|
|
|
|
|
|
|
|
|
||
обходимым грузом пунктов назначения). |
|
|
|
|
|||||||||||
3. |
п |
х - - |
^ |
а • |
|
(условие |
того, |
что общий планируемый |
|||||||
£ |
|
|
|||||||||||||
<»=». |
ч |
|
|
1 |
|
|
|
|
|
|
|
|
|
||
объем перевозок с данного склада не превосходит имеющегося на |
|||||||||||||||
нем запаса, |
- условие практической реализуемости плана). |
|
|||||||||||||
4. |
L |
m a x |
t L• |
ss t ’dot^Условие |
своевременной доставки грузов |
||||||||||
|
|
J . |
J оп |
|
|
|
|
|
|
|
|
||||
к пунктам назначения), |
где |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
к |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
£ |
1 |
( ^0 ~ ^ ’ ** = , |
|
||
|
|
|
|
|
|
|
|
|
|
Р=О |
|
|
|
|
|
|
|
|
fij) |
= |
| 1 > |
^ |
е |
|
» |
|
|
|
|
||
|
|
|
|
|
" |
[ о , |
|
|
|
|
|
|
|
|
|
5 . |
Функция |
|
т |
п |
Р ;; |
ал) |
|
|
- |
средние потери груза - |
|||||
|
£ |
D |
( О х ; . - |
||||||||||||
достигает |
|
|
|
|
L=f j = / |
|
Р — ( y J L^) |
- |
вероятность потери |
||||||
минимума. Здесь |
|||||||||||||||
груза на маршруте (ZJ). |
|
|
|
|
|
|
|
|
|
||||||
Совокупность чисел х - е Х я |
|
У |
» для которых выполня |
ются условия 1 - 4 , определяет допустимые маршруты и объемы
118
перевозок на них. План перевозок, для которого выполняется и условие 5, является оптимальным.
Заметим, что следствием условий 2 и 3 является неравенство
/77 |
П |
|
Е О; * |
Е % . |
(D |
При построении решения это не |
является принципиальным, так |
как введением фиктивного пункта назначения указанное неравен ство сводится к равенству.
Если суммарные запасы пунктов отправления меньше общей по требности пунктов назначения, то для приведения к равенству вводится фиктивный пункт отправления иш+ I й, в котором сосре доточено
л |
т |
а т + 1 ~ ^ |
^ |
<*=1 |
1=1 |
единиц однородного груза, |
а маршруты (m + i,j) |
имеют вероятно |
|
сти потери |
Рт+ц = 0 . |
|
|
Если в |
выражении (I) |
имеет место строгое |
неравенство, то |
для получения-равенства вводится фиктивный пункт назначения "п + Г о потребностью
Jп + 1 = |
Е <7; |
“ Е 6; . |
|
||
|
i=l |
L |
j =/ |
|
|
а маршруты ( i,r>+1) имеют вероятности потери PL n+J = 0 . |
В даль |
||||
нейшем будем считать, что |
в |
условии |
2 имеет место равенство |
||
|
£ |
О; = £ |
Ь: . |
(2) |
|
|
|
|
<и |
* |
|
Критериальная функция (условие 5) представляет собой нелиней ную функцию целочисленных искомых переменных • Таким образом, сформулированную задачу построения плана перевозок с минимумом потерь груза в пути при доставке в заданные сроки можно отнести к задачам нелинейного целочисленного программи рования.
В настоящее время не существует достаточно эффективных точ ных методов решения задач нелинейного целочисленного программи рования и обычно применяют следующие приближенные методы:
119
-увеличением числа переменных (для случая квадратичного программирования в два раза), преобразовывают задачу нелиней ного программирования в задачу линейного программирования;
-отказываясь от целочисленности, применяют различные гра диентные методы, округляя затем тем или иным способом получен ные значения результатов до целых чисел.
Известные методы решения преобразованной линейной задачи требуют значительного числа вычислений и обладают большой связ ностью. Решение подобных задач даже при использовании современ ных быстродействующих ЦВМ вызывает определенные трудности (боль шой объем вычислений, необходимость хранить большое количество исходной и промежуточной информации). Однако для решения по ставленной задачи не требуется привлечения математического ап парата, предназначенного для решения нелинейных целочисленных задач.
Анализ критериальной функции и условий 1 - 4 показывает,
что можно минимизировать сначала коэффициенты, зависящие от
(т#е. построить оптимальные в смысле предлагаемого кри
терия маршруты), а затем, зафиксировав переменные решить задачу линейного программирования. Это вытекает из следующего утверждения.
Пусть имеется набор функций fk {y) |
и Цк(зс)9 таких, |
что |
|||||
fk( y ) * 0 и |
Чк (х ) * ° - |
|
|
|
|||
Пусть далее |
у |
таковы, что при |
х = 5с и у = у T.fk(y)qk(x) |
||||
достигает минимума, т .е . |
|
|
|
||||
m tn |
T . f k ( y ) q k (x ) = ' Z f k { y ) q H { x ) . |
(3) |
|||||
у , |
о с |
|
к |
|
к |
|
|
Если в |
выражении |
(3) |
справа заменить fk (у) минимумом |
||||
m in ?к (у) |
~ ?к(У) |
9 т0 |
в СИЛУ неотрицательности функций по |
||||
лучим |
|
|
|
|
|
|
|
|
|
|
£ fk ( y ) 4 k ( x ) э £ fk ( y ) q k ( x ) , |
W |
|||
|
|
|
к |
|
|
к |
|
****** fk ( y ) *
Но, с другой стороны, всегда