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

книги / Статистический анализ временных рядов

..pdf
Скачиваний:
5
Добавлен:
19.11.2023
Размер:
47.02 Mб
Скачать

10 .2

ЭФФЕКТИВНОЕ ОЦЕНИВАНИЕ ФУНКЦИЙ ТРЕНДА

6 2 а

Условие 10.2.3. Предел при Т -+■ оо отношения

 

(61)

 

aif №

 

 

= rfj (h)

 

 

Va^WaJ^O)

 

существует для всех i, /,

h, принимающих значения i, ] = 1 ,

..., p

и h = 0 ,

± 1 , ... .

 

 

Пусть

 

 

 

(62) lim

rjj(h) = pij(h),

i, j = 1, ... ,/> , h = 0, ± 1 , . . . ,

 

T—►oo

И R (h) = [pij(h)].

Условие 10.2.4. Матрица R (0) не вырождена.

Первое предположение означает, что вклад каждой последова­ тельности г» неограниченно возрастает с ростом Т. Благодаря вто­ рому предположению исключается возможность заметного вкла­

да последних членов г?/ в сумму квадратов при больших Т. Третье предположение приводит к тому, что корреляции между независи­ мыми переменными при всех достаточно больших Т мало отличают­ ся от некоторых фиксированных значений. Наконец, четвертое ус­ ловие вводится из соображений удобства вычислений.

Последовательность {R (/г)} является положительно определен­ ной последовательностью эрмитовых матриц, т. е.

(63)

2! с'R (h— k) 1хьхк > 0 ,

 

h,k

где суммирование проводится по любому конечному множеству индексов, {xh} — произвольная последовательность действительных чисел, а с — произвольный комплексный вектор размерности р. Отсюда следует, что существует эрмитова матрица М (А,) с положи­ тельно полуопределенными приращениями, такая, что

Л

(64)

 

R(/t) =

j eMdM(%).

 

 

 

 

—Л

 

 

(См. § 7.3.)

 

 

 

 

 

Пусть

 

 

 

 

 

 

 

~ V атп (0)

0

. . .

0

(65)

Dr =

У' Й22 (6)

• • •

0

 

 

 

 

 

 

0

0

*• * У а рр(т0) _

*•€24 ЛИНЕЙНЫЕ ТРЕНДЫ И СЛУЧАЙНЫЕ СОСТАВЛЯЮЩИЕ Гл. 10.

Тогда

 

 

 

(66)

R(A) =

lim D f'A r (A) DF1,

где Ат (А) =

[а!/ (Л)1.

Т—►оо

 

 

 

Для того чтобы найти пределы величин

 

 

Г - 1

 

(67)

Df'Z'SZDF1=

D f1 s

a(A)Ar(A)DF‘

и

 

Л = —(Г — I)

 

 

 

 

(68)

D ^Z'Sr'ZD F1,

сделаем следующее предположение о последовательности ковари­ аций.

Условие 10.2.5. Имеет место представление

 

Я

<69)

а (А) = J eiMf(X)dX

 

—Я

с функцией f (Я), непрерывной на [—я, я1.

Указанные пределы найдем, аппроксимируя спектральную плот­ ность /(А) тригонометрическими полиномами /у(A) = fL(—А.) и /у (А) = / у ( —А), такими, что

*<70)

/у (А )< /(А )< /у (А), — я < А < я.

Поскольку 0 < / (А) < /у(А), то /у (А) является спектральной плот­ ностью. Ей соответствует ковариационная последовательность

Я

Оу (А) = £ еШ| fu(X) dX. Если 0 < /L(А), то и /у(А) будет спектраль-

—Я

ной плотностью, порождающей ковариационную последовательность

Я

d (А) = \ е,ХАfL (A) dX. Однако в ряде доказательств неотри-

—Я

дательность /у(А) не является необходимой.

Лемма 10.2.6.

Обозначим

2у = [оу(*— s)], 2у = [ay (/ — s)].

.Если выполнено соотношение (70), то для любого вектора х

г{71)

х'2ух <

х'2х < х'2ух.

Если, кроме того, fL (А) >

0,

то для любого вектора х

;(72)

х '2 у 'х

< х '2 ~ 'х < х ' Е Г ‘ х.

Доказательство. Прежде всего

 

 

г

(73)

х '2 х =

S

о (/ — s) xpct =

