
книги из ГПНТБ / Берзин Е.А. Оптимальное распределение ресурсов и элементы синтеза систем
.pdfдать достаточно высокой точностью, так как между задачами суще ствует определенная аналогия.
Рассмотрим частный случай, который подтверждает аналогию и может дать представление о возможной погрешности. Этот случай со
ответствует условию |
|
|
С0Ѵ(, |
если ] = і, |
( 2.200) |
Wjv/ — О, |
если ]Ф і. |
|
Предполагается также, что L = |
S. |
|
Физически это условие означает, что каждое активное средство |
||
при назначении в /-ю точку прицеливания (г = /) |
с вероятностью юѵг |
может поражать только і-й объект (точки прицеливания и объекты совмещены, а объекты находятся далеко друг от друга). По физичес кой постановке задача совпадает с задачей § 2.1(1) и остается прове рить, совпадут ли их алгоритмы. Для этого достаточно убедиться, что совпадут формальные выражения для приращений Ай и —А*).
То, что выражения для Ай совпадают, следует непосредственно из сравнения (2.188) и (2.7) при условии (2.200). Убедимся, что совпа дают и приращения — Айі.
Подставляя условие (2.200) |
в формулу в (2.192), |
имеем: |
|
|||||
- А й = - 3 ^ - 1)П в ѵг( - ^ - і ) * = |
|
|||||||
|
і |
V |
\ |
eki |
J |
|
|
|
— |
г ^ |
' - ' Ч |
П |
^ |
— . |
|
|
(2.201) |
|
І ф і |
' |
V |
/ |
&hi |
|
|
|
что и требовалось [см. (2.13)]. |
|
|
|
|
|
|
|
|
Таким образом, задача §2.5(1), |
как |
частный |
случай, |
включает |
||||
в себя задачу §2.1(1), |
и оценки точности, |
проведенные для |
алг. 2.1, |
|||||
справедливы и для алг. 2.5 при условии (2.200). |
|
|
|
|||||
Другим обстоятельством, |
увеличивающим |
аналогию, |
является |
тот факт, что и для условий данной задачи сохраняется ранее уста
новленная тенденция к слиянию алг. 80 и 8+ по мере увеличения раз мерности задачи. Такой вывод можно сделать на основе анализа фор
мулы |
(2.194). |
|
|
Действительно, при увеличении числа средств или их эффектив |
|||
ности |
второй член в круглых скобках (2.194) стремится к нулю, по |
||
этому |
|
|
|
min Aw да min |
1)еш ~ т а х ^ Л і і |
!>сош = max A ^ . (2.202) |
|
kJ |
k,l І |
kJ l |
kJ |
Требуемый для счета по алг. 2.5 объем вычислений характе ризуется следующими оценками количества элементарных операций:
умножений да SLN2,
сложений да SLN2,
сравнений да LN2.
* При і = I имеем (0 ;ъі = Шь/ и выражение в скобках обращается в нуль
(при іф і Bihi—1)> см. (2.200)).
90
Задача § 2.5(1) и алгоритм решения имеют много общего с неко торыми из предыдущих рассмотренных задач, что дает нам основание просить читателя обратиться к §2.1(4), где изложены возможные ва риации задачи и соответствующие им изменения в алгоритме, которые приемлемы, на наш взгляд, и для данной задачи.
Рассмотрим пример, исходные условия для которого заданы че тырьмя матрицами. Эти матрицы, по сути дела, являются одной из форм представления объемной матрицы ||e,-VJ- Цз/ѵП
{АІ} = {А1 = 50, |
|
Л* = |
20, |
||||
0,8 |
0,9 |
0,4 |
k |
0,4 |
0,7 |
0,6 |
|
k , V 0 ,4 | 0,9 |
0,9 |
0,9 |
0,9 |
0,1 |
|||
V |
|||||||
0,8 |
0,9 |
tO,2 |
|
0,9 |
0,2 |
0,8 |
|
|
j,l |
|
|
|
I |
|
Л = |
20, |
|
Л ,= 10}, |
|
|||
0,8 |
0,2 |
0,8 |
|
0,3 |
0,3 |
0,5 |
|
k 0,5 |
0,3 |
0,9 |
k |
0,2 |
0,1 |
0,9 |
. (2.203) |
V |
0,2 |
0,8 |
V |
0,1 |
0,4 |
0,5 |
|
0,9 |
|
|
|||||
|
I |
|
|
|
I |
|
|
Не будем приводить вычислительную схему, поскольку она до статочна проста и близка к тем, которые показаны в предыдущих за дачах, но остановимся на анализе полученного решения.
Ниже записано решение (матрица 1) и последовательность рас пределения средств, полученные в соответствии с алг. 2.5. Матрица 2 отображает оптимальное решение, а матрица 3 дает представление об одном из возможных худших вариантов решения. Соответствующие
этим вариантам значения целевой функции |
записаны под |
матрицами |
|||||||
М а т р и ц а 1 |
М а т р и ц а 2 |
М а т р и ц а 3 |
|
||||||
0 |
І2 |
0 |
0 |
0 |
12 |
1 |
0 |
0 |
(2.204) |
0 |
0 |
1з |
k, V h |
0 |
0 |
k , V 0 |
0 |
1 |
|
0 |
0 |
І ! |
0 |
13 |
0 |
1 |
0 |
0 |
|
|
I>/ |
|
|
It / |
|
|
І |
|
|
4>! = |
86,55, |
ф2= |
88,63, |
ф з= 57,28 |
|
Полученное решение отличается от оптимального примерно на 2,4%, причем других распределений, более близких к оптимальному,
нет. Интересно отметить, что и алг. 6+ дает в данном примере ту же
последовательность назначений, что и алг. 60.
Полученная погрешность в данном примере близка к своему мак симальному значению. Это следует из анализа матрицы |)Aw||wl д л я
первого шага процесса:
79,2 |
68,8 |
54,3 |
|
|
k 51,01 J |
72,5 |
94,7 |
— 1 ^hl |
(2.205) |
85,0 |
60,24 |
ф 48,2 |
|
|
Из матрицы видно, что А33 = 48,2 и лишь немного меньше элемен та Д21 = 51,0. Незначительное уменьшение параметра е121 = 0,4 или увеличение параметра е133 = 0,2 [см. (2.205) и (2.203), жирный шрифт] приведет еще к некоторому увеличению погрешности и затем к нару шению условия
^33 ^ ^21- |
(2.206) |
91
Нарушение условия (2.206) означает, что первое назначение бу дет теперь произведено в соответствии с оптимальным решением (62і = 1). Последующие назначения, как показала проверка, также соответствуют оптимальному решению, даже в случае использования алг. 6+.
Анализ исходных условий примеров, в которых возникает по грешность, и промежуточных результатов вычислений наводит на мысль, что одной из причин погрешности для данной задачи может быть многократный (L-кратный) учет каждого средства по каждой из целей. Это приводит к занижению абсолютной величины приращения —Ай [см. (2.192)3 и может явиться причиной сбоев алгоритма. Как один из возможных способов компенсации этого недостатка можно рассматри вать использование в качестве суммарного приращения линейной ком бинации приращений Ай и —Ай не с равными, как ранее, а с раз личными коэффициентами:
АЫ= Д ^ - | А 7 / , і > 0. |
(2.207) |
В условиях данной задачи можно принять | = L. Путем эквива лентных преобразований выражение для суммарного приращения при водится к виду:
П< |
( f - i ) |
= 2 л Г '> Hkl 1 + L |
(2.208) |
і = 1 |
atu |
В остальном алг. 2.5 остается без изменений.
Применительно к рассмотренному числовому примеру решение для первого шага процесса дает следующую матрицу || Akl [|, полу ченную с помощью (2.208).
75,5 |
73,4 |
82,6 |
|
|
54,5 |
77,2 |
81,4 |
6 3 3 — |
(2.209) |
90,3 |
70,0 |
51,8 |
|
|
Мы видим, что хотя решение и не изменилось, но относительная разница между конкурирующими элементами матрицы || Akl || умень шилась (с 5,5 до 5,0%). Это говорит о сокращении области возможных сбоев и об уменьшении погрешности алг. 2.5.
Заключительное замечание. Рассмотренный в данной главе метод ДФ по существу является способом, позволяющим путем квазиэквивалентных преобразований исходной задачи, свести ее к условиям, при которых она может быть эффективно решена методом МЭ. Пока зателем степени эквивалентности преобразований в большинстве слу чаев являлась оценка точности алгоритма.
Полученный опыт применения метода ДФ позволяет перейти к ре шению более сложных задач, для получения представлений о точнос ти решения которых мы будем иметь возможность пользоваться еще более ограниченными средствами. Основными из них будут являться аналогии, частные экспериментальные оценки (числовые примеры), предельные оценки и качественный анализ решений.
Г л а в а т р е т ь я
МЕТОД ВЕСОВЫХ КОЭФФИЦИЕНТОВ
И я больше всего дорожу Аналогиями, моими са мыми верными учителями, Они знают все секреты Природы и ими меньше всего следует пренебрегать в Геометрии.
Кеплер
По-видимому (а может быть, — тем более), не следует пренебре гать аналогиями и в современной математике. Аналогиями мы окру жены в повседневной жизни и часто пользуемся ими, даже не замечая этого. Однако при решении трудных задач нужно искать аналогии. Удачно найденная аналогия может явиться ключом к решению про блемы.
Решения задач § 3.2, 3.3, 3.4 получены на основе формальной ана логии с задачей §3.1(1). Проведенные оценки лишь подтверждают, что точность решения по сравнению с § 3.1(3) по крайней мере не ухуд шилась.
§ 3.1. Оптимизация системы с нелинейной
двухиндексной целевой функцией аддитивно-мультипликативного вида
при независимых элементах системы и метод ВК
1. Постановка задачи. Требуется |
определить |
матрицу у0 = |
||
— ІТп Цлілг, доставляющую минимум функции вида: |
|
|||
|
|
|
(3.1) |
|
при ограничениях на переменные |
|
|
|
|
N |
|
|
(3.2) |
|
2 уг] — 1 > |
г = 1. |
•••> м, |
||
и дополнительных условиях |
|
|
|
|
Уті € {0; 1}, |
г = |
1.......М, |
|
|
0 <(<7tf= 1—Р г іХ 1 |
І ~ <-£> сг + 1, • |
(3 .3 ) |
||
0 < ( с О ; = 1 - в ;) < 1 , |
||||
|
|
|
t'= l, ..., S.
93
В соответствии с принятыми обозначениями будем иметь в виду,
что
|
Сі |
= |
иг_i + l |
|
{с1 = 1, |
us = |
N), |
|
|
(3.4) |
||||
|
|
|
|
|
г |
— |
1, . . |
S. |
|
|
|
|
|
|
Из (3.4) следует, |
что вообще |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
/ = |
1......... N, |
|
|
|
|
|
(3.5) |
||
однако в силу того, |
что / = /'(г) для каждого фиксированного значения |
|||||||||||||
і переменная / (где значения |
prj отличны |
от нуля) будет изменяться в |
||||||||||||
более узких пределах, чем (3.5), а именно: |
|
|
|
|
|
|
||||||||
|
/(О = |
|
си |
ct + |
1, |
. . ., Ui, |
i |
= |
1, . . ., 5 . |
|
|
(3.6) |
||
При этом дополнительно обозначим: |
|
|
|
|
|
|
||||||||
|
rii = Ui — сг |
+ |
1, |
і = |
1, . . ., S. |
|
|
(3.7) |
||||||
Приведем одну из интерпретаций, которой в дальнейшем будем |
||||||||||||||
пользоваться для |
|
логического |
обоснования метода. Имеется S |
объ |
||||||||||
ектов. Вес (ценность) t-го объекта измеряется величиной A t (і |
= 1, . . . |
|||||||||||||
. . . , S). |
Объект с номером і |
(і = |
1, . . ., |
5) |
атакуется |
группой щ |
||||||||
целей, |
которые |
нумеруются последовательно возрастающими номе |
||||||||||||
рами / от Ui до Сі, начиная с |
= |
1. Таким образом, номер цели / связан |
||||||||||||
с номером объекта, |
|
который эта цель поражает с вероятностью сoj не |
||||||||||||
зависимо от других целей. Имеется |
|
М активных |
средств, |
|||||||||||
обороняющих объекты. |
Условная |
вероятность |
поражения |
г -м |
(г = |
|||||||||
= 1, |
. . ., М) средством /-й (/ |
= |
1, . . ., |
N) цели равна |
рг} (матрица |
\\Ргі\\лш задана).
Полагая, что вероятности <а^ известны, а также известно, какой объект атакует каждая цель, необходимо так распределить средства, для воздействия по целям (найти матрицу назначений у0). чтобы об щий ущерб объектам, выражаемый формулой (3.1), был минимален*. При S = 1 задача вырождается в задачу § 2.3(1). Рис. 4 поясняет при веденные обозначения.
Другая интерпретация может состоять в оптимизации надежной системы, состоящей из 5 подсистем. В этом случае под величиной Аі следует понимать условную вероятность выхода из строя всей системы, если прекращает функционировать ее і-я подсистема.
Оптимизация надежности состоит в оптимальном резервировании элементов подсистем.
Основное отличие данной задачи от задачи надежности, рассмот ренной в § 2.3(1), состоит в том, что здесь не учитывается в явном виде влияние структурной схемы системы на ее надежность. Такое влияние может учитываться через условные вероятности Аі? получаемые на ос нове моделирования системы. Меньшая степень зависимости матема тической формулировки данной задачи от структурной схемы оптими
* По смыслу задана строго противоположна задаче §2.1(1), где ущерб максимизируется путем оптимального распределения атакующих средств.
94
зируемой системы говорит о большей ее общности и пригодности для синтеза надежных систем со структурой более общего вида.
2.Метод весовых коэффициентов и алгоритм решения задачи
Существенное отличие функции (3.1) от ранее рассмотренных состоит в том, что она является выпукло-вогнутой. Непосредственное примене ние изложенных выше методов для ее оптимизации невозможно хо тя бы уже потому, что функция не может быть так просто (логарифми рованием) сведена к аддитивной, что позволило бы применить для ее
оптимизации ранее изложенные методы. |
Тем не |
менее, сведение ее |
|||||||
к аддитивному виду все-таки возможно. |
|
|
|
|
|
||||
:1 |
г . |
Uf, |
cz . . . uz .. • Ci . J (i) .. U{ . ..CS . . . N |
(3.8) |
|||||
Bf в2 . • • BUf? |
|
|
|
|
|
|
|||
................ |
Bq |
ü) Bu. .• ■Bcs Bfi |
|
||||||
tüf |
СО2 • • • |
. |
шсг. • • ■<^az” M q |
U)j |
|
ü)Uj ... Cüq |
|
||
1 |
1• • |
• iI |
І І ... І .... |
L |
i . |
, |
. ] . . |
|
|
|
|
|
Y ^ |
V |
У |
1 |
|
|
|
|
Y |
|
|
|
^ |
|
|
||
|
Af.; mi ßf |
Az ; nz ; B2 |
Аi ;j);пі ;ßj |
As;ns ;ß„ |
|
Рис. 4.
При этом мы будем использовать первую содержательную поста новку задачи, приведенную в § 3.1(1). За основу преобразования возь мем соотношение (2.149), которое процесс максимизации произведения функций позволяет заменить максимизацией суммы функций [(см. (2.149)]. Опираясь на физическое содержание задачи, можно сказать, что процесс минимизации ущерба объектам заменяется процессом мак симизации числа сбитых целей, атакующих эти объекты, с учетом важ ностей (весов) этих целей. Подставим (2.149) в (3.1):
min F (у) ~ |
|
S |
Аі |
и і / |
|
м |
qyy |
|
шах |
2 j |
П |
1 — |
П |
||||
У |
У i = 1 |
j — ci \ |
|
r = 1 |
1 |
|||
|
S |
_ |
u i |
I |
|
M |
|
\ |
^ m a x |
S |
A ß t |
S |
И/ |
1 - |
П qyry |
= |
|
V |
i = 1 |
|
i = ct |
\ |
|
r = 1 |
|
/ |
= |
max |
i f |
|
f 1 — П <7г/ф |
|
(3.9) |
||
|
V |
|
' |
r |
1 ' |
|
|
|
где Bj(/)= Ю/(,)= AiBi (üj(,)— коэффициент |
важности (весовой |
|||||||
коэффициент) /-й цели из группы пи атакующей і-й объект. |
||||||||
Примечание. Если |
компонента |
в>] |
коэффициента |
Bj |
определяет относи |
тельную важность /-й цели в составе группы пі, то коэффициент В% определяет относительный вес всей группы целей в.составе остальных группы (см. рис. 4).
С учетом (3.4)и (3.5)выражение (3.8) приводится к окончательному виду.
|
N |
/ |
М |
\ |
|
min F (у) |
шах 2 |
в / т ( 1~ |
ГІ |
W ) . |
(3.10) |
V |
V /= і |
V |
г= 1 П / |
|
95
Таким образом, исходная задача распределения М средств для воздействия по N целям из условия обеспечения минимума эффекта на S объектах заменена задачей определения матрицы 1у"/||млг = у0, обеспечивающей максимум эффекта на N целях с учетом их важно стей Bj. Остается теперь так определить значения коэффициентов В}
в зависимости от исходных параметров, чтобы для обеих задач реше ния (уо) совпадали. Формальными методами этого достигнуть невоз можно, так как, по сути дела, здесь мы касаемся вопроса выбора и обо снования непротиворечивости различных критериев для решения од ной и той же задачи. Обратимся к анализу физической стороны задачи.
Факторы, определяющие величину Вщ), должны зависеть по крайней мере:
от эффекта, который может обеспечить данная цель на «своем» объекте (Л ,©,(«);
суммарного эффекта (ущерба), достигаемого группой пь целей на
г-м объекте; |
|
|
|
|
|
|
|
|
|
|
|
|
|
количества п; |
целей в і-й |
группе. |
|
|
|
|
|
|
|
|
|||
Характер зависимости B^j) от первых двух факторов должен быть |
|||||||||||||
таким, чтобы по мере |
увеличения значения каждой из |
соответствую |
|||||||||||
щих величин увеличивалась бы и величина Вщ). |
|
|
|
|
|||||||||
Зависимость от величины nt должна носить противоположный |
|||||||||||||
характер, так чтобы при |
яг-> оо величина Вц^ |
асимптотически стре |
|||||||||||
милась к нулю. Можно видеть, что выражение (3.11) |
удовлетворяет |
||||||||||||
всем предъявленным требованиям: |
|
|
|
|
|
|
|
|
|||||
Л і |
1— П |
№%) |
A i |
|
1 — |
П |
СО; |
|
|
||||
в і (О |
|
ß= C; |
|
|
|
|
Ц = гг |
|
|
(3-11) |
|||
|
Уі Аі соц |
|
|
|
|
и і |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
S |
ffl4 |
|
|
|
|||
|
|
4= 0 |
|
|
|
|
4= 0 |
|
|
|
|
||
/ |
= |
1, |
..., |
N, |
i = 1, ..., |
S |
(i |
= |
i(j)). |
|
|
|
|
Зависимость |
B j |
от |
n t |
проявляется |
через |
количество |
слагаемых |
||||||
в знаменателе и может стать наглядной, |
если |
положить |
со,- = |
со = |
|||||||||
= const для всей группы щ (/(г) = си ..., |
ut). |
Тогда из (3.11) |
полу |
||||||||||
чим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В;ц) = |
АЛ1 - е |
п 0 |
|
|
|
|
|
|
(3.12) |
Таким образом, если величины Bj определены, то решение сво дится к максимизации функции (3.1 0 ), что и достигается с помощью ме
тода ДФ. Однако в ходе проведения вычислительного эксперимента было установлено, что в некоторых случаях использование формулы (3.11) приводит к значительным погрешностям. Анализ исходных дан ных таких задач показал, что погрешность возникает в тех случаях,
Подробнее о весовых коэффициентах см. в приложении II.
96
когда в составе одной или нескольких групп целей nt имеется хотя бы одна такая, для которой выполняется условие
(<°j |
qrj = |
1 |
(3.13) |
Физически условие (3.13) означает, что данная цель гарантиро ванно (bj —у1 ) достигает г-го объекта и обеспечивает максимум эффекта
на нем (со/-»-1). Но если сохранить і-й объект невозможно, то формула для вычисления весовых коэффициентов Вj должна автоматически обеспечивать важность, равную нулю для каждой цели из группы пи в составе которой находится та, для которой выполняется условие (3.13). Легко проверить, что формула (3.11) этого не гарантирует, так как она нечувствительна к величинам qTj. Новому требованию удовлет воряет формула:
|
|
|
( |
иі |
|
|
и і |
|
|
|
|
|
Ai(djl |
П (1—СОцбц)— п |
|
|
|
||||
|
в і U) |
|
= |
____________ Ѵ= сі |
|
|
|
|||
|
|
|
и і |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
2 |
°Ѵ |
|
|
|
|
|
|
|
|
|
Я = |
С ; |
|
|
|
|
|
A i |
CD; Гі |
|
0) ^ а д - а г ), |
/ = |
|
|
(3.14) |
||
|
- ^ - ( П ( 1 - |
|
|
|||||||
где |
|
|
|
|
|
|
|
|
|
|
|
М |
|
|
и і |
|
|
|
и і |
|
(3.15) |
Ьц = |
П |
qrß> |
Q-i~ |
П |
вjx (/), |
Qi ( / ) = |
2 |
й)ц. |
||
|
r = |
l |
|
n = Cf |
|
В = |
|
|
||
Действительно, если теперь со, Ьх= |
1 (е, = 0), |
/ |
£ ХПі **, то вы |
|||||||
ражение, заключенное |
в |
квадратных |
скобках |
(3.14), |
приводится |
|||||
к виду: |
|
|
|
|
|
|
|
|
|
|
ß / < o ~ 8 z ( |
гі |
( 1 — (ОцЬц)— |
П |
|
= 0, }: |
|
|
(3.16) |
||
|
|
|
|
|
и Ф I |
|
|
|
|
|
что и требовалось.
Формула (3.14) удовлетворяет всем ранее предъявленным требо ваниям, что следует из ее сводимости к виду (3.11), когда «угрожаю
щей» ситуации, |
связанной с условием (3.13), не существует (т. е. |
С О ,& ,->0 ДЛЯ I £ |
Х г ц ) . |
Перейдем к выводу основных соотношений, позволяющих сфор мулировать алгоритм решения задачи.
По аналогии с § 2.1(2) [см. (2.21)] запишем выражение для суммар
ного текущего приращения А к о т о р о е получает |
функция Ff на |
||
* |
Это |
возможно, если одновременно выполняются |
условия со,-->1 |
( E j 0 ) |
И |
Ь ] - > 1 . |
|
** Множества SSn _ составляют номера целей группы гц.
4 Зак. 1292 |
97 |
і -м шаге процесса в результате пробного назначения k-то средства по I-й цели:
Ai/, = ß{<“ ,) Pw - 2 |
1 a- |
Pk!~bi/ - 1), k e M W , |
(3.17) |
|
|
i*t |
Q h j |
|
|
|
неиспользованных к t-му шагу средств; |
|||
где МО) — множество |
||||
В« -О |
(у) ®/ |
|
|
1) |
У |
//V |
п |
( і - < » Г ,Ч ' - І)) - Л |
‘гл |
|
й< (у) ■ |
\ ^ = fi |
|
1 |
|
|
І = 1, ЛГ. |
(3.18) |
|
Выражения для текущих значений величин о//*, |
6 ,(>, ау^, 0,\() по |
аналогии можно получить на основе предыдущего, поэтому здесь мы запишем их без выводов в виде рекуррентных соотношений*:
|
|
|
a>f_1), |
если |
i=f=lt, |
** |
|
(3.19) |
||||
|
®/° = |
' |
|
l)4kt j, |
если |
j = lt, |
|
|
|
|||
|
.®/ |
|
|
|
|
|
||||||
|
b» = ^ L — |
' |
^ i - l , |
...,N |
(см.(3.15)), |
(3.20) |
||||||
|
|
ЧДу n |
|
|
|
|
|
|
|
|
|
|
|
|
если |
|
/( < с г, |
или |
/< > и г |
[см. |
(3.15)1, |
||||
J t ) |
- |
° |
V |
, |
„(У—1) |
|
|
|
|
|
|
(3.21) |
йі |
— |
|
|
|
|
|
|
|||||
|
1 |
- с о ^ |
- 1 ) |
|
‘ |
! _ |
“ ( ! - Г ) ’ c t |
< |
h |
< |
t t t , |
|
|
|
Н |
|
|
|
|
it |
|
|
|
|
|
о ( |
у ) _ [ 0 <( |
если |
lt<Ci> или |
/ у > « г |
[СМ. |
(3.15)], |
' ”
Начальные
( 0 )
Pkth, |
если |
сг< / ; < |
Мг. |
|
||
значения величин |
вычисляются |
по формулам: |
||||
З.о). |
®у. |
/ = |
1 , |
іѴ, |
|
(3.23) |
® Г = |
|
|||||
и і |
|
|
п |
|
|
/ = 1 , ..., S, |
і Иц, |
«І0) = |
( 1 - с о Г ) , |
||||
ц=сг |
м |
|
|
|
|
|
|
Р - = 1 ,...,уѴ. |
|
(3.23') |
|||
т= 1 |
|
|||||
|
|
|
|
|
Начальные значения весовых коэффициентов определяются по формуле
В(°) = W Ü Y |
f t |
( ! _ < * & « ) - « % ) . У =1 ......«• (3.24) |
Q i ( i ) |
V'i = f « |
/ |
Теперь можно сформулировать алгоритм решения задачи.
* Эти соотношения получаются из (3.15) путем сравнения (3.15) дои после назначения на t -м шаге процесса fty-го средства по /у-й цели.
** соО) = Q [°< i> j (2.177), 6<° = П ?г и т. д. rW^
98
|
Алгоритм 3.1. |
|
|
|
|
|
|
|
|
и |
1°. Вычислить начальные значения текущих величин по формулам (3.23) |
||||||||
(3.24). |
|
|
|
|
|
|
|
|
|
|
2°. |
Вычислить элементы матрицы || |
||Л^ |
по формуле (3.17). |
|||||
|
3°. |
Выбрать и запомнить пару индексов kt, |
It согласно условию |
||||||
|
|
|
Д^, =шах Д<*> , |
к^М(П |
(3.25) |
||||
|
|
|
|
k t i |
к . I |
k ‘ |
1 = 1 , . . . , |
N . |
|
и |
4°. Пересчитать текущие значения |
величин |
по |
формулам (3.19) —(3.22) |
|||||
(3.18). |
(< := .< + |
1). |
|
t < |
М: |
|
|
|
|
|
5°. Проверить |
условие |
|
|
|
||||
|
|
да — перейти |
к п. |
2°, |
|
|
|
|
|
|
6°. |
нет — перейти |
к п. |
6°. |
|
|
|
|
|
|
Вычислить значение функции F (у0) по формуле (3.1). |
||||||||
|
7°. |
Отпечатать |
результат (F (у0). (% к)м>’ |
прекратить вычисления. |
Полученное решение можно попытаться улучшить. Одним из воз можных способов является последовательное переназначение каждого из средств на N — 1 других целей с оценкой ожидаемого при этом эф фекта ДF. Следует оговориться, что невозможность улучшить решение таким способом не является достаточным условием оптимальности,
о чем свидетельствует следующий пример:
Вц |
1 |
1 |
|
|
|
|
|
Г= 1 0 . 8 |
0 . 6 |
0 |
1 |
I УriI • |
I I |
(3-26) |
|
|
0.5 |
\Ргі\\ |
|
|
|||
2 |
0.7 |
1 |
0 |
|
|
|
|
|
/ = |
1 2 |
|
/ |
|
|
|
Решение, заданное |
матрицей || yrj |, |
явно не оптимально и не мо |
жет быть улучшено указанным способом. Однако в силу определенной специфики алг. 3.1 решение, если оно таит в себе погрешность, как правило, удается улучшить указанным способом. Получим необ ходимые формальные выражения для вычисления эффекта от перена значения.
Изменение значения целевой функции (3.1) будет определяться увеличением ее г(/)-й составляющей (за счет снятия Іг-го средства с /-й цели из группы nt) на величину А/7* и уменьшением £ (/)-й составля ющей на величину АF^) (за счет закрепления /г-го средства за /-й це лью, входящей в группу л;). Указанные приращения на основании (3.1) можно получить в виде:
|
1 —« W - L «г |
|
АFm) = А { 1 — |
Чы П (і-®и |
|
|
ц= с. |
|
■А, |
1- П (1 -® Н |
(3.27) |
L |
п='г |
|
4 * |
|
99 |