Ayibazova Aminat Abdullakhovna
Karachai-Cherkess State Thechnological Academy.
E-mail: [email protected].
36, Stavropolskaya Street, Cherkessk, 357100, Russia.
Phone: +7782202387.
The Department of Mathematics; Competitor.
Zargarjan Jury Arturovich
Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”.
E-mail: [email protected].
44, Nekrasovskiy, Taganrog, 347928, Russia.
Phone: +78634371773.
The Department of Automatic Control Systems; Postgraduate Student.
УДК 681.518
E.B.Заргарян ФОРМАЛИЗАЦИЯ ПАРАМЕТРОВ ЗАДАЧ В УСЛОВИЯХ НЕОПРЕДЕЛЕННОСТИ С ПРИМЕНЕНИЕМ ТЕОРИИ РИСКОВ
Рассматриваются элементы теории рисков при задании параметров функционирования объекта и параметров задачи в условиях неопределенности в виде нечетких интер-
. .
Рассмотрены часто применяемые в математических моделях нечеткие случайные величины: биномиально распределенная нечеткая случайная величина, распределение Парето и его связь с теорией рисков.
Нечеткие интервалы; распределение Парето.
E.V. Zargarjan FORMALIZATION PARAMETER PROBLEMS UNDER UNCERTAINTY WITH APPLICATION OF THE RISKS THEORY
We consider the elements of the theory of risk when setting the parameters of the functioning of the object and the parameters of the problem under uncertainty in the form of fuzzy intervals. Lemma of representation of random variables as fuzzy intervals. Reviewed frequently used in mathematical models of fuzzy random variables: binomial distribution, a fuzzy random variable, the distribution of the pair-and his connection with the theory of risk.
Fuzzy intervals; Pareto distribution.
Достижения научно-технического прогресса реализуются на практике в виде создания производственных систем все более возрастающей сложности. Действующие производственные системы являются воплощением замыслов их создателей и естественным образом проходят традиционные стадии жизненного цикла: прикладные научные исследования, разработка технических требований и технического задания (фаза "Исследование"), эскизное, техническое и рабочее проектирование (фаза "Проектирование"), строительство и изготовление оборудования (фаза "Изготовление"), опытная и промышленная эксплуатация (фаза "Функцио-) , ,
( " ").
При этом для проектов любой степени сложности - от минимальных до особо крупных производственных систем - характерен риск не достичь исходных
требований, заложенных в техническом задании, а при реализации - риск нежелательного развития событий, приводящих к непредсказуемым последствиям.
Одним из возможных способов преодоления отмеченного недостатка может
,
производственной системы в целом, а не изолированно по отдельным фазам, пере.
Параметры исследуемой системы зачастую невозможно представить в виде , -лизуемой и зависит от большого количества параметров, поэтому параметры исследуемой производственной системы необходимо представить в виде нечетких .
условиях неопределенности с применением теории рисков.
Классические работы по моделям рисков дают следующую классификацию моделей риска [1]:
♦ модель индиви дуального риска;
♦ модель коллективного риска, (динамическая модель).
При анализе рисков, при построении математических моделей рисков важными проблемами являются:
♦ установление причинных и следственных взаимосвязей разных рисков;
♦ установление и построение модели вероятностных распределений и вероятностных характеристик рисков;
♦
.
Для формализации параметров рассматриваемой производственной системы необходимо использовать понятие нечеткого интервала.
Нечёткий интервал - это неопределенное множество А С Я со средним интервалом, чьи элементы обладают ф ункцией принадлежности цА(Х) = 1. Как и для , , -ной и по крайней мере кусочно-непрерывной.
В различных практических приложениях могут применяться разнообразные классы случайных величин как дискретных, так и (абсолютно) непрерывных случайных величин [2].
В прикладной математике [2,3] часто индивидуальный риск в условиях неопределенности X представляют в виде X = 1( А ) • У , в котором нечеткая случайная величина 1(А) ( А ) -
чение (1,1,0,0) при наступлении нечеткого случайного события А и нечеткое значение (0,0,0,0) при не наступлении А. Нечеткая случайная величина У описывает величину реально предъявленного производственной компании риска. Легко ,
Р(1 = (0,0,0,0)) = Р( X = (0,0,0,0)), Р(I = (1,1,0,0)) = Р(X > (0,0,0,0)), Р(У < у) = Р(X < у IX > (0,0,0,0)).
В частности, если нечеткая случайная величина X дискретного типа, то выполнятся равенства [4,5]:
Р(У = Хк) = Р(Х = Хк/Х > (0,0,0,0)) =------~Р(~Х ак)---------, к = 1,2,...
к к 1 — Р(Х = (0,0,0,0))
при каждом возможном нечетком значении а >0.
Нечеткий риск часто записывают в виде суммы нечеткого случайного числа
V случайных нечетких слагаемых Х1,...., Ху. Весьма типичны в этих моделях требования взаимной независимости и одинаковой функции принадлежности нечетких случайных величин Х1,...., Х]/.
В общем курсе теории вероятностей известен следующий результат о вычислении первых моментов для суммы случайного числа случайных величин.
1. Лемма (тождество А. Вальда) для представления случайных величин в виде нечетких интервалов.
Пусть дано нечеткое случайное число V взаимно независимых и одинаковой
функции принадлежности нечетких случайных величин Х1,....,Хх,, причем
Х = Х1 + - + Хх, и существуют нечеткое математическое ожидание В(Х) И
дисперсия Ъ(Х).
Тогда справедливы нечеткие равенства:
Е(Х) = Е( ~) • Е(Х1),ЩХ) = Е( ~) • ЩХ1) + 5( ~) • (ЩХ1))2; (1) :
Вычислим нечеткую функцию принадлежности случайной величины X с помощью формулы полной вероятности:
Р (X) = Р(Х < X) = £ Р(у = к)Р( Х < XIV = к) =
к =0
(2)
= ^ Р(У = к)Р(Х1 + • • • + Хк IV = к) =£ Р(у = к)Рк (X),
к =0 к =0
где
Бк(Х) = Р(Х1 + — + Хк < х/ V = к ) = Р( Х1 + — + Хк < х ), (3)
поэтому
^ ^
Е(Х) = | хйР (X) = 2 Р(V = к) | хйРк (X) =
—^ к=0 —^
= 2 Р (V = к) Е (X1 + • • •+X к) = 2 ^ (Х 1) Р (V = к) =
к=0 к=0
= Е (X1 )2 кР (V = к) = Е (^1) Е (V).
к =0
Аналогичные нечеткие преобразования позволяют выразить второй нечеткий начальный момент для Х в виде
Ё(Х2 ) = | х2ёР(х ) = X Р(~ = к ) |х2йРк(х ) =
—^ к =0 —го
го
= X Р(х = к)Ё((Х1 + + Хк)2) =
к =0
го Х ~ к х ххх
= X Р(х = к)( XЁ(Х2) + 2 X Ё(Х]Х>) =
к =0 г =1 1< ] <г < к
= X {кЁ( Х1) + к(2^( Ё(Х1))2}Р( V = к) =
к =0
= Ё( Х1 )Ё( V) + (Ё(Х1))2Ё( Х(Х — (1,1,0,0))), таким образом, для дисперсии следует цепочка равенств
Ъ(Х ) = Ё(Х2) — (Е(Х ))2 =
= Ё( X2 )Ё( V) + (Ё(Х1))2Ё{ у( V — (1,1,0,0))) — (Ё(Х )Ё( V ))2 =
= Е(Х2)Е(V ) + Ё2(Х ){Ё(V2 ) — Ё(V )} — Ё2(Х )Ё2(V ) =
= ^( V )ЦХ ) + (Ё(Х ))2 Щ V ).
Следовательно, лемма доказана.
Некоторые часто применяемые в математических моделях нечеткие случайные величины:
а) Биномиально распределенная нечеткая случайная величина Х - это целочисленная нечеткая случайная величина, функция принадлежности которой задается равенством [6]:
Р(Х = к) = скрк((1,1,0,0) — р)п—к,к = 0,1,...п,0 < р < 1
Этот факт часто для краткости записывают в форме Х Е Вп( п,р ) .
Отметим важное для разных приложений свойство нечеткой биномиальной :
Если независимые случайные величины имеют нечеткую биномиальную функцию принадлежности Х Е Вп( N ,р),У Е Вп( К,р ), то функция принадлежности суммы Х + У Е Вп( N + К, р ). Одно из применений этого свойства -возможность в простом выводе следующего свойства, если Х Е Вп( п,р ), то имеет место представление Х = Х1 +... + Хп, в котором все слагаемые - взаимно независимые нечеткие случайные величины, Хг Е Вп(( 1,1,0,0 ),р ), г = 1,2,...,п.
Ё (Х) = Ё (Х,) +... + Ё( Хп) = пЁ (Х,) = пр, щ Х) = й( Х,)+...+й( Хп) = пй( Х,) = пр((1,1, о, 0)—р).
б) Распределение Парето. Нечеткая случайная величина X имеет функцию
принадлежности Парето с параметрами > 0, Л > 0 ,если ее нечеткая функция принадлежности задается равенством [4]
f (*) = f (£Г^ x >°.
Легко вычислить среднее и дисперсию в условиях неопределенности:
го го ~ го
E(X) = \xf(x)dx = \х Л( -Л— )a+1dx = [(7^)“ dx = ^,
0 0 Л Л+х О
го
E(X2) = \x2f(x)dx = X Л( Л f+1dx = 2fx(Tfir dx =
0 0 Л A+x О
(2,2,О,О )12 I2 =__________
D (X)
(а-(2,2,О,О ))(«-(1,1,О,О ))2 {а-(1,1,О,О))2 {а-(2,2ДО))(а-(1,1,О,О)2)
Нечеткое математическое ожидание распределения Парето бесконечно при а > 1, а дисперсия бесконечна при <5 > 2 . Отметим, что поскольку нечеткая экспонента убывает быстрее степенной функции, то функции принадлежности Парето соответствует более частое появление больших рисков. Чтобы сравнить с экспоненциальным распределением, достаточно проверить, что при а = 3
Р(X >10£(X)) = (т=-\^~)& = С~(1ЙЩ0))2 = (0.006,0.006,0,0).
^+а-(1,1,0,0)
Если сравнивать гамма распределение и функцию принадлежности Парето, то иски большой величины с меньшей вероятностью будут появляться при гамма, .
О связи функции принадлежности Парето и гамма-распределения [5].
Хорошим примером этого является широко применяемый пример нечеткой рандомизации показательного распределения по смешивающему гамма-распределению [6]. Более конкретно, пусть
~(~’У) = = $е~х®,
У(Х) = ±х(х) = 1а>х а-1е~1х ,х > (0,0,0,0)
Тогда
го го
у(х) = | у(х,х )]~(х^у = | уе~ уа ■(Ш)’0) УЧу =
0 0
=_______аТ_______= а_ ( 1 )0-+(1,1,0,0)
(х+Х )“+(1Д,0,°) +!' '
Получившееся выражение определяет функцию принадлежности Парето. Та, ,
.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Саами Т.Л. Принятие решений. Метод анализа иерархий. - М.: Радио и связь, 1993.
- 320 с.
2. Gogus O., Boucher T.O. A consistency test for rational weights in multi-criterion decision analysis with fuzzy pair wise comparisons // Fuzzy sets and Systems. - 1977. - Vol. 86.
- P. 129-138.
3. Аверкин AM., Батыршин ИЗ., Блиншун А.Ф., Силаев Б.В., Тарасов Б.Н. Нечет кие множества в моделях управления и искусственного интеллекта. - М.: Наука, 1986. - 312 с.
4. . . -
риев // Известия ЮФУ. Технические науки. - 2009. - № 5 (94). - С. 117-122.
5. . ., . . -
ного сравнения // Известия ЮФУ. Технические науки. - 2010. - № 1 (102). - С. 100-104.
6. . ., . . -роса мнений // Известия ЮФУ. Технические науки. - 2010. - № 1 (102). - С. 104-110.
Статью рекомендовал к опубликованию д.т.н., профессор Я.Е. Ромм.
Заргарян Елена Валерьевна
Технологический институт федерального государственного автономного образовательного учреждения высшего профессионального образования «Южный федеральный университет» в г. Таганроге.
E-mail: [email protected]
347928, . , . , 44.
.: 88634371773.
Кафедра систем автоматического управления; доцент.
Zargarjan Elena Valerevna
Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”.
E-mail: [email protected].
44, Nekrasovskiy, Taganrog, 347928, Russia.
Phone: +78634371773.
The Department of Automatic Control Systems; Associate Professor.
УДК 681.142.34
B.B. Соловьев, В.И. Финаев МОДЕЛЬ РЕГУЛЯТОРА ДЛЯ ЛОКАЛЬНЫХ СИСТЕМ РЕГУЛИРОВАНИЯ
C
Выполнен анализ форм математических моделей. Приведена структура ПД-регуля-тора с нечеткой логикой. Разработана структурная схема локальной системы регулирования. Определены аналитические выражения для треугольных функций принадлежности. Разработана и скорректирована база правил нечеткого ПД-ре^лятора. Разработан алгоритм блока коррекции выходного сигнала регулятора. Выполнено исследование замкнутой системы с нечетким ПД-ре^лятором в среде MatLab Simulink. Выполнена оценка показателей качества управления.
Регулятор с элементами нечеткой логики; коррекция выходного сигнала.
V.V. Soloviev, V.I. Finaev
MODEL OF THE CONTROLLER FOR LOCAL SYSTEMS OF REGULATION WITH APPLICATION OF METHODS OF THE ARTIFICIAL INTELLECT
The analysis of forms of mathematical models is made. The structure of the PD-regulator with fuzzy logic is resulted. The block diagramme of local system of regulation is developed. Analytical expressions for the triangular membership functions are defined. The rule base of a fuzzy PD-controller is developed and corrected. The algorithm of the block of correction for output
1бб