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

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

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

60

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

Помимо рассмотренного способа дискретизации процесса х (t) существуют и другие, и ради полноты мы скажем о них несколько слов, хотя они, по-видпмому, ие столь важны с точки зрения прило­ жений. Пожалуй, наиболее часто встречающейся представляется «периодическая выборка с возмущениями», рассмотренная в работе Шапиро и Сильвермеиа [1960]. Предположим, что наблюдения произ­ водятся в моменты п ип, где ип — независимые одинаково распре­ деленные случайные величины, причем дисперсия ип сравнительно мала. Это соответствует ситуации, когда возмущения в некотором измерительном устройстве приводят к тому, что наблюдения произ­ водятся в моменты времени, слегка отличающиеся от намеченных. Полагая

У {п) = X (к 4- н.„),

легко получаем, что

00

Г (?г) = Е (у (т) у' (т 4- п)) = j Г (п 4- и) Р (du),

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

со

со

со

j

j

eK*+'MF(dh)P{du)= j einl | cp (X) \ZF (dX),

— 00

— 00

— CO

где cp (Л) — характеристическая фупкция un. Это равно

Я

j ginXp (^A.)t

—Я

где

со

F (X) = S 11Ф (Ä.■4- 2л/) |2F (X+ 2л/) - 1cp ((2/- 1 ) л) |2F ((2 /- 1) л)].

Функция F — это наибольшее, что мы можем получить из наблюде­ ний у (/г). В случае когда F (X) является постоянной вне интервала

[—я, я], мы, очевидно, можем найти F (А,) по F (X) при условии, что известно Р, а следовательно, и ф. В противном случае при интерпре­ тации наблюдаемого спектра следует соблюдать осторожность.

Шапиро и Сильвермен [1960] рассмотрели также случай, когда наблюдения производятся в точках tn, где (tn 4і-і) — независимые одинаково распределенные (неотрицательные) величины (иначе гово­ ря, tn является реализацией процесса восстановления) с характери­ стической функцией, которую мы вновь обозначим ф (X). Теперь,

3. Дискретизация процесса с непрерывным временем

(31

полагая у (п) — х (tn), имеем

 

 

со

 

Г (?г) = Е (у (т) у' (т +

п)} =

| Т (t) Рт (dt),

 

 

 

— со

 

где Р {П)— функция распределения

величины tn10. Как и раньте,

получаем

 

 

 

СО

 

 

 

Г(я )= J сf(X)F(dX),

/г>0,

(3.5)

00

 

 

 

так что у (п) — также стационарная последовательность. Зададимся теперь вопросом: может ли быть так, что последова­

тельность Г (п), содержащая в себе всю информацию, основанную на свойствах второго порядка, однозначно определяет функцию F (Я) (подчиненную обычным условиям нормировки). Другими сло­ вами, можно ли избавиться от наложения частот, используя опи­ санный способ выборки и подбирая надлежащим образом ср (Я). Ответ дается теоремой 5, принадлежащей Шапиро и Сильвермену [I960], доказательство которой мы опускаем.

Т е о р е м а 5. Для того чтобы спектральная функция F (Я) однозначно определялась по функции (3.5), достаточно, чтобы для любых двух точек (Я1, Я2): Я! Ф Я2, ЯІ5 Я2 6 (—°°, оо), выполнялось условие ср (Я!) Ф ср (Я2) (т. е. чтобы ср (Я) была унивалентна на

(—оо, оо)).

Представляет интерес случай, когда tn порождаются пуассонов­

ским

процессом с интенсивностью

появления событий

р. Тогда

ср (Я)

= (1 + гЯ/р)_1, и достаточное

условие теоремы 5,

очевидно,

выполняется. В этом пуассоновском случае утверждение теоремы легко доказать непосредственно г); в самом деле, полагая для просто­ ты обозначений р = 1, мы видим, что множество А комплексных линейных комбинаций функций вида (1 + іЯ)_п, п — О, 1, . . ., является не только линейным пространством над полем комплекс­ ных чисел, но и алгеброй, т. е. произведение двух таких комплексных комбинаций вновь является комплексной линейной комбинацией

