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

книги из ГПНТБ / Хеннан, Э. Многомерные временные ряды

.pdf
Скачиваний:
27
Добавлен:
23.10.2023
Размер:
21.96 Mб
Скачать

220 Гл. I I I . Прогнозирование и фильтрация случайных процессов

6. Вывести рекуррентные формулы, приведенные перед теоремой 10', кото­

рые выражают s(0>(га +

1) через s(0>(п) и

s^v"*"^(n) через s ^ ( n ) .

7. Показать, что если фильтр

 

 

т

т

X (г)

-J- ^ a (s) X (t s) ds, ^ I а (s) | ds < oo,

 

о

о

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

вперед любой многочлен степени d~.

то его частотная характеристика имеет

вид

 

d

 

 

h ( і г ) = £

(і т ),

 

о

 

где hд (z) аналитична.

Часть II

СТАТИСТИЧЕСКИЕ ВЫВОДЫ

Г л а в а I V

Законы больших чисел. Центральная предельная теорема

і. ВВЕДЕНИЕ

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

ны, мы не сможем оставаться на прежнем уровне общности, так как бессмысленно, например, рассматривать статистические проблемы, основываясь только на предположении интегрируемости / (Я). (С этой точки зрения некоторые результаты первой части являются излишне общими.) Это не означает, что допустимо какое-либо ослаб­ ление математической строгости, просто нам придется наложить некоторые дополнительные ограничения, которые, хотя и жестки с математической точки зрения, не являются таковыми, как мы надеемся, с практической точки зрения.'

Мы почти полностью ограничим нашу статистическую теорию случаем дискретного времени. Это ни в коей мере не является необ­ ходимым, однако представляется целесообразным, так как если бы мы этого не сделали, то потребовалось бы слишком много места для ясного изложения. Кроме того, почти все способы анализа времен­ ных рядов реализуются теперь на вычислительных машинах. В таких случаях дискретизация непрерывных данных необходима, и дискрет­ ное время отвечает существу дела. Это, конечно, не означает, что изучение в предыдущих главах процессов с непрерывным временем было ненужным, ибо обычно в основе анализируемых данных лежат непрерывные явления.

Мы не собираемся давать подробные инструкции, как следует производить те или иные расчеты (хотя часто будем обсуждать вычис­ лительные проблемы и подробности вычислений).

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

224 Гл. IV . Законы больших чисел. Центральная предельная теорема

2.ЭРГОДИЧЕСКАЯ ТЕОРИЯ ДЛЯ ПРОЦЕССОВ, СТАЦИОНАРНЫХ

ВУЗКОМ СМЫСЛЕ

Напомним,

что (векторный) процесс

называется стационарным

в узком

смысле, если

совместное

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

х (щ + т),

X (п2 + т),

. . ., X (nk +

т)

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

от пи п2,

. . ., «ь и не

зависит от

т для любых к,

щ, . . ., nk и т.

Мы сформулируем без

доказательства несколько фундаментальных теорем, касающихся усиленного закона больших чисел. Мы не даем их доказательств, потому что они длинны, сложны и не могут быть упрощены (по крайней мере автором). Эти доказательства легко найти в литера­ туре *).

Мы можем считать (см. § 1 гл. I), что случайные величины xj (п) определены на одном и том же вероятностном пространстве (Й, Л , Р). Класс Л всех измеримых подмножеств этого пространства содержит все цилиндрические множества

{® I (хі (Л) (гса, со), к =. I, . .., т) 6 В},

где' (зщю (nh, со), к = 1, . . ., т) — точка в Л"1 с указанными коор­ динатами, а В — борелевское подмножество Rm. На самом деле можно считать, что Л — наименьшая о-алгебра, содержащая все цилиндрические множества. На цилиндрических множествах можно

определить оператор сдвига

S ^ T S ,

который

на

множество

S = {С О I Xj(iy ( Щ , СО)

< 1 < Ц , . . . ,

Хцт) (llm t

0 ) )

й щ }

действует по формуле

 

 

 

 

TS = {a\Xj (^(Ні + І, со) < öl, . . ., .Tj(m)(nm-f 1, со) < am}.

Так определенный оператор продолжается на все множества из Л- Легко видеть, что Т сохраняет меру, так что S и TS всегда имеют одну и ту же вероятность. Очевидно, Т~г также корректно определен аналогичным образом и также сохраняет меру. Поэтому можно определить оператор сдвига, который мы также обозначим Т, на пространстве случайных величин (измеримых функций) на й по формуле Tf (со) = / (Г-1®). Если Tf — f, то мы будем говоритъ, что f инвариантна (аналогично для S).

Заметим, что Тх} (п, со) = xj (п-j- 1, со), так как для всех а

