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

книги из ГПНТБ / Судовые системы автоматического контроля (системный подход к проектированию)

..pdf
Скачиваний:
13
Добавлен:
21.10.2023
Размер:
9.77 Mб
Скачать

Из = 0 ,1 5 ; р4 =

0,35; р5 = 0,3; р6 = 0,1;

р7 = 0,05; р8= 0,02.

По значениям

весовых коэффициентов

построена гистограмма

(рис. 1.4), аппроксимируемая нормальным распределением с пара­ метрами М [Гер] = 600 ч, о [Гср] = 130 ч. Проверка этой гипо­ тезы может быть осуществлена известными методами математической статистики, например, с помощью критерия %2.

Как уже отмечалось, одной из важных задач информационного

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

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

средств судна, или САК на ранних этапах проектирования. В пре­

дыдущем параграфе рассматривались аналитические методы прогно­ зирования. Однако эти методы

дают

достоверные

результаты

р

 

/

 

 

при большом количестве систем-

 

 

 

прототипов. Если это требова­

 

 

 

 

ние не выполняется, для прогно­

0,3

М *

/ \М 5

 

зирования могут быть исполь­

 

 

1Л'

 

 

зованы

только

различные

ме­

 

\

 

 

 

тоды

экспертной оценки.

 

0,2

L

__

 

 

/

\

 

В

 

общем случае

прогнози­

 

 

руемое

значение

показателя

0,1

М з /

 

\

Мб

представляет

собой

случайную

 

/

 

\

Я17

величину х,

имеющую некото­

 

 

рую

функцию

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

 

М 2

 

 

Ms

 

200 т

600

800 t 1000

f {х).

Каждый

показатель

как

 

 

 

 

 

 

случайная величина имеет свое

 

Рис. 1.4.

 

 

специфическое

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

 

 

 

 

 

Задача

состоит

в построении такого

обобщенного распределения,

которое отражало бы наиболее

общие закономерности и могло слу­

жить

исходной

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

лей системы.

 

 

 

 

 

 

 

 

 

При выборе формы распределения необходимо учитывать следу­

ющие

 

соображения,

вытекающие из

особенностей

судовых

систем:

1) значения показателей ограничены сверху и снизу и, если рассматривать их как случайные величины, имеют ограниченную по абсциссе непрерывную, одномодальную функцию распределения; 2) значение показателя зависит от большого числа случайных факторов, каждый из которых в отдельности оказывает несуще­

ственное влияние. Вместе с тем, имеются факторы, число которых

невелико, а влияние существенно. Именно эти факторы приводят

к различию значений одноименных показателей в разных вариантах реализации системы. Влияние этих факторов вызывает асимметрию распределения.

Таким образом, типовое распределение значений показателей

судовых систем должно быть как правило непрерывным, одномо­

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

Среди практически применяемых распределений такого класса наи­

большей гибкостью обладает

бета-распределение [34, 57], которое

и целесообразно принять в

качестве типового.

31

Плотность бета-распрёделёния выражается следующим образом:

 

(1.4.1)

f (х) = 0

при х < а, х > Ь,

где В (р, q) — бета-функция; a,

b, р, q — параметры распределения,

причем а и b определяют левую и правую границы распределения

соответственно.

На рис. 1.5 представлены кривые бета-распределения при раз­ личных р и q. Из рисунка видно, что, меняя эти параметры, можно

получить практически все одномодальные распределения (при р ^ 1

и <7 1).

Выражения для математического ожидания М [х], моды Мо [х\

и дисперсии D [х] соответственно

имеют вид

 

 

 

М [х] = а ф a) j ^

,

 

 

Мо [х] — а -\- (Ь а) p + q 2

(1.4.2)

D

[х] ~ (b — а)2

PQ

 

 

 

 

