ТОПИН.Лекции, задания / ТОППиН_часть1 / глава 1234
.pdf
|
|
|
|
|
|
( y my )2 |
|
f ( y,x) |
1 |
|
|
22y |
d y, |
||
|
|
|
|||||
|
|
|
e |
|
|
||
|
|
|
|
|
|||
|
y |
2 |
|
|
|
||
где my – математическое ожидание; σy – среднеквадратическое отклонение. Случайную функцию можно также определить как функцию своего аргу-
мента, значение которой при любом значении ее аргумента является случайной величиной.
Характеристиками случайных функций, в отличие от числовых характеристик случайных величин, представляющих собой определенные числа, являются в общем случае не числа, а функции.
Математическое ожидание M[y(x)] случайной функции y(x) есть неслучайная функция my(x), которая при каждом значении аргумента x равна математическому ожиданию (МО) соответствующего сечения случайной функции
[ y(x)] my (x) y f ( y,x) dy,
где f(y, x) – одномерная плотность вероятностей семейства реализаций случайных функций.
По смыслу МО случайной функции есть некоторая средняя функция, около которой группируются и относительно которой колеблются все реализации случайных функций.
Дисперсия случайной функции y(x) есть неслучайная функция Dy(x), значение которой для каждого x равно дисперсии соответствующего сечения случайной функции
|
|
D[ y(x)] Dy (x) |
[ y my (x)]2 f ( y,x) dy. |
|
|
Дисперсия Dy(x) характеризует разброс реализаций в семействе относительно МО my(x).
Среднеквадратичное отклонение случайной функции
y (x) 
Dy (x).
Одномерная плотность распределения f(y, x) позволяет найти среднее значение случайной функции и ее дисперсию, но не содержит информации о поведении случайной функции. Например, две случайные функции обладают одинаковыми средними значениями, но одна из них представляет плавно меняющиеся случайные функции, а другая сильно перемешанные. У первой из них сильно связанные значения y для сечений x1 и x2 (рис. 4.3, а), у другой эти значения независимы (рис. 4.3, б).
51
а |
б |
Рис. 4.3
Вводится двумерная плотность распределения случайной функции, которая показывает насколько связаны между собой значения случайной функции при двух различных значениях аргумента. Для количественной оценки этой связи (изменения случайной функции с изменением аргумента x) используется
автокорреляционная функция (АКФ).
АКФ семейства реализаций случайной функции называется детерминированная функция двух аргументов, которая при каждой паре значений аргумента x1 и x2 равна корреляционному моменту соответствующих значений семейства реализаций случайной функции.
K y (x1,x2 ) |
|
|
|
|
|
|
||
|
|
|||||||
y(x1) my (x1) |
y(x2 ) my (x2 ) |
|||||||
|
|
|
|
|
|
|
||
f y(x1),x1; y(x2 ),x2 |
dyx |
dyx . |
|
|||||
|
|
|
|
|
1 |
2 |
|
|
Функция f [ y(x1),x1; y(x2 ),x2 ] в данном случае представляет собой дву-
мерную плотность вероятностей y(x).
АКФ содержит данные из двух сечений семейства реализаций случайной функции. Следовательно, она определяется двумерными свойствами семейства реализаций. Назначение АКФ – характеризовать степень зависимости между сечениями случайной функции, относящимся к различным значениям аргумента х.
Взаимосвязь вероятностных свойств двух случайных функций y(x1) и z(x2)
обычно характеризуется взаимной корреляционной функцией (ВКФ).
B(x1,x2 ) |
|
|
|
|
|
|
|
|
|||
y(x1) my (x1) z(x2 ) mz (x2 ) |
|||||
|
|
|
|
|
|
|
f y(x1),x1; z(x2 ),x2 dyx |
dzx , |
|||
где f y(x1), x1; z(x2 ), x2 |
|
|
|
1 |
2 |
– совместная плотность вероятностей значений |
|||||
y(x1) и z(x2).
Случайный процесс есть случайная функция от независимой переменной t- времени: y(t).
52
Следовательно, если в качестве аргумента используется время, то имеем дело со случайными процессами.
Реализации случайного процесса называются траекториями. Все что говорилось о случайных функциях полностью относится к случайным процессам с аргументом t. Характеристики случайных процессов аналогичны характеристикам случайных функций: МО, D, σ, АКФ, ВКФ случайных процессов.
В качестве дополнительной характеристики для стационарных случайных процессов (процессов, которые протекают во времени приблизительно однородно) используют интервал корреляции
k Ry ( ) d ,
0
где Ry(τ) – нормированная АКФ процесса.
Это позволяет более просто, хотя и более грубо, описать быстротечность процесса.
Величина τk дает ориентировочное представление о том, на каких интервалах времени в среднем имеет место существенная корреляция между значениями случайного процесса. Максимальный интервал корреляции τkмакс – такой интервал, за пределами которого корреляция между значениями случайного про-
цесса пренебрежимо мала, т.е. Ry k макс , где α= 0,1,…, 0,01 – некоторая малая, наперед заданная величина.
Втеории случайных процессов принято классифицировать их по тем или иным признакам, учитывая плавность или скачкообразность реализации, фиксированность или случайность моментов, в которые могут происходить скачки
ит.д., вид закона распределения отдельного сечения процесса или совокупности его сечений и т.д. Самая элементарная классификация случайных процессов
– «по времени» и «по состояниям».
Случайный процесс с дискретным временем – это процесс в системе, кото-
рая может менять свои состояния только в отдельные моменты времени t1, t2 и т.д.
Случайный процесс с непрерывным временем – если переходы из одного состояния в другое могут происходить в любой момент времени.
Случайный процесс с непрерывными состояниями, если его сечение в любой момент времени представляет собой непрерывную случайную величину и, значит, множество ее значений несчетно.
Случайный процесс с дискретными состояниями, если в любой момент времени множество его состояний конечно или счетно, т.е. случайная величина в сечении дискретна.
Всоответствии с вышеизложенным приведем одну из классификаций случайных процессов (четыре класса):
I. Параметр t дискретен и множество состояний дискретно. II. Параметр t непрерывен, множество состояний дискретно. III. Параметр t дискретен, множество состояний непрерывно.
IV. Параметр t непрерывен и множество состояний непрерывно.
53
К I типу можно отнести дискретные цепи Маркова, последовательности дискретных случайных величин. К II типу – Марковские цепи с непрерывным множеством состояний. Тип III представляет процесс Пуассона. В IV тип входит броуновское движение, сюда же относятся случайные процессы общего типа.
На рис. 4.4 представлены диаграммы:
а) дискретной случайной последовательности (х и t – дискретны);
б) дискретного случайного процесса (х – дискретно, t – непрерывно);
в) непрерывной случайной последовательности (х – непрерывно, t – дискретно);
г) непрерывного случайного процесса (х и t – непрерывны).
xk |
xk |
0 1 2 3 4 5 6 7 |
tk |
t |
а |
|
б |
х |
|
х |
|
|
|
|
|
|
|
|
|
tk |
tk |
|
|
|
|
|
|
|
|
|
||
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
|
|
|
|
|
|
|
|
|
в |
|
|
г |
|
|
|
|
|
|
|
|
|
Рис. 4.4 |
|
Все процессы можно также классифицировать и из других соображений, например:
1.стационарные – нестационарные;
2.гауссовские – негауссовские;
3.марковские – немарковские
ит.д.
Некоторые из случайных процессов (стационарный, Марковский и др.) широко используются в теории надежности, так как позволяют с достаточной точностью описать процессы, происходящие в элементах и системах.
4.2Марковский случайный процесс. Дифференциальные уравнения Колмогорова
Основными математическими методами анализа функционирования систем являются марковские процессы.
Марковские процессы могут быть процессами, как с дискретным, так и с непрерывным временем.
54
Важным частным случаем марковского процесса является процесс Пуассона. Основные виды марковских процессов:
–марковская цепь (дискретный марковский процесс с дискретным временем);
–марковская последовательность (непрерывный марковский процесс с дискретным временем;
–разрывной марковский процесс (дискретный марковский процесс с непрерывным временем);
–непрерывный марковский процесс (непрерывный марковский процесс с непрерывным временем);
–дискретный марковский процесс (дискретно-непрерывный марковский процесс).
Среди марковских процессов выделяются процессы, которые получили название «процессы гибели и размножения».
Марковский случайный процесс – это процесс, у которого в каждый момент времени вероятность любого состояния объекта в будущем зависит только от его состояния в настоящий момент времени и не зависит от того, каким образом объект пришел в это состояние. Следовательно, вероятность будущего состояния объекта не зависит от вероятности прошлого состояния, а только от вероятности настоящего состояния, состояния в данный момент времени. Таким образом, марковский процесс – это случайный процесс без последействия.
Рассмотрим систему, которая может находиться в трех состояниях (рис. 4.5). Это марковский процесс с дискретными состояниями и непрерывным временем.
L(t)
α1
α2
α3
t1 |
t2 |
t |
Рис. 4.5
Процесс L(t) может принимать только дискретные значения α1, α2, α3. Моменты времени t1, t2 – моменты переходов из одного состояния в другое. Переход из одного состояния в другое происходит за время t. Не может быть перехода из одного крайнего состояния в другое, крайнее состояние, например, из а1
в а3, и наоборот.
Наиболее подходящим способом представления процесса перехода системы в различные состояния является граф переходов (состояний) (рис. 4.6), в котором вершины (узлы) соответствуют состояниям системы, а дуги указывают все возможные переходы системы из одного состояния в другое.
55
Рис. 4.6
На рисунке обозначены: P |
( t),P |
ji |
( t) – вероятности переходов из одно- |
i j |
|
|
го состояния в другое за время t; Pii ( t) – вероятность сохранения i-го состояния за время t; i j , ji – интенсивности прямого и обратного переходов из
одного состояния в другое за время t.
Вероятность перехода – это важнейшая характеристика марковского процесса.
Возможно представление вероятностного процесса, описываемого графом переходов, матрицей вероятностей переходов (или просто переходов)
|
|
|
|
P00 P01 0 |
|
|
||
A |
|
|
|
P |
P |
P |
|
. |
|
|
|
|
10 |
11 |
12 |
|
|
|
|
|
|
0 P21 P22 |
|
|
||
Состояние системы в момент времени t определяется вектором P(t) с составляющими P0, P1, и P2 – вероятности состояний:
P0 (t)
P(t) P (t) .
1
P 2 (t)
Для момента времени t +Δt состояние системы будет определяться вектором P(t t) с такими же составляющими:
P0 (t t)
P(t t) P (t t) .
1
P 2 (t t)
Переход из состояния системы для момента времени t в состояние для момента t +Δt осуществляется преобразованием вектора P(t) в вектор P(t t) с
помощью транспонированной матрицы переходов AT:
P(t t) P(t) AT .
Академик А.Н. Колмогоров предложил систему дифференциальных уравнений для определения вероятностей каждого из состояний системы (дифференциальные уравнения типа массового обслуживания). Рассмотрим получение дифференциальных уравнений Колмогорова на примере системы, которая может находиться в двух состояниях. Граф переходов (состояний) будет иметь вид (рис. 4.7).
56
Рис. 4.7
Матрица переходов |
|
|
|
A |
|
|
|
|
P00 ( t) P01( t) |
|
|
|
|
. |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
P |
|
( t) P |
|
( t) |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
11 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Определим компоненты вектора P(t t) : |
|
|
|
|
|
||||||||||||||||||||||||||||
|
P (t t) |
|
|
|
|
P (t) |
|
AT |
|
|
P (t) |
|
|
|
|
P ( t) P ( t) |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
0 |
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
00 |
10 |
|
|
. |
|||||||||||||
|
P (t t) |
|
|
|
|
|
P (t) |
|
|
|
|
|
|
P (t) |
|
|
|
|
|
|
|
P ( t) P ( t) |
|
|
|
||||||||
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
01 |
11 |
|
|
|
|
P (t t) P (t) P ( t) P (t) P ( t) |
|
|||||||||||||||||||||||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
0 |
|
00 |
|
|
|
|
|
|
|
|
|
1 |
10 |
|
|
|
|||||
|
P (t t) P (t) P ( t) P (t) P ( t) . |
|
|||||||||||||||||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
0 |
01 |
|
|
|
|
|
|
|
|
|
1 |
11 |
|
|
|
||||||
Так как P ( t) 1 P |
( t) , то получаем уравнения для определения веро- |
||||||||||||||||||||||||||||||||
ii |
|
|
i j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ятностей состояний системы |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
P (t t) P (t)[1 P |
|
|
( t)] P (t) P ( t) |
|
||||||||||||||||||||||||||||
|
0 |
|
|
|
|
|
|
0 |
|
|
|
|
01 |
|
|
|
|
|
|
|
|
1 |
10 |
|
|
|
|||||||
|
P (t t) P (t) P |
|
( t) P (t)[1 P ( t)]. |
|
|||||||||||||||||||||||||||||
|
1 |
|
|
|
|
|
|
0 |
|
|
01 |
|
|
|
|
|
|
|
1 |
|
10 |
|
|
|
|||||||||
Данную систему уравнений можно записывать по графу переходов по следующему правилу (без предыдущих преобразований).
Вероятность того, что система в момент времени t +Δt будет находиться в 0-м состоянии равна вероятности того, что система в момент времени t находилась в 0-м состоянии и не перейдет из него за время t в первое состояние плюс вероятность того, что система в момент времени t находилась в первом состоянии и за время t перейдет из него в 0-е состояние (сумма вероятностей ситуаций).
Вероятность 1-й ситуации определяется произведением P0 (t)[1 P01( t)],
вероятность 2-й ситуации определяется произведением P (t) P |
( t) . |
|
1 |
10 |
|
Для второго уравнения аналогично.
Вероятности переходов с достаточной точностью можно выразить через интенсивности переходов (e-аt при аt 0,1, разлагаем в ряд 1 – аt).
P ( t) 1 e ij t |
1 1 |
ij |
t |
ij |
t, |
|
ij |
|
|
|
|
||
Pji ( t) 1 e |
ji t |
1 1 ji t ji t. |
||||
|
||||||
Тогда уравнения состояний примут вид:
P (t t) P (t) (1 |
01 |
t) P (t) |
|
t |
||||
0 |
0 |
|
|
1 |
10 |
|
||
P (t t) P (t) |
01 |
t P (t) (1 |
t) . |
|||||
1 |
0 |
|
|
1 |
10 |
|
|
|
57
В |
каждом |
уравнении |
Pi (t) |
переносим |
в |
|
левую часть, получаем |
|||||||||
P (t t) P (t) P (t) и уравнения можно записать: |
|
|
||||||||||||||
i |
i |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
P (t) P (t) |
01 |
t P (t) |
|
t |
||||||||||
|
|
|
0 |
|
|
0 |
|
|
1 |
|
10 |
|
||||
|
|
|
|
|
|
|
|
|
t P (t) |
|
t . |
|||||
|
|
P (t) P (t) |
01 |
|
||||||||||||
|
|
|
1 |
|
0 |
|
|
|
1 |
10 |
|
|
||||
Разделив левую и правую части уравнений на |
|
t и устремляя t к 0, имеем |
||||||||||||||
в левой части уравнений |
Pi (t) |
|
dPi (t) |
. |
|
|
|
|
|
|||||||
t |
|
d t |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
t0
Дифференциальные уравнения Колмогорова для однородного марковского процесса принимают вид:
dP0 (t)
dt
dP (t)
1
dt
01 P0 (t) 10 Pi (t)
01 P0 (t) 10 Pi (t) .
Чтобы решить дифференциальные уравнения Колмогорова и найти вероятности состояний, прежде всего надо задать начальные условия. Если мы точно знаем начальное состояние системы, то в начальный момент времени (при t = 0) Р0(0) = 1, а все остальные начальные вероятности равны нулю: Р1(0) = Р2(0) = 0.
Для любого момента времени выполняется нормировочное условие
n
Pi (t) 1 при t ≥0.
i0
Марковский процесс называется однородным, если закономерности его поведения на любом интервале времени t2 – t1, не зависят от расположения этого интервала на оси времени.
Дифференциальные уравнения можно составлять непосредственно по графу переходов по следующему правилу.
В левой части каждого уравнения записывается dPi (t) dt , а в правой части
столько членов сколько стрелок связано с данным состоянием. Каждый член равен произведению интенсивности перехода ( ij или ji ), переводящей си-
стему по данной стрелке, на вероятность того состояния, из которого исходит стрелка. Знак перед каждым членом зависит от направления перехода: если стрелка направлена в данное состояние (стрелка входит), то ставится знак плюс, если из данного состояния (стрелка выходит), то минус.
Можно пользоваться другим правилом: производная вероятности каждого состояния равна сумме всех потоков вероятности, идущих из другого состояния
вданное, минус сумма всех потоков вероятности, идущих из данного состояния
вдругое
dPi (t) |
n |
|
|
|
m |
|
|
|
|
P (t) P (t) |
|
, |
|||
|
|
|
|||||
dt |
j1 |
ji |
j |
i |
j1 |
ij |
|
58
где n – количество вершин графа (состояний), из которых стрелки (переходы) идут в данное (i) состояние; m – количество вершин графа (состояний), в которые направлены стрелки (переходы) из данного (i) состояния; Рi(t) – вероятность данного i-го состояния; Рj(t) – вероятности состояний, из которых стрелки (переходы) идут в данное (i) состояние.
Решение полученной системы уравнений, т.е. нахождение вероятностей состояний Р0(t) и Р1(t), производится с помощью преобразования Лапласа. При этом проводится замена
Р(t) → Р(z)
dP(t) z P(z) P(t 0). dt
Если число уравнений больше двух – трех, то обычно их решают численно
– вручную или на ЭВМ.
Решение существенно упрощается, если рассматриваемый процесс являет-
ся марковским стационарным процессом, для которого dPi (t) dt 0 , т.е. веро-
ятности состояний становятся постоянными Рi(t) = Рi = const. В этом случае система дифференциальных уравнений становится системой алгебраических уравнений.
Отметим, что описание функционирования системы можно выполнять с помощью интегральных уравнений [36].
4.3 Потоки событий
Поток событий – это последовательность событий, при которой они происходят одно за другим в случайные моменты времени. Случайный поток можно также назвать случайным точечным процессом, так как реализации такого процесса представляют собой случайную последовательность точек на временной оси.
Случайный точечный процесс, удовлетворяющий трем условиям: ординарности, стационарности и отсутствию последействий, называется простейшим или пуассоновским процессом или потоком.
Свойства простейшего потока событий.
Стационарность – свойство потока событий, заключающееся в том, что вероятность появления определенного числа событий за фиксированный промежуток времени не зависит от положения этого промежутка на оси времени, а зависит только он его длины, т.е. плотность потока появления событий постоянна во времени.
Ординарность – свойство потока, заключающееся в том, что события возникают поодиночке, а не группами, и практически невозможно одновременное появление двух и более событий.
Отсутствие последействия – свойство потока, заключающееся в том, что вероятность появления определенного числа событий в течение некоторого промежутка времени не зависит от числа и характера возникновения событий
59
до начала этого промежутка времени, т.е. события взаимно независимы и случайные промежутки времени между соседними событиями также взаимно независимы.
На практике не всегда одновременно выполняются все три свойства, но простейший поток очень удобен для решения различных задач.
Важным понятием является интенсивность потока (t) – это среднее число событий, приходящееся на единицу времени. Для стационарного потока событий его интенсивность есть величина постоянная, т.е. (t) = = соnst.
Простейший (стационарный пуассоновский) поток описывается распределением (законом) Пуассона
p(k) ak e a , где а = λ t, k!
т.е. p(k) ( t)k e t , (k 0,1, 2,...), k!
здесь k – число событий на интервале t; р(k) – вероятность появления k событий за интервал времени t.
Оценка параметра а выполняется следующим образом. Если п – число проведенных испытаний и ki – число событий, появившихся в i-м испытании, то
n |
k |
|
оценка максимального правдоподобия параметра а равна a |
i |
. |
|
||
i1 n
На рис. 4.8 приведено семейство распределений Пуассона для различных значений параметра а.
Для простейшего потока с интенсивностью интервал Т между соседними
событиями имеет показательное распределение с плотностью f (t) e t .
Если отсутствует свойство стационарности, то имеет место нестационарный пуассоновский поток ( – зависит от времени), у которого интервал между появлениями событий уже не подчиняются показательному закону распределения.
Стационарный поток с ограниченным последействием называется потоком Пальма (или рекуррентным). Для такого потока интервалы Т1, Т2, … между событиями представляют собой последовательности независимых случайных величин с одинаковым произвольным распределением. При показательном распределении интервалов между событиями поток Пальма становится простейшим потоком.
60
