Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лекции АСП 1-4

.pdf
Скачиваний:
31
Добавлен:
01.03.2023
Размер:
363.04 Кб
Скачать

АНАЛИЗ СЛУЧАЙНЫХ ПРОЦЕССОВ

Тема № 1 Основные понятия. Стационарные процессы

Рассмотрим произвольное вероятностное пространство ,A, . Вспомним

определение из курса теории вероятностей, призванное формально определить величины, подлежащие «измерению» в случайных экспериментах.

Определение. Случайной величиной называется функция w : ,

заданная на пространстве элементарных событий , принимающая значения в

и удовлетворяющая условию, что для любого вещественного х множество

w : w x является событием, т.е. w : w x A.

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

Основным понятием теории случайных процессов является следующее.

Определение. Случайным процессом (или случайной функцией) t,t T

называется семейство случайных величин t , зависящих от параметра t, пробегающего некоторое множество значений T . Таким образом, любому фиксированному t0 однозначно ставится в соответствие случайная величина t0

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

вероятностном пространстве ,A,P .

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

Определение. Реализацией, выборочной функцией или траекторией

случайного процесса

t,t T

называют действительную функцию

t 0 :T аргумента

t T , где

0 - некоторый фиксированный элемент из

.

На практике считается, что в результате некоторого эксперимента происходит элементарное событие 0 и наблюдается значение t t в

точке 0 . Поэтому траектория случайного процесса зависит от наблюдений при эксперименте. Иначе говоря, реализацией случайного процесса X(t) называется неслучайная функция x(t), в которую превращается случайный процесс X(t) в результате опыта; другими словами, конкретный вид, принятый случайным процессом X(t), который наблюдался на каком-то отрезке времени от 0 до τ.

Каждой сечению t

для фиксированного t T

однозначно ставится в

соответствие функция

распределения Ft x P t

x , которая является

полной и исчерпывающей характеристикой случайной величины t . Эта

функция

описывает поведение случайного процесса t,t T

в момент

времени

t

и называется одномерным законом распределения

случайного

процесса t

,t T .

 

Пример. Пусть T . По определению зададим t t t0 b, где t0,b

- некоторые фиксированные числа, - некоторая случайная величина, заданная на некотором вероятностном пространстве ,A,P , не зависящая от

параметров t, t0 , b.

 

 

 

 

 

 

 

 

При фиксированном t случайная величина t

будет обладать функцией

распределения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x b

 

 

 

 

 

F

 

 

 

,t t0

 

 

 

 

 

t0

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

x b

 

 

 

 

 

F x P

t

x

P (t t

) b x

P

 

 

 

,t t0 .

 

 

 

t

 

0

 

 

t t0

 

 

 

 

 

 

 

0,b x

 

 

 

 

 

 

 

 

 

 

 

P b x

 

 

,t t0

 

 

 

 

 

 

 

1,b x

 

 

 

 

 

 

 

 

 

 

 

 

x b

x b

 

x b

Заметим, что P

 

1 F

 

P

 

.

 

 

 

 

t t0

t t0

 

t t0

Если фиксировать элементарный исход 0

, то t 0 будет реализацией

случайного процесса t,t T

и будет представлять из себя прямую линию.

Тангенс угла наклона этой линии равен значению 0 . Так как случайная величина может принимать разные значения то все реализации (траектории) случайного процесса t будут принадлежать пучку прямых на плоскости. Каждая прямая этого пучка проходит через точку плоскости с координатамиt0,b . Эта точка в аналитической геометрии носит название центра пучка прямых на плоскости.

Случайный процесс t в данном примере носит называние веерного случайного процесса.

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

распределения случайного вектора t1 ,..., tn , координаты которого равны значениям процесса в моменты времени t1,...,tn .

Обозначим

Ft1,...,tn (x1,...,xn) F t1 ,..., tn (x1,...,xn) F (x1,...,xn) P t1 x1,..., tn xn .

Определение. Семейство функций

Ft1,...,tn (x1,...,xn), n , t1,...,tn Tn,ti tj

называется семейством конечномерных распределений случайного процесса

t,t T .

Это семейство удовлетворяет условиям согласованности:

 

 

 

 

 

 

i1,...,in -

 

 

 

 

 

 

 

 

 

 

1.

Пусть

 

перестановка

элементов

из

1,n,

тогда

Fti

,...,ti xi1

,...,xin Ft1,...,tn x1,...,xn .

 

 

 

 

 

 

 

 

 

1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

lim Ft ,...,t

 

,t

(x1,...,xn 1,xn) Ft ,...,t

(x1,...,xn 1).

 

 

 

 

 

 

 

 

xn

1

 

n 1

 

n

1

n 1

 

 

 

 

 

 

 

 

 

