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

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

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

170 Гл. I I I . Прогнозирование и фильтрация случайных процессов

сходится равномерно при | х | < 1 — б, б > 0. Имеем

N N со

 

 

2 1о° (! + а 2

.

4

2 { 2 ( -

} =

 

 

 

 

 

 

1

1

 

 

 

 

 

 

N

 

 

 

 

 

 

=

Hm

S ( - l

mN. P>

 

 

 

 

N-+OQ

 

 

 

 

и так

как

а < Л/-1, 2л/ (А,) ^

М,

то мы можем также разложить

log (1 + 2ла/ (А)) в логарифмический

ряд н проинтегрировать его

почленно,

откуда

 

 

 

 

 

 

 

со

 

 

со

 

 

Я

 

Ши 2 ( - i r i ^ - m w. p= 2 ( - 1)Р"1- у " ^ г J (2л/(А)}рс/А.

"*

і

 

 

1

 

 

—я

 

Но левая часть равномерно

ограничена величиной

 

 

 

V ( а М ) Р

V

1 1

г .

 

 

 

2 — —

< 2

y i n

< 00

 

 

 

 

 

 

 

 

CU

 

при а< (2М )~ 1. Последовательность

fN (a)= 2 ( — 1)р_1 — - mN, v

равномерно

ограничена при

| а | <

 

 

1

 

(2Л/)-1. Кроме того, для любого

вещественного а в интервале ( — 1/М, 1 /М)

функции fN (а)

сходятся к

 

 

со

 

 

 

л

 

 

 

 

/ ( с с ) = 2 ( - 1 Г 1- у ^ г j {2л/(А)}р <іА.

 

 

 

1

 

 

 

- Я

 

 

Тогда, согласно теореме Витали (см. Титчмарш, Теория функций, ГИТТЛ, М.— Л., 1951, стр. 194), последовательность сходится равномерно в круге | z | < (2М)-1 н, следовательно, / (а) является аналитической. Заметим, что коэффициенты последовательности ана­ литических функций, которая сходится к аналитической функции, сходятся к коэффициентам этой функции. Приравнивая коэффи­ циенты при одинаковых степенях а, получаем следующую теорему:

Т е о р е м а

7.

Если / (А) равномерно

ограничена и Ау-,Лг, / =

= 1, . . ., IV,— собственные

значения ГЛг,

то

 

N

 

 

п

 

lim IV1 2

 

к 2 — \ {2л/ (А)}р dA, р = 0, 1, ....

N-* оо

,

г п

•>

 

 

і

 

 

- я

 

Эта теорема показывает, что величины 2л/ (со7- — л),

с о = 2л//IV,

/ = 0,

1, . . ., N — 1, будут давать в среднем хорошую аппрокси­

мацию

для собственных значений ГЛпри большом N,

в чем можно

3. Общая теория для одномерных процессов с дискретным временем 171

убедиться, аппроксимируя интеграл конечными суммами. С другой стороны, выражение в правой части можно рассматривать как момент порядка р функции 2л/ (•) от случайной величины X, рас­ пределенной равномерно на интервале (—л, л]. Теорема утверж­ дает, что все моменты последовательности серий наблюдений {hj.Ni І = 1) • • м N} сходятся к моментам 2л/ (X ).

(ii) В

упражнении

4 гл. II мы видели, что процесс х (п) =

= cos п \

+ sin mp где

£ и г| — независимые величины с одинаковым

(симметричным) распределением вероятностей F, стационарен и имеет спектральную функцию F. Тогда если F несингулярна и log / интегрируем, то X (п) не может быть безошибочно предсказан с помощью линейного выражения от прошлых значений. Однако х (п), очевидно,

безошибочно предсказуем с помощью нелинейного выражения.

(iii) В заключение этого параграфа, следуя Робинсону [1962], обсудим смысл внутренних функций в теореме 4. Прежде всего для каждого из таких множителей групповая задержка (см. § 4 гл. II) неотрицательна. В самом деле, для ехр іа и zp это очевидно, а для других множителей имеем соответственно

 

i_ _d_

an

1 - K I 2

> 0 ,

 

i dX

1 —a ЛХ )}

I 1 a n e ~ ' I2}‘

 

 

1

d

2 {1 —cos (Ѳ-Х)}

(X(dQ) > 0.

i

dX { - і £ £ и * , ы

I J0_JKu

 

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

Рассмотрим теперь фильтр с частотной характеристикой

 

о

2 Р 0 ‘) еіД-

а2

2

ß O')eijK = /w -

 

Ѵ 2я

 

0

 

о

 

 

Тогда

 

 

 

 

 

S l ß O ) 12 = S « 0 ) 2,

2

Iß 0)12< S « ( / ) 2-

