Рассмотрим действительный стационарный в широком смысле случайный процесс
,
, с математическим ожиданием
,
, взаимной ковариационной функцией
, и взаимной спектральной плотностью
.
Предположим, имеются Т последовательных, полученных через равные промежутки времени наблюдений
за составляющей
, рассматриваемого процесса
. Как оценку взаимной спектральной плотности в точке
рассмотрим статистику
(2.1)
где
, - произвольная, не зависящая от наблюдений четная целочисленная функция,
для
, а
(2.2)
s – целое число,
- целая часть числа
.
Статистика
, называемая выборочной взаимной спектральной плотностью или периодограммой, задается соотношением
(2.3)
определено равенством (2.2).
Известно, если
рассматривать как оценку взаимной спектральной плотности
в точке
, то она является асимптотически несмещенной, но не состоятельной оценкой этой спектральной плотности. Заметим, что оценка (2.1) взаимной спектральной плотности
построена путем осреднения значений периодограммы в точках
некоторой весовой функцией
.
Лемма 3. Для любого действительного
, и любого
справедливо неравенство

где
- ядро Фейера, задаваемое равенством
(2.4)
, а
, (2.5)
Доказательство. Учитывая чётность функции
и элементарное неравенство
(2.6)
справедливое для всех x, таких, что
, имеем

Сделаем замену переменной интегрирования
тогда правая часть последнего неравенства примет вид

Применив для оценки первого интеграла, стоящего в квадратных скобках, неравенство
, а для оценки второго – неравенство
, получим

Лемма доказана.
Проведен численный анализ для соотношения (2.5) при Т=100 и при
, T
, где T - число наблюдений и получены следующие результаты
|
|
|
| 0,1 | 0.663138 | 2.13239 |
| 0,2 | 0.447986 | 1.48005 |
| 0,3 | 0.308154 | 1.04694 |
| 0,4 | 0.216092 | 0.7554 |
| 0,5 | 0.154768 | 0.556644 |
| 0,6 | 0.113483 | 0.41954 |
| 0,7 | 0.085422 | 0.323925 |
| 0,8 | 0.06619 | 0.256576 |
| 0,9 | 0.0529213 | 0.208718 |
| 0.0437283 | 0.348932 |
| α |
|
|
| 0,1 | 0.663138 | 1.63184 |
| 0,2 | 0.447986 | 1.10052 |
| 0,3 | 0.308154 | 0.755087 |
| 0,4 | 0.216092 | 0.527538 |
| 0,5 | 0.154768 | 0.375825 |
| 0,6 | 0.113483 | 0.273535 |
| 0,7 | 0.085422 | 0.203842 |
| 0,8 | 0.06619 | 0.155894 |
| 0,9 | 0.0529213 | 0.122613 |
| 0.0437283 | 0.0993358 |
ОКНА ПРОСМОТРА ДАННЫХ
Для выделения определенных характеристик спектральных оценок, нередко прибегают к сглаживанию значений на концах случайного временного ряда. Временное сглаживание представляет собой умножение ряда на «окно данных».
При определении расширенного конечного преобразования Фурье, задаваемого соотношением

введена функция
, называемая окном просмотра данных (множителем сходимости, коэффициентом сглаживания).
Функцию
(3.1)
называют частотным окном. Из соотношения (3.1) вытекает, что

Характерное поведение функции
состоит в том, что она становится все более сконцентрированной в окрестности нуля при
.
Примеры окон просмотра данных:
1.
1 – окно Дирихле;
2.
1-
– окно Фейера;
3.
;
4.
– окно Хэннинга;
5.
– окно Хэмминга;
6.
– окно Хэмминга;
7.
, где
– окно Хэмминга;
8.
1-
– окно Рисса.
ЗАКЛЮЧЕНИЕ
В данной работе исследована оценка спектральной плотности вида

где
, а периодограмма задана следующим соотношением

Оценивается смещение данной спектральной плотности. Построены графики этой оценки для различных окон данных на основании наблюдений за солнечной активностью по Вольфу с 1749 г. по 1901 г.
Также построены графики для центрированного случайного процесса.