CC BY
56
УДК 007:631.4
В.П. ДИМИТРОВ, Л.В. БОРИСОВА, И.Н. НУРУТДИНОВА
МЕТОДИКА ОЦЕНКИ СОГЛАСОВАННОСТИ МОДЕЛЕЙ НЕЧЕТКИХ ЭКСПЕРТНЫХ ЗНАНИЙ
Рассмотрены некоторые аспекты методики расчета показателей различия и согласованности экспертной информации. Получены аналитические выражения указанных показателей для 13-ти типовых вариантов парного расположения функций принадлежности (ФП). Приведен пример использования методики.
Ключевые слова: лингвистическая переменная, функция принадлежности, показатель различия, показатель согласованности.
Введение. При построении интеллектуальных информационных систем поддержки принятия решений в сфере эксплуатации сложных машин широко используются базы знаний, основанные на нечетких знаниях, т. е. нечеткие продукционные системы [1]. Поэтому одним из актуальных является вопрос представления нечеткой экспертной информации. Блок приобретения и корректировки знаний - один из основных блоков интеллектуальных информационных систем (экспертных систем).
Нечеткая экспертная информация трудноформализуема в рамках традиционных математических формализмов. При отображении лингвистических значений качественных признаков на числовые элементы порядковых шкал информация огрубляется, теряется та ее ценная составляющая, которая характеризует индивидуальный опыт и знания эксперта. Приближенное представление значений функций принадлежности термов семантических пространств может приводить к неадекватности нечетких моделей субъективным суждениям и исходным данным.
Например, для описания признаков предметной области «технологическая регулировка зерноуборочных комбайнов» эксперты могут применять разные множества их лингвистических значений. В одном случае возникают трудности в связи с недостаточностью значений, в другом в связи с их избыточностью. В результате чего следует ожидать увеличения нечеткости и рассогласованности поступающей от экспертов информации. Естественным вопросом при оценивании экспертом проявлений признаков является вопрос: «По каким критериям должен производиться выбор оптимального множества значений лингвистической шкалы, которая применяется при оценивании того или иного признака?».
Описание критериев оптимальности выбора значений лингвистических переменных в [2] содержит следующие требования:
- минимальная неопределенность для экспертов при описании реальных объектов;
- максимальная согласованность экспертной информации.
Важной практической задачей является выбор оптимального множества лингвистической шкалы, используемой для оценивания факторов внешней среды, регулируемых параметров машины и показателей качества работы.
Постановка задачи. Для адекватного представления нечеткой экспертной информации о предметной области необходимо установить оптимальное число термов (m0) лингвистической переменной (ЛП), учитывая ограничения снизу inf m0 и сверху sup m0. Сверху число термов ограничено соображениями точности измерения рассматриваемого параметра. Например, для лингвистической переменной «засоренность хлебной массы» эта точность не может быть высокой, а, следовательно, число термов не может быть большим. Если же выбрать число термов слишком малым, то затруднено выявление и описание взаимодействий данного фактора с показателями качества работы (такой случай приводит к вырождению данной зависимости). При решении рас-
сматриваемой задачи целесообразно проводить оценку согласованности нечетких экспертных знаний.
Рассмотрим некоторые аспекты вопроса о выборе оптимального числа термов лингвистической переменной на примере ЛП «засоренность хлебной массы». В данном случае пять экспертов дали оценки функций принадлежности для трех термов («низкая», «средняя», «высокая»), для четырех термов («низкая», «средняя», «высокая», «очень высокая») и пяти термов («очень низкая», «низкая», «средняя», «высокая», «очень высокая»). Вид функций принадлежности и значения параметров трех моделей представлены в табл.1-4.
Номер эксперта будем обозначать индексом I = 1, 2, ..., k, а номер терма - индексом I = 1, 2, ..., т. Функцию принадлежности, которую задал ий эксперт для 1-го терма, обозначим
Рп (х).
Таблица 1
Вид функций принадлежности
Вид терма Выражение для ФП
Левый терм ¡и (х, а, Ь) = - 1, если х < а Ь - х , если а < х < Ь Ь - а 0, если х > Ь
Центральные термы и (х, а, Ь, с, d) = < 0, если х < а х - а , если а < х < с с - а 1, если с < х < d Ь - х , , , если d < х < Ь Ь - d 0, если х > Ь
Правый терм ии (х, а, Ь) = - 0, если х < а х - а , если х < Ь Ь - а 1, если х > Ь
Область определения для х от 0 до 1 (нормированные значения). Область определения для а, Ь, с, d от 0 до 1.
Таблица 2
Значения коэффициентов а, Ь, с, d функций принадлежности для различных экспертов
5 термов Значения параметров ФП для различных экспертов
Вид терма 1-й эксперт 2-й эксперт 3-й эксперт 4-й эксперт 5-й эксперт
1 2 3 4 5 6
Левый а = 0,04 Ь = 0,10 а = 0,05 Ь = 0,12 а = 0,05 Ь = 0,14 а = 0,03 Ь = 0,10 а = 0,08 Ь = 0,16
Центральный а = 0,04 Ь = 0,20 с = 0,10 d = 0,14 а = 0,05 Ь = 0,20 с = 0,12 d = 0,15 а = 0,05 Ь = 0,22 с = 0,14 d = 0,16 а = 0,03 Ь = 0,20 с = 0,10 d = 0,14 а = 0,08 Ь = 0,26 с = 0,16 d = 0,20
Окончание табл.2
1 2 3 4 5 6
а = 0,14 а = 0,15 а = 0,16 а = 0,14 а = 0,20
Ь = 0,30 Ь = 0,32 Ь = 0,32 Ь = 0,30 Ь = 0,35
Центральный с = 0,20 с = 0,20 с = 0,22 с = 0,20 с = 0,26
d = 0,24 d = 0,26 d = 0,26 d = 0,25 d = 0,30
а = 0,24 а = 26 а = 0,26 а = 0,25 а = 0,30
Ь = 0,44 Ь = 0,44 Ь = 0,46 Ь = 0,42 Ь = 0,44
Центральный с = 0,30 с = 0,32 с = 0,32 с = 0,30 с = 0,35
d = 0,36 d = 0,38 d = 0,36 d = 0,36 d = 0,38
а = 0,36 а = 0,38 а = 0,36 а = 0,36 а = 0,38
Правый Ь = 0,44 Ь = 0,44 Ь = 0,46 Ь = 0,42 Ь = 0,44
Таблица 3
Значения коэффициентов а, Ь, с, d функций принадлежности для различных экспертов
4 терма Значения параметров ФП для различных экспертов
Вид терма 1-й эксперт 2-й эксперт 3-й эксперт 4-й эксперт 5-й эксперт
Левый а = 0,04 а = 0,05 а = 0,05 а = 0,03 а = 0,08
Ь = 0,12 Ь = 0,12 Ь = 0,14 Ь = 0,10 Ь = 0,16
а = 0,04 а = 0,05 а = 0,05 а = 0,03 а = 0,08
Центральный Ь = 0,22 с = 0,12 Ь = 0,20 с = 0,12 Ь = 0,24 с = 0,14 Ь = 0,22 с = 0,10 Ь = 0,26 с = 0,16
d = 0,16 d = 0,15 d = 0,16 d = 0,16 d = 0,20
а = 0,16 а = 0,15 а = 0,16 а = 0,16 а = 0,20
Центральный Ь = 0,44 с = 0,22 Ь = 0,40 с = 0,20 Ь = 0,36 с = 0,24 Ь = 0,36 с = 0,22 Ь = 0,44 с = 0,26
d = 0,34 d = 0,34 d = 0,28 d = 0,28 d = 0,32
Правый а = 0,34 Ь = 0,44 а = 0,34 Ь = 0,40 а = 0,28 Ь = 0,36 а = 0,28 Ь = 0,36 а = 0,32 Ь = 0,44
Таблица 4
Значения коэффициентов а, Ь, с, d функций принадлежности для различных экспертов
3 терма Значения параметров ФП для различных экспертов
Вид терма 1-й эксперт 2-й эксперт 3-й эксперт 4-й эксперт 5-й эксперт
Левый а = 0,04 а = 0,07 а = 0,04 а = 0,03 а = 0,08
Ь = 0,16 Ь = 0,20 Ь = 0,14 Ь = 0,15 Ь = 0,19
а = 0,04 а = 0,07 а = 0,04 а = 0,03 а = 0,08
Центральный Ь = 0,36 с = 0,16 Ь = 0,38 с = 0,20 Ь = 0,30 с = 0,14 Ь = 0,38 с = 0,15 Ь = 0,40 с = 0,19
d = 0,24 d = 0,24 d = 0,18 d = 0,24 d = 0,26
Правый а = 0,24 Ь = 0,36 а = 0,24 Ь = 0,38 а = 0,18 Ь = 0,30 а = 0,24 Ь = 0,38 а = 0,26 Ь = 0,40
При анализе согласованности нечеткой экспертной информации вначале вычисляются аддитивный и мультипликативный показатели общей согласованности, а потом по их величинам формулируется суждение о согласованности моделей экспертного оценивания. Затем строится матрица парной согласованности моделей X- и Х}- экспертов.
Общая согласованность множества моделей экспертного оценивания признака определяется аддитивным k и мультипликативным k показателями [3]:
Уг=1,2,..Д
k = —£т--------------; (1)
ш і=1 | тах ца (x)dx
Уї=1,2,...^
k = ш
1 тіп цп (х)аХ
УЇ=1,2,..., k
I
П Т------------71-Г ■ (2>
і=1 I тах ¡ла (х)ах
V;=1,2,..., k
0
Очевидно, что 0 < k < 1, 0 < k < 1.
Если все модели X, Х2, ..., Xk совпадают (все эксперты одинаково определяют функции принадлежности у всех термов), то k = k = 1. Если нет пересечений у функций принадлежности всех термов, то k = 0 . Если нет пересечений у функций принадлежности хотя бы одного терма, то ~ = 0.
Показатель различия между моделями двух экспертов /-го иу-го в рамках 1-го терма определяется как линейное расстояние (Хемминга) между нечеткими множествами с функциями принадлежности ¡лп (х) и Мц (х) [1].
1
d = Ц Ми(х) - Ма (х)| ^. (3)
0
Показатель согласованности между этими же моделями определяется величиной ^ [1]:
1
I т1п [ м~ (х X м~ (х)] ^ k у = т . (4)
Iтах [м~ (хX м~(х)]^
0
Затем строится матрица парной согласованности К моделей Xi и Ху экспертов. Очевидно, что на главной диагонали матрицы стоят единицы, и матрица симметрична.
На основе матрицы парной согласованности моделей для всех термов находим матрицу согласованности моделей Xi иХу по всем термам. Ее элементы определяются формулой [2]:
~ 1 т ~
~ = т£~ ' (5)
т I=1
где т - число термов.
Типичные случаи показателей различия. Рассмотрим возможные случаи определения выражений (3) и (4) для различных вариантов взаимного расположения функций принадлежности
Мя (х) и Мц (х).
1. Крайний левый терм
Носитель обеих функций принадлежности Мц (х) и Ма (х) включает точку х = 0. На тех интервалах, где Ми (х) = Мц (х) =1 или Мц (х) = Мц (х) =0, подинтегральная функция равна 0. Из интервала [0, 1] интегрирования (3) остаются только участки:
а) где одна из функций Ми (х) или Мц (х) тождественно равна 1, а другая является убывающей линейной функцией;
б) где обе функции Ми (х) и Мц (х) являются убывающими линейными функциями;
в) где одна из функций Мц (х) или Мц (х) тождественно равна 0, а другая является убывающей линейной функцией;
г) где одна из функций Мц (х) или Мц (х) равна 1, а другая равна 0.
2. Крайний правый терм
Носитель обеих функций принадлежности Мц (х) и Ма (х) включает точку х = 1. На тех интервалах, где Мц (х) = Ма (х) =1 или Мц (х) = Ма (х) =0, подинтегральная функция в формуле (3) равна 0. Из интервала [0, 1] интегрирования (3) остаются только участки:
а) где одна из функций Мц (х) или Ма (х) тождественно равна 1, а другая является линейно возрастающей функцией;
б) где обе функции Мц (х) и Ма (х) являются линейно возрастающими функциями;
в) где одна из функций Мц (х) или Ма (х) тождественно равна 0, а другая является линейно возрастающей функцией;
г) где одна из функций Мц (х) или Ма (х) равна 1, а другая равна 0.
3. Центральный терм
На тех интервалах, где Мц (х) = Ма (х) =1 или Мц (х) = Ма (х) =0, подинтегральная функция в формуле (3) равна 0. Из интервала [0, 1] интегрирования остаются следующие участки:
а) одна из функций Мц (х) или Ма (х) тождественно равна 1, а другая является линейно убывающей или линейно возрастающей;
б) обе функции Мц (х) и Ма (х) являются линейно возрастающими или линейно убывающими;
в) одна из функций Мц (х) или Ма (х) тождественно равна 0, а другая является линейно возрастающей или линейно убывающей;
г) одна из функций Мц (х) и Ма (х) является линейно возрастающей, а другая линейно убывающей;
д) одна из функций Мц (х) или Ма (х) тождественно равна 0, а другая равна 1.
Подразумевая в дальнейшем использование для расчетов показателей согласованности конкретной ЛП, приведем результаты для типичных случаев.
Толерантное нечеткое число рассматриваемого признака символически запишем в виде четвёрки чисел:
А = (а х,а 2,аь,ак) (6)
или
(х)= (ах, а2,аь,ак),
где а1, а2 - границы толерантности; аь,ак - левый и правый коэффициенты нечёткости. Аналогично
В = М 2,ьь,ьк)
или
МВ (х)= (Ь1Ь 2,ЬЪ,ЬК ).
Типичные случаи показателей различия. В дальнейшем используются следующие переобозначения параметров применяемых типовых функций (см. табл.1) соответственно:
- для левого терма - а есть а2 (для эксперта «1») или Ь2 (для эксперта «2»); Ь есть ая
или Ья;
- для правого терма - а есть а |_ (для эксперта «1») или bL (для эксперта «2»); Ь есть а1
или Ь1,
- для центрального терма - а есть а I (для эксперта «1») или bL (для эксперта «2»); Ь есть ац или ^; с есть а! (для эксперта «1») или Ь1 (для эксперта «2»); d есть а2 (для эксперта «1») или Ь2 (для эксперта «2»).
1. Крайний левый терм
Пусть переменные А и В описывают крайний левый терм, задаваемый разными экспертами. Тогда а1 = Ь= 0 и ах = ЬL = 0, (рис.1).
«ь=Ьі
аь=Ьь
а)
б)
Рис.1. Графическое изображение варианта
Рассмотрим случай, когда а2 <Ь2; aR <bR (см. рис.1,а), тогда:
(Ь2 + ЬЯ ) - (а2 + аЯ ).
d =
2
а 2 + а ^
к — __________________
У Ьп + ь
(7)
(8)
'2 + ЬЯ ■
Для случая, когда а2 >Ь2, aR >Ьц используются формулы (7) и (8) с учетом переобозначения а и Ь.
Рассмотрим случай, когда а2 <Ь2, aR > Ьц (см. рис.1,б), тогда:
d — (Ь2 ~ а2)2 + (аК ~ ЬК )2 , 2(Ь2 ~ а2 + аЯ ~ ЬК )
2
(9)
к —
(аД + а2 )~ Ь ЬК ) Ь
Ь2 ~а2 + аК ~ЬК
(аК + а2 )
+
(Ь2 ~а2)2
Ь2 ~а2 + аК ~ЬК
(10)
2. Крайний правый терм
Если ах <Ьх; а1 <Ьг (рис.2,а), тогда:
d =
Ь1 ~ а1 + Ьь ~ а£ 2
(11)
k =
2 - (Ь1 + Ьх )
2 - (а 1 + ах )
(12)
а) б)
Рис.2. Иллюстрация для крайнего правого терма
Для случая, если ах >Ьь а! >Ь!, путем переобозначения получаем формулы (11)-(12). Если а. > Ь; Эх <Ьх (рис.2,б), тогда:
d =
(Ь1 - а1) + (аь - Ьь ) _ 2(Ь1- а1+аь- Ьь) ;
2 - (Ь1 + Ьх ) -
(ах - ЬхУ
к =
(Ь1 - а 1 + ах - Ьх )
2 - (Ь1 + Ьх ) +
(Ь1 - а 1)2
(13)
(14)
(Ь1 - а 1 + ах - Ьх )
Если ах <Ьх, а! >Ь! используются формулы (13)-(14) с учетом переобозначения.
3. Центральный терм
Вариант 1 (рис.3). Для случая ах <Ьх, а! <Ь!, а2 >Ь2, aR > ЬR выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
(Ь1 + ЬХ ) + (а2 + аК ) - (а1 + аЬ ) - (Ь2 + ЬЯ )
а
2
к =
(Ь 2 + ЬR ) - (Ь1 + ЬЬ ) ( а 2 + aR ) - ( а 1 + аЬ )
(15)
(16)
Рис.3. Иллюстрация к варианту 1
Вариант 2 (рис.4). Для случая ах >Ьх, а! <Ь!, а2 <Ь2, ар> <Ьр> выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
а = (аЬ - Ьь )2 + (Ь1 - а1)2 + (Ь2 + ЬR ) - (а2 + aR ) _
2(Ь1- а1+аь - Ьь) 2 ;
(Ь - а О2
(17)
(ак ~ аь ) + (а 2 ~ а 1 ) ~
к —
(Ьі ~ а і + аь ~ Ьь )
(ЬК ~ Ьь ) + (Ь 2 ~ Ь1 ) +
(Ь1 ~ а 1)2
(18)
(Ьх - а! + ах - Ьх )
Вариант 3 (рис.5). Для случая ах >Ьх, а! >Ь!, а2 >Ь2, aR >ЬR выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
(а! + ах + а2 + aR ) - (Ь1 + Ьх + Ь2 + ЬR ) .
d —
2
к —
(Ь 2 + ЬД ) ~ (а 1 + аЬ ) (а 2 + ай ) ~ (Ь1 + ЬЬ )
(19)
(20)
Рис.4. Иллюстрация для варианта 2
Рис.5. Иллюстрация для варианта 3
Вариант 4 (рис.6). Для случая ах <Ьх, а! <Ь!, а2 <Ь2, aR <ЬR (Ь! > а2) выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
d —
(аЯ ~ аь ) + (а2 ~ а1) + (ЬЯ ~ Ьь ) + (Ь2 ~ Ь1)
(ая ~ Ьь )2
2
(Ь1 ~ а2) + (ай ~ Ь1 )
( ая Ьь)
к —
(Ь ~ а2 + ая ~ Ьь )
(аЯ ~ аь ) + (а 2 ~ а 1) + (ЬЯ ~ Ьь ) + (Ь 2 ~ Ь1) +
(аЕ ~Ьь)2 (Ь ~ а2 + ак ~ Ьь )
(21)
(22)
Вариант 5 (рис.7). Для случая ах <Ьх, а! >Ь!, а2 <Ь2, aR <Ьц выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
^ _ (Ьь ~ аь )2 + (а 1 ~ Ь1 )2 + (Ь 2 + ЬЯ ) ~ (а 2 + аЯ ) ■
2(а 1 ~ Ь1 + Ьь ~ аь )
2
(23)
(ая - ЬЯ + а 2 - Ь1 ) -
к =
(а 1 - Ь1 + Ьх - ах )
(ЬЯ - аЬ + Ь 2 - а 1 ) +
(а 1 - Ь1 + ЬЬ - ах )
(24)
Рис.6. Иллюстрация для варианта 4
Рис.7. Иллюстрация для варианта 5
Вариант 6 (рис.8). Для случая ах <Ьх, а! <Ь!, а2 >Ь2, aR <ЬR выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
а = (Ь1 + Ьь ) - (а1 + аь ) + (а2 - Ь2) + (ЬЯ - ак )2
2
(ЬК - Ьь ) + (Ь 2 - Ь1) -
2(а2 - Ь2 + ЬЯ - аД )
(Ьк - ал )'
к =
(а 2 - Ь 2 + ЬК - аЛ )
(ак - ах ) + (а 2 - а 1) +
(Ьк - ая)2
(а 2 - Ь 2 + ЬЛ - аК )
(25)
(26)
Рис.8. Иллюстрация для варианта 6 среднего терма
Рис.9. Иллюстрация для варианта 7 среднего терма
Вариант 7 (рис.9). Для случая ах >Ьх, а! <Ь!, а2 >Ь2, aR < ЬR выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
а = (аЬ - Ьь ) + (Ь1 - а1) + (а2 - Ь2) + (ЬК - аК )2
2(Ь1 - а1 + аь - Ьь ) 2(а2
(Ьх - ах)2
Ь2 + Ьй
(ЬR - ЬХ ) + (Ь2 - Ь1) -
~ая )
(ЬR - аR)
к =
(Ь1 - а1 + аь - Ьх ) (а 2 - Ь2 + ЬR - аR )
(aR - ) + (а2 - ах) +
+
(ЬR aR)
(27)
(28)
(Ь1 - а1)2
(Ь1 - + ах - Ьх) (а2 - Ь2 + ЬR - aR )
Вариант 8 (рис.10). Для случая ах >Ьх, а! <Ь!, а2 >Ь2, aR >ЬR выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
213
d■
к =
(aL ~ bL )2 + (Ь1 ~ а1)2 , (а2 + а2) ~ (Ь2 + bR )
2
(аь ~ Ьь )2
(Ьі ~ а 1 + аь ~ ьь )
2(Ь1 ~ аі + аL ~ bL )
(bR ~ ЬЬ ) + (Ь 2 ~ Ь1 ) ~
(aR ~ аі ) + (а 2 ~ а 1) +
(аь ~ Ьь )2
(Ьі ~ а і + ^ ~ Ьь )
(29)
(30)
Рис.10. Иллюстрация для варианта 8 среднего терма
Рис.11. Иллюстрация для варианта 9 среднего терма
Вариант 9 (рис.11). Для случая а>ЪЬ а1 >Ъ1, а2 >Ъ2, aR <ЪR выражения для расчета значения индекса нечеткости и показателя согласованности имеют вид:
d = (а1 + aL ) ~ (Ь1 + bL ) + (а2 ~ Ь2) + (bR ~ aR )2
2
2(а2
Ь2 + bR
aR )
(aR ~ аь ) + (а 2 ~ а 1 ) ~
( а 2 ~ Ь 2 ) 2
к=
(а 2 ~ Ь 2 + bR ~ aR )
(bR ~ ЬЬ ) + (Ь 2 ~ Ь1 ) +
( а 2 ~ Ь 2 ) 2
(31)
(32)
(а 2 - Ъ 2 + ^ ^ )
Для рассматриваемых функций принадлежности (см. табл.2-4) результаты расчетов аддитивных и мультипликативных показателей по формулам (1) и (2) представлены в табл.5.
Таблица 5
Результаты расчета показателей к и к
Модель к к
3-термовая 0,655 0,628
4-термовая 0,542 0,503
5-термовая 0,509 0,460
Анализ полученных показателей общей согласованности экспертной информации показывает, что наиболее согласованной является 3-термовая модель, менее согласована 4-термовая и еще менее - 5-термовая. При этом различие в степени согласованности 4- и 5-термовых моделей незначительно. Одна из причин этого заключается в том, что увеличение числа термов влечет за собой и усиление рассогласованности экспертной информации.
Сопоставим полученные результаты с результатами определения оптимального числа термов ЛП методом минимизации средневзвешенного квадратического отклонения Fm индивидуальных параметров, задаваемых экспертами, от усредненных значений этих параметров. Результаты расчетов Fm для лингвистической переменной «засоренность хлебной массы» представлены в [3], и оптимальное число термов согласно этим расчетам равно трем. Наблюдаемая корреляция
результатов, полученных методом минимизации и путем вычисления показателей общей согласованности, подтверждает, что такой выбор числа термов ЛП позволяет получить адекватное представление нечеткой экспертной информации о предметной области.
Однако, следует учитывать, что выбор оптимального числа термов для каждой ЛП индивидуален. Например, выбор 3-термовой модели относится только к ЛП «засоренность хлебной массы». В других случаях ограничение 3-термовой моделью может затруднить установление взаимосвязей данного фактора с регулируемыми параметрами машины и показателями качества ее работы.
После выбора в качестве оптимальной 3-термовой модели выполним расчеты показателей согласованности экспертных моделей для всех термов ФП. Результаты вычислений показателей
различия d(A,B) и показателей парной согласованности кАВ представим в виде соответствующих матриц, которые обозначим Dи К
Для крайних левых термов D и К имеют вид:
О =
( 0 0,035 0,01 0,01 0,035^ (0,09"
0,035 0 0,045 0,045 0,005 0,13
0,01 0,045 0 0,005 0,045 , 0,105
0,01 0,045 0,05 0 0,045 0,15
ч 0,035 0,005 0,045 0,045 0 V 1 0,13 )
К =
Для центральных термов D и К имеют вид:
D =
С 0
0,05
0,07
0,02
0,07
0,05
0
0,11
0,05
0,03
0,07
0,11
0
0,08
0,14
0,02
0,05
0,08
0
0,07
0,07 ^ 0,03 0,14 0,07 0
( 0,21^ (1
0,24 0,79
0,40 , К = 0,67
0,22 0,91
10,3Ъ 10,72
( 1 0,74 0,9 0,9
0,74 1 0,67 0,67
0,9 0,67 1 0,95
0,9 0,67 0,95 1
ч 0,74 0,96 0,67 0,67
0,79 0,67 0,91 0,72 ^
1 0,48 0,80 0,87
0,48 1 0,69 0,44
0,80 0,69 1 0,73
0,87 0,44 0,73 1 )
0,74 ^ (4,28л
0,96 4,04
0,67 4,20
0,67 4,19
1 ) 14,04 )
( 4,09 ^ 3,94 3,28 4,13 3,76
О =
Для крайних правых термов D и К имеют вид:
(4,852^
4.866
, К = 0,92 0,91 1 0,91 0,88 . 4,62
4.866 ч 4,78
Показатель общей согласованности моделей (5) попарно по всем термам имеет вид:
( 0 0,01 0,06 0,01 0,03^ (0,11л
0,01 0 0,07 0 0,02 0,1
0,06 0,07 0 0,07 0,09 , 0,29
0,01 0 0,07 0 0,02 0,1
V 0,03 0,02 0,09 0,02 0 ) 10,16 )
1
0,986
0,92
0,986
0,96
0,986
1
0,91
1
0,97
0,92
0,91
1
0,91
0,88
0,986
1
0,91
1
0,97
0,96 ^ 0,97 0,88 0,97 1
К =
1
0,84
0,83
0,93
0,81
0,84
1
0,69
0,82
0,94
0,87
0,69
1
0,85
0,66
0,93
0,82
0,85
1
0,79
0,81^ 0,94 0,66 0,79 1
( 4,45 ^ 4,29 4,03 4,39 4,2
В столбцах справа приведены соответствующие суммы элементов строк матриц.
Анализируя суммы элементов матриц О или К по строкам или столбцам можно выделить экспертные модели наиболее попарно согласованные с остальными (наименьшая построчная сумма элементов матрицы О или наибольшая матрицы К ) и наименее согласованные (наибольшая построчная сумма элементов матрицы О или наименьшая матрицы К ). Однако отметим, что значимость этого различия следует проверять в рамках статистического подхода, поскольку при построении обобщенной (усредненной) ФП исключение одной или более экспертных моделей влечет за собой потерю возможно ценной экспертной информации.
Выводы. В результате анализа различных вариантов расположения функций принадлежности выявлены типовые случаи, для которых получены выражения (7)-(32) для расчета коэффициентов нечеткости и согласованности. Полученные результаты будут использованы при разработке алгоритмов модификации в подсистеме ввода нечетких экспертных знаний интеллектуальной информационной системы для технологической регулировки машин.
Библиографический список
1. Борисова Л.В. К вопросу построения нечеткой экспертной системы продукционного типа для технологической регулировки машин / Л.В. Борисова, В.П. Димитров, А.К. Тугенгольд // Вестник ДГТУ. - 2008. -Т.8. - №3(38). - С.278-287.
2. Димитров В.П. Оценка параметров лингвистических переменных факторов внешней среды / В.П. Димитров, Л.В. Борисова // Искусственный интеллект в XXI веке. Решения в условиях неопределенности: сб. ст. V междунар. науч.-техн. конф.- Пенза, 2007.- С.30-32.
3. Обработка нечеткой информации в системах принятия решений/ А.Н. Борисов, А.В. Алексеев, Г.В. Меркурьев и др. - М.: Радио и связь, 1989. - 394 с.
Материал поступил в редакцию 26.02.2010.
V.P. DIMITROV, L.V. BORISSOVA, I.N. NURUTDINOVA
METHODS FOR ESTIMATING COORDINATION OF FUZZY EXPERT KNOWLEDGE MODELS
Some aspects of the methods for calculating indices of expert information difference and coordination are considered. Analytical expressions of the mentioned indices for 13 typical variants of belonging functions' pair arrangement have been obtained. The methods applying example is given.
Key words: linguistic variable, membership function, difference rate, cooperation index.
ДИМИТРОВ Валерий Петрович (р. 1953), декан факультета «Приборостроение и техническое регулирование», заведующий кафедрой «Управление качеством», доктор технических наук (2002), профессор (2003). Окончил РИСХМ (1975).
Сфера научных интересов: системы информационного обеспечения жизненного цикла продукции, экспертные системы, техническое обслуживание машин.
Имеет 349 публикаций (в том числе 32 учебных пособий и монографий).
БОРИСОВА Людмила Викторовна, заведующая кафедрой «Экономика и менеджмент машиностроения» Института энергетики и машиностроения (ИЭиМ) ДГТУ, доктор технических наук (2008), профессор (2010). Окончила РГАСХМ (1991).
Сфера научных интересов: статистика, системы информационного обеспечения жизненного цикла продукции, менеджмент качества.
Имеет 219 публикаций (в том числе 14 учебных пособий и монографий). borisova lv09@ma il.ru
НУРУТДИНОВА Инна Николаевна, доцент кафедры «Математика» Института энергетики и машиностроения (ИЭиМ) ДГТУ, кандидат физико-математических наук (1991). Окончила РГУ (1982).
Сфера научных интересов: математическое моделирование процессов и систем.
Имеет 70 публикаций.