(3.8)

о

о

 

о

о

 

В самом деле, мы можем представить х (іі) в виде

 

X (») =

2 ß (/) ’] (» — У'),

Е (V (,г) Л (т)) = бт ° 2-

 

 

о

 

 

 

 

Тогда a/Fn-i содержится в подпространстве, порождаемом величина­ ми ц (п — /), / > 1, так что I ß (0) I2 а2 — средний квадрат откло­

172 Гл. I I I . Прогнозирование и фильтрация случайных процессов

нения X (п) от его проекции на это подпространство — не превосхо­ дит о21 = а (О)2 о2. Утверждение в общем случае получается совер­ шенно аналогично. Если представить себе, что эти два фильтра имеют один и тот же вход, то оптимальный фильтр приписывает максимально возможный вес самым поздним входным сигналам (при фиксированном коэффициенте усиления). Произвольный односторон­ ний фильтр с заданным коэффициентом усиления получается путем последовательного соединения оптимального фильтра и фильтров, отвечающих внутренним множителям, которые, как мы видели, при­ водят к задержке выходных сигналов.

4 . О Б Щ А Я Т Е О Р И Я Д Л Я С Т А Ц И О Н А Р Н Ы Х О Д Н О М Е Р Н Ы Х П РО Ц ЕССО В С Н Е П Р Е Р Ы В Н Ы М В Р Е М Е Н Е М 1)

Мы не будем проводить здесь подробное рассмотрение, но укажем, каким образом результаты § 3 приводят к соответствующим резуль­ татам для непрерывного времени. Мы следуем Дубу [1953] и Гоф­ ману [1962]. Другой подход излагается в работе Хелсоиа [1964].

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

Это отображение переводит границу (z — exp іХ) в границу = іх), причем X = tg 1/2Х, откуда

1

d X

1

2 л

d x

л (1 + т2)

Основной результат предыдущего параграфа утверждает, что функция / (X) является квадратом модуля голоморфной в единичном круге функции тогда и только тогда, когда ее логарифм интегрируем.

Аналогичная теорема для

полуплоскости, принадлежащая

Пэли

и Винеру [1934], гласит, что спектральная плотность fix),

—сх><

< X < со, представима в

виде

 

/ ( т) = -^- jaOOe* ! ds

 

тогда и только тогда, когда

 

со

 

 

J

 

(4.1)

— со

 

 

Доказательство основывается’ на соответствии между пространства­ ми Н] ІІ = 1, 2) для единичного круга и для левой полуплоскости 2);

1)

Э тот

п ар агр аф м ож н о оп усти ть .

2)

См.

м атем ати ческ ое п р и л о ж ен и е.

4. Общая теория для одномерных процессов с непрерывным временем 173

в случае полуплоскости оии состоят из функций, голоморфных в полуплоскости и интегрируемых или соответственно квадратично интегрируемых на ее границе. Мы используем одни и те же символы Н j в случае круга и полуплоскости, так как аргументы (w или л и т. и.) позволят определить, о каком из пространств идет речь. Легко видеть, что g (w) £ Hj тогда и только тогда, когда *)