(Р +

<7)2(Р +

Ч +

О '

 

Эти выражения трудно использовать для решения рассматривае­

мой задачи, так как параметры р и q бета-распределения не имеют

ясной физической интерпретации. Поэтому проделаем некоторые преобразования.

Рассматривая выражения (1.4.2) как систему уравнений и исклю­

чая р и q, получаем

М 3

[х] (a -f Ь +

Мо [х]) М г [х] + (ab +

а Мо [х]

+

 

+

b Mo [х ] + D [х ] М [х ] + (aD

[х ] +

 

 

+ bD

\х] — ab Мо [х] — 3D [х] Мо [х] = 0 .

(1.4.3)

Для

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

расчетов целесообразно

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

это выражение

более

простым

 

 

 

 

М [х] ~ - - ± ^ М° [2Х| + b = О,

 

(1.4.4)

Из математической статистики [34, 57] известно, что для одно­ модальных распределений можно полагать

D

[*] = &(& — а)2.

(1.4.5)

Учитывая это, получим

 

 

(p +

q)4p + q + l)

= k = const.

(1.4.6)

 

 

32

Тогда

1

2.

(1.4.7)

6k

 

 

Для определения величины k можно использовать то известное из математической статистики обстоятельство [19], что по крайней

мере 89% площади под кривой плотности любого распределения

лежит внутри трех стандартных отклонений от среднего. Основы­

ваясь на этом, можно считать k = 1/36 [19]. Тогда

М м = Ж М о Ш ^ )

(1.4.8)

D[x] = ~ { b — a f.

Следовательно, чтобы найти значения математического ожидания

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

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

нить

граничные

значения

оцениваемого показателя

и Ь)

и

его

наиболее

ве­

роятное значение (моду Мо),

то,

используя

выражения

(1.4.8), можно вычислить ма­

тематическое ожидание и ди­

сперсию

значения

показа­

теля

х.

 

 

 

 

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

выраже­

ний (1.4.8) требует от экспер­

тов оценки значений трех

то­

чек

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

Обычно

для

вновь разрабатываемых

систем

оценка

моды некоторых показателей вызывает большие

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

нима при наличии информации только о его граничных значениях.

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

/ (х) = С (х а) (b — х)2.

(1.4.9)

Константа С в этом выражении определяется из условия нор­ мировки

ь

| f(x)dx— 1.

а3

3Зазка797

3 3

 

Тогда плотность распределения величины х в границах а и Ь

f(x) - {Ь^ ауг(х - а)(ь ~ * ) 2-

(1.4.10)

Для такого распределения математическое ожидание, мода и

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

равны

 

, .

г ,

За +

 

м

[X] =

— ± ----,

 

Мо [лг] =

--- g— ,

(1.4.11)

 

D [х] = 0,04 — а)2.

Для использования этих выражений от экспертов требуется

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

Рассмотрим пример использования последнего метода. В процессе проектирования судовой САК были предложены два варианта, отличающиеся структурой, методом отработки информации, номен­

клатурой блоков и загрузкой оператора.

На основе анализа особенностей контролируемых систем и тре­

бований к САК

были выбраны следующие основные показатели,

по совокупности

которых можно сравнивать возможные варианты:

1)надежность САК, оцениваемая как вероятность безотказной

работы за заданное время (х4);

2)универсальность, оцениваемая как процент основных узлов

САК, пригодных без переделок для проектируемых судов (х2);

3)достоверность информации, оцениваемая вероятностью пра­

вильного обнаружения неисправности (х3);

4)время обнаружения места неисправности в контролируемых системах (х4).

Остальные требования к системе, в том числе стоимость, были

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

В техническом задании были оговорены следующие значения

перечисленных показателей: хг ^ 0,95; х 2 ^ 0,9; х3 ^ 0,8; х4 ^ «с Ю с.

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

ково. По материалам предварительных разработок САК, с учетом сведений о контролируемых системах, эксперты дали оценки

показателей, представленные в табл. 1.1.

Для определения значения математического ожидания i-ro по­ казателя (i = 1, 2, 3, 4) по оценкам /'-го эксперта (/ = 1, 2, 3, 4)

использовалось выражение

М[х] = — iimax + 3Xi/mln .

34

 

 

 

 

 

Т абли ца 1.1

 

 

В а р и а н т I

В а р и а н т II

Э к с п е р т

П о к а з а т е л ь

м а к с и м а л ь ­

м и н и м а л ь н о е

м а к с и м а л ь ­

м и н и м а л ь н о е

 

 

 

 

н о е з н а ч е н и е

з н а ч е н и е

н о е з н а ч е н и е

з н а ч е н и е

 

 

0,98

0,95

0,96

0,94

 

* 2

0,95

0,9

0,9

0,85

 

х3

0,9

0,75

0,9

0 ,8

 

* 4

9,0

7,0

9,0

7,0

о

*1

0,97

0,94

0,98

0,96

* 2

0,95

0,85

0,95

0,9

 

х3

0,9

0 ,8

0 ,8

0,75

 

ч

8,5

8 ,0

9,0

7,5

Q

X1

0,97

0,93

0,98

0,94

х2

0,98

0,95

0,95

0,9

 

