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

ЕГОРОВ В.А.- Обработка траекторий случайных процессов

.pdf
Скачиваний:
83
Добавлен:
08.07.2017
Размер:
1.25 Mб
Скачать

6.2. Спектральное представление для стационарных процессов

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

Напомним, что M1 M – это множество всевозможных конечных линейных комбинаций случайных величин X (t) ; M H – замыкание M1 в метрике гильбертова пространства H. Обозначим через L1 множество все-

n it λ

возможных тригонометрических многочленов вида Q(λ) c je j , через L

j1

обозначим замыкание L1 в метрике L2 (dF ) , где F – спектральная функция из теоремы 6.3.

Определим взаимно однозначное линейное соответствие между множеством M1 и множеством L1 с помощью соотношения

n

n

 

c j X (t j )

c jeit j λ.

(6.3)

j1

j1

 

Отметим, что значению процесса X (t) в момент времени t

соответствует

функция eitλ , т. е.

 

 

X (t) eitλ.

(6.4)

Покажем, что соответствие « » является изоморфизмом, т. е. сохраняет скалярные произведения.

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

n

c j X (t j ),

j1

n

n n

 

 

d j X (t j )

c j d

j K (t j tk )

j1

j1k 1

31

n n

 

 

 

i(t j tk

 

 

 

n n

it j λ

 

 

 

it

 

λ

 

 

 

 

dF (λ)

 

c j d j e

 

c jd j

e

 

e

 

 

k

 

dF (λ)

j1k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j1k 1

 

 

 

 

 

 

 

 

 

 

n

it j λ

 

n

 

 

it

 

λ

 

 

n

it j λ

 

n

 

 

 

 

dk e

 

dF

c j e

,

d j

c j e

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j1

 

j1

 

k 1

 

 

 

 

j1

 

 

eit j λ .

Таким образом, соотношение (6.3) устанавливает взаимно однозначное соответствие между элементами линейных многообразий M1 и L1 , причем (6.3) сохраняет скалярные произведения. Это соответствие с помощью стандартных рассуждений и с сохранением скалярного произведения распространяется на соответствие между элементами M и L .

Распространение происходит по следующей схеме.

Выберем произвольное x M. Для него можно построить последова-

тельность xn , такую, что xn x, xn M1 в метрике пространства H . Затем, используя сохранение скалярного произведения при соответствии xn fn,

легко показать, что fn – последовательность Коши. Из полноты L2 (dF ) сле-

дует, что для некоторого f L2(dF ) справедливо соотношение fn n f . За-

тем полагаем x f . Изоморфизм при этом отображении сохраняется в силу непрерывности скалярного произведения.

Обозначим Z (a) случайную величину, соответствующую индикаторной функции I( , a](λ) , т. е. для которой I( , a](λ) Z (a) . Проверим, что процесс Z (a) является процессом с ортогональными приращениями. Пусть для чисел a, b, c, d справедливы неравенства a b c d. В силу того, что интервалы (a, b] и (c, d ] не пересекаются, получим

 

 

 

 

 

 

 

 

E (Z (b) Z (a)) (Z (d ) Z (c)) I(a,b](λ) I(c,d ](λ) dF (λ) 0 ,

 

 

 

 

поэтому случайный процесс Z (a)

имеет ортогональные приращения.

Для процесса Z (a) выполняется соотношение

 

 

 

 

E | Z (b) Z (a) |2

 

I(a,b](λ)dF (λ) F (b) F (a) .

32

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

f (a)dZ (a) f (λ).

Действительно, это соответствие верно для индикаторов скольку Z ( ) 0 и по определению процесса Z (a)

a

dZ (a) Z (a) I( ,a) (λ).

(6.5)

I(,a](λ) , по-

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

f из L является пределом в метрике L2 (dF ) последовательности функций fn из L1. Поскольку последовательность fn сходится, она является последо-

вательностью Коши. В силу сохранения скалярного произведения при соот-

ветствии « » последовательность интегралов fn (a) dZ (a) также являет-

 

 

ся последовательностью Коши. Из полноты пространства H и замкнутости

множества функций M следует, что эта последовательность сходится к неко-

 

 

торому элементу из M . По определению этот элемент равен

f (a) dZ (a).

 

 

Соотношение (6.5) доказано. Взяв в (6.5) f (λ) eitλ , получим

 

 

 

eitadZ (a) eitλ.

(6.6)

 

 

Сопоставляя (6.5) и (6.6), получим

 

 

 

X (t) eitadZ (a).

(6.7)

 

 

Это и есть спектральное представление процесса X (t) . ▼

 

В дискретном случае, используя аналогичную технику, можно получить

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

 

π

 

X n einadZ (a).

(6.8)

π

Упражнение. Вывести представление (6.8).

33

7. ОЦЕНИВАНИЕ МАТЕМАТИЧЕСКОГО ОЖИДАНИЯ СЛУЧАЙНОГО ПРОЦЕССА

Определение 7.1. Линейный оператор PI , заданный на замкнутом подпространстве M , называется оператором проектирования на подпространст-

во I , I M ,

если для любого

 

X M выполняются

соотношения PI X I ,

P

2 X P X

и для любого Y I

выполняется соотношение P Y Y . ▼

I

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

Теорема 7.1 (закон больших чисел). Пусть X n – стационарная в широ-

ком смысле последовательность. Тогда предел

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

n 1

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

 

X j Y

 

(7.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n m) n m j m

 

 

 

всегда существует.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если X

n

U n X

0

– представление процесса X

n

с помощью унитарного

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оператора U , то предел в соотношении (7.1) представляется в виде Y PI X0 ,

где PI оператор проектирования на множество I всех векторов, инвари-

антных относительно оператора U . ▼

 

 

 

 

 

 

 

 

 

Доказательство. Рассмотрим операторы

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

n1

 

 

 

 

 

 

 

 

Vm,n

 

 

 

 

 

 

U j .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

m j m

 

 

 

 

Очевидно, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n1

 

 

 

 

 

 

 

 

Vm,n X0 =

 

 

 

X j .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

m j m

 

 

 

 

Пусть I

– множество всех векторов, инвариантных относительно опера-

тора U , т. е.

I {Z M :UZ Z}. Если X I , то Vm,n X X , то предел (7.1)

существует и равен X .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть N – множество всех векторов, представимых в виде X UZ Z

для некоторого Z M , т. е. N {UZ Z : Z M}. Тогда

 

 

 

 

 

V

 

X

U nZ U mZ

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m,n

 

 

 

 

 

 

 

 

 

n m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

34

Поскольку U j унитарные операторы, то норма числителя ограничена, поэтому в рассматриваемом случае Vm,n X 0 при (n m) .

Докажем, что подпространства I , N являются ортогональными дополнениями друг друга.

Сначала докажем включение I N , где N – ортогональное допол-

нение N .

 

 

 

 

 

 

Рассмотрим произвольный элемент X из I

и произвольный элемент Z

из M. Поскольку X инвариантен относительно как оператора U , так и об-

ратного ему оператора U 1, то справедливы равенства

 

 

X , UZ Z

 

X , UZ X , Z

 

 

 

U 1X , Z X , Z

X , Z

X , Z

0.

Следовательно, X N , откуда следует, что I N .

Покажем, что N I .

Рассмотрим произвольный элемент X из N . Тогда для любого Z M справедливы равенства

0 X , UZ Z

 

X , UZ X , Z

U 1X X , Z .

Из этого соотношения

следует равенство

X U 1X . Следовательно,

X I , откуда следует, что N I .

Любой элемент X M может быть представлен в виде

 

X X I X N ,

(7.2)

где X I I , X N N, поэтому

 

Vm,n X Vm,n X I Vm,n X N X I Vm,n X N .

 

Соотношение (7.1) следует из этого равенства, поскольку Vm,n X N 0.

В представлении (7.2) X I PI X , X N PN X , поэтому

справедлива

вторая часть теоремы 7.1. ▼

Следствие. Предельная случайная величина Y в теореме 7.1 равна нулю тогда и только тогда, когда

35

 

 

 

 

1

N 1

 

 

 

 

 

 

 

 

K ( j) 0.

 

(7.3)

 

 

 

 

 

 

 

 

 

 

N j 0

 

 

 

 

Доказательство. Положив в (7.1) n m N и используя теорему 7.1,

получим

 

 

 

 

 

 

 

 

 

 

 

 

1

m N 1

 

 

1

N 1

Y , X m

lim

 

 

X j , X m

 

lim

K ( j) ,

 

 

 

 

N N

 

j m

 

 

N N j0

поэтому, если Y 0 , то выполняется условие (7.3).

Наоборот, если выполнено условие (7.3), то (Y , Xm) 0 для любого X m. Отсюда следует, что (Y , X ) 0 для всех X M , поэтому Y 0 . ▼

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

 

 

 

1

n 1

 

 

 

 

X j .

(7.4)

 

 

 

 

 

 

 

n j 0

 

 

 

Действительно, пусть EXn a . Предположим, что условия следствия из

теоремы 7.1 выполнены для случайного

процесса Yn Xn a. Тогда

1

n1

 

X j a .

 

 

 

 

 

n j0

 

Поскольку эта оценка H состоятельна, то она также и слабо состоятельна.

Теорема 7.2. Предельная случайная величина Y в теореме 7.1 равна нулю тогда и только тогда, когда спектральная функция F (x) непрерывна.

Доказательство. Определим обобщенную функцию δ(λ) равенствами:

 

 

 

 

 

 

 

π

 

 

δ(0) 1, δ(λ) = 0

при λ 0 и

δ(λ) dZ (λ) dZ ({0}) для любой допустимой

 

 

 

 

 

 

 

π

 

 

меры Стильтьеса dZ . Здесь dZ ({0}) обозначает меру точки ноль.

Сравним два предельных соотношения

 

1

N 1

 

1

N 1

π

π

1

N 1

 

X ( j)

 

 

eijλ dZ (λ)

eijλ dZ (λ)

 

 

 

 

 

N j0

 

N j0 π

π

N j0

36

 

 

 

π

 

 

 

 

 

 

 

 

δ(λ) dZ (λ) dZ ({0}),

(7.5)

 

 

 

π

 

 

 

 

 

1

N 1

1

N 1

π

π

1

N 1

 

K ( j)

 

eijλ dF (λ)

eijλ dF (λ)

 

N

N

N

 

j0

j0 π

π

j0

 

 

π

 

 

 

 

 

 

 

δ(λ) dF (λ) dF ({0}).

 

 

 

(7.6)

π

Поскольку E |dZ ({0})|2 dF({0}) , равенства нулю правых частей соотношений (7.5), (7.6) равносильны. Они выполняются только тогда, когда спектральная функция непрерывна. ▼

Пример 7.1. Рассмотрим прим. 5.2. В нем подпространство M гильбертова пространства H является n-мерным. Выберем в M базисные векторы

по формулам e j

 

Aj

 

 

 

 

 

 

 

Aj

, j 1,..., n. Тогда векторы из M однознач-

 

 

 

 

 

 

 

 

 

 

 

E | A

j

|2

 

 

Aj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

||A ||eiλ1t

 

eiλ1t

0

0

 

 

 

||A ||

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

||A2||eiλ2t

 

0

eiλ2t

 

 

 

 

||A2||

 

 

 

 

 

 

 

 

Ut X0.

Xt

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

||A ||e

iλ

n

t

 

0

 

iλ

n

t

 

||A ||

 

 

 

 

 

 

0 e

 

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

Если, например,

λ1 0, а остальные λ отличны от нуля, то инвариант-

ными будут векторы, у которых только первая компонента отлична от нуля. Рассматривая процесс только при t 0, 1, ... , получим

 

 

 

 

 

1

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

N 1

 

 

 

 

 

 

 

A1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

eiλ2 j

0

 

 

 

 

 

 

 

 

 

 

 

N 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

N j0

 

 

 

 

 

 

 

 

A2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

X

j

lim

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

N N

j0

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

N 1

iλ

 

j

 

 

 

 

An

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

e

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

A

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| An ||

 

 

0

 

 

 

 

 

Здесь использовано соотношение

 

1

N 1

1

 

iλNj

1

 

 

 

 

 

 

1,

если λ 0,

 

eiλj

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

iλj

 

 

 

 

 

 

 

 

 

 

 

 

 

N j0

 

 

1

 

 

N

0,

 

если λ 0.

 

N

e

 

 

 

 

 

 

 

Предел Y есть проекция вектора X

 

A

 

 

 

A

 

...

 

 

 

A

 

 

 

т на инвариант-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

1

 

 

 

2

 

 

 

 

 

n

 

 

 

 

ное подпространство.▼ Пример 7.2. Рассмотрим процесс скользящего среднего из прим. 5.7.

Пусть N – пространство, натянутое на случайные величины n}, N .

 

 

 

