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

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

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

Из (2.30) следует, что при N = 1 область сбоя отсутствует, так как А; = Аг+, что видно, например, из (2.29).

При N -*■оо алг. 60 и б+ также совпадают и дают точное решение. Наибольшая вероятность сбоя будет иметь место при N = 2, что мы и примем для каждого шага процесса, снова завышая оценку

погрешности. Тогда условие сбоя запишется так:

(2.31)

 

 

 

 

Рис.

2.

 

 

 

Если

обозначим

А =

Аі :

Аг и вспомним, что ег =

1 — и*

то

условие

сбоя после

некоторых

преобразований можно привести

к

виду:

 

 

 

 

 

 

 

 

Лсо/С](лг < 1

У I

— Лсог (2— о»;).

(2.32)

 

Для сбоя необходимы условия:

 

 

 

 

с о ; > 1 —

 

Y

и

(2.33)

Это следует из условия неотрицательности подкоренных выра­ жений (2.32) и (2.33).

Далее снова будем считать, что величины сог являются случай­ ными, равномерно распределенными на отрезке [0, 1].

На рис. 2 показаны области выполнения условий сбоя (2.31).

Так как А ^ 1 , а из (2.29) следует, что сог : <ог> А, то сбои возможны только в области со,.^сог.

При произвольном значении А ^ 1 вероятность сбоя Р (А) опре­ деляется как вероятность попадания случайной точки в область тре-

50

угольника ОВС. Р{А) равна отношению площадей S x = SobcuS = 1 . Исходя из геометрических соображений и зная координаты точек

В и С [определяются условием (2.32)], можно записать:

-

( 2 - 3 4 )

Вводя новое обозначение а (А)

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

Р (Л) = — ==

.

(2.36)

iS

2

 

Если сбой произошел, то это влечет за собой погрешность

= Аг+ - д Г = д ; (' 1 -

Л ^ W Д+ ß (Л),

(2.37)

V

сог ,/

 

где выражение в скобках — относительная погрешность ß (Л). Отношение параметров со / : мг, необходимое для определения

относительной погрешности ß (Л), примем равным его среднему зна­ чению в области сбоя (рис. 2):

Щ; _

1

■+ 1

 

=

-^-[2 — « ( ! + « ) ] •

(2.38)

со,.

2

 

 

 

 

 

 

После подстановки

и упрощений получим:

 

 

 

ß(3) = -% - =

1— Л - ^ = — - —

(2.39)

 

 

ѵ 2

Д7

сог

2(1 + а)

 

Суммарная погрешность ДF (80) равна сумме математических ожиданий погрешностей за N шагов процесса*. Учитывая (2.37),полу­ чаем:

Д/7(бо)= 2 Zt P(Ä)=p(Ä)P(Ä) 2 Аr\t).

(2-40>

t = 1

<=1

 

П р и м е ч а н и е . Следует заметить,

что величины ß (А)

и Р (А) также

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

ние А, при котором среднестатистическая погрешность наибольшая (рис. 3).

В работе [43, стр. 153] величина ДF определяется несколько иначе.

51

Поскольку максимальное значение функции F (б) получается

при алг. 6+ как сумма приращений

Аңо, то суммарную погрешность

ßs можно записать в таком виде:

 

 

 

АF (во).

р (л)Р(л) 2t А г+«)

 

(2.41)

ß s:

 

Ѵ(А)Р(А).

^(б+)

S a+(0

 

 

Учитывая (2.36) и (2.39), окончательно получим (рис. 3):

 

 

 

ßx =

а 2 (1 — а)

(2.42)

 

 

4 ( 1 ф а )

 

где

а = V 1— 1/А •

 

 

 

Проанализируем результаты.

 

 

Из (2.39)

следует, что

макси­

 

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

 

а - ѵ і

оо) и равна 25%, т. е. едва

 

раза меньше, чем у метода МЭ. Важно

 

отметить,

что максимальная погреш­

Рис. 3.

ность

метода

МЭ получается

при

А =

1 (а == 0)',

однако в этих

усло­

виях алг. б0 практически не дает сбоев [см. (2.36) при а

= 0], да и

сама погрешность мала [см.

(2.39) при а = 0].

 

Условие А^уоо (а -И ),

при котором возможна максимальная по­

грешность алг. б о, не совпадает с «опасной» областью для алг. б+. Сле­

довательно, алгоритмы б „

и б+ взаимно компенсируют

недостатки

друг друга и одновременный

счет по ним исключит с высокой гаран­

тией всякую погрешность.

 

 

В случае применения матриц больших размеров (10 X

10 и более)

в двойном счете нет необходимости, так как вступают в силу законы статистики, и если ММЭ имеет оценку погрешности сверху не выше 6,3%, то для МДФ она еще меньше — не более 3% (см. рис. 3).

Алг. б о

требует несколько большего

объема вычислений,

чем

алг. б+. Объем вычислений

определяется

примерными

оценками:

умножений

&SN (N + 1),

сложений

«

(S — 1) N (N +

1),

срав­

нений &SN (N + 1).

 

 

 

 

 

При разработке машинного алгоритма следует учитывать, что

максимальный элемент в каждой строке

матрицы ЦДыЦлга

соответст­

вует максимуму выражения

 

 

max Ij4,ö>w -

Ai ты П en \ *

(2.43)

1 < 7 < S

ehl

 

что следует из (2.17). Указанное обстоятельство, а также пересчет элементов матрицы от (t — 1)-го шага процесса к t-му позволяют создать экономные алгоритмы решения задачи.

* См. стр. 30.

52

4.

Некоторые

обобщения

и частные

случаи.

Полученный ре

зультат может быть расширен применительно к

следующим случаям.

а)

