Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Большие системы и управление

..pdf
Скачиваний:
10
Добавлен:
29.10.2023
Размер:
11.27 Mб
Скачать

 

 

 

 

 

 

 

 

 

н о

 

 

 

 

 

 

жением (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 ) *

Но, с другой стороны, всегда

Соседние файлы в папке книги из ГПНТБ