книги из ГПНТБ / Хеннан, Э. Многомерные временные ряды
.pdf40 |
Гл. |
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. Нелинейные фильтры и некоторые формы отклонений от стационарности, которые все еще позволяют исполь зовать спектральные методы, рассматриваются в § 6—1 1 , где также затрагиваются другие спектральные теории.
2.СПЕКТРАЛЬНЫЕ ТЕОРЕМЫ ДЛЯ СТАЦИОНАРНЫХ ПРОЦЕССОВ
СНЕПРЕРЫВНЫМ ВРЕМЕНЕМ
Докажем сначала следующую теорему.
Т е о р е м а 1. Если Г (s, t) = Г (t—s) — ковариационная матрица стационарного процесса второго порядка, то
оо со
|
Г (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 Э. Хеннан
