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

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

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

 

- A F W ) = A

1

 

пи

 

 

о

 

 

 

 

 

 

K=e£

 

 

 

 

- A t

i -

П

 

k = 1,..., M,

l Ф j.

После элементарных преобразований окончательно получим:

 

 

 

 

 

лШ) Ш(ІИ) п

д ( м > а (м> Du ■

AF(k, fl= A F /(0 - A F C(/) = - i —

 

 

----- (3.28)

 

где

 

 

 

( 1 - ® П * «

i - 4 M)

 

 

 

4

 

 

 

 

 

 

 

 

 

 

4"l = 4

П (i

 

<3'29>

 

4« =4 П (1 -<»<?>),

 

 

 

 

ll= cj

 

 

 

(<=С;

 

 

Средство r =

k

переназначается

c

/-й цели на /

= l'-ю соглас­

но условиям:

 

 

 

 

 

 

 

 

 

 

 

4F <*,/’) = min A F (k ,ir,\

/ =

1...... * ,

 

 

 

AF(k,l')<0,

 

J

k = l , . . . , M .

 

Полученное

значение функции

F

уменьшается

на величину

AF (k,

l’). Следует заметить, что

после

окончания

счета по алг. 3.1

(т. е.

t = М)

необходимые значения величин оз]М)

могут быть полу­

чены после выполнения п. 4° [см. (3.19)].

 

 

 

 

3.

Оценка эффективности алгоритма. В настоящее время не пре

ставляется возможным привести достаточно общую оценку точности алгоритма. Мы располагаем лишь качественными доводами и анало­ гиями, на которые обращали внимание читателя и которые говорят о логической целесообразности алгоритма. В дальнейшем, после зна­ комства с методом последовательных приращений, в § 4.1(5) нам пред­ ставится возможность получить метод точного решения данной зада­ чи для случая равноэффективных средств и сравнить его с алг. 3.1. Это позволит повысить уверенность в хорошей точности алг. 3.1. Однако, несмотря на большую по сравнению с предыдущими зада­ чами ограниченность в получении оценок погрешности, мы и сейчас имеем возможность рассмотреть один частный случай, который дает некоторое представление о точности алг. 3.1 или, по крайней мере, о порядке возможной погрешности. Этот случай соответствует ус­ ловиям (на каждый объект по одной цели):

л ,= 1, і = 1,.... S — N.

(3.31)

Индексы і и j при этом совпадают,

а целевая функция (3.1) прини­

мает вид:

м

 

s

 

f ( y ) =

г*=1

(3.32)

(■= 1

 

Для расчетов при j = I формула (3.28) непригодна (и не нужна).

100

Если обозначим Лгсог = A'h то можно заметить, что задача § 3.1 вырождается в задачу § 2.1, эффективно решаемую алг. 2.1. Остается сравнить решение, полученное на его основе, с решением, получаемым согласно алг. 3.1. Оценки будем вести при том предположении, что все средства равноэффективны и, следовательно, точное решение за­ дачи в этом случае обеспечивается методом МЭ, когда:

Аt(y+) = ÄiPi, і = 1,

(3.33)

Для оценки погрешности применим подход, использованный в § 2.1 (3). Оценим точность решения, обеспечиваемую каждой из формул (3.11) и (3.14). При условии (3.31) из (3.11) и (3.14) соответственно следует:

Д г- Л гсог = Л/,

г = 1,..., 5 ,

(3.34)

Bi = Ai a>t ( l - q f l) = A ; ( l - p * 1), i = l , . . . , S .

(3.35)

Подставляя (3.34) в формулу (3.17), видим, что оценка погреш­ ности в этом случае сводится к оценке среднестатистической погреш­ ности алг. 2.1 [см. § 2.1(3)] и, следовательно, она не превышает 3% от оптимального решения.

Для оценки погрешности при уточненном виде весовой функции

[см. (3.14)] подставим (3.35) в (3.17):

А; (Vo) ~ A'i (\ q¥)p.( \ -\-qf~1).

(3.36)

Условия сбоя по аналогии с (2.29) запишутся в виде системы не­ равенств:

[ Аг (у+) > А, (у+\,

іАг (у о )< А, (у0).

Сучетом формул (3.33) и (3.36) после элементарных преобразова­ ний условия (3.37) приводятся к виду, удобному для анализа:

Рг>РіА,

( l - 9 f ) ( l + 9f - ‘)

(3.38)

Р г <Р і А

где, как и прежде, через А обозначено А\ІА'Г.

Перемножая диагональные члены неравенства, видим, что сбои

возможны при условии

 

 

 

1 - q f

i + q f-

1 > 1.

(3.39)

і-Чг

i + ' ? " -

1

 

Поскольку каждая из дробей а и b может быть больше или мень­ ше единицы (в зависимости от соотношения значений величин рхи pr), то необходимо оценить погрешность в случаях, когда рТ< і р 1 и рг І> > Рі. Полагая, что величины рхи ртслучайны и равномерно распреде-

101

лены на отрезке 0 1 , указанные случаи можно считать равновероят­

ными*. Суммарная погрешность тогда будет равна

ßx = 0,5(ß(1>-f ß<2)).

(3.40)

где ß<’> и ß<2) — относительные погрешности в первом и втором слу­ чаях.

1. С л у ч а й P r <P i

{qr >

Рі)-

 

При любых значениях М ^ 1

имеем:

м

 

 

(3.41)

>

1 ,

 

 

 

1

1 + 9^ -*

Из (3.41) видно, что максимальная область сбоя (3.39) будет иметь место при Л4 = 1 , когда

(М 2 )

При М = 1 условия сбоя (3.38) запишутся в виде:

Рг >

Рі Л,

или j Pr>PiA, __

(3.43)

 

— р,

Рг <

Pi

\ P r < P l V Ä.

 

 

Pr

 

 

_ Из (3.43) видно, что при pr pt сбои возможны только в случае Л < 1. На рис. 5 показана область совместного выполнения условий

* В §2.1(3) рассматривался только один случай, так как во втором слу­ чае (рг < рі) погрешности не могло быть.

102

(3.43). Вероятность сбоя, рассматриваемая как отношение площади треугольника ОБЕ к площади OLC, будет равна

 

 

РА= ( \ - Ѵ Т ) Ѵ %

 

(3.44)

Возникающая при этом относительная погрешность как функция

от А запишется:

 

 

 

 

у,.,=

К рг7 К р,

= 1 _ ЛА = 1 _

м

=

(3.45)

 

А г р г

р г

Ѵ а + а

1 + Ѵ а

Опуская

промежуточные преобразования

[см.

§ 2.1(3)], можно

на основании (3.44) и (3.45) записать выражение среднестатистической

относительной погрешности для случая рг <

рг:

 

 

 

 

р ' у

 

I + V a

 

(3.46)

 

 

 

 

 

 

 

Из (3.46) видно, что в области А=

0,04 -у- 0,06 погрешность

значительна (ß < *)> 1 0 %), поэтому

необходимо уточнить

оценку по­

грешности по параметру А .

 

 

 

Положим,

что параметр А случаен и равномерно распределен на

отрезке О т 1.

Усредненная по параметру А оценка среднестатисти­

ческой погрешности будет равна

 

 

 

 

 

 

( 1 - Y ä T V ä dA= 81n 2 — - = 0 ,0 5 .

(3.47) -

 

 

о

I + V ä

 

2

 

 

 

 

 

 

 

 

Такое усреднение несложно проделать и при любом другом зако­

не распределения А.

 

 

 

 

2. С л у ч а й P r > pl(qr< qi).

 

 

 

При этом условии имеем [см. (3.39)]:

 

 

 

 

 

м '

 

 

 

 

 

 

1 ~ Я ‘І

< 1,

ь

 

(3.48)

 

 

 

1 - ч ?

 

 

 

 

 

 

 

 

Так как нас интересует верхняя граница погрешности, то, прини­

мая а =

1, на основании (3.38) получим условия сбоя алг.

3.1 в обла­

сти рг >

р

 

 

 

 

 

 

 

 

 

 

Рг>РіА,

м—1

 

(3.49)

 

 

 

 

 

 

 

 

 

 

 

Pr<PlÄ 1 + Яі

 

(3.50)

 

 

 

 

1+q™-1

 

 

* Отношение pi!pT заменено его средним значением в области сбоя, т. с. значением на линии OR треугольника ОБЕ (см. рис. 5).

103

Сравнивая (3.49), (3.50) с условиями (2.29), можно заметить, что они различаются только обозначениями. Это дает возможность не по­ вторять все выкладки, а ограничиться ссылкой на § 2 . 1 (3 ) и замеча­

нием о том, что полученную там погрешность следует удвоить учи­ тывая условие рт> рг. Тогда получим ß<2) « 0,06.

На основании рассмотренных двух случаев по формуле (3.40) мож­

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

ßs% = 0,5(0,05 + 0,06)

100«6.

(3.51)

Вероятность сбоев,

обусловливающих

погрешность,

минимальна

в первой половине

процесса распределения, когда

количество

средств М велико [см.

(3.38)]. Здесь даже можно использовать алг.

у+, так как Ahl ж Аti- Сбои в конце процесса меньше влияют на точ­

ность, так как в конце процесса распределяются оставшиеся менее эффективные средства.

Рассмотренный частный случай (цг = 1 ) интересен тем, что нет

необходимости учитывать возможное влияние

ситуации,

связанной

С^ СЛ0Вием

а значит> и усложнение формулы для весовых ко­

эффициентов до вида (3.14)

может вести только к увеличению возмож­

ной^ погрешности.

Именно

поэтому формула (3.11)

обеспечила луч­

ший результат (ß<2) « 3%). Однако в Общем случае

(tii ^

1 ) погреш­

ность, возникающая в связи с условием (3.13),

будет меньше при ис­

пользовании формулы (3.14), поэтому ее применение будет вполне оправданным, тем более что это не влечет за собой существенного увеличения требуемого объема вычислений.

В заключение следует отметить, что проведенную оценку погреш­ ности нельзя считать достаточно общей, так как условие п- — 1 ,

удобное с точки зрения оценки влияния фактора (3.13), игнорирует влияние факта нестрогой эквивалентности функции (2.149).

Представление о требуемом объеме вычислений можно получить на основании ориентировочных оценок числа элементарных операций,

получаемых из анализа вычислительной схемы алг. 3 .1 :

 

умножений да M2N2

+ MN,

 

сложений л? M2N +

MN, ’

/ 3 5 2 )

сравнений ä ; M2N +

MN.

V ‘

Использование формул (3.11) или (3.14) относительно мало влияет на объем вычислений, во всяком случае эта разница в объемах вычислении определяется с помощью коэффициентов пропорциональности при (3.52), но не показателями степеней при М и N. Некоторое увели­ чение объема вычислений (добавочные слагаемые в (3.52)) получается при попытке улучшить решение согласно условиям (3.30).

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

104

Чтобы уменьшить отрицательное влияние на точность решения накладываемых на переменные дополнительных ограничений, эти ограничения должны учитываться при вычислении и пересчете весо­ вых коэффициентов, в частности величин b^- Полученное решение должно быть проанализировано и по возможности улучшено извест­ ными методами [см. (3.30), § 2.1(4) п. г)].

5. Примеры расчета. В табл. 13 и 14 приведены исходные условия для двух примеров. Каждый из них решен двумя способами, с при­ менением весовых коэффициентов, рассчитанных по формулам (3.11)

и (3.14). Для первого примера оба способа дали одну и ту же последо­ вательность распределения средств:

(kt, l t) = (3; 1 -> 1; 2 -> 5 ; 6 -> 4 ; 3 -^ 2; 3)

(3.53)

и значение целевой функции F (у) = 140,0. В табл. 13 решение выделе­ но жирным шрифтом. Полученное решение согласно (3.30) улучшить не удалось.

 

 

 

М атрица

ИP r j И № 1

Т а б л и ц а

13

 

 

 

 

 

А

 

100

 

 

80

50

 

І

I

2

3

4

5

6

7

<°/

1, 0

0,8

0 , 6

0,8

0 , 3

0,6

0,2

г=1

0,1

0 ,7

0,1

0,1

0,7

0,1

0,8

2

0,4

0,6

0 ,5

0,1

0,5

0,9

0,5

3

0 ,7

0,2

0,6

0,2

0,8

0,4

0,2

4

0,1

0,6

0 ,4

0,2

0,5

0,1

0,6

5

0,2

0,4

0,5

0,1

0,7

0 ,8

0,2

 

 

 

М атрица

|| р г} || №

2

Т а б л и ц а 14

 

 

 

 

 

А

 

100

 

 

80

50

 

І

1

2.

3

4

5

6

7

“/

1, 0

0 , 8

0 , 6

0 , 8

0 , 3

0 , 6

0 , 2

Т = 1

0,1

0,1

0,1

0,1

0,7

0,1

0 ,8

2

0,4

0,1

0,5

0,1

0,5

0 ,9

0,5

3

0 ,7

0,1

0,6

0,2

0,2

0,4

0,2

4

0,1

0,1

0,4

0 ,2

0,5

0,1

0,6

5

0,2

0,1

0,5

0,1

0 ,7

0,8

0,2

105

шах

106

Продолжение табл. 15

5

СО

©

см

00

CM

oo

cg•—

СО

см

©

ю

со

Гм

со

о

мН

©

мН

чр

в

 

 

1

1

1

 

 

 

о

00

 

 

см

со

 

 

о

оо

 

О

 

 

 

 

со

о

 

 

О

СО

 

 

О

СМ

 

ю

Я

СО

 

 

о

Гм

 

О

 

 

•ч

•«*

о

о

00

©

ч ;

*

о

©

**м

 

со

о

со

 

СО

©

 

 

о

t-м

оо о

о00 СО

 

о

О

 

 

’~н

м.

о

00

со

О)

 

о

со

см

05

 

 

 

05

CM

Tp

 

см

со

о

о

оо

in

МН

00

см

о

Ч*

©

©

©

©

СО

о

г-

о

мН

оо

CM

©

мН

 

©

00

 

 

4f

^p

 

 

05

©

 

 

©

 

 

©

©

05

 

 

о

©

©

 

о

о

со

СО

05

©

CM

о

о

мН

о

см

 

мН

©

©

«—н

мН

Гм

о

см

СО

00

©

мН

со

со

Г"

 

см

со

 

 

о

о

со

 

 

о

со

 

 

 

со

Л

 

 

 

 

 

 

г -

Г -

 

00

чр

CM

©

CM

о

о о

 

8

о

CM

 

 

©

4P

СП

об

05

о

ю

со

I'M

©

CO

00

05

00

 

 

г -

о

 

 

©

со

©

 

 

со

©

 

 

о

СО

о

 

.

Гм

со

 

 

4 f

 

 

оо

Г -

CM

4p

©

«-.-*

СО

 

4t<

см

\ .у

©

CM

©

со

©

 

in

©

Гм

©

Cm.

со

о

со

 

 

00

Гм

 

 

со

г -

со

 

 

ю

CM

 

 

 

СО

ю

 

 

о

©

 

 

о

 

см

о

о

см

00

©

CO

 

' '

V- >

 

СО

©

CM

 

 

 

о

со

мН

 

 

о

мН

СО

о

M

©

CO

СО

чма

 

 

со

00

 

 

 

со”

 

 

со

чр

 

 

со

со

со

 

 

см

Гм.

 

 

00

00

оо

о

о

г -

©

8

Гм.

о

о

о

оо

о

©

CO

см

см

см

о

оо

см

Гм.

©

4p

05

О)

