Шалаев Ю.Н. каф. Информатики и проектирования систем. Институт кибернетики Теория случайных функций Случайной функцией называется случайная величина, зависящая от параметра t, т.е. X(t, ).

Если параметр t – время, то случайную функцию называют случайным процессом. Для дискретного случая – случайной последовательностью. Если зафиксировать элементарное событие = 0, то X(t, 0 ) будет неслучайной функцией аргумента t. Конкретный вид случайной функции при фиксированном в данном опыте называется реализацией случайной функции X( I,t). Если зафиксировать параметр случайной функции при t=t k, то она будет зависеть только от элементарного события и, следовательно, станет случайной величиной X(t к, ). При дальнейшем изложении аргумент для краткости опускается. Если зафиксировать элементарное событие = 0, то X(t, 0 ) будет неслучайной функцией аргумента t. Конкретный вид случайной функции при фиксированном в данном опыте называется реализацией случайной функции X( I,t). Если зафиксировать параметр случайной функции при t=t k, то она будет зависеть только от элементарного события и, следовательно, станет случайной величиной X(t к, ). При дальнейшем изложении аргумент для краткости опускается.

Законы распределения случайных функций Случайную функцию рассматривают как многомерную случайную величину. То есть X(t) можно представить как систему случайных величин: {X(t1), X(t2), …, X(tn)}, t1

F 2 (x 1, x 2 ; t 1,t 2 )=P{X(t 1 )

Для непрерывных случайных функций X(t) плотности распределения находятся как f 1 (x 1 ; t 1 )= F 1 (x 1 ; t 1 )/x 1, f 2 (x 1, x 2 ; t 1, t 2 )= F 2 (x 1, x 2 ;t 1,t 2 )/x 1x 2, f n (x 1, x 2,…x n ; t 1, t 2 …t n )= n F n (x 1, x 2,…x n ; t 1, t 2 …t n )/x 1…x n Их называют плотностями распределения первого, второго и n –го порядка.

Для плотностей распределения случайной функции X(t) имеет место интегральные соотношения: Для плотностей распределения случайной функции X(t) имеет место интегральные соотношения: Для независимых сечений X(t) n-мерная плотность вероятностей выразится через одномерную плотность распределения вероятностей формулой

Характеристики случайных функций Математическое ожидание случайной функции Математическим ожиданием случайной функции X(t) называется неслучайная функция m x (t), которая при каждом значении аргумента t равна математическому ожиданию соответствующего сечения случайной функции: Математическое ожидание случайной функции X(t) представляет собой некоторую среднюю функцию, около которой группируются и относительно которой колеблются все возможные реализации случайной функции.

Свойства математического ожидания Математическое ожидание неслучайной (детерминированной) функции равно самой этой функции MC(t)=C(t). Неслучайную функцию можно выносить за оператор математического ожидания MC(t)X(t)=C(t)MX(t). M(X(t)Y(t))=MX(t)MY(t). Для некоррелированных X(t) и Y(t) MX(t)Y(t)=MX(t)MY(t).

Дисперсия случайной функции Дисперсией случайной функции X(t) называется неслучайная функция D x (t), значение которой для каждого t равно дисперсии соответствующего сечения случайной функции Дисперсия характеризует разброс возможных реализаций случайной функции относительно среднего. Обладает свойствами дисперсии случайной величины.

Корреляционная функция случайной функции X(t) Корреляционной функцией случайной функции X(t) называется неслучайная функция двух аргументов t 1 и t 2, которая при каждой паре значений t 1 и t 2 равна корреляционному моменту соответствующих сечений случайной функции: После преобразования K x (t 1,t 2 ), получим K x (t 1,t 2 )=M(X(t 1 )X(t 2 ))-m x (t 1 )m x (t 2 ).

Корреляционная функция находится как Для независимых сечений случайной функции корреляционная функция равна нулю.

Свойства корреляционной функции Дисперсия случайной функции находится при равенстве аргументов t 1 =t 2 =t: DX(t)=K x (t,t). Корреляционная функция вещественной случайной функции симметричная функция: K x (t 1,t 2 )= K x (t 2,t 1 ). K x (t 1,t 2 ) убывает по мере увеличения длины интервала (t 1,t 2 ). Для вещественной случайной функции

Линейные преобразования случайных функций Прибавление неслучайного слагаемого Пусть X(t) – случайная функция, а C(t) – неслучайная функция: Y(t)=X(t)+C(t). Математическое ожидание: MY(t)=m x (t)+C(t). Корреляционная функция: K y (t 1,t 2 )= K x (t 1,t 2 ).

Умножение на неслучайный множитель Рассмотрим случайную функцию Y(t)=X(t)*C(t). Математическое ожидание: MY(t)=MX(t)*MC(t)=C(t)*MX(t). Корреляционная функция: K y (t 1,t 2 )=C(t 1 )C(t 2 )K x (t 1,t 2 ).

Дифференцирование случайной функции Пусть X(t) – случайная функция и заданы математическое ожидание m x (t) и корреляционная функция K x (t 1,t 2 ). Найдем характеристики случайной функции Y(t)=dX(t)/dt. Математическое ожидание: MY(t)=dm x (t)/dt. Функция корреляции: K y (t 1,t 2 )=K x (t 1,t 2 )/t 1t 2.

Интегрирование случайной функции Пусть X(t) – случайная функция и заданы математическое ожидание m x (t) и корреляционная функция K x (t 1,t 2 ). Найдем характеристики случайной функции Математическое ожидание Функция корреляции

Сложение случайных функций Рассмотрим сумму случайных функций: Z(t)=X(t)+Y(t). Найдем характеристики Z(t): По теореме сложения математических ожиданий получим: m z (t)=M((X(t)+Y(t))=m x (t)+m y (t)). Из определения корреляционной функции: После преобразования, получим:

Взаимная корреляционная функция: Корреляционная функция связи характеризует степень зависимости значения случайной функции X(t) взятого в момент t 1, от значения случайной функции Y(t) взятого в момент t 2. K x (t 1,t 2 ) – корреляционная функция связи одной случайной величины, поэтому иногда ее называют «автокорреляционной функцией».

Свойства взаимной корреляционной функции Для действительных случайных функций перестановка индексов при одновременной перестановке аргументов не меняет значения взаимной корреляционной функции: R xy (t 1,t 2 )=R yx (t 2,t 1 ). Взаимная корреляционная функция не изменяется при прибавлении любых неслучайных слагаемых. Если взаимная корреляционная функция равна нулю: R xy (t 1,t 2 )=0, то функции X(t) и Y(t)) называются некоррелированными (несвязанными).