Тогда для любого d N справедливо представление d

dnξn и выпол-

 

 

n

 

 

 

нено условие |dn|2 . Унитарный оператор U в пространстве N дей-

n

 

 

 

 

 

ствует следующим образом: Ud

dnξn 1, поэтому он в координатной

n

последовательности di тоже является оператором сдвига. В координатной последовательности инвариантное пространство состоит из последовательностей, состоящих из констант. Среди них квадратично-суммируемой последовательностью является только последовательность, состоящая из нулей. Таким образом, всегда Y 0 .

Это можно было установить также, используя следствие к теореме 7.1. Действительно, было показано, что для модели скользящего среднего

 

 

1

N 1

K(m) 0,

|m| , поэтому

K ( j) 0 .

 

 

 

N j 0

38

8. ОЦЕНИВАНИЕ КОВАРИАЦИОННОЙ ФУНКЦИИ

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

 

 

 

(m n),

 

 

 

 

K(n) EX (m) X

поэтому

выборочная ковариационная функция

ˆ

1 N 1

 

 

 

 

 

 

 

 

 

K (n)

 

X ( j n) X

( j) является несмещенной

оценкой

 

 

функции K (n).

 

 

 

 

N j0

 

 

 

 

 

 

 

Даже для постановки

задачи о

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