функций

(1 + іХ)~п. В самом деле,

(1 + іЯ)~х (1 — іЯ)' 1 =

= 1 / 2 ((1 +

гЯ)_1 4- (1 — іЯ)-1), откуда и

следует высказанное

утверждение. Таким образом, А является алгеброй непрерывных функций, обращающихся в нуль на бесконечности, замкнутой отно­ сительно комплексного сопряжения и обладающей тем свойством, что для любых двух точек Яі Ф Я2 найдется элемент А , принимающий в них различные значения.,Тогда (см. Иосида [1965], стр. 20) алгебра А плотна в смысле равномерной сходимости в пространстве С (— оо, оо)

J) Этот результат не понадобится в дальнейшем, н читатель может опустпть доказательство.

Ö2

Гл. I I . Спектральная теория векшорпы.ѵ процессов

всех непрерывных функций, обращающихся в нуль на бесконечно­ сти . Таким образом, если бы Г (п) не определяла однозначно F (dk) и существовали бы две разные функции, удовлетворяющие усло­ вию (3.5), то для матричиозначной функции М , которая равна раз­ ности этих функций, мы нмоли бы

оо

 

со

 

 

 

j

ф"(\) M(dk) =

j gFjkj М (dk) == Ü,

О,

— 00

 

— со

 

 

 

и из только что сказанного

следовало бы,

что

 

 

СЮ

 

 

 

 

 

j а (k) М (dk) —0,

а £ С ( — оо,

оо).

 

— со

 

 

 

 

Однако это возможно лишь в случае,

если М (dk) = 0 (Иоснда

[1965], стр.

170), так что

единственность доказана.

Мы не будем рассматривать здесь проблему явного выражения F

через у (п). Практическая ценность такого результата не ясна. Если имеется в виду планируемая выборка, то описанная выше про­ цедура представляется несколько громоздкой н использование перио­ дической выборки с очень малым интервалом выборки может ока­ заться более предпочтительным. Как видно пз формул (3.3) п (3.4), для абсолютно непрерывной F (А.) это будет эффективным средством избежать наложения частот (с (к) и д (к) очень малы при больших к).

Может встретиться н другая разновидность дискретизации: «периодическая неравномерная выборка» (Фримен [1965], стр. 77),

когда

моменты

выборки

имеют вид tj -j- иЛ,

/ = 1 , .

. ., N, п =

= 0,

± 1, . . .,

где 0 ^

tj < А. Такого рода

выборки

возникают,

например, в связи с некоторыми типами рынков, которые бывают открыты в какие-либо N нерегулярно выбранных дней в течение торгового периода (педелн, месяца, года) ѵ). Мы можем найти тогда Г (пД + tj th); j, к — 1, . . ., N. Достаточно рассмотреть лишь множество с¥ таких пар (/, к), для которых разности (tj th) отли­ чаются друг от друга на величины, не кратные А. Любая пара (tn, t0), для которой tn t0 отличается от разности, соответствующей какойлибо паре из <5°, на кА., очевидно, не несет новой информации. Зану­ меруем пары (/, к) из сf индексом т, пробегающим значения от 1

до М, и

обозначим через р т значение т-й разности tj th. Легко

видеть, что М ^

N', рассматривая для простоты скалярный абсолют­

но непрерывный

случай,

получаем

 

 

 

 

Я /Д

 

оо

 

 

 

У (пА +

рт) =

j

ЛпЛ*

2

ехР [йЛп ( к

) ] / ( к -1-

J dk,

 

-

Я /Д

U— -

оо

 

 

т = 1, . . ., М.

г) Примером может служить заключение фыочуреных сделок по шерсти

в Австралии.

4. Линейные фильтры

63

Чтобы разобраться в возникающем здесь наложении частот, пред­ положим, что / (А) равна нулю вне интервала [—МпІА, МпІА]. Тогда мы можем наблюдать функции

[М /2]

