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

книги из ГПНТБ / Оптимизация процессов грузовой работы

..pdf
Скачиваний:
10
Добавлен:
27.10.2023
Размер:
26.66 Mб
Скачать

ГЛАВА 4

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

В качестве иллюстрации метода рассмотрим приведенный в [10] пример системы, состоящ ей из четырех последовательно соединенных подсистем. Рассмотрение начнем с системы, для которой требуемое значение коэффициента готовности равно 0,99, а фактическое его значение 0,31 . В табл. 4 .17 приведены данные, необходимые для реше­ ния задачи. В последнем столбце указано число устройств внутри каждой подсистемы, которое позволяет получить требуемую вели­ чину коэффициента готовности подсистемы, равную 0,99 . Однако даж е если все четыре подсистемы будут иметь коэффициент готов­ ности, равный 0,99, то и тогда коэффициент готовности системы будет равен (0,99)4 = 0,96 . Поэтому значения, указанные в последнем столбце, являю тся лишь исходными для расчета требуемого значения коэффициента готовности системы, равного 0,99 .

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

т1 - l g ( l - f t r )

I g ( l — A rf) ’

где kvi — коэффициент готовности t'-й подсистемы.

И сходны е данны е для

расчета

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

алго р и тм а ,

Т а б л и ц а 4.17

 

позволяю щ его

с о кр а ти ть число

членов

опти м а льно й последовательности

 

Номер

Стоимость

Коэффициент

Коэффициент

Число

подсистемы

готовности

простоя

резервных

 

 

 

 

 

 

 

элементов

4—

2

-

0,8

0,2

3

 

2

 

3

 

0,7

0,3

4

 

3

 

1,5

 

0,8

0,2

3

 

4

 

1

 

0,7

0,3

4

.1 2 5

МЕТОДЫ ПОВЫШЕНИЯ готовности

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

kt-с— I — (1 - Q m,

где k'rj — коэффициент готовности /-го элемента і-й подсистемы. Алгоритм решения заклю чается в следующем. Сначала подсистемы

группируются в пары следующим образом:

Д ля первой пары подсистем составляется матрица, в верхней строке которой проставляю тся числовые значения стоимости и коэффициента простоя, получающиеся в результате увеличения кратности резер­ вирования в первой подсистеме. Аналогичным образом вычисляется левый боковой столбец матрицы, характеризующ ий вторую подси­ стему. Стоимость сц и коэффициент простоя кпц, полученные в ре­ зультате объединения членов і-й строки бокового столбца и /-го столбца верхней строки, соответственно равны

 

 

с і і = с и

+

с\/!

 

 

 

 

I?

. — Ь

_1_ ь

 

 

 

 

Кпіі

і I

 

 

 

 

где с 2і и kn2i

характеристики

второй подсистемы;

с 1;-

и knll-

характеристики

первой подсистемы.

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 4.18

Получение членов промежуточной

оптимальной последовательности

 

 

*112

 

0 , 0 0 8

 

0 , 0 0 1 6

0 , 0 0 0 3 2

 

0 , 0 0 0 0 6 4

 

6

 

 

8

1 0

 

12

с2

 

 

 

 

 

 

 

 

0,0081

0,0161

 

0,0097

0,0084

 

0,0082

12

 

18

 

 

20

22

 

24

 

 

 

 

 

 

1

2

 

0,00243

0,0104

 

0,004

0,0027

 

0,0024

15

 

21

 

 

23

25

4

27

 

 

 

 

 

 

3

 

0,0073

0,0087

 

0,0023

0,0011

 

0,00079

18

 

24

 

 

26

28

 

30

 

 

 

 

 

 

5

6

7

126

 

 

 

 

Г Л А В А 4

 

Так как заданный коэффициент готовности системы равен kr =

=

0,99, то ограничение по коэффициенту

простоя

равно kn огр =

=

0,0 1 . Ограничение по стоимости определяется для каждой матрицы

по

формуле

П

 

 

 

 

 

 

 

Согр k =~ ^orp п

Сlint, h

2, 3, . . . ,

tl,

 

 

C=k+ 1

 

 

 

Матрица, вычисленная

в соответствии

с этими правилами, пред­

ставлена табл. 4 .18 . Начинаем вычислять оптимальную последова­

тельность. Первый член ее равен { 0,2 о9 7 }. В соответствии с алго­

ритмом проверки записываем этот член в боковой столбец второй матрицы (табл. 4 .19) и последовательно проверяем все комбинации этого члена с членами последовательности, относящимися к третьей подсистеме, записанными в верхней строке. Характеристики комби­

наций вычисляю тся таким ж е образом, как

и у членов оптимальной

последовательности:

 

 

 

 

 

 

 

 

 

СЧ=

сз/ “Ь с,,

 

 

 

 

 

 

 

К ч= Кг, + Кі-

 

 

 

 

Здесь / =

1, 2,

. . . —

номера

членов оптимальной последователь­

ности; с,

и knl — характеристики проверяемого члена; с3,- и

kn3j

характеристики

третьей

подсистемы.

 

 

 

 

При проверке должны выполняться условия:

 

 

 

 

 

 

К і ] " К К . огр>

 

 

 

