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

книги / Математические методы принятия решений

..pdf
Скачиваний:
2
Добавлен:
13.11.2023
Размер:
22.94 Mб
Скачать

при условии, что

Y l Xj i - T l Xi j = d V ’

3

3

 

' % ij ^ U i j + V ij »

V ij ^ b i j ,

к

 

 

Xij > 0 ,

Vij ^

0 .

Описанная модель может быть использована для определения необходимой пропускной способности каналов, при которой задан­ ные требования выполняются с минимальными затратами. В этом

случае значения

обычно полагают равными нулю, если иц = 0,

а величины

выбирают достаточно большими. Другим примером

использования этой модели является задача минимизации стоимо­ сти передачи информации при фиксированных пропускных способ­

ностях дуг. В этом случае

= 0.

§ 5.7. Эвристический алгоритм решения задачи синтеза сети связи

Как уже отмечалось, сетевые задачи большой размерности ре­ шают, в основном, приближенными эвристическими методами. Рас­ смотрим один из эвристических методов на примере задачи синтеза некоторой сети связи.

Одна из основных задач при разработке схемы сети связи — это определение наивыгоднейшей конфигурации и распределения кана­ лов на сети. При этом могут решаться и другие задачи:

1)выбор типа линий связи (кабельные, радиорелейные, воздуш­ ные и т. п.) с соответствующими системами передачи;

2)последующее распределение стандартных каналов на сети. Структура сети —ее топология это совокупность пунктов

(узлов, станций и т. п.) и соединяющих их линий связи или каналов в их взаимном расположении. По структуре можно понять потенци­ альные возможности обеспечения связью отдельных пунктов этой сети.

Экономически оптимальный выбор сетки линий, их емкостей и плана распределения каналов при заданном наборе пунктов се­ ти называют задачей оптимального синтеза структуры сети. Ее

решение существенно зависит от того, насколько тщательно ото­ браны допустимые направления.

В качестве основных критериев оптимальности построения сети обычно рассматривают объем капитальных затрат, расход кабель­ ной продукции и т. д. Минимум расхода кабельной продукции обес­ печивается при минимуме общей протяженности трасс кабельных и радиорелейных линий связи. Иногда в качестве критерия опти­ мальности может рассматриваться и объем строительно-монтажных работ.

Реальная сеть должна отвечать одновременно в той или иной мере всем критериям оптимальности, т. е. задача определения наи­ выгоднейшей конфигурации и распределения каналов сети является многокритериальной. Однако в связи со сложностью решения мно­ гокритериальных задач обычно в качестве оптимальной принимают сеть, отвечающую одному, какому-либо заранее заданному, крите­ рию оптимальности. Очевидно, что сеть будет зависеть от выбран­ ного критерия оптимальности и в общем случае будет различной в зависимости от используемого критерия оптимальности.

Рассмотрим задачу синтеза сети при следующих требованиях.

1.Каждый пункт сети должен иметь два независимых пути (прямой и обходной) к центральному узлу (ЦУ). Прямой путь должен обеспечивать передачу по 70-75% каналов от их расчет­ ного числа, а обходной путь —по 25-30% каналов. Это соотно­ шение может измениться, а поэтому должно задаваться на входе задачи.

2.Должен существовать путь, соединяющий каждый пункт се­ ти с дополнительным узлом (ДУ) связи, не совпадающим с ЦУ. На некоторых участках этот путь может совпадать с прямым и об­ ходным путями.

Всетях связи налагают следующие дополнительные требо­

вания.

3.На линиях связи «пункт—ЦУ» число транзитов по высокой частоте (ВЧ) должно быть не более двух в случае применения ана­ логовых систем передачи.

4.Суммарная емкость (число каналов) на любом участке зоно­ вой сети должна быть кратна 12 для аналоговых систем передачи

и30 для цифровых систем.

5. Иногда требуется предусмотреть возможность того, чтобы че­ рез некоторые пункты проходило минимальное число путей. Эти пункты отмечают в исходных данных.

Задача синтеза структуры сети может быть наиболее адекватно представлена (в силу нелинейного ступенчатого изменения целе­ вых функций, в частности функции стоимости линии связи) как задача целочисленного нелинейного программирования в различ­ ных формах.