2 ехр \jprn ( А + - ^ р ) ] / ( х + ^ г ) » т = 1 , .... м,

[1/2(М +1)]

\

Д

7»

где [х] — наибольшее целое,

не превосходящее х. Это приводит

к системе из М уравнений

для

М величин / (А + 2пи/А),

и = — [Ѵ2 (М + 1)], . . ., [ѴгМ], с определителем, который на пере­ сечении т-й строки и н-го столбца имеет элемент ехр {ірт (А -j- 2пи/А)}. Легко видеть, что такой определитель равен выражению

мм

П еірт я Д

{£***.!* —

7П— 1

і

которое по предположению не равно нулю. Таким образом, мы оказываемся в той же ситуации, как если бы наблюдения производи­ лись через равные промежутки Д/М, хотя фактически точки наблю­ дений отстоят друг от друга в среднем на ДIN. Если t} расположены очень нерегулярно, то можно ожидать, что М велико по сравнению с N, и этот способ выборки в значительно большей мере способствует устранению наложения частот, чем чисто периодическая выборка. Это дает интуитивное объяснение явлению, которое было обнаружено выше для выборки, соответствующей пуассоновскому процессу.

4. ЛИНЕЙНЫЕ ФИЛЬТРЫ

Прежде чем приступить к изучению линейных фильтров, мы хотим ввести некоторые понятия, связанные с определением интегралов по неограниченному интервалу. Более полно эти понятия обсз^ж- даются в математическом приложении. Так, в скалярном случае

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

Фурье)

оо

 

j a(t)eitxdt.

(4-1)

— сю

 

Если I а (і)і интегрируема, то его можно рассматривать как инте­ грал Лебега. Однако для нас более важен другой случай (не охваты­ вающий полностью только что отмеченный), когда (4.1) определяется как предел при Т -*■ оо в среднеквадратичном с весом F (dX) (или,, коротко, (Е)-среднеквадратичный предел) интегралов

G4 Гл. I I . Спектральная теория векторных процессов

существование которых предполагается для всех 0 ^ Г < со. Таким ■образом, мы рассматриваем случай, когда функция от А, определяемая

выражением (4.1),

 

удовлетворяет

соотношению

 

со

 

СО

Т

 

 

Н т

(

 

j

a(t)eit}*cU— I

a(t)eil%dt

2 F (dK) ^ 0.

Т-*-со J

—со

—Т

 

—со

 

Это выполняется тогда и только тогда, когда

 

 

гс

 

т

 

s

 

lim

f

