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

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

.pdf
Скачиваний:
17
Добавлен:
10.02.2015
Размер:
993.95 Кб
Скачать

η(t)

g(v,t) ξ(v) dv.

 

ξ

Например, в качестве весовой функции в интеграле Дюамеля имеем импульсную характеристику g(v,t) h(t v).

Математическое ожидание и корреляционная функция процесса η(t) будет иметь вид:

 

 

 

 

ξ

 

 

ξ

(v) dv ,

(1.34)

 

M(η(t))

 

g(v,t) M(ξ(v)) dv

g(v,t) mξ

 

η(t1,t

2)

 

ξ

ξ

,t1) g(v2

,t2) Bξ (v1,v2) dv1 dv2

(1.35)

B

 

 

g(v1

1.6. Эргодические случайные процессы

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

ется одна какая-то реализация x(k)(t) случайного процесса ξ(t). Для изучения физического процесса возникает необходимость вычисления плотности распределения вероятности какого либо параметра процесса, функции распределения и моментных функций процесса по одной реализации на интервале наблюдения (0, Tн). Вероятностные характеристики случайного процесса, полученные по одной реализации процесса за ограниченное время наблюдения Tн, будут случайными. Следовательно, необходим критерий, по которому можно было бы отождествить вероятностные характеристики случайного процесса, полученные по одной реализации, с характеристиками случайного процесса, вычисленных усреднением по ансамблю. Для выбора критерия предварительно рассмотрим сходимость по вероятности и сходимость в среднеквадратическом.

Положим, имеется последовательность случайных величин ξ1, ,ξn и неслучайная величина m. Последовательность случайных величин ξ1, ,ξn сходится по вероятности к величине m, если для любого ε 0 выполняется соотношение

 

lim P(

ξn m

ε) 0,

(1.36)

или

n

 

 

ε) 1.

 

lim P(

 

ξn m

 

 

 

n

 

 

 

 

 

 

23

 

Последовательность случайных величин ξ1, ,ξn сходится в среднеквадратическом к величине m, если выполняется соотношение

lim M(ξn m)2

0.

(1.37)

n

 

 

Если m Mξ mξ, то сходимость в среднеквадратическом озна-

чает стремление дисперсии случайной величины ξnк нулю, то есть

lim Dξn 0.

(1.38)

n

 

Из сходимости в среднеквадратическом следует сходимость по вероятности. Действительно, согласно неравенству Чебышева

P ξn mξ ε Dεξ2n для любого ε 0.

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

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

Таблица 4.1.а

 

Средние по времени

 

 

Средние по ансамблю

ξ(k)(t)

 

1

Tн ξ(k)(t)dt

 

 

M(ξ(t))

x w(x)dx

 

Tн

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

2

 

 

1

 

T

2

 

 

 

 

(k)(t)

 

 

 

0 н (k)

(t) dt

 

 

M(ξ(t))2

x2 w(x)dx

 

 

 

 

 

 

 

 

 

Tн

 

 

 

 

 

ξ(k)(t)ξ(k)(t τ)

 

 

Bξ (τ)

 

 

1

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

0н ξ(k)(t)ξ(k)(t τ)dt

 

 

x1 x2 w(x1,x2,τ)dx1dx2

T

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В таблицах 4.1.а и 4.1.б приведены некоторые числовые характеристики непрерывного и дискретного по времени случайного процес-

24

са, вычисленные по времени и по ансамблю. Из приведенных таблиц видно, что при вычислении среднего по времени используется произвольная k-ая реализация случайного процесса и средние по времени

 

 

 

Средние по времени

 

 

Средние по ансамблю

 

 

 

 

 

 

 

1

 

 

N

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

(k)

 

 

 

x

(k)

 

 

 

M(ξ(t)) x[i]p(ξ(t) x[i])

 

ξ

 

 

(t)

N

 

 

[i]

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

1 N

 

 

 

 

2

M(ξ(t))2

 

 

 

 

 

 

 

 

(k)

 

 

 

 

 

 

(k)

 

 

N

 

 

 

 

 

 

 

 

ξ

 

 

(t)

 

 

 

 

x

 

 

[i]

 

 

(x[i])

2

p(ξ(t)

 

x[i])

 

 

 

 

 

 

 

 

 

N i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

ξ(k)(j)ξ(k)(j m)

 

 

 

 

 

N 1N 1

 

 

 

 

 

 

 

Bξ (m)

 

 

x[i]x[j]

 

 

 

 

 

 

N m

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

i 0 j 0,

i j

m

 

 

 

 

 

 

x(k)[i]x(k)[i m],

p(x[i], x[j]),

m 0,1, ,N 1

N m

 

i 1

 

 

 

 

 

 

 

 

j 1, , N m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, где означает усредняемую величину, не зависят от времени t.

Таблица 4.1.б

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

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

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

lim D 0.

(1.39)

Tн

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

Среднее по времени для непрерывного случайного процесса имеет вид (Таблица 4.1.а)

25

ξ(k)(t)

1

Tнξ(k)(t)dt.

(1.40)

T

 

0

 

 

н

 

 

Определим математическое ожидание и дисперсию среднего по времени:

 

 

1

T

 

 

 

M ξ(k) (t)

 

0

н

(k)(t)dt mξ ,

(1.41)

T

 

 

н

 

 

 

 

D ξ(k)(t)

M ξ(k) (t)

mξ 2 .

(1.42)

Формула (1.42) отражает среднеквадратическое отклонение среднего по времени от среднего по ансамблю. Если с увеличением вре-

мени наблюдения Tн дисперсия среднего по времени Dξ(k)(t)

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

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

(1.39)

lim Dξ(k)(t) 0

Tн

Преобразуем выражение (1.42)

Dξ(k)(t) M ξ(k)(t) mξ 2 M ξ(k)(t1)ξ(k)(t2) m2ξ

lim

1

Tн

Tн M(ξ

(k)(t

 

) ξ(k)(t

 

 

))dt

 

dt

 

m2

 

 

0

 

 

 

 

 

 

Tн T2

0

 

 

 

 

1

 

2

 

1

 

2

ξ

 

 

н

 

 

 

1

 

Tн

Tн Bξ (t2

 

 

 

 

 

 

 

 

 

 

lim

 

t1))dt1 dt2 mξ2

 

 

 