( , S = 1

10.2 ЭФФЕКТИВНОЕ ОЦЕНИВАНИЕ ФУНКЦИЙ ТРЕНДА 625

= 1

/.S-»

2

niXt f (X) dX.

Аналогичные выражения имеют место для х' 2^ х и х' Syx. (Если II (X) при некотором значении X отрицательно, то oL(u) уже не бу­ дет ковариационной последовательностью. Однако и в этом случае соотношения (71) и (74) выполняются.) Из (70) вытекает, что

{74) 2 / е (А) £&< ]Г 2 f(X)dX<

—Я

2

fu(X) dX.

Последнее равносильно (71). Что касается (72), то справедливость неравенства х'2х < х' Syx для всех х эквивалентна тому, что кор­

ни уравнения | 2 — 0 2 У| = 0 положительны и

не превосходят 1.

Они же будут являться и корнями уравнения |

2£/‘— 02-1 1 = 0,

а отсюда следует справедливость первого неравенства в (72) для всех х. Аналогично доказывается и второе неравенство, в

Теорема 10.2.6. При условиях 10.2.1,

10.2.2, 10.2.3 и 10.2.5

 

 

 

Я

 

(75)

lim Df'Z'SZQF1= 2л

f

/ (X) dtA (X).

 

Т_1____

4

**

 

Доказательство. Идея состоит в том, чтобы доказать, что левая

оо

часть (75) равна 2 а (Л) R (Л)- Последняя же сумма может быть

Л=—ОО

я

оо

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

fa (h) ethh с1ЩХ), а этой есть правая часть

Л = — оо

J

 

—Я

(75). Однако имеется некоторое затруднение, связанное с тем, что

число слагаемых в сумме (67) неограниченно возрастает,

и с тем,

что мы исходим из непрерывности f (X), а не из условия

OQ

Iа W I <

2

/Г-

—оо

< оо. [Предположение о непрерывности f (X) было сделано по той причине, что оно используется и при доказательстве теоремы 10.2.7.]

По теореме 7.5.2 для любого е >

0 найдутся такие два тригоно-

к

к

метрических полинома, /F (X) = 2

и fu(ty = 2

fc=-K

fc=-K

626 ЛИНЕЙНЫЕ ТРЕНДЫ И СЛУЧАЙНЫЕ СОСТАВЛЯЮЩИЕ Гл. 10.

=

CL-k и bk ~ b - k, k = 1, .... К, что выполняется (70) и

 

(76)

fu(l) — /ь(Я)<е,

— я < Я < я .

 

Далее,

 

 

Т

 

 

 

 

(77)

 

 

 

 

 

 

-

S /=12 xteiKi

f u(k) dX =

 

 

 

 

=

к

S

*<+**<•

 

 

2

 

 

 

 

k*=—K

t£Sk

— A, .... T},

k < 0.

где Sft

= {1, .... 7 - A },

£ > 0 ,

и

={1

Для произвольного вектора у размерности

p положим я* =

(zy, ...

.... г„<) D f‘ v. Тогда

 

 

к

 

 

 

(78)

 

 

 

 

 

 

Y'DF Z'SZDf'v <

2 ^ Y 'D F 'A r^ D F ’Y-

 

Поэтому

 

k——K

 

 

 

 

 

 

К

 

 

(79)

___

 

 

 

 

 

lim Y'DF’Z'SZDF'Y <

2 bky'R (k) у =

 

 

 

 

=

2лу'

J ]

j eiMdm(%) у =

 

 

 

 

 

k~ —K

 

 

 

 

 

 

 

л

 

 

 

 

 

=

2яу'

\fu{K)dfA(X)y.

 

Подобным же образом получаем, что

—71

 

 

 

 

 

 

Л

 

 

 

 

 

 

(80)

2яу' J fL(Я) dm (Я) у < limjy'Df'Z 'SZDFV

 

 

—Я

 

 

Г-coo

 

 

Поскольку е произвольно, из (79) и (80) вытекает

 

 

lim Y'DF'Z'SZDF'Y = 2яу'

Л

 

 

(81)

