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

книги / Статистический анализ временных рядов

..pdf
Скачиваний:
14
Добавлен:
19.11.2023
Размер:
47.02 Mб
Скачать

<Ш )

<iv)

(v)

 

 

 

УПРАЖНЕНИЯ

 

403

8Л2*

(2*)! Г ( т * + * ) _

 

 

 

 

И Г (4 Ч

 

 

 

 

 

_ 2г‘г(‘ + т)г( т " +‘)

Л-0, 1, ... ,

 

 

 

 

 

V 2 А - 1

Г (

т

М »

 

л - 1, 2.........

<0 ,

 

 

 

 

 

 

 

 

fc =

0 , 1, . . . .

 

г ( т ) г (

т

лг+ ‘ ) ’

 

 

 

(Разд. 6.8.1) Покажите, что

 

 

 

2я

 

 

я

 

 

 

j.

 

 

 

 

 

J (а + b cos х) dx =

2^ log (а +

b cos *) cf* =

 

 

 

=* 2я log —

(а + j/ а 262),

а > b > 0.

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

ляя cos х == cos2j

х — sin2

л; и полагая затем у =

6) tg

а262.)

70.

(Разд, 6.8.1) Докажите, что

 

 

 

2я

г

 

e - i2nt/T)k =

271

k »

0, 1. . . .

Г — 1.

Y

2

('ei2nt,T +

j (е?и + e~ iu)k du,

 

/=1

 

о

 

 

 

71.(Разд. 6.8.1) Покажите, что вплоть до порядка Т — 1 семиинварианты

Ох для циклического случая совпадают с соответствующими семиинвариантами случайной величины, имеющей производящую функцию моментов (20), сравнивая

разложения (18) и (19) в ряды по степеням 0. (Указание. Использовать упр. 70.)

72.(Разд. 6.8.1) Покажите, что если X — случайная величина с нулевым средним, дисперсией а2, значения которой с вероятностью 1 сосредоточены в ин­

тервале (a, b)t то а2< — ab. (Указание. Показать, что

ь

j (x — a)(b — x )f (х) dx > 0,

а

где / (я) — плотность распределения случайной величины X.)

73. (Разд. 6.8.1) Покажите, что нечетные моменты г\ до порядка Т — 2 вклю­

чительно и четные моменты до порядка Т —Л включительно совпадают с соот-

404

СЕРИАЛЬНАЯ КОРРЕЛЯЦИЯ

Гл. 6 .

ветствующими моментами функции

 

i 2

Т * 2 -)

^(г-а/гЛ

г +

1 ГЛ

Г ( ±

г ) г ( — )

V1

7^=~t

W *

(Заметьте, что эта функция не является плотностью, поскольку она отрицательна

при 1 - 2 Ц Т + 1) <

г\ <

1.)

74.

(Разд. 6.9.2)

Покажите, что

 

 

/ _

г3 ~ г \ Г3 "Г г \ *— ^Г1Г2 + г \г2

 

 

3

1 - 2г1* + 2г^ - г?

75.

(§ 6.10) Покажите,

что производная вероятности, дополнительной (до 1)

к (14), по переменной R* дается выражением (9) с g | (rj | 0), задаваемой выраже­

нием (53) § 6.7, с заменой г и г\ на R* и заменой pi на р.

76. (Разд. 6.11.1) Покажите, что матрица квадратичной формы в показателе

экспоненты равна — i-

= — 1/(2а2) Т'Т, где

“V I - р2 0 0 . .. 0 о"

р1 0 . .. 0 0

0

р 1 . .. 0 0

Т =

0 0 0 . .. 1 0

00 0 . .. р 1_

77.(Разд. 6.11.1) Убедитесь в том, что решения кубического уравнения (10)

не являются рациональными функциями полиномов P Q , P I и Р0.

78. (Разд. 6. 11. 1) Покажите, что

является решением уравнения (10) с точностью до о (1/Т ). (Указание. Положить в (10) р = —Pi/P 0+ v/T. Отбросить в полученном уравнении члены, имеющие

порядок малости выше 1/7\ и решить его относительно v. Заметьте, что при Т

оо

отношение Р^/Р0имеет (по вероятности) порядок 1/7\)

 