(4.33)

 

 

 

Cfj

С’огр-

 

 

 

 

 

 

 

Табл ица 4.19

 

 

 

 

 

 

Проверка членов промежуточной оптимальной последовательности

 

 

с членами последовательности третьей подсистемы

 

 

 

 

П1, 2

 

С г*пз

0 , 0 0 8

0 , 0 0 1 6

0 , 0

0 0 3 2

0 , 0

0 0 0 6 4

 

 

 

4 , 5

6

7

, 5

 

9

С 1, 2

 

 

 

 

 

 

 

 

 

0,0097

 

 

 

 

 

0,009764

 

20

 

 

 

 

 

 

29

 

 

 

 

 

 

 

 

1

 

0,0084

 

 

 

0,00872

0,008464

 

22

 

 

 

29,5

 

31

 

 

 

 

 

 

2

 

3

 

0,0011

 

 

 

 

 

 

 

 

28

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2 7

МЕТОДЫ ПОВЫШЕНИЯ готовности

П ервая, вторая и третья комбинации не удовлетворяю т первому

условию . Ч етвертая комбинация { ° ’°2 9 76} удовлетворяет обоим

условиям . Эту комбинацию заносим в левый боковой столбец сле­ дующей матрицы (табл. 4.20) и повторяем процесс проверки. Ч етвер­ тая и пятая комбинации в этой матрице удовлетворяю т обоим ука-

I 0,0097 \

занным условиям и, следовательно, член | 20 1 утверждается

в качестве первого члена оптимальной последовательности.

Т а б л и ц а 4.20

Проверка членов промежуточной оптимальной последовательности с членами последовательности четвертой подсистемы

 

* П 4

0 , 0 0 2 - 1 3

0 , 0 0 0 7 3

0 , 0 0 0 2 2

0 , 0 0 0 0 6 6

 

0 , 0 0 8 1

* П 1 - 3

4

5

6

7

8

С 1 — 3

 

 

 

 

 

0,00976

 

0,00998

0,00982

29

 

36

37

0,00872

0,00945

0,00894

 

29,5

35,5

36,5

 

0,008464

0,00919

 

 

31

37

 

 

[ . . .

В первой матрице (см. табл. 4.18) по алгоритму Кеттела определяем следующий член оптимальной последовательности. Подвергаем его аналогичной проверке. Дальнейш ее рассмотрение членов оптималь-

нои последовательности приведет нас к члену ( 28 |> проверка которого показывает, что ни одна из его комбинаций в матрице с членами, относящимися к третьей и четвертой подсистемам, не удовлетворяет условиям (4.33) и, значит, этот член не должен быть

занесен

в оптимальную

последовательность. Д л я оставш егося члена

Г Г }

тем более

не

будут удовлетворяться условия (4.33), и он

такж е

не заносится

в

оптимальную последовательность.

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

В табл. 4.21 приведена сравнительная характеристикарешения задач по алгоритму Кеттела с осуществлением проверки и без нее при различных значениях глубины проверки s.

1 2 8

ГЛАВА 4

 

 

 

 

Таблица

4.21

Сравнительная характеристика решения задач по алгоритму Кеттела

 

 

 

 

 

Максималь­

Сумма длин

Время реше­

Количество

Способ вычисления

ная длина

всех опти­

ния задачи

подсистем

оптимальной

мальных

на ЭВМ

 

 

 

последова­

последова­

«Минск

22»,

 

 

 

тельности

тельностей

м ин

 

 

Без проверки

56

268

1

 

8

С осуществлением

про-

 

 

 

верки:

1

30

154

1

 

 

s =

 

 

s = 6

24

125

1

 

 

Без проверки

69

419

1,5

10

С осуществлением

про-

 

 

 

верки:

1

39

237

1

 

 

s =

 

 

s =

8

26

177

10

 

 

Без проверки

227

2134

10

 

16

С осуществлением

про-

 

 

 

верки:

1

119

1295

5

 

 

s =

 

 

s= 8

67

385

6

 

ВЫБОР ОПТИМАЛЬНОГО ЧИСЛА

 

 

§ 4.6

КОНТРОЛИРУЕМЫХ ПАРАМЕТРОВ

 

 

 

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

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

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

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

9 А. Г. Варжапетян

129

МЕТОДЫ ПОВЫШЕНИЯ г о т о в н о с т и

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

Воспользуем ся информационным критерием объективности кон­ троля, под которым будем понимать отношение количества инфор­ мации /к, приобретенной при контроле определенного числа пара­ метров п, к энтропии системы # 0:

Д ля случая

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

 

 

 

П

 

 

В =

S Ьі>

(4.34)

 

 

/=1

 

д е bl — Ii/Ho,

/,• — количество

информации,

приобретенной при

контроле і-го параметра.

 

контроля D и опре­

Следуя [1 0 ],

введем понятие достоверности

делим его как произведение критерия объективности контроля на вероятность безотказной работы автоматической системы контроля Р:

D - ВР.

(4.35)

Д ля решения поставленной задачи необходимо максимизировать функцию D при заданных ограничениях на время проведения кон­ троля /к, на стоимость системы контроля С, на ее вес G и т. д.

Обозначим через Я. интенсивность отказов части автоматической системы контроля, необходимой для контроля t-го параметра. При этом рассмотрим два режима работы автоматической системы кон­ троля:

1) все части системы работают непрерывно в течение всего вре­