а д = (і+2г2/,я (^|)еяц

Если / (т) удовлетворяет условиям теоремы Пэли — Вииера, то

логарифм

функции ср (К) /

(tg 1/2Х) интегрируем,

откуда

/ (т) =

= I g (іт)

I2 и g (w) 6 Hz (так

как g (w), равная h

,

голо­

морфна в полуплоскости и квадратично интегрируема на границе, ибо / (т) интегрируема).

Теперь положим

оо

1' s ii^ e - ^ d x .

—оо

Для доказательства достаточности остается показать, что a (s) = О при s ■< 0. Можно предположить, что g (w) £ H t f| H 2, так как в противном случае вместо g можно рассмотреть функции

[1 — {w/(1 + w)}n] g (w),

которые принадлежат Hi f) П 2 при n > 0 и сходятся в среднеквадра­

тичном

к g. При этом предположении найдется функция Іг 6 / / ь

такая,

что

х)

Здесь,

очевидно, необходимо дать

точное определение классов Харди

H j для полуплоскости:

это пространства

голоморфных в левой полуплоскости

функций,

таких, что

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

J

1h (.т +

іу) р dy ^

const,

 

 

 

 

— ОО

 

 

 

 

 

Одна из теорем Пэли — Винера утверждает,

что h £ # 2 тогДа и только тогда,

когда

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h (ш) =

^ewsa (s) ds,

Re w

О,

 

 

 

 

 

О

 

 

 

ОО

 

 

 

 

 

 

 

 

где J

