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

книги из ГПНТБ / Берзин Е.А. Оптимальное распределение ресурсов и элементы синтеза систем

.pdf
Скачиваний:
34
Добавлен:
24.10.2023
Размер:
11.19 Mб
Скачать

В отношении возможности дальнейших обобщений путем введения раз­ личного рода дополнительных ограничений на переменные остаются справедли­ выми замечания, высказанные в §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

ы<-

 

 

 

 

 

 

 

 

 

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

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