{со I ТXj (п, со) а] = {со | Xj (п, Г-1со)^ а} = Т {со | Xj (п, со)^а} —

= {со I Xj (п -f 1, со)^а}.

Положим

_ N

X n = N ~ 1 2 х ( п ) ■

Имеет место следующая эргодическая теорема:

с-) См., например, Дуб [1953] и Биллингслей [1965].

2. Эргодическал теория для процессов, стационарных в узком смысле 225

Т е о р е м а 1. Если х (п) стационарный в узком смысле процесс, причем Е {| Xj (п) |} < оо, / = 1, . . р, то существует

вектор X инвариантных случайных величин, для которых Е {| Xj |} < < оо и

lim xN --- X п. н.,

Лт—>со

Е (п)) = Е (х).

Единственный минус этого факта в том, что х не обязан быть постоянным с вероятностью 1. Очевидно, что это требование будет выполнено, если инвариантными случайными величинами являются только постоянные с вероятностью 1, или, что то же самое, если всякое множество S , такое, что S = TS (с точностью до множества нулевой вероятности), имеет вероятность 0 или 1. В этом случае процесс X (п) называется эргодическим. (Такой процесс и оператор сдвига Т иногда называются также «метрически транзитивными».)

Свойство эргодичности не может быть установлено по наблюде­ ниям за единственной траекторией. Это можно показать на примере

скалярного процесс.а, имеющего скачок в спектре,

из-за

которого

в X {п) возникает

компонента

2 Re {z {dX0) е~іп Ь}.

Тогда, как мы

сейчас покажем,

| z (dX0) | —

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

величина.

Однако на основании одной траектории нельзя судить, является ли

она почти наверное постоянной. Чтобы показать, что

| z (dX0) |

инвариантна, рассмотрим точки непрерывности

, Х2 функции F (X),

такие, что Xj < Х0 <

Х2. Из теоремы 3" гл. II

следует,

что

 

Т (z (Х2) —z (А,і))= j e~i%z (dX),

 

 

так

что величина

hi

 

 

Xi

 

 

 

 

 

 

II

T (z (X2) - z (X,)) -

e-ao (Z (X2) - z (Xi)) II2 = j 4 sin2 V2 (X -

X0) F (dX)

стремится к нулю, когда Хь Х2 приближаются к Х0. Отсюда

Tz (dX0) = exp (—iX0) z (dX0)

и I z {dX0) I2 инвариантна, что и утверждалось. В гауссовском случае отсутствие скачков в спектре также и достаточно для эргодичности х). Несколько обобщая, имеем следующий результат. Рассмотрим про­ цесс

оо

Х { п ) =

2

ein X jZ j,

E ( z j Z k ) ^ = 8 j f j ,

 

j==—00

 

оо

 

 

 

 

 

 

Z - j — Z j ,

 

X_j —

X j ,

2 / j

° ° ■

J) Доказательства результатов, приводимых в атом абзаце, читатель может найти в книге Розанова [1963].

15 э. Хешіан

226 Гл. IV . Законы больших чисел. Центральная предельная теорема

Тогда для

эргодичности х (п) в первую очередь необходимо, чтобы

I Zj I были

почти наверное постоянны для любого / . (Это сразу же

исключает гауссовский случай.) Таким образом, случайность входит в X (п) только через фазы колебаний. Во-вторых, если S — любое конечное подмножество Х }-, такое, что для некоторых целых чисел n i f

2 mjQj = 0

(mod 2л),

то должно быть

 

2 nijQj = ср (mod 2л),

zj = | Zj I егвц

bj£S

 

где ф и. н. постоянна. Необходимость этого условия вновь довольно

очевидна, ибо если

 

Г[ = ехр г ( 2

,

/■j £ S

 

ТО

 

Ti] = exp і I 2 Щ (0J +

А.}) I = 11-

XjCS

 

Следующее условие, более сильное, чем эргодичность, смысл

которого

несколько более очевиден,— это

условие

 

 

lim Р (А П

Т~пВ) = Р (А)

Р (В),

А, В еШ.

(2.1)

 

П —► ОО

 

 

 

 

Если

это условие

выполнено, то

х (п)

называется

процессом

с перемешиванием. Чтобы показать, что перемешивание влечет за собой эргодичность, возьмем в качестве А = В инвариантное собы­ тие; тогда получаем Р (А)2 = Р (А), откуда Р (А) — 0 или 1. Грубо говоря, если события, определяемые значениями х (п), сильно раз­ деленными во времени, почти независимы, то выборочное среднее с вероятностью 1 сходится к среднему значению процесса. Более сильным условием, чем (2.1), является «условие сильного перемеши­ вания», которое мы сейчас введем. Рассмотрим все события (т. е. все множества из Л), определяемые только значениями х (п) при п ^ р.