Введем еще одно полезное

 

 

 

 

 

 

 

 

 

 

Определение. Случайные процессы t,t T и t,t T , заданные на одном

вероятностном пространстве ,A,P

и одинаковом множестве T , называются

стохастически эквивалентными, если t T

P t t 1.

 

 

 

 

 

 

Определим некоторые характеристики случайного процесса { t,t T}:

 

 

1. Пусть

 

m t E t для

всех t T , если это

математическое ожидание

определено

 

и конечно,

тогда

функция

m t :T

 

называется

математическим ожиданием случайного процесса.

2. Пусть D t D t =E t m t 2 для всех t T , если дисперсия определена и конечна, тогда D t :T - дисперсия случайного процесса.

Наряду с дисперсией в инженерных приложениях теории случайных процессов часто рассматривается функция среднеквадратичного отклонения

t D t :T .

3. Пусть m t , m s , тогда K s,t E s m s t m t -

ковариационная функция случайного процесса. Синонимами понятия ковариационная функция в технической литературе являются автоковариационная и автокорреляционная функции. Согласно мнению ряда ученых наиболее удачным названием является «автоковариационная функция», подчеркивающее, что речь идет о ковариации двух сечений одного и того же процесса (от латинского autos – сам, направленный на самого себя).

Докажите следующие свойства K s,t :

а) K s,t 2 D s D t ,

б) K t,t D t ,

в) K s,t K t,s .

4. R s,t

 

K s,t

 

- корреляционная (нормированная корреляционная,

 

 

 

D s D t

 

 

 

 

нормированная автоковариационная, нормированная автокорреляционная) функция случайного процесса.

Докажите следующие свойства R(s,t):

а) R s,t 1,

б) R t,t 1,

в) R s,t R t,s .

Пример.

Рассматривается случайная

функция X(t) Ut2 2, где

U

случайная

величина, распределенная по

нормальному закону N(2; 4).

Нам

требуется найти функцию распределения сечения этой функции, математическое ожидание mX (t), дисперсиюDX (t), среднеквадратичное отклонение X (t) и корреляционную функцию K X (t1,t2).

Согласно определению функции распределения случайной величины

FX (t) x P X(t) x P Ut2 2 x

 

 

 

 

 

 

x 2

 

x 2

P U

 

 

,t 0

FU

 

 

 

,t 0

t

2

 

t

2

 

 

 

 

 

 

 

 

 

 

 

0,x 2,t 0

 

0,x 2,t 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,x 2,t 0

1,x 2,t 0

 

 

 

.

Далее рассмотрим случай, когда t 0. Так как U N(2; 4),то

 

 

 

 

 

y

x 2

 

 

 

 

 

 

z 2 2

 

x 2

 

t2

 

1

 

 

 

 

FU y pU z dz

 

 

 

FU

 

 

 

 

 

 

 

 

e

 

8 dz

t

2

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

По определению

mX (t) E X t E Ut2 2 t2 EU 2,

 

 

 

DX (t) DX t D Ut2 2 t4 DU ,

 

 

 

 

X (t)

 

t2

 

,

 

 

 

 

 

 

 

 

 

DX t

 

 

 

 

 

 

 

 

 

DU

 

 

 

 

 

K X t1,t2 E X t1 mX t1 X t2 mX t2

E Ut12 2 t12 EU 2 Ut22 2 t22 EU 2 t12t22 DU .

Параметры нормального распределения являются его основными числовыми характеристиками (EU 2,DU 4). Следовательно,

mX (t) 2t2 2, DX (t) 4t4 ,X (t) 2t2 ,

K X t1,t2 4t12t22 .

Классификация случайных процессов

Основными признаками, по которым можно провести классификацию

случайных процессов, являются:

 

 

 

 

1). Пространство состояний E,

которому

принадлежат все

возможные

значения, принимаемые всеми случайными величинами t

по всем t T :

E Et

 

 

 

 

t

( ) .

 

 

t T

t T

 

 

 

Назовем множество E фазовым пространством.

 

 

Если E , то E – фазовое пространство целочисленного

случайного

процесса.

Если E – это некоторый непрерывный промежуток из , то E – фазовое пространство действительного случайного процесса.

2). Временной параметр T .

Определение. Если t T , то говорят, что процесс t,t T является

процессом с дискретным временем. Если параметр t принимает все значения из некоторого непрерывного промежутка на числовой прямой, то говорят, что

процесс t,t T

- случайный процесс с непрерывным временем.

 

3). Отношение зависимости между случайными величинами t .

 

Определение.

Рассмотрим

случайный

процесс t,t T , где

T . Для

всех n для

любых t1 t2

... tn T

рассмотрим случайные величины

t2 t1 ,..., tn tn 1 .

Случайный процесс называется процессом с независимыми приращениями,