О)

 

 

ІП

00

 

 

 

1

 

 

©

 

 

05

05

 

 

о

тр

 

 

 

о

см

о

8

©

©

 

4P

о

о

©

©

оо

оо

тр'

о

со

чр

Чр

мН

СО

СО

 

 

in

 

 

о

о

 

 

©

 

 

CO

4p

4P

©

©©

©

.

©

in

©

CM

 

©

oo

 

©

Гм

 

 

4 *

 

©

©4P

чр © мН ©

O0

©

CM

©

CO

©

Гм

©

CO

CO

©©

©©

©CO

©©

0-1

CN

 

© со

in -g* •m*

©

3

Wn

 

3 ОЭ

3 CQ

3

щ

-£й

 

 

 

 

 

-

I иg л е т и - е

I == "9 лвш и-t-

1

 

 

 

13?

Второй пример по исходным данным отличается от первого толь­ ко тем, что все элементы матрицы Цр^-Ц для второго столбца умень­ шены до 0,1. Этим было усилено влияние фактора, отмеченного усло­ вием (3.13).

Решения, полученные на основании формул для весовых коэффи­ циентов (3.11), (3.14), и соответствующие им значения целевой функ­

ции следующие:

 

 

 

 

 

 

 

 

 

 

(3.11):(3; 1-J-2;

6 -^ 5;

3 -^ 1;

5-^ 4;

4),

