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

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

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

Нетрудно убедиться, что подстановка условий независимости (2.66) в формулу (2.63) позволяет привести ее к виду (2.1).

Величину ал, обратную по смыслу величине ан , можно понимать как коэффициент снижения важности Л; і-го элемента в случае раз­ рушения /-го.

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

Пусть имеется система, состоящая из S элементов, между кото­ рыми установлены функциональные связи. Степень связей количест­

венно

характеризуется матрицей

flayjflss.

 

Количественным

показателем

эффективности функционирования

системы является критериальная функция F (6) или эквивалентная ей

функция F (6):

s

 

 

 

 

 

 

 

 

Р(6)= Ц Аі—F{b).

(2.67)

 

 

I=i

 

 

Задача состоит в таком распределении N единиц неоднородного

ресурса между элементами системы (матрица 6

= |6 v/||ws). которое

обеспечит минимум критерия F(б) [максимум

F( б)].

Рассмотрим подробнее функцию F (6). Из (2.63) видно, что функ­

ция

F (6) аддитивна

(суммарный

эффект складывается из эффектов

на элементах системы). Внешнее влияние (воздействие) на /-й элемент системы при некотором произвольно заданном распределении ресур­ сов (матрица 6) будем вычислять по формуле

°J>i=

Л

в -

(2.68)

1 — П

еѵ]>, / = 1,..., 5 ,

 

Ѵ =

1

 

полагая, что отдельные средства влияют на /-й элемент независимо

друг от друга. Подставив (2.68)

в формулу (2.63), получим:

 

s

/

Si

 

F(8)= 2

Лг

1 - П

(2.69)

 

V

і = 1

 

Из (2.69) видно, что полная вероятность а* выхода из строя (пре­ кращения функционирования) і - г о объекта как за счет внешнего воз­

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

сс,= 1— П ( 1 — SFjaH).

(2.70)

Полная вероятность невыхода из строя t-го объекта:

 

« t = 1 — a i = П ( 1 — SPjUji)-

(2.71)

i

 

«о

Тогда текущая важность (вес) і-то объекта

после распределения

(t — 1) единиц средств

запишется так:

 

 

А ? - 1) =

А ; • otf = A t п (1

а л)-

(2.72)

 

I

 

 

Следует обратить внимание на то, что нарушение функциональ­ ных связей в рассматриваемой схеме учитывается только между элементом, подвергшимся внешнему воздействию SPj, и остальными элементами. Здесь не учитываются более глубокие, ветвящиеся связи, которые могли бы характеризоваться цепями Маркова [44].

Из всей цепи

3bj - + a j i - + a il-+ai T-H

т. д.

(2.73)

в задаче учитываются только первые связи

(связи первого

порядка):

І=£ І ’ С / = 1 ,

 

(2.74)

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

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

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

2. Метод оптимизации и алгоритм решения. Для оптимизаци критерия F (б) примёним метод двух функций. Получим выражения для функций Ff, Ff и их приращений для ^-го шага процесса, если k-e средство назначается по 1-й цели*.

Пусть к моменту ^-го шага уже распределено t — 1 единиц ре­

сурса так, что эффективность воздействия

по /-й цели

равна

 

/ =

1,... , S .

(2.75)

Согласно матрицам (2.6) для функции выигрыша с учетом (2.63)

имеем:

 

 

 

F t i = І

At ( l — П(1 —

о.н)),

(2.76)

І = 1

J

 

 

2 Л * 0 - П ( 1 - ^ % г))-

(2.77)

Используется интерпретация задачи, данная на стр. 59.

61

Неотрицательное приращение функции выигрыша будет равно

Аt i = F t - F t v

2 Аі

П ( і

_gb(t- I )

« / / ) — Щ і -

-S6}0 « * )

.(2.78)

f= 1

/= і

 

 

 

 

 

С учетом (2.72)

формула (2.78) приводится к виду:

 

 

 

 

 

 

П(1

(Xj. )

 

 

 

 

 

 

П(1—

a}.)

.(2.79)

Назначение на t-м шаге процесса k-ro средства по /-й цели ведет

к изменению 1-й компоненты вектора

 

поэтому

 

 

 

