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

Nugmanov

.pdf
Скачиваний:
10
Добавлен:
10.02.2015
Размер:
1.6 Mб
Скачать

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

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

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

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

1. Ковариационная и корреляционная функции есть четные функции,

т.е.

Bξ (t1, t 2 ) Bξ (t 2 t1) Bξ (t1 t 2 ) Bξ (τ) Bξ ( τ).

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

2.Для стационарного процесса имеем (Рис. 1.8)

Bξ (τ) B(τ) mξ2 .

3.B(0) Dξ(t) .

Для доказательства этого свойства рассмотрим ковариационную функцию B(τ) при τ 0

lim B(τ) lim M(((ξ(t) mξ (t)) (ξ(t τ) mξ (t τ)))

τ 0

τ 0

M(ξ(t) mξ (t))2 Dξ(t) .

13

 

B ( )

 

4. Для любого значения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеем Bξ (0)

Bξ (τ)

.

 

 

 

 

 

 

 

5.

 

 

Нормированная

 

 

 

 

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

 

функции

 

 

 

 

B(τ)

 

 

относительно

 

 

 

 

B(0) Dξ(t) :

 

 

 

 

 

 

 

B(τ) B(0) ,

 

 

 

 

 

2

 

R ξ (τ)

 

Рис. 1.9

 

называется

 

коэффициентом

 

 

 

 

 

 

 

 

 

 

 

корреляции

 

 

 

случайного

процесса. Из пятого свойства следует, что

R ξ (τ) 1.

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

корреляции. Реализации, отстоящие по времени больше, чем интервал корреляции, считаются некоррелированными. Интервал корреляции может быть определен несколькими способами.

1) Площадь под кривой нормированной корреляционной функции

R ξ (τ) равна площади прямоугольника с основанием, равным 2 τ , и

высотой, равной единице (Рис. 1.9)

 

 

 

 

 

Bξ (τ) dτ

 

τ

 

(1.12)

 

2 Bξ (0)

 

 

 

14

2)Введем R* (τ) - огибающую нормированной корреляционной

R( )

2 ρ

- 0

Рис.1.10

 

функции и

ρ - допустимый

R*( )

коэффициент

корреляции,

 

0 ρ 1,

выбранный

из

практических соображений.

 

Будем

считать,

(Рис.

1.10),

0

значения

процесса

ξ(t)

 

 

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

на интервале

 

времени

τ , где τ выбрано из

решения равенства

 

R*

 

 

 

 

 

 

 

0

)

ρ :

τ

τ

0

.

Существуют и другие методы определения интервала корреляции.

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

τ t 2 t1 k t . Установим связь между значениями коэффициента

корреляции R ξ (τ) B(τ) B(0)

случайного процесса ξ(t) и

коэффициентом корреляции Ri j

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

ξ1, ξ1, , ξn , распределенных по одному и тому же закону с одинаковыми математическими ожиданиями m и одинаковыми дисперсиями σ2 .

Коэффициент корреляции Ri j является элементом корреляционной матрицы

15

 

1

R12

R1n

 

 

 

 

 

 

R 2n

 

 

R

R 21

1

 