х3

0,85

0 ,8

0,85

0 ,8

 

xi

9,0

8 ,0

9,0

7,0

 

Xl

0,96

0,93

0,98

0,95

4

х3

0,98

0,96

0,95

0,9

 

0,9

0,7

0,9

0,7

 

*4

1 0 ,0

9,0

9,5

8 ,0

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

приведены в табл. 1.2. Осредняя оценки экспертов, получаем сле­ дующие значения показателей (математические ожидания):

В а р и а н т

П о к а з а т е л ь

И

Xl

0,95

0,96

хг

0,92

0,91

* 3

0,81

0 ,8

ХА

8,14

7,95

Часто при использовании метода экспертных оценок исходят из предположения, что искомое значение показателя является сред­ ним на интервале [а, b ], т. е. распределение значений показателя принимается равновероятным. Как известно [34], для этого распре­

деления оценки математического ожидания и дисперсии опреде­

ляются выражениями

М[х] = ± ± ± , D [ x ] = - ^ ±

3!

35

 

 

 

Табли ца 1.2

 

 

 

В а р и а н т

Э к с п е р т

. П о к а з а т е л ь

 

 

 

 

I

II

 

X 1

0,96

0,96

 

*2

0,92

0,87

 

х 3

0,81

0,84

 

* 4

7,8

7,8

 

* 1

0,95

0,97

 

Х2

0,89

0,92

 

xl

0,84

0,77

 

8 ,2

8,1

 

Х1

0,95

0,96

q

х 2

0,96

0,92

 

хз

0,82

0,82

 

.

8,4

7,8

 

 

 

 

Х1

0,94

0,96

4

Х2

0,94

0,92

 

х3

0,78

0,78

 

* 4

8,16

8,06

Сравнивая эти выражения с (1.4.11), видим, что оценка диспер­

сии при использовании гипотезы о бета-распределении значений

показателя л: получается почти вдвое меньше, а оценка математи­

ческого ожидания меньше на величину 0,1 а). Таким образом, по выражениям (1.4.11) получается более осторожная оценка и с меньшей дисперсией, что является важным положительным фак­ тором.

При опросе N экспертов для оценки показателя х в качестве их коллективной оценки принимается составное распределение

N

 

Ы *) = S в Л М ;

(1.4.12)

 

/=1

 

 

 

N

 

 

 

§1*1 =

1,

(1.4.13)

где /г (х)

плотность распределения,

построенная

по оценкам t-ro

эксперта;

цг— «вес» t-ro эксперта (оценка его квалификации).

Выражения для математического ожидания и дисперсии вели­

чины х в случае использования

оценок N экспертов

имеют вид

М^[х] = § ц М [х];

N

(1.4.14)

Ds [х] = 2 p D [Х],

i—i

t=i

'

36

где Mi lx] и Dt [л:] — соответственно математическое ожидание и дисперсия распределения, полученного по оценкам г'-го эксперта.

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

циент р(- для каждого t-го эксперта. При этом можно считать, что

дисперсия оценки, даваемой экспертом, обратно пропорциональна

его квалификации р'. Тогда

щ : р 2:

: ^

=

~Ж\х] :

'■ ■ ■ ■ '• ТЩ /\ ■

О-4-15)