^ =

154,4,

 

 

(3.14): (2;

6 -+ 3 ;

1-+1;

7 -^ 4;

4 ^ 5 ;

5),

F2=

153,2.

(3.54)

В табл. 14 первое решение выделено прямоугольниками, второе—

жирным шрифтом. Вычислительная

схема одинакова для обоих слу­

чаев, поэтому в табл. 15 мы ее приводим только для основного

(второ­

го) метода расчета. Элементы табл.

15

заполняются

по такому же

принципу, как и в табл.

1 1 .

 

 

 

 

 

 

 

В

последней

строке табл.

15

записаны

компоненты

вектора

{соуЛ1)},

на основании которых можно попытаться улучшить решение.

Согласно формуле (3.28)

, имеем для k =

3 (/ =

1 , Г =

4).

 

 

АF ■

5,6-0,3-0,7

26,1-0,64-0,2

к „

n 0

 

(3.55)

 

0,7-0,3

0,36

- = 5,6 — 9 ,2 = — 3,6.

 

 

 

 

 

 

 

 

Переназначение средства k = 3 с первой на четвертую цель дает

дополнительный

эффект —3,6 ед. Полученное

решение оптимально

и ему соответствует значение целевой функции F (у0) =

149,6.

Прибли­

женное решение имеет погрешность 2,4% и укладывается в пределы полученной оценки погрешности. Решение с помощью формулы (3.11), дало, как и следовало ожидать, худший результат (ß2 == 4,55%).

