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

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

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

Сведение задачи к постановке § 4.1(1) осуществляется по схеме, аналогичной схеме (4.91):

Уі tziI

(ч Wi)

---------- * 9 ,

При этом следует иметь в виду, что в случае многозначности промежу­ точных функций 2 ; и уі необходимо вычислять значения gt для всех значений 2 ; и у% и из полученных значений gt выбирать максимальное,

если решается задача § 4.2(1), или минимальное для задачи § 4.1(1).

Других особенностей, существенного увеличения

объема счета

или ухудшения точности мы не ожидаем. Объем счета несколько воз­

растет,

если промежуточные функции не могут быть представлены

в явном

виде. Это потребует применения итерационного

процесса

для промежуточных вычислений.

 

 

 

 

 

 

 

5.

Пример

расчета.

Вычислительная схема достаточно ясна и

§ 4.1(5)

и 4.2(3j. Способ

избавления

от

нелинейности

ограничений

также может быть уяснен из § 4.3(2)

и не требует

новых иллюстра­

ций, поэтому

на данном примере мы продемонстрируем,

в основном,

возможность решения «непрерывной» задачи без существенного увели­

чения объема

вычислений по сравнению с решением дискретной.

Найдем вектор Х 0, доставляющий максимум

функции

 

 

 

 

f (x ) = (2x1-

x2)-^-x2 =

F1(x1) ^ F 2(x2)

 

 

(4.102)

при ограничении

 

 

 

 

 

 

 

 

 

 

 

 

2*i +

3 4 = gT (*x)+g2 (х2) < 6,

Xi, х2 > 0 .

 

(4.103)

Устраним нелинейность ограничения. В данном примере это можно сделать

аналитически, путем замены

переменных:

 

 

 

 

 

 

 

 

 

 

 

г — —

 

 

 

 

 

 

 

 

 

 

*1 = } / у

(£і = 2х?),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.104)

 

 

 

*2 = | / у

(& = Зх^).

 

 

 

Эквивалентная задача запишется

в виде:

 

 

 

 

 

 

Ф( D = ^ 2 f £

— y

j + ] / ^

 

= Фі (£і) Ф Фа (£*).

(4.105)

 

 

 

 

 

Яі + §2—6 »

 

 

 

 

(4.106)

 

 

 

 

 

gl, ё 2 > 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Неравенство заменено равенством,

так как <р2 fe) — неубывающая функ­

ция.

 

метод

ПН,

с

помощью ^совмещенной

матрицы

Цо/іМ ^/ііі

Используя

(табл. 31)

получаем первое решение [gi] =

(1;

5}

(Ag'j == 1).

 

 

170

/

 

 

Т а б л и ц а 31

/

 

І

1

2

 

1

0,914/1

0,577/1

2

0,086/1

0,239/1

3

0,184/1

4

0,155/1

5

0,138/1

Поскольку сопряженно-минимальные элементы не равны (t»x = 0,914 ф Ф 0,138), то решение может быть улучшено. Их уравнивания требует условие (4.95). Улучшение производится путем дробления шага Аg; и увеличения компо-

ненты g1 (Ѵ] > ѵ2) на величину, равную половине шага

Чтобы не нарушилось условие (4.105), необходимо на столько же уменьшить компоненту g2.

В табл. 32 приведены текущие значения величин, получаемые в процессе улучшения решения. Оказалось достаточно трех итераций, чтобы получить оп­ тимальное решение. Хорошая сходимость процесса улучшения позволяет не

беспокоиться

о начальном решении. Точность решения будет

ограничиваться

только возможностями вычислительных средств. При наличии

многих перемен­

ных (п > 2)

уравнивание элементов о* =

Ац>уAgt производится

согласно схеме:

 

* { ѵ к = т \ п ѵ і }

(шах ц; = иЛ*

(4.107)

 

І

і

 

т. е. переменная gk уменьшается, а gr увеличивается на величину Аgk.

 

 

 

 

 

 

 

 

 

Т а б л и ц а

32

 

 

1

2

3

4

5

6

1

7

8

9

t

Шаг

At)

 

*</>

4 °

 

Vi/Agz

 

Фі°

фі/>

фг

 

 

*1

 

 

 

1

1,00

