Розробляеться nidxid до побудови ModeMi ощнки ймов1рност{ виконання алгоритмгв функщонування дискретних систем "людина-машина"зр1знимир1внямиякост1. Розвиток видомихратше тдход1в полягае в розширент можливостей типових функщональних моделей i видмову в1д традищйног ттерпретащг поняття "помилка" i "безпомилков^ть". Працездаттсть тдходу демонструеться на nрикладi типовог функщональног структури "робоча-контроль функщонування"
Ключoвi слова: функщонально-структур-
на тeoрiя, рiвнi якoстi
□-□
Разрабатывается подход к построению модели оценки вероятности выполнения алгоритмов функционирования дискретных систем "человек-машина" с различными уровнями качества. Развитие известных ранее подходов состоит в расширении возможностей типовых функциональных моделей и отказе от традиционной интерпретации понятия "ошибка" и "безошибочность". Работоспособность подхода демонстрируется на примере типовой функциональной структуры "рабочая- контроль функционирования"
Ключевые слова: функционально-
структурная теория, уровни качества □-□
Developing an approach to building a model estimating the probability that the functional algorithms of discrete systems "man-machine" with different levels of quality. The development of the previously known approaches lies in the empowerment model functional models and the rejection of traditional interpretation conception of "error" and "infallibility". The efficiency of the approach is demonstrated by the example of the type functional structure "working control operation"
Keywords: functional-structural theory, levels of quality
УДК 658.5.011.56
подход к
вероятностной
оценке качества результатов функционирования систем «человек -машина»
Е.А. Лавров
Доктор технических наук, профессор Кафедра информационных систем в менеджменте Национальный университет биоресурсов и природопользования Украины г. Киев
Контактный тел.: 050-691-37-33 E-mail: [email protected]
Н.Л. Барченко
Аспирант
Кафедра кибернетики и информатики Сумской национальный университет г. Сумы, ул. Кирова, 160, 40020 Контактный тел.: 099-975-27-11 E-mail: [email protected]
1. Введение
Функционально-структурная теория (ФСТ) эрго-технических систем [1], разработанная школой проф. А. Губинского, предлагает ряд комплексных моделей моделирования процессов возникновения и устранения ошибок. В данных моделях принимается бинарная модель ошибки "есть ошибка - нет ошибки". Возможность наличия ошибок разного типа промоделирована в работах [2, 3]. Эти работы позволяют определить вероятность возникновения ошибок разных типов. Однако, часто задача моделирования результатов функ-
ционирования человеко-машинных систем (ЧМС) ставится таким образом, что необходимо определить не вероятность ошибки, а вероятность достижения того или иного уровня качества qi е й^сц,..^), где п-количество выделенных уровней качества. Примерами таких задач могут быть задача определения уровня качества выпускаемых изделий (1 сорт, 2 сорт, 3 сорт и брак) или задача прогнозирования оценки студента в диалоговой системе "студент-компьютер", где качество может быть описано следующими уровнями: "отлично", "хорошо", "удовлетворительно", "неудовлетворительно".
Целью настоящей работы является исследование возможности распространения методологии ФСТ проф. А. Губинского на задачу моделирования различных уровней качества. Для примера покажем возможность получения зависимости уровней выходного качества от характеристик выполнения операций процесса функционирования и характеристик контроля для одной из типовых функциональных структур (ТФС).
2. Постановка задач исследования
1. Известна структура взаимосвязи операций в алгоритме деятельности (например, ТФС "рабочая с контролем функционирования").
2. Для рабочей операции известны вероятности
её выполнения с разными уровнями качества: В>- ве-
роятность того, что рабочая операция выполнена с ]-ым уровнем качества, ] = 1,4 ; для операции контроля
функционирования К1 - вероятность того, что опера-
ция, выполненная с 1-м уровнем качества, признана
выполненной с ,-м уровнем качества, ] = 1,4 ,1 = 1,4 .
Необходимо оценить вероятность выполнения всей структуры с 1-м уровнем качества V; , 1 = 1,4 и математическое ожидание времени выполнения М(Т).
3. Разработка подхода к выводу расчётных зависимостей
3.1. Методология вывода расчетных зависимостей, принятая в ФСТ
В основе получения конечного набора формул для множества ТФС - процедура перехода от "графа работ", соответствующего ТФС, к "графу событий" [1]. Информация, необходимая и достаточная для построения моделей оценки АФ ЧМС, может быть задана на вероятностном графе G(X, Z), вершины которого (X - множество вершин) соответствуют событиям начала и окончания выполняемых операций, а дуги - возможным переходам ^ - множество дуг).
Задача анализа алгоритма функционирования ЧМС, заданного на вероятностном графе G(X, Z), сводится к преобразованию типа
G(X, Z) ^ G(X*, Z*),
где G(X*, Z*) - укрупненный вероятностный граф, в котором промежуточные вершины отсутствуют; веса дуг, направленных из начальной вершины, соответствующей принятию решения о начале процесса функционирования, в поглощающие вершины, представлены в виде функциональных зависимостей от весов дуг исходного графа G(X, Z), что позволяет вычислять вероятности и затраты, связанные с переходами из начальной в поглощающие вершины.
В основе метода укрупнения вероятностных графов лежит понятие Н-функции [2]. Пусть G(X, Z) - вероятностный граф. Каждую его дугу (;, ^ е Z, ; е X, ] е X взвесим значением ненулевой вероятности р, перехода ; ^ ^ которая соответствует уровню правильности (безошибочности) выполнения некоторого элемента процесса функционирования. При этом ^ р, = 1 для всех ; е X. Каждой дуге (;, ^ е Z поставим в соответствие случайную величину ^ затрат (времени, стоимости), связанных с переходом ; ^ . Величина ^ характери-
зуется функцией распределения Fij(t) с плотностью ^О, математическим ожиданием Му и дисперсией Dij. Предполагается, что случайные величины ^ являются независимыми для всех дуг (;, .) е Z.
Н - функцией дуги (;, .) е Z назовем величину, определяемую следующим образом:
Ну^) = р^у ( t ) ],
(1)
где L - условное обозначение применения оператора Лапласа к функции распределения
Fij(t): L[Fij(t)]= |е"8^)
Зная Н-функцию, получаем рц = Ну (8) | 8=0 ;
р^ ds
1 Н2 1
= —L"1[Hij(s)], . Рп j
(2)
(3)
(4)
(5)
где L -1 - условное обозначение обратного преобразования Лапласа.
3.2. Разработка модели для ТФС "рабочая - контроль функционирования"
Для примера возьмем ТФС "рабочая - контроль функционирования" без ограничения на допустимое число т циклов контроля (т.е. т = ~ ), схема уменьшения размерности графа и изменения матриц Н - функций для бинарных исходов операций рассмотрены в [1]. Граф-схема указанной ТФС изображена на рис. 1.
Рис. 1. ТФС «Рабочая — контроль функционирования».
Граф работ
Пусть рабочая операция может быть выполнена с 1, 2, 3 уровнем качества или с браком. Процесс функционирования и, естественно, результат, будут зависеть от установки на уровень качества. Имеется в виду минимально допустимый уровень качества, при котором не осуществляется возврат на повторное выполнение рабочей операции. В данном случае возможны три установки: а) установка на качество 1; б) установка на качество 1 или 2; в) установка на качество 1, 2 или 3. Чем ниже установка, тем меньше затраты на соответствующий контроль и меньше вероятность повторения рабочей операции. Итоговые результаты будут зависеть от исходных данных и от установки на уровень качества. В данной работе ограничимся примером получения расчётной установки для случая в).
Для построения математической модели обозначим события, связанные с выполнением операций, входящих в ТФС. Вершины вероятностного графа (рис. 4) соответствуют следующим событиям:
1. Начало выполнения структуры (вершина 13).
2. Окончание выполнения рабочей операции с 1-м, 2-м, 3-м качеством или браком и переход к операции контроля (вершины 14-17).
3. Окончание выполнения структуры с 1-м, 2-м, 3-м качеством (вершины 1-12).
(6)
Рис. 2. ТФС «Рабочая — контроль функционирования». Граф событий
Н-функции дуг графа, соответственно формуле (1), имеют вид:
Н^в^^)], Н^в^^)], Н^^В3^^)], Н13,17=(1-(В1+ В2+ В3)^^)],
Н^^К1^^)], ^5,13=^^)], Н16д3=К3^^ка)], ^7,13=^^)], Н^К1^^)], Н!4,2=К^ка)], Н!4,3=К^^ка)],
Нl5,4=K22L[Fk(t)], Н15,5=К2^^)],
Н^К2^^)],
Нl6,7=K33L[Fk(t)], Нl6,8=K31L[Fk(t)], Нl6,9=K32L[Fk(t)],
^7,10=^^)], Нl7,ll=K42L[Fk(t)], Нl7,l2=K41L[Fk(t)].
Здесь Fp(t), Fk(t) - законы распределения времени на выполнение рабочей операции и контроль функционирования, соответственно.
Для вероятностного графа ТФС "рабочая - контроль функционирования" без ограничения на допустимое число циклов контроля, изображенного на рис. 2, матрица Н-функций имеет вид:
Н =
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 Н13.14 Н13.15 Н13.16 Н13.17
Н14.1 Н14.2 Н14.3 0 0 0 0 0 0 0 0 0 Н14.13 0 0 0 0
0 0 0 Н15.4 Н15.5 Н15.6 0 0 0 0 0 0 ^15.13 0 0 0 0
0 0 0 0 0 0 Н16.7 Н16.8 Н16.9 0 0 0 Н16.13 0 0 0 0
0 0 0 0 0 0 0 0 0 Н17.10 Н17.11 Н17.12 Н17.13 0 0 0 0
Размерность матрицы и графа может быть сокращена за счет транзитивных вершин. При удалении
Н =
графе образуется неограниченная ( т = ^ ) ветствует матрица:
1 0 0 0 0 0 0 0 0 0 0 0
0 1 0 0 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0 0 0 0
0 0 0 0 1 0 0 0 0 0 0 0
0 0 0 0 0 1 0 0 0 0 0 0
0 0 0 0 0 0 1 0 0 0 0 0
0 0 0 0 0 0 0 1 0 0 0 0
0 0 0 0 0 0 0 0 1 0 0 0
0 0 0 0 0 0 0 0 0 1 0 0
0 0 0 0 0 0 0 0 0 0 1 0
0 0 0 0 0 0 0 0 0 0 0 1
Н13.1 Н13.2 Н13.3 Н13.4 Н13.5 Н13.6 Н13.7 Н13.8 Н13.9 Н13.10 Н13.11 Н13.12
петля. В результате получим граф, имеющий в вершине четыре неограниченные петли (рис. 3). Графу соот-
0 0 0 0 0 0 0 0 0 0 0 0
+ Н" + Н:
+ Н1У
13.3 т п13.3
Рис. 3. Вероятностный граф после удаления транзитивных вершин
Дальнейшее укрупнение графа, изображённого на рис. 3 , сводится к замене равноценных петель на одну однородную петлю [1]:
Н13.13= Н13.14Н14.13+ Н 13.15 Н 14.15+
+ Н13.16Н16.13+ Н13.17Н17.13 (7)
Пользуясь правилом удаления неограниченной петли (т.е. бесконечное число повторения цикла), получим итоговый граф (рис. 4), для которого Н - функции будут иметь вид:
Рис. 4. Итоговый вероятностный граф
H13.1 = H13.1
H13.3 = H13.3
H13.5 = H13.5
H13.7 = H13.7
H13.9 = H13.9
H13.11 = H13.17
1
1- H13.13
1
1- " H13.13
1
1- " H13.13
1
1- " H13.13
1
1- " H13.13
1
> H13.2 = H13.2
1
1 - H
> H13.4 = H13.4
13.13 1
1-H
> H13.6 = H13.6
1-H
> H13.8 = H13.8
1-H
(8)
> H13.10 = H13.10
1-H
1-H
' > H13.12 = H13.17
13.13 1
1-H
43.13 1 ^^13.13
На основе полученных итоговых Н-функций (8) и формул (2), (3), (4), выражений для Н-функций исходного вероятностного графа (формулы 6 и 7), получим вероятностные характеристики путем подстановки параметра 8=0:
Используя полученные исходы, определим вероятности выполнения структуры с качеством 1, 2, 3:
^=^3.1+ V13.5.+V13.8+V13.12 - вероятность выполнения структуры с качеством 1
V2=V13.2+ V13.4+V13.7+V13.11- вероятность выполнения структуры с качеством 2
V3=V13.3+ V13.6+V13.9 + ^3.10 - вероятность выполнения структуры с качеством 3.
Для упрощения записи формул обозначим: Ь= =Н1313(0). Тогда формулы для определения вероятностных характеристик ТФС "рабочая - контроль функционирования" примут вид:
V! = (Б'К11 + В2К21 + В3К31 + (1 - (В1 + В2 + В3))К41), V, = (В'К12 + В2К22 +В3К32 + (1 - (В1 + В2 + В3))^2)^,
V3 =
= (B'K13 +B2K23 +B3K33 +(1 -(B1 + B2 + B3))K43)
(9)
1 - b
где Ь=В1К14+ В2К24+ Б3K34+(1-(Б1+Б2+Б3))K44
Математическое ожидание времени выполнения структуры определяется по формуле (3). В итоге получим:
. М(Тр) + М(Тк)
M = -
1 - b
(10)
V13.1 = V13.2=
V13.3= V13.4 = V13.5= V13.6= V13.7= V13.8= V13.9= V13.10: V13.11: V13.12
H*13.1(0) H*13.2(0) H*13.3(0) H*13.4(0) H*13.5(0) H*13.6(0) H*13.7(0) H*13.8(0) H*13.9(0) =H*13.10(0) =H*13.11(0) =H*13.12(0)
-вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того, -вероятность того,
что структура что структура что структура что структура что структура что структура что структура что структура что структура что структура что структура что структура
выполнена с выполнена с выполнена с выполнена с выполнена с выполнена с выполнена с выполнена с выполнена с выполнена с выполнена с выполнена с
Эти модели соответствуют установке на допустимость любого из трёх уровней качества (1, 2 или 3). Аналогичным образом получены формулы, по которым определяются вероятностные характеристики итогового графа, отвечающего ТФС "рабочая
- контроль функционирования" для уровней качества 2 и 1.
3.3. Компьютерное моделирование
Рассмотрим моделирование процесса взаимодействия в диалоговой системе. Имеются следующие уровни качества: 1- выполнение этапа диалога без ошибок, 2 - выполнение этапа диалога с незначительными ошибками, 3 -выполнение этапа диалога со значительными ошибками, 4 - невыполнение этапа диалога.
Структура алгоритма функционирования процесса обработки данных задается ТФС "рабочая
- контроль функционирования". С использованием полученных моделей (9) и (10), определим V; - вероятность того, что структура выполнена с 1-м уровнем качества и М(Т)-математическое ожидание
времени выполнения
качеством 1, признана с качеством 1 качеством 1, признана с качеством 2 качеством 1, признана с качеством 3 качеством 2, признана с качеством 2 качеством 2, признана с качеством 1 качеством 2, признана с качеством 3 качеством 3, признана с качеством 3 качеством 3, признана с качеством 1 качеством 3, признана с качеством 2 браком, признана с качеством 3 браком, признана с качеством 2 браком, признана с качеством 1
структуры для трех установок на уровень качества. Вычисления выполняем в среде MS Excel.
Пример. Пусть имеются несколько операторов с такими показателями качества выполнения операций "рабочая" (табл. 1) и «контроль функционирования» (табл. 2):
1
1
1
П Г.! 1иЛ7л'1Ш,
Таблица 1
Характеристики операции «рабочая»
Вероятность выполнения рабочей
Операторы операции с j-м качеством, В' (j = 1,4 ) М(Т3)
В1 В2 В3 В4
1 0,1100 0,5021 0,3781 0,0098
2 0,4470 0,3380 0,2100 0,0050
3 0,5220 0,3750 0,0941 0,0089 20
4 0,7310 0,2531 0,0110 0,0049
5 0,8940 0,0570 0,0430 0,0060
Таблица 2
Характеристики операции «контроль функционирования» (одинаковые для всех операторов'
К11 К22 К33 К41 К42 К43 К14 К24 К34
0,9741 0,9741 0,9741 0,0009 0,0013 0,0778 0,0129 0,0139 0,0149
К44 К12 К13 К21 К23 К31 К31 К32 МОк)
0,92 0,01 0,002 0,01 0,003 0,003 0,001 0,01 7,7
Результаты моделирования процесса функционирования с исходными данными для разных операторов представлены на рис. 5.
а) б)
Рис. 5. Результаты моделирования процесса функционирования (фрагмент) а) вероятность выполнения структуры с качеством 1; б) математическое ожидание времени выполнения структуры
4. направления дальнейших исследований
В последующих исследованиях необходимо: 1) получить расчётные зависимости для других ТФС; 2) разработать процедуру редукции функциональной сети таким образом, чтобы получать модели качества алгоритмов функционирования произвольной структуры.
5. Выводы
Продемонстрирована возможность расширения известного подхода к выводу расчётных зависимостей для моделирования безошибочности и времени выполнения ТФС, развитого в [2,3], на случай, когда необходимо определить вероятности достижения различных уровеней качества завершения алгоритмического про-
цесса. Подход может быть использован как в промышленности, АПК, так и, например, при прогнозировании результатов обучения в диалоговых системах.
Литература
1. Губинский А.И. Надежность и качество функционирования эргатических систем. Л.: Наука, 1982.- 270с.
2. Козачко О.М. Моделювання надшност алгоритмiчних процеав, як виконуються з помилками рiзних тишв, Автореферат на здобуття наукового ступеня кандидата техшчних наук, Вшниця - 2006. - 19 с.
3. Лавров Е.А., Пасько Н.Б. моделирование надёжности че-
ловеко-машинных систем: учет ошибок разных типов//
Восточно-Европейский журнал передовых технологий.
Сер. Системы управления. - Харьков, 2007 - 2/2 (32) -
2007. - С. 63-70.