в среднем квадратичном

оценки

ˆ

 

 

 

ˆ

 

2

. Это усло-

 

 

 

 

 

 

 

 

 

K (n) потребуется дополнительное условие E

K (n)

 

 

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

Теорема 8.1. Пусть X (t) − стационарный вещественный гауссовский процесс с дискретным временем и с математическим ожиданием, равным ну-

лю. Для состоятельности оценки ˆ необходимым и достаточным является

K (n)

условие

1N 1

K 2 (n) 0.

N n0

Это условие равносильно непрерывности спектральной функции. Доказательство. Для простоты будем считать исходный процесс веще-

ственным. В силу несмещенности оценки ˆ условием, равносильным со-

K (n)

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

 

 

 

 

 

 

 

 

ˆ

(8.1)

 

 

 

 

 

 

 

 

DK(n) 0.

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

во тождество (см. [1])

 

 

 

 

 

 

 

 

E(X1X2 X3X4 ) E(X1X2 ) E(X3X 4 )

 

 

E( X1X3) E( X 2 X 4 ) E( X1X 4 ) E( X 2 X3).

Используя это равенство, получим

 

E X n k X k K (n) X n X 0 K (n) EX n k X k X n X 0 K 2 (n)

EX n k X k EX n X0 EX n k X n EX k X 0

 