1,00

5,00

0,707

1,291 0,914/1,0 0,138/1,0

 

0,914

1,291

2,205J

2

0,50

1,50

4,50

0,866

1,225 0,136/0,5 0,140/0,5

 

0,982

1,225

2,207 '

3

0,25

1,25

4,75

0,791

1,258 0,168/0,25 0,132/0,25]

0,956

1,258

2,214

 

Рассмотренный пример может быть решен и чисто аналитическим

путем. Для этого

необходимо неравенство (4.103)

заменить равенст­

вом. Достаточным условием для такой замены, как отмечалось ранее, является наличие хотя бы одной неубывающей функции [<p2(g)]. Иск­

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

В работе [20] рассмотренный пример решается методом аппрокси­ мации (см. [20], стр. 123).

* Имеется в виду, что функции ф* (g{) Еогнуты (заменены их выпуклыми

кверху оболочками).

171

§ 4.4. Задача синтеза оптимальной системы

при ограничении на ресурс и возможный состав

ееэлементов

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

Требуется определить матрицу Х0 = ЦлщЦ, доставляющую мак­ симум функции вида [см. (3.148)]:

 

t w =

 

 

П B* / j

(4.108)

при линейном ограничении на переменные

 

 

т

S

ан хп ^ Ь

*

 

а

а

(4.109)

 

1=1

і = 1

 

 

 

и дополнительных условиях:

 

 

 

 

Хц € {0 ;

1 ,

 

1

1 .......S

(4.110)

0 < ( е ,,=

1 - с о , , ) < 1 ,

Аі> 0 ,

ап > 0 ,

 

/ =

1 , ..., m.

 

 

 

 

 

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

только видом ограничения (4.109).

По физическому содержанию задача может быть истолкована

как распределение

активных средств для воздействий по объектам

из множества і С Is

с целью обеспечения максимума эффекта (4.108).

Отличие задачи от § 1.3(1) состоит в том, что

предстоит не только

оптимально распределить активные средства,

но и оптимально вы­

брать состав средств из заданной совокупности т типов, не выходя за пределы ограничения по бюджету (4.109). Величина ад в этом слу­ чае рассматривается как стоимость одного образца /-го типа при сна­ ряжении его для воздействия по t-му объекту. Через величины aj и Xj будем обозначать соответственно полные затраты и общее коли­ чество задействованных средств /-го типа, вычисляя их по формулам:

S

 

S

 

 

 

öj* =

ад, Хд, Xj

Хці

j ~ 1,

m.

(4.111)

i= I

(= 1

 

 

 

Состав отобранных средств при этом может быть представлен в

виде вектора X =

{xj}m.

 

 

 

 

Общий ресурс, выделенный

для организации

воздействия по

t-му объекту, обозначим через Ьр

 

 

 

 

т

/

\

t = l,...,m .

(4.112)

 

\ b = ' 2 j b l J,

172

Поскольку целевая функция не убывает, а ограничение линейно, то решение будет выполняться при знаке = в ограничении (4.109)

сточностью, определяемой дискретностью величин аjt.

2.Метод и алгоритм решения. Первым интуитивным движен ем, которое возникает на основе знаний материала гл. I, является попытка решить задачу методом МЭ, предварительно вычислив мат­ рицу I ид I!, где

 

=

‘ 1>

- * S*

(4.113)

11 O j i

a.ji

/=1,

...,m .

 

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

легко понять на следующем примере.

 

 

 

Имеются два (/

=

1, / =

2) средства, которым соответствуют сле­

дующие значения

величин

 

и ан :

 

 

 

для

/ =

1

<ön =

1 ,0 ,

au =

1 0 ,

А = 1 0 0

 

для

/ =

2

со21 =

0,5,

а21=

2,

(4.114)

1

 

Согласно (4.113)

имеем:

 

 

 

 

 

 

 

 

 

ѵп =

10,

 

ѵ21 = 25.

 

Применяя метод МЭ, первую и все последующие выборки элемен­

тов Vji следует делать в пользу второго средства, так как ѵ21 >

ѵІЬ

однако ясно, что абсолютной

(1 0 0 %-ной)

эффективности можно

до­

биться только за счет первого средства, так как неравенство

 

 

 

 

(Л =

1 - 8 ^ ) < 1

(4.115)

сохраняется при любом конечном значении х21, т. е. Ь{. Следовательно, начиная с некоторого критического значения величины ресурса bt > >- bікР, становится выгоднее использовать по 1-му объекту средство первого типа (/ = 1). Формально условие целесообразности замены средства k-то типа средством /-го типа можно записать в виде:

Ліыл

1hi

 

ъЫі

(4.116)

hi

ahi xhi

 

 

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

Аналогичные требования будут предъявляться и к синтезу на­ дежной системы, при этом условие (4.116) можно истолковать следую­ щим образом: резервирование малонадежных, но дешевых элементов

173

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

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

Средство типа k может быть исключено из рассмотрения для возможности его использования по г-му объекту, если

[vhi ^ vji>

[ahi > а р ,

/ = 1, ... , от, j =j= k.

(4.117)

<

и тем более если ^

Ісодг < ( ä ji

lcöfti < (ü ji

 

 

Средство типа k может быть вообще исключено из рассмотрения, если условия (4.117) для него выполняются при всех і.

Таким образом, если известны условия оптимальной замены средств,то может быть определена зависимость изменения оптимального состава средств и максимального уровня достигаемого эффекта от вели­ чины бюджета Ьі, выделяемого для организации воздействия по t-му объекту — фіі). При наличии таких функций исходная задача мо­ жет быть записана как задача максимизации функции

 

 

ф $ =

2

 

Ч>і Фі)

 

 

 

 

 

(4.118)

 

 

 

7= 1

 

 

 

 

 

 

 

путем определения вектора b =

{&;}s, компоненты которого удовлет­

воряют условиям:

 

 

 

 

 

 

 

 

 

 

 

i= 1

bt >

0,

t =

l , ... ,S .

 

 

(4.119)

 

 

 

 

 

 

 

 

 

 

 

Исходная задача сведена к задаче § 4.1(1). В целом порядок ре­

шения следующий.

 

 

 

 

 

 

 

 

 

 

Алгоритм 4.4.

 

 

 

 

 

 

 

 

 

 

1°.

Вычислить элементы матрицы

 

по формуле (4.113).

рассмотрения

2°.

Согласно принципу доминирования

исключить из

«лишние» элементы матрицы Цу ц Ц.

 

бюджет

bi <

b

и

применяя условие

3°.

Последовательно

наращивая

(4.116),

определить зависимости cpі(Ьі) для г'= 1, ..., S.

 

 

 

4°.

Согласно алг. 4.1а определить вектор Ь0 =

{&?},

доставляющий мак­

симум функции (4.118).

 

 

 

 

 

 

 

 

 

 

5°. Записать решение, прекратить вычисления.

 

 

 

 

 

3.

Оценка эффективности алгоритма. Решение

исходной задач

расчленено на три этапа, эффективность решения

каждого

из кото­

рых определит и общую эффективность решения задачи.

 

 

Первый этап (проверка на

доминирование)

состоит

в

последо­

вательном сравнении каждого (/-го) элемента t'-ro

столбца

матрицы

Vji |, начиная сверху,

с каждым из нижележащих элементов. Одно­

временно такая же процедура проделывается с элементами t-го столбца

матрицы I (£>л|, чего требуют условия

(4.117). Число операций сравне-

174

нйя для проверки всех элементов одного столбца в худшем случае (если не будет «отсеивания» элементов) составит величину

(/п— 1) + (/п— 2) + - + 1 = <Я(” ~ 1)-

(4.120)

Если учесть, что подобная процедура проделывается и для мат­ рицы И(itjj |, а в каждой из них по п столбцов, то общее число опера­ ций оценится сверху величиной

т(т— 1)л да « т 2.

(4.121)

При больших размерах матрицы такая проверка может стать бо­ лее обременительной, чем в случае ее игнорирования*, поэтому п. 2° алг. 4.4 может быть опущен. Достаточно очевидно, что отсеивание явно невыгодных элементов не сможет повлиять на точность решения задачи.

Второй этап (п. 3° алг. 4.4) расчленен на S циклов, в каждом из которых решается одна и та же промежуточная задача с размер­ ностью, в S раз меньшей, чем размерность исходной задачи. Какихлибо трудностей или неожиданностей на этом этапе ожидать, по-ви­ димому, не следует.

Третий этап (п. 4° алг. 4.4 — собственно метод ПП) подробно рассмотрен в § 4.1. Поскольку характер функций ф;(6 г) известен

(выпуклы кверху), то объем вычислений при составлении совмещенной

матрицы будет несколько меньше, чем в общем случае.

 

 

4. Некоторые обобщения. Так

как

формально

задача сведена

к^условиям^ задачи (4.1), (4.3),

тол все

обобщения, л рассмотренные

в^§ 4.1(4),

остаются справедливыми и здесь. Некоторую опасность

смысле

появления погрешности)

может представлять ограничение

на

имеющееся количество средств каждого типа (х}

х"р), и осо­

бенно если

 

 

 

 

 

 

 

^ € { 0 ; 1},

/ = 1,

.... т.

(4.122)

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

V т а хап > Ь .

(4.123)

/= 1 ‘ £ 's

Впротивном случае ограничение по ресурсу (4.109) становится несущественным, и задача вырождается в задачу § 2 . 1 (1 ).

££■ Однако если условие (4.123) и соблюдено, то все же при ограни­ чении (4.122) придется отойти от алг. 4.4., так как кривые срг(&г) в си­ лу жесткой ограниченности количества средств каждого типа не

могут быть определены заранее. В этом случае более удобным ока­

* Подробнее о подходе к оценке целесообразности проверки на доминиро­ вание см. в приложении III.

175

жется алг. 2 .1 ,

если приращение Д « = А*/ — Д й вычислять

не в

расчете на единицу средств, а в расчете на единицу стоимости'.

 

Да/

= Д&; •

k = 1, . . т, I= 1, ..., 5.

(4.124)

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

ных

5. Пример расчета.

Решим задачу §4.4(1) при исходных условиях, задан­

матрицей эффективности Цсо^Ц, матрицей стоимостей

|| ayt || и вектором

{Л;}.

Определим такой

набор средств поражения ||лу-/||,

не превышающий

по стоимости b ед. (b пока не фиксировано), чтобы при оптимальном их распре­

делении по объектам воздействия

Л;

(/' =

1,

...,

3) обеспечивался наибольший

ущерб этим

объектам:

 

{Лг} = {100;

80;

40)3,

 

 

 

 

 

 

 

 

 

і

 

 

 

і

 

 

 

 

і

 

1,0

0,2

0,2

 

80

5

5

 

 

1,25

3,20

1,60

0,8

0,9

0,6

 

40

40

8

 

/

2,00

1,80

3,00

0,5

0,7

0,9

 

І 20

20

16

 

2,50

2,10

2,24

0,2

0,8

0,4

4,3

20

45

8

4,3

 

1,00

1,42

2,00

 

 

 

 

 

H

 

i

l

 

f

,

Согласно

(4.113)

вычислены

элементы

матрицы

||t»^||. В ходе

первого этапа

решения

из рассмотрения

выбывают

элементы

матрицы

||у/;||, выделенные

жирным шрифтом, так как для них выполняются условия (4.117).

 

Результаты

второго этапа

расчетов для точного решения

представлены

в табл. 33,

в соответствии с которой на рис.

15 построены кривые ф* (Ьр.

Пунк-

 

 

 

 

 

 

 

 

 

Т а б л и ц а 33

Вариант

 

г =

1

 

 

 

і —2

 

 

i = 3

 

Ьі

Фі

{/}*

b2

 

Фг

{/}

 

Фз

{/}

 

 

t>3

1

20

50

3

 

5

 

16,0

1

8

24,0

2

2

40

80

2

 

20

 

56,0

3

16

36,0

3

3

80

100

1

 

40

 

72 ,8

3— 3

32

39,6

3 — 3

4

 

 

 

 

50

 

7 5 ,5

3— 3— 1

48

39,9

3 - 3 - 3

5

 

 

 

 

80

 

79,6

3 X 4

 

 

 

 

 

 

 

 

jXXj 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

120

79,96

3 X 6

 

 

 

* {/}—множество

номеров оптимального

набора средств,

соответствующего

бюджету Ьі.

 

 

 

 

 

 

 

 

 

 

176

тиром

показаны

кривые

ф* (bj)',

полученные для приближенного

решения

(см.

начало

§ 4.4(2)).

 

 

что средство

/ =

2 по

второму

(г =

На

втором

этапе расчетов выяснилось,

2)

объекту

никогда

не использовалось,

несмотря на

то,

что

величина

(сі)22

=

0,9)

максимальна

во втором

столбце

матрицы ||(<)7-;||.

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

показал,

что средство j — 2 (г = 2) не может быть использовано для получения

оптимальных решений, поскольку две единицы средства третьего типа обеспечи-

5 0

О

 

 

%

 

 

5 0

 

 

О

50

Ьі

 

 

Рис. 15.

 

вают величину ш = 1 — е§2, ббльшую чем одно средство / = 2 = 2), при рав­ ных затратах. Это позволяет обобщить условия доминирования (4.117), приме­ няя их к комбинациям средств. Однако это и увеличивает объем счета на первом этапе, что будет оправданным, если увеличение счета компенсируется со­ кращением вычислений на втором этапе.

9

Точные уіешения

200

Приближенные

решения

100

t. ■

 

 

 

 

0

50:

100

/50

6

 

Рис.

16.

 

 

Из табл. 33 видно, что по мере увеличения ресурса изменяется и качест­ венный состав средств (см. третий столбец — {/}).

Решение третьего этапа представлено табл. 34* (точное решение) и 35 (при­ ближенное решение). В соответствии с этими таблицами на рис. 16 построены кривые решений (зависимости максимальных значений целевой функции от общего ресурса 6).

* Табл. 34 содержит элементы совмещенной матрицы, полученной на осно­ вании табл. 33, например, с21 = (80—50): (40—20)= 1,5 = Дф2і-Д62і.

То же относится к табл. 35.

177

k

1

2

3

4

5

 

Т а б л и ц а 34

 

І

 

1

2

3

 

Vki/Abki

 

2,50

3,20

3,00

20

5

8

 

 

1

1,50

2,66

1,50

20

15

8

 

 

2

0,50

0,84

0,21

40

20

16

 

 

3

 

0,27

0,01

 

10

32

 

 

4

0,14

 

30

5

 

 

6

Т а б л и ц а 35

І

1

2

3

3

/(«')

 

>

2

 

2,50

3,20

3,00

20

5

a

1,25

2,56

1,20

20

5

8

0,62

1,85

0,48

20

10

8

0,31

1,18

0,19

20

10

8

0,16

0,75

0,08

20

10

8

 

0,48

 

 

10

 

Выделенным в табл. 33 вариантам соответствуют следующие решения для

ЪX 136 ед.:

6 =

136:

{6г} = {80; 40; 16},

(р = 209,

 

6' =

134:

і} = {60; 50; 24},

q/ = 193.

(4.126)

Распределение и состав средств представлены матрицами (см. также табл. 33 — жирный шрифт):

1

0

0

 

 

I0

10

0

 

0

0

0

,

I! Хц И=

0

0

2

(4.127)

0

2

1

1

1 J*4

3

0

0

 

0

0

0

 

 

0

0

0

 

(X0) =209

 

ф (X) = 193

 

Из рис. 16 видно, что максимальные расхождения между точным и при­ ближенным решениями для условий данной задачи невелики и практически не превышают 8—9%. Однако при определении минимальных затрат для обес­ печения требуемого значения целевой функции (обратная задача) относительная погрешность существенно увеличивается. Так, например, для фзад = 204 ед.

178

необходимо затратить в первом случае 124 ед., а во втором — 200 ед., что состав­ ляет погрешность в определении затрат 61% (всегда в сторону завышения). Уменьшение требуемого значения функции ф3 всего на 2% снижает погрешность почти в 2 раза (до 33%).

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

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

Заключительное замечание. Метод последовательных приращений явился логическим следствием и обобщением известных по отече­ ственной и зарубежной специальной литературе методов решения задачи, записанной в § 4.1(1). Задача §.4.1 может быть решена, в част­ ности, методом динамического программирования. При прочих равных условиях объем счета будет примерно в b раз больше, чем у метода ПП (шаг дискретности одинаков и равен единице).

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

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