Диплом / suslov_ibragimov_ekonometrika
.pdf414 |
Глава13.Спектральный и гармонический анализ |
Исследованияпоказывают,чтоналичиенепериодическоготренда(трендасбесконечным периодом)дает скачок на нулевой частоте,т.е.в начале координат спектральной функции.При наличии циклических составляющих в соответствующих частотах имеется всплеск;если ряд слишкомÇзазубренÈ,мощность спектра перемещается в высокие частоты.
Типичнымдлябольшинстваэкономическихпроцессовявляетсяубываниеспектральной плотности по мере того,как возрастает частота.
Процесс выделения существенных гармоникÑитеративный.При изучении периодограммы выделяется две-три гармоники с максимальной интенсивностью. Находятся оценкипараметров этих наиболее существенных гармоник,и они удаляются из временного ряда с соответствующими весами.Затем остатки временного ряда,получающиеся после исключения значимых гармоник,снова изучаются в той жепоследовательности,т.е.строится периодограмма для этих остатков,ипроявляются те гармоники,которые на начальном этапе былинезаметны,и т.д.Количество итераций определяется задаваемой точностью аппроксимации модели процесса, которая представляется в виде линейной комбинации основных гармоник.
Понятие спектра,являясь основополагающим в спектральном анализе,для экономистов играет важную роль еще и потому,что существует функциональная связь выборочного спектра и оценок автоковариационной функции.
13.4.Связь выборочного спектра с автоковариационной функцией
Покажем,чтовыборочныйспектрпредставляетсобойкосинус-преобразование Фурье выборочной автоковариационной функции.
Теорема ВинераÑХинчина:
|
|
|
|
0 |
|
|
T −1 |
1 |
|
|
0 |
|
T −1 |
1 |
|
||||
|
|
|
|
|
|
k! |
|
|
|
! |
|
||||||||
|
|
p (f ) = 2 c0 + 2 |
ck cos 2πf k |
= 2s2 |
|
1 + 2 |
rk cos 2πf k , |
(13.25) |
|||||||||||
|
|
|
|
|
|
|
=1 |
|
|
|
|
|
|
k=1 |
|
|
|
|
|
где r |
k |
= c /c |
= c /s2 |
Ñвыборочные автокорреляции. |
|
|
|
|
|
||||||||||
|
k |
0 |
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Доказательство. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Объединим коэффициенты Фурье αf , βf |
в комплексное число df = αf − iβf , где |
|||||||||||||||||
|
i Ñмнимая единица.Тогда |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
T |
0αf2 + βf21 = |
T |
|
|
|
|
|
T |
|
|
|
|||
|
|
|
p (f ) = |
|
|
(αf − iβf ) (αf + iβf ) = |
|
df df , |
(13.26) |
||||||||||
|
|
|
2 |
2 |
2 |
где df комплексно сопряжено с df .
416 |
Глава13.Спектральный и гармонический анализ |
Теперь допустим,что выборочный спек тр,характеризующийся эмпирическими значениями частоты,амплитуды и фазы,вычислен для ряда из T наблюдений и мы можем неоднократно повторить этот эксперимент,соответственно собрав множество значений αj , βj и p (f ) по повторным реализациям.Тогда среднее значение p (f ) будет равно:
0 |
T −1 |
1 |
|
k! |
|
||
E(p (f )) = 2 |
E(c0) + 2 |
E(ck ) cos 2πf k , |
(13.28) |
|
=1 |
|
|
где c0 и ck Ñэмпирические значения автоковариации. Cучетом того,что |
E(ck ) |
при больших T стремится к теоретической автоковариации γk ,получим,переходя
к пределу,так называемую теоретическую |
спектральную плотность,или спектр |
||||||||||
мощности: |
|
|
|
|
|
|
|
|
|
|
|
p(f ) = lim E(p (f )), 0 |
|
f |
|
1 |
|
|
|||||
|
|
|
|
|
|
||||||
! |
! 2 |
|
|
||||||||
|
T →∞ |
|
|
|
|
|
|
||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
∞ |
|
|
|
1 |
|
|
|
|
|
|
k! |
|
|
|
|
|
|
|
||
p(f ) = 2 |
γ0 + 2 |
γk cos(2πf k ) |
|
= |
|
|
|||||
|
|
|
=1 |
∞ |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
1 |
|
||||
|
|
k! |
|
|
|
|
|||||
|
= 2σ2 |
|
1 + 2 |
ρk cos(2πf k ) , |
(13.29) |
||||||
|
|
|
|
=1 |
|
|
|
|
|
|
|
где ρk = γk /γ0 = γk /σ2 Ñтеоретические автокорреляции.
Итак,это соотношение связывает функцию спектральной плотности с теоретическими автоковариациями.
Иногда более удобно использовать автокорреляции:разделим обе части p(f ) на дисперсию процесса, σ2,и получим нормированный спектр g(f ):
|
∞ |
|
1 |
|
||
0 |
k! |
1 |
|
|
|
|
! f ! 2 . |
||||||
g(f ) = 2 1 + 2 |
=1 |
ρk cos 2πf k , 0 |
||||
|
|
|
|
|
Если процесс представляет собой белый шум,то,согласно приведенным формулам, p(f ) = 2σ2,т.е.спектральная плотность белого шума постоянна.Этим объясняется терминÇбелый шумÈ.Подобно тому,как в белом цвете смешаны в одинаковых объемах все цвета,так и белый шум содержит все частоты,и ни одна из них не выделяется.
13.5.Оценка функции спектральной плотности |
417 |
13.5.Оценка функции спектральной плотности
На первый взгляд,выборочный спектр,определенный как
|
|
T −1 |
|
|
|
|
|
p (f ) = |
2 |
! ck e−i2πf k = |
|
|
|
(13.30) |
|
|
k=−(T −1) |
|
|
|
|
|
|
|
0 |
T −1 |
1 |
|
1 |
|
|
|
k! |
! f ! |
|
|
|
||
= |
2 |
c0 + 2 |
ck cos 2πf k , 0 |
2 |
, |
(13.31) |
|
|
|
=1 |
|
|
|
|
|
является естественной и правильной оценкой функции спектральной плотности:
p(f ) = 2 |
+∞ |
γk e−i2πf k = 2 γ0 + 2 |
∞ |
γk cos 2πf k , 0 ! f ! |
1 |
. (13.32) |
|
! |
k! |
|
|||||
|
0 |
1 |
|
|
|||
|
|
−∞ |
|
|
|
||
|
k= |
|
|
=1 |
2 |
|
|
|
|
|
|
|
|
Известно,что выборочная автоковариация ck Ñэто асимптотически несмещенная оценка параметра γk ,так как
lim E(ck ) = γk ,
T →∞
и поэтому выборочный спектр есть также асимптотически несмещенная оценка функции спектральной плотности.
Однакодисперсияоценки(выборочногоспектра)неуменьшаетсяпомерероста размера выборки.Это означает,что рассмат риваемая оценка несостоятельна.В то время как график функции теоретической спектральной плотности стационарного стохастического процессаÇгладкийÈ, Ñ график выборочного спектра,построенный на основе эмпирических данных, ÇнеровныйÈ.Использование данной оценки в качестве оценки функции спектральной плотности может привести к ложным выводам.
Поведение выборочного спектраиллюстрируютспектрограммы нарис. 13.1а), б),в).Гладкая жирная линия соответствует теоретической спектральной плотности случайного процесса,а неровная тонкая линияÑоценке по формуле(13.24). Видно,что с увеличением длины ряда оценка не становится более точной,а только увеличивает частоту флуктуаций.
Существует два подхода к решению проблемы несостоятельности выборочного спектра как оценки теоретического спектра.Оба заключаются в том,что выборочный спектр сглаживается,так что за счет некоторого увеличения смещения этой оценки достигается существенное снижение дисперсии(см.рис. 13.1г) ).Один подход сглаживает оценку спектра вÇчастотной областиÈ,видоизменяя формулу (13.24),другой же подход видоизменяет формулу(13.25).
418 |
Глава13.Спектральный и гармонический анализ |
|
|
$) |
-) |
|
T = 100 |
T = 300 |
|
,) |
") |
|
T = 1000 |
T = 1000, |
|
|
!"#$%&''$( |
|
|
)*&'+$ |
Рис. 13.1
1)Взвешивание ординат периодограммы.
Первый способ сглаживания выборочного спектра применяет метод скользящего среднего,рассмотренный в предыду щей главе,к значениям периодограммы. Сглаживающая оценка определяется в форме
|
#T $ |
|
k= M |
k |
|
# |
T |
$ |
k= M |
k j −k |
|
||
|
|
|
|
|
M |
|
|
|
|
|
M |
|
|
ps |
|
j |
= |
!− |
|
p |
|
j − k |
= |
!− |
|
(13.33) |
|
|
|
|
µ |
|
|
|
µ I . |
||||||
Здесь {µ−M , . . . , |
µ−1, µ0, µ1, . . . , µM } Ñ 2M +1 коэффициентов скользящего |
среднего,которые в сумме должны давать единицу,а также должны быть симметричными,в смысле µ−k = µk .Как правило,коэффициент µ0 максимальный, а остальные коэффициенты снижаются в обе стороны.Таким образом,наибольший вес в этой оценке имеют значения выборочного спектра,ближайшие к данной
частоте Tj ,в связи с чем данный набор коэффициентов называют спектральным окном.Через это окно мы как быÇсмотримÈна функцию спектральной плотности.
Формула сглаженной спектральной оценки определяется только для значений j = M, . . . , [T /2] − M .Для гармонических частот с номерами j = 0, . . . , M − 1 и j = [T /2] − M + 1, . . . , [T /2] оценка не определена,поскольку при данных значениях j величина( j − k)может принимать значения −M, . . . , −1; [T /2] + 1, . . . , [T /2] + M .Однако проблема краевых эффектов легко решается,
13.5.Оценка функции спектральной плотности |
419 |
если доопределить функцию выборочного спектра(и,соответственно,периодограмму),сделав ее периодической.Формула(13.24)дает такую возможность,поскольку основана на синусоидах и косинусоидах.Таким образом,будем считать,что выборочный спектр определен формулой(13.24)для всех частот f (−∞, +∞). При этом ясно,что выборочный спектр будет периодической функцией с периодом 1,зеркально-симметричной относительно нуля(и относительно любой частоты вида i/2,где i Ñцелое число):
p (f − i) = p (f ), i = . . . , −1, 0, 1, . . .
и
p (−f ) = p (f ).
Формулу(13.33)несложнообобщитьтак,чтобыможнобылопроизводитьсглаживание не только в гармонических частотах.Кроме того,в качестве расстояния междуÇусредняемымиÈчастотами можно брать не 1/T ,а произвольное > 0. Таким образом приходим к следующей более общей формуле:
k |
M |
|
|
!− |
|
|
|
ps (f ) = |
µk p (f − k |
). |
(13.34) |
= M
Просматривая поочередно значения выборочного спектра и придавая наибольший вес текущему значению,можно сгладить спектр в каждой интересующей нас точке.
2)Взвешивание автоковариационных функций.
Известно,что при увеличении лага k выборочные автоковариации ck являются все более неточными оценками.Это связано с тем,что k-я автоковариация вычисляетсякаксреднеепо T −k наблюдениям.Отсюдавозникаетидеявформуле (13.25)ослабить влияние дальних авток овариаций за счет применения понижающих весов так,чтобы сростом k происходило уменьшение весового коэффициента при ck .
Если ряд весов,связа нных с автоковариациями c0, c1, . . . , cT −1,обозначить как m0, m1, . . . , mT −1,оценка спектра будет иметь вид:
|
0 |
T −1 |
1 |
1 |
|
|
|
|
k! |
где 0 ! f ! |
|
|
|
||
pc(f ) = 2 |
m0c0 |
+ 2 |
mk ck cos 2πf k , |
2 |
. |
(13.35) |
|
|
|
=1 |
|
|
|
|
|
Набор весов mk ,используемый для получения та кой оценки,получил название
корреляционное,или лаговое окно.
420 Глава13.Спектральный и гармонический анализ
Прииспользованиикорреляционногоокнадляуменьшенияобъемавычислений
удобно брать такую систему весов,что |
mk |
= 0 при k " K ,где |
K < T .Тогда |
||
формула(13.35)приобретает вид |
|
|
|
|
|
0 |
|
K −1 |
1 |
|
|
|
k! |
|
|
||
pc (f ) = 2 m0c0 |
+ 2 |
|
|
mk ck cos 2πf k . |
(13.36) |
|
|
=1 |
|
|
|
Корреляционное окно удобно задавать с помощью функции m(á),заданной на интервале [0; 1],такой,что m(0) = 1, m(1) = 0.Обычно функцию выбирают так, чтобы между нулем и единицей эта функция плавно убывала.Тогда понижающие веса mk при k = 0, . . . , K вычисляются по формуле
mk = m(k/K ).
Ясно,что при этом m0 = 1 (это величина,с помощью которой мы взвешиваем дисперсию в формуле(13.35))и mK = 0.
Наиболее распространенные корреляционные окна,удовлетворяющие перечисленным свойствамÑэто окна Парзе на и ТьюкиÑХэннинга(см.рис. 13.2).
Окно ТьюкиÑХэннинга: |
|
|
|
|
|
|
|
|
|
|
|
|
m(z) = |
|
1 |
(1 + cos(πz)) . |
|
|
|
|
|||
|
2 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
||
Окно Парзена: |
|
|
|
|
|
|
|
|
|
|
|
m(z) = |
|
1 − 6z2 + 6z3, |
если z [0; 21 ]; |
||||||||
|
|
2(1 |
|
z) |
3 |
если z |
|
[ |
1 |
; 1]. |
|
|
|
− |
, |
|
2 |
||||||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
Можно доказать,что сглаживание в частотной области эквивалентно понижающему взвешиванию автоковариаций.Ч тобы это сделать,подставим в(13.34) выборочный спектр,выраженный через комплексные экспоненты(13.30):
ps(f ) = 2 !M µj T!−1 ck e−i2π(f −j )k . j =−M k=−(T −1)
Меняя здесь порядок суммирования,получим
ps (f ) = 2 T!−1 ck e−i2πf k !M µj ei2πj k . k=−(T −1) j =−M
13.5.Оценка функции спектральной плотности |
|
421 |
|||
|
m(z) |
|
|
|
|
1 |
|
|
|
|
|
|
|
окно Парзена |
|
|
|
0.8 |
|
|
|
|
|
0.6 |
|
|
окно ТьюкиÑХэннинга |
|
|
|
|
|
|
|
|
0.4 |
|
|
|
|
|
0.2 |
|
|
|
|
|
|
|
|
|
|
z |
0 |
0.2 |
0.4 |
0.6 |
0.8 |
1 |
|
Рис. 13.2.Наиболее популярные корреляционные окна |
|
Введя обозначение
|
M |
|
|
M |
(ei2πj k |
+ e−i2πj k ) |
|
!− |
|
|
! |
||
mk = |
j = M |
µj ei2πj k = µ0 + j =1 µj |
||||
|
|
|
M |
|
|
|
|
|
|
j! |
|
|
|
= |
µ0 + 2 |
µj cos(2πjk |
), |
|
|
|
|
|
|
=1 |
|
|
|
придем к формуле
ps (f ) = 2 T!−1 mk ck e−i2πf k = k=−(T −1)
0 |
T −1 |
1 |
|
k! |
|
||
= 2 m0c0 |
+ 2 |
mk ck cos 2πf k , |
|
|
=1 |
|
|
которая,очевидно,совпадает с(13.35). |
|
|
|
Кроме того,без доказательства отметим,что подбором шага |
и весов µj мы |
всегдаможемсымитироватьприменениекорреляционногоокна(13.35)спомощью (13.34).Существует бесконечно много способов это сделать.
В частности,как несложно проверить,частотное окно ТьюкиÑХэннинга полу- |
|||||||||||
чим,если возьмем M = 1, µ |
0 |
= 1 |
, µ |
1 |
= µ |
−1 |
= 1 |
и |
= |
1 |
. |
|
|||||||||||
|
2 |
|
|
4 |
|
|
2K |