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

Tom_2

.pdf
Скачиваний:
59
Добавлен:
18.02.2016
Размер:
3.2 Mб
Скачать

Закон распределения случайной функции X (t) , наряду с

плотностью, можно задать n-мерной функцией распределения.

Например, двумерной функцией распределения случайной функции X (t)

называется неслучайная функция F (x1, x2 ,t1,t2 ) = P( X (t1) < x1, X (t2 ) < x2 ) , которая при фиксированных значениях t1, t2 дает функцию совместного распределения соответствующих сечений X (t1) и X (t2 ) .

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

Действительно, при фиксированном t получим случайную величину X (t) . Математическое ожидание этой величины будет

mX (t) = M [ X (t)] .

Когда t – переменная величина, то mX (t) – также переменная

величина. Таким образом, математическим ожиданием случайной функции X (t) называется такая неслучайная функция, значение которой

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

M éX (t)ù =

 

(x,t)dx = m

 

(t) , t T ,

 

ò

xf

X

(4)

ë

û

1

 

 

 

−∞

здесь f1(x,t) – одномерная плотность распределения.

Из (4) следует, что математическое ожидание случайной функции X (t) есть какая-то средняя функция mX (t) , около которой

разными способами размещаются конкретные реализации случайной функции (см. рис. 1).

Дисперсией случайной функции X (t) называется неслучайная функция DX (t) , значение которой при каждом значении аргумента t

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

DX (t) = D[X (t)] = ò (x - mX (t))2 f1(x,t)dx , t T ,

−∞

404

где f1(x,t) – одномерная плотность распределения случайной функции X (t) , mX (t) – ее математическое ожидание. Дисперсия

случайной функции характеризует степень разброса возможных реализаций случайной функции относительно средней функции mX (t) . Среднее квадратичное отклонение случайной функции равно

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

σ X (t) = DX (t) , t T .

Пример 2. Найти математическое ожидание и дисперсию

случайной

функции X (t) = ξ sinωt , ω число, ξ – случайная

величина:

Mξ = a, Dξ = σ 2 .

Решение. Используя свойства математического ожидания и дисперсии, находим:

mX (t) = MX (t) = M (ξ sinωt) = sinωt × Mξ = a sinωt , DX (t) = DX (t) = D(ξ sinωt) = sin2 ωt × Dξ = σ 2 ×sin2 ωt .

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

корреляционной функцией или функцией связи между сечениями случайной функции.

Корреляционной функцией случайной функции X (t) называется такая неслучайная функция двух аргументов KX (t1,t2 ) , значение которой при каждой паре значений t1 и t2 равно корреляционному моменту соответствующих сечений случайной функции, это значит

K

X

(t ,t

2

) = M é( X (t ) - m

X

(t ))( X (t

2

) - m

X

(t

2

))ù

(5)

или

1

 

ë

1

1

 

 

û

 

 

+∞ +∞

 

 

 

 

 

 

 

 

 

 

 

 

 

(x1 - mX (t1))(x2 - mX (t2 )) f2 (x1, x2 ,t1,t2 )dx1dx2 ,

 

KX (t1,t2 ) = ò ò

 

−∞ −∞

где f2 – двумерная плотность распределения.

Свойства корреляционной функции.

1) Функция KX (t1,t2 ) симметрична относительно своих

аргументов, т.е. KX (t1,t2 ) = KX (t2 ,t1) . Это следует из (5). 2) KX (t1,t1) = DX (t1) .

Доказательство. Действительно, при t1 = t2 формула (5) имеет

вид K

X

(t ,t ) =

= M ( X (t ) - m

X

(t ))2

= D

X

(t ) .

 

1 1

1

1

 

1

405

Это значит, что при t1 = t2 корреляционная функция случайной

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

Вместо корреляционной функции часто пользуются

нормированной корреляционной функцией

 

 

 

 

rX (t1,t2 ) =

 

KX (t1,t2 )

 

 

.

 

 

 

 

 

 

 

(6)

 

 

 

 

σ X