Пусть р(, р',

. . ., р^

положительные числа, удовлетворяю­

щие (1.4.15). Чтобы учесть нормирующее условие (1.4.13), доста­

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

2 > ;

i=\

 

Из (1.4.15) следует, что p iD a [х] = • • • = рд,DN

[х] — D, где

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

полученного

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

получить оценку Мг [х] с наименьшей дисперсией.

В ряде случаев вместо выражений (1.4.14) удобнее пользоваться

усредненными значениями оценок границ распределения а и Ь. Тогда при использовании мнений N экспертов для оценки математи­

ческого ожидания и дисперсии получаем соотношения:

 

 

N

 

N

 

 

 

 

3 2

М-Шг + 2

Рibi

 

 

 

М2 [х] = -

^

------------------ ,

 

 

 

/

N

N

\ 2

 

 

Ds [X] = 0,04

2 РА- 2 j V-fii

>

 

 

\ i = 1

1=1

1

 

где cti,

bt — оценки границ

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

i

экспертом; р; —

весовые

коэффициенты экспертов.

 

 

 

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

устройств регистрации судовой САК.

В качестве характеристики надежности принято среднее время

исправной работы между отказами ГСр. Была сформирована группа

из пяти специалистов — работников НИИ, КБ и пароходства, ко­ торых ознакомили со всеми предполагаемыми отличиями разрабаты­

ваемых вариантов от известных. Им предложили оценить предпола­ гаемые границы значений Тср, минимальное /гоШ и максималь­

37

ное tm!a. По выражениям (1.4.11) определены значения математиче­

ского ожидания и дисперсии среднего времени исправной работы. Результаты приведены в следующей таблице:

Эксперт

^mln*

^тах*

м [Ар]’

D [Ар]-

 

ч

ч

Ч

Ч

 

I

вариант

 

 

1

5500

6000

5700

10

2

5500

6500

5900

40

3

6000

6500

6200

10

4

5000

6000

5400

40

5

6000

6500

6200

10

 

 

II вариант

 

 

1

6000

6500

6200

10

2

6000

7000

6400

40

3

6500

7000

6700

10

4

6000

7500

6600

60

5

6500

7500

6900

40

Используя выражения (1.4.14) в предположении, что квалифика­

ция экспертов одинакова, получим средние результаты соответственно для I и II вариантов:

M iz[Т’ср] =

5880

ч,

Dlz [Tcp} =

22

ч;

M2 z [Тср] =

6560

ч,

D2Z [Тср] =

32

ч.

Разница в оценках математического ожидания Тср вариантов составляет ~ 1 0 % , что весьма существенно для дальнейшего сравне­ ния вариантов САК по совокупности показателей.

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

практики известно, что именно такие сведения эксперты дают наи­

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

нием, получить значения математического ожидания и дисперсии.

Вкачестве аппроксимирующего распределения можно использовать

нормальное распределение. Объясняется это тем, что процесс оценки

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

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

38

факторов И, согласно центральной предельной теореме теории вероят­

ностей [19], погрешность измерения является случайной величиной,

имеющей нормальное распределение.

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

тоды же определения характеристик законов распределения на ос­

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

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

нее дают информацию о границах значений показателей, а не точеч­

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

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

ляет интуитивную информацию, имеющуюся у специалистов и осно­

ванную на их знаниях и опыте, выразить в вероятностных числовых

оценках.

§ 1.6

ВОПРОСЫ ДОСТОВЕРНОСТИ ИНФОРМАЦИОННОГО ОБЕСПЕЧЕНИЯ

Достоверность информации, принимаемой в качестве исходной

при разработке судовых САК, во многом определяет ее эффектив­

ность, сроки и стоимость проектирования и освоения. Поэтому оценка

и повышение достоверности исходной информации являются важ ­

нейшей задачей при проектировании судовых систем.

Способы оценки и повышения достоверности информации зависят

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

систем недостаточность информации при априорной оценке показате­

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

при проектировании, получается путем статистической обработки

результатов испытаний и эксплуатации. К ним относятся сведения

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

сведения о влиянии САК на технико-эксплуатационные характери­

стики контролируемого оборудования и судна в целом. Для возмож­

ности использования статистических данных об этих показателях

при проектировании необходимы их обобщенные статистические

оценки в виде законов распределения, средних значений, дисперсий, показателей корреляции и др. Упоминавшиеся уже особенности

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

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

виям состоятельности, эффективности, достаточности и несмещен­

ности [57]. Имеющиеся отдельные реализации в этом смысле яв­

ляются малоинформативными за исключением, конечно, экстра­

39

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

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

любыми методами.

Особый интерес представляют вопросы оценки и повышения до­

стоверности информации, получаемой эвристическими методами. При

использовании этих методов появляются дополнительные источники

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

группы.

1. Случайные ошибки у. Это типичные ошибки любых измерений

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

N дающих согласованные оценки экспертов и

усредненная при

N —> о о , является достоверной. Поскольку всегда

возможен опрос

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

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

(школы) и больше всего сказываются при оценках показателей, не

имеющих общепринятых единиц измерения.

3.Психологические ошибки ф. Специфические при использовании

эвристических методов они вызываются как чисто психологическими

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

оценки зависит от желания эксперта дать добросовестную оценку

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

наоборот, занижены у разработок организаций с низкой репута­

цией. Возможно стремление эксперта не давать крайних оценок

(ошибка центральной тенденции) или давать одинаковые оценки ло­

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

Средней ошибкой оценки при использовании эвристических мето­ дов является случайная величина

б = М [у] ± М [е] ± М [ф],

где М [у], М [е] и М [ф] — математические ожидания случайной,

систематической и психологической ошибок соответственно.

40

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