Рассмотрим задачу выбора оптимального варианта построения сети с двумя независимыми путями к единственному узлу (цен­ тральному узлу), к которому передается информация со всех других узлов сети. В качестве целевой функции рассматривается функция стоимости сети.

Задан неориентированный граф (сеть) G с N + 1 вершинами. Все вершины, кроме (N + 1)-й, являются источниками информации. Вершина N + 1 является стоком (центральным узлом). Путем ytk от вершины к к стоку N + 1 будем называть последовательность

вершин (к,i , j , . . . , N +

1), через которые проходит данный путь.

Для каждой вершины

к, к = 1,2, ... , N, необходимо выбрать два

пути к стоку, не имеющих общих вершин, кроме к и N + 1. За­

даны потоки информации Pi и Рг, которые требуется передать от вершины к к стоку соответственно по первому и второму путям. Дугам графа G приписаны веса, равные их длинам. Для общно­ сти рассмотрения заданный граф G дополним до полного дугами с достаточно большим весом.

Цель задачи — выбор топологии сети и маршрутов передачи ин­ формации на ней, минимизирующих общую стоимость сети. Об­ щая стоимость сети определяется как сумма стоимостей входящих

в нее дуг. Стоимость

дуги является функцией от длины дуги

(г, j) и суммарного потока информации,

передаваемой по

ней,—

суммы потоков от вершины i к вершине j

и от вершины j

к вер­

шине г.

 

 

 

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

получена по одному и тому же алгоритму при соответствующем выборе величины Су.