Р е ш е н и е

з а д а ч и б о л ь ш о й

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

Пусть

все N средств составляют N типов.

Все

N}

средств

однородны

и

потому внутри своей группы не нумеруются.

Тогда задача (2.1) —

(2.3) запишется в следующем виде.

 

\\х°і\\,

 

 

 

Требуется

найти

матрицу

 

Х 0 =

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

симум

функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

'

N

 

\

 

 

 

 

 

F(X)

 

 

1

 

 

 

(2.44)

 

 

 

= 2

 

4

1 — П

8/ / г ]

 

 

 

 

 

 

'V

/ = 1

1

1

 

 

 

 

 

 

 

i=i

 

 

 

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

на переменные

 

 

 

 

 

 

 

 

S

Xj

 

 

 

1....... , N,

 

~N

Nj =N*

(2.45)

 

S

 

 

/ =

 

2

 

і= 1

 

 

 

 

 

 

/=i

 

 

 

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

условиях:

 

 

 

 

 

 

 

 

 

 

Хп € {0,

1, ••• > N},

 

 

 

N,

 

 

 

 

° < ( ел = 1— ©л)<

1. Г.~ .

 

(2.46)

 

 

 

S.

 

 

*

 

 

 

 

 

і = \

 

 

 

Для решения задачи можно использовать алг. 1.2 со следующими изменениями:

элементы матрицы || Д 1 1| 2vs будут вычисляться по формуле

(2.47), вместо (2.21):

N

(2.47)

Ш

ек і /= 1 '

элемент xjt определяется как суммарное количество средств, назначенное от /-й группы (источника средств) для воздействия по і-й цели за все N шагов процесса распределения:

N

_

 

Х ц = Ъ Ь и і ,

N.

(2.48)

Так как N ^ N , то возможностьразделения средств по группам позволяет существенно увеличитьразмерность (по j) решаемой зада­ чи при том же объеме ОЗУ. Объем вычислений также сократится. Оценка по числу элементарных операций показывает, что для реше­ ния потребуется:

умножений «