Сложение случайной функции со случайной величиной Пусть X(t)-случайная функция, ξ-случайная величина; они некоррелированы. Получим случайную функцию Z(t)=X(t)+ ξ и определим ее характеристики при известных m x (t) и K x (t 1,t 2 ): MZ(t)=m z (t)=m x (t)+M ξ. Для корреляционной функции получим: K z (t 1,t 2 )=K x (t 1,t 2 )+K ξξ = K x (t 1,t 2 )+Dξ.

Нормированная взаимная корреляционная функция связи Это безразмерная характеристика связи между случайными функциями:

Оценка характеристик случайной функции Пусть имеется n реализаций случайной функции X(t): x 1 (t),x 2 (t),…, x n (t). Требуется найти оценки характеристик случайной функции: m x (t), D x (t) и K x (t 1,t 2 ). Для этого рассмотрим ряд сечений X(t) для моментов t 1,t 2,…, t m. Каждому из этих моментов будет соответствовать n значений случайной функции. Моменты задаются обычно равноотстоящими или из технических условий. Значения X(t) заносятся в таблицу:

t X(t) t 1 t 2 t m X 1 (t) x 1 (t 1 ) x 1 (t 2 ) … x 1 (t m ) X 2 (t) x 2 (t 1 ) x 2 (t 2 ) … x 2 (t m ) … … … … … X n (t) x n (t 1 ) x n (t 2 ) … x n (t m )

Каждое сечение t k есть n значений случайной величины и оценка математического ожидания находится по известному соотношению: Получим m точек оценки функции математического ожидания. Это усреднение по множеству реализаций. Оценка не смещена, так как Mm(t)=m(t). Для дисперсии:

Для корреляционной функции: По полученным значениям можно построить функцию математического ожидания и дисперсии по точкам, а функция корреляции двух аргументов воспроизводится по ее значениям в прямоугольной сетке. Эти функции можно аппроксимировать аналитическими выражениями.

Стационарные случайные функции Различают стационарность случайной функции в узком и в широком смысле. Стационарность в узком смысле случайной функции называется такая случайная функция X(t), для которой n-мерная плотность распределения вероятностей f n (x 1, x 2,…x n ; t 1, t 2 …t n ) при любом n зависит только от величины интервалов t 2 -t 1, t 3 -t 1,…, t n -t 1 и не зависит от положения этих интервалов в области изменения аргумента t.

Так, одномерная плотность распределения вероятностей f 1 (x;t) стационарной в узком смысле не будет зависеть от t – f(x). Двумерная будет зависеть от разности t 2 -t 1 = то есть f 2 (x 1,x 2 ;). N-мерная будет зависеть только от их разностей 1 = t 2 -t 1, 2 = t 3 - t 1,…, n-1 = t n -t 1, то есть f n (x 1, x 2,…x n ; 1, 2 … n-1 ).

Стационарная функция называется стационарной в широком смысле, если ее математическое ожидание постоянно, а корреляционная функция зависит только от разности аргументов =t 2 -t 1 : m x (t)=const; K x (t 1,t 2 )=K x ().

Стационарная функция в широком смысле может быть нестационарной в узком смысле. Наоборот, случайная функция стационарная в узком смысле, является стационарной в широком смысле:

Свойства K x () Если X(t) вещественная и стационарная, то ее корреляционная функция является четной функцией: K x ()= K x (-). Если X(t) стохастически непрерывна, то ее корреляционная функция K x () есть функция непрерывная. Если X(t) вещественная и стационарная, то имеет место неравенство |K x ()|K x (0).