Из табл. 15 для второго шага процесса, на котором происходит неоптимальное назначение средства k = 3 , видно, что отношение зна­

чений соответствующих элементов матрицы ||AftZ|| равно, примерно, двум, поэтому если предположить, что изменение погрешности про­ порционально изменению элементов Ahl, то в момент выполнения ра­ венства А 31 = Аз* погрешность не превысила бы 5%.

Рассмотренные примеры являются частным подтверждением спра­ ведливости полученных представлений о точности алг. 3 .1 .

§ 3.2. Оптимизация системы с усложненной двухиндексной целевой функцией

аддитивно-мультипликативного вида

при независимых элементах системы

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

Требуется найти матрицу ||у°/||л<лг, доставляющую минимум функции вида:

F ( y )= 2

V

fW))

(3-56)

i = 1

r/= = 1Л ' J

 

* При Дзі <С Д34 решение становится точным.

108

при ограничениях на переменные

 

N

 

 

 

(3.57)

 

2 j Y t f = l .

r = l , . . . ,

Af,

 

 

/ = 1

 

 

 

 

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

 

 

 

 

Yrj 6

{0 ; 1 },

 

1,

М,

 

0 ^

ІЯгі — 1 — P t j X

 

 

U

1 , •••. Л’,

(3.58)

 

 