(1.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R n2

 

1

 

 

 

R n1

 

 

и вычисляется как

Ri j

M (ξi m) (ξ j m)

.

(1.14)

 

 

σ2

 

 

 

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

ξ(t)

в дискретные моменты

времени ti и пусть i 0,1, , n 1. Так как ковариационная функция процесса ξ(t) рассматривается в дискретные моменты времени и на конечном интервале времени, равном n t , коэффициент корреляции случайного процесса ξ(t) удобно представить через корреляционную матрицу

R( t)

 

 

1

 

R( t)

 

R(2 t)

R((n 1) t)

 

 

 

R( t)

 

 

 

R( t)

R((n

 

 

 

 

 

 

1

 

2) t)

(1.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R((n 2) t)

 

R((n 3) t)

1

 

 

 

 

R((n 1) t)

 

 

 

 

где R(k t)

M[(ξ(t2 ) m) (ξ(t1) m)]

,

t2 t1 k .

 

 

 

 

 

σ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравнивая элементы Ri j и

 

R(k t)

корреляционных матриц R и

R( t), видно,

что R i j R(

 

i j

 

t) ,

т.е. значения коэффициента

 

 

 

 

 

 

 

 

 

 

 

 

 

корреляции случайного процесса ξ(t)

в дискретные моменты времени

совпадают со значениями коэффициента корреляции системы случайных величин ξ1, ξ1, , ξn , если τ i j t .

16

1.4. Характеристическая функция случайного процесса

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

вероятности, но и с помощью характеристической функции.

Рассмотрим случайную величину ξ , имеющую плотность

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

wξ (x) для

непрерывной

случайной

величины и распределение

вероятности

P(ξ xk ) pk ,

k 1, 2,

для дискретной случайной величины. Характеристической функцией

θξ (v) случайной

 

величины e jv ξ

называется математическое

ожидание случайной величины e jv ξ :

 

 

 

 

 

 

 

 

 

 

 

 

 

w ξ (x) e jv x dx для непрерывных сл. в.

θ (v) M(e

jv ξ

 

 

 

 

 

)

 

 

 

 

 

ξ

 

k

jv xk

 

 

 

 

 

 

p e

 

для дискретных сл. в.,

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.16)

 

 

 

 

 

 

 

где v - произвольный вещественный параметр, j 1 .

Обратное преобразование имеет вид

 

1

 

w ξ (x)

θξ (v) e j v x dv .

2 π

 

 

 

 

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

17

1.

 

 

 

 

 

 

 

θξ (v)

1.

 

 

 

 

Действительно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θξ (v)

 

 

wξ (x) e jv x dx

 

 

 

wξ (x) e jv x dx

 

 

wξ (x) dx 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

θξ (0) 1, следует из определения характеристической функции.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

θξ (v)

w ξ (x) cos (v x) dx j w ξ (x) sin (v x) dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θξ ( v)

 

 

 

 

 

 

4.

θ

ξ (v) следует из пункта 3,

 

где

θξ (v) - функция ,

комплексно сопряжённая функции θξ (v) .

 

 

 

 

 

 

5.

lim θξ (v) 0. Доказательство этого свойства приведено в [5]

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.

 

 

Если

 

плотность

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

 

 

вероятности

w ξ (x)

симметричная функция, то характеристической

функции

θξ (v) -

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(v) w (x) cos v x dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

Если характеристическая функция

θξ (v)

 

- действительная

четная функция, то плотность распределения вероятности w ξ (x) -

симметричная функция.

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

Продифференцируем θξ (v) :

dk θξ (v)

 

 

 

 

 

jk

x k w ξ (x) e jv x dx .

(1.17)

dv

k

 

 

 

 

 

 

 

 

18

Из этого выражения получим моменты k-го порядка:

 

k

 

1

dkθ

ξ

(v)

 

M ξ

 

 

 

 

 

mξ k .

 

j

k

dv

k

 

 

 

 

 

 

v 0

 

Если существуют моменты любого порядка, то разложение характеристической функции в ряд Тейлора в точке v 0 имеет вид

 

 

 

 

m

ξ k

 

 

θ

 

(v) 1

k 1

 

( j v)k

(1.18)

ξ

 

 

 

 

k!

 

 

 

 

 

На практике часто приходится сталкиваться с поиском функции и плотности распределения вероятности суммы независимых случайных величин ξi с известными характеристическими функциями θξi (v) . Согласно определению (1.16) характеристическая функция суммы независимых случайных величин равна произведению характеристических функций, т.е.

n

θ (v) M(e jv ξ ) M(e jv in 1ξi ) M(e

ξ

i 1

n

jv ξ i ) θξi (v) . (1.19)

i 1

Формула (1.15) позволяет вычислить плотность распределения вероятности суммы независимых случайных величин.

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

ψ(v) ln θξ (v) .

(1.20)

19

Кумулянтная функция может быть разложена ряд по степеням v в

точке v0 0 . Для этого рассмотрим разложение функции ln(1 z) в

ряд Тейлора в точке z0 0 :

 

 

 

 

 

 

 

 

 

 

 

 

ln(1 z) z

z2

 

z3

 

z4

(1)k 1

zk

. (1.21)

 

 

 

 

 

 

2

 

3

 

 

4

 

 

k 1

 

 

k

 

Сделаем подстановку z θξ (v) 1 и получим, используя (1.18) и

(1.21),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mξ k

 

 

k

 

 

1

 

 

mξ k

 

2

 

 

 

ψξ (v)

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

( j v)

 

 

 

 

 

( j v)

 

 

 

k 1

k!

 

 

 

 

 

2

 

k 1

k!

 

 

 

 

 

 

1

 

 

mξ k

 

3

 

1

 

 

mξ k

 

4

 

 

 

 

 

 

 

 

( j v)

k

 

 

 

 

 

( j v)

k

 

3

 

k!

 

4

 

k!

 

 

 

k 1

 

 

 

 

k 1

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

χ k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(v)

 

 

( j v)k ,

 

 

 

 

 

 

(1.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1 k!

 

 

 

 

 

 

 

 

 

где

коэффициенты разложения

 

χ k ,

называемые

кумулянтами

(семиинвариантами),

 

 

являются

 

полиномами

от

моментов

mξ k , k 1, 2, 3,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o

 

 

 

 

 

 

 

χ

1

m

ξ1

, χ

2

m

ξ2

m2

 

M(ξ)2 Dξ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ1

 

 

 

 

 

 

 

m

 

3 m

 

 

 

 

 

 

 

 

 

 

 

o

 

 

 

 

o

o

χ

3

ξ3

ξ1

m

ξ2

2 m3 M(ξ)3

,

 

χ

4

M(ξ)4 3(M(ξ)2 )2 (1.23)

 

 

 

 

 

 

 

 

 

 

ξ1

 

 

 

 

 

 

 

 

20

 

Пример 1. Пусть случайная величина ξ распределена

по

нормальному закону с математическим ожиданием mξ

и дисперсией

σ 2

, N m

 

, 2 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

(x mξ )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w(x)

 

 

 

 

 

 

 

exp

 

 

,

 

 

x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 σ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 π σ

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

Характеристическая функция случайной величины ξ имеет вид

 

 

 

 

 

 

 

 

 

 

 

jv mξ

1

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θξ (v) e

 

 

σξ v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а кумулянтная функция –

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

(v) j v m

 

 

 

1

σ2 v2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2. Рассмотрим сумму трех независимых нормально

распределенных случайных величин ξ ξ1

ξ2

ξ3 , математические

ожидания и дисперсии которых соответственно равны mξ1, mξ2 ,

mξ3

, σξ21, σξ22 ,

σξ23 .

Определим плотность распределения вероятности

w ξ (x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Характеристические функции случайных величин ξi

равны

соответственно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jv mξi

 

σξi v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θξ

 

(v) e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По формуле (1.19) вычислим характеристическую функцию

случайной величины ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

2

 

 

 

 

 

 

 

v2

 

2

θξ (v) exp i3 1

( j v mξi

 

 

 

σ

ξi

v

 

)

exp j v i3 1 mξi

 

i3 1

σξi .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

21

Введем обозначения m 3 m , σ2 3 σ2 .

ξ i 1 ξi ξ i 1 ξi

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

равна

 

jv mξ

1

2

2

θξ (v) e

 

σξ v

 

2

.

 

 

 

 

 

Но полученная формула представляет характеристическую функцию нормально распределенной случайной величины с параметрами mξ , σξ2 .

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

θξ (v1, , vn ) M e

n

j vi ξi

i 1

 

 

n

 

 

j vi xi

 

 

 

 

 

 

w ξ (x1, , x n ) e i 1

dx1 dx n

(1.24)

 

 

 

 

 

Соответственно можно записать совместную плотность распределения вероятности

w ξ (x1

, , x n )

1

 

2 π n

 

 

θξ (v1, , vn ) e

n

j vi xi

i 1 dv1 dvn (1.25)

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

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

 

 

 

θξ (v, t) M(e jv ξ(t) )

w ξ (x, t) e jv x dx

(1.26)

22

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