SN(N-}- 1),

 

сложений «

(5 — 1) N (N + 1),

(2.49)

сравнений «

S N (N -f- 1),

 

* См. § 1.3(1).

53

При большом количестве средств N объем расчетов может быть еще более сокращен. Это вытекает из того, что в этом случае О и поэтому [см. (2.47)]

Ahl = Mti — Д« « Аti = А\*~ 0

(2.50)

Это значит, что вначале распределение можно производить соглас­

но алг. б+ до тех пор, пока не выполнится равенство

 

2 x n = Nr

(2.51)

І

 

для г-го (г = 1, ..., N) источника средств.

В этом случае средство г-го источника, которое привело к выпол­ нению равенства (2.51) возвращается в тот же источник и дальнейшее распределение ведется методом двух функций (алг. 2.1), причем сум­ марное приращение Aki вычисляется по формуле (2.21) в предполо­ жении, что в каждом источнике только единица средств.

Если после перехода к алг. 2.1 все же будет назначено средство

г-го источника, то можно снова перейти к алг. 6+ до выполнения условия (2.51) для какого-либо дугого источника.

Если после перехода к алг. 2.1 последнее средство г-го источника пока не назначается, то распределение ведется по алг. 2.1. Переход к этому алгоритму должен осуществляться с тем расчетом, чтобы ис­ ключить возможность распределения последнего из остающихся в каж­

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

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

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

б) Д о п о л н и т е л ь н ы е

о г р а н и ч е н и я

вида:

 

N

 

 

 

хі = ' 2 1хіі =1 ,

/ = 1, ..., S

(N —S).

(2.52)

;= 1

 

 

 

При таких ограничениях задача сводится к задаче о назначениях

[ 10].

Алгоритм решения изменится лишь в том, что после каждого шага процесса уменьшается на единицу не только количество средств, но исключается из дальнейших расчетов и цель, на которую назна­

чается средство.

 

 

в) Д о п о л н и т е л ь н ы е

о г р а н и ч е н и я вида:

 

Хі = Ъхи £ { 0 , 1 } ,

f = l ....... S ( N < S ) .

(2.53)

/

 

 

Задача может быть названа задачей о назначениях с отбором целей

(незамкнутая задача о назначениях).

54

По сравнению с предыдущим случаем область оптимизации рас­ ширена, зависимость между назначениями ослаблена и вероятность

получить погрешность меньше.

 

г) Д о п о л н и т е л ь н ы е

о г р а н и ч е н и я вида:

 

Xt = HxJt= l ,

/ = 1 ....... S (N > S ).

(2.54)

/

 

 

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

Формальным условием для замены &г го средства г-м (не вошед­ шим в решение) является неравенство:

<ükt it <®rit , /t= l , . . . , S .

(2.55)

Номер г наиболее эффективного из множества остающихся не­ использованными N' средств определяется из условия

и>гі

= max töji

(2.56)

*

/елг

1

 

д) Д о п о л н и т е л ь н ы е

о г р а н и ч е н и я

вида:

*і = 2 * я € { 0, 1,...,

üj},

г = 1,...,

S (и, > 0). (2.57)

і

 

 

 

Это более общий случай, включающий все предыдущие (кроме а).

Если ^ v ^ N , то в алг. 2.1 необходимо проверять дополнитель-

і

ные условия (2.57) и по мере их выполнения исключать из дальней­ шего рассмотрения соответствующие им цели.

Если 2 к г < N, то задача решается с «отбором средств» и требует

І

последующей проверки условий (2.55), как и в случае г).

е) Алг. 2.1 позволяет решать обратную задачу, которая формули­ руется следующим образом. Выбрать минимально необходимое число

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

min N = 2 Xi = 2 2

Хи

(2.58)

X

I

i i

 

 

при ограничениях [см. (2.44)]:

 

 

 

F ^.F3,

F i^ F ig

для

некоторых i.