Эти события образуют ст-алгебру $ ѵ. Аналогично, пусть ерч есть

a-алгебра событий, определяемых х (п) при п

q. Пустъ существует

положительная функция g от целочисленного аргумента, такая, что

g (п) ->- 0

при

I п I —>- оо и

 

I Р (В П

F) -

Р (В) Р (F) I < g (q - р), В 6 .«V F е

(2-2)

Тогда говорят, что х (п) удовлетворяет условию сильного перемеши­ вания. Теперь мы утверждаем, что события, разделенные во времени

промежутком q р, приближаются к

независимости

равномерно

в том смысле, что

отклонение от нуля

левой части

(2.2)

зависит

только от

q — р,

но не от

конкретных

рассдіатриваемых

В и F.

Из этого

условия,

очевидно,

вытекает (2.1), а значит, и

эргодич­

2. Эргодическая теория для процессов, стационарных в узком смысле 227

ность {я (тг)}. Условие сильного перемешивания является довольно жестким, но, с другой стороны, и более привлекательным. Например, есть все основания полагать, ято, скажем, поведение морских волн, которые омывали берега много лет назад, совершенно не связано с поведением волн, омывающих берега сегодня, и что это верно для всех явлений, определяемых этими волнами. Правы мы или нет в таком предположении, в настоящей книге нам придется существен­ но использовать это условие, введенное в статистические исследова­

ния М. Розенблаттом. Впрочем, оно

понадобится нам не сразу,

и мы ввели его сейчас лишь ддя того,

чтобы указать на его связь

с перемешиванием и, следовательно,

с эргодичностью х).

Пусть X (п) — стационарный в узком смысле эргодический про­

цесс; рассмотрим величины х (п) х (п +

т)'. Если зафиксировать т

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

смысле,

так что

определена Г (т) = Е (х (п) х (п + т)'). Если

(7г)}

эргодичен,

то эргодичен и

процесс

 

 

(п) X (п +

то)'},

ибо оператор сдвига для этого процесса совпадает с оператором сдви­ га для {а: (п)} на своей области определения (т. е. на множествах, определяемых условиями на xj (п) xh (п + т)), а значит, любое инвариантное множество имеет меру 0 или 1.

Т е о р е м а 2. Если х (п) стационарный в узком смысле эрго­ дический процесс, такой, что Е {| Xj (п) |} < оо, то

ІѴ

lim jjr 2 х in) = Е (п)) п. н.

IV—>-оо •

7 1 = 1

Если, кроме того,

Е (xj (п)2) <

оо,

то

 

 

N

 

 

lim -jr- 2

х (т ) х +

п)' = Г (п) п. и.

іѴ-ѵоо

1

 

772=

 

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

X {т) X (т + п)'.

*) Розанов [1963] обсуждает смысл этого условия для гауссовских процес­ сов и устанавливает ряд интересных соотношений. Мы используем это условие главным образом в случаях, когда предположение о гауссовости не делается.

15*

228 Гл. IV . Законы больших чисел. Центральная предельная теорема

Если X (п) — процесс с перемешиванием, то любой новый процесс, определенный соотношением

уЦп. со) = / (Т~па).

где / — векторнозначная измеримая функция на (й, А , Р), также будет процессом с перемешиванием. В самом деле, если А и В — множества вида

I

(&)і

®)і

k

• • •» гі) £ С} —

 

 

=

{со

I (fkh

k =

1. , r .,

m) 6 C),

где С — борелевское

множество

в Rm, то, как

легко

видеть, усло­

вие (2.1) выполняется для класса измеримых множеств J-.v, порож­

даемого множествами вида А

и В, в вероятностном пространстве

У, А Ч,Р), соответствующем

процессу у (п). Нетрудно показать,

что последовательность независимых, одинаково распределенных случайных векторов е (п) с нулевым средним и конечной матрицей

ковариаций

обладает

свойством перемешивания

(мы оставляем это

в качестве

упражнения). Отсюда вытекает,

что

если

 

X (п ) = 2

- 4 (у) е (га — / ) ,

2

II

-4 (/) II2 < оо,

то X (га) — также процесс с перемешиванием и, следовательно, эргодичен. Процессы вида (2.3) иногда называются «линейными». Мы, однако, будем использовать этот термин в более узком смысле. (Мы знаем, что если процесс х (п) гауссов с абсолютно непрерывным спектром, то он может быть представлен в виде (2.3).)

Т е о р е м а 3. Если х (га) дается соотношением (2.3), где е (га) независимы и одинаково распределены, то х (п) процесс с переме­ шиванием и, следовательно, эргодичен.

3.ЭРГОДІІЧЕСКАЯ ТЕОРИЯ ДЛЯ ПРОЦЕССОВ, СТАЦИОНАРНЫХ