I [a (t) ei a dt

f a (t) eia dt

2 F (dK) = 0.

s, r-rc

J

I

-L

.L

 

 

-rc

 

-T

 

 

Аналогично можно определить интегралы типа

оо

^eii%m (dt)

СО

(преобразование Фурье — Стилтьеса), где т — функция, имеющая ограниченную вариацию на любом конечном подинтервале (— оо, оо). В векторном случае поступим следующим образом. Сначала введем обозначеипе

X (А) F (dK) X* (К)

для матрицы интегралов, типичный элемент которой имеет вид

ОО

2 ( xih(K)^T(K)Fhl(dk).

h , I — оо

Мы рассматриваем только тот случай, когда F (К) имеет эрмитовы неотрицательные приращения; будем говорить, что X (К) (^-квадра­ тично интегрируема (или принадлежит Ьг (F)), если

 

СО

 

 

Tr j

j X (K)F (dX) X* (A)} < oo.

 

 

— CO

 

 

'Теперь определим

 

 

 

оо

 

со

 

j

A(t )e ^dt ,

j eiim (d t) ,

(4.2)

— oo

 

— oo

 

шодобло тому как это было сделано выше. Именно, если для

тs

XS, T(X)= [ A ( t ) e ^ d t — \ A ( t ) e i a dt -т -s

4. Линейные фильтры

65

имеет место соотношение

lim Тг ^ X s, т(Ä.) F (dty

X(X)j- = 0.

S, Г-оо

 

то мы будем говорить, что первый из интегралов (4.2) есть (^-средне­ квадратичный предел выражений

т

j A{i) еІІХ dt.

Второй интеграл (4.2) определяется аналогично.

Всюду в этом параграфе мы используем именно такие определе­ ния преобразований Фурье (или Фурье — Стилтьеса) функций а (t), А {t), m (t) и М (t) и рассматриваем только такие функции, для которых преобразование Фурье (Фурье — Стилтьеса) в этом смысле существует. В этом параграфе процесс х (£) стационарный второго порядка.

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

(і) Иптегральные операторы

Рассмотрим сначала операцию, которая преобразует х (t) в

СО

 

y{t)= j

A(s)x(t — s)ds,

где функция А (s) такова, что интеграл

 

со

В Д =

J A (t )e ^d t

 

—оо

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

 

СО

У (t) — j M(ds)x(t —s),

где у (t) — вектор с компонентами

 

СО

2

j xh(t —s)mjh(ds),

k

—ОО

оЭ. Хеннан

66

Гл. I I .

Спектральная

теория

векторных процессов

а

М (s) — матрица

функций,

имеющих ограниченную вариацию'

ыа каждом конечном интервале,

такая,

что

 

 

 

СО

 

/і(А) = j e^M(ds)

— со

существует в смысле определенной выше среднеквадратичной сходимости относительно F.

Т е о р е м а 6. Для того чтобы случайный процесс у (t) был корректно определен как среднеквадратичный предел интегралов по конечным промежуткам, необходимо и достаточно, чтобы h (X) суще­ ствовала как (Р)-среднеквадратичный предел. Если это выполнено, то у (t) стационарен и имеет ковариационную функцию

СО

Г „(0 =

j

eu4i(X)F(dX)h(X)*.

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

Типичная компонента вектора y ( t )

имеет вид

оо

 

 

 

2

f

xh(t — s)mjk(ds),

k

— со

 

поэтому мы должны рассмотреть среднеквадратичную сходимость интегралов

ь

j хц (t s) mjh (ds)

—a

при a, b ->■ оо. Согласно теореме 3 гл. I, необходимым и достаточным условием среднеквадратичной сходимости является существование интеграла

оо

j j yk (и ѵ) mJk (du) mjk (dv);

но

 

b

CO

ь

 

j j" ук (п — V) mjh (du) mjh (dv) =

j

j

j

e^u- v^Fk (dX) mjk (du) mjk (dv) =--

 

—a

—oo

 

CO

b

eiuhnjh (du) Fh (dk),

=

j

I

J

так что для среднеквадратичной сходимости необходимо и достаточ­ но, чтобы функция hjk (А.) существовала как (^-среднеквадратичный

4. Линейные фильтры

67

предел. Это условие выполняется для всех ]', к тогда и только тогда, когда второй интеграл (4.2) существует как (.^-среднеквадратичный предел. Ковариационная функция равна

 

 

 

 

СО

 

Гу (t — s) = Е (у (s) у (і)') =

^ j М (du) Е (х (s и) х' (t — v)) M ’ (dv) =

 

 

 

 

— o o

 

 

CO

 

CO

 

 

=

j j

M (du)

j

еЫ-Ч-ѵ-МР (dk) M' (dv) =

 

— ОЭ

 

— CO

 

 

 

o o

 

 

 

 

=

j

 

(к) F (dk) h (k)*,

 

 

— CO

 

 

 

 

откуда видно,

что у (t) — стационарная

последовательность, имею­

щая указанный спектр, и теорема доказана.

Таким

образом,

действие фильтра

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

в спектральных терминах. Это в особенности относится к скалярному

абсолютно непрерывному случаю, когда / (А,)

заменяется просто

на

I h (к) |2/(Х).

В векторном

абсолютно

непрерывном случае f

(к)

переходит в

h (к) f (к) h (к)*.

Ситуация,

в

которой

 

оо

y(t) = 2 Ahx(t — tk),

— с о

охватывается случаем, когда М (t) возрастает только скачком, так как А h можно считать величиной скачка функции М (t) в точке th. В случае когда М (t) абсолютно непрерывна н М (ds) можно записать как A(s)ds, функцию Л (s) иногда называют матричной импульсной переходной функцией фильтра, ибо если вместо х (t) в формуле для

у (t) подставить импульсную функцию (т. е. дельта-функцию с осо­ бенностью в начале координат; см. § 3 математического приложения), то на выходе будет именно А (і). Функция h (к) называется матрич­ ной частотной характеристикой фильтра. Смысл термина «частот­ ная характеристика» станет более попятным, если рассмотреть спек­

тральное

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

у (t). В

общем случае

имеем х)

 

СО

СО

с о

с о

 

у (I) =

j

М (ds) j e-i(f-s)X.z (dk) =

^ e~itl I

j eis%M (ds) j

z (dk) =

— CO

— CO

— OO

— CO

 

 

 

 

00

 

 

=j e-il4(k)z(dk).1

1)Полное обоснование этих преобразований мы даем в приложении к этой главе. (См. доказательство теоремы 9.)

5*

68

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

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

(іі) Дифференциальные операторы Рассмотрим теперь операторы вида

о

где Ak — матрицы. Предположим, что

СО

 

(4.3)

— со

 

Положим теперь

 

р

(4,4)

}нХ) = ^ А ІІ( - і Х ) к.

Ü

Т е о р е м а 7. Для того чтобы х (t) был р раз среднеквадратично дифференцируем, необходимо и достаточно, чтобы выполнялось условие (4.3). В этом случае процесс

О

является стационарным и имеет ковариационную функцию

со

Г„(*)= I е>Чг{Х)В((1Х)!1(Х)*

оо

где h (X) определяется соотношением (4.4).

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

lim У(62— 6])— у ( — б,) — ѵ(6г)-И(0)

öl, Ö2-»-0

6j62

 

00

=

1б2)-1 {<?І(02-б,)Х_е-іб1?у_ еібЛ+ typ (dX) =

 

со

=

ЬД (е~і6іХ — і) б; 1 ів2к— 1 ) F (dl).

00

4. Линейные фильтры

69

Подинтегралытоо выражение оценивается по абсолютной величине как

4

sin 1/2 ÖiX-sin !/2 62Х

< я 2,

 

6j62

 

поэтому, в силу мажорированной сходимости1), условие является достаточным. С другой стороны, полагая 6j = ö2, получаем

lim

j 46-2 sin2 і/26ЯF (dX) = lim

J

WF (dX) >

а-о

 

 

 

 

> > ! “

J w

№ w = J

 

 

—а

—а

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

Если теперь х (t) среднеквадратично дифференцируем р раз, то, согласно теореме 2 гл. I,

V

 

E(y(s)y(ty)=

^

A{k ) E { ( ^ x ( S)) ( J L x{i) ) ' ) A ( iy =

 

 

 

 

 

О

 

 

 

 

 

 

 

p

 

 

 

 

 

 

=

2

А ^ і ^

ж T(t- s)A(iy .

Это равно

 

 

 

ft, г=о

 

 

 

 

 

 

 

СО

 

Р

 

 

 

 

со

j

ei(t-s)x

2

A{k)(iX)l { - iX)kF{dX)A{l)' = j e ^ - s^h{X)F{dX)h{X)*.

— со

h,

1= О

 

 

 

 

—со

Появление

множителя

( — i)k в

члене ( — i)hXh можно объяснить,

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

случае2)

 

 

 

 

 

 

 

 

 

со

р

 

 

со

 

У(і)=

J

 

 

 

je~^h(X)z(dX

 

 

 

— СО

О

 

 

—со

1)Здесь и далее имеется в виду известная теорема Лебега о возможности продельного перехода под знаком интеграла в случае, когда последовательность подинтегральных функций мажорируется фиксированной интегрируемой функ­ цией.— Прм. персе.

2)Подробное обоснование спектрального представления мы даем в прило­ жении к этой главе (см. доказательство теоремы 9).

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