Введем следующие обозначения: x\j — булева переменная, рав­ ная 1, если вершина j непосредственно следует за вершиной г на первом пути [if от вершины к к стоку, и равная 0 в против­ ном случае, г = 1 ,2,..., N, j = 1,2,..., N + 1, г ф j, к = 1 ,2,..., N; x\j булева переменная, равная 1, если вершина j непосредствен­ но следует за вершиной г на втором пути (if от вершины к к стоку, и равная 0 в противном случае, г = 1,2,..., N, j = 1,2,..., N + 1,

г Ф j,

к = 1,2, ... , N; Q ÿ суммарный поток информации по дуге

(г, j),

определяемый формулой

N

2

 

 

 

Q ii = У

У

Pik(Xlh +

Xl£ ),

(5-15)

fct i

г=1

;

3

г = 1,2,

j = 2, 3, . . . , N

+ l,

i < j .

Тогда рассматриваемую задачу можно сформулировать следую­ щим образом. Необходимо найти значения булевых переменных х £1 , минимизирующих общую стоимость сети:

N

N + 1

 

5 = 2

2 C i j Q i j ^ n a i n ,

(5.16)

i — 1 j=i+\

где Qij вычисляют по формуле (5.15), при следующих ограниче-

НИЯХ:

 

N+1

 

N

 

 

 

 

2 Х%=

2 Хкг = °>

 

(5.17)

 

 

j

= 1

 

i — 1

 

 

N

j^k

 

 

 

 

 

 

 

 

 

 

2

*i5v+i = 1.

k = l , 2 , . . . , N ,

1 = 1 ,2 ,

(5.18)

2=1

 

 

 

 

 

 

N

N + l

 

 

 

 

 

 

2 ^ =

2

a;S - ^ 1’

r ,k

= \ , 2 , . . . , N ,

1 = 1,2,

(5.19)

2=1

j= 1

 

 

 

 

 

 

iФг

j ^ r

 

 

 

 

 

 

 

2

N

 

 

 

 

(5.20)

 

2

2

ar^

1’

r ,k = 1,2,

 

г=12=i 2фТ

Ограничение (5.17) означает, что /г —начало как первого, так и второго путей к стоку. Ограничение (5.18) предопределяет, что все пути должны оканчиваться в вершине N + 1. Ограничение (5.19) указывает на условия сохранения потока (равенство) и на запрет существования петель в путях от вершины к стоку (неравенство). Заметим, что если функция стоимости монотонно возрастает с ро­ стом суммарного потока, то ограничение (5.19) в виде неравенства можно исключить. Ограничения (5.20) реализуют условие непересечения в вершинах первого и второго путей от каждой вершины к стоку.

Данная математическая постановка предполагает наличие 27V3 булевых переменных и 5N 2 + 6N ограничений, что в практически интересных случаях (N > 20) не позволяет рассчитывать на точное решение.

Исключение составляет случай, когда функция стоимости ли­ нейна относительно суммарного потока. Поскольку целевая функ­ ция (5.16) становится линейной по переменным х 1£ , это позволяет осуществить декомпозицию исходной задачи на N независимых за­ дач булевого программирования с целевой функцией

 

N

N+1

 

 

 

S f = E

£

ï j P f x ÿ

(5.21)

 

2=1 .7=2+1

 

 

и ограничениями (5.17)-(5.20)

при

фиксированных

значениях к,

к = 1,2,..., N, где

— стоимость

единицы потока,

проходящего

по дуге (г, j).

 

 

 

 

Последняя задача (при фиксированном к) может быть сформу­ лирована достаточно просто: найти два пути от вершины к к вер­ шине N + 1, не имеющих общих вершин, кроме к и N + 1, и дви­ жение по которым обеспечивает минимум суммарной стоимости путей. При этом стоимость пути равна стоимости входящих в него дуг, а стоимость дуги (г, j ) равна для 1-го пути, 1 = 1 ,2 .

Математическую постановку задачи с двумя независимыми пу­ тями к центральному узлу и одним путем к дополнительному узлу можно записать в следующем виде.

Пусть вершина с номером N является дополнительным узлом. Через [ij обозначим путь от вершины к к ДУ, к = 1 ,2,..., N — 1.

При этом для любого k = \ , 2 , — I путь р3 не должен содер­ жать вершину N + 1. Пусть xjj — булева переменная, равная 1, если вершина j связана с вершиной г на пути [i3 от вершины к к ДУ, и равная 0 в противном случае. Через Р* обозначим требуемые ве­ личины потоков информации, передаваемой от вершины к к ДУ, к= 1 ,2 , . . . , Л Г - 1.

Тогда, придерживаясь прежних обозначений, суммарный поток на дуге (г, j ) определяем по формуле

Qu= 2 S

+xïù+ Ё рзН-+4-)•

(5-22)

fc=i;=i

fc=i

 

Для дуг (г, N + 1), г = 1,2,..., N , суммарный поток, как и ранее, вычисляют по формуле (5.15).

Рассматриваемую задачу можно сформулировать следующим образом. Найти значения булевых переменных x l£ , I = 1,2,3, ми­ нимизирующих общую стоимость сети, вычисляемую по форму­ ле (5.16) при ограничениях (5.17)-(5.20) и ограничениях

N

 

 

N-

1

 

 

 

2 4

$ = 1 ,

 

£

> ^

= 0,

fc = 1 ,2,..., JV — 1,

(5.23)

j = 1

 

 

2=1

 

 

 

 

j ^ k

 

 

1фк

 

 

 

 

 

TV— 1

 

 

 

 

 

 

^

1

^

=

1, Ac =

1,2, . . . , N — 1,

(5.24)

 

2=1

 

 

 

 

 

T V -1

TV- 1

 

 

 

 

 

 

£ 4 *

= £

a $

^

1,

r, fc =

1,2,..., TV —1, г ф к .

(5.25)

i= l

i = 2

 

 

 

 

 

 

Ограничения (5.23)-(5.25). аналогичны ограничениям (5.17)- (5.19). Аналог ограничения (5.20) отсутствует, поскольку к ДУ необходимо построить только один путь от каждой вершины.

Как легко видеть, требование существования путей к ДУ уве­ личивает число переменных на (N —I)2 и число ограничений на (N — 1 )(N + 2). Однако, как и ранее, в случае, когда функция стои­ мости дуги от величины суммарного потока линейна, задача распа­ дается на две независимые задачи:

1)определение двух непересекающихся путей, направленны

кцентральному узлу;

9— 4077

2) определение путей к вершине N , имеющих минимальную суммарную стоимость.

Математическая постановка первой задачи уже рассматрива­ лась, а вторая, как легко видеть, приводит к задаче построения кратчайших путей от каждой к-й вершины к N -й. В этих двух по­ становках случайными величинами являются значения Р к, так как они получаются в результате прогнозов и, в крайнем случае, долж­ ны задаваться в виде интервалов с соответствующими значениями вероятностей. Другой случайной величиной может быть доля ин­ формации, передаваемой по прямому (или обходному) пути.

Полученная нелинейная целочисленная задача может быть ре­ шена на современных ЭВМ лишь для сети малых размеров: число пунктов равно 10 при 20-30 допустимых трассах. Размеры реаль­ ных сетей, как правило, больше. Поэтому иногда проводят лине­ аризацию данной задачи, приводя ее к целочисленной линейной или просто к задаче линейного программирования. Однако в этом случае для получаемого решения не гарантируется определенная степень приближения критерия затрат к минимуму.

Большая размерность описанных задач и стохастическая приро­ да исходных данных приводят к необходимости разработки эври­ стических алгоритмов решения. Кроме того, сама специфика струк­ туры сетевых задач позволяет находить более эффективные алго­ ритмы поиска решения, чем алгоритмы в точных методах решения.

Известен ряд эвристических методов, в которых учитывается специфика поставленной задачи, но в этих алгоритмах для реше­ ния задачи используют дуги минимальной длины. Следовательно, алгоритмы ориентированы на расчет сети минимальной длины, что не всегда позволяет найти наилучший результат, тем более этот недостаток будет сказываться при построении коротких сетей, в ко­ торых расстояния между узлами малы, так как в этом случае стои­ мость сети определяется в основном стоимостью аппаратуры и ма­ ло зависит от длины трасс.

Алгоритм синтеза сети без Д У

Алгоритм для решения общей задачи (5.11)—(5.16) состоит из трех этапов:

I — построение первых путей от всех вершин к стоку;

= (к, N +1), p f(r)

II — построение вторых путей от всех вершин к стоку; III — улучшение сети, полученной на первых двух этапах. I. Построение первых путей от всех вершин к стоку.

Ш аг 1. Строят начальную звездообразную сеть, в шторой каж­ дая вершина непосредственно, без промежуточных узлов, соедине­ на со стоком, и считают стоимость такой сети согласно заданной функции стоимости. Если для какой-то вершины связь со стоком отсутствует в исходном графе, то стоимость такой несуществую­ щей дуги полагают равной достаточно большому числу.

Ш аг 2. Для каждой вершины к, путь от которой является ра­ диальным путем:

t f = ( k , N + I),

вычисляют экономию стоимости Т£, которая будет получена при замене пути р f на путь

=(к, г, р?),

иравна разности их стоимостей. При этом промежуточная верши­ на г принимает значения всех вершин, инцидентных вершине к (т. е. связанных с вершиной к разрешенными дугами), и вершин пути р*. Пути от вершины до стока не содержат несуществующих дуг. Затем

из всех экономий стоимости

выбирают максимальную.

На рис. 5.6 приведен пример такой замены. Экономия дости­

гается за счет ликвидации

дуги (к, N + 1), но при этом может

увеличиться стоимость дуг пути р р где pj = (г, гь гг, N + 1), р* =

= (к,г,ц, i2, N + 1).

Шаг 3. Выбирают максимальное значение Tfc* из всех эконо­ мий стоимости Т/с, к = 1,2,..., N . Если Tfc* ^ 0, то этап I считают завершенным.

Шаг 4. Для выбранного Tfc* осу­

ществляют замену пути pf* = (к*, N

+ 1)

 

на путь

pk*(r) = (к*, г, р^), на

котором

 

достигалась максимальная экономия стои­

 

мости Tfc*. Осуществляют замену путей

 

для всех

вершин

из множества V(k*),

 

пути от

которых проходили

по

дуге

 

( k \ N +

1), т.е.

вместо

пути

р^*** =

Рис. 5.6. Пример заме-

= (V(k*),..., к*, N + 1)

формируют

но-

ны дуги

вый

путь

=

На

рис.

5.6 такими

вершинами являются вершины j\

и j 2- Новый

путь, например,

для

вершины j\

будет ^ ’(г) = (ji, к, г, i\ ,i 2, N + I)

вместо р]1=

