CC BY
14
Доклады БГУИР
2009 № 4 (42)
УДК 62-506.29:519.21
АНАЛИТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ГРАНИЧНЫХ РЕЖИМОВ РАБОТЫ СТОХАСТИЧЕСКОЙ СИСТЕМЫ
А.А. ЛОБАТЫЙ, Ж.М. САИД
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 3 марта 2009
На основе теории марковских процессов случайной структуры решается задача вероятностного анализа граничных режимов работы стохастической динамической системы. Аналитически получены выражения для вычисления интенсивностей поглощения реализаций случайного процесса. Это позволяет по известным статистическим характеристикам фазовых координат оценивать вероятность безотказной работы системы управления.
Ключевые слова: стохастическая система, случайная структура, плотность вероятности, срыв управления, функция поглощения.
Введение
Существует широкий класс автоматических систем, у которых переход из одного состояния в другое происходит при достижении одной или несколькими фазовыми координатами определенных границ. Задачи вероятностного анализа таких систем рассматриваются в ряде работ на основе различных теоретических подходов при применении двух классов моделей: имитационных и аналитических. Сложность имитационного моделирования состоит в необходимости обеспечения адекватности математической модели реальному процессу. Иногда необходимо моделировать процессы, происходящие на молекулярном уровне, это требует больших экспериментальных исследований на реальных образцах техники. Предсказательные возможности имитационного моделирования значительно меньше, чем у аналитических моделей. Для проведения предварительного проектирования и оптимизации системы, исследования ее характеристик используются аналитические модели, отражающие основные, наиболее важные свойства системы, позволяющие решать широкий круг задач анализа и синтеза. Аналитическая модель дает решение поставленной задачи в законченной форме в виде функциональной зависимости исследуемых характеристик от параметров модели и обеспечивает глубину анализа.
При аналитическом моделировании граничных режимов работы системы управления целесообразно использовать теорию динамических систем случайной структуры [1], в основе которой лежит изучение разрывных марковских процессов с поглощением и восстановлением реализаций. Рассмотрим распространенный класс задач анализа систем случайной структуры — задачи анализа стохастических динамических систем с учетом срыва управления, когда при достижении фазовыми координатами (параметрами) системы определенных границ она переходит в принципиально иное состояние, например неисправное или неработоспособное.
Теоретический анализ граничных режимов
Достаточно общая математическая модель стохастической динамической системы может быть представлена векторно-матричным уравнением в форме Ланжевена вида
X(t) = D(X, t) + F(X, Ш), X(t0 ) = X0,
где Х0 — случайный вектор начального состояния; X ) — «-мерный вектор фазовых координат системы; 0( X, t), Е (X, t) — в общем случае произвольные нелинейная векторная и матричная функции; £,(7) — вектор белого шума с нулевым математическим ожиданием и матрицей интенсивности G(1;) .
Пусть X^) может находиться в различных состояниях (структурах). При этом 1-я структура характеризует основное (рабочее) состояние системы. Выход процесса X(1;) за границы некоторой известной области О , ограниченной гиперповерхностью Г 1г X, I , означает
переход системы в одну из г-х (г = 1,л ) структур (состояний). Порядок смены структур (чередование индексов I, г, к, ...) случаен, представляет собой дискретный процесс ) и подчинен определенным статистическим закономерностям. Нахождение процесса в 1-й структуре обозначим X (I, t) .
Наиболее полной характеристикой процесса 1-й структуры X(I, t) является плотность распределения вероятности (ПРВ) /(X, I, t) , которая при отсутствии восстановления реализаций процесса X(I, t) удовлетворяет обобщенному уравнению Фоккера-Планка-Колмогорова вида [1]
/(X,/,0 = Утхп(Х,/,0 - ЦХ,/,0 . (1)
Здесь VI =
— векторный оператор дифференцирования; п(Х,/,1)
д д д дх, ' дх, дх,
1 а }
— вектор плотности потока вероятности, равный
1 г
п(х,/, о = А(х, /, 0/(х, 1,0-- V* в(х, и 0/(х, /, о ,
где А( X, 1,{) — вектор сноса с составляющими, характеризующими средние значения локальных скоростей компонент непрерывного векторного марковского процесса; В( X, I, t) — диффузионная матрица, каждый элемент которой характеризует скорость изменения условного момента связи компонент векторного марковского процесса. При математическом описании локальных характеристик А(X, 1,() и В(X, I, t) будем рассматривать стохастические интегралы в смысле Ито как более простые для вычислений [2].
Уравнение Фоккера-Планка-Колмогорова является уравнением в частных производных параболического типа. Нахождение решения этого уравнения представляет собой значительную трудность, особенно для многомерных систем. Однако в большинстве случаев этого решения не требуется, так как хорошие результаты позволяет получить гауссова двухмоментная аппроксимация ПРВ /(X, I, t) .
Основная проблема при аналитическом исследовании систем случайной структуры заключается в определении функций поглощения /, I). В зависимости от условий поглощения реализаций процесса X(1;) различают процессы с сосредоточенными и распределенными переходами системы из одного состояния в другое.
В процессах с сосредоточенными переходами смена состояний связана с достижением реализациями X ^) границ некоторой области в фазовом пространстве. Эти границы перехода
процесса из области I в область г представляются детерминированными функциями у 1г (I). Функция поглощения /, I) для многомерного процесса при этом имеет вид:
£
= Х К л Х,1,1 5 Г1г Х,1 -у1г , (2)
г=1
где 5(...) — дельта-функция; п°1г — орт внешней нормали к гиперповерхности I' 1г X,! ; пи. 71 X,/, / — скалярное произведение векторов п/г и п X, /, / .
В процессах с распределенными переходами многомерная функция поглощения для 1-й структуры имеет вид:
£
3(Х,1,1) = у,(Х,1)/(Х,1,1) = . (3)
г=1
Функция гг(Х, — интенсивность поглощения реализаций 1-й структуры. Функции \'1г (X, I) в выражении (3) называются интенсивностями смены структур (переходов из области I в г-ые области). Они могут быть отличными от нуля во всей или части области существования процесса. Если процессы ) и X ) независимы, то интенсивности переходов не зависят от
координатXи становятся функциями только индексов /, г и времени / \> 1г(Х,/) = \>/г(I).
Интеграл от ПРВ / (X, I, ^) по всей бесконечной открытой области без учета границ и стоков вероятности представляет собой вероятность существования реализаций 1-й структуры в каждый текущий момент времени
3(0 =)/(*,/,/)£ЙГ.
Проинтегрировав уравнение (1) по Хв бесконечных пределах, получим дифференциальное уравнение для Р1 (^) .
Если функция поглощения &(Х,/,/) рассматривается в выражении (2), то вследствие свойств дельта-функции задача вероятностного анализа сводится к отысканию значения плотности потока вероятности на поглощающей границе я(угг,0 . В работе [3] проведен подробный
анализ такой задачи. Однако точного решения, особенно для многомерных систем, получить не удается. При этом приходится делать ряд допущений и громоздких преобразований.
Простое решение получается при рассмотрении функции поглощения в уравнении (3) при Э(Х,I) — V, (/)/(Х,/,/). В этом случае при интегрировании уравнения (1) по X получаем простое для решения дифференциальное уравнение:
Таким образом, при решении задач анализа и особенно синтеза удобнее рассматривать функцию поглощения 9{Х, I) в виде (3). Однако основная сложность состоит в приведении
функции V,,. (X, 0 к виду \>1г (0 . В общем случае вопрос аналитического определения интен-сивностей \>1г (0 авторами теории динамических систем случайной структуры подробно не рассматривается, предполагая, что можно определять экспериментальным путем.
Методика определения интенсивности поглощения реализаций
Воспользуемся для определения \>1г (/) эргодическим свойством случайного процесса X (^), считая, что математическое ожидание (корреляционную функцию) случайного процесса можно определить по достаточно длинной реализации X(^). На практике большинство рассматриваемых процессов, происходящих в технических системах, являются стационарными и
таким свойством обладают. Кроме того, гиперповерхность I' 1г X,! сведем к одной поглощающей границе для конкретной фазовой координаты, что также достаточно характерно для реальных условий.
Рассмотрим реализацию компоненты ) процесса X(1) . Пусть переход X(1) из состояния I в состояние (структуру) г проявляется в превышении реализацией х(1) некоторого допустимого уровня Д . Это превышение называется выбросом случайного процесса х{!) [4].
Среднее число положительных выбросов Х(1) за уровень А на интервале [0, 7] случайно и определяется следующим выражением
N.
1
(Т)= ¡х(0 8[х(0-Д] \[х(1)]Ж .
(4)
Подынтегральная функция в (4) вследствие свойств дельта-функции 8[л'(/) — А] и единичной ступенчатой функции 1[л*(7)] равна нулю всюду кроме тех точек, где случайный процесс х(1) пересекает уровень А. В точках пересечения процессом х(1) уровня А интеграл скачком возрастает на единицу. Следовательно, интеграл (4) равен числу положительных пересечений случайным процессом х(/) уровня А на интервале [0,Т].
Математическое ожидание случайной величины МА(Т) вычисляется по формуле:
М N.
«и ои
(5)
где /{х, х, — двумерная плотность вероятности случайного процесса Х(1) и его производной х(0 .
Подставляя формулу (4) в (5) и пользуясь известным правилом интегрирования произведения произвольной функции на дельта-функцию, проинтегрируем (5) по переменной х :
1 сс
м ЛГД(Г) =
(6)
Для случайного процесса х(1), стационарного в узком смысле, для которого любые сечения по временной оси ПРВ /(х, Л", /) зависят лишь от разности аргументов ,..., 1п
/(л",Л",/) = /(\*,Л",/( — / ,) (/, / = \,п), внутренний интеграл в выражении (6) не зависит от времени и выражение (6) можно записать в следующем виде:
М N.
(Т) =г|х(0/(А,х>й.
(7)
Разделив среднее число выбросов на интервале [0, Т] (7) на длину интервала, получим интенсивность числа выбросов:
(Г) = ±М[МА (Г)] = )х /(А, х)с!х . 1 -1
Для гауссова процесса, обозначив /(х,х) = /(х,х,Г), имеем
(8)
Дх,х) =
2тфхД(1 -г2)
ехр
1
2(1 -г2)
(х-тх)2 Д.
2 г(х - тх )(х -т^) + (х-т±)2
А
(9)
о
о
о
где тх , тк — математические ожидания; /),.. /X — дисперсии процессов и Х(1); г —
коэффициент корреляции Х(/) и Х(1) соответственно.
Я*
г = г. =
л/^А
где — корреляционный момент (момент связи) Х(/) и Х(1).
Подставляя в (8) выражение для плотности вероятности /(Х,Х) (9), в котором заменен параметр х на Д, получим
1 00
ехр
I
2(1-г2)
(А-тх)2 2г(А-тх)(х-тх) (х-т,)
Д.
+ -
Д,
йх, (10)
В установившемся режиме при тх = тк = г = 0 интеграл (10) легко вычисляется и интенсивность выбросов определяется по формуле:
1
А д.
ехр
,2 Л
2Д
(11)
X у
Таким образом, задача вычисления интенсивности \'Л (/ ) сводится к определению установившейся дисперсии компоненты х() и дисперсии ее производной.
При достаточно большом значении А с (А > Зах, ох = ) выбросы стационарного процесса X () становятся редкими явлениями, а интервалы между выбросами будут настолько велики по сравнению с длительностью выбросов, что сечения случайного процесса, разделенные такими интервалами, будут практически независимыми. При таких предположениях закон распределения числа выбросов будет близок к пуассоновскому закону, для которого
т
(12)
где Рт — вероятность того, что число положительных выбросов за уровень А случайного процесса на интервале [0, /| с Травно числу т.
Вероятность отсутствия выбросов Р0 и вероятность хотя бы одного выброса 1\ на интервале [0, на основании (12) определяются выражениями
Р0 = ехр , Рх =1-ехр . (13)
Если условием работоспособности элемента системы является нахождение X() в диапазоне \Н, А ], где Н — допустимый уровень отрицательных выбросов Л"(/), то кроме вычисления Уд (Т) необходимо определять интенсивность отрицательных выбросов Ун(7). При этом для вычисления среднего числа отрицательных выбросов М Ын (Т) по (7) следует интервал интегрирования брать (-да, 0].
Выражения (11)—(13) позволяют оценить вероятность бессрывной работы системы управления на основе известных статистических характеристиках случайных процессов, описывающих ее работу. Эти характеристики могут быть получены путем математического моделирования или экспериментально.
о
Экспериментальная часть
В качестве примера рассмотрим систему второго порядка, описываемую уравнением х + ах + Ьх = £}, х(0) = 0, (14)
где Е, — белый шум с нулевым математическим ожиданием и заданной интенсивностью С. Система прекращает функционировать (происходит срыв управления) при достижении фазовой координатой х1 границы Д. Начальные значения х1 и х2 считаем нулевыми. Требуется определить вероятность бессрывного управления в течение заданного времени г.
На рис. 1 изображен график процесса х(1) до момента срыва управления . На рис. 2 представлен график ПРВ не поглощенных реализаций, находящихся в области (/) и не перешедших в область (г).
х {$ )
▲ ('")
-►
Л
г
Рис. 1. График процесса х(1) до момента срыва управления
Рис. 2. График ПРВ не поглощенных реализаций
Представим уравнение (14) в виде системы, обозначив х = хх, х = х2.
=х2, Х1(0) = Х10, (15)
х2 = -ах2 —Ьх1+'%,т± 0 = т±0. (16)
Для процесса, описываемого уравнениями (15)-(16), дифференциальные уравнения для входящих в выражение (11) дисперсий имеют вид:
Д1(О = ё11= 20125£)я О = £>я0,
Дет =612 =е22-2^®0е12-Ю0еП' Кк 0 " КхО '
(17)
(18)
А = 022 = -4^ш0е22 - 2ю02е12+о, е22 (0) - е220. (19)
В установившемся режиме тх—т± — = 0 и входящие в формулу (11) значения /),. и на основании уравнений (17)—(19) вычисляются по формулам:
А = еп = —-т' А = е22 = —— •
Например, при значениях параметров системы а=4, 6=9000, (7= 1 и времени работы системы /=20 с, при значениях А =3 а,.. А =4 а,.. А =5 а,, в соответствии с формулой (11), интенсивности срывов управления соответственно равны уд (Т) =0,168, уд (Т) =0,00507, Уд (Г) =0,0000563. Вероятность срыва управления (хотя бы одного выброса), вычисленная
по формуле (12) для приведенных выше трех уровней, равна Р\( Д=3 <зх )=0,965, Л( А =4ох )=0,0960, Л( А =5 ах )=0,00113 соответственно.
Таким образом, очевидно, что уровень 5 ах в данном случае является вполне подходящим, так как при нем вероятность срыва управления является допустимо малой.
Заключение
Предложенная аналитическая математическая модель вероятностного анализа стохастической системы позволяет на основе экспериментально или аналитически определенных статистических характеристик исследуемого процесса и заданных эксплуатационных параметров элементов системы (границ работоспособности), подверженных случайным воздействиям, определить диапазон работоспособности и вероятностные характеристики безотказной работы системы.
ANALYTICAL MODELING BOUNDARY REGIMES OF STOCHASTIC SYSTEMS
A.A. LOBATY, G.M. SAID Abstract
On the basis of the Marcov processes theory of casual structure the problem of the likelihood analysis of boundary operating modes of stochastic dynamic system dares. Expressions for intensity of absorption of realisations of casual process are analytically received. It allows to estimate probability of non-failure operation of a control system under known statistical characteristics of phase coordinates.
Литература
1. Казаков И.Е., Артемьев В.М., Бухалев В.А.. Анализ систем случайной структуры. М., 1993.
2. Пугачев В.С., Синицин И.Н. Теория стохастических систем. М., 2004.
3. Казаков И.Е., Мальчиков С.В. Анализ стохастических систем в пространстве состояний. М., 1983.
4. Тихонов В.И. Выбросы случайных процессов. М., 1970.