ВШИРОКОМ СМЫСЛЕ

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

3. Эргодическая теория для процессов, стационарных в широком смысле 229

Т ео р ем а 4. Пустъ х (п) — стационарный процесс1) и

 

 

Л

 

 

 

 

 

 

 

Л

 

X (п) = j

e~in%z (dl),

Г (л) = Е (х (т) х(т-\- ?г)') = | einX F (dl).

 

 

- л

 

 

 

 

 

 

 

Л

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l.i.m . xN = z (0) — z (0—),

 

 

 

ІѴ-І

 

 

 

N->со

 

 

 

 

 

 

Г ( n ) = F ( 0 ) - F (0 - ) = Е {(2 (0)-2 (0 -)) (* (0 )-« (0-))*}.

l i m ^

2

И-ЮО

 

 

 

 

 

 

 

 

 

=

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

 

Положим ф (X) =

0 при X Ф 0; ф (0) =

1.

Тогда,

в силу

мажорированной сходимости,

 

 

 

 

 

я

 

N

 

 

 

 

 

 

 

lim

\

| - ^

2

е- 1пХ- Ф М Г ^ ( ^ ) =

 

 

 

 

N —юо J

I ІѴ

 

 

1

 

 

 

 

 

Л

 

 

 

 

 

 

 

 

 

=

j { + ■ » ( “ - з д ч } * т - о.

Поэтому

 

 

 

 

 

N

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l.i.m . xN = l.i.m . -jjr- у.

f e -inXz(dl) =

 

 

 

 

N -> o o

 

N - ю о

1S ■“ 1 -JЯ

 

 

 

 

 

 

 

 

 

Я

N

 

 

 

 

 

 

 

Л

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= j Ф (X) 2 (dX) = z (0) —z (0 —),

а

эта

величина

имеет

ковариационную

матрицу, указанную

в формулировке теоремы. Чтобы доказать^'второе утверждение, заметим, что

N - 1

я

N —1

e™F(d1)=

JV—►*ітсхэ4 2 Г (”) =ІѴііт—^COJJ

ПТ■* 2

О

-Jt

0

 

Jt

=j <p(^)E(dX) = E ( 0 ) - E ( 0 - ) ,

-Л

также в силу мажорированной сходимости.

*) В с ю д у в дал ь н ей ш ем , к а к

и в н ач ал е к н и ги , «стационарность)) б у д ет

озн ачать «стац и он арн ость второго

п ор я дк а » .

230 Гл. IV . Законы больших чисел. Центральная предельная теорема

Из теоремы 4. сразу следует, что если х (п) — стационарный

процесс со средним р, то вектор выборочных средних x N сходится в среднеквадратичном к р тогда и только тогда, когда спектр процесса не имеет скачка в нуле. Это утверждение, часто называемое, в отли­ чие от «индивидуальной эргоднческой теоремы» предыдущего пара­ графа, «статистической эргоднческой теоремой», принадлежит Дж. фон Непмаиу.

Мы знаем, что если х (п) стационарен в узком смысле, то x N имеет

предел почти наверное. Последняя теорема утверждает, что x N сходится в среднеквадратичном к случайной величине z (0) — z (0—). Так как из сходимости в среднеквадратичном и сходимости почти

наверное вытекает сходимость по вероятности, то х к почти наверное должно сходиться именно к z (0) — z (0—).

С л е д с т в и е 1. Если х (п) стационарный в узком смысле

(а также в широком смысле) процесс с нулевым средним, то х N схо­ дится п. н. к 0 тогда и только тогда, когда спектр процесса не имеет скачка в нуле.

Совершенно аналогичные рассуждения, примененные к величинам

Л'

-jr 2 х («)еіп%.

показывают, что они сходятся в среднеквадратичном и и. н. к слу­ чайной величине z (X) z (X—); при этом они сходятся к нулю тогда

итолько тогда, когда спектр процесса не имеет скачка в точке X. Сейчас мы докажем теорему Дуба [1953], которая дает условия,

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

Т е о р е м а

5. Пустъ х (п) удовлетворяет условиями теоремы 4

и Е (гс)) =

0.

Тогда, если существуют константы К >

0 и а > 0,

такие, что

 

 

 

N - 1 N - i

N - 1

 

 

І і У 1- 0 і - т ) = ± 2

(3.1)

0

0

-Л’+і

 

то xN сходится к нулю с вероятностью 1.

Очевидно, эту теорему достаточно доказать для одномерного про­ цесса. В этом случае левая часть формулы (3.1) есть дисперсия х^. Найдем ß > 1, такое, что ßa >-1. Тогда при N >- имеем Е (х%) ЦІ Если е > 0 и N (М) — наименьшее целое, большее или

Соседние файлы в папке книги из ГПНТБ