79.(Разд. 6.11.1) Получите уравнения для оценок максимального правдо­ подобия, когда %yt = р.

80.(Разд. 6. 11. 1) Покажите для случая стационарного гауссовского процес­

са первого порядка (%yt = р, $ (yt — р) (ys — р) = а (0) р*

1), что оценкой

максимального правдоподобия для р при известном р (0< р <

1) является

 

Г -1

 

У\ + (1 — р) У yt Ут

р =

t=2

— 2) (1 — р) + 2

УПРАЖНЕНИЯ

40S

и проверьте, что ее дисперсия равна

 

 

Var ц = а (0)

1 + Р

 

— 2) (1 — р) + 2

 

 

 

(Указание. Использовать результаты § 2.4.)

 

81. (Разд. 6. 11. 1) Найдите дисперсию величины у для

стационарного гаус-

 

 

Л

совского процесса первого порядка и сравните ее с дисперсией для р,.

82. (Разд. 6.11.2) Пусть у стационарного гауссовского

процесса первого по­

рядка р = е~Ут (0 < k). Покажите, что при Т -* оо

 

Нш Var |х = а (0)

Г-*-ОО

lim Var у = а (0)

T-+OQ

2

2+ k

P~k 1 + k k2

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

lim Var р < lim Var у < о (0).

Г - * оо

7 -+ 0 0

83.(Разд. 6.11.2) Проверьте, что якобиан преобразования (12) действитель­ но равен 1 — (—Р)г .

».А

84.(Разд. 6. 11.2) Существует ли решение (—Р) уравнения (20), являющееся

монотонно возрастающей функцией от г?

л

85. (Разд. 6.11.2) Пусть (—р7) — корень уравнения (20), расположенный

между —1 и 1, и пусть г = гт. Покажите, что Tk (гт+ Рг) -> 0 по вероятности;

при Т

сю для любого k >

0. (Указание. Покажите, что для достаточно больших

Т со сколь угодно большой

вероятностью выполняются неравенства | гт+

Л

Рг | <с

<

 

А

 

< а <

е и а < ~ P j < Ъдля любых заранее заданных значений е, а и b, — 1

<

& <

1. Затем покажите, что Tk (—РГ)Г—1 -> 0 по вероятности.)

 

* 86. (Разд. 6.11.2) Убедитесь в том, что (21) соответствует критерию отноше­ ния правдоподобия для проверки гипотезы Р = 0.

87. (Разд. 6.11.2) Соответствует ли указанный критерий отношения правдо­ подобия отвержению гипотезы р = 0 при г > п и г < г2?

88. (Разд. 6.11.2) Найдите уравнения максимального правдоподобия, когда

Ш = Р-

89. (Разд. 6.11.2) Найдите уравнения максимального правдоподобия для <т* и р в модели с последовательными разностями, т. е. в модели разд. 6.5.3 с у0=

=

(1 +

Ра)/а2 и

yi = 2р/а2. (Указание. Использовать лемму 6.7.11 с

0=

=

—р/(1 + Р2) для вычисления определителя в соответствующем распределе­

нии.)

 

 

 

 

90. (Разд. 6.11.2) Выпишите уравнения максимального правдоподобия для

оценок параметров а2, Р и ц модели с последовательными разностями, когда

v0

=

(1 +

Р2)/а а , Vi =

2р/о* и gyt = JA.

 

Глава 7

СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ

ПРОЦЕССЫ

7.1. ВВЕДЕНИЕ

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

у0, уъ ... называется случайным процессом с дискретным параметром времени.

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

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

В гл. 5 нас интересовали процессы, описываемые конечным чис­ лом параметров. Теперь мы будем рассматривать процессы, в кото-

7.2. ОПРЕДЕЛЕНИЯ И ПРИМЕРЫ 407

рых число параметров может быть бесконечным. Например, гаус­ совский (нормальный) стационарный процесс имеет конкретные среднее значение и дисперсию, но бесконечное число коэффициентов ковариации. Поставим вопрос: какую информацию о них можно получить исходя из конечного числа наблюдений?

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