(2.59)

Следует иметь в виду, что область возможных решений существует только при условии

Ъ ? 1г > Рд .

(2.60)

і

 

Задача будет окончательно определена после указания совокуп­ ности средств, из которой производится выборка. Если указаны только

55

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

в § 1.3 (1). Она решается на основе.алг. б+.

При ограничении числа средств каждого типа следует исполь­ зовать алг. 2.1 или его обобщение а). После достижения требуемого эффекта (2.59) полезно убедиться, что оставшиеся средства не могут улучшить решение путем замены согласно условию (2.55).

Рассмотренные выше видоизменения задачи не влекут за собой снижения точности решения, так как по отношению к ним практически

целиком справедливы доводы,

приведенные в § 1.1 (4).

5. Примеры расчета. Задача

(2.1)—(2.3) решается для исходных значе­

ний величин, заданных вектором

{4j}s

и матрицей II й/г || NS:

{10,0

8,00

6,00

2,00}

0,40

0,40

0,20

0,50

1,00

0,20

0,80

1,00

0,80

0,10

0,60

1,00

1,00

0,60

0,00

(2.61)

1,00

0,50

0,00

0,40

0,70

0,20

0,30

0,40

0,50

Вычислительная схема, соответствующая алг. 2.1, представлена табл. 11, элементы которой записаны в виде пар чисел (число без скобок и число в скоб­ ках). Каждая пара чисел вычисляется по формулам:

At (i>hl

N

I A i(ü k i

П ец;

I “ TT11

Ъы

 

/ = 1

e&z

первое

число

второе число

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

[см. (2.21)], которое записывается в предпоследнем столбце таблицы. Максималь­ ный элемент этого столбца определяет второй индекс — индекс средства (&<), распределяемого на t-м шаге процесса. Максимальные элементы столбцов также выделены жирным шрифтом.

Характер изменений функций F f, Fp и F f уже был показан на рис. 1. Матрицами б0 и 6+ представлены решения, полученные с помощью алгоритмов б0 и 6+:

0

1 0

0

 

 

0

1 0

0

1 0

0

0

 

 

0

0

1 0

0

0

1 0

,

Ö+ -

0

0

0

1

0

1 0

0

1 0

0

0

 

 

0

0

0

1

 

 

0

0

1 0

0

0

1 0

 

 

0

1 0

0

алг. б0, F0 =

22,04;

алг. б + ,

F $ = 21,92.

56

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

 

11

 

 

 

 

 

 

 

І,

1

 

 

•ч

 

 

 

 

 

1

 

2

1

3

1

 

 

 

 

 

 

 

4

<

 

 

t

k, і

 

 

 

 

л (° )

 

 

 

II

 

 

 

 

 

 

 

 

 

<

 

 

 

 

 

 

 

 

А1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

10,00

 

8,00

1

6,00

1

2,00

Л _

*вГ

 

 

 

В

 

1

4 , 0 0

(0.00)

 

3,84 (0.64)

 

1,24(0.04)

 

1,00(0.00)

3,32

 

 

 

2

10,00

(0.00)

 

1,84 (0.24)

 

5,35 (0.55)

 

2,00(0.00)

9,21

 

 

 

3

8 , 0 0

(0.00)

 

0,91 (0.11)

 

3,81 (0.21)

 

2,00(0.00)

7,68

2;

1

 

4

10,00

(0.00)

 

6,25(1.45)

 

0,00(0.00)

 

2,00(0.00)

8,55

1

5

5 , 0 0

(0.00)

 

0,00 (0.00)

 

2,49 (0.09)

 

1,40 (0.00)

4,91

 

 

 

6

2,00

(0.00)

 

2,81

(0.41)

 

2,49 (0.09)

 

1,00(0.00)

2,31

 

 

2

О

4

5

6

 

0,00

8,00

6,00

2,00

 

1

_

4 , 0 0

(0.80)

1,37(0.17)

1,00 (0.00)

3,03

3

0,93 (0.13)

6 , 3 6 (2.76)

2,00 (0.00)

3,47

4

6,61

(1.81)

0,00 (0.00)

2,00 (0.00)

4 ,8 0

5

0,00 (0.00)

2,86 (0.46)

1,40 (0.00)

2,40

6

2,91

(0.51)

2,86 (0.46)

1,00 (0.00)

1,8У

 

 

 

Л<2>

 

 

 

0,00

3,20

6,00

2,00

 

1

2,08(0.80)

1,37(0.17)

1,00(0.00)

1,11

3

0,45 (0.13)

6,36(2.76)

2,15(0.15)

3,3 2

5

0,00(0.00)

2,86 (0.46)

1,40(0.00)

2,40

6

1,47(0.51)

2,86 (0.46)

1,00 (0.00)

1,89

 

 

 

л<3>

 

 

 

 

 

 

1

 

 

 

0,00

3,20

2,40

2,00

 

1

_

3,17(0.89)

0,78 (0.17)

1,15(0.15)

0,96

5

0,00 (0.00)

1,42 (0.46)

1,75 (0.35)

0,94

6

1,53 (0.57)

0,42 (0.46)

1,15(0.15)

0,35

 

 

 

л<4>

 

 

 

0,00

 

1,92

2,40

2,00

 

5

0,00 (0.00)

1,54 (0.58)

2 ,1 0 (0.70)

0,8 2

6

1,13(0.57)

1,54 (0.58)

1,30 (0.30)

0,09

 

 

 

іл<5>

 

 

 

0,00

 

1,92

2,40

0,60

 

6

0,58

0,96

0,30

0,08

 

 

 

л<6>

 

 

 

0,00

 

1,92

1,44

0,60

 

4; 2

1; 2

£ 4.

6; 3

57

Текущие матрицы назначений, соответствующие алг. б0, представлены

ниже.

 

для P t

 

для Ft

 

 

для P t

 

для Ft

 

 

0

0

0

0

1

1

1

1

 

 

0

0

0

0

1 1 1 1

 

0

0

0

0

1

1

1

1

 

со

1 0

0

0

1 0 0 0

 

0 0

0 0

 

1 1 1 1

0 0 1 0

0 0 1 0

 

 

 

 

 

 

 

 

 

 

II

 

 

 

 

 

 

 

0

0

0

0

 

1

1

1

1

633— 1

0

1 0

0

0 1 0

0

 

0

0

0

0

1

1

1

1

 

 

0

0

0

0

1 1 1 1

 

0

0

0

0

1

1

1

1

 

 

0

0

0

0

1 1 1 1

 

0

0

0

0

1

1

1

1

 

 

0

1 0

0

0 1 0 0

 

1 0

0

0

1 0

0

 

0

 

 

1 0

0

0 1 0 0 0

t = l

0

0 0 0

1

1

1 1

 

II

0 0

 

0

0 0 1 0

б21 = 1 0 0 0 0

 

1 1 1 1 бі 2 — 1 0 1 0 0

0 1 0 0

 

0

0

0

0

1

1

1

1

 

 

0

0

0

0

1 1 1 1

 

0

0

0

0

1

1

1

1

 

 

0 0 0 0

1 1 1 1

 

0

0

0

0

1

1

1

1

 

 

0

1 0

0

0 1 0 0

 

1 0

0

0

 

1 0

0

0

1 0

0

0

1 0 0 0

t = 2 0 0 0 0

 

1 1 1 1

t = 5

0 0 1 0

0 0 1 0

б42= 1

0

1

0

0

 

0

1 0

 

06 5 4 = 1

0

1 0

0

0 1 0 0

 

0

0

0

0

1

1

1

1

 

 

0

0

0

1

0 0 0

1

 

0

0

0

0

 

1

1

1

1

 

0

0

0

0

1 1 1 1

 

 

для /=)+

 

 

0

1 0

0

 

 

10

0

0

t =

6

0 0

1 0

Ö«3 =

1

0 1 0

0

 

 

0 0

0

1

 

 

0 0

1 0

для

0 1

1

0

0 0

0 1

0 0

0 0

Ft

0 0

0 0

1 0

0 0

0 1

1 0

Пример хорошо иллюстрирует «чувствительность» алг. б0, которым «замечена» разница решений в 0,55%:

Ро~Р£

22,04—21,92

 

 

 

Ft

 

= 0,0055

 

Решение б0 оптимально.

 

22,04

 

 

получено

решение и при введении дополнительных

С помощью табл. 11

ограничений (случай «г»

обобщений)

вида:

 

 

2 6 7., = л:г= 1 ,

t = 1____ s,

N > S .

(2.62)

і

 

 

 

 

 

Решение также оптимально. Последовательность назначений и отбора средств при этом следующая:

(kt, /() = (2; 1 4-, 2-*3; 3 ^ 5 ; 4).

Средства / = 1 и / = 6 при этом оказались забракованными. Улучшения решения путем замены средств достигнуть не удалось.

58

§ 2.2. Обобщенная задача оптимального распределения неоднородного ресурса между взаимозависимыми элементами системы

 

1. Постановка обобщенной задачи.

Формальная постановка зада­

чи

следующая.

 

 

 

 

 

 

 

 

 

 

 

Требуется определить

матрицу 8 0 =

|| б“/ || n s , доставляющую

максимум функции

 

 

 

 

 

 

 

 

 

 

 

s

 

г

s

 

 

 

N

 

 

 

F(8)=

2

А;

1 - П

1а}і I 1—

П

8.ву/

(2.63)

 

 

г=і

 

{

і= 1

 

 

 

 

VI

 

 

 

 

 

 

 

 

 

 

при

ограничениях

на

переменные

 

 

 

 

 

 

 

 

 

2

б ѵ, = 1,

Ѵ =

1,

,

N,

 

(2.64)

 

 

1=1

 

 

 

 

 

 

 

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

условиях:

 

 

 

 

 

 

 

6VJ-6 ( 1 ,0 },

 

 

 

 

 

 

 

 

0 ^ (еѵ/ =

1 — (öVJ) sc: 1,

i, / — l,***,

 

 

0 < ( ан =

1 — O jiX

 

(2.65)

 

 

V =

1........

N.

 

 

 

 

 

 

 

 

®U= !•

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

Имеется N активных средств (ракет).

Вероятности поражения

ѵ-й ракетой (ѵ = 1, ..., N) j-го объекта (/ =

1, ..., S) заданы матрицей

®v/!wsОтносительные значимости (веса)

объектов заданы векто­

ром { Л;}5. Если в результате обстрела /-й объект (цель) уничтожен, то

с вероятностью

выходит из строя (прекращает функционирование)

і= 1, ...,

S)

объект.

 

 

Необходимо так распределить ракеты по объектам (найти матрицу

6 о = I Iбѵ/І),

чтобы

максимизировать

функцию ущерба (2.63),

учитывающую

связи

между

объектами

(коэффициенты связи аji).

Отличие от постановки задачи

§ 2.1 (1) обусловливается учетом зави­

симостей, существующих между объектами. Это требует рассмотрения

и объектов, имеющих собственно нулевые веса (Аи = 0),

если с дру­

гими

объектами

они

связаны

ненулевыми

коэффициентами

связи

(ahi ^

0)-

 

 

 

 

 

 

Если /-й объект уничтожен, то он прекращает функционирование

и, значит, ап — 1

для всех Если і-й объект не связан с /-м,

то

0.

Следовательно, в случае независимых объектов имеем:

 

 

 

аН= 0,

<%}}= 1,

/, г'= 1,...,

S, ! ф і .

 

(2.66)

59

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