[ a

(s)

I2 ds <; оо.— Прим,

перев.

 

 

 

----ОО

174 Гл. I I I . Прогнозирование и фильтрация случайных процессов

И при S < 0

СО СО

— f

g(ix) e~üxd% = — f

(1— it)2 xfr~T^~ d% =

л .1

0 ' '

Я J 1 — IT

'

' 1-f- X-

— CO

 

— CO

 

 

Я

=-^- [ k L(e^)e^d.X = 0,

-Л

так как

 

 

 

hi —h(z) exp ( —s i ^ i ) £ //j.

 

 

Необходимость доказывается легко: пусть / (т) имеет указанный

вид;

тогда

/ (т) = | g (іх) |2, где

g (w)

определенно

принадлежит

II2 для полуплоскости;

но тогда

/ (т)

переходит в

cp (А) =

| h р,

h р # 2, так что

log ср (А.) интегрируем,

откуда сразу

следует

необ­

ходимость.

 

 

 

 

 

(4.1),

 

Далее,

если

log / (т)

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

соотношению

то

 

 

 

 

со

 

 

 

 

 

 

 

* ( 0 = j a{t + s)l(ds),

 

 

 

 

 

 

о

 

 

 

 

где I

(s) — как

обычно,

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

такой,

что

Е (Е (tfe)2) =

ds.

 

 

 

 

Теперь,

следуя Дубу [1953], мы построим процесс у (п) с дискрет­

ным параметром, соответствующий процессу х (і) с непрерывным

параметром,

так что подпространство Л 0, порождаемое величинами

X (t), t ^ 0,

в пространстве

(порождаемом всеми х (t)) отобра­

зится В ПОДПрОСТраПСТВО Жп.

ПОООЖТЯРМПР литгпшиачя у (/»), и о,

в пространстве $£■ у (порождаемом всеми у (п)). Пусть спектральная

функция у (п) равна

G (А) =

F (tg Ѵ2А),

где F (т) — спектральная

функция процесса а: (1). Заметим теперь,

что аII0 изоморфно подпро­

странству в

Ь2 (F),

порождаемому

функциями exp (—itx),

t ^ 0,

а J T оизоморфно подпространству в Ь2(G), порождаемому функциями

ехр (—іпХ),

п ^ 0.

Чтобы установить нужное соответствие, рассмо­

трим

отображение

Ь.г (F) на

Ь2 (G), задаваемое формулой

 

 

 

 

 

j ^ ^ e x p t A .

 

(4.2)

Имеем

 

 

 

 

 

 

 

 

 

4 —р -г—=

— 1 - f - ~7~.— =

— 1 “Ь 2 ( е ~ ите 1 d t,

 

 

 

1 + іх

1 1-рix

 

J

 

 

 

 

 

 

 

 

— СО

 

поэтому функция

ехр (—іХ),

а значит,

и ехр (—іпХ) при любом

п ^

0, соответствует

некоторому элементу из wJ(0 (последний инте­

грал

может

быть

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

в среднеквадратичном с

весом

F {dx) липейными комбинациями функций ехр (—Ux), t ^ 0, и поэто­

4. Общая теория для одномерных процессов с непрерывным временем 175

му принадлежит Л 0). Итак, образ jfT0содержится в М й. Обратно, функция

 

 

exp wt = ехр{£ (г — 1 )/(z + 1)}

голоморфна при

| z | >- 1

и не превосходит по модулю единицы при

t <

0. Она разлагается

в ряд по неположительным степеням z.

На единичной

окружности эта функция

 

 

exp it

(tgV2Ä,) = exp itx, t < 0,

n. в.

равна предельному

значению ограниченной функции exp wt,

t <

0. Таким образом, функция exp itx, t <; 0, может быть аппрокси­

мирована ограниченно, а значит, н в среднеквадратичном (по мере F (dx)) образами функций из следовательно, образ j(Fо содержит­ ся в ojHо. Итак, эти два подпространства отображаются одно в другое при соответствии (4.2). Поскольку процесс х (t) чисто детерминирован

(по определению)

тогда и только

тогда, когда у//0 =

е#_те, и а//'.«,

соответствует jp''

 

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

выше)1),

то имеет

место2)

 

 

Т е о р е м а

3'.

Процесс х (t)

чисто недетерминирован в том

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

место (4.1).

 

 

 

 

(Мы

называем

х (t) чисто недетерминированным,

если & # _ —

нулевое

подпространство.)

аналог теоремы 5.

Мы опускаем

Далее проще всего установить

доказательство, так как оно является очевидным видоизменением

доказательства

теоремы

5.

 

 

 

 

 

 

Т е о р е м а

5'. Пустъ х (t)

недетерминированный процесс со

спектральной

плотностью

 

 

 

 

 

 

 

f (х)=

j a

(s) eUxds = | g (ix) |2.

 

(4.3)

 

 

 

о

 

 

 

 

 

 

Tогда для w = a -f- іт

имеем

 

 

 

 

 

 

g (w) = eiaB (w) C (w) Q (tr),_

 

 

 

В (w) = (4 + ^ - )p П {

I 1 — ßr, I

u;+ ß„

I Pn,

 

 

K 1

1 \—

w )

-А-I

\

/

 

 

 

 

 

71

 

 

 

 

 

 

 

 

 

 

CO

 

 

 

 

C (w) = e~pwexp [ —

j

l ~ ™

V (гіт)]

,

 

 

 

 

 

 

—00

 

 

 

 

 

 

 

oo

 

 

 

 

 

Q (w) = exp [-g^-

j log

 

 

(4.4)

*)

См. упражнения к этой главе.

2)

Номера теорем в этом параграфе соответствуют номерам аналогичных

теорем

предыдущего параграфа.

176 Гл. I I I . Прогнозирование и фильтрация случайных процессов

где ßn — нули функции g (w) в левой полуплоскости, р — неположи­ тельное число, а V — сингулярная мера.

З а м е ч а н и е . Множитель ехр (—рш) появляется из-за (воз­ можного) скачка ц (dX) в точках ± я . Как и прежде, функция Q определяется однозначно, если потребовать, чтобы она была голо­ морфна в левой полуплоскости и удовлетворяла соотношению

 

со

е ( - 1 ) = ехр-^

j log/(x)T^ r .

— со

Мы можем также положить

 

со

 

Q{w) = j

а (s) eswds

о

 

для w, лежащих в левой полуплоскости, причем на границе интеграл сходится в среднеквадратичном.

Рассмотрим теперь задачу нахождения минимума выражения

OQ

j I е~іи — h{ix) ^ F (dr),

где функции h (w) должны быть голоморфными в левой полуплоско­ сти. Это есть среднеквадратичная ошибка прогноза х (t) по сУ/0, записанная в спектральной форме. Если F абсолютно непрерывна и отвечает недетерминированному процессу, то это выражение можно переписать в виде

j|е _і(т()(п:)— h {ix) Q {ix) \z dx',

СО

выполнив преобразования, аналогичные использованным в § 1, мы получим, что минимум достигается при

h[{ix) = [e~iixQ (іт)]+ Q (it)“1,

где

СО

с о

[e~iixQ (іт)]+ = I

eisxe~itxa{s)ds— | e’sxa {s 1) ds.

Ошибка прогноза при этом равна

j |[е -і,т№ )]_ |Ч х ,

где

[e~itxQ (іт)]_ =

о

4. Общая теория для одномерных процессов с непрерывным временем. 177

j а1{и) du.

О

Доказательство того, что сингулярная часть F не дает вклада в ошиб­

ку прогноза, почти не отличается от соответствующего доказатель­ ства для дискретного времени, данного в § 3, и мы его опускаем. Теперь мы докажем теорему 2', после чего сформулируем теорему 4' и закончим параграф теоремой 6'.

Т е о р е м а 2'.

г

 

X (t) = и (t) -f V (t) = j a(t — s)£,{ds) + v(t),

(4.5)

где

 

(i) £ {s) —■процесс с ортогональными приращениями,

такой, что

Е(I (dsf) = ds;

(ii)E (ц (s) V (t)) == 0;

(iii)V (t) чисто детерминированный, и (t) чисто недетерми­ нированный процессы.

Разложим в сумму &//_«> ® о/ll, где ыіі — ортогональное допол­ нение е'//_оо в afë. Обозначим через u(t) проекцию х (t) на о/М, а через

и (t) — его проекцию на ы Н Оператор сдвига U (т), U (т) х (t) =

= X (t + т),

коммутирует с проектором на ©//_<», так как, очевидно,

U (т) сИ_оо —

о'і_со, следовательно, он коммутирует и с проектором

на oft, а значит, и (t) и ѵ (t) стационарны. Кроме того, величины и (і) порождают ale, а u(t) порождают сЖ_«>, и (іі), очевидно, выполняется.

Пусть

f “о — пространство, порождаемое

величинами ѵ (t), t

t0.

Тогда

7 'оer о/ІІц, іі поэтому

Т 0 = а#..»,

ибо аМ-оо er J tio,

Т 0 с=

сг с// _оо и о.Иі0 = 7г о +

Тіо,

где *2/о порождается величинам

и (t),

t ^ t0.

Таким образом,

v {t)

— детерминированный процесс. С дру­

гой стороны, обозначая бесконечно далекое прошлое процесса и (t) через ЧІ-оо, имеем 4L-оо с: о/И- « и 4L-оо с: <Л, так что ЧІ-оо есть нулевое подпространство и и (t) — чисто недетерминированный про­ цесс. Поэтому и (t) имеет абсолютно непрерывный спектр и спек­ тральную плотность, удовлетворяющую соотношению (4.1); пред­ ставление (4.5) вытекает из теоремы 15 приложения к гл. II.

Ошибка

прогноза,

очевидно,

равна

 

 

t-j-S

 

f

 

 

 

j

a2(su)du^= j

a2 (u) du,

 

 

s

 

0

 

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

a (s)

действительно

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

Фурье функции Q (£т).

место

 

Так же

как в

§ 3, имеет

 

12 э. Хеннан

178 Гл. I I I . Прогнозирование и фильтрация случайных процессов

Т е о р е м а 4'. Пустъ х (t) не является чисто детерминирован­ ным, и пустъ F = Flli + jFt2>+ F'3>разложение Лебега спек­ тральной функции F\ тогда производная абсолютно непрерывной ком­ поненты F'u равна спектральной плотности процесса u(t) из теоре­ мы 1 a F(2>(X) + F(3' (X) есть спектральная функция процесса v(t).

Т е о р е м а 6'. Пустъ х (і) чисто недетерминированный про­ цесс и Мг>(іх) частотная характеристика фильтра, дающего наилучший прогноз величины х (s -f- it) по наблюдениям до момента s. Тогда

/i(t> (іх)

[е xUQ(ix}]+

Q(it) ’

где Q (w) определяется формулой (4.4). Имеем также

2

j e-sska{s)ds=exp | a 0 + 2 ^ a j

} >

lmu><0,

0 * 0

1

 

 

 

CO

 

 

 

a>= ш j (1о^ 7 м ) ( т т ё ) п- [ ^

-

 

 

— CO

 

 

Пример. Пусть / (т) = (1/2л) | р (іх) |~2, где р (w) — (веществен­ ный) многочлен, не имеющий корней в левой полуплоскости. Тогда, очевидно, Q (w) = р (и;)-1. Имеем

ОО

=і т w e' UxdT==

ОО

hj

впредположении, что корни ßj многочлена р (w) в правой полупло­

скости простые. Таким образом,

i НФз

и, значит, наилучший прогноз х (s -f- t) по наблюдениям до момента s дается выражением

- [ 2 « - v n h ^ r ( £ - h ) } ] * < * > -

іЬФз

5. Прогнозирование векторных процессов

179

5.ПРОГНОЗИРОВАНИЕ ВЕКТОРНЫХ ПРОЦЕССОВ

СДИСКРЕТНЫМ ВРЕМЕНЕМ 1)

Многомерный случай общей теории прогнозирования значительно более сложен, чем одномерный, в силу того обстоятельства, что / (%) может быть вырожденной, но ненулевой матрицей. Тем не менее, благодаря работам Винера и Мазани [1957], [1958], Хелсона и Лауденслагера [1958] и Розанова [1963], существенная часть теории была распространена иа векторный случай. Мы ограничимся дискрет­ ным временем. Относительно случая непрерывного времени см. Робертсон [1968].

Введем гильбертово пространство Ж , порождаемое

величинами

Xj (гг), і = 1, . . ., р; п — 0, ± 1,

. . ., и обозначим через Ж п замкну­

тое подпространство, порождаемое xj (т), /

=

1,

. . .,

р; т $7 п.

Мы получим наилучший ѵ-шаговый линейный прогноз для х (п

ѵ),

если спроектируем компоненты этого вектора

иа

вІ/п.

Ошибку про­

гноза для V = 1 обозначим е (п + 1).

Тогда,

как

и

прежде,

 

Е {е (тгг) е (га)'} =

8^G.

 

 

 

 

 

(5.1)

Полагая

 

 

 

 

 

 

 

 

А (]) = Е (гг) е (гг — /)'} G-1,

j

>

0

 

 

(5.2)

(где G-1 — обобщенная обратная

матрица для

G),

имеем

 

А (j)G = E{ z

(n ) e ( n - j ) '},

 

 

 

 

 

ибо вектор е (п j) с вероятностью 1

ортогонален нулевому

под­

пространству матрицы G. Следовательно, А (/) е (п — /) есть проек­

ция X (п) на пространство, порождаемое е (п — /),

так

как

 

Е [{ж (7г) — А (/) s (п — /)} е (тг — ;)'] = 0.

Таким образом, мы можем построить процесс

со

со

 

u ( n ) = ^ A ( i ) e ( n — i),

SX(/)GM (/)* < °°.

(5.3)

о

о

 

Последнее соотношение вытекает из неравенства Бесселя. Полагая V (п) — X (п) и (тг), имеем следующую теорему:

Т е о р е м а 2" (разложение Вольда). Если х (?г) — стационар­ ный векторный процесс с р компонентами, то х (;г) — и {п) -1- г; (тг), где и (гг) определяется соотношениями (5.1), (5.2) и (5.3), и (п) — детерминированный процесс и Е (гг) е (ггг)') = 0.

Доказательство в основном такое же, как в скалярном случае. Если присутствует только компонента и (гг), то мы говорим, что

процесс X (гг) ч и с т о н е д е т е р м и н и р о в а н н ы й .

х) Этот параграф может быть опущен.

12*

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