Случайный процесс у (f) с непрерывным параметром времени можно определить для 0 < / < оо или — оо < t < оо и рассматри­ вать с привлечением вероятностной меры на пространстве функций у (/). Выборка из такого процесса состоит из наблюдений в конечном числе точек времени, или из непрерывных наблюдений в интервале времени. Пусть, например, процесс представляет собой теорети­ ческую модель наблюдений температуры в некоторой местности. Тогда последовательность ежечасных отсчетов или диаграмма непре­ рывных показаний является выборкой из случайного процесса. Случайный процесс с дискретным временем часто понимают как выборку через равные промежутки времени из случайного процес­ са с непрерывным временем.

7.2.СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ, ОПРЕДЕЛЕНИЯ И ПРИМЕРЫ

7.2.1.Определения

Случайным процессом с дискретным временем называется бес­ конечная последовательность случайных величин уъ у2, ... или ..., У-1, уп, ух.......Наблюдение процесса, часто называемое реализацией* есть точка в соответствующем бесконечномерном пространстве, где определена вероятностная мера. Вероятность определяется на некоторых множествах, называемых измеримыми. Этот класс мно­ жеств включает вместе с любым множеством его дополнение, а так­ же объединение и пересечение счетного числа множеств этого клас­ са; вероятностная мера на этом классе множеств определяется таким образом, что вероятность объединения непересекающихся множеств равна сумме вероятностей отдельных множеств.

4 08

СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ

Г л . 7.

