![](/user_photo/_userpic.png)
книги / Статистический анализ временных рядов
..pdf10 .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) _ |
![](/html/65386/197/html_fXvLgHGOl5.4JhA/htmlconvd-SYthD1625x1.jpg)
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 . .. ъх Ьц _ |
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ц (ы) неотрицательны, 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 (Я) положительно полуопределена не более чем в р точках.и
Следует отметить, что если действительная часть М (Я) не воз растает на некотором множестве значений Я, то на этом множестве не возрастает также и мнимая часть М (Я). Поэтому указанный критерий можно формулировать как с использованием самой М (Я), так и с использованием лишь ее действительной части.