(t-i)

если

j =f= I,

 

(2.80)

 

 

 

 

 

 

 

1 q V l) ehl, если

/'=/ .

 

(2.81)

После сокращения формула (2.79) приводится к окончательному

виду:

 

 

 

 

 

 

 

АЙ= І

---- ^

т Г<<) = ^ ~

1> VQ<C -I r n

• (2-82)

Сравнивая средние части равенств (2.82) и (2.79), можно получить рекуррентную формулу для вычисления значений А\1\ которые не­ обходимы для формирования алгоритма:

а Р =

а \‘- 1)

1—

ац

 

.83)

 

І =

1, ... , S.

Для учета потерь, которые могут возникнуть при использова­

нии алгоритма б+,

получим выражение для

неположительного при­

ращения — Аы функции

потерь.

 

 

 

По аналогии с

§ 2.1

суммарное воздействие по /-му

элементу до

распределения /-го средства будет включать действие уже распределен­

ных (/ — 1) единиц ресурса

и влияние

оставшихся Nt средств

Ісм. (2.19)]. С учетом матриц (2.6) можно записать:

 

/Т - > = І М г( і - П ( 1 - ( 1

 

(2.84)

где

 

 

 

 

а/і - І ) =

П

Ëvj, } = 1...... N;

(2.85)

— множество еще не

распределенных

к моменту

/-го шага

средств.

 

 

 

 

62

Для функции F t в соответствии с (2.6) и на основе (2.63) по­ лучим*:

Ft = ,?1Л ,| 1—

. /.

Х

X п 1

 

8Ь£

аи

 

 

 

Q(.'-')a(<-1:

 

1-

1

( 2. 86)

^kj

/= 1

 

 

 

 

 

Ha основе формул (2.84) и (2.86) можно получить выражение для неположительного приращения функции потерь при 6hi= 1:

 

- а Г і = 2 М ,

п [1—(1 - ( &

*

- "

 

 

і

 

 

 

 

 

 

l - ( l

- Q ^ - ' ’а ^ - 1)) ап

п

1

и

 

I) . (2.87)

 

 

 

4 j

 

 

 

 

I 1

I

V

8ki

1

 

 

 

 

 

Поскольку положительный первый член в формуле (2.87) не зависит от индексов k и /, то он может быть опущен. Если учесть (2.72) и (2.82), то после некоторых преобразований получим оконча­

тельное выражение для максимизируемого элемента Л*/* ~ Д ы —

= А*/ — А

 

 

Ым^

1>(Ѵ

 

i - ( i -

q</-» «</-■>)«„

X

 

 

 

 

 

 

 

1—â!>\t~ 1)ali

i - і 1

а . .

 

 

 

 

 

ll

 

п 1 -

 

 

 

 

 

eģ

 

1

Bkj

ан

,

k £ N W , j, l = l , . . . , S .

(2.88)

X

 

 

 

Теперь можно сформулировать алгоритм решения задачи (2.63)

(2.65).

 

 

 

 

 

 

 

Алгоритм

2.2.

 

 

 

 

величин по формулам:

1°. Вычислить начальные значения текущих

 

 

 

(Q(°> = l),

 

 

 

 

 

 

N

 

 

і, j = 2 ........S.

(2.89)

 

а(0)= П

еѵ,

 

 

/

А *

ѴМ

 

 

 

 

 

1

Ѵ = 1

 

 

 

 

 

А{0) = Аі,

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

63

2°. Вычислить элементы текущей матрицы

|| Л*/ ||jvs 110 Ф°РмУле (2.88)

для k £

N ^ \ 1 = 1 ,

..., S.

 

определяемых

из условия:

3°.

Выбрать и запомнить пару индексов kt, It,

 

Akt \ t = max &ki ,

/ = 1 , . . . , S .

(2.90)

4°.

Пересчитать

текущие значения величин по формулам:

 

 

 

 

если

ІФ lt>

(2.91)

 

 

[ l —Qf

0 s k j, если j = l t ,

 

 

 

 

 

,V>= -L

 

S ,s

(2.91')

 

 

4

/

 

 

 

_ л « - і > _____ 1__ ILL ,

S.

(2.9Г)

чч 1

5°. Вычислить текущее значение функции F^~;

 

 

 

s

 

(ok ,

а,

,

p +

p + А Д +

c +

'V

A^~ П

—-

-_______-

(2.92)

 

 

F + = 0 ,

<: =

/ +

!.

 

 

6°.

Проверить условие

t < N:

 

 

 

 

 

 

да — перейти

к п.

2°;

 

 

 

 

 

 

нет — перейти

к п.

7°.

 

 

 

 

 

7°.

Отпечатать результат

(F (60) =

F ^ ;

{kf,

lt },

/ = 1 , . . . ,

N; или ||б°(.|),

прекратить вычисления.

3. Оценка эффективности алгоритма. Для оценки среднеста­ тистической погрешности рассмотрим вначале два предельных по сте­ пени зависимости между элементами системы случая:

1)

ал =

0

(а]}= \) , j,i= 1.......

S,

(2.93)

2)

<Х]і=

1,

j,i= l, ... ,S .

 

(2.94)

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

1) Н е з а в и с и м ы е э л е м е н т ы . Совпадение решений, по­ лучаемых с помощью алг. 2.1 и 2.2, при условиях (2.93) гарантирует-

* См. сноску на стр. 48.

64

ся, если совпадут соответствующие формулы для вычисления величин

АР( , Аkt и Аki- Такое совпадение имеет место. Действительно, при условии (2.93) из (2.91") имеем:

 

 

 

Ä( t ) [ A i{

 

 

 

если і=£1,

 

(2.95)

 

 

 

 

I A t ~ X) гкі,

если /==/.

 

 

 

Сравнивая

(2.95) и

(2.24),

видим,

что

необходимое

совпадение

в этом случае достигается.

 

 

 

 

 

 

 

 

 

Аналогично

из

формулы (2.82)

при условиях (2.93)

получаем:

 

 

 

 

 

Л Й = Л (/ - ' Ч

г,

 

 

(2-96)

что

и требовалось [см.

(2.7)].

 

 

 

 

 

 

 

Из формулы (2.87) при условиях (2.93) получаем:

 

 

 

 

 

 

 

2 л г Q-Г

&ki

 

 

 

(2.97)

 

 

 

 

 

i¥=l

 

 

 

 

 

 

Учитывая,

что

А fit*

 

!)

= А\‘

и объединяя два

первых

члена

в один, получим

окончательно:

 

 

 

 

 

 

- 4 5 = 2 4 : ‘ - 1І4 ' - І,- 2 - 4 Г п

а( t - 1)

 

 

 

 

&hi

 

 

 

 

 

іф I

 

2

 

 

і ФІ

 

 

 

 

 

 

 

-

-

Л ' - 1».a ki

П

6vi

 

(2.98)

 

 

 

 

ІФ

 

eki

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что

и

требовалось

[см.

(2.13)].

 

 

 

 

 

 

2)

З а в и с и м ы е

э л е м е н т ы .

При условиях (2.94)

целевая

функция (2.63)

принимает вид:

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

s

N

e6Y/

 

 

 

 

 

F(8)

2

А, I

П

П

 

(2.99)

 

 

 

VI

 

 

 

 

 

 

i=

1

 

/ =

1 V =

1

 

 

 

Из (2.99) следует, что

 

 

 

 

 

 

 

 

 

max F (б) =

г

— (2^ г ) min (П П e®p)~max2öv./®v.7-*

(2.100)

 

 

6

 

і

 

 

б

/ ѵ

 

б /,ѵ ■

 

 

Таким образом, из (2.100) следует, что оптимальный алгоритм состоит в выборе N наибольших элементов матрицы [[ соѵ/ 1] л/’s — по одному из каждой строки, т. е. на каждом шаге процесса назначение следует производить согласно условию

^ = maxtöw,(6ftt lt= 1), k^NW, Z = 1 ,...,S . (2.101) k,I

* Произведению N наименьших элементов Evj (по одному из строки мат­ рицы II Еѵі II) соответствует минимум их суммы (т. е. максимум суммы элемен­ тов соѵу=1—£vj)

3 Зак. 1292

65

Покажем, что алг. 2.2 обеспечивает такое же решение, как и опти­ мальный алгоритм, основанный на требовании (2.100).

Из формулы (2.88) при условии (2.94) получаем:

 

І

\

і

®Äehj7

/

 

 

 

 

(0hl'

П

(2.102)

Из (2.102) следует, что

 

^ k j i

 

&kj

 

 

 

 

 

 

шах

шах a>hl—cow ,

1, ..., TV7.

(2.103)

1 </ < s

/

 

г

 

 

 

Таким образом,

каждый

шаг

соответствует

назначению

одного

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

Выясним теперь, как зависит погрешность алг. б+, возникающая при решении задачи § 2.2 (1) от значений исходных величин (в том числе и коэффициентов связи а^). Это необходимо для того, чтобы показать монотонный характер изменения погрешности между двумя рассмотренными предельными случаями. Тогда оценка погреш­ ности определится одним из этих случаев. По аналогии с § 1.3 (3) рассмотрим задачу 2 x 2 , заданную матрицами

 

I I

®

f

® i

®

> .

И

v

1

2

1

а

2 =

1

2

 

f

 

2

 

j1 а

f

1 l

 

1(2.104)

и вектором

[Аъ А 2).

 

 

 

 

 

 

 

 

 

 

 

 

 

Будем полагать, что оптимальным решением является решение

612 = 1 и

б21 =

1. Элемент

 

 

со22

матрицы || соѵ/ [|

приравняем

 

нулю.

При оптимальном процессе решения значение целевой функции

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

исходных

параметров будет равно

 

 

 

 

/%>= 2

i A Ä (0 = A ti(l) +

At2(2)*,

 

 

(2.105)

где для данных условий задачи в соответствии с формулой (2.82) имеем:

 

 

 

Д21 (1) = ^4і ®21 Н"~ А 2®21^12>

 

 

 

Aj2 (2) — А 2 (1

®21®1г) ®12 “Ь ^4 1 (1

®2і) ®12®^21-

(2.106)

Пусть

 

алг. 6+ обеспечивает решение

8И = 1,

б21 =

1.

Для

этого необходимо и достаточно, чтобы соблюдались условия:

 

 

для t= 1:

 

Ап ( 1 ) > АД (1) и Ап ( 1 ) > АД(1)**

 

(2.107)

* Запись

Д+; (t)

заменяет выражение Д £ ^

можно

найти

и по

фор­

муле (2.63)

для 6і 2 =

1, б * і=

1.

 

 

 

 

** Дц (1)

>

Д22 (1),

так как

со22 = 0.

 

 

 

 

66

 

Тогда

для

t = 2 имеем ö21 =

1 и решение,

полученное на основе

Ö+,

в зависимости от исходных

параметров запишется так:

 

 

 

П = ДГі (1) + Д2+і (2),

(2.108)

где

Дц

(1) — Аі(Иц -f- А 2 ®ц К^2і

(2.109)

 

Даі (2) =

А 1 (1 — ®іі) ®2! + А 2 (1 — ЮцЯіа) «2іа і2-

Теперь задача определения максимальной относительной по­ грешности ß формально будет состоять в максимизации функции

ß _ J F t _ I

A l С0П

Ю11 и 12~|^^1 е 11 03 21

(1 — ЮХ1 а і г ) Ю21 а 12

Al (Ö21 "ф- А2 (Ö21 С6і2"Ф“Л2 О — ^21 аіг) ^12Ф Лі 821 ^12 &21

(2.110)

при условиях (2.107), которые соответственно примут вид:

А і(Оц + Л 2соиа 12 > А хю21 + А 2(о21а 12,

Лі<Йі + А2І0цКі2 !> Л2 (Оі 2 Ах(і)і2® 21-

(2.111)

Следует отметить, что хотя неравенства (2.111) строгие (это гаранти­ рует однозначность решения), для оценки максимальной погрешности, переходя к пределу, можно заменить их равенствами (см. § 1.3 (3))*, ко­ торые после элементарных преобразований соответственно примут вид;

Ю2і = 0ц . « і2 = «и

, г д е Л = ^ - .

(2.112)

1 -ф. Л а 2і

Л 2

 

Производя подстановку, приводим формулу (2.110) к окончатель­ ному виду:

 

 

Лвц-ф- (1 —t°ii иіг) иі2

 

 

ß = l

Л rj-aі2

 

 

(2.113)

 

Л811 G&21 "ф 1

0)ц а12

 

1Л-

 

 

 

1 -Ф" Л а 2і

 

 

 

 

 

 

 

В случае полностью зависимых элементов (а12

a 21 1)

имеем:

 

 

 

 

 

 

ß '= l -

8ii -ф О

 

1-f е11 (Л-^-l)

\ __ Q

(2.114)

 

Л-f 1

 

Л +

1

)

 

Впрочем этот результат не нов. Мы его предвидели при анализе случая (2.94). Именно поэтому ß = 0 является необходимым усло­ вием, которому должно удовлетворять выражение (2.113), если в ло­ гике его вывода и в преобразованиях не было допущено ошибки. Вы­

* В § 1.3 (3) стр. 30 мы вводили малые величины а > 0 и переходили к ра­ венствам.

3*

67

деляя в знаменателе второго члена выражения (2.113) целую относи­

тельно

а 21 часть, получим:

 

 

Лец -Ml — Mn а12) и12

 

Ä^-a12 ______

 

ß = l

(2.115)

 

1 -рeu ^-

о)ц (1 —а 12)

 

 

1 "4" а21 А

Из

(2.115) видно, что по мере

увеличения коэффициента связи

а 21 погрешность монотонно убывает (или остается постоянной при

а 12 =

1) независимо от значений остальных параметров.

Если к числителю второго члена выражения

(2.113) добавить и

вычесть единицу, а А + а 12

заменить, используя

для этого условие

(2.112),

то формула погрешности примет вид:

 

 

 

2

Mil

 

 

 

Wl2 1 ■ф’^4 а2і

 

 

ß = l

 

(2.116)

 

j .

А вц а2і^ 1 —Мц а12

 

 

 

 

 

 

1 -f- А сс2і

 

Из

(2.116) видно, что и

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

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

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

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

В силу аналогии структур алг. 6+ и 6 0 полученный вывод рас­ пространяется и на алг. 2.2.

Монотонный характер изменения погрешностей в зависимости от коэффициентов связи ап позволяет для оценок точности алг. 2.2 ограничиться худшим из двух рассмотренных предельных случаев

(2.93) и (2.94).

Согласно случаю (2.93) максимально возможная и среднестати­ стическая погрешности оцениваются сверху величинами 25 и 3%.

Для реализации алг. 2.2 требуется следующий оценочный объем элементарных операций:

умножений

S3N3,

сложений

t t S 3N3,

сравнений

&SN2.

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

Небольшого снижения объема вычислений можно достичь за счет некоторого ухудшения точности решения путем перехода к алг.

68

б+, пренебрегая величиной Дкі и полагая, что

АІі [см. (2.82).

(2.88)].

 

Использование обоих алгоритмов (б0 и 6+) дает надежную гаран­ тию в получении точного решения при любых сочетаниях исходных параметров.

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

Реальные системы часто могут иметь один основной элемент,

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

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

Тогда исходная задача выродится в задачу минимизации функции

F (б) = 1 — F (б):

 

N

 

 

 

F(ö) = П

П г6Ѵ

V j S -

(2.117)

/= і

ѵ =

VI

 

 

1

 

 

при тех же исходных ограничениях и условиях (2.64) и (2.65). Матрицу связей lo /ils s при этом можно свести к вектору-столбцу [| ajs [| s.

Алгоритм решения несколько упрощается и уменьшается объем вычислений (за счет сокращения числа операций при вычислении суммарного приращения). На основании (2.88) будем иметь:

l_ (l_ Q j* -l> aj*-l))

a i S

X

щ

 

 

 

\

Zhl

/

a is

 

 

 

(2.118)*

Объем вычислений можно несколько сократить еще и за счет

того, что произведение П (в числителе), как видно из (2.118), не зависит

/

* В § 2.4 будет дана более удобная формула для

учитывающая муль­

типликативность функции (2.117).

 

69

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