Практически мы интересуемся вероятностями, которые свя­ заны с конечным числом случайных величин. Эти вероятности вклю­ чают в себя функцию совместного распределения Рг {ух < ь и ...

..., Ут < Ьт}, а также вероятности, полученные из нее. Потребуем, чтобы в класс измеримых множеств входили множества {ух <; Ьъ ...

..., Ут < Ьт\. Множества, которые определяются в терминах конеч­ ного числа координат уи, ..., ytn, называются цилиндрическими

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

Наоборот, пусть заданы конечномерные функции совместного распределения (ф.с.р.), тогда на бесконечномерном простран­ стве вероятностную меру можно определить таким образом, что вероятностная мера цилиндрических множеств согласуется с ф. с. р. в том смысле, что может быть вычислена по соответствующим ф. с. р. Это расширение единственно. Согласованность семейства ф. с. р. означает: (i) любая маргинальная (совместная) ф. с. р., полученная из заданной ф.с.р., приводит к соответствующей ф. с. р. меньшего числа переменных, (И) ф.с.р., соответствующие наборам точек параметра t, написанным в различном порядке, согласованны. Точ­

нее, пусть F (alt

..., ап\ tv ...,

tn) означает ф. с. р.

от сц,

...,

ап, лю­

бых tu ...,

tn (из

множества ..., —1, 0 ,

1, ... или

1,

2 , ...) и

п = 1,

2,

.... Потребуем выполнения равенств

 

 

 

 

 

■(1)

F(av

.. . ,

ат, оо, . . . ,

оо; tlt

. . . , tm,

 

 

. . . ,

/„) =

 

 

 

= F(alt . . . , Umt

... »

^m),

m =

1,

. . . , ftt

(2)

F (ax,

. . . .

an; tlt . . . .

tn) = F (a<„ . . . .

aln;

tlt, . . . .

tin),

тде (iv ...,

in) есть произвольная перестановка чисел (1,

...,

п). Со­

гласованное семейство ф. с. р. определяет случайный процесс yt

в

том смысле, что F (аъ ..., а„; tv ..., /„) = Рг {уи < ах, . . . , ytn <

^

а„} для всех а,., t{ и п. Это утверждение принадлежит Колмогоро­

ву (1933), доказательства его мы не приводим.

В качестве примера зададим ф. с. р. как ф. с. р. многомерных нормальных распределений, определив средние значения и диспер­

сии

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

(3)

%yt — т (t),

<4)

%[yt — m (0) [ys — т (s)] = Соv (yt, у,) = о (/, в).

Последовательность т (t) произвольна. Двойная последовательность <т (t, s) должна удовлетворять обычным свойствам ковариаций, а именно а (/, s) = а (s, t) для каждой пары s, t и любая матрица

7.2.

ОПРЕДЕЛЕНИЯ И ПРИМЕРЫ

[сг(/*, </)], /, / = 1, ..., п> положительно полуопределена (или, эк­ вивалентно, неотрицательно определена.— Перев.). Эти условия оп­ ределяют случайный процесс. Такой процесс называется гауссов­ ским. Случайный процесс с дискретным временем называется ста­ ционарным (или стационарным в узком смысле), если распределение величин уи% ..., ytn совпадает с распределением уи+и ..., y tn+t

для любого конечного множества целых чисел {/ь ..., /„} и любого целого /. Это определение эквивалентно требованию того, что веро­ ятностные меры последовательностей {ys} и {f/s+*} совпадают для любого целого /.

Пусть существуют первые моменты, тогда из стационарности

следует

 

 

(5)

%ys = %ys+t, s, t = . . . ,

— 1, 0, 1, . . . .

или

 

 

(6 )

т (s) = т (s +

t) — т

для любых s, t. Поскольку распределения (ytl, yt,) и (ytl+t, Уи+t> совпадают, то из существования моментов второго порядка и стацио­ нарности следует

(7)

о (flt

*4)

= о

+ t; t2+ t).

Подставляя t = 12, получаем

 

(8 )

a (flt t j

= a

 

0 ) = a (tlt — 1%).

В случае нормальных ф. с. р. свойства (5) и (7) эквивалентны ста­ ционарности случайного процесса.

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

если функция средних т (t) и ковариационная функция о (t, s), определенные в (3) и (4), существуют и удовлетворяют соотно­ шениям (5) и (7), т. е. средние постоянны и не зависят от времени, а ковариация любых двух значений зависит только от их расстоя­ ния во времени. Очевидно, любой процесс, стационарный в узком смысле и имеющий конечную дисперсию, является также стацио­ нарным и в широком смысле. (В рассмотренном выше нормальном случае стационарности в узком и широком смысле эквивалентны.)

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

Рассмотрим несколько примеров стационарных случайных про­ цессов.

Пример 7,1. Пусть yt независимы и одинаково распределены с

(9)

%yt = т,

Var */, = ста.

410

СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ

ГЛ. 7.

Тогда

 

 

t — s,

 

<Ю)

o(t, s) = о*

при

 

 

0

при

t=£s.

 

Этот процесс стационарен в узком смысле. Если отбросить условие одинаковой распределенности величин yt (в то время как (9) и (10) сохраняются), то результирующий процесс будет все же стацио­ нарен в широком смысле.

Пример 7.2. Пусть все yt тождественно равны случайной вели­ чине у. Тогда, если существуют два первых момента у,

( 11)

 

 

 

 

 

 

%yt = %у = т,

 

 

 

 

(12)

 

 

 

 

 

 

o(t, s) — Var у =

о2.

 

 

 

Этот процесс стационарен в узком смысле.

 

 

 

 

Пример 7.3.

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

{yt}

следующим

образом:

 

 

 

 

 

 

 

 

 

 

 

 

(13)

У( — 2

И / cos V

 

+

В/ siri Я/0,

t =

• • » - 1 , 0,

1, •.. ,

 

/=i

 

 

 

 

 

 

 

 

 

 

 

 

где Я/ — постоянные, a A t......

5:1> •••»

 

 

- случайные вели

чины, такие, что

 

 

 

 

 

 

 

 

 

 

 

(14)

 

* А , =

ИВ/ =

0,

 

/■=

1 ...

Я,

 

 

 

 

 

 

ю

 

to

 

 

 

 

 

 

(15)

 

 

 

йо

 

Q

 

/ =

1,..

<7,

 

 

 

II ОБ

 

II •j.

 

 

(16)

% AtA j =

g

 

 

= 0, 1 Ф! >

*. / =

1.

• • » q,

(17)

g^B/ =

о,

 

 

 

 

/ =

1,

• • *q-

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

(18)

 

 

 

 

 

 

%yt =

0 ,

 

 

 

 

 

(19)

%ytyt =

Я

 

(Л cos V + Bi sin V ) И / c°s Я/ s + Bf sin kts) =*

2

8

 

 

/ , / s s l

 

 

 

 

 

 

 

 

 

 

 

=

q

[%>A2j cos Xft cos KfS +

gВ/ sin Xft sin A,/$J =

 

2

 

 

/=1

 

 

 

 

 

 

 

 

 

 

 

 

=

*

2

IcosV COS XfS +

sin %it sin Л/s] =*

 

 

2

°f

 

 

 

/«i

 

 

 

 

 

 

 

 

 

 

 

 

 

q

O/COSkf{t — s).

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

/='

 

 

 

 

 

 

 

 

 

 

 

Следовательно, процесс стационарен в широком смысле. Если Л/ и Bj нормально распределены, то yt будут также распределены нор-

7.2.

ОПРЕДЕЛЕНИЯ И ПРИМЕРЫ

411

Bj

Рис. 7.1.

Определение параметров 0/ и /?/.

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

Этот процесс можно также представить следующим образом:

(20)

 

у , = 2

Я/cos (V — в/).

<

 

 

1, о,

1, .. . ,

 

 

 

м

7.1) R}

 

 

 

 

 

 

где

(как

на

рис.

0,

 

 

 

 

 

(21)

 

 

 

$ = А а, +

$ ,

1 = 1 ,

. . . , q ,

 

(2 2 )

 

 

tg 0/ = Bj/Aj,

) =

1,

. . . , q,

 

0 <

0 / <

я,

если

B j> 0 , и

я <

0 / <

2 я,

если

£ / <

0 .

Если ^4/ и 5/ распределены нормально, то R/

пропорционально Xs

с двумя степенями свободы и 0 /

равномерно распределено между О

и 2 я (ввиду симметрии нормального распределения) и не зависит от Я/. Для доказательства стационарности процесса заметим, что

(23)yt+s = 2 Rj cos [У — (0/ — М ],

/=i

 

где (0; — X/s) равномерно

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

интервале длины 2 я.

Таким образом, yt и yt+s

имеют одно и то же

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

А/ и В) имеют распределение, отличное от нормального, то последо­ вательность {yt} не обязательно стационарна в строгом смысле.

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

Пример 7.4. Процесс скользящего среднего.

Пусть ..., о_ь v0, vu ... — последовательность независимых и одинаково распределенных случайных величин, а а 0, а ь ..., aq

412

СТАЦИОНАРНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ

Гл 7.

последовательность q + 1 чисел. Определим процесс скользящего среднего

(24) yt = cc0vt + ajVt-i + • • • +

1, О, 1,... .

Это стационарный случайный процесс. Если $vt = v и Var vt = а2* то

(25)

8 % S5=v(a0 + a 1+

•••

 

 

 

<26)

Cov (yt, yl+s)

I a2(a0as +

+ Oq-aX

s =

0 , ...

<7»

= o

s =

9 ■+■1»

 

 

 

 

 

и последовательность \yt) стационарна в широком смысле. Для того чтобы последовательность {yt\, определенная формулой (24), была стационарна в широком смысле, необходимо только, чтобы величины vt имели одни и те же среднее значение и дисперсию и бьши взаимно не коррелированы.

Мы можем также определить процесс скользящего среднего с бесконечными пределами суммирования. Бесконечная сумма

(27)

2

 

5=0

означает случайную величину уи если она существует, такую, что

(28)

/

п

\ 2

limg[у,— 2

aso<_3 = 0 .

Если выполнено

"-*’00 \

s = 0

/

 

 

 

(29)

 

2 а 3<

оо,

 

 

s = 0

 

такая случайная величина существует (следствие 7 .6 .1). Говорят, что сумма (27) сходится в среднем или среднеквадратичном. Беско­ нечная сумма

(30)

2

asVt~s

определяется как случайная величина yt, такая, что

(31)

lim

в ( yt —

2

a ^ - s V - 0 .

 

n,m -*oo

\

s==t_ m

J

Достаточным условием существования случайной величины yt яв­ ляется

(32)

2 а‘<°°-

 

S = — о о

Условия (29) и (32) являются достаточными, если vs не коррели­ рованы и имеют общие среднее и дисперсию. Но, вообще говоря, эти случайные величины могут быть зависимыми.