( / (X) dm (Я) у.

 

 

 

 

 

 

 

 

Так как последнее соотношение выполняется для любого вектора у. отсюда и следует утверждение теоремы, щ

Условие 10.2.6. Функция / (Я), определяемая в (69), удовлетворяет соотношению f (Я) > 0, —я < Я < я.

Теорема 10.2.7. При условиях 10.2.1,10.2.2, 10.2.3,10.2.5 и 10.2.6

Л

(82)

10 .2

ЭФФЕКТИВНОЕ ОЦЕНИВАНИЕ ФУНКЦИЙ ТРЕНДА

627

Доказательство. В силу следствия 7.5.2 для любого е > 0 най­ дется такая пара спектральных плотностей, а именно, Д. (Я) =

к

к

= [2я 12 аь

|2р ' иДг (Я,)= [2я | 2 bk e‘XkI г]-1> соответствующих

й=0

k—Q

некоторой паре процессов авторегрессии, что при этом выполняются соотношения (70) и

(83) 7 I W ^ 1 ^ W < 8 ’ ~ л < к < п-

Функция fu (Я) является спектральной плотностью процесса {до,}, удовлетворяющего стохастическому разностному уравнению

(84)

2 bkwt- k = vt,

* = . . .

, — 1, 0, 1, . . . ,

 

fe=0

 

 

в котором {о,} — последовательность некоррелированных случай­ ных величин с нулевыми средними и единичными дисперсиями.

Предположим, что Т >■ 2К. Элементы Оти матрицы Sy1 равны

 

 

 

 

min(f,s)—

 

 

 

 

 

 

 

(85)

 

 

оти<3 и-—

2

bjcbk+\t—sl

 

 

 

 

 

 

 

fc=0

 

 

 

 

 

 

 

 

 

 

=

о Тт Г - * 'Т+1- ° ,

 

* , / -

1 .............. к .

Для К <. шах (t, s) и min (t, s)

T

К они равны

 

(86)

 

Оти ~

2

bkbk-\-\t-s\t

0 < | f —■s\

К,

 

fc=.0

 

 

 

 

 

 

 

 

 

 

 

 

0,

 

 

 

 

 

 

K < \ t - S \ .

(См. § 6.2.) Пусть

0

 

 

 

 

0

 

0 . ..

 

0 ”

 

 

~Ьп

0 .. .

 

0

 

0

 

 

bn

bn

0

..,.

 

0

0

 

0 . ..

0

0

 

 

bn

bn

bn

• •

 

0

0

 

0 . ..

0

0

(87)

В

Ьк\

Ьк2

Ькз

•.

ькк 0

 

0 . ..

0

0

 

 

Ьк

bK-i

 

,,

 

ьх

ъй

0 . ..

0

0

 

 

Ьк-2 • ■

 

 

 

 

0

 

 

0

Ьк

Ьк- 1 ..

 

Ь2

ьх Ьо ..

0

 

 

6

0

0

.... 6 6 6 . • * ь0 0

 

 

0

0

0 . ..

 

0 0 0 . .. ъх Ьц _

628

ЛИНЕЙНЫЕ ТРЕНДЫ И СЛУЧАЙНЫЕ СОСТАВЛЯЮЩИЕ

Гл. 10,

где bn , bn , bi2, ..., Ькк выбраны так, что

 

(88)

V = В'В.

 

[Условием (88) треугольная матрица В определяется однозначно. Последние Т К строк правой части (87) непосредственно восста­ навливаются по имеющимся представлениям (85) и (86).J При этом

(89)х'2у1х = х'В'Вх =

 

 

 

К

/ i

 

\ 2

 

Т

/

К

 

\ 2

 

 

 

= 2( S b ijX j) -f-

2

( 2 b jX t- A

 

 

 

<-1

\ / = l

 

/

teK +1 \/= o

 

/

\2

 

 

 

К

т

 

 

 

 

К

/ {

 

 

 

=

2

2

 

bjb^xt-ixt-k +

2

2

Ьих,- .

 

 

 

/.ft-0 <-K+l

 

 

 

 

«=1\/=l

 

/

Для произвольного вектора Y размерности

р

положим х,

= (zit, ...

-.Zp/JDr1V- Тогда

 

 

 

 

 

 

 

 

 

 

 

 

(90)

Y'DF'Z'S^ZDF’Y >

Y'DF’Z 'V Z DF'Y =

х'В'Вх >

 

 

 

 

 

к

 

т

bjbkx,-jXt-k

 

 

 

 

 

 

>

2

2

 

 

 

 

и

 

 

 

j,k—0 Ы К + 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(91)

lim -y'D r^'S

*ZD

 

 

 

 

 

 

 

 

 

Т^оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

т

 

 

 

 

 

 

 

 

 

 

bfik

 

 

 

 

 

 

 

 

>

Иm ' 2

2

 

v 'D ^ z ^ /Z ^ D ? 1?

 

 

 

^

Г-иэо /.Д -0

 

*=/С+1

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

=

Iim

2

bjbk4 ' DTXk T (k — / ) D r 1Y =

 

 

 

 

T-+OQ j,k=0

 

 

 

 

 

 

 

 

 

 

=2 ^ A v'R (^ — /)v =

j,k«*0

кn

=Y' 2 bjb* f eiUk~n dm (X) у =

«Л/ dM (ty Y =

- W l i *

i r v ' 5

—Я

10.2

ЭФФЕКТИВНОЕ ОЦЕНИВАНИЕ ФУНКЦИЙ ТРЕНДА

62»

Аналогично получаем

 

Я

 

 

 

 

 

 

 

(92)

Шп v'DF'Z'S-'ZDF'v < ~

у'

j

М V-

 

Поскольку е выбирается произвольно, то из (91) и (92) вытекает,

что

 

 

 

я

 

 

(93)

ton Y'D F 'Z 'S ^Z D ^ Y = ^

y

f

(X) у.

 

 

 

 

—Я

 

 

Это соотношение выполняется для каждого вектора у, что и доказы­ вает теорему.

Следствие 10.2.1. При выполнении условий 10.2.1—10.2.5

 

Иш DriS (b*Р) (b* — р)'Dr «=2nR_1 (0)

Я

(94)

ff (X) dM (X) R_1 (0);

 

Т-юо

_jt

при выполнении условий 10.2.1—10.2.6

 

(95)

lim Dr S (b — p) (b — P)' Dr =

 

 

T-+oo

 

Здесь использован тот факт, что lim D f1 Z'ZDf1= R (0).

T-*о©

Пусть

(96)S(u) = {Х|2я/(Х)<и},

(97)

T (u) = j dM (X),

 

S(u)

где m — 2я min / (X) и M — 2n max f (X), —n <. X < я. Поскольку / (X) = f (—X), то множество S (и) симметрично относительно нуля,

так

что матрицы Т («) и Т (и2) — Т (%), т < мх < ы2 <

М, дей­

ствительны и положительно полуопределены. При этом

 

 

 

 

м

 

(98)

 

lim D f’Z'ZDf1== f dT («),

 

 

 

 

M

 

(99)

 

lim DF'Z'SZDF1=

f udT («),

 

 

 

Г -к »

OT

 

(100)

lim D F'Z'S-’ZDF1=

f — dT (a).

 

 

 

Г-*>оо

Jmw

 

Существует

такая невырожденная

матрица P, что P'R(0) Р —

*=I,

Р'* lim DF* Z'SZDF* P = D, где

D — диагональная

матрица

 

T-*oo

 

 

 

630

ЛИНЕЙНЫЕ ТРЕНДЫ И СЛУЧАЙНЫЕ СОСТАВЛЯЮЩИЕ

Гл. 10.

с

элементами dn > ... > dpp. Положим

L {и) = Р'Т (и) Р,

т <

<

а <

Af. Тогда матричная функция L (и) и ее приращения поло­

жительно полуопределены. Кроме того,

м

 

 

 

 

 

<101)

1 = Р' lim D7lZ'ZDF‘P =

[dh (и),

 

 

 

г-*.»

щ

 

 

 

 

 

м

 

<102)

D = Р' lim DF'Z'SZDf'P =

f udL (и),

 

 

 

T-+°°

 

m

 

 

 

 

М

 

 

<103)

Р' lim DF'Z'S-'ZDF'P =

f — dL (и).

 

 

 

T-+OQ

J

U

 

 

 

 

m

 

 

Две интересующие нас предельные нормированные ковариационные матрицы будут совпадать тогда и только тогда, когда (103) равно

D -1. Диагональные элементы матриц (101), (102) и (103) равны

м

<Ю4)

1 = j d/t4. («),

t =

l .......... р,

 

ТП

 

 

 

М

 

 

<105)

dH = J udlu (и),

i =

1, . . . . р,

 

m

 

 

 

М

 

 

<106)

^ ± d l (((U),

i = i ..........р.

 

m

 

 

В силу того, что приращения функций (ы) неотрицательны, 1ц (и) для каждого i является функцией распределения вероятностей. Выражение (105) можно рассматривать как математическое ожидание случайной величины, имеющей такое распределение, а (106) —. как математическое ожидание обратной ей по значениям случайной

величины. Если (103) есть D-1, то (106) равно l/du, и из леммы 10.2.1 вытекает, что для каждого i функция 1ц (и) имеет только одну точку роста, причем в этой точке она имеет скачок величины 1.

Обратно, пусть каждая из функций /« (и) имеет только одну точку роста. В соответствии с (104), в этой точке будет скачок ве­ личины 1. Пусть скачки диагональных элементов матричной функ­ ции L («) соответствуют значениям и, среди которых различными являются «х > ... > UQ. Поскольку приращения L (и) положительяо полуопределены, возрастание недиагональных элементов может происходить только путем скачков при тех же значениях и — ul t ...

UQ. Пусть L/ обозначает приращение L (и) в точке и щ, / =

10 .2

ЭФФЕКТИВНОЕ ОЦЕНИВАНИЕ ФУНКЦИЯ ТРЕНДА

6 3 t

r= 1,

G. Тогда Ц01) принимает вид

 

(107)

/ i L/ = i-

 

Матрицы L/ соответствуют отличным от нуля матрицам С<Л) из дока­ зательства теоремы 10.2.2. При этом /-й диагональный блок матрицы L является единичной матрицей, порядок которой равен числу диагональных элементов L (и), возрастающих в точке щ. Все дру­ гие блоки L/ являются нулевыми матрицами. Ранг матрицы L/ совпадает с порядком единичного блока. Соотношения (101) и (102) принимают вид

(Ю8) 2 U/Lj = D, /=1

в

Теорема 10.2.8. Для асимптотической эффективности оценок наименьших квадратов при условиях 10.2.1—10.2.6 необходимо и достаточно, чтобы каждый диагональный элемент матричной функции L (и) возрастал только в одной точке.

Отметим, что число диагональных элементов матрицы L/, рав­ ных 1, совпадает с рангом L/.

Следствие 10.2.2. Для асимптотической эффективности оценок наименьших квадратов необходимо и достаточно, чтобы матрич­ ная функция L (и) имела G, G <■ р, точек роста, где р — сумма рангов приращений функции L (и).

Множества

S/

= (Я. 12я/ (Я) =

щ, —я < Я < я}, / = 1,

..., G,

являются множествами значений Я, на которых возрастает

М (Я),

Поэтому

 

 

 

 

(ПО)

Р'

У^М(Я)Р = Ь/(

/ = 1 ..........О.

 

 

 

si

 

 

Теорема 10.2.9. Для асимптотической эффективности оценок наименьших квадратов при условиях 10.2.1—10.2.6 необходимо и. достаточно, чтобы на множестве тех Я, при которых возрастает М (Я), спектральная плотность / (Я) принимала всего G, G < р, значений, где р равно сумме рангов интегралов JdM (Я) по тем G множествам значений %, на которых f (Я) принимает эти G значе­ ний.

Равные единице диагональные элементы матрицы Lj должны располагаться рядом, а диагональные блоки, порождаемые ими,,

'632 ЛИНЕЙНЫЕ ТРЕНДЫ И СЛУЧАЙНЫЕ СОСТАВЛЯЮЩИЕ Гл. 10.

должны составлять единичные подматрицы. Учитывая это, соотно­ шение (108) можно переписать в виде

 

 

u j

0

.

О

( 111)

D =

0

ы21

.

О

 

 

 

 

 

 

 

 

 

О0 ... UQ \J

где матрица и/1 имеет порядок, равный рангу L/. В силу того, что единичные матрицы в Lx, ..., L<? не перекрываются,

(112)

 

L,L/ = 0,

1 ф и

Множества 5/, /

=

1, ... G, значений Я, при которых возрастает

М (Я), называются

элементами спектра М(X). Объединение этих

множеств называют

спектром

М(Я). Поскольку множество S/

симметрично, интеграл (ПО) содержит только действительную часть М (Я). Фактически соотношение (110) имеет вид

(113)

2Р' j

йШ(К)Р = 1/,

0 g S y,

<114)

Р'Л [ М ( 0 + ) - М ( 0 Н 1 Р + 2Р'

J йЛш (Я)Р = L/,

 

 

 

5 ; П (0 ,Л ]

 

 

 

 

 

 

о e s , .

Отметим, что JdM (Я,)=

R (0),

где интегрирование

производится по

объединению всех G указанных множеств значений Я, и что последо­

вательность G матриц

JdM (Я) является

ортогональной в метрике

R (0)

в силу (112).

Sy

 

 

 

Теорема 10.2.10. При условиях 10.2.1—10,2.4 оценки наименьших квадратов асимптотически эффективны для всех стационарных процессов с непрерывными положительными спектральными плот­ ностями тогда и только тогда, когда М (Я) возрастает не более, чем при р значениях Я, 0 < Я я, где р сумма рангов прираще­ ний функции М (Я).

Доказательство. Теорема 10.2.9 может оставаться справедливой для всех указанных спектральных плотностей только в том случае, если матрица JrfM (Я) положительно полуопределена не более чем в р точках.и

Следует отметить, что если действительная часть М (Я) не воз­ растает на некотором множестве значений Я, то на этом множестве не возрастает также и мнимая часть М (Я). Поэтому указанный критерий можно формулировать как с использованием самой М (Я), так и с использованием лишь ее действительной части.