2

 

 

 

Tн T

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

Tн

Tн (Bξ (t2

 

 

 

 

 

 

 

 

 

lim

 

1

 

t1) mξ2)dt1 dt2

 

 

 

 

2

 

 

 

Tн T

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

Tн

Tн B(t

 

 

 

 

 

 

 

 

 

 

lim

 

1

 

2

t1)dt1 dt2.

 

 

 

 

2

 

 

 

 

Tн T

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

Согласно критерию эргодичности должно соблюдаться

26

D ξ(k)(t) lim

1

Tн

Tн B(t2

t1)dt1 dt2

0. (1.43)

2

Tн T

0

0

 

 

 

н

 

 

 

 

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

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

τ t2 t1,

t1 t0 τ/2 ,

0 t1

Tн,

t0 (t2 t1)/2,

t2 t0 τ/2.

 

(1.44)

0 t2 Tн.

Определим область интегрирования для переменных τ и t0 (Рис.

1.11). Из условия (1.44) имеем:

1.

0 t0 τ/2 Tн.

(1.45)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно левой части неравенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

2.2

 

 

 

 

 

 

(1.45)

на плоскости ( , t0 ) имеем

 

 

 

 

 

1.2

 

 

 

 

 

 

 

 

 

 

прямую 1.1. Соответственно для пра-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вой части неравенства (1.45) на плос-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

кости ( ,

t0 ) имеем прямую 1.2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

2. Точно также согласно условию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.44) должно выполняться

 

 

-T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.1

 

 

 

 

 

 

 

 

 

0 t0 τ/2 Tн.

(1.46)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Левая часть неравенства (1.46) со-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1.11

ответствует прямой 2.1 на плоскости

 

 

 

 

 

 

 

(τ, t0

).

Правая часть

неравенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.46) соответствует прямой 2.2 на плоскости (τ, t0 ). Область, ограниченная прямыми 1.1, 1.2, 2.1 ,2.2, будет областью интегрирования для новых переменных τ, t0 . Якобиан преобразования

 

 

t1

 

t1

 

 

J

 

t0

 

 

 

1

 

t

2

 

t2

 

 

 

t

0

 

 

 

 

Вычислим двойной интеграл (1.43)

27

1

0

 

T / 2

 

 

T τ / 2

 

 

 

 

 

 

 

1

 

T

 

 

 

 

 

T / 2

 

T τ / 2

 

 

 

 

 

н

 

 

 

н

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

н

 

н

 

 

 

 

B(τ)

 

dt0

 

dt0

 

 

 

 

B(τ)

dt0

dt0

 

T2

T

2

 

 

н Tн

 

τ / 2

 

 

Tн / 2

 

 

 

 

 

 

 

н 0

 

 

 

 

 

τ / 2

 

Tн / 2

 

 

 

 

 

1

0

 

 

 

(τ) T τ dτ

1

Tн

 

(τ) T

τ dτ

 

 

 

 

 

 

 

 

B

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T2

 

 

 

н

 

 

 

 

 

 

 

 

T

2

 

 

 

н

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Tн

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

1

Tн

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B(τ) dτ.

 

(1.47)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

B(τ) 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н T

 

 

 

 