/ =

о < К - і = i — e j i X i . Л > о ,

Данная формальная запись отличается от задачи § 3.1(1) наличи­ ем второго индекса (г) у величины со7-г. Запись со^г> означала, что соj является функцией от і, так что соj Ф 0 только при условии

сі < / (

0 < “г-

(3.59)

Новая запись (сол) говорит о

независимости

индексов і и /. Ве­

личина (Oji может принимать любые значения на отрезке 0 -f- 1 при

каждом сочетании і и /. Для задачи

§ 3.1(1)

величины

№,-(*) также

могли быть

заданы

в виде

матрицы

|| со;-; Ц^,

однако

эта матрица

тогда должна была бы удовлетворять дополнительным условиям:

 

 

 

/= 1

=

“ /(/),

 

j =

 

(3.60)

 

 

 

 

 

 

 

 

 

или,

точнее:

 

 

 

 

 

 

 

 

 

 

1

0 ,

если

/ < с г, или

j > u u

(3.61)

 

 

 

I сод,-),

если

ci ^s j ^ . u i.

 

 

 

 

 

 

Условие (3.61) отражает

собой факт блочной структуры матри­

цы I

Цл/s

для задачи

§ 3.1(1):

 

 

 

 

 

 

Й!

 

і,

I

 

 

 

 

 

 

0

0

0

0

0

 

 

 

 

щ

 

 

. . .

 

 

 

 

со«с

 

о

 

 

 

 

 

 

 

о

 

со..

 

 

 

 

 

 

 

сг

 

 

ajtbs-

(3-62)

 

 

/

 

®/<о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

соcS

 

 

 

 

0

 

0

 

 

 

 

 

 

 

0

о

о

0

0

соN У; S

 

 

Обобщение задачи влечет за сбой увеличение круга практических

задач и физических интерпретаций, соответствующих новой формаль­ но

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