Эргодические свойства стационарных случайных функций Так как X(t) вещественная и стационарная и процесс протекает однородно по времени, то по одной реализации достаточной продолжительности можно оценить характеристики случайной функции. Эргодическое свойство состоит в том, что каждая отдельная реализация X k (t) несет как бы информацию всей совокупности возможных реализаций X(t), т.е. одна реализация достаточной продолжительности может заменить при обработке множество реализаций.

Оценка характеристик стационарной случайной функции Если X(t) обладает эргодическим свойством, то для нее среднее по времени приближенно равно среднему по множеству наблюдений. Для m x (t)=const:

Для корреляционной функции: где Тогда Вычислив интеграл для ряда, можно приближенно воспроизвести по точкам весь ход корреляционной функции.

На практике интегралы заменяют конечными суммами. Для этого интервал Т разбивают на n равных частей длиной t=T/n.

Обозначим середины полученных участков t 1, t 2,…, t n, тогда получим: Для корреляционной функции введем: =mt=mT/n. Для интервала интегрирования: T-=T-mT/n=(n-m)/n·T, 1/(T-)=n/(n-m)·T. Тогда получим оценку корреляционной функции:

m =0, 1, 2,…. Вычисления проводятся до тех m при которых корреляционная функция становится равной нулю.

Марковский случайный процесс Случайный процесс называется Марковским, если все вероятностные характеристики процесса в будущем зависят лишь от того, каким образом этот процесс протекал в прошлом. То есть будущее зависит от прошлого только через настоящее.

Пусть A={a 1, a 2,…, a n } – пространство исходов эксперимента или пространство состояний некоторой системы, одинаковые для каждого шага случайного процесса. Тогда по определению Марковского процесса вероятность того, что система переходит из состояния a i в состояние a k зависит только от состояния a i из которого она исходит в процессе рассматриваемого перехода и не зависит от того, каким образом система пришла в это состояние.

Случайный процесс, протекающий в физической системе называется цепью Маркова, если переходы системы из одного состояния в другое возможны только в определенные дискретные моменты времени t 1, t 2,…. Обозначим через P ik вероятности перехода системы из состояния a i в состояние a к. Марковская цепь характеризуется тем, что вероятности P ik определяются для всех упорядоченных пар состояний и задано исходное состояние.

Вероятности перехода Р 11, Р 22, Р 33 означают, что система остается в состоянии а 1 с вероятностью Р 11, в состоянии а 2 Р 22, а 3 Р 33. Равенство нулю вероятности P ik означает невозможность соответствующих переходов. Вероятности перехода из одного состояния в другое можно представить двумя способами.

Первый способ состоит в том, что вероятности перехода записываются в виде квадратной матрицы. Для Марковской цепи с тремя состояниями а 1, а 2, а 3 матрица имеет вид: Матрица вероятностей перехода (стохастическая матрица) должна обладать условием, что сумма элементов каждой строки равна единице.

Второй способ представления вероятностей перехода состоит в построении диаграммы перехода, когда возможные состояния системы S наглядно изображают с помощью графа состояний. Возможные состояния системы на графе изображаются окружностями. Вероятности перехода изображаются ребрами с соответствующими числами вероятностей перехода. Сумма вероятностей для ребер, выходящих из любой вершины графа должна равняться единице. По матрице можно построить граф и наоборот.

При изучении Марковских цепей иногда возникает задача: найти вероятности того, что через n шагов процесс перейдет из состояния a i в состояние a k. Зная матрицу P n можно найти матрицу P n+1 по соотношению: P n+1 =P·P n.

Марковская цепь называется регулярной, если какая – либо степень ее матрицы вероятностей перехода не содержит нулевых элементов. Любая стохастическая матрица, не содержащая нулей, определяет регулярную Марковскую цепь.

Марковская цепь называется эргодической, если из каждого ее состояния можем попасть в любое другое состояние.

Моделирование случайных величин Случайные величины моделируют с помощью преобразований одного Или нескольких независимых Значений случайной величины L, равномерно распределенной в интервале (0,1).

Моделирование дискретных случайных величин Общий метод моделирования основан на следующем очевидном равенстве которые связаны рекурсивными формулами и моделирование производится по схеме

да нет M=L, m=0, P=P0 M=M-P M

Для биномиального распределения с параметрами (P,n)

Для распределения Пуассона с параметром a

Моделирование непрерывных случайных величин Случайная величина моделируется по формуле вида: ξ=φ(L), где φ(L) – строго монотонная и непрерывная функция на интервале (0,1). Для случайной величины ξ задана плотность распределения f(x) на интервале (а, в). Для монотонно возрастающей функции φ(L) моделирующая формула для непрерывной случайной ξ запишется в виде

Для монотонно убывающей функции φ(L) моделирующая формула для непрерывной случайной ξ запишется как которая эквивалентна предыдущей формуле т. к. случайные величины 1-L и L одинаково распределены.

Для экспоненциального закона Плотность вероятностей экспоненциального закона имеет вид Функция распределения имеет вид Из этих соотношений определяется формула для Моделирования экспоненциального закона