Быстрый поиск

Учебное пособие: Анализ временных рядов

yt-1 < yt > yt+1 yt-1 > yt < yt+1

yt-1 yt yt+1 yt-1 yt yt+1

Рис. Поворотные точки.

Для определения поворотной точки требуются три последовательных значения. Начальное и конечное значения не могут быть поворотными точками, т. к. неизвестно y0 и yN+1. Если ряд случаен, то эти три значения могут следовать в любом из шести возможных порядков с равной вероятностью. Только в четырёх из них будет поворотная точка, а именно, когда наибольшее или наименьшее из трёх значений находится в середине. Следовательно, вероятность обнаружения поворотной точки в любой группе из трёх значений равна 2/3.

с с c c c c

b b b b b b

а а a a a a

Рис. Варианты взаимного расположения трёх точек.

Для группы из N величин определим счётную переменную Х.

ì 1, если yt-1 < yt > yt+1 или yt-1 > yt < yt+1

Х = í

î 0, в противном случае.

Тогда число поворотных точек р в ряде есть просто , а их математическое ожидание есть М[p]=2/3(N-2). Дисперсия числа поворотных точек вычисляется по формуле D[p]=(16N-29)/90, а само распределение близко к нормальному.

2. Критерий, основанный на определении длины фазы

Интервал между двумя поворотными точками называется фазой. Для того , чтобы установить наличие фазы длины d (например, восходящей) , нужно обнаружить d+3 членов, содержащих падение от первого члена ко второму ,затем последовательный подъем до (d+2)-го члена и падение к (d+3)-ему члену.

1 2 3 4 d+1 d+2 d+3 N

рис. 3. Фаза длины d.

Рассмотрим группу из d+3 чисел, расположенных в порядке возрастания. Если, не трогая двух крайних членов, извлечь пару чисел из оставшихся d+1 и одно из них поставить в начало, а другое в конец, получим фазу длины d. Существует способов такого выбора пары чисел и каждый член пары может быть поставлен в любой конец, следовательно число восходящих фаз равно d(d+1).

Кроме того, поворотные точки будут иметь место, если первый член последовательности поставить в конец, а любой из оставшихся, за исключением второго, поместить в начало. Число таких последовательностей составит (d+1) . Еще столько же последовательностей получиться если последний член в исходной, возрастающей, последовательности поставить в начало, а любой другой, кроме последнего, в конец. Во избежании двойного счета следует исключить случай, когда первый член ставится на последнее место, а последний на первое. Таким образом, в последовательности из (d+3) чисел с фазой длиной d число случаев роста составит

d(d+1)+2(d+1)-1 =+3d+1 .

Число возможных последовательностей из (d+3) чисел равняется числу перестановок (d+3)!, так что вероятность либо восходящей, либо нисходящей фазы равна

В ряде длины N последовательно можно выделить N-2-d групп по d+3 членов. Т.о. математическое ожидание числа фаз длины d

Можно показать, что математическое ожидание общего числа фаз длины от 1 до N-3

3.Критерий , основанный на знаках разностей

Данный критерий состоит в подсчете числа положительных разностей первого порядка в ряде , иначе говоря, числа точек возрастания ряда. Для ряда из N членов получаем N-1 разностей . Определим счетную переменную как

Если теперь обозначить через с число точек возрастания случайного ряда , то

Распределение довольно быстро стремится к нормальному с дисперсией

В основном данный критерий рекомендуется для проверки наличия линейного тренда . С другой стороны, критерий, основанный на поворотных точках , плохо подходит для обнаружения тренда, т.к. наложение заметных случайных колебаний на умеренный тренд приводит примерно к тому же множеству поворотных точек, что и при отсутствии тренда. .Более совершенным, но более сложным критерием для обнаружения линейного тренда являются регрессия y на t и проверка значимости регрессионного коэффициента.

4.Критерий, основанный на ранговых сравнениях

Идею сравнения соседних значений ряда можно развить до сравнения всех значений. Для данного ряда подсчитаем число случаев, когда очередной член ряда превышает все последующие. Всего для сравнения имеется N(N-1) пар. Пусть n общее число случаев превышения. Подсчитывают ранговый коэффициент корреляции Кендэла