если величины t2 t1 ,..., tn tn 1 независимы в совокупности.

Случайный процесс называется процессом с некоррелированным

приращением, если величины t

2

t

,..., t

t

некоррелированные

 

1

n

n 1

 

 

случайные величины.

Случайный процесс называется процессом с ортогональным приращением,

если t

2

t

,..., t

t

некоррелированы, и m t 0

для всех t T .

 

1

n

n 1

 

 

Случайный процесс называется марковским, если при всех n для любых

t1 t2 ... tn tn 1 T ,

для

любых x1,...,xn,xn 1 E,

где

E

фазовое

пространство случайного процесса t,t T , выполняется

 

 

P tn 1 xn 1 | tn

xn, tn 1 xn 1,..., t1 x1 P tn 1

xn 1

| tn xn .

Определение. Функция

P t,x,s,y P s y| t x

называется

функцией

переходных вероятностей марковского процесса.

Данная функция играет важную роль при изучении подобных процессов.

Стационарные процессы

Рассмотрим случайный процесс t,t T .

Определение. Случайный процесс называется стационарным в узком смысле, если при всех n для любых t1 t2 ... tn T для всех h таких, что t1 h, … , tn h T , совпадают конечномерные распределения

Ft1,...,tn (x1,...,xn) Ft1 h,...,tn h(x1,...,xn).

Определение. Случайный процесс называется стационарным в широком смысле, если:

1)

существует a такое, что

m t a для всех t T ,

2)

существует b такое, что

D t b для всех t T ,

3)

для любых s,t T , для всех

h таких, что t h,s h T , выполняется

K s,t K s h,t h , т.е. K s,t

зависит только от разности t и s.

Проверьте, всегда ли из стационарности в узком смысле следует стационарность в широком смысле.

Для стационарного в широком смысле процесса корреляционная функция

R s,t

зависит только

от разности t s

и может быть

обозначена

R s,t R t s R . Это верно и для ковариационной функции

 

 

 

K s,t K t s K

 

Пример. Рассмотрим

случайный процесс

t

1 cos wt 2 sin wt , где

w

- константа; 1

и 2 - случайные

величины, заданные

на одном

вероятностном пространстве и не зависящие от t,

t .

 

Ответим на вопрос, когда данный случайный процесс является стационарным в широком смысле и найдем корреляционную функцию R(s,t).

m(t) E 1 cos(wt) E 2 sin(wt), т.е. для того, чтобы m(t) не зависело от t, необходимо, чтобы E 1 E 2 0.

Пусть s t . Тогда

K s,t E t s E 12 cos wt cos ws 22 sin wt sin ws

1 2 cos wt sin ws 1 2 sin wt cos ws E 12 cos wt cos ws

E 22 sin wt sin ws E 1 2 sin t s w .

Функция K s,t

не будет зависеть от t, если E

2

E 2

2

и E

2

0.

 

 

 

 

 

 

1

2

 

1

 

Из того, что E

E

2

0

следует, что D D

2

2 . Так как

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

,

 

 

E 1

E 1 2 E 2

 

E 1 2

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

D 1 D 2

 

 

2

 

 

 

 

то, следовательно, случайные величины 1 и 2 не коррелированны. Более того, мы будем считать, что они независимы.

Следовательно,

K(s,t) 2 cos wt cos ws sin wt sin ws 2 cos w t s 2 cos w .

Мы получили, что K(s,t) зависит только

от разности s t . Теперь

воспользуемся независимостью случайных величин 1

и 2 :

D t D 1 cos wt 2 sin wt D 1 cos wt D 2 sin wt

cos2 wt D 1 sin2 wt D 2

2 cos2 wt sin2 wt 2 ,

т.е. дисперсия случайного процесса D(t)

не зависит от t.

Отсюда получаем

 

 

E

 

 

 

R s,t R

 

 

t

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D t D t

 

2

cos w

cos w .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 2

 

 

 

 

 

 

Итак, если случайные величины 1

и

2

независимые с одинаковыми

математическими ожиданиями, равными 0, и одинаковыми ненулевыми

дисперсиями 2 , то t ,t -

стационарный в широком смысле случайном

процесс с корреляционной функцией R cos w .

 

 

 

Случайный процесс t ,t

представляет собой гармоническое колебание

 

 

 

 

 

arctg

1

,

со случайной амплитудой

 

2

2

и случайной фазой

 

 

 

1

2

 

 

2

 

 

 

 

 

 

распределения которых не зависят от времени t.

Линейные преобразования случайных процессов, дифференцирование и интегрирование случайных процессов

Определение. Назовем оператором отображение, ставящее в соответствие одной функции другую функцию. Примеры оператора являются:

Оператор дифференцирования

y dy A dx

