книги из ГПНТБ / Берзин Е.А. Оптимальное распределение ресурсов и элементы синтеза систем
.pdfВ отношении возможности дальнейших обобщений путем введения раз личного рода дополнительных ограничений на переменные остаются справедли выми замечания, высказанные в §3.1(3).
По сравнению со случаем § 3.2(3) объем вычислений увеличится примерно в S раз (за счет усложнения формулы весовых коэффициентов) и будет характери зоваться следующими оценками числа элементарных операций:
сложений ж S2N2M2,
умножений ж S2N2M2, сравнений ж NM2.
§ 3.4. Оптимизация системы с обобщенной
двухиндексной целевой функцией аддитивно-мультипликативного вида при наличии
связей между элементами системы
1.Постановка задачи. Запишем формальную постановку задачи.
Требуется |
определить |
матрицу ||уЛ] II лгм, |
доставляющую |
мини |
|||||||||||
мум функции вида: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
S |
/ |
S |
п |
|
|
|
N |
|
|
|
|
|||
Р(У)= 2 At 1— |
|
|
■ П |
|
1 ©іх/ X |
|
|||||||||
|
» = 1 |
\ |
/ = 1 |
м |
|
|
ц. = 1 |
|
|
|
|
||||
|
N |
|
1 - |
|
|
|
|
|
|
|
|
|
|||
|
X П |
1 — |
п аугц |
|
|
|
|
|
|
(3.118) |
|||||
|
Т ) = 1 |
|
|
|
г = 1 |
ГЦ |
|
|
|
|
|
|
|
||
при ограничениях на переменные |
|
|
|
|
|
|
|
|
|
||||||
|
N |
|
|
|
|
|
|
|
М, |
|
|
|
(3.119) |
||
|
2 |
Утт) — 1 1 |
|
Г = |
1 , |
|
|
|
|
||||||
|
Л= 1 |
|
|
|
|
|
|
|
|
|
|
|
|
||
и дополнительных условиях: |
|
|
|
|
|
|
|
|
|
|
|
||||
Угц € {0; 1}, |
|
|
|
|
r = |
1, .... |
М, |
|
|||||||
0 |
(Угц — 1 |
Ргц) |
|
1, |
|
|
|||||||||
|
|
/, |
і |
= |
1, |
.>., |
S, |
(3.120) |
|||||||
0 йС (Юц/ = 1 |
Bp,/) sjC 1, |
||||||||||||||
р, |
т)= 1, |
|
..., |
N, |
|
||||||||||
|
|
(“ л = |
*)> |
|
|
||||||||||
|
|
ц, р = |
1, |
..., |
N, |
|
|||||||||
|
0 :;С ßr)n |
1 |
(ßrin— |
1 )> |
|
||||||||||
|
|
Лг> 0, |
|
г= 1, .... |
S. |
|
|
|
|
|
|||||
Сравнивая функции (3.118) и (3.109), можно заметить, что они |
|||||||||||||||
различаются только сомножителем при со^/. Функция |
|
||||||||||||||
|
|
|
|
|
г 6 |
п |
ЯгцѴгц |
|
|
|
|
|
|
(3.121) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
заменена функцией |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
п |
|
1 - |
п |
ГГЦ |
|
JT1H |
|
|
(3.122) |
|||||
|
|
1гц |
|
|
|
|
|||||||||
|
Г) = 1 |
|
|
|
г 6 М(П |
|
|
|
|
|
|
|
120
При условии |
|
|
ßri|j,= 0 для всех г| и [X, |
кроме |
[х = rj (ßtlt) = l), (3.123) |
оба эти выражения совпадают и |
задача |
§ 3.4(1) сводится к задаче |
§ 3.3(4). С функцией вида (3.122) мы уже |
знакомы по § 2.2(1). Вели |
|
чины |ЦЙ являются коэффициентами связи между N целями аналогич |
но тому, как величины ajt характеризуют связи между 5 элементами системы (объектами).
Поскольку задача, по существу, является объединяющим итогом всего предыдущего материала*, в заключение дадим ее обобщенную физическую трактовку.
Имеется система, состоящая из 5 элементов. Составляющие ее элементы имеют веса А t ^ 0 и взаимосвязаны между собой, что коли чественно отражается матрицей коэффициентов связи Ца^Щз.
Каждый элемент системы испытывает внешнее воздействие со стороны N целей. Эффект воздействия р-й цели (р = 1, ..., Л') по /-му элемен ту системы количественно определяется матрицей || toц/1| jvs.
Действия целей на элементы системы независимы, так что суммар ное воздействие, например, двух целей на /-й элемент системы будет равно величине 1 — ( 1 — ®ю) ( 1 — «щ).
Цели взаимосвязаны между собой так, что состояние любой из них зависит от состояния остальных. Эта зависимость задается матри цей коэффициентов связи ЩпиЦлглг-
Имеется М единиц ресурса (средств), которые необходимо так распределить для воздействия по N целям (найти матрицу ||y“ri| atw), чтобы в максимальной степени (с учетом связей между элементами системы и связей между целями) уменьшить отрицательное влияние общего воздействия целей на элементы системы. Степень снижения влияния трй цели на систему при назначении для ее обстрела r-то сред ства задается матрицей Ц^ЛлшСуммарный (главный) показатель (критерий) качества функционирования системы определяется функ
цией F (у) = |
от параметров системы с учетом связей между ними, |
І |
что система функционирует тем лучше, чем ближе |
причем считается, |
|
к нулю значение целевой функции F (у). Если F (у) = 2 А ь то систе- |
|
|
І |
ма не функционирует. |
|
2. Метод и алгоритм решения. Метод и подход к решению остают |
ся прежними, однако при формировании алгоритма для новых усло
вий задачи мы должны теперь ориентироваться на алг. 2 .2 , |
который |
и возьмем за основу. |
|
Центральным становится вопрос о том, каким образом |
следует |
обобщить функцию для весовых коэффициентов (3.110), чтобы она, |
удовлетворяя всем прежним условиям, соответствующим образом учи тывала и новый фактор. В результате анализа задачи оказалось воз
можным влияние нового фактора учесть через замену величины Ь(£}
* В зависимости от тех или иных условий [см., например, (3.123), (3.117), /ij = 1 и др.[ задача вырождается в любую из рассмотренных ранее задач, за исключением трехиндексных.
121
в формуле (3.110) выражением (3.122) и приращения (3.101) прираще нием, вычисляемым по формуле (см. (2 .8 8 )):
&ki - 2 |
ы<- |
|
|
|
|
|
Dß |
|
|
|
Iß |
|
|
u=i |
|
|
|
|
|
|
|
|
N |
|
|
|
|
i - ( i - д {Г Х)а\1- |
Х)ж |
, |
5 |
t |
1)a^ ) ß^ |
(3.124) |
1 (i Q\i~ i) a/<>) ß/n |
|
N |
|
|
||
|
|
|
|
|||
|
|
|
|
|
||
|
|
|
T]=l |
|
|
|
k e M w , |
/ = |
1........N. |
|
|
Формула весовых коэффициентов после указанной замены в фор муле (3.110) примет следующий вид:
|
|
|
|
S |
|
а«) |
S |
|
|
|
|
|
|||
|
|
|
|
V— |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
в\і] = |
|
S |
|
|
2 |
|
|
|
|
|
|||
|
|
|
|
Е = 1 |
2 |
«ЧW == 1 |
|
|
|
|
|
||||
|
|
|
|
/ |
|
п = 1 |
|
|
|
|
|
|
|
|
|
|
|
- |
N |
1- 4 |
7 |
N |
|
|
|
|
|
|
|
||
X |
|
п |
п {1— |
|
Jim) |
ап |
|
||||||||
|
|
|
|
\ |
|
|
|
Т| = |
1 |
|
|
|
|
|
|
|
|
1\(і--[1 —а}0] ал) |
|
1 = 1 ....... |
N. |
(3.125) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сформулируем алгоритм решения задачи. |
|
|
|
|
|||||||||||
Алгоритм 3.4. |
начальные значения |
текущих величин |
по |
формулам: |
|||||||||||
1°. |
Вычислить |
||||||||||||||
“и/) = ®Д/ |
рі, |
Т) = |
1........N, |
|
а^0)= пN |
ей /, |
|
I, |
п = |
|
|||||
|
м |
|
|
|
|
|
|
|
|
и = і |
|
|
|
(3.126) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
< ' = |
П <7ГТ), |
/ = 1 .......s > |
|
ßL0)= |
2 |
«V, |
= |
1, . . . , S, |
|
||||||
|
г—1 |
|
|
|
|
|
|
|
|
И= 1 |
|
|
|
|
|
|
|
|
|
< С = 1 . |
л= |
1 .......N |
|
|
|
|
|||||
ß |0> вычисляются по формуле (3.125) |
при t = |
0. |
|
|
|
|
|
||||||||
2°. Вычислить элементы матрицы ЦД^Ілш по формуле (3.124). |
|
||||||||||||||
3°. Выбрать и запомнить пару |
индексов |
kt, |
It согласно |
условию |
|
||||||||||
|
Д£°г = т а х |
Д ^ , |
к£М<{), |
/ = |
1........ |
N. |
|
|
(3.127) |
||||||
|
* |
f |
k ,I |
|
|
|
|
|
|
|
|
|
|
|
|
122
4°. Пересчитать текущие значения величин по формулам:
|
|
|
|
|
( |
|
|
|
|
|
если г1 у |
/(, |
|
|
|
||
|
п ( 0 |
— 1__а о ( 0 __I |
11 |
|
’ |
|
|
если |
|
|
|
|
(3.128) |
||||
|
|
|
* 4 |
1 |
о Ц - І ) |
( Г|’ |
r\ — lt, |
|
|
|
|||||||
|
|
|
|
|
"т) |
^ |
|
|
|
|
|
|
|
||||
|
|
|
|
1—( і |
—q S'O Рі |
u |
* І = |
1 ..........s , |
|
(3.129) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
p = |
l ........ |
N. |
|
|||
|
w |
|
i - 0 - ^ - , ) ) P / |№ ’ |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||
|
< |
) = |
n |
( ! - < |
; ) , |
/ - 1 , |
|
|
5 . |
|
|
|
(3.130) |
||||
|
|
|
M= 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
... |
N, |
|
|
|
(3.131) |
|
|
|
|
|
ц = |
1 |
|
|
/ = |
і . - |
|
S, t: |
= |
t - \ - |
1. |
(3.132) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
5°. Проверить условие і < М: |
нет — перейти |
к |
п. |
6°. |
|
|
|
||||||||||
6°. |
да — перейти |
к |
п. |
2°, |
|
|
или (3.118). |
||||||||||
Вычислить значение целевой функции по формуле (3.116) |
|||||||||||||||||
7°. |
Отпечатать результаты |
(F (у0), |
||ѵ?т]|л<л? или (Мг) |
Для |
t = |
I, ..., М), |
прекратить вычисления.
В заключение может быть предпринята попытка улучшить реше ние путем проверки М (N — 1) дополнительных вариантов, возникаю щих за счет поочередного переназначения каждого из средств.
3. Оценка эффективности алгоритма. Чтобы выяснить влияние нового фактора на точность решения, оценим величину погрешности для двух предельных случаев:
1 ) |
ßm i = 0 |
для всех |
т) |
и |
р, |
кроме т] = р (ßnn = 1 ), |
2 ) |
ßru.1 = |
1 Для всех |
т) |
и |
р. |
|
Уже отмечалось, что формула для весовых коэффициентов в пер вом случае сводится к виду (3.110). В свою очередь, выражение для приращения (3.124) вырождается в формулу (3.101). Это было показано в § 2.2(3). Следовательно, алг. 3.4 и 3.3а совпадут.
Во втором случае для (3.122) будем иметь:
П |
( і — |
( і — |
П |
= П П |
(3.133) |
Г)=1 |
\ |
\ |
Г = 1 |
/ / Г) г |
|
На основании этого для целевой функции (3.118) имеет место экви валентность соотношений (стр. 85):
min/*'(у) ~ |
min П П qj^ ~ m a x |
SPrrpYm. |
(3.134) |
||
V |
V Ч Г |
1 |
у |
Г1 Г |
|
что определяет точный алгоритм для данного случая: назначение каж дого средства в соответствии с максимальным элементом рГТ) матрицы ІІРгпІ (или минимальным элементом матрицы ІдУчІІлш).
* |
Формула |
(3.129) — аналог формулы (2.83), а (3.130) и (3.132) более |
общи, |
чем (3.98) |
и (3.99). |
123
К такому же решению приведет и алг. 3.4, причем в этом предель ном случае значения весовых коэффициентов не оказывают никакого влияния на решение. Это было показано в § 2.2(3).
Таким образом, после введения фактора, учитываемого величи нами ß,,^, точность алг. 3.4 будет, по крайней мере, не ниже, чем у лю бого из предыдущих алгоритмов, рассмотренных в данной главе*.
Объем вычислений характеризуют следующие оценочные величи ны числа элементарных операций:
сложений « |
N3M (S2 + |
М2), |
|
умножений |
N3M (S 2 + |
М2), |
(3.135) |
сравнений » |
NМ2 |
|
|
Некоторого уменьшения объема вычислений (примерно в N раз)
можно добиться за счет перехода от алг. \ 0 к алг. у+. При этом поря док расчетов изменится лишь в той части, что при вычислении элемен тов матрицы I А*} I необходимо будет использовать формулу (3.124) без ее второго отрицательного члена. Большой потери точности при этом ожидать не следует.
§ 3.5. Некоторые нерешенные целочисленные
задачи в игровой постановке
1. Игровая постановка задачи § 3.1(1). Анализируя ранее рассмо тренные задачи, можно заметить, что если, определяя матрицу ||б;-г||, мы максимизировали целевую функцию F (б) (см. § 2.1), то ее миними зацию осуществляли нахождением соответствующей матрицы ||уГц| (см. § 3.1). Можно также заметить, что с ростом величин N и <д7-г зна чение функции F (б) все более увеличивается и стремится к пределу, равному 2 А, но с увеличением величин М и рГ11 (см. (3.1) ) значение
І
целевой функции все более приближается к нулю.
Таким образом, мы имеем две противоположные тенденции изме нения целевой функции F (б, у) в зависимости от параметров б и у, что говорит о наличии конфликтной ситуации и целесообразности пои ска гарантированного решения.
Условимся считать, что существуют две соревнующиеся сторо ны — А и Б. Сторона А управляет параметром (матрицей) у и заинте ресована в минимизации функции F (б, у), а сторона Б управляет параметром б и максимизирует значение целевой функции F (б, у).
Итак, наметилось основное отличие новой постановки задачи от
всех предыдущих: взаимное влияние и необходимость |
учета |
влия |
ния действий одной стороны на действия другой, т. |
е. б0 = |
б (у) |
и Уо = у (б). На формальной постановке задачи это отразится следую щим образом.
* Данный алгоритм в силу своей общности пригоден и для решения задач гл. 1 и 2, однако использовать его там явно невыгодно, так как частные алгоритмы экономнее, а в иных случаях и точнее.
124
Требуется найти |
матрицу |
8 0 |
= |
| | 6 / ) j j , доставляющую |
максимум |
||||
функции min F (б, у), |
т. |
е. |
8 0 |
и у0, |
которые обеспечивают |
|
|||
V |
|
|
|
|
|
|
|
|
|
|
F (б0 Yo) = |
max min F (б, у) = |
|
|
|||||
|
|
|
|
|
б |
|
у |
|
|
|
|
S |
|
[ |
|
N |
/ |
М |
|
= max min |
2 |
А |
1 |
— |
П |
1 — 8Jta>Jt |
П qy/i |
(3.136) |
|
б |
у |
І = 1 |
|
\ |
|
/ = 1 \ |
г — 1 1 |
|
|
при ограничениях |
на переменные |
|
|
|
|
1 8 „ = 1 ,
»= 1
!1
/ = 1 ....... |
N, |
г = 1 , ..., м, (3.137) |
и дополнительных условиях: |
|
|
|
|
..., м, |
|
|
1 , .... |
N, |
Trj С {0 ; |
i>, 1 |
r = \ , |
(3.138) |
||
0 < о ) л < 1 , 1 і = 1 , .... |
1 , |
|
1 |
/ = 1 , |
..., |
м. |
|
|
|
||||||
Следует заметить, что вместо (3.136) может использоваться экви |
|||||||
валентная форма записи целевой функции: |
|
|
|
|
|
||
F Фо, То) = max min 2 |
А |
^ 1 — П |
1 — |
П qjfi'j |
|
. |
|
Физическая интерпретация задачи может состоять в следующем. |
|||||||
Сторона А, управляя обороной (матрица у), |
минимизирует ущерб |
объектам при любом варианте удара стороны Б (матрица б). Требует ся найти такой вариант действий стороны Б (у0), чтобы ущерб был ма ксимален.
Если матрица б фиксирована и не зависит от у, то будем иметь-
|
I |
I |
Ö |
« |
® |
н |
Матрица |
|| со7-г || вырождается |
в матрицу (3.62) и задача фактиче |
||||
ски сводится к задаче § 3.1(1). |
|
|
|
|
||
Если М = |
0 или Pr, — 0 для всех г и /, то будем иметь дело с дру |
|||||
гим предельным случаем, |
который |
формулируется |
в виде |
задачи |
§2. 1( 1).
2.Игровая постановка задачи § 2.5 (1). По аналогии с 3.5(1) за
пишем формальную постановку игровой задачи для |
§ 2.5(1). |
|||||
Требуется найти |
матрицы 8 0 |
и у0, доставляющие значение max |
||||
min функции F (б, у): |
|
|
|
|
|
6 |
у |
|
|
|
|
|
|
S |
( |
L |
N / |
М |
\ è i h \ |
|
F(6,y)= 2 |
Л, 1 - |
п |
П ( і — |
п |
qy ) |
(3.139) |
i= 1 |
\ |
А=1 /= 1 \ |
r= 1 |
' / |
) |
при ограничениях на переменные
тN
f t = l |
= 1, / = |
1, ....N, r = l , ..., м, |
i = l |
|
|
Матрица ||согу-|| |
характеризует |
возможности /-й единицы при |
действии по t'-му объекту.
125
и дополнительных |
условиях: |
|
|
|
|
|
||
fyfte { o ; i } , |
\ |
|
|
|
I |
г = 1 ’ |
М' |
|
0 < ® ;>г< |
1 , / |
....... |
’ |
0 < 7 rJ. < l , J |
j = l , . . . , N , |
(3.140) |
||
|
I— I) • • •>О» |
|
|
|
|
|
||
|
|
Аг> |
0, |
і = 1, |
S. |
|
|
|
При фиксированной матрице Цб^Ц получим: |
|
|
||||||
|
|
|
|
|
(ІФІ(і)), |
|
(3.141) |
|
и задача вырождается в задачу § 3.2(1). |
|
|
|
|||||
В случае М = |
0 (или prj = |
0) для всех г и / задача сводится по |
||||||
своей постановке к задаче § 2.5(1). |
|
|
|
|
||||
На основании |
указанных |
задач |
можно |
представить себе и фи |
зическую интерпретацию данной задачи. Задача существенно услож нится, если величины qrj поставить в зависимость от индекса k (qrih).
3. Игровая постановка задачи § 3.3(1). Требуется найти матриц б0 и у0, доставляющие значение max min функции F (б, у):
бу
S |
. |
|
S ( |
Г |
N |
NI / |
М |
|
|
|
ан |
>(3.142) |
|
F ( 6,ѵ)= |
|
1— п |
1 - |
1 - П I |
|
|
|
|
|||||
і= 1 |
\ |
|
/= 1 \ |
L |
• |
( 1— бціСОці П qjn* |
|
|
|||||
|
іх= I \ |
|
|
|
|
|
|
||||||
при ограничениях на переменные |
|
|
|
|
|
|
|
||||||
|
|
|
|
s |
|
|
N |
|
|
|
|
|
|
|
|
|
і= 1 |
|
|
2 |
7гц — 1 > |
|
|
|
|
(3.143) |
|
|
|
|
|
|
ц = 1 |
|
|
|
|
|
|||
|
|
|
р,== 1, ..., N, |
1, .... м, |
|
|
|
|
|
||||
и дополнительных |
условиях: |
|
|
|
|
|
|
|
|
||||
бцг€{0; |
1}, |
I |
ц = 1 , |
N, |
7 г ц 6 { ° ; 1}. |
1 |
|
г = 1, |
..., M, |
||||
0<соц/<1, |
I |
і= 1 , |
.... S, |
0 < p rtl< l , |
I |
|
И- = 1 .......N, |
||||||
|
А > 0 , |
|
|
|
I . |
..., |
S. |
|
(3.144) |
||||
|
0 < а л < 1 |
(aw=,l), j |
h t |
|
|||||||||
|
|
|
|
|
|
||||||||
При фиксировании |
матрицы Цб^Ц |
мы имеем |
дело |
с задачей |
|||||||||
§ 3.3(1), но при Л4 = |
0 или рГІ1 = |
0 для всех г и р, |
данная задача вы |
||||||||||
рождается в задачу § |
2.2(1). Матрица Ца^Ц вводит |
зависимость меж |
|||||||||||
ду объектами Лг (г = |
I, |
..., |
S). |
|
|
|
|
|
|
|
|||
4. |
Игровая |
постановка задачи |
§ 3.4(1). |
По |
аналогии с предыд |
щими задачами запишем формальную постановку и данной задачи. Требуется определить матрицы б0 и у0, доставляющие значение
max min функции F (б, у):
бY
s |
|
L N / |
|
F(б, 7) 2 |
л ; |
П ГТ і і |
X |
|
|
' - Д М 1- k—1 ц= 1 \ |
|
|
N |
|
|
X |
П |
|
(3.145) |
|
Т) = 1 |
|
|
126
при ограничениях на переменные
т N
|
|
|
2 |
угц |
- |
1 у |
|
(3.146) |
|
|
N, |
Т) = 1 |
|
|
|
||
ц = |
1, |
r = |
1, .... М |
|
|
|||
и дополнительных условиях: |
|
|
|
|
|
|
|
|
1= 1 , •••. а,- |
|
|
|
|
|
|
||
0 < а л < |
1 |
(а})= |
1), |
/, і = |
1, |
|
S, |
(3.147) |
О ^ ßrin ^ |
1 ) |
(ßrin= |
1 )> |
"П»М'= |
|
•••> |
|
|
|
|
|
||||||
|
Лг> |
0 , і = 1 ....... S. |
|
|
|
|
||
Как и ранее, при фиксировании матрицы 6 |
имеем дело с уже рас |
|||||||
смотренной задачей § 3.4(1), |
а в случае М = |
0 или р гц = 0 для всех |
/- и ті задача сводится к постановке § 2.5(1), но при наличии связей меж ду элементами системы. Отдельно такая задача нами не рассматрива лась.
Если величины рг7] поставить в зависимость от индекса k (prilfe), то задача усложнится, хотя физический смысл изменится несущественно.
Матрицей" I РццЦіѵлг вводятся зависимости между атакующими целями.
Заключительное замечание. В гл. 3 рассмотрен ряд задач оптими зации целевых функций специального вида, отвечающих некоторым физическим системам при наличии ограничений типа (3.146). По физи ческому смыслу эти ограничения означают, что каждая из N (или М) единиц разнородного ресурса независимо от его физической природы может быть задействована только один раз. Однако часто для практи ки представляют интерес не только задачи оптимального распределе ния уже имеющегося набора N' (или М ’) некоторых средств, но и зада
чи оптимального |
выбора |
и |
распределения некоторого |
количества |
|
средств из N (или М) возможных их типов, различающихся эффектив |
|||||
ностями и стоимостями при |
общем ограничении на бюджет G. Фор |
||||
мально отличие в постановке такой задачи состоит только |
в |
измене |
|||
нии вида ограничения. Вместо (3.146) будем иметь: |
|
|
|||
2 |
2 |
|
2о ' ц У 2г ц < о ' , |
|
( з . і 4 8 ) |
где Xp,ft, Gllk и г/ГТ1, С/ч — |
количество и стоимость образца |
средств |
р,-го (r-го) типа, снаряжаемых по k-му (т]-му) варианту или назначае мых на k-й (т]-й) объект, цель.
Решение задач с такими ограничениями на переменные требует применения методов оптимизации аддитивных функций более общего вида с заранее нефиксированной величиной приращения аргумента на каждом шаге процесса. Один из таких методов, изложенный в гл. 4 (§ 4.4), в сочетании с ранее рассмотренными, позволит получить реше ние задачи при ограничениях вида (3.148).
Г л а в а ч е т в е р т а я
МЕТОД ПОСЛЕДОВАТЕЛЬНЫХ ПРИРАЩЕНИЙ
(МЕТОД ВЫПУКЛЫХ ОБОЛОЧЕК)
Максимизация достаточно хорошо ведущей себя функции представляется сравнительно простой зада чей...
Однако здесь приходится считаться с некоторыми осложняющими дело обстоятельствами. Прежде всего, эффективное аналитическое решение систем, состоящих из большого числа даже простых уравнений (например, линейных), является весьма нелегким де лом. Омрачая наши перспективы еще более, напомним, что численное решение даже таких систем обычно на талкивается на целый ряд трудностей — как техниче ских, так и принципиальных. Следовательно, когда число переменных велико, определение этого максиму ма отнюдь не является шаблонным.
Все это дает нам право говорить о «проклятии многомерности». И хотя это проклятие много лет тяготело над головами физиков и астрономов, все же не следует терять надежду на получение вопреки ему важных результатов.
Р. Веллман
Однако и получение таких результатов не снимает проблему мно гомерности. Это и понятно, если она возникает, то всегда является препятствием при получении этих результатов. Особо важное зна чение проблема многомерности приобретает тогда, когда приходится иметь дело с оптимизацией систем, состоящих из большого числа эле ментов, а потери, связанные с неоптимальностью построения или функ ционирования системы, могут быть очень велики. В качестве примера можно привести систему народного хозяйства страны или даже отдель ных ее отраслей. Никакая сколь-либо разумная степень централиза ции управления такой системой с применением экономико-математи ческих методов не может не столкнуться с проблемой многомерности. Борьба за размерность — это одна из задач, которая ставится в данной главе.
Другая задача будет состоять в обобщении вида целевой и ограни чивающей функций при сохранении, однако, условия аддитивности. Необходимость такого обобщения на практике определяется сложно стью реальных процессов, протекающих в системах, в связи с чем кри териальные и ограничивающие функции также могут иметь довольно произвольный вид.
* В гл. 4 нам будет удобнее обозначать вектор в виде X в отличие от
X = (xj}n, как это делалось раньше. Теперь X — сумма компонент вектора X •
128
§4.1, Максимизация аддитивной целевой функции произвольного вида при линейном ограничении на переменные и метод ПП
1. Постановка задачи. |
|
Требуется найти |
вектор |
Х0 = {х°і}п, до- |
||
ставляющий максимум аддитивной функции: |
|
|
||||
F( X) = |
і |
Ft (x,) |
|
(4.1) |
||
|
|
|
i = |
1 |
|
|
при линейном ограничении на переменные |
|
|
||||
Х = |
2 |
* ,) |
<Ь, 6 > 0, |
целое |
(4.2) |
|
і |
= |
1 / |
|
|
|
|
при условии их целочисленности и неотрицательности:
|
|
€ {0, 1, .... Ь), |
і = \ , ... , п. |
(4.3) |
||
Если |
хотя |
бы |
одна |
из |
функций |
Ft (хг) не убы |
вает в области о, |
Ь, то |
этого |
условия |
достаточно для замены нера |
венства (4.2) равенством. (Это условие не является необходимым.) Каждая из слагаемых функций Ft (xt) является однозначной функ
цией своего аргумента хг. В случае многозначности берется ее наиболь шее значение, что соответствует требованию максимизации. В отноше нии функций Fi (Хі) будем предполагать, что они проходят через на
чало координат: |
|
|
|
|
|
|
F i i 0 )= 0 , і = |
1, |
.. . , |
п, |
(4.4) |
и каждая из них |
неотрицательна: |
|
|
|
|
F i ( X |
i ) > 0 для х*£{0, |
... , |
b}, |
і = 1 , ... , п. |
(4.5) |
В дальнейшем (§ 4.2) от условия (4.5) откажемся.
Условие (4.4) не сужает общности задачи, что видно из следующих эквивалентных преобразований:
ш а х 2 |
Фі (*і) ~ max 2 (фг (*г)—Фг (0)) = тах 2 Fi (хі)> (4.6) |
|
X * |
к 1 |
X ‘ |
где Ф г(0)|0, |
но Ег(0) = фг(0) — ф,(0) = |
0. |
К рассмотренной формальной схеме могут сводиться многие задачи, встречающиеся на практике и связанные с оптимальным, в не котором смысле, распределением ограниченного количества однородно го ресурса между различными возможными областями его исполь зования. Такие задачи могут встречаться как в технико-экономической (распределение сырья между различными предприятиями; оптималь ное обеспечение запасов; вопросы надежности средств и др.), так и в военной областях (целесообразное применение боевых и обеспечи вающих средств, их размещение ит. д.). Не будем останавливаться на возможных сферах приложения данной задачи и на ее возможных
5 Зак. 1292 |
129 |