=( j u k , N + l).

Шаг 5. Экономию стоимости Т* полагают равной нулю для

всех к = 1,2, ... , N . Переход к шагу 2.

Шаги 2-5 выполняют до тех пор, пока может быть достигнута экономия при замене какого-либо радиального пути на путь, прохо­ дящий через промежуточную вершину. Пересчет экономии стоимо­ сти T*, осуществляется после каждой замены радиального пути.

Поскольку каждый раз происходит уменьшение на единицу чис­ ла радиальных путей, то этап I конечен.

И. Построение вторых путей от всех вершин к стоку. Идея ал­ горитма на данном этапе та же — замена радиального второго пути на путь, проходящий через промежуточную вершину.

Ш аг 1. Этот шаг аналогичен шагу 1 этапа I. Отличие состо­ ит в следующем: если первый путь от некоторой вершины после завершения этапа I оставался радиальным, то радиальный второй путь прокладывают по несуществующей второй дуге к стоку и его стоимость полагают вдвое большей, чем стоимость просто несуще­ ствующей дуги. Смысл этого состоит в том, что именно такие пути выгоднее заменять, а значит, они будут заменены в первую очередь. Кроме того, при расчете стоимости сети учитывают первые пути от всех вершин к стоку.

Ш аг 2. Этот шаг аналогичен шагу 2 этапа I. Только теперь существует две возможности замены радиального второго пути от вершины к на путь, проходящий через промежуточную верши­ ну г, а именно:

Р2Л (г) = (к, г, pï), р£’2(г) = г>Нг)>

т. е. новый путь идет от вершины к к вершине г, а затем либо по первому, либо по второму пути от вершины г к стоку (конеч­ но, при этом второй путь р2 должен быть «реальным»).

Обозначим через V(k) множество вершин, вторые пути которых проходят по дуге (к, N + 1), при этом keV (k) . Тогда путь и мож-

к \

к2

но заменить на путь |i2’

или [i2* только в том случае, если для

любого j е V(k)

n u * ’1e ( k , N + 1) или [i{ n (ij’2 e (k, N + 1)

соответственно. Сказанное выше схематично представлено на рис. 5.7.

Ш а г и

3-5

аналогичны

 

тем же шагам I этапа. Они вы­

 

полняются до тех пор, пока

 

замена какого-либо радиального

 

пути на путь, проходящий через

 

промежуточную вершину, может

 

дать положительную экономию

 

стоимости.

 

 

 

После

выполнения этапов I

 

и II мы получим

сеть с раз­

Рис. 5.7. Пример замены пути

меченными на ней маршрутами

 

первых и

вторых

путей от каждой вершины к стоку. Эта сеть

с маршрутами удовлетворяет всем условиям задачи, т. е. является допустимой. Кроме того, она обладает следующими особенностями.

1. Каждая вершина к = 1,2, . . . , N имеет не менее двух инци­ дентных вершин сети, поскольку осуществляет два непересекающихся по вершинам, а следовательно, и по дугам пути. Инцидент­ ными вершинами называют вершины, связанные разрешенными ду­ гами.

2. По каждой дуге (г, j) сети согласно выбранным маршрутам обоих путей происходит передача информации либо только в одну сторону, например от вершины г к вершине j, либо в обе стороны. В первом случае дугу (г, j) ориентируем от вершины г к вершине j, а во втором случае дуга (г, j) будет ориентирована в обе стороны. После ориентации сети описанным образом получим, что каждая вершина сети имеет ровно две выходящих из нее дуги (кроме, ко­ нечно, стока). Входящих дуг может быть произвольное число. Это следует из алгоритма замены радиального пути на путь, проходя­ щий через промежуточную вершину.

3. Если по дуге (г, j ) от вершины i к вершине j идет поток ин­ формации Pij, то в дальнейшем он полностью проходит по первому

Соседние файлы в папке книги