н

 

 

 

 

 

 

 

 

н T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

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

lim

D ξk (t)

0 дос-

таточно выполнения условия

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Tн

 

 

 

 

 

 

 

 

Tн

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

(τ) dτ 0,

 

 

(1.48)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

н T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

которое называется условием эргодичности Слуцкого.

Если условие Слуцкого выполняется, то процесс ξ(t) можно считать эргодическим по отношению к среднему по времени и использо-

вать среднее по времени ξ(k)(t) для оценки математического ожи-

дания M(ξ(t)). Время наблюдения необходимо взять достаточно большим, но не настолько большим, чтобы нарушились условия, обеспечивающие стационарность исследуемого процесса.

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

η(τ) ξ(k)(t)ξ(k)(t τ), то для применения условия Слуцкого необ-

ходимо вычислить соответствующую ковариационную функцию и проверить выполнение условия (1.39). Для процесса η(τ) необходимо вычислить B(ν). При выполнении условия

lim

1

TнB

 

(ν) dν 0

 

T T

 

 

н

н T

 

 

 

 

н

 

 

можно говорить, что процесс ξ(t) - эргодический по отношению к ковариационной функции процесса η(τ).

28

Пример 4.1. Пусть процесс ξ(t) - стационарный случайный процесс с ковариационной функцией B(τ) m сos ω0τ. Доказать, что

процесс ξ(t) - эргодический процесс.

 

 

 

 

 

 

 

 

 

 

Ковариационная функция B(τ) B(t2 t1) при

 

t2 t1

 

 

 

 

не стремится к нулю. Поэтому используем условие Слуцкого

 

 

1

Tн

sinω

0

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B(τ) dτ 2m

 

 

 

н

 

0.

 

 

 

 

 

T

ω

T

 

 

 

 

 

 

 

 

н

Tн

 

0 н

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

 

 

Условие Слуцкого выполняется. Следовательно, процесс ξ(t) - эргодический по отношению к среднему по времени.

Пример 4.2. Найти условие эргодичности по отношению к среднему по времени для случайного процесса η(t) k ξ(t) Y b, где k

– постоянная величина, не равная нулю, ξ(t) - стационарный эргоди-

ческий случайный процесс с Mξ(t) mξ, Dξ(t) σξ2, Y – случайная

величина с MY mY,

DY σY2 , процесс ξ(t)

и величина Y - неза-

висимы, b – постоянная величина.

 

 

о

 

о

о

1.Вычислим η(t) k (ξ(t) mξ ) (Y mY) k ξ(t) Y,

оо

2.Вычислим B(t1,t2) M(η(t1) η(t2))

k2 B(t1,t2) σ2Y k2 B(t2 t1) σ2Y B(t2 t1).

3. Из

полученной формулы для

ковариационной функции

B(t1,t2)

видно, что ковариационная

функция с увеличением раз-

ности t2 t1 стремится к постоянной величине σ2Y, а не к нулю, как

того требует условие эргодичности. Более того, если произвести проверку эргодичности процесса η(t) по условию Слуцкого, то увидим, что оно нарушается.

Следовательно, для того, чтобы процесс η(t) был эргодическим необходимо отсутствие случайной величины Y. Тогда в выражении для B(t1,t2) будет отсутствовать величина σ2Y и процесс η(t)бу-

дет эргодическим. Постоянные величины k и b не влияют на условие эргодичности исследуемого процесса.

29

Пример 4.3. Рассмотрим, как оценивается одномерная плотность распределения вероятности wξ (x) случайного эргодического про-

цесса ξ(t) по одной единственной k-ой реализации ξ(k)(t), (Рис. 1.12).

Преобразуем процесс ξ(t) в новый процесс η(t), который принимает значения y0 0, если процесс ξ(t) не принадлежит интервалу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(vi, vi v), и

y1 1,

 

 

 