EX

n k

X

0

EX

k

X

n

K 2 (n) K (n)2 K 2

(k)

 

 

 

 

 

 

39

K(n k)K (n k) K 2(n) K 2(k) K (n k)K (n k).

Следовательно, в гауссовском случае достаточным условием для состоя-

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

ˆ

 

 

 

 

 

 

 

 

 

K (n) является условие: при любом n

 

 

 

ˆ

 

1

N 1

 

 

2

 

 

 

 

DK (n)

 

 

 

K

 

(k) K (n k)K (n k)

 

0.

(8.2)

 

 

N k 0

Из неравенства

K (n k) K (n k) K 2(n k) K 2(n k)

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

N 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

K 2 (n) 0.

(8.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

N n0

 

 

 

 

 

 

 

 

С другой стороны, если оценка

ˆ

 

 

 

 

 

 

 

 

K (n) состоятельна, то, используя (8.2),

 

 

 

 

 

ˆ

3

 

N 1

2

(k) 0.

 

 

 

 

 

 

получим DK (0)

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N k 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь докажем второе утверждение теоремы 8.1. Справедливы равен-

ства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

N 1

 

 

1

N 1

π

 

 

 

 

2

π π

1

N 1

 

K 2 (k) =

 

eikλdF (λ)

 

 

ei(λ – ν)k dF (λ) dF (ν)

 

 

 

 

 

 

 

 

 

 

N k 0

 

 

 

 

 

π

 

 

 

 

 

π π

 

N k 0

 

 

 

 

N k 0

 

 

 

 

 

 

 

 

 

π

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fN (λ, ν) dF (λ) dF (ν),

 

 

 

 

 

 

 

 

 

 

 

 

π π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,

 

 

 

 

если λ = ν,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i(λ– ν)N

 

 

 

 

 

 

 

 

 

1, если λ = ν,

 

 

 

 

fN (λ, ν)

 

1 e

 

 

 

, если

λ

ν

δ(λ – ν)

 

 

 

 

 

 

 

 

 

 

 

i(λ– ν)

 

 

 

 

0, если λ ν.

 

 

 

 

 

N 1 e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что обобщенная функция δ(λ) уже была использована в соот-

ношениях (6.5), (6.6). Таким образом,

 

 

 

 

 

 

 

1

N 1

 

 

π

π

 

 

 

 

 

 

 

 

π

 

 

 

 

K 2 (n)

 

 

 

δ(λ ν) dF (λ) dF (ν) F ({λ}) dF ({λ}) dF 2({λ}).

 

 

 

 

 

N n 0

π π

 

 

 

 

 

 

 

 

π

 

 

λ

40