/
Текст
ТЕОРИЯ ВЕРОЯТНОСТЕЙ
МАТЕМАТИЧЕСКАЯ СТАТИСТИКА
Ю.А. РОЗАНОВ
ТЕОРИЯ
ОБНОВЛЯЮЩИХ
ПРОЦЕССОВ
ТЕОРИЯ ВЕРОЯТНОСТЕЙ
И МАТЕМАТИЧЕСКАЯ СТАТИСТИКА
Ю. А. РОЗАНОВ
ТЕОРИЯ
ОБНОВЛЯЮЩИХ
ПРОЦЕССОВ
ИЗДАТЕЛЬСТВО «НАУКА»
ГЛАВНАЯ РЕДАКЦИЯ
ФИЗИКО-МАТЕМАТИЧЕСКОЙ ЛИТЕРАТУРЫ
МОСКВА 1974
517.8
Р64
УДК 519.21
Теория обновляющих процессов. Ю. А. Роза-
нов, Главная редакция физико-математической
литературы изд-во «Наука», 1974.
В книге изучаются закономерности обновления
данных о «наблюдаемом» случайном процессе в зада-
чах линейного оценивания, прогнозирования и филь-
трации, которые приводят к проблеме факторизации
корреляционного оператора. Она рассчитана на
специалистов по теории вероятностей и функцио-
нальному анализу.
© Издательство «Наука», 1974 г.
Розанов Юрий Анатольевич
ТЕОРИЯ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ
М., 1974 г., 128 стр.
Редактор М. П. Ершов
Техн, редактор 4. П. Колесникова Корректор В. П. Сорокина
Сдано в набор 15/Х 1973 г. Подписано к печати 8/IV 1974 г. Бумага
84Х108|/з2, тип. № 1. Физ. печ. л. 4. Услов. печ. л. 6,72. Уч.-изд, л. 5,75.
Тираж 6800 экз. Т-20864. Цена книги 39 коп. Зак. № 848
Издательство «Наука»
Главная редакция физико-математической литературы
117071, Москва, В-71, Ленинский проспект, 15
Ордена Трудового Красного Знамени
Ленинградская типография № 2 имени Евгении Соколовой
Союзполиграфпрома при Государственном комитете Совета Министров
СССР по делам издательств, полиграфии и книжной торговли
198052, Ленинград, Л-52, Измайловский проспект, 29.
_20203_-_058_
053(01)-74
ОГЛАВЛЕНИЕ
От автора............................................... 4
Глава I. Обновляющие процессы и канонические пред-
ставления ............................................. 5
§ 1. Введение.................................. 5
§ 2. Структурные типы и подчиненные процессы ... 14
§ 3. Некоторые примеры обновляющих процессов ... 29
Глава II. Регулярные процессы.................. 39
§ 1. Изометричные семейства и некоторые примеры . . 39
§ 2. Некоторые модели случайных процессов. Понятие
регулярности и проблема факторизации........... 47
§ 3. Одна теорема о факторизации................... 56
Глава III. Регулярные стационарные процессы............ 66
§ 1. Структурный тип регулярного стационарного про-
цесса ............................................ 66
§ 2. Представление Вольда и факторизация спектраль-
ной плотности.................................. 71
§ 3. Кратность регулярного стационарного процесса . . 79
§ 4. Условия регулярности........................ 83
Глава IV. Эквивалентные случайные процессы.............100
§ 1. Понятие эквивалентности. Вероятностная интер-
претация в случае гауссовских распределений . . 100
§ 2. Эквивалентность стационарных процессов .... 102
§ 3. Случайные процессы, эквивалентные винеровскому
процессу....................................... . 124
ОТ АВТОРА
В этой небольшой книге излагаются новые во-
просы общей теории случайных процессов второго
порядка, в рамках которой случайный (одномерный)
процесс рассматривается как функция со значениями
в гильбертовом пространстве случайных величин,
имеющих конечный второй момент. Основные про-
блемы касаются «обновления» данных о случайном
процессе с течением времени. Типичной в этом смысле
является проблема о полной недетерминированности
(бесконечномерного) стационарного в широком смысле
процесса. Как оказалось, эти проблемы (и особенно
общий вопрос о регулярности случайного процесса)
тесно связаны с вопросом о факторизации опера-
торов в гильбертовом пространстве относительно за-
данной цепочки расширяющихся подпространств, и
в силу этого обстоятельства в книге дается также по-
дробное доказательство некоторых, сравнительно не-
давних, результатов в этой области функционального
анализа.
ГЛАВА I
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
И КАНОНИЧЕСКИЕ ПРЕДСТАВЛЕНИЯ
§ 1. Введение
Пусть £(/) = {£/ (0)Г, tQ < t < Г,—многомерный слу-
чайный процесс с компонентами h (/), i — 1, 2, ..т
(число пг может быть равным оо). Предположим, что
все значения & (/) имеют нулевые средние E£f(/) = 0
и конечные вторые моменты Е| h (t) |2 < оо. Обо-
значим замкнутую (в среднем квадратичном)
линейную оболочку всех значений g/(s), и
будем рассматривать Ht(l) как подпространство гиль-
бертова пространства всех комплекснозначных слу-
чайных величин т), Е| л |2 < оо, со скалярным произ-
ведением
01b Пг) = ЁП1П2
(не делая различия между величинами, совпадаю-
щими с вероятностью 1).
В линейных задачах теории случайных процессов
таких, например, как линейное прогнозирование и
фильтрация, соответствующее подпространство
представляет собой «совокупность данных», которыми
располагает «наблюдатель» к моменту времени t
оценивания той или иной величины т], Е|т]|2<оо
(за наилучшую оценку принимается проекция этой
величины на подпространство /7Д£)).
Обозначим Pt оператор ортогонального проектиро-
вания на и введем пространство
и
h<t<T
6
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
Семейство проекционных операторов Pt, tQ<t < Т,
в Я(£), а также семейство соответствующих подпро-
странств
растущих с течением времени /:
= при
будут главным объектом нашего исследования. Основ-
ной вопрос, который нас будет интересовать, заклю-
чается в следующем. Как охарактеризовать эволюцию
подпространств Ht(g) с течением времени /?
Чтобы пояснить общий подход к решению этого
вопроса, обратимся к одномерному случайному про-
цессу I (/), t0<t < Т, с некоррелированными
приращениями. Будем предполагать, что этот
процесс непрерывен слева:
lim g(/-A) = g(/),
h->+0
и положим F (t) = Е| £ (/) |2. Назовем монотонно не-
убывающую, непрерывную слева функцию F(f), tQ<
< t< Г, структурной функцией. Соответствующее под-
пространство //(g) состоит из всех величин т), пред-
ставимых в виде стохастического интеграла
т
п = / С (t)d£ (О,
*0
где комплекснозначная функция c(t), удо-
влетворяет условию
т
J|c(/)|2dF(0<oo.
/о
Если рассмотреть гильбертово пространство С всех
таких функций c(t), < Т, со скалярным произ-
ведением
т
§ П
ВВЕДЕНИЕ
7
и в нем подпространства Ct всех функций, обращаю-
щихся в 0 вне интервала [/0, 0> то окажется, что инте-
ресующие нас подпространства будут унитарно
изоморфны соответствующим подпространствам Ct,
поскольку состоит из всех величин вида *)
t
П= J c(«)^(s).
t
где j* | c(s) pdF(s) < оо. Грубо говоря, семейство под-
^0
пространств Ht (g), <0</< Т, устроено точно так же, как
и семейство подпространств Ct> t0<.t < Т, эволюция
которых с ростом t представляется достаточно нагляд-
ной и полностью характеризуется соответствующей
функцией F(t).
Рассматривая общий случай, мы будем предпола-
гать пространство Н(g) сепарабельным, а семей-
ство Ht(g), t0<t < Т, непрерывным слева:
lira = (1.1)
Л-»+0
Эти условия будут выполнены, например, если исход-
ный случайный процесс £ (/) = {£//0 < t < Г,
является непрерывным слева:
lim — A) = gj(f), /=1, ..., m.
Л->+0
t
*) Здесь и далее в интегралах вида J интегрирование
S
ведется по интервалу [s, /), замкнутому слева и открытому
справа. При этом
t t
j = pg(«) = g(/)-g(S)
t9 s
И
t t
j dF (s) = F (t), J dF (u) = F (t)-F (s), t0<s<t<T.
8
8
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. 1
Вопрос об эволюции подпространств Ht (£), /0< /< Г,
для произвольного случайного процесса £(/) = {£/
tQ<t<T, можно было бы считать решенным, если бы
удалось найти некоррелированные между собой про-
цессы Xj (/), j = 1, ..., М, t0< t < Tf с некоррелиро-
ванными приращениями, такие, что отвечающие много-
мерному процессу X(/) = {X/tQ<t<T, подпро-
странства Ht(X) совпадают с интересующими нас
подпространствами
Ht(X) = Ht®, t0<t<T. (1.2)
Действительно, тогда можно было бы представить
подпространства Ht(g) в' виде ортогональной суммы
где каждое из ортогональных друг другу семейств
Ht(Xj), tQ<t< Т, описывается, как указывалось выше,
с помощью соответствующей структурной функции
Fj(t) = E\Xi(t)\\ tQ<t<T- (1.3)
Случайный процесс X (/) = {Xj (t)}™ с ортогональ-
ными компонентами Xj (/), tQ < t < Т (Xf (s) ± Xk (/)
для всех 5, t при / У= /?), каждая из которых пред-
ставляет собой процесс с некоррелированными прира-
щениями, удовлетворяющий условию (1.2), будем назы-
вать обновляющим процессом для случайного процесса
U0 = •&(<> k<t<T.
Отметим, что при нашем предположении (1.1),
согласно которому
lim Ht_h{X) — Ht(X),
h~>+0
обновляющий процесс является непрерывным слева.
Очевидно, обновляющий процесс всегда существует.
Его компоненты Xj (f), tQ< t < Т, можно построить,
например, следующим образом. Выберем какой-либо
элемент х1 е Н (I) и определим процесс с некоррели-
рованными приращениями ХД/) = PtxXi t0<t<T (на-
§ И
ВВЕДЕНИЕ
9
помним, что Pt обозначает оператор проектирования
на подпространство Если Н (Л\) У= Н (g), то
выберем затем какой-либо элемент х2еЯ(£), орто-
гональный подпространству Я(^), и определим про-
цесс с некоррелированными приращениями Х2(7) =
= Ptx2, t0<t<T. Поскольку подпространство Н (Л\)
и его ортогональное дополнение инвариантны относи-
тельно проекционного семейства Pt, tQ<t < Т, про-
цесс Х2(7), tQ<t<Tf будет некоррелирован с про-
цессом Jfi(7), t0<t<T. Если Я(Л'1)фЯ (Х2) У= #(g),
то выберем следующий элемент х3, ортогональный
подпространствам Н (Xj) и Я(Х2), и определим про-
цесс с некоррелированными приращениями Jf3(7) =
= PfXa, t0<t < Т, который будет некоррелирован с (7)
и Х2 (7), t0 < t < Т. Аналогичным образом выбирая по-
следующие элементы х^\
k-\
xkA. ф H(Xk),
мы в конце концов исчерпаем все пространство 77(g),
получив некоррелированные между собой процессы
= 70<7<Г; / = (1.4)
с некоррелированными приращениями (их число М
может быть бесконечным). Очевидно, X (7)=={Xy (7)}f\
t0 < t < Т, является обновляющим процессом для g (7),
70 < t < Т, поскольку Pt+hXj — PtXj _L Ht (g) при h > 0
и система величин Psx^ tQ< Pt+hXj—Ptx]f h> 0;
/= 1, ..., Л7, полна в H (g), так что система величин
PsXj, to< j = 1, ..., M, полна в Ht(g), и, таким
образом,
= Ш t.<t<T,
Если ввести гильбертово пространство С всех век-
торных функций с (7) = {с}- (7)}^, 7о<7<Г, с число-
выми компонентами ^(7), 70^7<Т, удовлетворяю-
щих условию
т м
j (!-5)
/=1
10
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
со скалярным произведением
т м
(С1/(С, {С2/(СеС.
t, 1=1
где
F/(0 = E|Xz(0P, t0<t<T,
— соответствующие структурные функции обновляю-
щего процесса X(t)== {Xz(/))^( t0<t<T (см. (1.3)),
то можно описать пространство Н (g) как совокуп-
ность всех величин вида
г м
(1.6)
*. i=i
где (cjf1 е С, и, очевидно,
t м
Pt^ ^ct(s)dX1(s), (1.7)
to /=1
поскольку Н (X) = //(£) и Ht (АЭ = Ht (£). Видно, что
подпространство Н,(£) = PtH (X) состоит из всех вели-
чин вида
t м
П = J ^icj(s)dXl(S), {ct(S)}“<=C,
t0 /=1
и, имея в виду унитарный изоморфизм ц <-> \Cj (s)}f*,
можно сказать, что семейство Н( (g), tQ<t < Т, того же
типа, что и семейство подпространств Cti t$<t <Т, ка-
ждое из которых образовано функциями {Cj(s)]M <= С,
обращающимися в 0 при s > t (для соответствую-
щего /).
Формула (1.6) в применении к значениям рассмат-
риваемого процесса £(/) = {g/ (/)}{”, tQ < t < Т9 дает так
называемое каноническое представление
t м
^(/)== J s)dX/(s)’ (1.8)
§ 1]
ВВЕДЕНИЕ
11
которое позволяет описать проекционные опера-
торы Pt, tQ<t < Т, в пространстве Н (£):
и М
/ ^c4(t, s)dXf(s), i=l, tn, (1.9)
to /=1
при tQ^u < t < T.
С точки зрения приложений, по-видимому, наи-
более важным следует считать вопрос об эффектив-
ных методах построения линейных преобразований,
связывающих процессы £(/) и X(f),—в частности,
эффективных методов построения канонического пред-
ставления (1.8). В теоретическом плане большой
интерес представляют и вопросы о том, как опреде-
лить тип обновляющего процесса X(f) = {Xf (/)}^,
tQ < t < Т, или как сравнить его с некоторым задан-
ным типом, зная корреляционную функцию исходного
процесса £(/) == [Ь (/))™, tQ<t<T.
Особый интерес к этим вопросам возник после
одного примера Г. Крамера*), показавшего, что обно-
вляющий процесс для одномерного непрерыв-
ного процесса g(/), /0 < t < Г, может быть, вообще
говоря, произвольным. Отправляясь от заданного
обновляющего процесса X (t) — {Xj (/)}f, tQ<t < Т,
произвольной структуры (с произвольными структур-
ными функциями Fj (t) = Е| Xj (/) |2, М)
соответствующий процесс £ (f), /0 < t < Г, может быть
построен следующим образом.
Пусть Ду, / = 1, ..., М, — непересекающиеся между
собой измеримые множества на интервале (/0, Г)»
м
(J Ду = (/0, Г) (mod 0), обладающие тем свойством, что
/«1
пересечение каждого из них с любым интервалом (s, f),
tQ< s < t < T, имеет положительную меру. Определим
интегрируемый процесс £(/), t0 < t < Т, над-
*) Н. Cramer, Stochastic processes as curves in a Hilbert
space, Теория вероят. и ее прямей. IX, 2 (1964), 193—204.
12
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. 1
лежаще выбрав функцию а (0 > 0 и положив | (0 =
= при /еДр М. Рассмотрим
t
g (/) = J i ($) ds. Очевидно, Ht (g) s Ht (X) при всех t,
^0
причем A’y(s) совпадает с производной g'(s) для почти
всех seAy, и поскольку при
а множества Ду всюду плотны в интервале (tQ, Т) и
процессы Xj(s) непрерывны слева, то JQ (s) s Ht(g)
при всех так что Ht (X) s Ht (£) и, следова-
тельно, X (/) является обновляющим процессом для
непрерывного процесса g (/), /0 < t < Г.
Приведем один из примеров такого типа множеств
Al, Д2, ...
Пусть Д’— канторово множество интервала / =
= (/0, Г), получающееся последовательным выбрасы-
ванием «смежных» интервалов /1/2; /1/4, /3/4’, /1/8, /з/8,
/5/8, Л/8 и т. д., которые выбираются так, чтобы мера
Н (Ю = Н U) — 5 Н (4) = т,
а
где т>0 — заданное число (m<Jp(/)).
На каждом смежном интервале/О1 «первого класса»
возьмем свое канторово множество Ка1 меры mai,
которое получается выбрасыванием из Iai смежных
интервалов /ai, a2 (a2 пробегает двоично-рациональную
последовательность 1/2, 1/4, 3/4,...). Положим
/С(1) = (J /Ох,- Далее, на каждом смежном интервале ZaiCt2
at
«второго класса» возьмем снова канторово множе-
ство меры maia?, получающееся выбрасыванием
из /aja2 смежных интервалов /aia2a3. Положим 7<(2) =
= U ZCata2. Продолжая эту процедуру последова-
(*1, а2
тельно для всех смежных интервалов 1а{ ... и пола-
гая (J ... ал, мы получим последова-
ar ...» а„
тельность непересекающихся множеств Д(°)=Д,
Д(1), /С(2), • • •> /С(п), • • •• Предположим, кроме того, что
§ I] ВВЕДЕНИЕ 13
00 oo
объединение (J KWl имеет полную меру: 2 И (К(п>) “
п—О п=0
— у. (Г) (этого можно добиться, выбрав надлежащим
образом меры та1 ... а„ > 0 отдельных канторовых
компонент 7Ц ... ап)-
Напомним, что канторово множество К.а{ ... ап
нигде не плотно в 1а{ ... ап и всякий интервал
slai...an имеет непустое пересечение с некоторым
смежным интервалом следующего класса 7а1 ... апа„+1.
Возьмем любой интервал (a, b)sl и рассмотрим
вложенный в него интервал (а', Ь'), а' = а^-г;
b' = b — е. Как было отмечено, (а', Ь') имеет непустое
пересечение с некоторым смежным интервалом 1а1‘,
положим б1 = (а', Z/)f|Zai. Интервал 6^/а, имеет не-
пустое пересечение с некоторым смежным интерва-
лом /а,/а2; положим 62 = 61 fl /а1а2. Уже ясно, что интер-
вал (а', Ь') при любом п имеет непустое пересечение
с некоторым смежным интервалом ...«„• Но
ц(At] ... ага)->0, и при достаточно больших п (для
которых ... а„) < е) соответствующий интервал
Iat ... ап, пересекающийся с интервалом (а4~е, b — е),
целиком входит в первоначальный интервал (а, Ь).
Итак, всякий интервал (a, b)^.I для всех достаточно
больших п содержит некоторый смежный интервал
«n-го класса» /Я1 ... ап, а следовательно, содержит и
множество К(п) с:/а1 ... ап положительной меры
ц (K(n)) — tnw- Разобьем теперь последовательность
целых положительных чисел п — О, 1, ... на счетное
число непересекающихся подпоследовательностей пп,
п12, ...; /г21, /г22, и положим
Р=1>2, ...
Ф=1
Очевидно, всякий интервал (a, b) s I содержит при
достаточно больших q множества К\пря> положитель-
14
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
ной меры так что
ц{(а, Ь) Л Др} >0
для всех р = 1, 2, ...
§ 2. Структурные типы и подчиненные процессы
1. Некоторые вспомогательные предложения *).
Пусть Ht (g), tQ < t < Г, — семейство монотонно воз-
растающих непрерывных слева подпространств
в сепарабельном гильбертовом пространстве Н (|),
и Pt, tQ<t < Т, — соответствующее семейство про-
екционных операторов (Pt — операторы ортогональ-
ного проектирования на Ht (g), tQ < t < T).
Возьмем произвольный элемент х^Н(£) и рас-
смотрим замкнутую линейную оболочку всех величин
tQ<t<T, обозначив ее Н(х). Подпространство
Н (х) состоит из всех величин вида
т
т)= J ф(0<*Ф(0, (2.1)
t0
где = t0 < t < Т, а функция <р(£) удовлетво-
ряет условию
т
J | ф(0 fdF(t) < ОО
t0
(F(t) = Е| Ф(/) I2 — структурная функция процесса с не-
коррелированными приращениями Ф (/), tQ < t < Т).
Пусть у^.Н(х). Согласно общей формуле (2.1)
величина у представима в виде стохастического
интеграла
т
У = / <fyx(t)d<I>(t),
to
*) По существу, здесь в основном излагаются хорошо
известные факты теории «спектральных типов» (см., например,
обзорную статью А. И. Плеснера и В. А. Рохлина,
Спектральная теория линейных операторов, Успехи матем.
наук I, 1 (11), (1946)).
§ 2] СТРУКТУРНЫЕ ТИПЫ И ПОДЧИНЕННЫЕ ПРОЦЕССЫ 15
так что
t
Р{У— J <Pyx(s)dO(s), t0<t<T,
и структурной функцией процесса Т (/) = Pty, t0 <
<t <Т, будет
t
G(O = E|^(OP=f I Фг/Х (s) Ws), t0<t<T.
Как мы знаем, пространство Н (у) состоит из всех
величин вида
п = / ф (/) dw (/) = / ф (0 <рух (О ЙФ (О,
fo ^0
где
т т
J | ф (0 I2 dG (/) = J | ф (/) <р^ (0 I2 dF (/) < оо.
t0
Обозначим А носитель меры dG(t). Очевидно,
т
J Ф (/) dT (t) = J ф (0 ач (/) = J ф (/) <?ух (0 йф (0.
f« Л Л
Поскольку dG(t)=] cpyX(t) \2dF(f), то почти
всюду на множестве Д, и, положив <р (/) = ф (/) (рух (0
при /еД, получим, что пространство Н(у) состоит
из всех величин вида
Т) = / <р (0 d(b (0, (2.2)
А
где 11 <р(0 fdF(t) < оо. Полезно отметить, что в слу-
Л
чае, когда меры dG(/) и dF(t) эквивалентны
dG ~ dF, Д есть одновременно носитель меры dF(t) и
имеет место равенство
Н(у) = Н(х).
16
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
Обратимся теперь к некоторому обновляющему
процессу X (/) = {Xj tQ<t < Т, определенному
соотношениями типа (1.4): Xj (/) = PtXj, М,
где элементы х; выбраны таким образом, что
Xy(s) ± Xk(t) для всех s, t при / =/= k, так что
м
Н®=® Я(х7).
/=i
Поскольку каждое из ортогональных между собой
подпространств Нинвариантно относительно
проекционных операторов Pt, t0<t<T, для всякого
элемента у^Н(£) имеем
м
Pty = Ф Р{Уи, t0<t<T,
a=i
где yk — проекция у на соответствующее подпростран-
ство Н (xk), k = 1, ..., М. Мы знаем, что структурные
функции Gk (0 = Е | Ptyk I2 абсолютно непрерывны от-
носительно соответствующих функций Fk — Е | Ptxk |2,
/0<^<Г, и, следовательно, структурная функция
G(t) = Е| Pty I2,
м
G{f)=^Gk(f}, t0<t<T,
k=l
очевидно, F* (/) = Е| Ptx* |2 при
. Мы видим, что все возможные «струк-
абсолютно непрерывна относительно функции
м
О) =2^(0, t0<t<T, (2.3)
k=l
где постоянные с^, сг2> ••• выбраны так, чтобы схо-
м
дился ряд 2 оЦхкII2
м ~~
* V
X = 2j GkXk
k=\
турные типы» dG(t) абсолютно непрерывны относи-
тельно некоторого максимального структурного типа
dF (/).
Пусть Н — некоторое подпространство, инва-
риантное по отношению к семейству проекционных
операторов Pt, t0<t<T, Назовем элемент х^Н
максимальным в Н, если для любого у^Н струк-
§ 2] СТРУКТУРНЫЕ ТИПЫ И ПОДЧИНЕННЫЕ ПРОЦЕССЫ 17
турная функция G (/) = Е | Pty (2, tQ<t < Т, абсолютно
непрерывна относительно структурной функции F(t) =
= Е | Ptx |2, /0 < t < Т (G (t) F (/)). Как было факти-
чески показано выше, максимальный элемент х^Н
с максимальным структурным типом dF (/) всегда
существует (см. (2.3)).
Пусть хх^Н(£) — максимальный элемент в про-
странстве и F{ (/) = Е| PtXi |2, tQ<t<T. Пусть х2—
максимальный элемент в подпространстве Н (%) QH (xj
и F2(t) — Е) Ptx212, t0 < t < Т. Вообще, пусть xk —
максимальный элемент в подпространстве
W)e[®W/)
и
/М0 = Е| Ptxk |2, to<t<T\
Очевидно, обновляющий процесс X (/) = {X; (/)}^,
tQ < t < Т, с компонентами Х^ — Ррсф t0<t<T-,
j = 1, ..., М, будет иметь структурные функции Fj(t),
tQ<t<M, упорядоченные в том смысле, что
dFi(t)>dF2(f)> ... >dFM(tY (2.4)
(т. е. каждая из мер dFj(f) абсолютно непрерывна
относительно предшествующих dFi(t), ..., dFy-JO).
Ниже мы покажем, что упорядоченные структурные
типы (2.4) определяются однозначно по семейству
tQ<t<T\ точнее, для любого обновляющего
процесса У (/) = {/*(< (Ht(Y) = t0<t<T)
с упорядоченными структурными функциями Gk (t) =
= Е| Yk(t)\2, iQ<t<T (т. е. такими, что dGJ/)^
>йЮ2(/)> ... >^(/)),
N = M, (2.5)
dGj(t)~ / = 1, M
(соответствующее число М обычно называют крат-
ностью).
Будем для краткости говорить, что мера dG(t)
подчинена dF(t), если она абсолютно непрерывна
относительно меры dF (/): dG (t) dF (/), и dG (t) орто-
18
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
гональна dF(t), если эти меры имеют непересекаю-
щиеся носители: dG(t) ± dF(t).
Пусть dF (/)— произвольная положительная огра-
ниченная мера на отрезке [/0, Г), подчиненная макси-
мальному структурному типу dF* (t) ~ dF{ (I). Суще-
ствует элемент x^H(g) со структурным типом,
в точности равным dF(t), т. е. такой, что
t
Е|Лх|2= j dF (s'), t0<t<T,
to
например, можно взять
*0
где, напомним, (/) = Е| Х1 (/) |2 — максимальная
структурная функция (см. 2.4)). Очевидно, все эле-
менты у е //(g) со структурными типами dG(f), под-
чиненными dF (t), образуют линейное подпрост-
ранство, которое инвариантно относительно проекто-
ров Pt, tQ<t <Т.
Рассмотрим систему элементов у{, у2, .77(g)
с одним и тем же структурным типом dG(t), такую,
что подпространства H(yk), порождаемые величи-
нами Ptyk, tQ<t < Т, ортогональны при различных
k=l, ..., п. Систему z/b ..., уп назовем макси-
мальной, если ее нельзя расширить; точнее, если не
существует элемента z/^7/(g), ортогонального под-
пространствам H(yk), k= 1, ..., п, и имеющего
своим структурным типом dG(t\ Очевидно, макси-
мальная система существует.
Лемма 1. Если у19 ..., ym — максимальная си-
стема, а у'{, ..., у'п — некоторая система с тем же
структурным типом dG(t), то
п<т.
(2.6)
Доказательство. Предположим сначала, что
m
подпространство L— ф Н (у^ совпадает со всем про-
§ 2] СТРУКТУРНЫЕ ТИПЫ И ПОДЧИНЕННЫЕ ПРОЦЕССЫ
19
странством Н (£). Если ввести У, (/) = PtyfJ tQ < t < Т
/==1,..., m, то элементы y'v ..y'n^L можно
представить в виде
т щ ,
и 1=1
где функции ckj(t) удовлетворяют условию
t, f=i
Поскольку все элементы y'v .... у'п имеют одну и
ту же структурную функцию
G'(O==J ^\ckj(s)?dG[s), k = \,...,n,
*0 /=1
а при условии ортогональности подпространств
Н (у'), Н(у'^ должно быть
t m
(РМ) = J S Cii (s> dG (S) = °’ еСЛИ Z’ k’
Г. /=1
TO
tn tn _______
Si ck{(t) |2= 1, 2сц(0М) = ° при i=£k
для почти всех t относительно меры dG(t). Видно,
что векторы ck (/) == {cki (t)}™, k=l, n, с компо-
нентами /=1, ..., tn, в m-мерном векторном
пространстве (с обычным скалярным произведением)
образуют ортонормированную систему. Как известно,
число элементов такой системы не превосходит раз-
мерности пространства, т. е.
Общий случай сводится к только что рассмотрен-
ному с помощью следующего приема. Введем про-
странство Н всех элементов, структурные типы ко-
торых подчинены dG (/), а в нем подпространство
20
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
т
Е—@Н(у]) и его ортогональное дополнение LQ =
— HQL. По определению максимальной системы
yl9 ут в подпространстве Lo нет ни одного эле-
мента со структурным типом dG(f). Поэтому, если
dG0(t)—максимальный структурный тип в инвариант-
ном подпространстве £0, то dG (/) = dG$ (/) ф dG (f),
где dG(t) — некоторая ненулевая мера, ортого-
нальная dGQ(t). Рассмотрим новое пространство
Н — подпространство в Н из всех элементов со струк-
турными типами, подчиненными dG(t). Легко видеть,
что поскольку меры rfGo(O и dG(f) ортогональны,
подпространства £0 и Н с максимальными структур-
ными типами dG0(/) и dG(t) будут также ортого-
~ т
нальны, и потому H^L, где, напомним, L= ф Н(у^.
Если До и Д— непересекающиеся носители орто-
гональных Mep.dGo(Z) и dG(Z), то всякая величина
x\^H{yk) может быть представлена в виде
т
п = j (t) dW, (0 = J (0 dv} (О Ф / ф (0 dwl (/)
^9 До Д
(Т/(0 = р^, t0<t<T),
где
По = / Ф (0 I (0 <= Ль Й = j Ф (0 dWf <= Н.
До А
Ясно, что элементы 2/ — | j(f), j—l, .... т,
д
с одним и тем структурным типом dG(t) образуют
в Н максимальную систему, причем подпространство
т ~
£=ф//(3у) совпадает с пространством Н. Таким
образом, вместо исходных систем ylt ..., ут и
у\, ..., у'п со структурным типом dG(t) мы можем
рассматривать системы элементов zjt j=l, ..., mt
§ 2] СТРУКТУРНЫЕ ТИПЫ И ПОДЧИНЕННЫЕ ПРОЦЕССЫ
21
И z' = p4"fc(/), k = \..n (Tft(0 = W> t,<t<T)
д
в пространстве И со структурным типом dG(t). Как
было установлено выше (в случае L = Н), имеет
место неравенство п^т.
Вернемся к обновляющему процессу X(t)—{Xj (/)}“,
t0 < t < T, с упорядоченными структурными функ-
циями F!(t) = E\X!(,t)?, t0<i<T (см. (2.4)).
Пусть Y (/) = {Уй tQ<t < Т,— некоторый про-
цесс с некоррелированными компонентами вида
Yk (t) = Ptt/k, t0<t <T, структурные функции кото-
рого Gk(t) — Е| Yb(t) [2, t0 < t < Т, также упорядочены:
dG\(t)>dG2> ... >dGN(t).
Лемма 2. Имеют место следующие соотношения:
N^M,
dGk(t)<dFk(f), k —
(2.7)
Доказательство. Предположим, что неко-
торая мера dGn (/) не является подчиненной соответ-
ствующей мере dFn(t~). Тогда существует нену-
левая мера do (/), подчиненная dGn (7) и ортого-
нальная dFn(t). Поскольку меры dP}(t) при j^n
подчинены dFn(f), все они будут ортогональны da(t).
С другой стороны, daft) подчинена максимальному
структурному типу dF\(f) и, вообще, каким-то dFj(f),
1 т, где число т всех таких структурных типов
из последовательности dF{ (t) )> dF2 (t) )> ... строго
меньше п. Если xh х2, ..., хм — элементы, поро-
ждающие обновляющий процесс X (/) = {Xt (Z)}^
(Xj(t) = Ptxlt то, учитывая структуру
ортогональных подпространств И (xft) (см. (2.1), (2.2)),
легко понять, что элементы 2/~ J ^/(0> j— 1>..-> tn,
д
где А — носитель меры do(t), образуют максимальную
систему в Н (g) со спектральным типом da{t). Соот-
22
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
ветствующие элементы j dYk(t), k = l, .п,
д
также образуют систему со спектральным типом do(/),
причем п > т, что противоречит неравенству (2.6)
леммы 1. Следовательно, на самом деле мера dGn(t)
подчинена dFn(t) при любом п. Фактически доказано
также, что N М.
Добавим здесь, что если рассматриваемый про-
цесс Y (/) = {Yk (/)}/, t0<t<T, сам является обно-
вляющим: Ht(Y) = tQ<t<T, то процессы
Y (/) и X (/), фигурирующие в лемме 2, можно поме-
нять местами, и в этом случае вместо «неравенств»
^2.7) будем иметь N = M, dGj(t)~dF/==1,..., N
(см. указанные ранее соотношения (2.5)).
Пример. Пусть g (/) = (Z)}7, /0 < ^ < Л — про-
извольный процесс с некоррелированными прираще-
ниями. Предположим, что заданы
^(0 = EU0W •••>'"•
Найдем соответствующие структурные типы dFj (/)>*...
...^dFM(t). Очевидно, максимальным типом будет
т
dF4t)~2fidFkk(t),
где строго "положительные .... <т^ выбраны так,
чтобы функция 2i&kFkk(t) была ограниченной. Вместо
самих структурных типов dFj (t) можно указать плот-
ности fi(t) — dFj^dF* (/) или множества
А/={Л 7/(0 > 0}
— носители мер dFj(t), j = l, ..., М. В рассматри-
ваемом случае носители \} можно описать, обра-
тившись к матричной функции f (/)}, t0 < t < T,
с компонентами
dFij (t)
fij (t) = о,. dF,-^ <Jlt i, / = 1, ..., tn.
§ 21 СТРУКТУРНЫЕ ТИПЫ И ПОДЧИНЕННЫЕ ПРОЦЕССЫ 23
Именно, если ... —упорядо-
ченная система всех ненулевых собственных зна-
чений матричной функции f(t), то
Д, = {t: Kj (0 > 0}, /=1.............М.
Соответствующий обновляющий процесс X (/) =
= {X,(7)}|М, t0 <t<T, с упорядоченными структурными
типами может быть получен как
t tn
t0 Z = 1
где и (/) = {щ/ (/)) — унитарное преобразование, при-
водящее положительную матрицу к диа-
гональному виду (с элементами
kk(f)==O при М < k^n, на главной диагонали).
2. Подчиненные процессы. Рассмотрим пару про-
цессов: £(0 —UiWir И п(0= Говорят, что
процесс т](/) может быть получен из процесса £(/)
неупреждающим (линейным) преобразованием, если
(2.8)
Обозначим Ht (г])1 ортогональное дополнение в п р о-
странстве Н(ц) к подпространству Н((т\У, ЯДт))1 —
= Н (n) © Ht (-и). Аналогично, ЯД^)1 = Я(£)@ЯД£).
Будем говорить, что процесс ц(/) подчинен процессу
£(/), если
ЯДт))еЯД£) и ЯДт])х S ЯД|)х, t0<t<T. (2.9)
Второе из соотношений (2.9) означает, что всякая
величина уеЯ(т|), ортогональная подпространству
ЯДт]), должна быть также ортогональна и более
широкому подпространству ЯД£).
Рассмотрим два важных примера подчиненных
процессов.
24
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
Пример (каноническое представление). Пусть
I(0= {Ь (0)Г» to<t < — процесс с некоррелиро-
ванными приращениями и
t m
n* = J S Cki V’ s) d^l (s)’ to < Z < T’ (2 10)
k = 1, ..., n.
Обозначим Qt операторы ортогонального проектиро-
вания на подпространства Ht (г|), /0 < t < Т. Выраже-
ние (2.10) называется каноническим представлением,
если проекции значений на подпространства
#иСп)> u^t, получаются по формулам
Q«nfe(0= f У ck}{t, s)dZj(s), t0<t <Т,
/tf (2.Н)
k = 1, . .., n.
(Например, так будет в случае, когда £(/), t0 < t < Т,
является обновляющим процессом для ц(/), /0 < / > Т.)
Очевидно, представление процесса '*1(0= {лИО}”
в виде стохастического интеграла (2.10) возможно
тогда и только тогда, когда выполнено соотношение
(2.8). При условии же (2.11) имеем
nfe (0 — (0 = / У Ckj (t, s) d% j(s);
U /—1
видно, что величины (t) —Q^i (/), и; k = 1,n,
порождающие подпространство //„(ц)1, ортогональны
Яи(£), и, таким образом, имеет место соотношение
(2.9), т. е. процесс г)(/) = {тр. (/))", t0 < t < Т, подчи-
нен процессу с некоррелированными приращениями
Пример. Пусть g (/), t0 < t < Т, — случайный про-
цесс со стохастическим дифференциалом
§ 2] СТРУКТУРНЫЕ ТИПЫ И ПОДЧИНЕННЫЕ ПРОЦЕССЫ 25
точнее,
t
Ц0=/ a(s)ds + W(t), t0<t<T, (2.12)
/о
где а (/), /0 < t < Л — некоторая случайная функция
со средним Еа (/) = 0 и конечным вторым моментом
Е| a(t) |2,
t
j Е | a (s) |2 ds < oo, tQ < t < T,
A)
и W (/), t0 < t < T, — процесс с некоррелированными
приращениями, такими, что при h > 0
+ tQ<t<T
(величины W (t + h) — W (/), h > 0, некоррелированы
co значениями g(s), s^Z).
Рассмотрим процесс ц(/), t0 < t < T, получающийся
из процесса g(/), /0 < t < Г, линейным неупреж-
дающим преобразованием вида
t
п (0 = ЦО — / «(«)ds =
t,
t
= J[a(s) —4(s)]Js + IF(0, tQ<t<T,
где d(s) — Psa(s) обозначает проекцию величины a (s)
на подпространство Hs (g). Поскольку a (s) — d(s)LHt (g)
при приращения
t+h
П(/ + ^)-П(0= J {a(s)-a(s)]4/s + [F(/ + A)-r(/)]
t
обладают тем же свойством, что и приращения
IF(/ +A) —а именно, при А>0
+ — т](0 ± t0<t<T.
При этом
26
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
так что ц (0, to<t < Г, — процесс с некоррели-
рованными приращениями. Для такого типа
процесса подпространства (т))1 порождаются при-
ращениями т](/ + А)—т](/), h > 0, а мы видели, что
они ортогональны соответствующим пространствам
/7Д£), и, таким образом,
t0<t<T,
т. е. процесс ц (/) подчинен процессуg(/),t0 < t < Т.
Обратимся к произвольной паре процессов |(/) =
= & « и Т) (0 = {Т]А (<, /о < t < т. Пусть dFx (/) >
>...> dF м (t) и dG{ (/)>...> dGN (/) — упорядо-
ченные структурные типы их обновляющих процессов,
Теорема 1. Если процесс т](0, t0 < t < Т, под-
чинен процессу g (t), t0<t < Т, то
dGk(t)<dFk(t), k = l, .... N. (2,13)
Доказательство. Обозначим Pt и Qt опера-
торы ортогонального проектирования на подпростран-
ства Ht (£) и Ht (ц), t0 < t < Т. В силу соотношений
(2.9) имеем
ЗД(ц) = ДДт1), ад(п)х = 0
И
PtH (Т)) = Pt [Ht (п) е Ht (п)Ч = PtHt (n) =
= = (2.14)
Для проекторов Pt и Qt отсюда вытекает, что опера-
тор Pt в инвариантном подпространстве Н (т]) s Н (£)
совпадает с оператором Qt) tQ<t < Т. Следовательно,
обновляющий процесс Y (/) = {}% для процесса
т](/), tQ<t<T, получается по формуле
Yk(t)=-=Ptyki k = \y ..., N.
где z/i, ..., yN — некоторая система величин в под-
пространстве Н (ц) Н (£), и теорема 1 является,
таким образом, простым следствием доказанной ранее
леммы 2.
§ 2] СТРУКТУРНЫЕ ТИПЫ И ПОДЧИНЕННЫЕ ПРОЦЕССЫ
27
Попутно отметим здесь, что обновляющий про-
цесс Х(/)== {%,(/)} f, tQ<t<T, для процесса £(/),
/0 < t < Т, обладает следующим свойством. Инте-
гральное представление
t м
П* (О = / S ck! (t, s) dX^s), t0<t<T,
t, /=i
k = 1, ..., n,
подчиненного процесса я (t) — {т)& (/)}", t0 < t < T, яв-
ляется каноническим:
и M
Qu^k (i) = Pa^]k w = f S Ckl S) dX№’ u^(-
При рассмотрении процесса t](0, tQ < t < T, полу-
чающегося из процесса g (t), t0<t <T, некоторым не-
упреждающим линейным преобразованием (ЯДт)) s
s Ht (g), t0 < t < T) может возникнуть вопрос о том,
является ли это преобразование «обратимым», точ-
нее, вопрос о том, выполняется ли условие
= tQ<t<T. (2.15)
Если рассматриваемое преобразование таково, что
процесс т](7) подчинен процессу g (/), /0 < t < Т, то
условие (2.15) равносильно (вообще говоря, более
слабому) условию
н®=щц), (2.16)
поскольку в этом случае
= PtH ® = =
при всех t, t0< t < Т (см. (2.14)).
Теорема 2. Если процесс т](/) подчинен про-
цессу Z (f), to<t<T, и они имеют обновляющие про-
цессы одного и того же типа кратности
М < <х>
(dGk(t)~ dFk(t), k = i.......М),
то имеет место соотношение (2.15).
28
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
Прежде чем доказывать эту теорему, отметим,
что она неверна в случае бесконечной кратности
М — оо. Например, если £(/) = &(/)}“» to<t <Т,—
процесс с ортогональными компонентами, каждая из
которых является процессом с ортогональными при-
ращениями, причем структурная функция F(t) =
= E|gz(/)|2, t0<t<T, одна и та же для всех
/=1,2,..., то процесс г] (/) = {т)й (/)}~, t0 < t < Т,
с компонентами (О — b+i (/), 6 = 1,2... будет
того же типа, что и процесс £ (0 = {£/(0)Г> будет
подчинен g(/), t0 < t < Т, но, очевидно, Н(г)) не со-
держит значений g( (t), t0<t < Т, так что /У (r))=#//(g).
Доказательство теоремы 2. Предположим,
что Н (т|) #= Н (g) и возьмем ненулевой элемент
Пусть F (/) = Е | Ptx |2, t0 < t < Т.
Среди упорядоченных структурных типов dFx (/) )> ...
.. .>dFM(t) процесса g(/), t0 < t < T, выделим те,
которым подчинен тип dF(fy, пусть это будут
dFt (/) >- ... >> dFp (/), где dFр (t) — самый «младший»
из них (р^М < оо). Если xlt ..., хм — система эле-
ментов в пространстве Н (g), порождающая обно-
вляющий процесс Х(/) = {X/ (/)}“, t0<t< Т; Xj(f) =
= PtXj, и F, (/) = Е | Ptx, |2, to < t < T:
то система элементов
Xj = j dXjit), j=l, p,
Д
где A — носитель меры dF(t)y очевидно, буде^ макси-
мальной ортогональной системой со структурным
типом dF (t) (см. определение на стр. 18).
Не ограничивая общности, можно считать dGj(t) =
= dF j(t), взяв соответствующие элементы уъ ..., ум^
порождающие обновляющий процесс Y (/) =
= {Y1 (С, to < t < Т (У, (/) = Ptyi, Gj (t) = Е | Ptyi |2,
tQ < t < Г). Ортогональная система элементов
/=!>•••> P>
Д
§ 3] НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ 29
с тем же структурным типом dF (t) и той же крат-
ности р, что и максимальная система xJf j = 1, ..., р,
может быть дополнена выбранным ранее элементом х,
и, таким образом, в Н (g) найдена система у[9 ..., ур, х
кратности р +1, большей, чем максимальная крат-
ность р, что невозможно (см. неравенство (2.6)).
Следовательно, предположение о том, что Н (ч\)^Н (g),
является неверным, и тем самым теорема доказана.
§ 3. Некоторые примеры обновляющих процессов
1. Марковские процессы. Мы рассмотрим ниже
два важных класса процессов — марковских и ста-
ционарных в широком смысле.
Многомерный случайный процесс g (/) = {gz (/)}™
называют марковским в широком смысле, если для
любого t проекции Pfa (t + h) величин gt- (t + h), 0,
на все «прошлое» //z(g) принадлежат замкнутой
линейной оболочке значений gz(Z), Z = 1, ..., tn.
Ограничимся расмотрением лишь конечномерных
процессов. В этом случае для любого
tn
Psb(t) = ^c4(t, (3.1)
где коэффициенты cZ/(Z, s) легко определяются по
корреляционным функциям
s) = Egz(OiT(s), Z, /=1, т.
Например, в случае линейно независимых
значений gy (s), / = 1, ..., tn, когда матрица В (s, s) =
= {Bij(s, s)} является невырожденной и существует
обратная матрица В (s, s)“! = {bki (s, $)}; коэффициенты
в выражении (3.1) можно найти по формуле
Q/(Z, «)= 2 Bih(t, s)bkl(s, s), i, j=\...m.
k
Согласно общему предположению (1.1)
lim =
s->t
30
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
и в случае линейно независимых значений
/==1, т, имеется некоторый интервал /j s t
такой, что величины
= /=1, т,
образуют базис в линейной оболочке значений Е/(s),
Z=l,..., tn. Поэтому ортогональное дополнение
Ht (g) @ Htx (g) к подпространству (£) в простран-
стве порождаемое всеми величинами
(s) — Pt& (s), /1 < s < /; i = 1, ..., m,
совпадает с замкнутой линейной оболочкой величин
Л/ (з)—Л/ (О = Л/(з)—^.Л«(Д z = l,...,m.
Случайный процесс г] (s) = {т)г ($)}{", tt s t, является
процессом с некоррелированными приращениями и
обладает тем свойством, что
^2(g)e^,(g)=^(n)e^.h) (з.2)
при любых t\ S] s2 t, ив этом смысле является
обновляющим процессом для g (s) = {gz (з)}™,
(Можно затем перейти к обновляющему
процессу с некоррелированными компонентами и
упорядоченными структурными типами, как это сде-
лано в примере на стр. 22.)
Простейшим примером марковского в широком
смысле процесса может служить многомерный про-
цесс g (/) = {£; (/)}”, t > t0, удовлетворяющий системе
стохастических дифференциальных уравнений вида
т
^z(0=2 + z = l, ..., tn, (3.3)
/=i
или, в матричной форме,
dl(t) = A(t)Z(t)dt + dx](t),
где A (t) = {Atj (i)} — непрерывная матричная функция
(неслучайная), а т](0 — (Л/ Ю}Г> —процесс с не-
коррелированными приращениями.
§ 3] НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ
31
Решение системы (3.3) при заданных начальных
значениях £ (/0) = (4)}”, точнее, решение интеграль-
ного уравнения
t
&(/) = £ (/о) + J A (s) I (s) ds + [n (0 - T) (Ш t > /о- (3.3')
t<s
может быть представлено в виде
t
= Ш Ш(/о)+/Я(Л (3.4)
tr\
где m X-матричная функция R(t, s) — {Rl{(t, s)},
t^s,—так называемая резольвента—есть ре-
шение обыкновенного дифференциального уравнения
d^ ^’-s—= A(t) R(t, s), t> s,
с начальным условием R(s,s) = / (I — единичная
матрица). Действительно, по самому определению
имеем
t
R(t, s) = 1 + j A (и) R (и, s) du
s
И
t
s
t
t
j R(t, s) dv\ (s) = j I + J A (w) 7? (u, s) du d^is) —
to to L
t
= [»] (/) — П (f0)] + J J A (u) R (u, s) du dx\(s) =
to
t Г и
= h (/) —11 (f0)] + J A («) J R(u, s) df](s) du.
Видно, что
to
0>
to
32
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
есть решение интегрального уравнения (3.3') с на-
чальным значением |(/о) = О; разность же х(/) =
— W)~W), как решение однородного диф-
ференциального уравнения dx (t)/dt = А (/) х (/) с на-
чальным условием x(tQ) = имеет вид
x(t) = R (t, to) g (t0).
Предположим, что приращения T)z (/)—т|£- (s),
/=1, т, некоррелирован ы с начальными
значениями (/0), j = 1, ..т. Тогда, если переопреде-
лить процесс т](/)={т]г- (/)}™, t < tQ, точнее, рассматривая
вместо первоначального процесса т)(/), t < tOf процесс
с некоррелированными приращениями W)— Wo)+
+ £(/0), можно (не меняя обозначений) считать
Пг (0 = ^(0, 1=1,..., т (3.5)
и переписать выражение (3.4) в виде
t
R(t, s) dt\(s). (3.6)
^0
Очевидно, подпространства ЯДт]) порождаются
значениями (/0) и приращениями тр ($) — T]z (/о),
/0< s^t. Используя тот факт, что резольвента R(t, s)
удовлетворяет условию
+ s) = R(t + h,t)- R(t, s)
при всех t^s и /г^О, из формулы (3.6) получаем,
что
t
R(t + h,t)Z(f)=A RU + h, t) R (t, s) (h}(s)
= j R(t + h, s) dx\(s)
h
и
t+b
l(t-^h)-R(t + h, J R(t + h, s)dr\(s).
t
Видно, что разности
litt+V-^Rtjit + h, t)li(t), i^l, ..., m,
§ 3] НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ 33
ортогональны подпространству т. е. величины
т t т
+ 1)^(1) = J h, s)dr\j(s)
/=i и. /=1
являются проекциями значений g; (t + Л), i = 1.т,
на подпространства s //z(r|).
Таким образом, процесс £(/) = {М0}Г, t^t0,
является марковским в широком смысле, а фор-
мула (3.6) дает нам его каноническое предста-
вление с помощью описанного выше (см. (3.5)) про-
цесса г] (/) = {т](- (Z)}j , t^tOt с некоррелированными
приращениями.
Очевидно, Ht (t]) S Ht (g) (см. (3.3')), т. е. этот про.
цесс является обновляющим для §(/) = {?;(/)}”,
а именно,
при всех t^t0.
2. Стационарные процессы. Рассмотрим стацио-
нарный в широком смысле процесс g (t) — {gj. (/)}{",
— oo < t < oo, т. e. такой, что корреляционные функ-
ции Bkl (t, s) = Eg* (/) gy (s) зависят лишь от разности
t — s:
M,s) = M-s)’ fe,/ = l,...,m.
Стационарный процесс g (/) = {gi (f)}", — oo < t < oo,
называют регулярным, если
f|^(g) = °. (3.7)
t
Как известно *), регулярный процесс имеет спек-
тральную плотность / (Л) = {^/(Л)}; это положительно
определенная матричная функция от X, — оо < А, < оо,
компоненты которой связаны с корреляцион-
ными функциями соотношением
оо
Bfe/(0= J e^fkji^dKj, m.
*) См., например, Ю. А. Розанов, Стационарные слу-
чайные процессы, М., Физматгиз, 1963,
2 Ю. А. Розанов
34
обновляющие процессы
[ГЛ. 1
Более того, условие (3.7) равносильно тому, что
спектральная плотность f (X) почти всюду имеет один
и тот же ранг и представима в виде
/(А) = ^ф(А).ф(А)* п.в., (3.8)
где tn X «-матричная функция <р (Z) = {фА/-(А.)} удовле-
творяет условиям
оо
/ || ф (А) ||2 dK < оо,
—-оо
(3.9)
J ешф (A) dA = 0 при t < О,
—-оо
а ф(А)* обозначает матрицу, сопряженную к ф(А).
Обозначим И2 класс всех матричных функций-
ф(А), удовлетворяющих условиям (3.9). Преобразо-
вание Фурье
оо
= i pw<p(W (3.10)
—оо
функций из Н2 обращается в 0 при t < 0 и
Ф (X) = [ е~шс (/) dt
о
аналитически продолжается в нижнюю полуплоскость
ImX<0 комплексного переменного X.
Отправляясь от спектрального представления
= feiMd<Pk(K), k = (3.11)
—оо
где Ф (X) = {Ф& (Х)}^, —оо < X < оо,—процесс с не-
коррелированными приращениями, такой, что
н
ЕФ*(н)Фу(н)= ^fkl(y)dk, k,/=1, ..., m,
§ 3] НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ 35
можно получить следующее выражение для рассма’
триваемого стационарного процесса:
t п
lk(t)= f (3.12)
— oo /=!
где m X «-матричная функция c (t) = {ckj (/)}, />0,
есть преобразование Фурье матричной функции <р(Л) =
= {фй;(А)} класса Я2, удовлетворяющей условию (3.8),
а сам стохастический интеграл в правой части вы-
ражения (3.12) отвечает некоррелированным
приращениям
оо
Г piKt — ,
n/W-n/(s)= J -e ./......ф(А)~‘ЛФ(Л), (3.13)
—сю
j = 1, ..., n.
Уточним здесь, что, во-первых, ср (А,)-”1 есть п X w-
матричная функция, обратная к (р(Л) (т. е. про-
изведение ср (Л)-”1 ср (Л) есть п X ^-единичная матрица),
и, во-вторых, для непересекающихся интервалов (sb
и ($2, Q величины nz (^1) — П/и П/ОУ — П/(^2) не-
коррелированы, причем
E|n/W-n/(s)|2 = /-s, /=1, п, (3.14)
а величины (У — i}k (У и П/(У~П/(52) некоррели-
рованы, для любых интегралов, если k =/= /. Можно
рассматривать выражение (3.12) как неупреждающее
линейное преобразование я-мерного обобщенного про-
цесса л (0 = {'П/(0}? типа «белого шума», переписав
(3.12) в виде
Ц0 = J c(f —s)r)(s)ds. (3.15)
— 00
Обозначим ЯДц) подпространство, порожденное
всеми величинами цД/)—t]/(s)> / = 1......
и Qt — оператор проектирования на Поскольку
Ht(i])=>Ht(l), —00 </<00,
2*
36
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
имеют место следующие неравенства: для любого
Л>0
h т п т
I ckj(s) |2ds= 2 El b{t + h) -Q£k+
0 fel >1
m
C 2E|^(/ + A)-P^(/ + A)f,
k—i
где Pt, как и раньше, обозначает оператор проекти-
рования на подпространство Если допустить,
что среди матричных функций ф(А) класса Я2, удо-
влетворяющих условию (3.8), найдется функция ф°(А),
для которой соответствующие приращения (/) — (s),
j = 1....п, дают равенство
= (3.16)
оо
то матричная функция cQ(f) = J ешф° (Л) tZZ по от-
ношению ко всем матричным функциям c(t) —
оо
— J еш<р (A.) dK будет максимальной в том
— оо
смысле, что
h т п т
о k=\ /=1 k=\
h т п
I S S' Ck>^'2rfs
о £=1 /=1
при всех Такая функция ф°(Л)^Я2 действи-
тельно существует и определяется тем условием, что
она максимальна в прямом смысле этого слова:
пт пт
II <Р° (А) II2 = 2 21 Ф°А/ (А) I2 > 2 2 I ф6/ (А) I2 = (| ф (А) II2
т k 1 II 1 1 (3.17)
при всех комплексных А, 1тА<0.
Фактически мы уже описали обновляющий про-
цесс, задав соответствующие приращения tjJ(/) —
§ 31 НЕКОТОРЫЕ ПРИМЕРЫ ОБНОВЛЯЮЩИХ ПРОЦЕССОВ 37
/=1, т. Как показывает общая формула (3.14),
кратность М обновляющего процесса равна п — рангу
спектральной плотности f(K):
М = п,
а в упорядоченной системе структурных типов
dF{ (/)>...>• dFM (t) (см. (2.4)) все dFj(t) эквива-
лентны лебеговской мере
dFj(t)^dt, /₽1.......М. (3.18)
Несколько иначе обстоит дело, если стационарный
процесс £(0 = {|i (t)}™, — оо < / < оо, рассматривается
лишь на конечном интервале t0 < t < Т (или при
Поясним это на примере одномерного стационар-
ного процесса £(/) со спектральной плотностью вида
п
vjifi P(z)=^i pkZk — полином степени п, все корни
fe=0
которого лежат в левой полуплоскости Re z < О
(в частности, р0 0). В этом случае максимальной
функцией, удовлетворяющей условию (3.8), будет
а формула (3.13) дает
f № _
П (0 - n(s) = ----.т— Р (Л) ЙФ (А) =
J Ifv
с piM_JKs
= а -Ро^Ф(М +
J L'v
— ОО
00 п-1
+ / - eiKs) 2 Pfe+i W =
— оо k—Q
* П — 1
=ро J । («)<*«+S pk+i (*)]• (3-20>
38
ОБНОВЛЯЮЩИЕ ПРОЦЕССЫ
[ГЛ. I
Если это соотношение переписать в дифференциаль-
ной форме
dr\ (t) = 2 pkl^ (t) dt + pn (t),
k*=0
то легко заметить, что многомерный процесс (^ (/)}”,
с компонентами
k = \,...,n,
удовлетворяет системе стохастических дифференциаль-
ных уравнений типа (3.3), а именно,
= k = \, ..., п-1,
= + (3’21)
\ Рп / Рп
Кроме того, как и в случае общего стационарного
процесса, определяемые выражением (3.20) величины
т](/) — л^о), ортогональны к замкнутой линей-
ной оболочке всех значений g(s), и> в част-
ности, к величинам
= k = l,..., п. (3.22)
Таким образом, многомерный процесс (/)}”,
принадлежит к тому типу, что был уже рас-
смотрен нами, и из полученных выше результатов
(см. стр. 33) следует, что в рассматриваемом случае
обновляющим для {^(/)}р будет процесс
{лИО}?» t^tQ, с компонентами
= k = 1, ..., п — 1,
nn(O = n(O-n(^o) + ^(^ ( }
Очевидно, этот же процесс будет обновляющим
и для исходного стационарного процесса £(/), t^tQ.
ГЛАВА II
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
§ 1. Изометричные семейства и некоторые примеры
Пусть Ht, t$<t < Т, —семейство подпространств
в гильбертовом пространстве Я, удовлетворяющее
условиям (1.1) гл. I. Пусть Ut, tQ<t<T,—другое
семейство (в гильбертовом пространстве Я). Будем
называть эти семейства изо мет ранными, если суще-
ствует изометрический оператор X из U в Н такой,
что
Ht = XUt, tQ<t<T. (1.1)
Чтобы не вводить новых обозначений, будем счи-
тать, что Н и U являются замыканием объединения
всех подпространств Ht, tQ<t<T, и Ut, tQ<t <Т
соответственно; тогда X — унитарный оператор,
отображающий U на Н.
Обозначим Pt оператор проектирования на Ht,
Qt — оператор проектирования на Ut, Оче-
видно, при условии (1.1)
Pt = XQtX~\ t«<t<T, (1.2)
поскольку для унитарного оператора X вместе
с условием (1.1) выполняется также условие Ht~ =
= xuj- (где = U^UQUt), и потому
при h^Ht
X'lh<=Ut, QtX~'h = X~xh, XQtX~'h = h,
а при h ± Ht
X'‘hLUlt QiX~ih = G, XQtX~lh = Q.
40
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. И
Возьмем произвольные элементы и[9 um^U,
порождающие «обновляющий» процесс /0 <
< t < Т, для семейства Ut, tQ<t<T (т. е. подпро-
странства Uky порождаемые всеми величинами Qtuki
tQ<t < Т, ортогональны между собой при различных
k = 1, ..., N, и U = ф Uk). Положим Fk (/) = || Qtuk ||2,
k=\
tQ<t<T; 6 = 1, N. Очевидно, соответствующие
элементы xk = Xuk, k = \, ..., N, порождают в про-
странстве Н обновляющий процесс [Ptxk]^9 tQ <t < Т,
точно такого же типа, поскольку
Ptxk = XQtuk, t0<t<T,
откуда видно, что подпространства Hk = XUk (поро-
ждаемые величинами Ptxk, tQ<t<T) ортогональны
между собой и
II Ррсь II2 = II XQtuk II2 = II Qtuk II2 = Fk (/),
t0<t<T-, k = \, ..., N.
Таким образом, для изометрических семейств Ufi
tQ<t < Т, и Ht, tQ<t < Т, обновляющие процессы
имеют один и тот же тип.
Верно и обратное утверждение. Действительно,
рассматривая «циклические» подпространства Uk и Hk
с эквивалентными структурными типами, можно вы-
брать порождающие элементы uk е U k и xk е Hk так,
чтобы их структурные типы в точности совпадали:
IIQ/«dl2 = ll^ll = ^W. * = 1, .... N.
N X
Пространства U — ф Uk и Н=® Hk унитарно изо-
k=\
морфны пространству С всех векторных функций
c(t)=^ [ck (Z)}^, tQ^t< T, с компонентами, удовлетво-
т N
ряющими условию J ck(t) ^dFk(t) < оо, в котором
/0
§ И
ИЗОМЕТРИЧНЫЕ СЕМЕЙСТВА
41
скалярное произведение элементов с, (/) = [clk (Z))^ и
c2(Z)= {с2И0)Г определено как
Т N
J (1.3)
причем при унитарном отображении Н, U-> С под-
пространствам Ht и Ut соответствует подпростран-
ство Ct всех функций c(s), из простран-
ства1 С, таких, что с($) = 0 при s>t (см. § 1 гл. I,
стр. 10).
Представьте теперь, что имеются два семейства
подпространств: Ht> t0<t<T, и Ut, to<t <Т, свя-
занных друг с другом следующим образом:
Ht = AUt, t0<t<T,
(1.4)
где А — некоторый линейный оператор из гильбертова
пространства U в гильбертово пространство Н. Или
представьте еще, что соотношение типа (1.4) связы-
вает (незамкнутые) подпространства Я? и (7?:
Я° = ЛС/?, (1.4')
замыкания которых есть
Ht^Hi Ut = lft, to<t<T.
Спрашивается, для каких операторов А можно утвер-
ждать, что соответствующие семейства Ht и Ut,
tQ<t < Т, будут изометричны?
Чтобы пояснить, какое отношение этот вопрос
имеет к теории обновляющих процессов, рассмотрим
несколько примеров.
Пусть г](/) = {ц£ (/)}Г, tQ<t < Т,—какой-либо «стан-
дартный» процесс, для которого известен тип обнов-
ляющего процесса или даже определены проекторы Qt
на подпространства ЯДт])> ^<1<Т>
Пусть £(/) = {£, (/)}{", tQ< t< Г, —другой процесс,
который мы желаем сравнить со стандартным про-
цессом п(/) = {т]х (/)}?, tQ<t<T, Введем оператор Л:
= tQ<t<T, (1.5)
42
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
и линейно продолжим его на линейную оболочку Я0 (ц)
всех значений г)г- (/), i = 1, ..., m, tQ < t < T (это воз-
можно, например, если указанные значения (t)
линейно независимы). Если обозначить Я?(т]) ли-
нейную оболочку значений щ (s), *’=1, •••> k <
< s /, то, очевидно,
И ____ ______________
Мы видим, что поставленный выше вопрос в отно-
шении оператора А и #, = //,(£), Ut = Ht(y^ озна-
чает следующее: при каких условиях на оператор А
обновляющий процесс для £(/), /0 < t < Г, имеет тот
же тип, что и обновляющий процесс для т](/), tQ <
< t< Т?
Отметим здесь, что для процессов одного и того
же типа важной является задача отыскания соот-
ветствующего изометрического оператора X*.
= t0<t<T
(см. (1.1)), с помощью которого проекторы Pt на
подпространства могут быть определены по
формуле (1.2):
Pt — XQtX~\ t0<t<T.
Рассмотрим несколько примеров, показывающих,
что структурный тип может меняться самым неожи-
данным образом при переходе от семейства подпро-
странств Ut, t0<t<T, к семейству Ht = AUt, /0 <
< t < Г, где А — ограниченный линейный опе-
ратор.
Пусть ц(/), 0 < t < 1, — стандартный винеровский
процесс. Как известно,
п(0 = 2паФаЮ, о</<1,
£=0
1
где па = J П (0 Фа (t) dt, а ср* (/) = sin (k + 1/2) nt — соб-
О
ственные функции ядра В (s, f) = min (s, /), 0 s, t 1.
§ И
ИЗОМЕТРИЧНЫЕ СЕМЕЙСТВА
43
Пусть А — оператор проектирования на конечномер-
ное подпространство L, порожденное величинами
(отметим, что т)^., £ = 0, 1, ..., есть орто-
гональный базис в пространстве £/ = //(£)). Рассмот-
рим процесс
п— 1
g (/) = аТ] (/) = 5 ад* (/), 0 < t < 1.
/г=0
Для любого t > 0 при надлежащем выборе точек
О < tj < ... < tn < t матрица {<р* (//)} будет невырож-
денной, и из уравнений
п—1
n*<Pfe (tj) = I (if), j=l, • • •, «,
находим, что величины т]э, ..., являются линей-
ными комбинациями значений gft), ..., l(tn) и при-
надлежат подпространству Очевидно,
0</< 1,
и процесс X (/) = {Xk (t)}” с компонентами Xk (/) = ти-н
k = 1, ..., п, будет обновляющим для g (/), 0 < t < 1.
Таким образом, если исходное семейство U t =
имело кратность 2V=1 (обновляющим
является сам процесс ц(/), 0 < / < 1, со структурным
типом Л), то семейство Ht = AUt, 0</<1, будет
иметь кратность М — п, а структурным типом (крат-
ности М) будет мера, целиком сосредоточенная
в точке t— 0.
Как показывает этот пример, даже для ограни-
ченного оператора А на месте точек непрерывности
семейства Ut (U/+0= Ut) у нового семейства Ht = AUt
могут появиться «скачки» (Нм =# Ht).
Вообще, для простейшего стандартного семейства
Ut = ЯДц), 0</< 1, где ц(/) — винеровский процесс,
можно указать такой проектор Л, что семей-
ство Ht = AUt, 0</<1, будет иметь произвольный
44
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
заранее заданный структурный тип. Именно, если
взять процесс вида
f / 1
О 'О
для которого обновляющий процесс имеет заданный
структурный тип (о существовании такого процесса g (/),
0<t < 1, с любым наперед заданным структурным
типом говорилось ранее на стр. 11), и взять орто-
гональный ему стандартный винеровский процесс
W(t), 0</<1 (Я (Г) ±//(g)), положив
n(/) = g(/) + lF(/), 0</<1,
то семейство Ht(г|), 0 < t < 1, будет изометрично
семейству Ut = 0</<1 (см. стр. 25). При
этом
£(/) = ЛШ 0</< 1,
где А — оператор проектирования на под-
пространство Н (g).
Вернемся к произвольным семействам подпрост-
ранств Ut и Ht = AUt, tQ<t<T (Л — ограниченный
линейный оператор).
Будем называть оператор А обратимым, если
существует ограниченный обратный оператор Л-1.
Для операторов такого типа вопрос об изометрич-
ности семейств Ut и Ht — AUt, tQ<t<T, легко ре-
шается в случае, когда Ut, tQ < t < Т, есть дискретная
цепочка подпространств, а именно,
Ut= ® (1.6)
tk<t
где tk, k = 1, 2, ..., — конечное или счетное мно-
жество точек интервала [/0, Г), в которых
= Utk + о © (Utk+h © Utk) ¥= 0.
Структурные типы
dFl(t)>dF2(t)>...>dFN(t)
§ И
ИЗОМЕТРИЧНЫЕ СЕМЕЙСТВА
45
такого семейства являются чисто дискретными, со
скачками в точках tk, k = l, 2, ... Напомним, что
Fj (t) = || Qtuj ||2, tQ < t < T, где Qt — оператор проек-
тирования на подпространство Ut, a щ, ..., uN —
полная система «циклических» векторов в гильбер-
товом пространстве U, для которых порождаемые
элементами Qtu^ tQ< t < Т, подпространства Uj орто-
гональны при различных и U=®Uf
/=i
(см. стр. 40). В рассматриваемом случае цикличе-
ские векторы щ, ..., uN могут быт выбраны следую-
щим образом:
Uj = 2 Ckuk}, j=l, N,
где ukj, j=\, , Nk, — ортонормированный базис
в соответствующем подпространстве Ак== Utk+0Q Utk
размерности Nk (Nk^N), ukj = Q при j>Nk и
2l Ck I2 <
k
Очевидно, семейства Ut и Ht — AUt, tQ<t<T,
будут изометричными тогда и только тогда, когда Ht,
t0<t<T, будет чисто дискретной цепочкой со скач-
ками в тех же самых точках tk, k = 1, 2, ..., причем
dim (Я^+о © Htk) = dim (Utk+o © (1.7)
при этом условии, выбрав любой унитарный
оператор X:
х(Uik+0k = 1, 2.......
будем иметь
Ht = XUt, t0<t<T.
Пусть A — ограниченный обратимый оператор
и Ht = AUt, Поскольку AUt+0 S AUt+h =
= Ht+h при всех h > 0, имеем AUt+b s f~) Ht+h = Ht+0.
h>0
Учитывая, что Ut = A~xHt, tQ< t < T, имеем A~'Hf+0 S
s Ut+0> Ht+Os AUt+0 и в итоге получаем
Hi+0 = AUt+0, t0<t<T. (1.8)
46
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. И
Равенство (1.8) показывает, что семейство Ht,
tQ< t < Т, является дискретной цепочкой подпро-
странств и может иметь скачки лишь в тех же точ-
ках tk, что и семейство Ut, /0 < / < Г, поскольку при
Ut+o = Ut имеем Ht+0 = AUt+0= AUt = Ht.
По поводу этого уместно заметить здесь, что без
условия обратимости оператора А, вообще говоря,
равенство (1.8) может не выполняться; более того,
что кажется несколько неожиданным для ограни-
ченного оператора Л, может оказаться, что f/z+0 =
= lim [//+л==0, тогда как Я/+о = lim AUt+h =# О
h-»+0 /i->+0
(ср. с примером на стр. 42).
Покажем, что выполняется условие (1.7). Дей-
ствительно, обратимый оператор А прямую сумму
подпространств Utk+o = UtkA~ переводит в пря-
мую же сумму Я^ + о = Htk + ЛА* (точнее, ни при
каком и е А* элемент х = Аи не принадлежит Ht^
откуда следует, что подпространство
Htk+oeHtk=(i~ptk) лд,
(Pt означает оператор проектирования на имеет
ту же размерность, что и подпространство ЛА*, ко-
торая совпадает для обратимого оператора Л с раз-
мерностью подпространства
Ай = Utk>
Таким образом, справедливо следующее предло-
жение.
Пусть Ut, —дискретная цепочка под-
пространств в гильбертовом пространстве U и Л —
ограниченный обратимый оператор из U в гильбер-
тово пространство Н\ тогда семейства Ut и Ht = AUh
t0<t < Т, изометричны *).
*) Существует гипотеза о том, что это верно для обра-
тимого оператора Лив случае произвольного семейства £//,
<t <Т,
§ 2] НЕКОТОРЫЕ МОДЕЛИ СЛУЧАЙНЫХ ПРОЦЕССОВ 47
§ 2. Некоторые модели случайных процессов.
Понятие регулярности и проблема факторизации
В дальнейшем нам удобнее будет рассматривать
обобщенные случайные процессы. Например, сравни-
вая «обычный» процесс g (/) = {gz (Z)}™, t0 < t < T,
с некоторым стандартным процессом ц (/) = {щ (/)}?*,
/0 < t < Т, можно ввести обобщенный процесс (g, и),
и е £7°, на предгильбертовом пространстве 47° — Н° (13)—
линейной оболочке всех значений т]/ (/), положив
(g, и) = Аи, (2.1)
где А— линейный оператор из 47° в Н (g), опреде-
ленный формулой (1.5) (а именно, Ax\t (/) = gz (/), i —
= 1, ..., m, /0 < t < T).
Представьте, что «наблюдается» некоторый обоб-
щенный случайный процесс (g, и), и е 47°, опреде-
ленный на некотором подпространстве 47° в гиль-
бертовом пространстве U со скалярным произведе-
нием (и, v), u, v U, в котором задано некоторое
семейство подпространств 4/?, /о < t < Г, такое, что
к моменту времени t «наблюдатель» располагает
всеми величинами (g, и), и е 47?. Предположим, что 47°
плотно в U '.U = UQ, и в гильбертовом пространстве
задан корреляционный оператор В:
Е (£, и) (g, и) = (Ви, v), и, v(= UQ.
Обозначим Ht замыкание подпространства всех
величин (g, и), и е 47?. Спрашивается, при каких
условиях на корреляционный оператор В семейства Ht,
tQ<t<T, и Ut = UQt, to < t < T, имеют один и тот
же тип, точнее, являются изометричными?
Если ввести оператор А, удовлетворяющий условию
А*А = В (2.2)
(например, можно взять Л = В1/2 — положительный
квадратный корень из положительного оператора В),
то будем иметь
Е (g, и) (g, v) = (Ви, v) = (Au, Av), и, U°,
48
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
откуда видно, что изометричными являются семей-
ства и Alft, to<t<T, и поставленный выше
вопрос можно переформулировать следующим обра-
зом: при каких условиях на оператор А будут изо-
метричными семейства
tft = AU° и Ut = Ul h<t<T
(2.3)
(ср. (1.4) и далее)?
Вообще, рассматривая вопрос о типе семейства
подпространств //,(£), tQ < t < Г, возникающих при
«наблюдении» обычного случайного процесса £(/),
/0 < t < Т, или сравнивая /7z(g), t$<t<T, с некото-
рым другим «стандартным» семейством (7Z, t0 < t < Т,
в том или ином гильбертовом пространстве £7, как
правило, можно перейти к описанной выше схеме
обобщенного процесса (£, и), и е t/°, гильбертовом
пространстве U с заданным семейством подпространств
£7?, to < t < Г, таких, что совпадает с замы-
канием подпространства всех величин (g, и), и s [7?,
a Ut = UQt, tQ<t< Т.
Проведем эту простую редукцию для беско-
нечномерного процесса, заданного таким обра-
зом, что его компоненты, обозначаемые {£(/)> х},
to<t < Т, отмечены «индексом» хе /?, где 7?— сепа-
рабельное гильбертово пространство со скалярным
произведением {х, у}, х, у R. Будем предполагать
при этом, что корреляционная функция B(t, s), tQ <
< t, s < T, такого процесса, определяемая из соот-
ношения
ВО), *}U(s), = s)x, у}, x,y<=R, (2.4)
является слабо непрерывной операторной функцией
в R.
Введем линейное пространство £7° всех функций
со значениями в /?, являющихся линейными комби-
нациями вида
и (/) = 2 Xkck (/), to < t < Т,
k
§ 2]
НЕКОТОРЫЕ МОДЕЛИ СЛУЧАЙНЫХ ПРОЦЕССОВ
49
где коэффициенты xk есть элементы из гильбертова
пространства R, а числовые функции ck(f) финитны
(обращаются в 0 вне конечного интервала) и удо-
т
влетворяют условию J | щ,(/) |2Л < оо. Обозначим [7?
^0
подпространство, образованное всеми функциями
u = u(s), tQ<s<T, указанного типа, удовлетворяю-
щими дополнительному условию u(s) = 0 при 5>Л
Введем также гильбертово пространство U = L2(R)
всех измеримых функций u = u(t), со зна-
т
пениями в /?, для которых J || u(t)\\2dt < оо, опреде-
лив скалярное произведение элементов и, v е U как
т
(и, у) = | {и (/), v (/)} dt. (2.5)
to
Очевидно, замыкание if} в U представляет собой
подпространство Ut всех функций w = w(s),
из U — UQ, удовлетворяющих условию
и ($) = 0 при s>t, (2.6)
и структура семейства подпространств U(, t^<t<T,
может быть принята за определенный стандарт.
Введем обобщенный случайный процесс (g, и),
u^U, положив
т т
(I, ы) = / а (0.«(/)} dt=2 J а ю, xk) мо dt (2.7)
to k to
при и (0 = 2 xkck (f) e U°.
k
Очевидно, замкнутая линейная оболочка //z(g)
всех значений {£($), х}, tQ< s^t, совпадает
с замыканием подпространства //?(£) всех величин
(g, u), (напомним, что мы считаем корреля-
ционную функцию B(s, t) слабо непрерывной, а это
50
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
равносильно непрерывности в среднем квадратичном
всех компонент {g(/), х}, t0<t< Т, исходного процесса).
Далее, из формул (2.4), (2.7) легко вывести, что
Т Т
= j’ J {В (s, f) и (s), v (/)} ds dt, u,veUr\ (2.8)
to to
и если предположить выполненным условие
т т
J J || В (s, t)ll2dsdt< оо,
to t0
(2.9)
то обобщенный процесс (g, и), и е (/°, по непре-
рывности может быть продолжен на все гильбертово
пространство U = L2(R), а ограниченный положи-
тельный оператор В с ядром B(s, /):
т
Ви (t)= j В (/, s) u(s)ds, u^U, (2.10)
to
будет корреляционным для % (и), u^U.
Как уже отмечалось выше (см. (2.3)), вопрос о том,
будет ли семейство Я/(£), tQ<t< Г, изометрично стан-
дартному семейству Uti tQ<t < Т, сводится к анало-
гичному вопросу для семейства Ht = AUh tQ<t<T,
где A = Bi/2 — квадратный корень из положительного
оператора В.
Остановимся еще на одном важном случае —
бесконечномерном стационарном процессе g(/),
— оо </<оо, с компонентами {£(/), х}, x^R, которые
являются стационарно связанными в обычном
смысле процессами. Условие стационарности состоит
в том, что корреляционная функция В(/, $)— см. (2.4) —
зависит лишь от разности t — s.
Предположим, что имеется спектральная плот-
ность fKi —оо<Л<оо,—измеримая операторная
функция в гильбертовом пространстве R (при фикси-
§ 2] НЕКОТОРЫЕ МОДЕЛИ СЛУЧАЙНЫХ ПРОЦЕССОВ 51
рованном А представляющая собой положительный
оператор fK в R), такая, что
Е {£ (0, х} {g (s), у} = {В (t — s) х, у} =
= J elK»-s>{fKx, y}dk, x,yt=R. (2.11)
—оо
Предположим для простоты, что функция fK является
ограниченной (||/%|]< С). Обратившись к соот-
ветствующему обобщенному процессу (g, и),
(в формуле (2.7) нужно положить /0 =— оо и Т == оо),
будем иметь
оо
E(g, н)(П0= /О(Л), 6(Z)}dA, (2.12)
—оо
где й (Л) обозначает преобразование Фурье соответ-
ствующей функции и е £7°, а именно,
й (Л) = J еши (t) dt =
—оо
оо
= xk j* eiUck(i)dty —оо < Л < оо, (2.13)
k — ОО
При и (/) = S ХкСк (t), —оо <t < оо.
k
Удобнее перейти к новому параметрическому про-
странству, взяв вместо £/° пространство UQ всех функ-
ций й(к), которые связаны с функциями u(f)^UQ
с преобразованием Фурье, и положив
(g, u) = (g, и), йе=й°. (2.14)
Из формулы (2.12) видно, что обобщенный процесс
(g, й), й е £7°, имеет своим корреляционным операто-
ром «оператор умножения» на спектральную плот-
ность fK:
[Вй](к) = /\й(Л), —оо<Л<оо. ' (2.15)
52
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
Обозначим U°t подпространство в пространстве 0й,
образованное всеми элементами й е U0, которые по-
лучаются преобразованием Фурье (см. 2.13)) функций
u = u(t) из соответствующего подпространства U°t,
-- ОО < t < оо.
Как известно, преобразование Фурье есть унитар-
ный оператор в L2(7?), и, следовательно, семейства С/?
и Я?, —оо <t< оо, унитарно изоморфны. Как пока-
зывает формула (2.12), унитарно изоморфными
будут также семейства подпространств Я?(£) и ЛЯ?,
— оо < t < оо, где А — Вх'2 есть «оператор умноже-
ния» на fj/2:
[лвд=^2ад,
— оо < Л < оо.
(2.16)
Таким образом, вопрос о том, при каких условиях
на спектральную плотность fK семейство ЯД£) будет
того же типа, что и описанное выше стандартное
семейство U t, — оо < / < оо (см. (2.6)), является част-
ным случаем общего вопроса об изометричности
семейств вида ЛЯ?, t0 < t < Г, и Я t = UQv tQ<t<T,
обобщенного процесса (g, и), и е Я0, на пред-
гильбертовом пространстве Я0 с корреляционным
оператором В = Л*Л.
Обратимся к произвольному обобщенному про-
цессу (g, и), UQ, с корреляционным оператором В
в гильбертовом пространстве Я = Я0 с заданным
семейством подпространств Я?, /о < t < Г, и рассмо-
трим соответствующее семейство ЯД£), /0 < t < Т
(напомним, что ЯД£) обозначает замкнутую линей-
ную оболочку всех величин g (и), и е Я?). Положим
Ut = lPt и будем считать, что пространство Я совпа-
дает с замыканием объединения всех подпро-
странств UtJ tQ<t <Т.
Назовем процесс (g, и), и е Я0, регулярным, если
семейства Я,^) и Ut, tQ<t<T, изометричны:
tQ<t<T,
(2.17)
§ 2]
НЕКОТОРЫЕ МОДЕЛИ СЛУЧАЙНЫХ ПРОЦЕССОВ
53
где X — некоторый унитарный оператор из U
в пространство Н (g) — замкнутую линейную оболочку
всех величин g(u), и е 77°.
Будем говорить, что оператор В допускает соб-
ственную факторизацию, если он представим в виде
В = С*С, (2.18)
где С—некоторый линейный оператор из U° в U
такой, что
CU*=Ut, tQ<t<T.
(2.19)
Теорема. Процесс (g, и), u^U°, является регу-
лярным тогда и только тогда, когда его корреляцион-
ный оператор В допускает факторизацию.
Доказательство. Как уже отмечалось, семей-
ство Ht(Q, tQ<t<T, изометрично любому семейству
подпространств вида AU°f, tQ<t<T, где А — линей-
ный оператор, определенный на UQ и удовлетворяю-
щий условию А*А = В. Поэтому если процесс (§, и),
и е 77°, регулярен, то семейство подпространств АЩ,
tQ<t < Т, где А = В1/2 — положительный квадратный
корень из корреляционного оператора В, изометрично
семейству Ut = UQt, t^<t < Т, т. е. существует уни-
тарный оператор X из U = UQ в AUQ, такой, что
AUQt = XUt, tQ<t<T. Очевидно, оператор С = Х~1А
удовлетворяет условиям (2.18), (2.19).
С другой стороны, если корреляционный опера-
тор В допускает факторизацию (т. е. существует
некоторый оператор С, удовлетворяющий условиям
(2.18), (2.19)), то подпространства AUQt при А = С
просто совпадают с Ut, а как только что было
доказано, регулярность равносильна изометричности
семейств AUQt и Ut,tQ<t<T (для какого-либо А,
А*А —В). Следовательно, если корреляционный опе-
ратор В допускает собственную факторизацию^ то
процесс (£, и), и е £7°, является регулярным.
54
регулярные Процессы
[ГЛ. п
Как показывают приведенные ранее примеры,
отнюдь не всякий оператор В допускает фактори-
зацию.
По поводу факторизации отметим также следую-
щее.
Предположим, что положительный оператор В
в гильбертовом пространстве U представим в виде
В = С*С. Такое представление, конечно, не един-
ственно; например, имеем также В — С\-С\ для
С1 = ХС, где X — произвольный унитарный опе-
ратор в U. С другой стороны, если для каких-либо С
и Cj выполняется условие В = С*С = C\Ci, то опе-
ратор X, определенный на подпространстве CU° ра-
венством ХСи = Сщ, u^UQ, будет и з о м е тр и че-
с к и м:
(ХСи, XCu) = (CJu, С^)=:
= (С1Сщ, v) = (С*Си, v) = (Cu, Cv), и, v^U.
При условии же, когда Си] = Схи^= Ut для всех t
(см. (2.19)) имеем CU° — C}UQ = U, и, следовательно,
оператор X однозначно продолжается до унитарного
оператора на всем пространстве U. Кроме того, по-
скольку для любого u^Ut найдется последователь-
ность Un^U°t, п=\, 2, ..., такая, что и = lim Сип, то
Хи = lim ХСип = lim Сщп е Ut.
Таким образом, при условии (2.19) факторизую-
щий оператор С в соотношении (2.18) определен одно-
значно с точностью до унитарного множителя X,
удовлетворяющего условию
xut=ut, t0<t<T.
(2.20)
Для заданного семейства t/0, t0 < t < Г, со струк-
турными типами dFx (/)>... ";>dFN (t) можно исполь-
зовать следующую стандартную модель: U—простран-
ство векторных функций и (/) = {uk (/)}^, /0 < t < Г,
с компонентами uk(t),
uk (0 = 0 при dFk (t)/dF{ (/) = 0, k = 1, ..., N,
§ 2]
НЕКОТОРЫЕ МОДЕЛИ СЛУЧАЙНЫХ ПРОЦЕССОВ
55
в котором скалярное произведение элементов u(t) =
= {uA,(/)}^ и v (0 = {vk (01Г определено формулой
(«, <0 =
Uk (О Vk (О
dFx(fy,
семейство Ut, tQ<.t<T, в функциональном простран-
стве U реализуется как семейство подпространств
Ut = {u(s), t0 < s < Т: &($) = 0 при s > t}
(см. по этому поводу § 2 гл. I). Обозначим Е TV-мерное
векторное пространство с обычным скалярным произ-
ведением S xktfk векторов х = {xjf, у = {ykW е Е и
k=\
£z — подпространство в нем, образованное векторами
х = с компонентами
Xk — 0 при dFk (t)ldFx (/) = 0.
Пусть X — унитарный оператор в U, удовлетво-
ряющий условию (2.20). Тогда при всех t
t N
J '^i(Xu)k(s)-(XvK(sjdFl (s) =
to k~l
t
= j uk (s) • Vk ($) dF{ (s), u, V (= U.
t0
Видно, что при фиксированном s отображение
{^а<->{(хи)и<,
определенное для почти всех s, является унитарным
в подпространстве Es Е. Отсюда уже легко вывести
формулу, описывающую общий вид унитарного опе-
ратора . X, который удовлетворяет условию (2.20)
в функциональном пространстве U:
(Xu)(f) = Xtu(t), tQ<t<T, (2.21)
где Xt при фиксированном t есть произвольный уни-
тарный оператор в подпространстве Et векторного
пространства Е (необходимо только, чтобы оператор-
ная функция Хр t$<t <Т, была измеримой),
56
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
§ 3. Одна теорема о факторизации
Нашей целью будет здесь доказательство следую-
щего предложения *).
Теорема. Положительный (обратимый) опера-
тор В в гильбертовом пространстве U вида В = 1 -— G,
где G — оператор Гильберта — Шмидта, допускает
собственную факторизацию
В = С*С (CUt = Ut, tQ<t<T) (3.1)
относительно любого семейства подпространств Ut,
t0<t<T.
Предварительно отметим, что оператор В допускает
факторизацию (3.1), если он представим в виде
В = С_ • С+, (3.2)
где операторы С_ и С+ являются обратимыми и
удовлетворяют условиям
C+Ut = Ut, C-Ut- S Uh, to < t < T;
здесь «сопряженное» семейство подпространств
u^-=ueut,
to < t < T,
инвариантно относительно оператора C_ тогда и
только тогда, когда само семейство Ut, tQ<t<T,
инвариантно относительно сопряженного оператора
CL : CLUt <= Ut.
Действительно, в случае положительного опера-
тора В положительным будет также оператор
d = (с;1)’ вс+1 = (с’+) 1 вс+1 =
=(с’+)’|с,+с:с;,=с’_с;1,
для которого, очевидно, DUt^Ut. Поскольку опе-
ратор D является обратимым, на самом деле имеем
*) Ср. с изложенным по этому поводу в книге И. Ц. Г о х-
берга и М. Г. Крейна, Теория вольтерровых операторов
в гильбертовом пространстве и ее приложения, М., «Наука»,
1967 (теоремы 6.2 гл. IV и 10.1 гл. I).
ОДНА ТЕОРЕМА О ФАКТОРИЗАЦИИ
57
§ 3]
DUt=Ut. Если ввести оператор £>1/2—положитель-
ный квадратный корень из D, то получим, что
В = С*+£>С+ = С*С,
где оператор С — Dil2C+ удовлетворяет условию
< cut — ut, t0<t< т.
Мы покажем далее, что обратимый оператор
В = I — G можно факторизовать с помощью операто-
ров С+ и с_ = вс;‘ (см. (3.2)), взяв С+= (/+ G+)~‘,
где оператор G+ определен как своего рода интеграл
G+ = / / (/ - QtGQt-)-' dQsG dQt (3.3)
s<t
по «операторной мере» dQsG dQt (здесь Qt — проек-
торы на подпространства Ut, t0<t<T).
Отметим’ сразу, что поскольку В — положительный
обратимый оператор, то
sup (Gu, и) = 1 — inf (Ви, u) — r< 1
II «11=1 II «11=1
и для любого проектора Q
inf ([/ — QGQ] и, и) — 1 — sup (QGQu, и)^ 1 — г > О,
Цы||=1 И«11=1
и следовательно, I — QGQ является положительным
обратимым оператором, причем
II7 — QGQ 1Г‘<. (3.4)
Выберем конечное разбиение у = {/о<Л<
... < tn < tn+\ — T}. Положив AQ/fe = Q/fe+1—Qtk, всякий
оператор А можно представить в виде
Л =( 2 AQfjX ( 5 AQf,) =
U=0 \f=o 4
= 2 AQ/ XAQ/ + 2 AQ/.XAQ/.,
k<i R 1 R 1
где
5(Л)= s AQ^AQz.
k<j R '
58
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
обладает тем свойством, что
п
Utk =
Аг-1
= 2 (QtjA &Qtj) Utk S utk_t, a = 1,..., n 4-1, (3.5)
и
S(X) Ut, = 0,
а дополнительная часть S'(X) = X — S(X),
S'M)= 2 AQi/AQ/.,
такова, что
n
S'WUi=^Q AQ =
n
= 2 (AQtAQf\ ut = ut, j = O,...,n. (3.6)
R.=J x К K~r I / J ]
Посмотрим, что дает «операция усечения» S(A)
для оператора
X = G(/ + G+),
где G = I — В и
G+ = 2 (/ - QtjGQtjy' &QtkG SQtj. (3.7)
Имеем
5(G-G+)=2Q//(G-G+)AQ,/ =
n n
= 2 Q,,c 2 (l - ()-’ Q,fi AQ,t AQ,t =
= 2 Q,tC II - QitGQlty’ Q,fi AQ,t,
а поскольку каждое из подпространств t/^ инвариантно
относительно операторов I — QtfGQtk, (/ — QtkGQtkyl
и проектор Qtk перестановочен с ними, то
s(G • G+) = 2 QiaGQ<a(/ -QtkGQtkY' Qt/G\Qt
& I fv /V \ /V /V/ /V гм
§ 3] ОДНА ТЕОРЕМА О ФАКТОРИЗАЦИИ 59
В итоге получаем
S[G(Z4-G+)] = S(G)4-S(G-G+) =
п
= 2 [' + Q'fiQ',. (' - <?<,««<,)I Q'fi -
= 2 у -QifiQltr' Gift ЭД, = <;+
। \ /V К/ К K>
Таким образом, оператор G+ получается «операцией
усечения» оператора А — G (Z + G+):
S[G(/ + G+)] = G+.
Легко проверить, что
5(ЛГ = 0 (3.8)
(ср. с (3.5)), и следовательно, спектр усеченного опе-
ратора G+ = S (Л) состоит из единственной точки
Л = 0, так что оператор (/ + G+) является обратимым.
Если взять G~ = S'(Л) для X = G(/ + G+), то будем
иметь
G+ + G^ = G(/ + GJ
и
(7-G)(/+GJ = / + G+--G(/ + G+)==
= /+ G+ —G+ —G_ = / —G_,
откуда находим
(Z —G) = (Z —G_)(Z+G+)-1. (3.9)
Подпространства Utk, k = 0, 1, ..., n + 1, инвариантны
относительно усеченного оператора G+ = S [G (I + G+)]
(cm. (3.5)), а следовательно, и относительно опера-
торов Z4-G+, (Z + G+)-1. В свою очередь, подпро-
странства Ut~k, k = 0, 1, ..., п+1, инвариантны
относительно операторов G_ = S' [G (Z + G+)] и Z — G_
(см. (3.6)). Таким образом, соотношение (3.9) дает
нам факторизацию оператора B — I — G с помощью
С+ — (I + G+)-1 и С_=1—G_ относительно дискрет-
ной цепочки Utk, t0 < tt < ... < tk+i — T (см. (3.2)).
60
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. и
Выражение (3.7) для оператора G+ = G+* напоми-
нает своего рода интегральную сумму, и если пред-
ставить себе, что при все более мелком разбиении
V —{/о<Л < ••• <tn<tn+i — T} существует пре-
дельный оператор G+ == lim G(+\ то его естественно
П-»оо
было бы обозначить так, как это сделано в фор-
муле (3.3). Очевидно, если существует предельный
(в смысле сильной сходимости) ограниченный опера-
тор G+, такой, что Z4-G+ имеет обратный оператор
C+ = (/ + G+)~1’ то существует также предельный
оператор
С_ = lim В (/ + G(;0 = В (Z + G+)
П->оо
и равенство В = С_ • С+ дает факторизацию отно-
сительно полного семейства Ut, t0< t < Т.
В самом деле, всякий элемент u<=Ut есть предел
элементов Uk^Utk при tk->t — 0, и потому
G+uk = lim G{t^uk е Ut. s Ut, G+u = lim G+uk s U(,
П->оо t^t—O
так что каждое из подпространств Ut инвариантно
относительно операторов I-\-G+, (/4-G+)-1. Анало-
гичные рассуждения применимы к оператору С_ =
= lim С{-, где С(- = В (/ + G+*) удовлетворяет уело-
П->оо
вию C{-Ut~k — UДля любых tk из соответствующего
разбиения = {/0 < tx < ... < tn < Т}.
Покажем теперь, что существует равномерный
предел limG(+) = G+ в случае, когда
G = l — В
есть оператор Гильберта — Шмидта с конечной
«следовой нормой»:
| G |2 = Sр G*G = 21| Gup |р = 51 (Gup, uq) |2 < oo (3.10)
P P,q
(где щ, u2, ... —какой-либо ортонормированный базис
в гильбертовом пространстве G), причем, так же как
§ 31
ОДНА ТЕОРЕМА О ФАКТОРИЗАЦИИ
61
и операторы G(+, предельный оператор G+ имеет
лишь единственную точку спектра Л = 0, так что
оператор (Z4-G+) будет обратимым.
Для любых разбиений —{/0 < s1 < ... < sm <
< sm+1 < Т} и y(rt) = {/0 < . <tn<tn+l==T}, таких
что г у(п) (каждый интервал разбиения [$г, s/+I)
совпадает с объединением некоторых интервалов
k/> ^/+i)), разность G{+—G^1’ можно представить
в виде
G(n) _ = 2 _ Fs Q
где Ft = (l — QtGQt) a $(/) обозначает левый конец
интервала из разбиения у(,п\ содержащего точку
i/<=yM (Si<tz< su.! при s (j)=- st), —см. фор-
мулу (3.7).
Заметим, что для любых операторов Л,
12 Afi \Qt. |2 = 2 | Afi\Q{j |2. (3.11)
Это равенство легко получить, если выбрать полные
ортонормированные системы в ортогональных
подпространствах = UtJ+1 Q Utj; j = 0, 1, ..., n,
и взять их объединение щ, u2j ... за базис во всем
пространстве U. Действительно, A}G \Qt.up = 0 при
j ф j (р), если ир е AQZ/ (р), откуда следует, что
[2 Afi \Qtjup||2 = 21| Afi \Qt.up||2
и
A,G &Qtj |2 = 2 j 2 A}G \Qtjup |2 =
= 2 2 И/G &Qtjup ||2 = 2 [21| AG aq<7«p II2] =
= 21 AjG &Qt] I2.
62
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
Кроме того, для любых ограниченных Alt А2 спра-
ведливы следующие соотношения:
] A,GA2]2 = H\(A1GA2up, uq)? =
= 21 (A'2GAiUq, ир)? = [ A*2GAi |2
P> q
|Д,ОД2|2=2|| Д^Д^ДС
< IIЛ Ip • 2II GA,up IP = II A! II2 • IGД212 =
= ||Д1|р.|Д5О|2<|| Д, |p-|| АИР -1 G|2 =
= || AJp-H A2[p-|G|2. (3.12)
Используя соотношения (3.11), (3.12), легко полу-
чаем, что
I G<;' - G(r г = 2 I (F4Qtj - Fs(i)QS(i}) G AQt{12 =
= 2 | (Ft -Fs (1))QS (I)G bQt/+Ftj (Qt-QS(n) G |2<
<2 2 (|| (Ftf - Fs (>>) J)21| Qs (n II21G \Qt/12 +
+ ||FZ/||2|AQs(/)GAQZ/|2,
где AQH/) = Qtj—Qs (,)(s (j)<*/)- Поскольку || Qs (/) ||= 1
и \\Ft 5^ t 2. r (см. неравенство (3.4)), имеем
I G(;> - G(+m) I2 < 2 2 (I Ftj - Fs (/) ||21G AQZ/12 Д-
/=1
+ 7r^7F SlA^</)GAQ//|2- (3.13)
/=1
Согласно общей формуле (3.11) функция Л(Д/) =
= |GAQJ2 на полуинтервалах Д/ = [/, s) является
аддитивной. Выбрав собственный для оператора G
§3]
ОДНА ТЕОРЕМА О ФАКТОРИЗАЦИИ
63
ортонормированный базис щ, и2, ... \Gup = kpup,
JjAp < °°У имеем
р /
| G \Qt |2 = | \QtG |2 = 2II &QtGup IP =
p
= ^lK2p\\NQ(up\^Q при /->s, (3.14)
p
поскольку проекторы Qs непрерывны слева: AQZ =
= Q5 — Qt~>0 при t->s, и для любого фиксирован-
ного элемента u^U его проекция AQziz на под-
пространство АС/1= USQUt такова, что || &Qsu || —> 0.
Таким образом,
Л (АО = | G \Qt |2
является не только аддитивной, но и непрерывной
функцией множеств А/ = [/, s) и, следовательно, может
быть продолжена до ограниченной борелевской меры
на интервале [/0, Т).
Очевидно, такими же свойствами обладает и функ-
ция множеств
|i(AsXA/) = |AQsGAQJ2,
определенная на прямоугольниках AsXA^ = [si> s2)X
X [G, G)‘, именно, из общей формулы (3.11) вытекает,
что ц (As ХА/) аддитивна, а из оценки (3.14) — что
она непрерывна, и следовательно, продолжается до
ограниченной борелевской меры на квадрате [/0, Г)Х
Xko, Т).
Положим f(t) = ц {[/0, t) X [Лъ t)}> tQ<t<T. Для
операторов Ft = (I — QtGQt)~l имеем
F;'-Frl = QsGQs-QtGQt,
|F71-FrI|2 = f(0-f(s),
I Ft - Fs |2 = | F, (F71 - Ff1) Fs |2 <
< (147)2-1 FT1 - FJ2 = -^4^ [f (0 - f (s)].
Видно, что первый член в правой части неравенства
64
РЕГУЛЯРНЫЕ ПРОЦЕССЫ
[ГЛ. II
(3.13) удовлетворяет условию
п п
SII Ftj -Ps(n ||21G \Qt. |2 < If (//) - f (s (/))] =
/=1 /=1
~ — г)2 Im)’
m n
где lm = 2 f (s (i)) A, In = 2 f (^/) к (МЛ есть «лево-
го /=0
интегральные суммы» для неубывающей функции
f(t),tQ<t<T. Очевидно, последовательность интег-
ральных сумм In, п = 1,2, ..., является ограниченной
и неубывающей (при все более мелких разбиениях),
так что существует lim 1п и, следовательно, при
П->оо
т, п -> оо
п
2 и<, - А I’ I са<г,, |= < (/„ -1 „) -» 0.
/=1
Далее, введем множество
Ал ==={«, /: s < t < s + h}.
п
Очевидно, объединение (J [s (/), О Х[0> О+i) содер-
/=i
жится в Алпри/г = max] sz+1 — s, |, когда/,-+1 — з(/)</г.
Следовательно, для второй суммы в правой части не-
равенства (3.13) справедлива оценка
2 I AQS (f)G AQt I2 < ц (АЛ) -> 0 при Л->0.
/==11 71
Из полученных оценок вытекает, что для последо-
вательности вложенных друг в друга разбиений
п = 1, 2, ..соответствующие операторы G{± таковы,
что при п, т->оо, и следова-
тельно, существует предельный оператор О+:
(О^ — О+|2->0 при п->оо. (3.15)
Нам осталось показать, что предельный опера-
тор G+ имеет, как и операторы G™ = S(G(l + <?'+))
(см. (3.8)), единственную точку спектра Л = 0.
§3]
ОДНА ТЕОРЕМА О ФАКТОРИЗАЦИИ
65
Являясь вполне непрерывным оператором (и даже
оператором Гильберта — Шмидта), G+ может иметь
лишь не более счетного числа точек спектра, причем
предельной для них может быть лишь точка Л = 0.
Таким образом, всякая точка Z 0 имеет окрест-
ность, за вычетом X целиком принадлежащую резоль-
вентному множеству оператора G+, и на любом
замкнутом контуре Г этой окрестности определена
резольвента RZ = (G+— zl)~\ которая непрерывна и
ограничена на Г, При условии равномерной схо-
димости G^~>G для всех достаточно больших п
№ = (G(;’ - z/)"‘ = Rz[l - Rz (G+ - =
= Rz J Rkz(G+-G^k
k—Q
и R(z}->Rz при n->oo равномерно по геГ, так что
zkRz dz = lim f
ГС~>ОО
zkR{z} dz = Q при всех k 0,
поскольку единственной особой точкой аналитических
функций R{z} является z = 0 и j zkRzn) dz = 0. Таким
г
образом, J zkRzdz = 0, что означает аналитичность Rz
г
в точке Л. Следовательно, единственной точкой спек-
тра оператора G+ является /1 = 0.
В итоге мы установили, что положительный обра-
тимый оператор B — I — G, где G есть оператор
Гильберта — Шмидта (см. условие (3.10)), допускает
представление В = С~ • С+ типа (3.2); в качестве С+
можно взять оператор С+ = (/4-С+)-’1, где G+ = lim G^
есть предел по «следовой норме» операторов Гиль-
берта — Шмидта G+\ определенных формулой (3.7).
Отсюда, как уже отмечалось, вытекает возможность
факторизации относительного любого семейства под-
пространств U(,
3 Ю. А. Розанов
ГЛАВА III
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
§ 1. Структурный тип регулярного
стационарного процесса
Обратимся к описанной ранее модели стационар-
ного процесса £(/), — оо < t < сю, со стационарно
связанными компонентами {| (/), х}, x^R, каждая
из которых соответствует некоторому х е R, где
R — сепарабельное гильбертово пространство со ска-
лярным произведением {х, у}, х, у е R, Будем счи-
тать, что существует спектральная плотность fK,
— оо <Д < оо, —- измеримая операторная функция
в гильбертовом пространстве R (при фиксирован-
ном Л представляющая собой положительный опера-
тор в А?) такая, что
со
Е {£ (0» • {s («), У} = J (i~s) {f\x, у} dk =
iK(t-s) fifty} X,IJ<=R, (1.1)
где f^2 обозначает «квадратный корень» из точ-
нее, /J/2, — оо < Л < оо, — операторная функция, удо-
влетворяющая условию (f}/2)2==/л-
Введем гильбертово пространство L2(R) интегри-
руемых в квадрате функций х(Л), —оо < Л < оо, со
СТРУКТУРНЫЙ тип
67
§ 1]
значениями в /?, определив скалярное произведение
в L2(R) равенством
(*, у) = / {X (Л), у (Л)} dK, х, у ^L2 (/?).
—оо
Формула (1.1) показывает, что отображение
{g (/), х} еа71/2*, — оо</<оо, (1.2)
является изометрией из //(g) в L2(R). Поэтому,
рассматривая такие вопросы, как, например, вопрос
о структуре семейства подпространств 77z(g), — оо <
< t < оо (Z/z(g) есть замкнутая линейная оболочка
всех значений {g (s), х}, х е /?; можно вместо
самого стационарного процесса
{g(0, х}, хе/?; — оо</<оо,
использовать изометричную ему функциональную
модель
^Ш^/2Х, — ОО < Л < оо}, хе /?; — оо < / < оо,
в гильбертовом пространстве А2(/?). В этом смысле
семейство //Jg), — оо</<оо, имеет ту же струк-
туру, что и семейство подпространств
Ht = V e^f^R, - оо < / < оо (1.3)
s
(Ht есть замкнутая линейная оболочка подпространств
eiKsfl/2j^ Е [2 (fy, И Т. П. Положим
Я = V e^f^R. (1.4)
—оо < s < ОО
Обозначим Pt оператор проектирования в про-
странстве Н на подпространство Hf и Et — унитар-
ный оператор умножения на функцию е1М от Л,
— оо < А, < оо.
Очевидно,
V eiKis+tV2X = Ht+a
s^t
(1.5)
3*
68
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. Ш
при всех и, t. Поскольку оператор Еи унитарный,
также имеем EuHi = Ht+U (где Hi = HQ Ht обозна-
чает ортогональное дополнение к подпространству Н^.
Таким образом, для любого элемента хе Ht+u имеем
Е—цХ s Н0 Н/Е—цХ — Е—цХ, ЕиРtE_ux — х,
а для хе Hi+u
E_ux^Hi, PtE-ux = 0, EuPtE-ux = 0,
откуда видно, что произведение EuPtE_u совпадает
с проектором Pt+u- Учитывая, что Е_и = Ей\ полу-
чаем следующее равенство:
Pt+u = EuPtEul, —oo<t<oo. (1.6)
Оно, по существу, равносильно условию (1.5), по-
скольку при всех t, и
Pt+u(EuHt-) = Eu (PtHi) = 0
И
EuHi S Hi+U, Hi = Е-и (EuHi) S E.uHi+u,
откуда (заменив параметры t, и на t-\-u, —и) полу-
чаем также, что Hi+U Е EuHi, и следовательно,
EuHt — Ht+u, EuHt — Ht+u.
Соотношение (1.6) показывает, что семейства Н,
и Ht+u> — оо < / < оо, изометричны при любом ut
— ОО < U < оо (см. § 1 ГЛ. II).
Ясно, что если Fj (t) > F2 (f) > ... >FM (/) — упоря-
доченная система структурных функций для семей-
ства проекторов Pt, —оо </<оо, то F\(t—и)>-
>-F2(/—-м))> ... — «) будет аналогичной си-
стемой для семейства Pt+U, — 00 < t < °°, и поскольку
семейства Ht и Ht+u, —оо < / < оо, изометричны,
структурные функции Ff(t) и Ff(t — и) должны быть
одного и того же типа:
dFt(t)~ dF^t — u), /=1, .... Al, (1.7)
при всех и, — оо < и < оо.
СТРУКТУРНЫЙ тип
69
§ 11
Свойство (1.7) обычно называют квазинвариант*
ностью относительно сдвигов. Если рассматривать
наши меры dF/(t) лишь на прямой —оо </< оо (без
точки t0 —— оо), то, как известно, при условии ква-
зиинвариантности
dFjify'^dt (1.8)
(отметим, что значение F, (/0) = lim Fj(t) может быть
/-> —00
положительным).
Для удобства читателя приведем здесь простое
доказательство того, что всякая квазиинвариантная
мера и эквивалентна лебеговой мере. Если взять
свертку ограниченной квазиинвариантной меры ц
с гауссовским распределением Р, имеющим плотность
p(t) = JLr е~/2/2, то окажется, что
У 2л
Ц * р (Д) = | р, (Д — t) р (/) dt = О
тогда и только тогда, когда р (Д) = 0, т. е. что р * Р~ р.
Но свертка р * Р имеет преобразование Фурье вида
ср (X) e~v/2, где ср(^) — преобразование Фурье ограни-
ченной меры р, а следовательно, р * Р имеет гладкую
плотность, которая совпадает с обратным преобразо-
ванием Фурье функции ф(Л)в“Л2/2, —оо < X, < оо. По-
скольку мера р*Р является также квазиинвариантной,
то ее (гладкая) плотность не должна обращаться в О
ни на каком интервале, откуда следует, что эта
плотность положительна почти всюду.
Напомним, что в свое время мы приняли за стан-
дарт структурный тип семейства подпространств
Ut = Lt (Rm), — оо < t < оо, в функциональном про-
странстве U = L2(RM), где RM — некоторое гильбер-
тово пространство размерности М (Lt (Rm) есть под-
пространство всех функций u(s), —оо<$<оо, из
L2(PM), обращающихся в 0 при s>t).
Остановимся подробнее на структурном типе
семейства подпространств Ut = L2i(Rm), —оо </< оо.
Обозначим Q/ проекторы на Ut. Функция u(s),
— оо < $< оо, как элемент гильбертова пространства
70
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
t/z = L2(/?M), имеет своей структурной функцией
t
F(t) — J || tz(s)||2ds, — оо < t < co. Таким образом,
т-ОО
всякий элемент и е U такой, что и (s) #= 0 почти
всюду, имеет максимальный лебеговский тип:
dF(t)~dt, Очевидно, элементы
Uj (s) == c(s) —оо < s < oo; j= 1, ..., M
(где числовая функция c(s), —оо < s< оо, интегри-
руема в квадрате и почти всюду отлична от 0,
а х{, ..., хм — ортонормированный базис в гильбер-
товом пространстве RM), образуют максимальную
систему в L2(RM) с максимальным структурным ти-
пом (см. § 2 гл. I). В частности, соответствующая
кратность М совпадает с размерностью простран-
ства RM.
В силу соотношения (1.8) семейство Ht, —оо <
</<оо, имеет тот же структурный тип, что и се-
мейство Ut = Lt (Rm), — оо < / < оо, тогда и только
тогда, когда
0
(1.9)
(что равносильно условию Ff(to) = O в точке /0 = — оо
для всех j = 1, . . ., М).
Обратимся теперь к самому стационарному про-
цессу £(/), — оо < t < оо. Следуя ранее предложен-
ному определению регулярности (см. стр. 52), назовем
процесс £(/), —оо</<оо, регулярным, если семей-
ство Ht(Q, — оо < / < оо, изометрично «стандартному»
семейству Ut — Lt (Rm), — оо < t < оо. Выше мы полу-
чили фактически следующий результат.
Теорема. Стационарный процесс % (t), — оо <
</<оо, регулярен тогда и только тогда, когда
выполняется условие
Г\Н№ = Ъ. (1.10)
§ 2]
ПРЕДСТАВЛЕНИЕ ВОЛЬДА
71
Отметим, что условие (1.10) было положено в основу
хорошо известного понятия регулярности, введенного
А. Н. Колмогоровым при рассмотрении одномерных
стационарных процессов (dim/?=l) с дискретным
временем *).
§ 2. Представление Вольда и факторизация
спектральной плотности
1. Группы унитарных операторов, связанных со
стационарным процессом. Стационарный процесс £(/),
— оо</<оо, можно представлять себе следующим
образом: в гильбертовом пространстве Н (£) задана
группа унитарных операторов Et, — оо < t < оо, и
некоторое подпространство Я0, «сдвиги» которого
Н* = порождают все Н (g); совокупность стацио-
нарно связанных процессов, определенных соотноше-
ниями {£(/), х} = £zx, — оо < t < оо, где х е Я0, будет
стационарным процессом рассмотренного ранее вида.
С другой стороны, если заданы стационарно связан-
ные компоненты {£(/), х}, —оо < t < оо (где параметр
х е /?, а гильбертово пространство R может быть
произвольным), то соотношение
£/{£($), x} = {g(s + /), х}, — oo <s<oo; хе/?, (2.1)
определяет в H(g) унитарные операторы Et, — оо <
</< оо, которые образуют (непрерывную) группу:
EsEt = Es+(. (2.2)
Взяв за Я0 замыкание подпространства величин
,{g(0), х), хе/?, мы получим пару (Bz, Я0) описанного
выше типа. (Ранее, в § 1 мы фактически использовали
функциональную модель стационарного процесса как
пары (/?„ Я0), рассматривая подпространство HQ = f^2R
в функциональном пространстве Л2(/?), и группу уни-
тарных операторов Et — операторов умножения на еш.)
Предположим, что стационарный процесс £(/)>
— оо < t < оо, удовлетворяет условию (1.10) (или, что
*) А. Н Колмогоров, Стационарные последовательности
в гильбертовом пространстве, Бюлл. МГУ 2, № 6 (1941), 1—40.
72
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ [ГЛ. III
то же, его спектральная плотность удовлетворяет
условию (1.9)). Как мы знаем, в этом случае семей-
ство подпространств Ht(%), — оо < /< оо, изометрично
семейству Ut = Lt (Rm), — оо < t < оо, где Rm — гиль-
бертово пространство размерности М (М — кратность
обновляющего процесса для £(/), —оо < t < оо,
a Lt (Rm)—подпространства в функциональном про-
странстве L2(RM), состоящие из функций u(s), — оо <
< s < оо, равных 0 при s > t).
Покажем, что существует обновляющий процесс
с некоррелированными приращениями Х](Ь)—Х;(а),
— оо < а, b < оо; /=1,..., обладающий сле-
дующим свойством:
Et [X, (b) - X, (а)] = Х} (b + t)- Xi (a +1), (2.3)
—.00
t
oo,
где Et— унитарные операторы в пространстве Н (g),
связанные со стационарным процессом £(/) соотноше-
ниями (2.1). Обновляющий процесс, обладающий ука-
занным свойством, будем называть стационарно
связанным с процессом £(/), — оо < t < оо.
Сначала покажем, что группа унитарных операто-
ров Et, — оо < / < оо, изометрична группе унитар-
ных операторов в функциональном пространстве
U = L2(RM), имеющих вид VtTt, где Г,—оператор
сдвига:
Ttu (s) = u(s — t), — оо < 5 < оо,
a Vti —оо </ < оо, —некоторая группа унитарных
операторов в гильбертовом пространстве RM, дейст-
вующая в L2(RM) по правилу:
(Vtu) (s) = Vtu (s), — oo < s < oo.
Действительно, если обозначить Qt проекторы в про-
странстве U на подпространства Ut = L2t(RM),—оо <
< t < оо, то будем иметь
Qt +$—TsQtTs , Pt+s — EsPtEs ,
— oo < s, t < OO
(2.4)
ПРЕДСТАВЛЕНИЕ ВОЛЬДА
73
§ 2]
(ср. с (1.5), (1.6)). Учитывая, что
Pt = XQtX~\ — оо < t < ОО,
где X — некоторый унитарный оператор из U —
= L2(RM) в из второго соотношения (2.4) по-
лучаем
Q<+S = (X-'£SX) Qt(X~'EsX)~l,
что равносильно условию
(X~lEsX) Ut = Ut+s.
Очевидно, унитарный оператор (X~'ESX) T-s = Vs об-
ладает тем свойством, что
VsUt = Ut, -оо</<оо, (2.5)
а как мы знаем, этим свойством обладают лишь
«операторы умножения» на унитарный множитель
(в гильбертовом пространстве RM); см. формулу (2.21)
на стр. 55. В итоге получаем, что
X~'ESX=VSTS, — oo<t<eo.
Очевидно, для любых s', s" операторы Vsf и Ts„
перестановочны.
Возьмем теперь ортонормированный базис xlf ...
.. ., хм в гильбертовом пространстве RM и рассмот-
рим функции интервала \ — [а,Ь) со значениями
в функциональном пространстве U = L2(RM), опреде-
ленные следующим образом:
(х/ —%“)(/) = VfX/• %д (/), — оо</<оо, (2.6)
где хл(0 обозначает индикатор интервала Д = [а, Ь):
%д(/)=1 при /еД, хд(/) = 0 при /<£Д. Очевидно,
w? - «я </-»)]=
= vtxr Х4+,(<) = ^+‘--'“+'-
Очевидно также, что «приращения» хЬ—х“1, хь.г—xaf
ортогональны в L2(RM), если интервалы Д1 = [а1, &2),
Д2 = [о!2, 6) не пересекаются, а при различных j. они
74
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
ортогональны для любых А] и Д2, поскольку при
различных j ортогональны в RM элементы Vtxr
Кроме того, поскольку при любом s элементы Vsx^
/ = 1, А4, образуют ортонормированную систему
в гильбертовом пространстве Z?(/?M), функции
— оо < s < оо, Д = [a, b), ~оо < a; b < t; ]=== 1М,
порождают все подпространство U t = [^(R^. Поло-
жим
%.(&) — X. (а) = Х(хъ. — х“\ /=1, М, (2.7)
где, напомним, X есть унитарный оператор из U =
=L2(RM) в Н (g), удовлетворяющий условию XUt=Hz(£),
— оо < t < оо. Ясно, что определенный выше процесс
с некоррелированными приращениями {Х^Ь)~ХДа)}^,
— оо < а\ Ь < оо, является обновляющим процессом
для g (/), — оо < t < оо. При этом
Et [J, (b) - (a)j = EtX (х> - х*) =
= W/W - - */) ==
= Х/(& + 0-^(а + /)>
так как по определению оператора Vt
X~iElX = VtTt,
а кроме того,
||Х/(гО-*Да)||2 = &-а, /=1,...,Л1. (2.8)
Таким образом, построенный нами обновляю-
щий процесс со «стохастическим дифференциалом»
{dXj(t)}M является стационарно связанным со ста-
ционарным процессом £(/), —оо < t < оо *).
*) Другой метод построения обновляющего процесса был
предложен Ханнером в одномерном случае, dim 7? = 1 (О. Han-
ner, Deterministic and non-deterministic stationary random
processes, Arkiv for Matematik 1, 44 (1950), 161 — 177. См. также
G. К a 1 1 i a n p u г, V. Mandrekar, Multiplicity and represen-
tation theory of purely non-deterministic stochastic processes,
Теория вероят. и ее примен. X, 4 (1965), 614 — 644, где метод
Ханнера применяется к многомерным процессам).
§ 2]
ПРЕДСТАВЛЕНИЕ ВОЛЬДА
75
Теорема. Имеет место следующее канониче-
ское представление:
t м
U(0,4 = j^CjO-OdX^s), -со</<оо. (2.9)
— оо /=1
(Обычно его называют представлением Вольда.)
Действительно, величину {g (0), х} е Яо (g) можно
представить в виде
о м
й(0),4 = ^c/(s)dXl (с),
— оо /=1
где о м
J I Cj (s) |2 ds < OO (2.10)
— oo /=1
(см. формулу (1.8) гл. I), и в случае обновляющего
процесса, который стационарно связан с процессом
|(0, имеем
{g(/),4 = Ez{g(0),4 =
ом t м
= J 2 <4 (s) dxs (s + t)--= / £ с, (s -1) dxs (s).
— oo f==l — oo /==1
2. Факторизация спектральной плотности. Как мы
знаем, стационарный процесс со спектральной плот-
ностью /(Л) регулярен тогда и только тогда, когда
выполнено условие (1.9), а именно,
Л [sV e‘^K2R] = 0.
Докажем следующее предложение *).
*) По-видимому, впервые это предложение появилось для
случая дискретного времени в статье Ю. А. Розанова,
Спектральная теория многомерных стационарных случайных про-
цессов с дискретным временем, Успехи матем. наук ХП1, 2 (80),
(1958). Одновременно для .невырожденной спектральной
плотности проблема факторизации рассматривалась в работах
Винера и Мазани, а также Хелсона и Лоуденслегера (N. W i е-
п е г, Р. Masani, The prediction theory of multivariate stoc-
hastic processes. II, Acta Math. 99 (1958); H. H e 1 s о n, D. Low-
den s 1 ager, Prediction theory and Fourier series in several
variables, Acta Math. 99 (1958)).
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
Теорема. Для регулярности стационарного про-
цесса необходимо и достаточно, чтобы спектральная
плотность допускала факторизацию
— оо<Л<оо,
(2.Н)
где срл при почти каждом К есть линейный оператор
из гильбертова пространства Д в М-мерное гильбертово
пространство Дм, а как функция от К, — оо < А < оо,
Флх е L2 (Д) для всех х е Д, причем
оо
J е~ш<рхх dK = 0 при / > 0.
—оо
(2.12)
(Ср. с факторизацией корреляционного оператора В, one-
ратора умножения на fK, рассмотренной в § 2 гл. II.)
Доказательство. Рассмотрим определенную»
почти всюду функцию
м
их ($) = 2 С/ (^) xh —оо < s < оо (их (s) = 0 при $ > 0),
/=1
где Xi, ..., хм — ортонормированный базис в Дм,
а коэффициенты с{ (s), ..., см (s) те же, что и в пред-
ставлении Вольда (2.9) для компоненты {£(/), х}
(ОО м М
напомним, что J I cj (s) I2 ds < оо, V| су (s) I2 < оо
-оо /==1 /=1
при почти всех s и функция их (s)eZ.2 (А?л1)^. При усло-
вии (2.8) из формулы (2.9) следует, что
оо
EOi)> 4-Ш у}= J {«%($ —иу(s — 4)}ds.
—oo
Положим
oo
фх (Л) = J eiKsux(s)ds, —oo < Л < oo.
—oo
Легко видеть, что соответствие элементов (/), х
§ 2] ПРЕДСТАВЛЕНИЕ ВОЛЬДА 77
из НЦ) и ешфх(Л), — оо < Л < оо, из L2 (Rm) является
изометрическим:
Е{Ш, М У} =
= j eiK{ti~t2) {fKx, y}dk = J {фх (Л), qy (Л)} dk
— OO —00
при всех tlt t2 и x, у e R. Поэтому
{fax, i/} = (<Px(X), <рг/(Л)}
при почти всех Л для любых х, у R. Соответствие
х -> фх (X) е L2 (R^ (такого же типа, как x->f\[2x^
е L2 (/?)) является линейным, иначе говоря, ф есть
линейный оператор из гильбертова пространства R
в функциональное пространство А2(/?м). При этом
мы можем определить ф таким образом, чтобы на
некотором всюду плотном в R подпространстве R
каждое значение фх(Л) было бы линейным оператором
фх (Л) = ф^х, х е R,
при почти всех Л (например, можно взять надлежащие
п
ФХ*(Л), 6 = 1, 2, и положить Фа* == 2 (Ч
П
при х== S ckxk, где х}, х2,... —базис в /?). Из ра-
k=^\
венства
{/\*, У} = {<IV, У R,
вытекает, что при почти всех Л оператор фЛ в R
является ограниченным и может быть однозначно
продолжен с подпространства R на все пространство R
(напомним, что самосопряженный положительный опе-
ратор fK, определенный на всем пространстве R,
является ограниченным при почти всех Л). В итоге
получаем
78
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. Ill
где — оператор из RM в 7?, сопряженный к фл;
Фь удовлетворяет условию (2.12), поскольку
о
= "Йл- J* e‘KSUx dS’ Ux S L2
В свою очередь, факторизация спектральной плот-
ности fK влечет за собой справедливость условия
регулярности (1.9)—(1.10). Действительно, взяв в функ-
циональном пространстве U = L2(RM) подпространство
{7°^£о(^м)> порожденное функциями
оо
= f e~iKt^xdk, xeeR
у 2л J
—оо
(их (/) = 0 при t > 0),
будем иметь
оо
Е{Ш- = / e^~*{fKx, y}dk =
— оо
оо оо
= / el^b~a} {qKx, <px£/}dZ,= J {Taiix(t), Tbuy(t)}dt,
—оо — оо
где Ts — операторы сдвига в L2(/?M). Видно, что ото-
бражение {g(a), х} <-> Таих есть изометрия из Н (g)
в L2(RM), и поэтому семейство ЯД£), — оо < t < оо,
изометрично семейству подпространств вида
Ht=\J TSU\ — оо </< оо.
Поскольку Ht ^L2 (Rm) и Q Ht £ Q L2 (RM) = 0, имеем
t
В заключение этого пункта введем аналити-
ческие классы Я6, 6=1, 2, операторных функций Гх,
— оо < Z < оо, значения которых суть линейные опе-
раторы в гильбертовом пространстве (I\: R'R").
Будем говорить, что функция 1\, —оо < Z < оо, при-
§ 3]
КРАТНОСТЬ
79
надлежит классу Я6, если для всех х е /?', у е /?"
скалярное произведение {1\х, у} как функция от А,
принадлежит пространству L6 (на прямой) и
J e~~iKt {1\х, у} dk = 0 при t < 0 *). (2.13)
—оо
Например, если операторная функция <рА удовле-
творяет условию (2.12), то сопряженная функция
фА = ФА принадлежит классу Н2, поскольку для любых
№ R', у е R" (где R' — RM, R" = R)
{ФлХ, у} = {х, qKy}(=L2 (на прямой)
и
| е-ш{фАх, y}dk = J etM{<pKy, x}dK — 0 при t < 0.
—оо —оо
Отметим, кстати, что если
/л = "Фл ' "Ф!» — оо<Л.<оо,
где операторная функция е Я2, то функция
будет удовлетворять условию (2.12).
§ 3. Кратность регулярного стационарного процесса
Рассмотрим регулярный стационарный процесс
g(Z), —оо </ < оо, с компонентами {£(/), х}, —оо<
</<оо, где хе/? (/?— «параметрическое» гиль-
бертово пространство). Пусть fK — спектральная плот-
ность. Как мы знаем, тип семейства ЯД£), — оо <
</<оо, и, в частности, кратность М обновляющего
*) По поводу определения ‘ классов Я6 и их свойств
см., например, книгу К. Г о ф м а н а, Банаховы пространства ана-
литических функций, М., ИЛ, 1963, и обзор В. И. Крылова,
О функциях, регулярных в полуплоскости, Матем. сб. 4, 46
(1938), 9—30. См. также: И. И. Привалов, Граничные свой-
ства аналитических функций, М.—Л., Гостехиздат, 1950; Б. Секе-
фаль ви-Надь и Ч. Ф о я ш, Гармонический анализ опера-
торов в гильбертовом пространстве» М»» изд-во «Мир», 1970,
80
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
процесса для £(/), —оо</<оо, целиком опреде-
ляются спектральной плотностью /х(см. § 1). Возни-
кает вопрос, как по спектральной плотности fK найти
кратность М обновляющего процесса.
Из условия факторизации fx = qp*.(pA, где qpx —
операторная функция из гильбертова пространства R
в гильбертово пространство RM (см. формулу (2.11),
вытекает, что
(3.1)
Действительно, положив фл=(р*, получим
М = (3.2)
поскольку fKx — Q тогда и только тогда, когда
{/\х, х} = {<Рлх> Фх-Н = 0, и если — подпространство
всех элементов х е R, для которых fKx = <ркх = 0, то
Заметим, кстати, что из условия (3.2) вытекает
следующее свойство спектральной плотности fK:
dim fKR = const. (3.3)
Действительно, если взять любые конечномерные
подпространства /?' RM, R" R (выбрав в них
ортонормированные базисы х', ..., х'т и х", ..., х")
и рассмотреть проекции элементов фхх, xg R', на R",
то окажется, что соответствующее подпространство
п
элементов вида S х'£\ х'£, х' е /?' (обозначим
k—\
его qpx^7^,z) будет иметь постоянную почти всюду
размерность, поскольку она совпадает с рангом ана-
литической матрицы с компонентами (фдХ<, х^},
/ = 1, ..., m; k=l, ..., п, каждый из миноров
которой как функция от Л, —оо < Л < оо, принад-
лежит одному из хорошо известных аналитических
классов Н6 при надлежащем 6 > 0 и либо почти
всюду отличен от 0, либо равен 0 тождественно *)♦
*) См. сноску на стр. 79,
КРАТНОСТЬ
81
§ 3]
Теперь уже ясно, что размерность
dim^7?z==max(dim tyKR'/R")
R”^R
постоянна почти всюду для любого конечномерного
подпространства R' RM, а следовательно, постоянна
и размерность
dim dim = dim fKR.
Покажем, что
М < dim fKR. (3.4)
Наше определение кратности М случайного процесса
£(/), —оо < / < оо j связано с семейством проекто-
ров Pt на подпространства ЯД£), —оо</<оо.
Будем исходить из того, что М равно минималь-
ному числу циклических подпространств, замкнутая
линейная оболочка которых совпадает со всем про-
странством Н (£) — см. § 2 гл. I — (напомним, что
«циклическим» по отношению к проекторам Pif
— оо < t < оо, мы называем подпространство в Н (£),
порождаемое величинами вида Ptx, — оо < / < оо,
где х— некоторый элемент из Я(£)).
Рассмотрим циклические подпространства ЯД£),
/ = 1, ..., Л4, порождаемые соответствующими при-
ращениями Xj (/) — Xj (s), — оо < s, t < oo, обновляю-
щего процесса в представлении Вольда (2.9) (поро-
ждающим элементом является, например,
оо
&/ = / С ($) dXf(S),
— ОО
где подынтегральная функция интегрируема в квад-
рате и отлична от 0 почти всюду). Выберем такую
функцию c(s), — оо < s < оо, чтобы ее сдвиги
c(s — t), —oo<s<oo, где параметр t меняется
в пределах —оо</<оо, порождали все простран-
ство L2 (на прямой). Ясно, что тогда величины
оо
= / c(s — f)dXf(s), —оо<К<х>,
порождают все подпространство ЯД£).
82
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
Обратимся теперь к группе унитарных операто-
ров Et, —оо</<оо, связанных со стационарным
процессом g(/), —оо</<оо, соотношением (2.1), и
используем известное представление Стоуна *):
Et = \e^dEK, (3.5)
— оо
где Ек, — оо < Z < оо, —семейство проекционных
операторов в пространстве Н (g). Очевидно,
МО ==£/£/> — оо</<оо,
откуда следует, что ортогональные подпространства
(g), /= 1, ..., М, являются циклическими также
и по отношению к семейству проекторов Ек,
—оо < Л < оо (Нj (g) порождаются величинами £xgz,
—оо < А < оо). При этом все «циклические векторы» gz
имеют один и тот же «лебеговский тип»: функции
Fj (X) = gy) имеют положительную почти всюду
плотность
оо 2
J eiKs с (s) ds .
—оо
-- ОО < Л < ОО,
и, таким образом, семейство проекторов — оо <
<Х<оо, в пространстве Н (g) имеет М-кратный
«лебеговский тип». Поэтому, М. равно минимальному
числу циклических (по отношению к семейству
—оо < Л < оо) подпространств, замкнутая линейная
оболочка которых совпадает со всем пространством
/7(g). Если же обратиться к функциональной модели
(1.2) рассматриваемого стационарного процесса и
соответствующим операторам Et, — оо < / < оо, умно-
жения на функцию еш, —оо < Л < оо, в гильберто-
вом пространстве L2(R), то легко указать N цикли-
ческих подпространств (N = dim fKR) вида
нк = V еМ№хк,
— oo<Zt < оо
*) См., например: Ф. Рисе, Б. Секефальв и-Н а д ь, Лек-
ции по функциональному анализу, М., ИЛ. 1954.
§ 4]
УСЛОВИЯ РЕГУЛЯРНОСТИ
83
замкнутая линейная оболочка которых дает все
пространство Н= V etKtf}!2R; именно, можно
— оо < t < оо
взять функции f^Xk, k = 1, ..N, значения которых
образуют базис в подпространстве f[/2R s 7?
(см. далее лемму на стр. 87). Отсюда следует, что
dim f^R(= dim fKR).
Соединяя неравенства (3.1) и (3.4), получаем сле-
дующий результат.
Теорема. При условии регулярности крат-
ность М обновляющего процесса равна размерности
подпространства fKR е R:
М = dim f\!2R п. в. (3.6)
§ 4. Условия регулярности
1. Общий критерий регулярности. Нашим исход-
ным пунктом будет теорема о факторизации, со-
гласно которой стационарный процесс g(Z), — оо <
< f < оо, со спектральной плотностью —оо<Л<оо
(в гильбертовом пространстве /?) является регуляр-
ным тогда и только тогда, когда существует опера-
торная функция <рх, —оо < К < со, класса Н2
такая, что
^ = Фл-Фх п-В- <4Л)
(определение класса Н2 и аналогичного класса Я1,
который нам понадобится в дальнейшем, дано на
стр. 78). По сравнению с ранее предложенным усло-
вием (2.11) мы поменяли местами фА и ф*, так что
теперь qp^ — линейный ограниченный оператор из гиль-
бертова пространства R в Л4-мерное гильбертово
пространство RM. Поскольку всегда dim/?^M
(Af = dim/\/? п. в.), то можно считать, что RM S /?,
и доопределить оператор срл = (ф^)* на всем простран-
стве R, положив Ф^(7?@ RM) = 0.
Введем оператор Ух, положив
уЖ/2*) = ф1х> х6Е/?-
84
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
(ГЛ. III
Из равенства (4.1) видно, что
и, таким образом, изометрично отображает под-
пространство f\>2R на подпространство Следо-
вательно, сопряженный к нему оператор
(<₽1Х)= Я/2х’ х *= также является изометриче-
ским; при этом
Поскольку то для сопряженного опера-
тора = (<₽’)’ получаем, что <рку ±= f%2VKy при y<=q%R,
и если___доопределить оператор VK так, чтобы
= 0, то будем иметь
Фх = П/2^
поскольку ортогональное дополнение RQVkR к под-
пространству ф^Т? состоит из «нулей» оператора ф^,
сопряженного к ф^. Очевидно,
f^KR=i\R = ffR
(здесь и в дальнейшем /\1/2 обозначает обратный
оператор к сужению на подпространстве
Ясно, что функция
(4-2)
удовлетворяет условиям
а кроме того,
/ | ||2 dK < оо, xt=R,
§ 4]
условия . регулярности
85
поскольку
= 1v^x' 11 v>x ч 11 х и-
При этом функция ф^, — оо < Л < оо, принадлежит
классу Н1 П Я2, потому что для любых х, у е R ска-
лярное произведение {«р^х, у} как функция класса Я2
есть преобразование Фурье некоторой функции c(t)
из L2 (на прямой — оо < t < оо), с (/) = 0 при t < 0, а
1
1 — г'Л
J etMe~fdtt
о
и следовательно,
оо
{Фх*, у} = д = J еМс * e~tdt *= Н' Н2,
О
где с*е~*> t > 0, обозначает свертку указанных
функций.
Докажем теперь следующее предложение *).
Теорема. Стационарный процесс со спектраль-
ной плотностью fK является регулярным тогда и только
тогда, когда существует операторная функция
класса Н] такая, что
f;l,\KR = f^R п. в. (4.3)
и при любом xeR функция Д12ф^х, — оо < А < оо,
принадлежит пространству L2{R):
00
У ||frI/4xl2rfx< °° •
*) При доказательстве мы следуем методу, который был
фактически предложен в одной нашей работе о стационарных
процессах с дискретным временем (Ю. А. Розанов, О ли-
нейном интерполировании стационарных процессов с дискретным
временем, ДАН СССР 116 (1957), 923—926); см. также Yu. A. Ro-
zanov. Some Approximation Problems in the Theory of Sta-
tionary Processes, J. of Multivariate Analysis 2, 2 (1972)„
135-144).
86
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. Ш
(При этом необходимым является существование опе-
раторной функции удовлетворяющей условиям
(4.3)—(4.4), которая принадлежит не только классу Я1,
но одновременно также и классу Я2.)
Мы уже убедились выше, что если стационарный
процесс регулярен и, следовательно, его спектральная
плотность fK допускает факторизацию Д —
с помощью операторной функции срх класса Я2, то
определенная формулой (4.2) функция г|^ принадлежит
обоим классам Я1, Я2 и удовлетворяет всем условиям
нашей теоремы. Таким образом, эти условия являются
необходимыми. Докажем, что они являются доста-
точными.
Напомним, что регулярность имеет место тогда и
только тогда, когда выполняется условие (1.9):
»--n[V «Тем.
t Ls < t J
Рассмотрим в пространстве
Н = v еа7'/2/? <= L2 (/?)
—оо < S < оо
ортогональное дополнение
\ = HQH0
к подпространству Яо = V Очевидно, все
подпространства eiKt\ (состоящие из функций вида
ешх(А), — оо < Л < оо, где х (Л), — оо < К < оо, е Д)
ортогональны определенному выше подпростран-
ству Я-оо, поскольку eiKt& ортогональны соответствую-
щим подпространствам
Я,= V easf^f> = eiMHQ'
Ясно, что если
V ешд = я,
—оо < t < оо
(4.8)
то Я-оо — 0.
УСЛОВИЯ РЕГУЛЯРНОСТИ
87
§ 4]
Определим пространство-функцию Л (Л), —оо <
< Л < оо, для любого подпространства ЛьА2(/?),
выбрав в А полную систему функций {^1(Л),а2(Л),...} = S
и положив А (Л) = Л5 (Л), где (Л) есть замкнутая ли-
нейная оболочка в гильбертовом пространстве У? всех
значений а{ (Л), а2(Л), ... Поскольку для любой функ-
ции a(Z)e4 найдется последовательность линейных
комбинаций вида ^ckak(K), сходящаяся в Л2(7?)ка(Л),
k
то некоторая подпоследовательность S ckak (Л) схо-
k
дится к а (К) п. в., так что a(A)e/ls(Z) при почти
всех Л, й поэтому для любой другой полной системы
Sz= [а! (Л), 6/2 (^), ’ • •} в A^L2(R) имеем
Д5,(Л) = Л5(Л) п. в.
В этом смысле рассматриваемая п. в. пространство-
функция__А (Л), — оо < Л < оо, определяется равен-
ством Л(Л) = Л5(Л) однозначно.
Лемма. Подпространство
L = V ~^А L2 (/?) (4.6)
—оо < t < оо
состоит из всех функций х(Л), —оо < Л < оо, е Z,2(/?),
значения которых удовлетворяют условию
х(Л)еЛД) п. в. (4.7)
Доказательство. По определению подпро-
странства L^L2{R) всякая функция .r(Zjs7 есть
предел функций вида ^el tkak(K), где а^(л)еЛ, и,
k
очевидно, значения х(Л) удовлетворяют условию (4.7).
Далее заметим, что для любой функции а (Л) е А
/ оо
и скалярной измеримой функции с(Х) I J | с (Л) I2 X
' —оо
X || а (Л) II2 < °°) произведение х (Л) = с (Л) а (Л) есть
функция из подпространства L, поскольку функция
с (Л) может быть сколь угодно точно аппроксимиро-
вана в среднем (с весом g (Л) —1| а (Л.) |р) линейными
88
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
комбинациями вида S cketUk и для некоторой после-
k
довательности вида ^etKtkcka(X) в L2(R) имеем
k
оо
/ |х(А) — =
— ОО k II
2
с (Л) — ^ckeiM* 2 g (X) <Д->0.
k
Аналогично, для любых функций (Л), ..., аг (Л) е А
и векторной измеримой функции с(А)= (С[(Л), , сг(А)},
ОО
удовлетворяющей условию J с (Л) g (Л) с (Л)* dK < оо
(где g(K)— положительная матричная функция с ком-
понентами gpq (Л) = {ар(Л), aq(h)}, р, q = 1, ..г),
г
функция х (Л) = 2 СР (А) аР (А/) принадлежит подпро-
р=1
странству А, поскольку векторная функция с (Л) может
быть сколь угодно точно аппроксимирована в среднем
(с весом g(A)) линейными комбинациями вида
^etKtkck (с векторными коэффициентами ck =
k
— [ск1, ..., ckr}), и для некоторой последовательности
функций 2 e'Uk S ck[>aD(k) в L2(R) имеем
k р—1
2
— oo || k P=1
oo
Р=1
k
_ eiMkck
k
Возьмем теперь произвольную функцию х(Л)е£2(/?),
удовлетворяющую условию (4.7), и полную систему
функций aJA.), а2(Х), ... s А, значения которых поро-
ждают подпространство А (Л) s Я при п. в. К. Проек-
ция (в гильбертовом пространстве. R) величины.
УСЛОВИЯ РЕГУЛЯРНОСТИ
89
§ 4]
х(А)еЛ(А) на подпространство, порожденное конеч-
ным числом величин (Л), ..., аг(Л), имеет вид
г
хг (Л.) = 2 с₽(^) apW> гДе коэффициенты (Л), ..., сг(К)
p=i
как функции от Л являются измеримыми в силу того,
что для любых функций х(Л), у (Л) е Л2 (7?) скалярное
произведение (х(Л), у(Л)) является измеримой функ-
цией от Л. Мы показали выше, что .ггР.)е£. Но при
каждом фиксированном Л
|| х (Л) — xr (X) ||2-» 0 при г—>оо,
причем || х (X)—хг(К) ||2 ^|| х(Х) ||2, и следовательно,
оо
j || х (Л) — xr (X) II2 dk -> 0 при г-> оо,
—оо
так что Лемма доказана.
Согласно этой лемме подпространство Н —
= V eiKtfb2R^ L2(R) состоит из всех тех функ-
ций х(л)еА2(/?), значения которых удовлетворяют
условию
x(V^2R п. в., (4.8)
и условие регулярности (4.5) можно выразить сле-
дующим образом:_____ ___
________ А(Л) = ^п.в., (4.9)
где А (Л), — оо < Л < оо, обозначает пространство-
функцию, порождаемую подпространством А =
= HQH^L^.
Предположим теперь, что существует некоторая
операторная функция % класса Я1:
сю
J e~iU {фдх, у} dk = 0 при t < О
—оо
для всех х, у s R, удовлетворяющая требованиям
(4.3), (4.4). Очевидно, функция х(Л) = f^ll2tyKx е L2(R)
удовлетворяет условию (4.8) и, следовательно, при-
надлежит подпространству Н= V elMf'fcR. Кроме
go
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. Ш
того,
оо оо
J [х(Х), f'fcy] dk = | e~iU Цлх, Ц dk—Q при t <0,
— оо ’ —СО
т. е. функция х(Л), — оо < А < оо, ортогональна
подпространству HQ— V eMflJ2R-, иначе говоря, функ-
t < о
ция х{К) = f^i/2фдХ, —оо < Л< оо, принадлежит под-
пространству k = HQHQ. Поэтому
П. В.,
и в силу соотношений (4.3) получаем, что
ДД) = п. в.,
откуда следует условие регулярности (4.5).
Теорема полностью доказана.
2. Некоторые выводы и примеры. Остановимся
подробнее на предложенных выше условиях регуляр-
ности (4.3)—(4.4).
Рассмотрим конечномерный случай, dim/? = JV.
Будем считать, что в унитарном пространстве R
выбран некоторый ортонормированный базис xlf ...
..., xN и спектральная плотность задана соответ-
ствующей матрицей /(Л) = {fkJ (Л)} с компонентами
fkj<X) = {h.Xk, х^, k, / = 1, N.
В одномерном случае, dim/? = l, условия (4.3) и
(4.4), очевидно, означают, что скалярная спек-
тральная плотность /(Л) должна быть отлична от 0
почти всюду и для некоторой скалярной функции ф(Л)
класса Я1:
J J -Ф(Л) Г/(Л)"1 tZX < оо.
—оо
Отсюда следует, что
J ------Г+v--------
—oo
log 1w |2 , f 10g^W/7A<'oo
—1~dK ~ J -Г+Л^Й'Л< °0’
§ 41
УСЛОВИЯ РЕГУЛЯРНОСТИ
91.
и поскольку для функции ф(Л) класса Я1
I log I -ф (Л) I21 „ .
--------^<00,
мы приходим к хорошо известному условию Колмо-
горова— Крейна *):
logfW
1 +V
dK
— оо.
(4.Ю)
В случае произвольного конечномерного процесса
с обратимой п. в. спектральной плотностью f(K)
условия (4.3)—(4.4) означают, очевидно, что суще-
ствует невырожденная NXN матричная функ-
ция Ф(Л) = {ф^у(Л)} с компонентами из класса Я1,
для которой
оо
J 8р[ф(Л)Г‘(Л)ф(Л)‘]</Л<оо.
—оо
Отсюда следует, что
log det [i|> (Л) f (Л) ip (Л)*]
det i|> (Л) |2 _
1 + V
? det f (Л)
1 4-V
и поскольку для матричной функции ф(Л) класса Я1
оо
Г I log I det ф (Л) |21 ,,
J ----Г+V ^<оо>
*) В цитируемой ранее работе Колмогорова (см. сноску на
стр, 71) получено условие регулярности для одномерных ста-
ционарных процессов с дискретным временем; на случай непре-
рывного времени оно было обобщено Крейном (М. Г. Крейн,
Об одной экстраполяционной проблеме А. Н. Колмогорова,
ДАН СССР 46 (1944), 306-309).
92
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
мы приходим к хорошо известному условию регу-
лярности *)
log detf(X)
1 + Л2
dK > — оо.
(4.11)
Отметим, что это условие равносильно следующему:
Iogll/л ‘II
1 + А2
> — оо,
поскольку
log det f (А) = 2 log mk (A), log || fK 11 ‘ = log m, (A),
k—\
где m, (A) ... mN (A) —собственные значения
спектральной плотности fK.
Для регулярного стационарного процесса ранг
спектральной плотности f (А) п. в. равен постоян-
ному М (см. теорему § 3). Следовательно, хотя бы
один из ее главных миноров порядка М отличен от О
на множестве положительной меры. Пусть это будет
определитель М X Л/-матричной функции ?м(Ь) =
= {^р^(Л,)} с элементами fkpk(l (А) = {/Лр, xkq], р,
<7=1, М. Обозначим RM подпространство в R,
порожденное элементами xkp, р=1,..., М. Оче-
видно, является (матричной) спектральной плот-
ностью 7И-мерного регулярного процесса с ком-
понентами (I (0, х}, х е RM, и, следовательно,
det fM (Л) У= 0 не только на множестве положительной
меры, но и почти всюду. Как мы только что отме-
чали, в случае обратимой спектральной плотности fM (Л)
*) Для случая дискретного времени это условие было пред-
ложено в заметке В. Н. 3 а с у х и н а, К теории многомерных
стационарных процессов, ДАН СССР 33 (1941), 435—437. Слу-
чай непрерывного времени рассмотрен в работе Е. Г. Глад ы-
ш е в а, О многомерных стационарных случайных процессах,
Теория вероят. и ее примен. III, 4 (1958), 458—462.
УСЛОВИЯ регулярности
93
ее наименьшее собственное значение tn (К) удовлет-
воряет условию
При этом, очевидно,
м м
inf 5 lifkpk.Wcp <
М q=\ р=1 р 4
?ыг=*
< inf 2 2 fjkp W ср
м /=»1 р=1 Р
S I с„ |2= 1
1 I р I
inf ||fxx||.
x^Rm, 11*11=1
Рассмотрим теперь подпространство fKR = RQR^,
где R% — нулевое подпространство для положитель-
ного оператора fK. Поскольку dim/\/?M = dim А 7?,
имеем f\RM = f\R п. в., откуда следует, что проек-
ция хк любого элемента х е RM, х =£ О на подпро-
странство fKR отлична от 0 п. в. Кроме того, f^x —
= fKxK для № R, и подпространство всех элемен-
тов хк, xs R, совпадает с fKR. Очевидно,
inf II A* IK inf II А* 11 =
где, напомним, fK 1 есть обратный оператор к поло-
жительному оператору fK в унитарном простран-
стве fKR. Таким образом, т (Л) С | f\l |] п. в., и
следовательно, в регулярном случае
(4.12)
Как мы знаем, кроме этого условия для регуляр-
ности необходимо, чтобы существовала операторная
94
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. III
функция qpx класса Я2, обладающая следующим
свойством:
П. В.
Очевидно, в конечномерном случае, dim R < оо>
это равносильно условию
fi.R = <PkR п- в. (4.13)
Из предложенного нами общего критерия регу-
лярности легко вытекает, что
совокупность условий (4.12) — (4.13) не только не-
обходима, но и достаточна для регулярности конеч-
номерного стационарного процесса *)•
Действительно, как известно, при условии (4.12)
интегрируемая скалярная функция
^||Д||) допускает факторизацию Ц^Ц 1 = | 0 (Z) |2 с
помощью некоторой скалярной функции класса Я2
(см. сноску на стр. 79), и если взять функцию
(Л) = 0 (A,) qp (Л), принадлежащую классу Я1, то ока-
жется, что она удовлетворяет всем условиям (4.3) —
— (4.4), поскольку
II fr1/2%x К РГ1/2II • 10 (ЧI • || qv II < II4V ||> *65 Я.
Стоит отметить, что условие (4.12) равносильно
следующему:
log det fM (Л)
1 + V
dK > — оо
(4.129
для некоторого главного минора det fM (%) порядка М
матричной спектральной плотности (ранга М).
*) В случае дискретного времени аналогичные условия
были предложены Хелсоном и Лоуденслегером (Н. Н е 1 s о п,
D. L о u d е n-slager, Vector-valued processes, Proc. Fourth Ber-
keley Symposium Math. Statist. Probability, Berkeley, 1961).
Условия несколько другого типа были ранее найдены Матвее-
вым (Р. Ф. М а т в е е в, О регулярности многомерных стацио-
нарных процессов, ДАН СССР 126 (1959)).
§ 41
УСЛОВИЯ РЕГУЛЯРНОСТИ
95
Вместе с условием (4.13) это позволяет, например,
легко получить следующий важный результат: вся-
кий стационарный процесс, имеющий спектральную
плотность f (Л) = {fkj (Л)} с рациональными элемента-
ми fkjWy /—1, ..., является регулярным.
Остановимся теперь на некоторых особенностях
бесконечномерного случая, dim/?==oo.
Вообще говоря, условие (4.12) уже не будет не-
обходимым *). Например, если взять стационарный
процесс со спектральной плотностью вида —
= fk(h)xk, k=l, 2, ..., где х2, ...—ортонорми-
рованный базис в R и скалярные спектральные
плотности fk (Л),
Г Jog fk W
J 1 + X2
dK > — oo,
монотонно убывают таким образом, что
f lim =
J /г->оо
то окажется, что, хотя процесс регулярен, условие
(4.12) будет нарушено.
В то же время условие (4.12) в случае обрати-
мой спектральной плотности fK является достаточным
для регулярности**).
*) Некоторый класс спектральных плотностей (обладаю-
щих «скалярным кратным»), для которых условие типа (4.12)
является необходимым условием регулярности, указан в книге
Б. Секефальви-Надя и Ч. Фояша, Гармонический
анализ операторов в гильбертовом пространстве, 1970,
**) В случае дискретного времени это фактически было до-
казано в работе Ю. А. Розанова, Спектральная теория
многомерных стационарных случайных процессов с дискретным
временем, Успехи матем. наук XIII, 2(80), (1958), 93— 142; см.
также: A. D е v i n a t z, The factorization of operator-valued
functions, Ann. of Math 73 (1961), 458 — 495.
В связи с условием (4.12) следует упомянуть интересный
пример Лакса, в котором (при обобщении его на случай непре-
96
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ. Ill
В самом деле, пусть fK равномерно интегрируема:
оо
J Ilf (Л) ИЛ < оо,
—оо
и при почти всех Л существует ограниченный
обратный оператор f^1 в R (7? = удовлетворя-
ющий условию
dK > — оо.
Пусть
। О W г
где О (Л) — скалярная функция класса Я2, фактори-
зующая интегрируемую функцию Ц^Ц 1 (Ц^Ц
<||/\||), а —единичный оператор в /?. Очевидно,
операторная функция класса Я!ПЯ2 будет отве-
чать требованиям (4.3) — (4.4), и следовательно, ста-
ционарный процесс со спектральной плотностью fK
будет регулярным.
В случае, когда спектральная плотность fK при
почти всех Л является обратимой на замкнутом под-
пространстве fKR (fKRc:R), регулярность будет иметь
место, если кроме условия (4.12) потребовать еще,
чтобы пространство-функция fKR, —оо < Л < оо, бы-
ла «аналитической» в том смысле, что
«• в. (4.14)
рывного времени)
log fk
1 +Л2
dK > — cl,
где с — некоторая постоянная, а I — единичный оператор в
гильбертовом пространстве R, и тем не менее, соответствующий
стационарный процесс не является регулярным (Р. Lax, On
the Regularity of Spectral Densities, Теория вероят, и ее при-
мен. VIII, 3 (1963)).
§ 4] УСЛОВИЯ РЕГУЛЯРНОСТИ 97
(ср.с (4.13)) для некоторой операторной функции <рх
из класса Н2.
Действительно, в качестве функции удовлет-
воряющей условиям регулярности (4.3)—(4.4), можно
взять
= 6 (А.) • <р (А,),
где скалярная функция 6 (X) класса Я2, как и прежде,
факторизует функцию Ц^Ц (напомним, что —
ограниченный обратный оператор к fK на замк-
нутом подпространстве fKR = f^R, так что =
= f'^R, если = fKR).
В заключение выведем из наших условий регуляр-
ности (4.3)—(4.4) одну теорему сравнения.
Для краткости назовем спектральную плотность
регулярной, если регулярным является соответствую-
щий стационарный процесс.
Предположим, что спектральные плотности Д и
gK при почти всех Л имеют в гильбертовом прост-
ранстве R одни и те же «нули» и, более того, для
элементов хп е R
h.xn->0~ gKxn-+0. (4.15)
Если fK регулярна, а
gx>fk п. в., (4.16)
то спектральная плотность gK также регулярна.
Действительно, при условии (4.15)
причем
где С\ — обратимый оператор в подпространстве
R^ (непосредственно из соотношения (4.15) вытекает,
что определенные равенствами CJ^2x = g^2x и
== /д/2х, х е R, операторы Ск и являются
4 Ю. А. Розанов
98
РЕГУЛЯРНЫЕ СТАЦИОНАРНЫЕ ПРОЦЕССЫ
[ГЛ, III
непрерывными в R}). Поэтому g1/2 = f’/2C[ и
^/2=^/2(С1)_1’ откуда получаем, что^Я = =
= fx2cK^K ~ I'kR’ и> аналогично, f^R<= g^R, т. е.
О = п-в-
Для регулярной спектральной плотности fK сущест-
вует операторная функция класса Я!ПЯ2, удо-
влетворяющая условиям (4.3) — (4.4). Учитывая, что
^Г1/2== гДе — обратимый оператор (со-
пряженный к С\), приходим к заключению, что усло-
вие (4.3) выполняется и в отношении спектральной
плотности g(Z):
(“ftM
Кроме того, поскольку gK fK, имеем п. в.,
и следовательно, в отношении g’(A) выполняется
также условие (4Л)
оо со
— оо —оо
при всех хе/?. Таким образом, спектральная плот*
ность g(Z) является регулярной»
В частности, для конечномерного случая отсюда
вытекает, что если невырожденная спектральная плот-
ность Д является регулярной, то регулярной будет
и всякая спектральная плотность gK, gK^fK- Для
бесконечномерного случая без ограничений типа (4.15)
это, вообще говоря, уже неверно. (Именно, для всякой
функции х(^)еЯ0, х(Х)^0, где Я0 есть замыкание
в L2(R) подпространства всех функций g^2x, — оо <
< Z < оо, x^R, в случае регулярности gK одномер-
ный «стационарный процесс» {eiKtt х(А)} также регу-
§ fl
УСЛОВИЯ РЕГУЛЯРНОСТИ
99
лярен и должно быть выполнено условие
f log IIX WII
J 1 4-х2
—ОО
dK > — оо;
в частности, х (Л) =# 0 п. в. Если же взять строго по-
ложительные (постоянные) операторы Р и Q в R та-
кие, что Q Р и Qxn —> х0 #= О, Рхп -> 0 для некото-
рой последовательности хп е R, взять скалярную
функцию 0(Л) класса Я’ПЯ2, |0(Л) ^1, то ока-
жется, что операторная функция fK= 0(Л) |2Р2 регу-
лярна, /Л = ФЛ-<р;, фа = 0(Л)Р, а функция gKi
_( 10 (Л) I2Р2 при |Л| < 1,
10(Л) I2Q2 при |Л | > 1,
хотя и мажорирует fK, не будет регулярной, поскольку
я1/2^„->^(х)={
О при | Л, | I,
| 0 (Л)| х0 при | X | > 1.
Уточним, что этот контрпример существенно беско-
нечномерный. В качестве Р и Q можно взять сле-
дующие операторы: Q — положительный невырожден-
ный оператор, для которого QR R, a P — QnQ,
где л — оператор проектирования на ортогональное
дополнение к элементу х0 е R, х0 ф QR *).
Имея в виду условия регулярности (4.3) — (4.4),
заметим, что для определенных выше плотностей
и gK, IЛ | > 1, g^2fKR = (QnQ)R = nR^^R=R,
и это указывает на невозможность подобрать опе-
раторную функцию <= Н' П Н2, удовлетворяющую
условиям (4.3) одновременно в отношении fK и gK
(как это было сделано выше при доказательстве
нашей теоремы сравнения).
*) См. R. G. D о u g 1 a s , On factoring positive operator
functions, J. Math. Meeh. 16(1966), 119— 126.
4*
Г Л А В A IV
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
§ 1. Понятие эквивалентности.
Вероятностная интерпретация в случае
гауссовских распределений
Пусть £(/) = {£ (О, x}xl-R и t](O = {n(O> 4^ —два
случайных процесса на интервале tQ < t < Т, где па-
раметр х пробегает некоторое множество /?. Имея
в виду вопрос о том, при каких условиях их обно-
вляющие процессы будут одного и того же типа,
рассмотрим отображение
Л: {г|(/), *}->{£(/), Д x^R- tQ<t<T. (I.I)
Будем говорить, что процесс £(/) = {£(/),
эквивалентен процессу ц(/) = {ц(0, x)x<=r на интеР"
вале tQ < t < Г, если это отображение продолжается
до линейного ограниченного обратимого оператора А
из гильбертова пространства Н (т]) в гильбертово
пространство Н (g) и, кроме того, если разность
/—-Л*Л будет оператором Гильберта — Шмидта (на-
помним, что Я(ц) и Н (%)— замкнутые линейные обо-
лочки соответствующих значений {т](0, 4 и {£(0, х}>
хе/?; /0 < t < Г).
Как было предложено ранее (см. § 2 гл. II), при
условии, что А — ограниченный линейный опе-
ратор, вместо исходных процессов g(f) и T](f), /0 <
< t < Tt можно рассмотреть обобщенный процесс
g(u) = Ли, и е U,
(1.2)
§ и
ПОНЯТИЕ ЭКВИВАЛЕНТНОСТИ
101
на гильбертовом пространстве U — Н (rf), в котором
выделено семейство подпространств Ut = Ht (r|), to <
< t < Т. Именно, вопрос об однотипности обновля-
ющих процессов для £ (/) и т) (/) (на интервале t0 <
< t < Т) есть, по существу, вопрос о том, при каких
условиях семейство подпространств
Ht(£) = AHt(vi), t0<t<T,
(1-3)
будет изометрично семейству ЯДц), /0 < t < Т:
= /0</<Л (1.4)
для некоторого изометрического оператора X.
Было доказано (см. теорему § 3 гл. II), что если
корреляционный оператор В = Л*Д является обрати-
мым и, кроме того, разность 1 — В является опера-
тором Гильберта — Шмидта, то рассматриваемые се-
мейства Ht(g) и ЯДц), tQ<t<T, будут изометрич-
ными, а следовательно, обновляющие процессы для
исходных случайных процессов £(/) и ц(/), /0 < t < Т,
будут одного и того же типа. Таким образом, экви-
валентные процессы £(/) и т) (£), /0 < t < Т, имеют
обновляющие процессы одного и того же типа.
Отметим, что поскольку А = ХВХ1\ где X — уни-
тарный оператор, и
/ — (Л"1)* (Д-1) = I — хвх-' =
= Х (/ — В) X~l = X(I- AM) X-',
то в случае, когда разность I — А*А является опера-
тором Гильберта — Шмидта, оператором такого же
типа будет и разность / — (Л-1)*Л-1. Таким обра-
зом, определение эквивалентности |(£) ~ "П (£) явля-
ется симметричным. Легко видеть, что оно также
транзитивно и, конечно, рефлексивно.
Введенное выше условие эквивалентности хорошо
известно для гауссовских случайных процессов
%(t) и т|(0> t0<t<T, когда оно означает эквива-
лентность (взаимную абсолютную непрерывность) их
102
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
распределений вероятностей и в том или ином
функциональном пространстве *).
В связи с этим стоит сказать, что рассматривае-
мый нами вопрос о структуре обновляющих процес-
сов касается лишь тех свойств случайных процессов,
которые целиком определяются их вторыми момен-
тами (корреляционной функцией), так что можно
было бы, не ограничивая общности, считать эти про-
цессы гауссовскими.
Итак, имеет место следующее предложение.
Теорема. Эквивалентные случайные процессы
имеют обновляющие процессы одного и того же типа.
§ 2. Эквивалентность стационарных процессов
1. Постановка вопроса. Операторы Гильберта —
Шмидта в функциональном пространстве L2(jR).
Пусть £(/) = {£(/), x}x^R — стационарный процесс с ком-
понентами {£(/), х}, хе R (/? — параметрическое гиль-
бертово пространство), который мы рассмотрим на
некотором интервале t$<l <Т.
Обозначим fK его спектральную плотность. На-
помним, что fK — положительная операторная функция
в гильбертовом пространстве R такая, что при всех
х, у е= R
оо
Е(т х} • = J e^^{fKx,y]dK.
—-оо
Наряду со спектральной плотностью мы будем
использовать и другую характеристику стационар-
*) Это условие эквивалентности гауссовских распределений
было предложено Фелдманом (J. Feldman, Equivalence and
perpendicularity of Gaussian processes, Pacif. J. Math. 8 (1958),
699—708). Подробное изложение разных вопросов, связанных
с эквивалентностью гауссовских распределений, имеется в мо-
нографии Ю. А. Розанова, Гауссовские бесконечномерные
распределения, Труды МИАН 108, (1968).
§ 2] ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ ЮЗ
ного процесса £(/)— его корреляционную функцию
оо
В6(0 = J etMfKdh,
—оо
-- ОО < t < оо,
являющуюся (слабым) преобразованием Фурье (слабо)
интегрируемой операторной функции /А, — оо < Л< оо.
Пусть ц(/) = {ц (/), — Другой стационарный
процесс со спектральной плотностью gK и корреля-
ционной функцией Вп(/) в параметрическом гильбер-
товом пространстве R. Спрашивается, при каких
условиях на спектральные плотности fK и gK (или
на корреляционные функции (/) и Вп(/)) процессы
£(/) и т](/) будут эквивалентными на интервале
/о < t < Г?
Имея в виду изометрию
I (/) е{М^2, т] (/) eiMgW
(см. § 1 гл. III), рассмотрим в гильбертовом про-
странстве L2 (7?) подпространства Н (f) и Н (g), ка-
ждое из которых есть замыкание соответствующих
линейных подпространств всех функций вида ?А/2х(Л)
И g'^xft), где
х(Л) = 2 e‘KtkXk, — оо < Л < оо; (2.1)
*0<^<Г
именно,
H(f) = \/ eiKtf^2R,
h<t<T
H(g) = V (2.2)
Эквивалентность процессов и т](/) на интервале
tQ<t<T означает, что отображение
Л:
(2.3)
есть линейный ограниченный обратимый оператор
из гильбертова пространства Н (g) в гильбертово
104
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
пространство H(f), а разность / — А* А— оператор
Гильберта — Шмидта *).
Прежде чем говорить о тех или иных условиях
эквивалентности, покажем, что всякий оператор Гиль-
берта— Шмидта в функциональном пространстве
L2 (7?) в некотором смысле есть «ядро» — факт, хорошо
известный в конечномерном случае dim7?<oo.
Обратимся к гильбертову пространству S2(/?) всех
операторов Гильберта — Шмидта в гильбертовом про-
странстве R, определив скалярное произведение ф,
ф е S2 (7?) как
{(р, ф} = Sp {<реА, =
k
= 3 {ф^, в/} et}, (2.4)
i
где eh e2t ••• —ортонормированный базис в /?. Хо-
рошо известно, что стоящие в правой части равен-
ства (2.4) суммы не зависят от выбора ортонорми-
рованного базиса в 7? и что пространство S2(T?) яв-
ляется полным относительно «следовой нормы»:
I <рI2 = Sp <p*<p = S II <f>ek II2 = Si {ф*?ь Sj} I2-
k k, /
Рассмотрим пространство всех измеримых опера-
торных функций ф(Л, ц), — оо < Л, ц < оо, в гиль-
*) Отметим, что излагаемые ниже методы и результаты
легко распространяются на операторные функции и «мед-
ленного роста»:
оо оо
J 1^Г2"||МИХ<оо, j |X|-2n||gJ|dA< оо
— оо —оо
при каком-либо п > 0, что соответствует случаю обобщенных
стационарных процессов и процессов со стационарными при-
ращениями. Вместо исходных функций х (Л) типа (2.1) нужно
лишь взять функции
х (Л) = JeiUc (t) dt,
h
являющиеся преобразованием Фурье достаточно гладких функ-
ций с (/), — оо < t < оо, обращающихся в 0 вне интервала (/q, Т).
§ 2] ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ Ц)5
бертовом пространстве /?, удовлетворяющих условию
ф (Л, ц) е S2 (/?) и J* J |ф(Л, р)|2 dK d\x < оо, (2.5)
— 00 —оо
определив в нем скалярное произведение формулой
(ф, ф) = / J {ф (X, ц), ф (X, ц)} dX (2.6)
— оо —оо
Это будет полное гильбертово пространство; обозна-
чим его L2 (R X R)-
Примером операторной функции ф(Х, [i)eA2(/?X^)
может служить функция
ф (Л, ц) = х0 (Л) X Уо (н)> — °° < X, ц < оо,
определенная соотношениями
ф(Х, ц)х = х0(Х) -{t/o(H)> х}, xe=R, (2.7)
где х0 (Л), у0 (%) s L2 (R). Действительно,
| ф (X, и) I2 = II х0 (X) II2 21 {//0 (ц), ek} |2 = II х0(Х) ||2 • || Уо(ц) II2
k
и
J j Iф (X, ji) I2 dX dp — j II Х(Д) II2 dK- | II г/0(р) Ipdp.
— oo —oo —oo —oo
ЕслихДХ), x2(X), ... —некоторая ортонормирован-
ная система в функциональном пространстве L2(R),
то система операторных функций
Фй/(Х, Н) =-Vfe (X) X (ц), k, ] = 1,2.... (2.8)
будет ортонормированной в пространстве LPtRXR),
поскольку
{ф//.(^> Н), ФР<7 (Л, р)} =
— {Xi.(X); хр (X)} s {Xj (р), ek} • {xq (ц), ek} =
= [xt (X), хр (X)} • {xt (ц), xq (р.)1
106
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
И
оо оо
(ф</> фр9) = J J {ф;/ (X, и), фр, (X, ц)} d'k du =
= J {xz (X), (X)} </Х • j {Xj (ц), xq (р.)} dp = 0
при i р или j У= q.
Более того, используя тот факт, что всякий опе-
ратор Гильберта — Шмидта есть предел по следовой
норме «конечномерных» операторов *), можно пока-
зать, что если ортонормированная система функций
Xi (X), х2(^)> ••• полна в £2(/?), то соответствующая
система операторных функций qki (Л, ц) = xk (X) X
Х^(ц), £,/=1,2, ..., будет полной ортонорми-
рованной системой в пространстве L2(RXR) (это
предложение нам не понадобится, и его доказатель-
ство мы опускаем).
Теорема. Всякий оператор Гильберта — Шмидта
в функциональном пространстве L2 (R) задается не-
которым ядром ф(Х, ц) е L2 (/? X R)-
Доказательство. Пусть F—некоторый опе-
ратор Гильберта — Шмидта (в некотором подпрост-
ранстве Н L2 (R)). Это значит, что
Sl (Fxk, х}) |2 < оо (2.9)
k, j
для полной ортонормированной системы функций xk —
= хЛ(Х), £ = 1,2,.,., в пространстве H^L2(R).
Возьмем ортонормированную систему операторных
функций
ФА/ (X, ц) = xk (X) X Xj (ц), k, j= 1, 2, ....
и положим
Ф (X, Н) = 3 bkjfpkj (X, И),
k,i
*) См., например, книгу Н. Данфорда и Дж. Шварца,
Линейные операторы. Спектральная теория, М., изд-во «Мир»,
J966,
§ 2]
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
107
где коэффициентами служат
ьк{ = х}), Si bkl |2 < оо.
/
Ясно, что <р(Л, ц) есть элемент (полного) функцио-
нального пространства Л2(/?Х^)- При всех k, j и р
Фй/ н) Хр (ц) = Xk (AJ • {хj (ц), Хр (ц)} еЕ /7
и
Фй/ (*> Ц) хр (ц) dp =
при / = р,
при j=£p.
xk W
О
Легко видеть, что операторная функция ф(Л, р)
обладает следующим свойством: для всякой функ-
ции х (Л) = S архр (Л) (Sl«pl2<°°) из пространства
р \ р )
Hs=L2(R)
оо
J ф (х, ц) х (ц) dp = ckxk (Л),
— оо k
где
Сй = 2«рбРй, Si cki2<Siарi2-Siм2<00•
р k р k, f
Можно сказать, что оператор Ф, определенный в про-
странстве L2(R) равенством
G>x(A.) = S<Vft W
h
задается ядром ср (Л, ц) е А2 (7? X Я), и этот опера-
тор Ф совпадает с исходным оператором F, поскольку
(Л) = 2 W = 2^/2 xk W =
i i k
= 2 xk = 2 Ckxk (^)«
k \ j ) k
Итак, мы показали, что всякий оператор F типа
Гильберта — Шмидта в функциональном пространстве
108
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
Н L2 (7?) задается некоторым ядром ф (Л, р) е
eL2(/?X^)> а именно,
Fx(k) = J ф (Л, p)x(p)dp. (2.10)
— оо
В свою очередь, как легко видеть, всякое ядро
ф(Л, ц) е L2(R X R) задает некоторый оператор F
типа Гильберта — Шмидта в пространстве А2(/?)>
причем
|F|2=Spr’F= J J |<р(Л, (2.II)
— оо —оо
В дальнейшем нам понадобятся еще некоторые
свойства преобразования Фурье операторных функ-
ций ф (Л, ц) е L2(/? X RY определенного равенством
b(s,f) = J J е-'<Л5-^>(р(Л, ix)dXdn, (2.12);
— оо —оо
— оо < s; / < оо.
В случае интегрируемой функции |ф(Л, ди J f это
обычный интеграл от функции со значениями в се-
парабельном гильбертовом пространстве S2(/?). В об-
щем случае преобразование (2.12) можно определить
в смысле слабой сходимости. Именно, если выбрать
какой-либо ортонормированный базис eh е2, ... в гиль-
бертовом пространстве 7? и рассмотреть обычное
преобразование Фурье интегрируемых в квадрате
функций {ф(Л, ^)е^, ef}, обозначив его {b(s, f)ek, е;},
то окажется, что для любых линейных комбинаций
х = 2 c'kek> 7/ = 2 c"ej функция
{b (s, t) X, у\ = 2 c'kc" {b (s, /) ек,
совпадает с преобразованием Фурье {ф(Л, //},
и потому мы можем рассматривать это преобра-
зование как билинейный функционал на всюду плот-
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
109
§ 2]
ном в 7? подпространстве 7? всех линейных ком-
бинаций вида x — ^jckek. При этом
J | 5\{b(s, f)ek, ej}\2dsdt =
— оо — оо
оо оо
= S / J |{^(s> t)ek, et}\2dsdt =
k, f —oo —00
oo oo
= 4jl2S I J |{ф(Л, p)eA, еДРrfXcfp. =
k, j
= 4n2 J J | {<p (Л, ц)ек, e/}l2d^dp =
— oo —oo k, i
= 4n2 J J |<р(Л, p) p dk dp < oo,
— 00 —oo
откуда видно, что
2 |{6(s, t)ek, в/}|2< oo п. в.,
k, i
и следовательно, билинейный функционал {b(s, t)x, у]
задается некоторым оператором b(s, t) типа. Гиль-
берта— Шмидта (в пространстве 7?). Как оператор-
ная функция b(s, t) принадлежит функциональному
пространству L2(R%R):
J j \b(s, t)\2 ds dt = 4n2 j J |ф(Л, р.)|2(Д dp < oo.
— OO —OO —00 —oo
(2.13)
Операторную функцию b(s, t), — oo < $, z < oo, мы
и называем преобразованием Фурье от <р(Л, ц) е
е£2(/?Х^).
В свою очередь, операторная функция <р (X, ц)
может быть получена из операторной функции
по
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕСС^
[ГЛ. IV
#($, t) <= L2(R X R) с помощью обратного преобра-
зования Фуръе\
оо оо
<р(Л, Н) = т”? / / {Ks~'xt) b(s, t)ds dt, (2.14)
— оо < Л, [х < оо.
Используя хорошо известные свойства преобразова-
ния Фурье скалярных функций, легко убедиться,
например, что если для некоторого натурального п
оо оо
J J (1+ЛТ(1+нТ|ф(^, ц)|2^н< оо, (2.15)
— оо —оо
то операторная функция b(s, f) имеет все (слабые),
производные до порядка п—1 по каждому перемен-
ному, причем
-h^b(s’z)=
dsRdtm
= J J (— iX)k (щ)т e"“l(Ks^t} ф (Л, ц) dh
— oo —oo
^2(rt-l)
и, более того, функция c(s, t) = b (s, ^аб-
солютно непрерывна относительно лебеговой меры
ds^dt в том смысле, что аддитивная функция
множеств
m(A) = c(s2, t2)~ c(s2, tl) — c(sl, t2) + c(sl, tt),
A = (s1, s2]X(/i, t2],
представима в виде интеграла
m(A)== [ f [-?^(М)
J J ds dt
Si ti
ds dt
от функции, которую мы обозначили
д2п
dsndtn
b (s,
О,
являющейся преобразованием Фурье операторной
§ 2]
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
111
функции ф(Л, ц) = (— /V) (/ц)" <р (Л, ц) е L2(R\ R) и,
следовательно, удовлетворяющей условию
d2n I2
dsdKoo. (2.16)
С другой стороны, если операторная функция
b(s, /), скажем, является финитной: b(s, /) = 0
вне некоторого квадрата tQ< s*, t < Т и имеет про-
изводную b(s, О, удовлетворяющую условию
л2"
(2.16), т. е. ----- b (s, f) Л2 (7? X 7?), то все млад-
шие производные (включая саму функцию b(s, /))
также принадлежат пространству L2(R\R), На-
пример,
т т
I2
ds dt
T T
I J dsndtn
2
&(s, t) dsdt.
При этом обратные преобразования Фурье оператор-
^2/г
ных функций ""d~ndtn ^)и /) связаны равенством
ф(А, ц) = (— /Л/ХфГф^, ц),
которое получается, как и для скалярных функций,
с помощью повторного интегрирования и интегриро-
вания по частям тождества
(/А)-/г(— ц) =
— OQ
д2п
—-r—^b(s, f)dsdt.
ds dt
112
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
Таким образом, операторная функция ф(Л, р), опре-
деленная формулой (2.14), удовлетворяет условию
оо оо
f / |лПн12п|<р(^ н)12^н =
—оо —00
оо оо
= J J | г|5 (X, р) |2 dk dp, < оо, (2.15z)
— 00 —оо
а значит, и условию (2.15).
В заключение приведем еще один простой факт,
касающийся операторных функций ф (Л, ц) е £2(/?Х^)«
Именно, пусть измеримые операторные функции аА
и Pv — оо < Л < оо, в пространстве R удовлетво-
ряют условию
оо ОО
J || ajl2 сА < оо, j || ||2 dk<оо.
—ОО —00
Тогда произведение (Л, р) ац также принадлежит
пространству L2 (R X R)> и для всех х, у е R
ОО
I ф(Л, p^xdp, frKy
—оо
dk —
оо оо
= J J {₽1<Р (Л, ц) z/} dk dp. (2.17)
—оо —оо
2. Общие условия эквивалентности. Как было
отмечено, для стационарных процессов g (/) и т](/) со
спектральными плотностями fK и gK (в гильбертовом
пространстве /?) эквивалентность на интервале t0 <
< t < Т означает, что определенный соотношением
(2.3) оператор Л из пространства Н (g) в простран-
ство Н (f) является ограниченным и обратимым,
а разность / А*А является оператором Гильберта —
Шмидта в Н (g). Уточним, что оператор А:
Аи(к) = ^2х(Л), — оо < Л < оо?
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
113
§ 2]
первоначально определен на всех функциях
«(Л) = ^’х(Х)
/х(Л) = 2 elUkXk, — оо < Л < ooY
\ *<4<т /
причем сами пространства Н (g) и Н ([) совпадают
с замыканием в L2 (R) подпространств из соответст-
вующих функций £'(!х(Л) и (Л) (см. (2.1) — (2.3)).
Будем предполагать в дальнейшем, что спект-
ральная плотность g(K) равномерно интегрируема:
сю
/ UWII<a<oo. (2.18)
—оо
Положим В —Л* Л и F — I — В. Ясно, что усло-
вие ограниченности и обратимости оператора Л равно-
сильно тому, что теми же свойствами обладает опе-
ратор В. Далее, если F = I — В есть оператор Гиль-
берта— Шмидта, то В = 1 — F будет ограниченным
оператором, причем условие обратимости положи-
тельного оператора В = Л*Л:
inf (Ви, и) = 1 — sup (Fa, и) > О
М=1
равносильно тому, что самосопряженный оператор F
(типа Гильберта — Шмидта) не имеет собственных
значений, равных 1 — его максимальное собственное
значение меньше 1,
sup (Fu, и) < 1. (2Д9)
Пи [1=1
Уместно отметить, что ограниченность и обратимость
исходного оператора Л означает следующее *):
J {fjxW.jeWJdXX J {g(X)x(A), x(A)}rfA, (2.20)
—оо —оо
*) Здесь и далее запись аХР для переменных а, Р озна-
чает, что
0 < Ci Су <с2 < оо,
где с 1, с2 — некоторые постоянные.
114
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
где х (Л) пробегает все функции вида х (Л) =
= 2 e^xk.
tk<t<T
Как мы показали выше (см. теорему п. 1), опе-
ратор Гильберта — Шмидта в подпространстве Н (g) =
^L2(R) задается некоторым ядром <р (Л, ц)еЛ2(7?Х^)>
что в отношении оператора F — I — B равносильно
следующему: для полной в Н (g) системы функций
и (Л) = e'Mg^x, tQ<t < Т, хе R,
и некоторого ядра <р (X, ц) е L2 (R X R) должно иметь
место равенство
оо
(Fu, f) = J
оо
J (Л, ц) g'fx dp,, eiKsg'jVy
—<x>
dK =
oo oo
= J J e-!'(Xs-|lZ'{4/2cp(Z, y] dXdp.
—oo —oo
(см. формулу (2.17)). С другой стороны, для u(X) =
= elKsg^2x, v(K) = eiKtg^2y имеем
(Fu, v) — (u, v) — (Bu, v) —
— J eiK(s-t) (g^x, y] dK— j elKi-s~t} (fKx, y] dK =
—oo —oo
oo
= J ~ x> y} dK = b(s — t),
— oo
где b(s — t) есть разность корреляционных функций:
b (s — t) = (s — t) — B% (s — /), tQ < s, t <T.
Сравнивая приведенные выше выражения, при-
ходим к следующему результату.
Лемма. Разность I — В является оператором
Гильберта — Шмидта тогда и только тогда,. когда
§ 21
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
115
разность корреляционных функций допускает пред-
ставление
b(s — t)=\ / p)gWdkdp, (2.21)
—оо —оо
Zo < /, s < Г,
с помощью некоторого ядра q? (Л, ц) е L2 (R X R)»
Здесь подынтегральная функция
Ф(^. P-) = gl/2<₽(^ ^g’/2
принадлежит пространству L2 (R X R), поскольку
| ||gj,/2|f2^ < 00 (см. условие (2.18)), и формула (2.21)
—оо
показывает, что операторная функция b(s — /),
tQ < s, t < T, может быть продолжена в некоторую
функцию
оо оо
b(s, t) = J J e-z ф (A, p)dAdp,
—oo —oo
— oo < s, t < oo,
преобразование Фурье которой
oo oo
ф(А, н) = "4^2’ J J е1{Кз~»*> b(s, t)dsdt,
— 00 — oo
— OO < л, Ц < oo,
удовлетворяет условию
g^ (*. H) g;1'2 e L2 (R x R). (2.22)
Очевидно также, что если функция b(s — /), /0 <
< s, t < Т, допускает подобное продолжение, то она
допускает и представление типа (2.21) с ядром
<р (А-, и) = grI/24> ц) g^1/2-
116
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
Таким образом, имеет место следующее пред-
ложение *)
Теорема. Разность I — В является оператором
Гильберта — Шмидта тогда и только тогда, когда
разность корреляционных функций b(s — f), tQ<s,
t < Г, допускает продолжение в некоторую функцию
b(s, f), — оо < s, /<оо, из L2(/?X^?) с преобразо-
ванием Фурье ф(Л, ц), удовлетворяющим условию
(2.22).
Как это ни кажется странным с первого взгляда,
в случае конечного интервала (/0, Т) при весьма
широких ограничениях на спектральные плотности
fK и gK из условия, что разность I — В является
оператором Гильберта — Шмидта (и даже просто
самосопряженным вполне непрерывным оператором),
вытекает обратимость оператора В, и в этом
случае предложенная выше теорема дает, как мы
убедимся в дальнейшем, весьма эффективный кри-
терий эквивалентности.
Предположим, что спектральная плотность fK
является невырожденной хотя бы на каком-либо
множестве положительной меры. Предположим, что
спектральная плотность gK является регулярной,
точнее,
J iOg1ll+J dX>-oo (2.23)
(ср. с условием (4.12) гл. III). Тогда имеет место
следующее предложение.
Теорема. При условии, что разность I — В
является вполне непрерывным оператором, оператор
В —А*А будет обратимым.
*) Эта и последующие теоремы являются обобщением ре-
зультатов для одномерных процессов, подробное доказательство
Которых имеется в книге И. А. Ибрагимова и Ю. А. Ро-
занова, Гауссовские случайные процессы, М., изд-во «На-
ука», 1970.
§21 ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
117
Доказательство этого предложения основано на
том, что при условии (2.23) всякий элемент u(k)^H (g)
представим в виде
и (Л) = §^/2х(Л), — оо < Л < оо, (2.24)
где х (К) —целая аналитическая функция, являю-
щаяся среднеквадратичным пределом с весом ш(Л) =
===||£Г1|Г1 аналитических' функций типа (2.1) :
х (Л) = lim S емЬхь.
Действительно, всякий элемент ц(А)е//(§) по
определению пространства Н (g) есть предел в L2(R)
некоторой последовательности ип (Л) = g]K12 хп (Л), где
хп (Л) — аналитические функции типа (2.1). Последо-
вательность ц„(Л) является фундаментальной в про-
странстве L2 (R), так что
сю
JH/2 - Xm
— oo
oo
> J IK — xm (в ||2 ® (k)dk -> 0 при n,
— oo
поскольку для любого x^ R
II 0-1/2 II2
OO
и следовательно, существует функция х (Л), J || х(Л.) Ц2 X
X (Л) dk < оо, такая, что
Оо
lim [ || хл (Л) — х (Л) |F w (Л) dk = 0.
п -> 00
— оо
Таким образом, для некоторой подпоследовательности
имеем
xn^)~>x{k), gl/2xn(k)-^g>/2x(k) п. в.,
118
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ ,[ГЛ. IV
что приводит нас к формуле (2.24). Поскольку при
любом х е R
оо
lim [ | {хп (А), х}— {х(А), х} |2да (А)</А = О,
П > оо J
— оо
скалярная функция {х(А), х} как предел функций вида
{хп(А),х}=
(0 < Ч < т
является целой аналитической функцией*). Например,
при условии
(2.25)
где п — некоторое целое положительное число (что
несколько сильнее, чем условие регулярности (2.23)),
для соответствующих функций х(Л) можно дать сле-
дующее представление**):
п-1 т
X (А) = (iK)keiU'xk + (14- г'А)” J емс (t) dt,
k^O О
где функция c(f) принадлежит пространству L2(R).
Итак, пусть пространство Н (g) состоит из эле-
ментов вида и (Л) = g^2x (Л), где х(Л) — целые анали-
тические функции.
*) Подробное доказательство этого имеется в книге, цити’
рованной на стр. 95. Описание аналитических функций х (Л)
в скалярном случае дается в заметке М. Г. К р е й н а, Об основ-
ной аппроксимационной задаче теории экстраполяции стационар-
ных случайных процессов, ДАН СССР 94 (1954), 13—16; см.
также: N. Levinson, Н. Mckean, Weighted trigonometrical
approximation in with application to the germ field of stationary
Gaussian noise, Acta Math. 112, 1—2 (1964), 99—143.
**) Отметим, что это представление может быть использо-
вано для нахождения проекторов Qt на подпространства
^U)= v eiKsg^R, t0<t<T
<C t
(см. задачу о прогнозе на конечном интервале в книге Ю. А. Р о-
з а н о в а, Стационарные случайные процессы, М., Физматгиз,
1963).
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
119
§ 2]
Предположим, что оператор В = не являет-
ся обратимым. Тогда, как мы показали выше (см.
(2.19)), самосопряженный вполне непрерывный опе-
ратор F — I— В должен иметь равное 1 собственное
значение. Обозначим и (Л) е/7(g) соответствующую
ему собственную функцию. Возьмем последователь-
ность аналитических функций хп(Л) типа (2.1) таких,
что и (Л) = g^2 х (Л) есть предел элементов ип (Л) =
= g^xn (Л) и х(Л) — lim хп (Л) п. в. Получается, что
fl!2хп (^) /V2 х (М п.в. и одновременно последователь-
ность элементов fH2xn (К) = Аип (Л) такова, что
(Аип, Аип) === (Вип, ип) (цп, ^п) (Fun, ип) > О,
поскольку ип->и, Fun->u. Таким образом, в Л2(/?)
последовательность f\!2xn(h) сходится к 0, так что
Д/2х(Л) = 0 п. в.,
и, следовательно, х(А,) = 0 на множестве положи-
тельной меры (где ^/2не вырождена), а для анали-
тической функции х(Л) это возможно лишь в случае,
когда х(Л) = 0. Но это противоречит тому, чго и(Л) =
= ё\!2х G4— собственная функция с единичным
собственным значением, и значит оператор В = А'А
является обратимым.
3. Некоторые выводы и примеры. Покажем, что
предложенные в п. 2 общие условия эквивалентности
стационарных процессов £ (/) и г](/), /0 < t < Т, со
спектральными плотностями tK и gK в гильбертовом
пространстве В являются достаточно эффективными.
Для- краткости назовем спектральные плотности
Д и g^ эквивалентными, если эквивалентны соответ-
ствующие процессы £(/) и г) (/), /0 < t < Т.
Предположим, что
(2.26)
Тогда, как это следует из условия (2.22) (см. тео-
рему на стр. 116), для эквивалентности необходимо,
120
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
чтобы разность корреляционных функций
оо
b(s—f) = B^fs—f) J elK (s ~ ° (^ - A) dK
—oo
на квадрате /0 < s, t < T совпадала с некоторой опе-
раторной функцией &(s, /), — oo < s, t < oo, из про-
странства L2 (R X Я) с преобразованием Фурье q) (Л, p),
удовлетворяющим условию (2.15):
oo oo
| J (1 + A.2)"(l + ц2)"|ф (X, n)|2dXdp. < oo,
— oo —oo
поскольку
и)г-1'»Г>+1Т(1 +|?)"|ф(1, Н)|!.
Действительно, в рассмотренном случае
у х ||2 = || g'fg^x II2 < к (1 + кТп К1/2*||2
|£Г,/2ф(*. н)£Л2 = SК1/2фн)^'Ч||2 >
> 2 А-1 (1 + IIФ (*, Н) g»1/2 хр |(2 =
= /Г'(1+л2)га|ф(^ н)£7'/2|2 =
ц)’Г=
=А-1 (1 4- к2)п 2II г;1/2ф (К, Ц)’ хр II2 >
> /С~2 (1 + Л2)" (1 + ц2)" 2IIФ (л, I*)’ ХРII2 =
иП! =
= к~Ч1 +>.!)’(1+й’1ф(1.н)|2.
Так как b (s, /) = &($ — /) при /0 < s,? t < Т, из усло-
вия (2.15) вытекает, что для эквивалентных спек-
тральных плотнрс.т.еИ и gK операторная, функция.
§ 2]
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
121
b (t) на интервале (— т, т), где т =Т — tQ, должна
иметь все производные до порядка 2п, причем 2п-я
производная (/), —- т < / < т, должна удовлетво-
рять условию
т т
j J | b{2n> (s — /) I2 ds dt < ОО. (2.27)
Предположим теперь, что выполняется соотноше-
ние (2.25), иначе говоря,
+А2)""/. (2.28)
Если операторная функция £(/), /0 < t < Г, обладает
описанными выше свойствами, а именно, существует
производная -^nbif), удовлетворяющая условию
(2.27), то операторную функцию b(s, t} = b (s — t)
можно продолжить с квадрата /0 < s, t < Т на всю
плоскость — оо < s, / < оо в некоторую финитную
гладкую функцию b(s, t) такую, что*)
оо оо
/ J I ds^dtn b(S’ dsdt < 00 •
—oo —oo
Преобразование Фурье if (Л, p.) такой функции b(s, t)
будет удовлетворять условиям (2.15') и (2.15), а сле-
довательно, и условию (2.22)
оо оо
/ / | (*-> н) ^1/212 dk dp. < оо,
*) В ортонормированием базисе еь e2i ... R операторная
функция b (s, t) со значениями в S2 (R) задается некоторой
матрицей bkj (s, О, S I (5> О |2< 00 и речь идет об очевид-
fc, /.
ном продолжении каждой скалярной компоненты bkj ($, /)
с квадрата tQ<s, t<T на всю плоскость — oo<s, /<оо в фи-
нитные функции bkj (s, t) с сохранением свойств гладкости и
оо оо
интегрируемости типа J J | bkj (s, t) |2 ds dt < oo.
— 00 —OO kt j
122
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
поскольку в нашем случае
|^(Х, н)^-,/2|2<Ьг,/2||-к71/21|-^(^ н)12<
<Я2(1+Л2)П(1 +fxT|№ н)12.
Поэтому разность / — Л*Л будет оператором Гиль-
берта— Шмидта, что при условии (2.28) равносильно
эквивалентности спектральных плотностей gK и fK
(если только fK не является вырожденной п. в.).
Таким образом, если разность корреляционных функ-
ций на интервале (—т, т), т=Т — /0, имеет произ-
водную Ь(2п) (/), удовлетворяющую условию (2.27), то
спектральные плотности fK ugK6ydyT эквивалентными.
В частности, спектральные плотности fK и gK будут
эквивалентны, каково бы ни было п в условии (2.28)
и каков бы ни был конечный интервал (tOf Т), если
fk = ёх +
где ^ — произвольная финитная функция с интегри-
руемой в квадрате следовой нормой’.
а
Дх = 0 при Р|> a, J |Дх|2<а< ОО. (2.29)
Действительно, в этом случае операторная функция
оо
b (f) = — J eiU Дх d)., t0< s, t < T,
“oo
будет иметь производные всех порядков, которые
удовлетворяют условию (2.27).
В конечномерном случае, dim/?<oo, получаем,
например, что если спектральная плотность fK как
угодно отличается от спектральной плотности gK на
каком-либо конечном интервале (в частности, может
быть Д=0 при | X |*Са), то fK и gK эквивалентны.
Отметим также, что для спектральной плотно-
сти gK типа
+ +V)-"/
ЭКВИВАЛЕНТНОСТЬ СТАЦИОНАРНЫХ ПРОЦЕССОВ
123
§ 2]
вместо условия (2.29) можно потребовать, чтобы
оо
J |^-,/2W/2|2^<°°- (2.29')
—ОО
Сформулируем наш основной результат в виде
отдельной теоремы.
Теорема. Если gK К (1 + /, то для экви-
валентности fK и gK на интервале (/0, /0 + т) необхо-
димо, чтобы операторная функция
оо
b (/) = Въ (/) - Bl (t) = J е™ (gK - fK) dX
—oo
на интервале — т < t < т имела все производные до
порядка 2п, причем
т т
J 11 b(2n] (s — /) |2 ds dt < оо.
о о
Если же gK k (1 + А,2)”"/, то указанные условия
достаточны, для эквивалентности.
В заключение, как следствие полученных выше
результатов, дадим описание обновляющего процесса
для класса стационарных процессов g(f), /0</<Г,
с корреляционной функцией В (/) следующего типа:
на открытом интервале (0, т), х = Т — /0, функция В(/)
имеет все производные до порядка 2п, причем функция
B(2rt“b(0 непрерывна всюду на интервале (—т, т),
кроме точки t = 0, где она имеет разрыв первого рода,
а функция В(2п) (/) удовлетворяет условию
т %
J j | В(2п) (t — s) |2dsdt<<x>.
о о
Каждый такой процесс £(/) эквивалентен стационар-
ному процессу г)(/), tQ<t<T, со спектральной плот-
ностью вида
g W = । р (а) |2 -
где P(iX)— полином степени п с соответствующим
124
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
коэффициентом при старшей степени (г’Л)п. Следова-
тельно, обновляющий процесс X(t) ~{Xk (/)}"+1, 4 <
< t < T, для g (/) будет такого же типа, как и для
процесса т](/), а именно, первые п компонент будут
постоянными:
£=1, и;
они порождают «начальное» подпространство *)
П Ht.+h (g),
h >0
а одна компонента Xn+i (/), описывающая дальней-
шее обновление подпространств
является процессом типа броуновского движения
(см. § 3 гл. I).
§ 3. Случайные процессы,
эквивалентные винеровскому процессу
С точки зрения нашего подхода к вопросу об эво-
люции подпространств tQ<t<Ty связанных
с соответствующим случайным процессом £(/), tQ <
< t < Т, винеровский процесс г] (/), tQ < t < Т (с компо-
нентами {ц(/), х}, t$<t<Ty где параметр х пробегает
некоторое гильбертово пространство R) может быть
взят в качестве определенного стандарта. Именно,
естественная изометрия
{ц(0, x}«->&(s)x (3.1)
между величинами {ц (/), х} е Н (ц) и элементами
(s) х s L2(R) позволяет отождествить ЯДц) с под-
пространствами L2(R), t0<t<T\ здесь скалярная
функция %i(s), есть индикатор интервала
(/0, f], L2(R) — гильбертово пространство во всех изме-
римых /?-значных функций u(s), t0<s<T, с интегри-
руемым квадратом || и (s) II2, со скалярным произве-
дением
т
(и, о) = J {«(s), v ($)} ds, и, vsL2 (/?),
^0
*) При указанных выше условиях процесс J (t) имеет
ровно п — 1 производных.
§ 3]
ПРОЦЕССЫ, ЭКВИВАЛЕНТНЫЕ ВИНЕРОВСКОМУ
125
a Lt (/?) — подпространство всех функций и (s) е L2(R)
таких, что «(з) = 0 при s>t. Мы считаем, что стан-
дартный винеровский процесс i\(t), t0<t<.T, имеет
нормированную корреляционную функцию
Вп(s, t) — min(s — s0, t —t0) • I, t0<s,t<T. (3.2)
Пусть l(t), t0<t<T,— какой-либо случайный про-
цесс с компонентами {£(/), х}, t0<t<T, хе R, и пусть
В^ (s, t) — его корреляционная (операторная) функция
в гильбертовом пространстве R:
Е {g (s), х} {g (/), у} = {В* (s, t) х, у}, X, yt=R. (3.3)
Как мы знаем, семейство подпространств Ht(g),
t0<i<T, будет изометричным семейству Lt (R),
to<t<T, если, скажем, процесс g(/), /0</<7’, экви-
валентен винеровскому процессу tj(/), ta<t<T. В дан-
ном случае эквивалентность означает, что оператор А:
Л : &(s)x->{g(i!), х}, (3.4)
определенный на полной в L2(R) системе элементов
и(s) == %<(s)х (/0<s<T; хе/?; t0<t<T) (3.5)
может быть продолжен до линейного ограниченного
обратимого оператора А такого, что разность / — А*А
есть оператор Гильберта—Шмидта в функциональном
пространстве L2(R).
Согласно общей теореме п. 1 § 2 всякий оператор
Гильберта—Шмидта в функциональном пространстве
L2(R) задается некоторым ядром К (s, ^)е/.2(/?Х
X R) — операторной функцией в гильбертовом про-
странстве R со значениями К (s, f) е S2 (R) и интегри-
руемой в квадрате следовой нормой:
т т
/ J |K(s, f)l2dsdt< ОО, (3.6)
где
K(S, /)|2 = Sp[K(s, ty . Ш /)].
126
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
Следовательно, при условии эквивалентности
т т
(и, v) - (Au, Av) = J j {К. (s, t) и (s), V (/)} ds dt, (3.7)
^0 ^0
где /С (s, t) — некоторое ядро из L2 (RR).
Очевидно, соотношение (3.7) будет выполнено
для всех элементов и, уе£2(/?), если оно выпол-
няется для некоторой полной системы элементов
(скажем, для системы (3.5)). При и (s) = (s) х,
x)(t) = ^h(t)y соотношение (3.7) можно представить
в операторной форме:
ЗД, h)-BK(tx, f2)=| J K(s, t)dsdt. (3.8)
to to
Очевидно также, что если представление (3.8) имеет
место, то соответствующий оператор / — А* А будет
оператором Гильберта—Шмидта.
Далее, если разность / — В(В —Л*Л) есть опера-
тор Гильберта—Шмидта, то условие обратимости огра-
ниченного оператора В можно выразить следующим
образом: оператор / — В не имеет собственного зна-
чения, равного 1 (см. по этому поводу стр. 113).
В итоге мы приходим к следующему результату.
Теорема. Для эквивалентности случайных про-
цессов и т](/), t0<t<T9 необходимо и достаточно,
чтобы разность b(s, /) = ВП($, /) — B$(s, t) была абсо-
лютно непрерывна относительно dsy^dt, точнее,
s t
b (S, t) = j J К (s', f) ds' dt', t0<s, t<T, (3.9)
to to
d2
где производная Д (s, t) — ds b (s, t) имеет интегри-
руемую в квадрате следовую норму*.
т т
I3-10)
to to
§ 3]
ПРОЦЕССЫ, ЭКВИВАЛЕНТНЫЕ ВИНЕРОВСКОМУ
127
а кроме того, у ядра К (s, t) как оператора в L2 (/?)
нет единичного собственного значения *).
Этот критерий эквивалентности сохраняется и
в случае более сложного «стандартного» процесса
т>(/), tQ<t<T, имеющего корреляционную функцию
вида (s, t) = min (F (s), F (/)) • /, где F (/), t0 < t < T, —
непрерывная слева монотонно неубывающая скаляр-
ная функция. Нужно только вместо лебеговой меры
dsy^dt взять соответствующую меру dF (s) X dF (t).
Пример. Пусть т](0, /0 </< Т, — одномерный
случайный процесс с некоррелированными прираще-
ниями и F (t) = Е |ц(/) |2, tQ < t < Т, — его структур-
ная функция. Пусть a (t), t0<t < Т, — случайный
процесс такой, что
т
J Е\ a(f) \2dF (t) < <*>.
to
Пусть, наконец,
т
Ш = J a(sW(s) + n(0. k<t<T. (3.11)
t,
Учитывая, что проекция величины а (и) на под-
пространство Н t (т\) представима в виде стохастиче-
ского интеграла J с (ц, где подынтеграль-
to
ная функция удовлетворяет условию
t
11 с (и, ti) |2 dF (v) Е | а (и) |2,
to
*) По-видимому, первые публикации об условиях такого рода
для скалярного случая даны Шеппом (L. A. S h е р р, Ra-
don— Nykodym derivatives of Gaussian measures, Ann. Math. Stat.
37, 2 (1966), 321—354). Отметим, что эти условия легко полу-
чаются из условий эквивалентности стационарному процессу
белого шума с простейшей спектральной плотностью g(X) = /.
128
ЭКВИВАЛЕНТНЫЕ СЛУЧАЙНЫЕ ПРОЦЕССЫ
[ГЛ. IV
будем иметь
t
Еа (и) т] (0 =
' t
= Е | с (и, f) dx\ (v) • J dr) (v) = | c (u> f) dF (v).
л
Ч,
Элементарные выкладки приводят к равенству:
s t
(s, 0 — (s, t) = J j к (U, v) dF («) dF (»),
где /((w, v)~—c(u, v) — c(v, u)— E [a (u) a (y)].
Таким образом, при дополнительном условии (3.10)
с заменой dsdt на dF (s)dF (/), случайный процесс £(/),
tQ<t<T, будет эквивалентен процессу ц(/), /0</<Г,
и, следовательно, будет иметь обновляющий процесс
такого же типа, как ц (/), /0 < t < Т. Это условие,
означающее, что оператор
является ограниченным и обратимым, автоматически
выполняется, например, если процессы a(t) и r|(f)
н е к о р р е л и р о в а н ы; Еа (s) ц (/) = 0, /0 < s, КТ.
Действительно, в этом случае ядро К (s, t) =
== — E[a(s)a(/)] является отрицательно опре-
деленным:
2
< О
т т
J j /C(s, t)u(s)int)dF(s)dF(t) =
Т
= —Е J а (/) и (t) dF (t)
т
для любой ненулевой функции u(Z), J| u(t) |2dF(f)<oo.