(t1X (t2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Величина rX (t1,t2 )

 

является коэффициентом корреляции сечений

X (t1) и X (t2 ) . При t1 = t2

из формулы (6) имеем

 

 

 

 

 

 

 

 

 

r (t ,t ) =

 

KX (t1,t1)

 

=

DX (t1)

=1.

 

 

 

 

 

 

 

 

σ X (t1X (t1)

 

 

 

 

 

 

 

 

 

X 1 1

 

 

 

 

 

 

σ

2

(t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

1

 

 

 

 

 

 

 

 

3)

 

KX (t1,t2 )

 

£

 

 

 

= σ X (t1X (t2 ) .

 

 

(7)

 

 

 

 

DX (t1)DX (t2 )

 

 

Неравенство (7) следует из (6), так как

 

rX (t1,t2 )

 

£ 1 .

 

 

 

 

 

 

4) Если к случайной функции прибавить неслучайную, то ее

корреляционная функция не изменится.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Пусть

имеем случайную

функцию

X (t) и

неслучайную

ϕ(t) .

 

Рассмотрим

 

 

сумму

 

%

 

 

 

 

 

 

 

 

X (t) = X (t) +ϕ(t) .

Математическое

ожидание

M [X (t) +ϕ(t)] = mX (t) +ϕ(t) .

Согласно

(5),

 

(t ,t

 

) = M

é(

 

 

 

 

 

 

 

(m

 

(t ) +ϕ(t )))´ =

 

K %

2

X (t ) +ϕ(t ) -

X

 

X

1

 

 

 

ë

 

 

1

 

 

1

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

´(X (t2 ) +ϕ(t2 ) - (mX (t2 ) +ϕ(t2 )))ûù =

 

 

 

 

= M é( X (t ) - m

X

(t ))( X (t

2

) - m

X

(t

2

))ù = K

X

(t ,t

2

),

 

 

 

ë

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

û

 

 

1

 

 

что и требовалось доказать.

5)Если случайную функцию X (t) умножить на неслучайную

ϕ(t) , то корреляционная функция случайной функции умножится

на ϕ(t1) ×ϕ(t2 ) .

 

а ϕ(t) –

Доказательство. Пусть X (t)

– случайная функция,

неслучайная функция. Рассмотрим

%

 

произведение X (t) = X (t) ×ϕ(t) .

Математическое ожидание M ( X (t)ϕ(t)) = ϕ(t)M ( X (t)) .

Согласно

формуле (5),

KX% (t1 ,t 2) = M[(X (t1)ϕ(t1) -ϕ(t1)MX (t1))(X (t 2)ϕ(t 2) -ϕ(t 2)MX (t 2 ))] =

=ϕ(t1)ϕ(t 2)M[(X (t1) - mX (t1)) × (X (t 2) - m X (t 2))] = ϕ(t1)ϕ(t 2)K X (t1 ,t 2),

что и требовалось доказать.

406

Пример 3. Найти корреляционную функцию случайной

функции

 

X (t) = ξ sinωt ,

где

ω –

число, t

время,

ξ –

 

случайная

величина:

Mξ = a,

Dξ = σ 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Согласно формуле (5),

 

 

 

 

 

 

 

 

 

 

K

X

(t ,t

2

) = M é(ξ sinωt - m

X

(t )) ×(ξ sinωt

2

- m

X

(t

2

))ù .

 

 

 

 

1

 

 

ë

 

 

1

1

 

 

 

 

û

Ввиду примера 2, mX (t) = a sinωt , поэтому

 

 

 

 

 

 

 

 

.

K

X

(t ,t

2

) = sinωt ×sinωt

2

× M - a)2 = σ 2

×sinωt ×sinωt

2

 

 

1

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

Если заданы две случайные функции X (t)

и Y (t) ,

то для них

вводится понятие взаимной корреляционной функции. Взаимной корреляционной функцией или корреляционной функцией связи двух случайных функций X (t) и Y (t) называется корреляционный момент

сечений этих функций при любых взятых значениях аргументов t и s:

K

XY

(t, s) = M é( X (t) - m

X

(t))(Y (s) - m

(s))ù

(8)

или

 

ë

 

 

Y

û

 

 

 

 

 

 

 

 

 

 

 

 

 

 

é

0

0

ù

 

 

 

 

 

KXY (t, s) = M ëX (t),Y (s)û ,

 

 

0

 

 

0

 

 

 

 

 

 

где X (t) = X (t) - mX (t), Y (s)

= Y (s) - mY (s) .

 

 

Формулу (8) можно записать

 

 

 

 

 

 

 

+∞ +∞

 

 

 

 

 

 

KXY (t, s) = ò

ò (x - mX (t))( y - mY (s)) f2 (x, y,t, s)dxdy ,

 

 

−∞ −∞

 

 

 

 

 

 

где f2 – совместный

закон

(плотность)

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

случайных

функций X (t) и Y (t) . Случайные функции X (t) , Y (t) называются коррелированными, когда их взаимная корреляционная функция отлична от тождественного нуля: KXY (t, s) ¹ 0 для некоторых значений t и s. В случае KXY (t, s) º 0 случайные функции X (t) , Y (t) называются

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

Нормированная взаимная корреляционная функция имеет вид:

rXY (t, s) =

KXY (t, s)

.

σ X (tY (s)

 

 

§ 2. Операции над случайными функциями

407

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

10. Сложение случайных функций. Пусть X (t) и Y (t) две случайные функции одного и того же аргумента с математическими ожиданиями mX (t) , mY (t) . Рассмотрим их сумму

 

 

 

 

 

 

 

Z(t) = X (t) + Y (t) .

 

 

 

 

 

 

(1)

 

Найдем математическое ожидание функции Z(t):

 

 

 

 

 

M[Z (t)] = M[X (t) + Y (t)] = M[X (t)] + M[Y (t)] = mX (t) + mY (t) .

(2)

 

Это означает, что при сложении случайных функций их

математические ожидания складываются.

 

 

 

 

 

 

 

 

 

Вычтем из выражения (1) выражение (2):

 

 

 

 

 

 

 

 

 

 

 

 

Z(t) - mZ (t) = X (t) - mX (t) + Y (t) - mY (t) ,

 

 

 

 

0

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. Z(t) = X (t) + Y (t) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем корреляционную функцию случайной функции Z(t):

 

K

 

(t ,t

 

 

 

0

0

 

)ù

= M

0

 

0

0

 

 

0

 

))úù =

 