(k)(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

если процесс ξ(t) на-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ходится

в

интервале

vi+ v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(vi, vi v).

В ре-

vi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зультате

 

образуется

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

последовательность

(t)

 

 

 

 

 

 

 

 

 

 

единичных импульсов

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(k)(t),

появляющих-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ся в случайные момен-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ты времени tig и со

 

 

 

 

 

tig

 

tig+ ij

 

 

 

 

 

 

 

 

случайной

длительно-

 

 

 

 

Рис.1.12

 

стью τij,

(индекс i

 

 

 

 

 

при моментах времени

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tig и τij означает, что вычисляется распределение плотности вероят-

ности для v vi):

 

 

 

если v

 

ξ

(k)

(t) v

 

v

 

 

 

1,

i

 

 

i

 

y(k)(t)

 

 

 

 

(k)

 

 

(k)

 

 

 

 

 

 

если vi ξ

(t) или ξ

(t) vi v.

 

 

0,

 

 

 

Вычислим среднее по времени процесса η(t):

 

 

 

T

 

 

 

 

 

 

 

tig τi j

 

 

 

 

 

 

 

1

н

 

 

1

 

 

 

 

 

 

1

τi j .

η(k)(t)

y(k)(t)dt

y(k)(t)dt

T

T

T

 

н

0

 

 

н

 

j

 

tig

 

 

 

 

н

j

Среднее по ансамблю процесса η(t) равно

 

 

 

 

Mη(t) 1P(η(t) 1) 0 P(η(t) 0) P(vi

ξ(t) vi v)

 

 

vi v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(v) dv w(vi) v.

 

 

 

 

 

 

 

 

vi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

Докажем, что процесс η(t) эргодический. Рассмотрим ковационнуюриационнуюфункциюфункциюпроцессапроцессаη(t).η(t).

B(t1,t2) M(η(t1),η(t2)) Mη(t1)Mη(t2).

Корреляционная функция процесса η(t)имеет вид

1 1

Bη(t1,t2) M(η(t1),η(t2)) yi yj P(η(t1) yi,η(t2) yj)=

i0j 0

vv v v

P(η(t1) 1,η(t2) 1)

 

wξ (v1,v2,t2 t1)dv1 dv2.

 

 

 

v

v

Подставим значения Bη(t1,t2) и Mη(t) в B(t1,t2)

 

v v v v

 

 

B(t1,t2)

 

(wξ(v1,v2,t2 t1) wξ(v1,t1)wξ (v2,t2))dv1 dv2

vv

B(t2 t1).

Как видно из предыдущего выражения, поведение ковариационной функции зависит от поведения совместной плотности распреде-

ления wξ (v1,v2,t2 t1) при

 

t2 t1

. Если статистическая связь

между

 

значениями

v1 и v2случайного

процесса ξ(t) исчезает при

 

t2 t1

 

, то

 

 

 

 

 

 

 

 

 

 

 

 

 

wξ (v1,t1)wξ (v2,t2),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wξ (v1,v2,t2 t1)

 

 

t

2

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

и B(t2 t1) 0 при

 

t2 t1

 

 

, что является достаточным усло-

 

 

вием эргодичности процесса η(t).

 

 

Учитывая, что процесс η(t) - эргодический, при достаточно

 

большом времени наблюдения Tн можно записать

 

 

 

 

1

τij w(vj) v

или

 

 

 

 

T

 

 

 

 

н

j

 

 

 

 

 

 

 

 

 

 

 

 

w(vi)

1

 

 

 

τij .

(1.49)

 

 

 

v T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

н

 

 

j

 

Разделив возможные значения v случайного процесса ξ(t) на интервалы с шагом v, определим величину плотности распределения вероятности для каждого i-го интервала при помощи (1.49).

31

1.7. Спектральная функция стационарного случайного процесса

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

Однако применить непосредственно преобразование Фурье к анализу случайных процессов невозможно, т.к. случайному процессу

ξ(t) соответствует множество возможных реализаций x(t) ξ(k)(t) и для каждого из них можно применить преобразование Фурье. В результате нельзя однозначно отождествить обратное преобразование Фурье со случайным процессом ξ(t). Поэтому для описания процесса ξ(t) пользуются усредненной спектральной характеристикой.

Рассмотрим стационарный случайный процесс ξ(t). Для каждой реализации ξ(k)(t) на интервале времени ( Tн /2, Tн /2) найдем преобразование Фурье и рассчитаем квадрат модуля преобразования Фурье FξTн (jω), приходящуюся на единицу полосы частот. Опреде-

лим математическое ожидание квадрата модуля преобразования Фу-

рье FξT (jω). Предел этой величины при

Tн называют спек-

н

 

 

 

 

 

 

 

 

 

 

тральной плотностью Fξ (ω) случайного процесса ξ(t).

 

Покажем связь спектральной плотности

Fξ (ω) случайного про-

цесса (t) с корреляционной функцией B ( ).

 

Пусть известна реализация ξ(k)(t) на интервале ( T

/2, T /2).

 

 

 

 

 

 

 

 

 

н

н

Преобразование Фурье этой реализации имеет вид

 

FξT

(jω)

Tн / 2

ξ(k)(t) exp( jω t)dt.

 

 

 

 

н

 

Tн / 2

 

 

 

 

 

 

 

 

 

 

Рассмотрим квадрат модуля преобразования Фурье FξT

(jω) дан-

 

 

 

 

 

 

 

 

 

н

ной реализации, приходящуюся на единицу полосы частот 1/Tн

F

(ω)

1

 

 

F

(jω)

 

2.

 

 

 

 

 

 

 

 

 

 

ξT

 

Tн

 

ξT

 

 

 

 

н

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

32