мени

контроля;

 

 

 

 

 

2)

отдельные части

системы

вклю чаю тся только

на время кон­

троля

определенного

параметра.

 

 

 

Введем переменные xt таким образом, что xt — 1, если і-й пара­

метр

принят к контролю, и Х[ =

0, если он не принят. Тогда

 

 

 

П

 

 

 

 

 

д =

2

ьіхі;

 

(4.36)

 

 

 

1 = 1

 

 

 

 

 

 

 

первого

режима

работы;

 

 

 

 

 

 

(4.37)

 

 

для

второго

режима

работы.

1 3 0

ГЛАВА 4

П одставляя (4.36) и (4.37) в (4.35), получим значения достовер­ ности контроля для первого и второго режимов работы соответ­ ственно:

О

Ограничения, имеют вид

п

2

(= 1

п

|]

(

\ х і

V

 

 

ехр

t У

) \

 

 

 

\

<•=1

 

і=і

(4.38)

 

 

п

 

\ л

п

 

 

— 2

%itlxl

2 Ьіхі-

 

 

 

i=i

 

/J i=i

 

накладываемые

на переменные х и х 2,

XПI

пП

ctX[ ^ С;

2 g,-

G,

(4.39)

1=1

І = 1

 

 

где tt — время контроля і-го параметра; ct и g{ — стоимость и вес элементов автоматической системы контроля, обеспечивающих кон­ троль t-го параметра.

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

рования

[7 ].

 

 

 

 

 

 

 

 

 

 

Ф ункциональное уравнение для первого режима работы автома­

тической системы контроля имеет вид

 

 

 

 

 

 

 

 

 

fk (t, с,

g) = max ((& L i + h ) exp

[— t (Я,*_і + A *)]),

 

(4.40)

где A,*_p

bk*_ j —

интенсивность отказов

и значение критерия объек­

тивности

части

системы,

необходимой

для

контроля

А —

1

пара­

метров; А,*, bk — интенсивность отказов

и значение критерия объек­

тивности части системы, необходимой для

контроля А-го параметра;

t — время, необходимое

для контроля

А

параметров.

 

 

 

Функциональное уравнение для второго режима работы автома­

тической системы контроля можно записать в виде

 

 

 

fk (t, с, g) = max ((6fe_i -+- bk) exp [—

(^ _ i +

h)

(A,LI +

A *)]),

(4.41)

где A’ j ,

b\_v

fk_ x — интенсивность

отказов,

значение

критерия

объективности и

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

контроля k — -1

параметров;

6А, 4 — интенсивность отказов,

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

Д л я

иллюстрации рассмотрим численный пример [1 0 ].

П ример 4.3. Пусть имеется пять выходных параметров радиоэлектронной аппа­

ратуры 2 і,

гг, z3, z4, z5, характеристики системы контроля которых даны в табл. 4.22.

Требуется выбрать для контроля набор таких параметров, чтобы функция (4.38)

имела максимум, а вес системы контроля и время контроля в относительных едини­ цах (о. е.) не превысили значений G = 38, (к = 18.

9*

13 1

МЕТОДЫ ПОВЫШЕНИЯ готовности

Т а б л и ц а 4.22 Характеристики системы контроля выходных параметров________________________

Контролируемые выходные

 

 

 

и

ьі

параметры

 

 

 

гі

0,08

2

5

0,055

г2

0

, 0

1

4

8

0,0335

г3

0

, 0

2

5

6

0,0704

Z4

0,03

1

5

0,136

26

0

, 1

 

8

2 0

0,428

Предположим, что стоимость не является лимитирующим параметром. Решение задачи сводится к нахождению мажорирующей последовательности согласно функ­ циональным уравнениям (4.40), (4.41). Табл. 4.23—4.26 отображают процесс выбора

Т а б л и ц а 4.23

Получение членов промежуточной оптимальной последовательности при рассмотрении параметров гх и г2

О,

іг

0

0,0469

0

2

0

5

0

0

0,0469

0

0

2

0

0

5

 

 

1

0,0322

0,0322

0,0724

4

4

6

8

8

13

 

 

2

Т а б л и ц а 4.24

Получение членов промежуточной оптимальной последовательности при рассмотрении параметров гг—г3____________________________________________

1 , а

0

*lf 2 '

0

0

0

0

0469

0,0724

■>

6

5

13

0,0469

0,0724

0

0

2

6

0

0

5

13

 

 

1

3

0,0637

0637

0967

1177

5

5

7

1 1

6

6

И

19

 

2

4

5

1 3 2

ГЛАВА 4

Т а б л и ц а 4.25

Получение членов промежуточной оптимальной последовательности при рассмотрении параметров гг— z4

Di—3

Dt

Л - Э

0

0 , 0 4 6 9

0 , 0 6 3 7

0 , 0 7 2 4

0 , 0 9 6 7

0 , 1 1 7 7

 

й і - з

0

2

5

6

7

1 1

14

 

0

5

6

1 3

1 1

1 9

 

 

Üi

 

 

 

 

 

 

 

 

0

0

0,0469

0,0637

0,0724

0,0967

0,1177

 

0

0

2

5

6

7

1 1

 

0

0

5

6

13

11

19

 

1320

1320

1579

0,1812

1784

1956

2120

 

1

1

3

6

7

8

12

 

5

5

10

11

18

16

24

 

 

1

2

3

 

4

5

Т а б л и ц а 4.26

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

Di-i

\ u -t - gl-1

h Въ

0

0

0

0

0 , 1 3 2 0

0 , 1 5 7 9

0 , 1 8 1 2

0 , 1 9 5 6

0 , 2 1 2 0

0

1

3

6

8

1 2

0

5

1 0

1 1

1 6

2

4

0

0,1320

0,1579

0,1812

0,1956

0,2120

0

1

3

6

8

12

0

5

10

11

16

24

 

 

1

2

3

4

 

5

0,1923

0,1923

0,2459

0,2300

0,2503

0,2317

 

 

8

8

9

11

14

16

20

"

20

20

25

30

31

36

44

 

 

 

6

 

7

 

 

 

МЕТОДЫ ПОВЫШЕНИЯ г о т о в н о с т и

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

Для определения оптимального числа контролируемых параметров необхолимо произвести анализ всех мажорирующих последовательностей, начиная с последней (четвертой). Находят член этой последовательности с максимальной достоверностью и проверяют, не входит ли он в состав предыдущих последовательностей. Если не входит, то параметр гъ подвергается контролю. Анализ третьей мажорирующей после­ довательности начинают с члена, имеющего характеристики 0 = 0,1812; / = 6; g = 11, так как с помощью этих характеристик и характеристик, соответствующих параметру z5, образован член с максимальной достоверностью четвертой последова­ тельности. Найденный член третьей мажорирующей последовательности подвер­ гают проверке, аналогичной проводившейся для члена четвертой последовательности

Т а б л и ц а 4.27

Получение членов промежуточной оптимальной последовательности при рассмотрении параметров гх и z2

Т а б л и ц а 4.28

Получение членов промежуточной оптимальной последовательности при рассмотрении параметров zx—z3

\Д1,2

N .

1, 2

0

0,0-169

0,0515

N.

gl. 2

0

2

6

Dr, \

N.

0

5

13

Із

 

 

 

Яэ

N^

 

 

 

 

 

0

0

0,0469

0,0515

 

 

0

0

2

6

 

 

0

0

5

13

 

 

 

 

1

 

 

 

0,0637

0,0637

0,0623

0,0476

 

 

5

5

7

11

 

 

6

6

11

19

 

 

 

2

 

 

 

 

 

 

Т а б л и ц а 4.29

Получение членов промежуточной оптимальной последовательности

 

при рассмотрении параметров гх—z4

 

 

 

 

 

D i —з

 

 

 

 

 

 

0

0 , 0 4 6 9

0 , 0 5 1 5

 

S 1 - 3

0

2

 

6

 

0

5

 

1 3

Ëi

~~~~___

 

 

 

 

 

0

0

0,0469

0,0515

 

0

0

2

 

6

 

0

0

5

 

13

 

0,1320

0,1320

0,1373

0,1529

 

1

1

3

 

6

 

5

5

10

 

11

 

 

1

2

 

3

134

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