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

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

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

40

Гл.

I . Вводные сведения

 

 

Рассмотрим теперь функции ср (?) ішда

 

 

Г

0,

? iß*

 

 

ср (С) =

<ц>.

/ = 1.

• ••>

(3)

 

0,

?jy<?,

 

 

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

 

 

 

 

 

 

N

 

 

cp (г) 4_>cp= j Ф ( < Ш < ? 0 = 2

(4)

 

 

1

 

 

В левой части мы рассматриваем ср (г) как элементы полного линейного про­ странства L 2 (/■') функций, квадратично интегрируемых относительно Г (dt),

тогда как в правой мы имеем случайные величины, принадлежащие гильбер­

тову пространству <й?,

определяемому

процессом Е (/). Конечно,

Ь 2 (F) само

является гильбертовым

пространством, и, как легко проверить,

 

 

j Фі (0 Фг (0

F {dt) = (Фі, фз)

(5)

для любых двух функций вида (3).

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

относительно F, ибо если ср„

(?) — последовательность функций указа иного

типа, такая, что

 

 

lim

f

I фт (0 — фп(0 I2 F {dt)— 0,

т> п -+ со

J

 

то фт (?) сходится по норме L->

{Г) к функции ф (?), квадратично интегрируемой

относительно F, а (5) гарантирует, что соответствующая последовательность

Ф,п сходится к элементу пространства &(!. Кроме того, ступенчатые функции (3) «плотны» в Ь 2 (F) в том смысле, что все квадратично пптегрируемые функции

могут быть получены из сходящихся в среднеквадратичном последовательностей таких фупкцпй. (На самом деле мы можем считать, например, что точки разрывов tj рациональны, так что (/’) содержит счетное плотное множество, т. е. L 2 (F)

«сепарабельпо».) Таким образом, соответствие (4) позволяет для любой ф (?),

принадлежащей

Ь г (/•’), определить случайную

величину, которую

мы обозна­

чим символом

 

 

 

 

 

 

 

 

 

(?) S

(dt).

 

 

 

Такие же

рассуждения

позволяют

определить

интеграл

 

 

 

j ® ( ? ) i( d « ) ,

 

 

(6)

где I (?) теперь — векторный

процесс с

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

Е (II «а) - Е (s2)l ІЕ (П) -

Е (*і)]*) =

0,

і* > s2 > *і >

Sl.

Здесь, конечно, Ф (?) — матрица, у которой число столбцов совпадает с размер­

ностью Е, а элемент

(?)

матрицы

Ф (?)

должен

принадлежать L z (Fj), где

Fj

(t) -

Fj (s) = E

(I Ej

(l) - Ej

(s) I2}.

В самом деле, (6) определяется как вектор с компонентами

5 [ ФП(0 Er И ) .

о

а каждое отдельное слагаемое уже было определено.

Приложение

41

 

 

 

 

 

 

2.

Доказательство теоремы 4'

 

 

 

 

 

 

 

 

 

В §

3 мы ввели м одел ь (3 .1 0 ),

к отор ая пр ивела н ас к р ассм отр ен и ю

реш ен ий

систем ы у р а в н ен и й

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[2В(/)*9-і]

Ь = 0.

 

 

 

 

 

 

 

 

 

 

(1)

П усть. z a — кор ен ь у р а в н ен и я

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d e t[2 ß (/)4 ^ ] = 0;

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обозн ач и м

соответствую щ и е

р еш ен и я

систем ы

(1)

ч ер ез

b

(г, и ) ,

г =

1, . .

 

s u

С ущ ествован и е т а к и х

р еш ен и й

и и х к р атн ости

к ак собств ен н ы х

век тор ов

у с т а ­

навл и ваю тся т ео р и ей эк ви вал ен тн ости м атр и ц ,

элем енты к отор ы х я в л я ю тся

м н о ­

гоч лен ам и от п ер ем ен н ой z

(см .,

н ап р и м ер ,

 

М акдаф ф и

[1 9 5 6 ], стр .

4 0 ,

п л и

Б о х е р

[1 9 0 7 ],

гл .

X X ) .

Д о к а зы в а ется ,

что для

м атрицы

в л ев о й

части

(1), к о то р у ю

мы

о бозн ач и м

С (г),

су щ ест в у ю т м атрицы

Р

и

 

Q ,

элем енты

к отор ы х

я в л я ю тся

м н о ­

гоч лен ам и

от

z,

с о п р едел и тел я м и ,

равны м и

н ен ул евы м к он стан там ,

так и е,

что

 

 

 

 

 

 

 

 

 

 

Г>4 (*)

к 2 (z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P C

( z ) Q =

 

 

 

 

 

 

 

 

 

 

 

 

 

= #,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

h p

(z) J

 

 

 

 

 

 

 

 

где

н еди агон ал ьн ы е элем ен ты

равны

н у л ю ,

 

(z)

— м н огоч л ен

от z

степ ен и

p it

п р ич ем

Аг д ел и т

Лі+І

дл я

л ю бого

і.

(В общ ем

слушав

в н п зу на

д и агон ал и

м огут

бы ть н ул ев ы е

элем ен ты ,

н о

д л я С

(г)

это

н ев о зм о ж н о ,

поскольку7 мы п р ед п о л о ­

ж и л и ,

что

В

(0)

— еди н и ч н ая

м атр и ц а .)

М ногочлены

h j

(z) н азы ваю тся

и н в а ­

риантны м и

м н ож и тел я м и .

П у ст ь

z u

к о р ен ь Іц (z),

н е

явл я ю щ и й ся

 

кор н ем

/гг _!

(г);

тогд а , уч иты вая ,

что Р - 1

— т а к ж е м атри ц а м н огочлен ов от z,

мы види м ,

что

су щ ест в у ет s u =

( р г - f

1)

реш ен и й

 

Ъ (г,

и ) ,

і =

1, . .

.,

su ,

ур ав н ен и я

С {zu ) b =

0;

и м ен н о ,

 

п одставим зн ач ен и е

z

z u

в Р ,

С ж Q ж возьм ем

в

к ач е­

стве

b (г, и )

п о сл ед н и е

( р

і

+

1)

стол бц ов

п ол уч аю щ ей ся

м атрицы

 

Q .

Эти

р еш ен и я ,

оч ев и дн о ,

ли н ей н о

 

незав и си м ы ,

поскольку7

оп р ед ел и тел ь

Q

 

равен

н ен у л ев о й

п остоя н н ой .

Т ак и м

о б р а зо м ,

каж дом у7

гц мы со п о ста в и л и s u

реш ен и й

у р а в н ен и я

С (zu ) b =

 

0 ,

а

и м ен н о ,

ук азан н ы е

вы ш е

Ъ (г, и).

П у ст ь

р

„ —

степ ен ь м н ож и тел я

(г — z u )

 

в

и н вар и ан тн ом

м н о ж и тел е ,

соответствую щ ем

Ъ (і,

и ) .

П остр ои м тогда

си стем у

векторов :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ь (і, u ) z £ a i ,

 

и = 1,

 

. . . , s,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

І

= 1,

 

. . •,

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7 = 0,

. . . ,

і . ц — 1),

 

 

 

 

 

 

 

 

где

s —

ч и сло

р азл и ч н ы х к о р н ей

z u .

О тм етим ,

что ч и сло

т а к и х

векторов

со в п а ­

д ает

со

степенью

p q оп р ед ел и тел я

С

(z)

к ак

м н огоч л ен а

от z.

О бр а зу ем

теп ер ь

вы р аж ен и е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с (п)= У] 2

2

с (*• й “) zunib (й “)>

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

и

г

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

с (і, ) ,

и )

дол ж н ы

быть

вы браны так ,

 

чтобы этот

вектор

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

начальны м

зн ач ен и я м

 

п р и

и

=

— 1,

. . .,

 

q. В ек тор

с

( п ) , оч еви дн о ,

у д о в л е ­

тв ор я ет

ур ав н ен и ю 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2О В (/) с (п~ 7)= °-

(4>

42

Гл. / . Вводные сведения

(То, что z^n’b (і,

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

(4) при j

> 0, можно проверить, дифферен­

цируя і раз обе части соотношения

С (z) Q =

Р~гП п полагая z = zu .) Остается

только доказать существование чисел с (і, j,

и) для заданных начальных усло­

вий, а это сводится к проверке того, что матрица, у которой в к-й группе из р строк (А- = 1, . . q) столбец с индексом (і, /, и) равен

zuh( ~ k)ib (*.

имеет ненулевой определитель.

Но если бы этот определитель был равен нулю,

то можно было бы найти ненулевые с (/, /, и),

такие, что с (л) =

0, при п =

= —1, . . ., —q, а это, очевидно, невозможно,

так как компоненты

вектора

с (п) являются многочленами степени не выше (q — 1).

 

 

Теперь мы можем завершить доказательство теоремы 4'. Пусть

(J~i(z_1) =

— ^ s Q j z ~ i u P ( z ~ l) — y^ ! P]z~i.

Тогда, переходя к новому процессу

і/ (л) ==

= 2 Qix (n — /) 11 заменяя правую часть (3.10)

выражением

 

 

У і Р} ^ і А ( к ) г ( п - к - і ) ,

 

 

j

ft

 

 

 

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

вторая часть доказательства, и мы предположим теперь, что ни один из корней уравнения (2) не лежит на единичной окружности. Тогда для ( ^ В (/) z1)”1

имеет место разложение Лорана в кольце, содержащем единичную окружность; записывая его в виде

СО

00

мы приходим к выражению

со

2Р (п и) А (к) е {и к).

и= — со

Это выражение определено корректно как предел в среднеквадратичном, иосколь ку элементы F (п) экслопепцнально стремятся к пулю при возрастании и; но

тогда, в силу соотношения

 

 

/ ) ^ öo7p-

q

со

 

É s x/){

2

и — /) А (к) в к)} = А (к) в (п — к)

0

j = — сю

 

u

 

 

 

S

о о

22 F (n~ и) А (к) в (а— к)

к= 0 и — — оо

является решением (3.10). Полагая

S

л (;■)= 2 р (І—к)А (к'і< ft= 0

мы можем записать это решение в виде

оо

 

2 Л(/)е(я —/).

(5)

—ОО

 

Комбинируя (3) и (5), получаем теорему 4'.

Г л а в а II

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

1. ВВЕДЕНИЕ

Эта глава посвящена анализу Фурье для временных рядов. Методы Фурье тесно связаны с понятием стационарности, так что случайные процессы, обладающие этим свойством, будут находиться в центре внимания, хотя мы рассмотрим также и некоторые отклонения от стационарности. Пожалуй, сначала стоит пояснить, почему методы Фурье играют такую существенную роль в рассматриваемой теории. Наложенное нами фундаментальное условие стационарности для скалярных процессов имеет вид у (s, t) = у (t s). Таким образом, мы рассматриваем ограниченный класс ковариационных функций, которые инвариантны по отношению к группе сдвигов вещественной прямой, т. е. у (s, t) = у (s+ т, t + т). В связи с этим возникает вопрос: нельзя ли заменить функцию времени х (t) какой-нибудь

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

диагональную, форму: у (s, it) = 0, s =£t. (Это может показаться

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

00

(1.1)

о

справедливое для некоторого специального класса стационарных процессов. Здесь процесс х (t) выражен линейно через новый слу­ чайный процесс £ Ц), хотя и не стационарный, но имеющий стацио­ нарные и, более того, некоррелированные приращения; таким обра­ зом, рассуждая очень нестрого, можно рассматривать £ (dt)/dt как

i (it).

Однако представление (1.1) не дает того, что нам нужно, так как оно не может быть реализовано, если нам известно лишь, что у (s, t) = = у (t — s), но функция у (t) в целом не известна. Мы же ищем преобразование, которое диагонализирует у (s, t) и задано a priori, если только известно, что у (s, t) = у (t s) (но сама функция у (t)

44

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

не известна). Не удивительно, что такое преобразование выполняется методами Фурье, ибо в некотором смысле, который мы здесь не уточняем 1), exp it% является собственной функцией оператора U (т), действующего по формуле U (т) / (I) = / (t -f- т) и, следовательно описывающего стационарность процесса.

Прежде чем продолжить обсуждение намеченной процедуры, мы рассмотрим пример, в некоторых отношениях более простой и более общий, который иллюстрирует сказанное выше. Рассмотрим непре­ рывный в среднеквадратичном скалярный случайный процесс х (t), ковариационная функция которого у (s, t) известна, и ограничимся сначала случаем интервала [О, Г]. Тогда по теореме Мерсера (Рисе и Надь [1956], стр. 265)

со

у {s, 0= 2 н-іФі(Оф;(0.

U

где ф; (О — собственные функции ядра y(s, О,

Т

j V (s, 0 Фі (i)* = WPi(s),

J Фі (s) Фj (s) ds = 6].

и

Двойной ряд для у (s, 0

сходится к своему пределу равномерно

по s и I. Образуем теперь ряд

СО

Т

00

2 ф/ (о J х (о Ф< ( о ds= 2 а /Фі (О-

о

и

о

Коэффициенты а г являются

случайными величинами (определение

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

в п. 1 приложения к гл. I).

Они удовлетворяют соотношениям

Е (а-іо.]) =

г

E{x(s)z (0) Ф; (s) q>j (t) ds dl —

j j

 

о

 

 

 

T

 

 

=

[ j

У(s, 0 Фі О) Фі (0 ds dt =

M i •

 

о

 

 

Таким образом, они «ортогональны». Кроме того,

 

 

 

т

 

Е (х (/) а,-) =

j у (s, I) фг (s) ds =

(О,

 

 

о

 

х) См. § 10.

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

45

откуда

следует, что

 

 

N

N

 

Е [{я (*)— 2«іФі (О)2] =У (t, t)— S М-гФГ (^).

 

Ü

о

.а это выражение, как нам известно,

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

Итак,

имеем

 

оо

*(0 = 2о аіЧ>і (і)

всмысле среднеквадратичной сходимости частичных сумм1). Этот результат остается верным, например, в случае, когда [0, Т] заме­

няется на (— оо, оо), у (s, t) квадратично интегрируема на плоскости и интегралы по [0, Т] заменяются на интегралы по (—оо, оо). Двой­ ной ряд для у (s, t) будет сходиться теперь только в среднеквадра­ тичном, т. е.

ооN

lim f

^

[ y (s, t)— 2

PiCPi(s)cpi(i!)]2rfs*

= 0>

N-+ ОО V

со

*-

-I

 

-

0

 

 

хотя, если у (s, г) непрерывна, сходимость будет

равномерной на

любом конечном

интервале.

 

 

Мы получили здесь некоторую диагоиализацшо ковариационной матрицы (новые переменные — это величины ссг), которая имеет, однако, ограниченную применимость, поскольку требует полного знания у (s, t). Некоторое обобщение этих пдей читатель может

найти

в работе Парзепа [1961].

В §

2 мы сформулируем и докажем основные теоремы, касающие­

ся анализа Фурье для х (t). В § 3 мы обсудим связь с аналогичным анализом Фурье для временных рядов с дискретным временем. В § 4 обсуждаются линейные фильтры; некоторые специальные моде­ ли рассматриваются в § 5. Нелинейные фильтры и некоторые формы отклонений от стационарности, которые все еще позволяют исполь­ зовать спектральные методы, рассматриваются в § 61 1 , где также затрагиваются другие спектральные теории.

2.СПЕКТРАЛЬНЫЕ ТЕОРЕМЫ ДЛЯ СТАЦИОНАРНЫХ ПРОЦЕССОВ

СНЕПРЕРЫВНЫМ ВРЕМЕНЕМ

Докажем сначала следующую теорему.

Т е о р е м а 1. Если Г (s, t) = Г (ts) ковариационная матрица стационарного процесса второго порядка, то

оо со

 

Г (0 = j

e^F (dX) = j

{cos tXC (dX) + sin tXQ (dX)},

(2.1)

 

—CO

0

 

 

Э

Полученное

разложение

называется разложением Каруиена—Лоэва.

Прим, перев.

46

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

где F (А) — матрица, приращения которой F (А,,) — F (Х2), Аі ^ Х2г эрмитовы и неотрицательны. Функция F (X) определяется однозначно,

если мы дополнительно потребуем, что (і) lim F\X) = 0, (іі) F (А)

Я->—OO

непрерывна справа. Матрицы С (А) и Q (А) вещественны и являются соответственно симметричной и кососимметричной, причем С (Aj) —

С (А2) неотрицательна при Aj ^ А2.

Матрица F (А)

называется матричной спектральной функцией,

а С (А) іі Q (А) —

соответственно коспектральной матричной функ­

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

При доказательстве теоремы 1 мы используем классические «теоремы единственности и непрерывности для характеристических функций» (см., например, Крамер [1946]), утверждающие, что после­ довательность функций распределения сходится к собственной функ­ ции распределения тогда и только тогда, когда соответствующая последовательность характеристических функций сходится поточеч­ но к функции, непрерывной в пуле, которая оказывается характе­ ристической функцией предельного распределения, и что эта харак­ теристическая функция однозначно определяет функцию распреде­ ления, нормированную, например, условием непрерывности справа. Мы будем использовать эти теоремы в случаях, когда верхний предел всех функций распределения равен одному и тому же положительно­ му числу (не обязательно равному единице); очевидно, что указан­ ные теоремы применимы и в этом случае.

Чтобы доказать теорему 1, построим сначала скалярный процесс

ха (t) = а*х (I),

где а — фиксированный вектор из' р комплексных постоянных. Положим

Уа (0 = Е (ха (s) ха (s+ і)) = а*Г (t) а.

Тогда имеет место основная

Л е м м а 1. Функция уа (t) неотрицательно определенна в том смысле, что для любых п, . . ., tn и комплексных постоянных

ß„ • • К

У

0.

( 2. 2)

і =1/1=1

 

 

 

Это следует из простых преобразований

У У

ßjß/,ya ( t j - h ) = y y РАЕ(Яа(tk) УУЩ )) =

і = 1=1

з =1/і== і

 

= E (|S ß ;* « (4 )|2)> 0 -

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

47

Для продолжения доказательства теоремы 1 положим

•у(Т) (/) _

1(і

V")Ѵа (0 >

 

 

ІО,

\ t \ >T,

 

 

СО

 

fa\X) =

j e - ^ ) ( t ) d l .

Поскольку

 

 

 

T

 

 

T

f P ( X ) = j

 

ф - ) T« ( 0 * = Y j J

-T

 

 

0

то /<Р(Я )> 0. ибо

правая

часть может быть с любой точностью

аппроксимирована

суммами

вида (2 .2) (где ß7= exp (itjX) 8t, а tj

равноотстоящие точки). Образуем теперь чезаровское среднее функ­

ции /£°(Я):

 

м

 

■ш ( О - Т т К " » ) « ' “ “»-'

(2-3)

 

которое сходится к (t) при М — оо. (См. § 3 математического приложения.) Так как (2.3) является характеристической функцией

(подинтегральный множитель перед exp ИХ положителен), то у(<Р (t) как предел сходящейся последовательности характеристических функций, непрерывный при t = О, также является характеристиче­ ской функцией. Теперь, устремляя Т к бесконечности, мы видим, что функция уа (t) также является пределом последовательности характеристических функций, и так как она непрерывна при t = О, то является характеристической функцией:

СО

Та 00= j eiixFa {dX).

— оо

Единственность Fa (Я) при выполнении условий (і), (іі) теоремы 1 следует из теоремы единственности для характеристических функций.

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

СО

Тj(t)= j e^Fj(dl).

СО

Если мы положим

а\х(і> (t)=xj (t) + xh (0 , xam (t)=x} (t) + ixh 00,

то, учитывая, что

у (Tad) (0 — Ту00— Тл 001 + Y ГТа<2>00 — Т.7 (0 — Та(01 = Т;. ft00.

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

имеем

 

 

со

 

V j . k ( t ) =

j e ^ F j h(dX),

 

где

 

 

 

Fjk (Я) I [^аш ( \ ) - F j (X)~FU(Я)] + 1 [^e(8) (Я) - F J (Я)-

Fh (Я)].

Обозначим через F (Я) матрицы

с элементами Fjh(Я). Из

того что

(*) = Ѵ/и (— 0. следует,

что F (Х) = F* (X), ибо

 

СО

оо

 

 

[ e^F*{dX)=:

J e-'^F ' (йЯ) = Г' (—і) = Г(*).

 

— со

— 00

 

 

Из доказательства что F (Я) эрмитова и ния. Так как

видно,

что a*F (X) а = Fa (X), откуда

следует,

имеет неотрицательно определенные

прираще­

СО

00

 

Г (

- 1) = j e-«\F (dX) = Г' (/) =

j

e ^ F ' (dX),

 

 

— СО

 

— со

 

 

то мы видим,

что F' ( — Я) =

Г(0) — F (Я)

во

всех точках непрерыв­

ности. или. что равносильно,

 

 

 

 

F(Xt) - F ( X 2) = F' ( ~ X 2- ) ~ F ' ( - X i - ) , Я ^ Я 2,

(2.4)

где. например,

F'( — Я!— )=

lim F' (Я). Мы будем сокращенно запи-

сывать (2.4) в виде F (dX) = F ’ (d(—X)), что равно, конечно, F(d{—Я)). Используя обозначения Re для вещественной части, Im для мни­

мой части, положим

f

2Re {Fjk {dX)},

Я > 0 ,

Cjh{dl)- \

Fjh(dX),

Я=0,

(^Я) = -

2Іхп {Fя, (dX)},

Я > 0.

Эти соотношения определяют две вещественные функции ограни­ ченной вариации Cjh (Я) и (?j7i (Я). Как и в § 1, гл. I, мы опускаем один из индексов, если эти индексы совпадают. Тогда Qj (Я), очевид­ но, равны нулю, а Cj (Я) имеют неотрицательные приращения. Соот­ ношения F (Я) = F* (Я), F (dX) — F' (d (—Я)) показывают, что С (dX) — симметричная матрица и ее можно определить так, чтобы она была четной функцией от Я, тогда как Q (dX) — кососимметрич­ ная матрица и ее можно определить так, чтобы опа была иечетиой функцией от Я. Имеем

ОО

СО

оо

T(t) = j j cos iXC (dX) -f Y

j sin tXC (dX) —^

j cos tXQ (dX) +

оо

со

те

+ ' j sin tXQ (dX)

j cos IXC (dX) — у j sin tXC (dX) -f

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

49

 

ОО

СО

 

 

f y

^ cos IkQ {dk) + у

j s*n

(^ )

 

 

о

0

 

 

oo

CO

 

 

 

^ cos tkC {dk) -|- j sin ikQ {dk).

о

о

 

 

Это завершает доказательство теоремы 1.

 

Для матрицы

F (Я) имеет место

разложение Лебега

 

F (К) = Fll>(k) + F<2>(к) + Б<3>(к).

 

Здесь F^^ {к) абсолютно непрерывна

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

меры Лебега

на прямой), так

что

 

 

 

— ОО

 

 

где / (к) — матрица с элементами

к (к), называемая матричной

спектральной плотностью. Пусть / (к) = Ѵ2 (с (Я.) —

iq (Я)); тогда

с (Я) и q (Я) называются соответственно матричными коспектралъной

и квадратурной спектральной плотностями. Функции h (Я) для ] ф к называются взаимными спектральными плотностями. Компо­ нента FW (Я) может возрастать только скачками в конечном или счетном множестве точек интервала (— оо, оо), не имеющем предель­ ных точек (помимо ±оо).

Элементы матрицы F ^ (Я) непрерывны и имеют почти всюду равные нулю производные относительно лебеговской меры. Было бы трудно придать им какой-либо физический смысл, и мы часто будем предполагать, что эта третья компонента отсутствует. На самом деле, как мы увидим дальше, в случаях, когда стационарная модель адекватно отражает всю сложность реальной ситуации, часто (но, возможно, не всегда) можно считать, что F(1) (Я) является единствен­ ной присутствующей компонентой. Ниже мы обсудим это более подробно.

Теперь мы должны пояснить физический смысл утверждения теоремы 1. Для этого рассмотрим сначала несколько искусственный

частный случай. Предположим, что Xj (/) = xj {t -f- 2/гл), / = 1, . . .

. . ., p, k = 0, ± 1 , . . ., T .

e. что наблюдается периодическая слу­

чайная функция. Например,

Xj (t) при —я < t ^ я может соответ­

ствовать некоторому измерению, производимому иа некоторой фикси­ рованной широте и долготе t, причем результаты измерений про­

должаются

просто по периодичности. Тогда также

yjh (t) =

= Vjh (t +

2kn). Мы будем писать х} (ф) и yjk (ф)

в случае,

когда

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

<С ф ^

я.

 

4 Э. Хеннан

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