Z

2

) = M éZ(t ) Z(t

2

êé(X (t )

+ Y (t ))(X (t

2

) + Y (t

2

 

 

1

 

ë

1

 

û

 

ë

1

1

 

 

 

û

 

 

 

 

 

é 0

 

0

0

 

0

 

0

0

0

 

0

ù

 

 

 

 

= M

ëX (t1 ) X (t2 ) + X (t1)Y (t2 ) +

X (t2 )Y (t1) + Y (t1)Y (t2 )û

=

 

é 0

 

0

ù

 

é 0

 

0

ù

 

é 0

0

ù

 

 

é 0

 

0

ù

= M ëX (t1) X (t2 )û + M

ëX (t1)Y (t2 )û

+ M ëX (t2 )Y (t1 )û + M

ëY (t1 )Y (t2 )û =

= KX (t1,t2 ) + KXY (t1,t2 ) + KXY (t2 ,t1 ) + KY (t1,t2 ).

Если X (t) и Y (t) – некоррелированны, то KXY º 0 и

KZ (t1,t2 ) = KX (t1,t2 ) + KY (t1,t2 ) .

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

n

Z(t) = åXi (t) .

i=1

Для математического ожидания имеем

n

M [Z(t)] = M åXi

i=1

n

(t) = åmXi (t) .

i=1

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

408

n

n

 

KZ (t1,t2 ) = åKXi (t1,t2 ) +

å KXi X j (t1,t2 ) .

(3)

i=1

i, j=1 (i¹ j)

 

Если X1(t), X2 (t),..., Xn (t) некоррелированы, то из (3)

имеем

n

KZ (t1,t2 ) = åKXi (t1,t2 ) .

i=1

20. Понятие о дифференцировании и интегрировании случайных функций. Если реализации функции X(t)

дифференцируемы, то производной случайной функции X(t) назовем случайную функцию Y(t), состоящую из производных x′(t) всех

реализаций функции X(t). Аналогично, если реализации функции X(t)

интегрируемы, то интегралом от случайной функции X(t) назовем

t

случайную функцию Y(t), состоящую из интегралов y(t) = ò x(τ )dτ от

0

всех реализаций функции X(t).

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

Если Y (t) = dX (t) , то m

(t) = M éY (t)ù =

dmX (t)

,

(4)

 

dt

Y

ë

û

dt

 

 

где mX (t) = M [X (t)] . Это означает, что математическое

ожидание

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

Корреляционная функция

 

производной Y (t) случайной

функции X (t) равна второй

смешанной производной

ее

корреляционной функции, взятой по обеим переменным:

 

K (t ,t

 

) =

2 K

X

(t ,t

2

)

.

(5)

2

 

 

1

 

Y 1

 

 

t1t2

 

 

 

 

Если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

t

 

 

 

Z(t) = ò X (τ )dτ , то

mZ (t) = òmX (τ )dτ ,

(6)

0

 

 

 

 

 

 

0

 

 

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

Корреляционная функция интеграла Z(t) от случайной функции X (t) равна двукратному интегралу от ее корреляционной ф у н к ц и и :

409

t1 t2

 

KZ (t1,t2 ) = ò ò KX 12 )dτ1dτ2 .

(7)

0 0

Формулы (4) – (7) доказываются в теории случайных функций. Часто действие изучаемой технической системы (процесса)

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

Пример 1. На вход интегратора подается случайная функция X (t) с математическим ожиданием mX (t) = 4t + 5 и корреляционной

функцией KX (t1,t2 ) = cost1 cost2 . Найти характеристики на выходе системы-интегратора.

 

 

 

 

t

 

 

 

Решение. Пусть Y (t) = ò X (τ )dτ , тогда, согласно (6),

 

 

 

 

0

 

 

 

t

t

 

 

 

mY (t ) = òmX (τ )dτ = ò(+ 5)dτ = 2t2 + 5t.

 

 

 

0

0

 

 

 

Согласно выражению (7), имеем

 

 

 

 

t1 t2

 

 

 

KY (t1,t2 ) = ò òcosτ1 cosτ2dτ1dτ2 = sin t1 ×sin t2 .

 

 

 

 

0 0

 

 

 

Чтобы получить

DY (t) , в выражении для KY (t1,t2 ) полагаем

t

= t

2

= t . Тогда D (t) = sin t ×sin t = sin2 t .

1

 

Y

 

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

Элементарной случайной функцией называется случайная функция вида X (t) = X × f (t) , где X – случайная величина, математическое

ожидание которой равно нулю, а f (t) – неслучайная функция.

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

410

Для элементарной случайной функции M ( X (t)) = 0 . Найдем ее

дисперсию

 

 

 

 

 

 

 

D éX (t)ù = M[X × f (t)]2

= f 2 (t)M (X 2 ) = f 2 (t)D

X

,

D

X

= DX.

ë

û

 

 

 

 

Это означает, что дисперсия элементарной случайной функции равна дисперсии случайной величины, умноженной на квадрат неслучайной функции. Найдем корреляционную функцию KX (t1,t2 )

элементарной случайной функции:

KX (t1,t2 ) = M [Xf (t1)Xf (t2 )] = f (t1) f (t2 )M (X 2 ) = f (t1) f (t2 )DX ,

где DX – дисперсия случайной величины X.

Производная элементарной случайной функции X (t) равна

[X (t)]¢ = [Xf (t)]¢ = Xf ¢(t) .

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

Найдем интеграл от элементарной случайной функции:

t

t

t

ò X (t1) dt1 = ò Xf (t1)dt1 = X ò f (t1) dt1 .

0

0

0

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

Каноническим разложением случайной функции X(t) называется представление ее в виде суммы математического ожидания и взаимно некоррелированных элементарных случайных функций:

n

 

 

X (t) = mX (t) + å Xiϕi (t),

(8)

i=1

X (t) , Xi

 

где mX (t) – математическое ожидание

взаимно

некоррелированные (M [Xi X j ]= 0,i ¹ j)

случайные

величины,

которые называются коэффициентами, ϕi (t) – неслучайные функции,

называемые координатными функциями.

Выражение (8) можно записать так:

n

0

n

 

X (t) - mX (t) = å Xiϕi (t)

или X (t) = å Xiϕi (t) .

(9)

i=1

 

i=1

 

Теорема 1. Если случайная функция задана каноническим разложением (8), то ее корреляционная функция имеет следующее каноническое разложение:

411

n

 

KX (t1,t2 ) = åDXi ϕi (t1i (t2 ),

(10)

i=1

 

где DXi дисперсия коэффициента Xi .

Доказательство. На основе определения корреляционной

функции имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

(t ,t

 

 

 

é

0

 

 

0

(t

 

ù

 

 

 

 

 

 

X

2

) = M êX (t

 

) X

2

)ú .

 

 

 

 

 

 

 

1

 

 

ë

1

 

 

 

û

 

 

 

Согласно (9), можно написать

 

 

 

 

 

 

 

 

 

 

K

 

(t ,t

 

) = M

é

n

X ϕ (t )

n

X

ϕ (t

 

ù

 

X

2

ê

å

å

2

)ú =

 

 

 

1

 

 

ê

i i 1

 

 

 

i i

ú

 

 

 

 

 

 

 

 

ëi=1

 

 

 

 

i=1

 

 

 

 

û

 

n

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

= åM [Xiϕi (t1)Xiϕi (t2 )] = åDXi ϕi (t1i (t2 ),

i=1

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

что и требовалось доказать.

 

 

 

 

 

 

 

 

 

 

 

 

Обратное утверждение тоже верно.

 

 

 

 

 

 

Теорема 2.

Если

 

известно

 

 

каноническое

разложение

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

корреляционной функции KX (t1,t2 ) = åDXi ϕi (t1i (t2 ) ,

то случайная

 

 

 

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

 

функция X(t) допускает каноническое разложение вида (8), где коэффициенты Xi имеют дисперсии DXi .

Примем это утверждение без доказательства.

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

Пример 2. На вход динамической системы, работа которой

описывается

оператором вида

Y(t) = 2t

dX (t)

 

+ 3t2

,

подается случайная

 

 

 

 

 

 

 

dt

 

 

 

функция,

заданная

 

каноническим

 

разложением

X (t) = t + X1 cos2t + X2 sin 2t ,

DX

= DX

2

= 2.

 

 

 

 

 

1

 

 

 

 

 

 

 

Найти каноническое разложение случайной функции Y(t) и ее характеристики на выходе системы.

Решение. Вычисляя 2t

dX (t)

+ 3t2 , получаем каноническое

dt

 

 

разложение для Y (t) :

 

Y(t) = 2t (1- 2X1 sin 2t + 2X2 cos2t) + 3t2 =

=2t + 3t2 - 4tX1 sin 2t + 4tX2 cos2t.

412

Отсюда M [Y (t)] = 2t + 3t2 . Из формулы (10) получаем

KY (t1,t2 ) = DX1 (−4t1 sin 2t1)(−4t2 sin 2t2 ) +

+DX2 (4t1 cos2t1)(4t2 cos2t2 ) = 32t1 ×t2 cos2(t2 - t1) .

При t1 = t2 = t имеем

D[Y (t)] = KY (t,t) = 32t2 .

Из примера 2 видно, что функция Y (t) , получаемая на выходе в

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

§ 3. Стационарные случайные функции

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

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

Случайная функция X (t) называется стационарной, если ее

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

t1

,

ч

т

о

о

з

н

а

ч

а

е

т

 

 

mX (t) = c = const,

KX (t1,t2 ) = KX (τ ),

где τ = t2 - t1 .

 

 

Таким образом, корреляционная функция стационарной случайной функции зависит только от одного аргумента. Найдем дисперсию стационарной случайной функции. Как известно, D[X (t)] = KX (t,t) .

Для данного случая D éX (t)ù

= K

X

(0) = D

X

.

ë

û

 

 

 

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

корреляционная

функция

стационарной

случайной функции является

ч

е

т

н

о

й

.

413

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]