Например:sinх cosх,х2 2х

A

A

Оператор интегрирования

 

х

 

 

y ydx

 

 

A 0

 

 

Например: x2

хy2dy x3

A

0

3

Определение. Оператор A называется линейным, если для любых элементов x1 и x2 пространства и любого λ из произвольного числового поля (к примеру, поля вещественных чисел) выполняются соотношения

1.A(x1+x2)=Ax1+Ax2.

2.A(λx)=λAx.

Определение. Производная случайного процесса определяется как предел lim X t t X t X ' t dX

t 0

t

dt

Так как X(t) случайный процесс, то предел нужно понимать в вероятностном смысле: это может быть, к примеру, сходимость с вероятностью единица либо сходимость в среднем порядка 2 (в среднем квадратическом).

Рассмотрим случайный процесс X(t) с математическим ожиданием mX(t) и автоковариационной функцией KX(t1,t2). Обозначим производную случайного процесса X(t) как Y(t) :

Y t dX X ' t dt

Требуется найти математическое ожидание и автоковариационную функцию производной случайного процесса, т.е. mY(t) и KY(t1,t2).

Найдем математическое ожидание левой и правой частей равенства из определения производной случайного процесса, предполагая, что соответствующие переходы возможны:

m

t EY t lim

E X t t E X t

m'

X

t

dmX t

 

 

Y

t 0

t

 

 

dt

 

 

 

 

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

Для определения KY(t1,t2) перейдем к центрированным случайным функциям

 

 

 

 

 

 

 

 

 

 

 

 

 

Y t Y t mY t и

X t X t mX t ,

т.е. отнимем от рассматриваемых

случайных величин их математические ожидания. Тогда получим:

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

 

t d X t .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

По определению автокорреляционной функции:

 

 

K

 

t ,t

 

E

 

 

 

 

 

.

 

Y

2

Y t

Y t

2

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим вместо Y t1 и Y t2 их выражения:

KY t1,t2 E d X t1 d X t2 .

1 dt2dt

Можно представить выражение под знаком математического ожидания в виде второй смешанной частной производной:

 

 

 

 

d X t1 d X t2 2 X t1 X t2 .

dt1

dt2

t1 t2

Аналогично тому как выше показано, что математическое ожидание производной случайной функции равно производной от математического ожидания, т. е. знаки дифференцирования и математического ожидания можно менять местами, можно показать, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

X

t1 X t2

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

K

 

t ,t

 

E

 

 

 

 

 

E

X t

X t

 

 

 

 

 

K

 

t ,t

 

.

Y

2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

X

2

 

1

 

 

 

 

t t

 

 

 

 

t t

 

 

 

1

 

 

 

t t

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

1

 

2

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

t ,t

 

 

2

 

K

 

t ,t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

1

2

 

2

 

X

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

кматематическому ожиданию исходного случайного процесса; для нахождения автоковариационной функции тот же линейный оператор применяется дважды

кавтоковариационной функции исходного случайного процесса.

Аналогично можно показать, что математическое ожидание интеграла от случайного процесса равно интегралу от ее математического ожидания. Т.е. операцию интегрирования можно менять местами с нахождением математического ожидания случайного процесса.

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

Примеры. 1. Уровень радиоактивности зависит от типа изотопа и количества вещества:

X t Ae αt ,

где α – const, случайная величина A - распределена по нормальному закону N a, 2 . Найдем математическое ожидание, автокорреляционную функцию,

дисперсию и нормированную автоковариационную функцию производной случайного процесса.

Легко можно найти, что:

E X t ae t Kx t1,t2 2e t1 t2

Rx t1,t2 1

Найдем математическое ожидание производной случайного процесса, используя выведенную выше формулу

EX t ae t ' a e t

Найдем автоковариационную функцию производной случайного процесса по выведенной выше формуле:

 

 

 

 

 

 

2Kx t ,t

2

 

 

 

 

 

 

2

 

t t

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1

 

t ,t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

Kx

 

 

t t

 

 

 

 

 

 

t t

 

 

 

 

 

1

 

2

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2

e

 

t1

t2

 

 

2 2

e

t1 t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для нахождения дисперсии случайного процесса x t подставим t1 t2 t в

выражение для автоковариационной функции. Получим

2 x t 2 2e 2 t

Нормированная автоковариационная функция равна:

R

, t ,t

 

 

Kx t1,t2

 

2 2e t1 t2

1

 

 

e t1 e t2

x

1

2

t1 t2

 

 

Можно провести проверку. К примеру, непосредственное нахождение характеристик случайного процесса после нахождения производной исходного случайного процесса:

X t A e t

EX t E A e t e t E A ae t

2.Пусть случайный процесс представляет собой гармоничное колебание со случайной начальной фазой

X t Acos wt