Оценивание смещения статистики взаимной спектральной плотности многомерного временного ряда

Главная задача спектрального анализа временных рядов. Параметрические и непараметрические методы спектрального анализа. Сущность понятия "временный ряд". График оценки спектральной плотности для окна Дирихле, при центрированном случайном процессе.

Рубрика Математика
Вид курсовая работа
Язык русский
Дата добавления 17.09.2009
Размер файла 332,8 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

23

Математический факультет

Кафедра информатики и прикладной математики

Курсовая работа НА ТЕМУ:

«ОЦЕНИВАНИЕ СМЕЩЕНИЯ СТАТИСТИКИ ВЗАИМНОЙ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ МНОГОМЕРНОГО ВРЕМЕННОГО РЯДА»

Брест 2009

СОДЕРЖАНИЕ

  • ВВЕДЕНИЕ
  • 1. ОСНОВНЫЕ ОПРЕДЕЛЕНИЯ, ИСПОЛЬЗУЕМЫЕ В РАБОТЕ
  • 2. ОЦЕНИВАНИЕ СМЕЩЕНИЯ СТАТИСТИКИ ВЗАИМНОЙ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ
  • 3. ОКНА ПРОСМОТРА ДАННЫХ
  • ЗАКЛЮЧЕНИЕ
  • СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
  • ПРИЛОЖЕНИЕ

ВВЕДЕНИЕ

Современный этап развития теории вероятностей и математической статистики характеризуется значительным расширением теоретических исследований по статистическому спектральному анализу (анализу в частной области) временных рядов и их практическим применением во многих областях человеческой деятельности, таких, как экономика, спектроскопия, медицина, биология, страхование, финансы, социология, радиоэлектроника, электротехника, геофизика, геология и многие другие. Цели изучения временных рядов могут быть различными. Можно, например, стремиться предсказывать будущее на основании знания прошлого, управлять процессом, порождающим ряд, выяснить механизм, порождающий ряд, или просто сжато описать характерные особенности ряда. Поэтому под статистическим спектральным анализом временных рядов понимают статистический спектральный анализ стационарных случайных процессов.

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

Существуют параметрические и непараметрические методы спектрального анализа. Среди непараметрических методов выделяют метод, в котором для построения оценки спектральной плотности производится осреднение периодограмм, построенных по непересекающимся интервалам исходной последовательности наблюдений и вводятся окна просмотра данных для уменьшения смещения оценок.

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

Также построены графики для центрированного случайного процесса.

1. ОСНОВНЫЕ ОПРЕДЕЛЕНИЯ, ИСПОЛЬЗУЕМЫЕ В РАБОТЕ

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

Совокупность функций вида

назовём r-компонентным векторным временным рядом (r-мерным временным рядом).

Переменная t обычно соответствует времени выполнения или регистрации наблюдений и измерений.

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

Если , или подмножество из , то говорят, что , - случайный процесс с дискретным временем.

Если , или подмножество из , то говорят, что , - случайный процесс с непрерывным временем.

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

При каждом фиксированном , , - множество случайных величин.

Если в определении случайного процесса , , , то =называется -мерным случайным полем.

Введем характеристики случайного процесса , , во временной области.

Математическим ожиданием случайного процесса , , называется функция вида

, .

Дисперсией случайного процесса , , называется функция вида

, .

Корреляционной функцией случайного процесса , , называется функция вида

, .

Ковариационной функцией случайного процесса , , называется функция вида

.

Заметим, что если , то , .

Смешанным моментом го порядка, , случайного процесса , , называется функция вида

, , .

Заметим, что

, .

Пусть - значения случайного процесса в точках .

Смешанный момент го порядка, , можно также определить как

, , .

Смешанным семиинвариантом ( кумулянтом ) го порядка, , случайного процесса , , называется функция вида

, , ,

которую также будем обозначать как .

Между смешанными моментами и смешанными семиинвариантами го порядка, , существуют связывающие их соотношения, которые имеют вид

(1.1)

(1.2)

суммирование по всевозможным разбиениям множества .

Спектральной плотностью случайного процесса , , называется функция вида

=, ,

при условии, что

.

Из определения видно, что спектральная плотность непрерывная, периодическая функция с периодом, равным по каждому из аргументов.

Семиинвариантной спектральной плотностью го порядка, , случайного процесса , , называется функция вида

=, ,

при условии, что

.

Лемма 1. Для любого целого справедливо соотношение

(1.3)

Теорема 1. Для смешанного семиинварианта го порядка, , случайного процесса справедливы представления

, (1.4)

Доказательство. Домножая обе части соотношения (1.1) на

, ,

и интегрируя обе части полученного неравенства по на , получим

.

Используя лемму 1, получим при требуемый результат. Теорема доказана.

Лемма 2. Если функция интегрируема и периодична с периодом , то для любого действительного имеет место соотношение

Доказательство. Предположим, что >0. Можно записать

В третьем слагаемом правой части последнего равенства сделаем замену переменных интегрирования и, учитывая периодичность с периодом функции , получаем требуемое. Случай, когда <0, доказывается аналогично. Лемма доказана.

Спектральной плотностью случайного процесса , , называется функция вида

=, ,

при условии, что

.

Из определения видно, что спектральная плотность непрерывная, периодическая функция с периодом, равным по каждому из аргументов.

2. ОЦЕНИВАНИЕ СМЕЩЕНИЯ СТАТИСТИКИ ВЗАИМНОЙ СПЕКТРАЛЬНОЙ ПЛОТНОСТИ

Рассмотрим действительный стационарный в широком смысле случайный процесс,, с математическим ожиданием , , взаимной ковариационной функцией , и взаимной спектральной плотностью .

Предположим, имеются Т последовательных, полученных через равные промежутки времени наблюдений за составляющей , рассматриваемого процесса . Как оценку взаимной спектральной плотности в точке рассмотрим статистику

(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

1

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

1

0.0437283

0.0993358

3. ОКНА ПРОСМОТРА ДАННЫХ

Для выделения определенных характеристик спектральных оценок, нередко прибегают к сглаживанию значений на концах случайного временного ряда. Временное сглаживание представляет собой умножение ряда на «окно данных».

При определении расширенного конечного преобразования Фурье, задаваемого соотношением

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

Функцию

(3.1)

называют частотным окном. Из соотношения (3.1) вытекает, что

Характерное поведение функции состоит в том, что она становится все более сконцентрированной в окрестности нуля при .

Примеры окон просмотра данных:

1. 1 - окно Дирихле;

2. 1- - окно Фейера;

3. ;

4. - окно Хэннинга;

5. - окно Хэмминга;

6. - окно Хэмминга;

7. , где - окно Хэмминга;

8. 1- - окно Рисса.

ЗАКЛЮЧЕНИЕ

В данной работе исследована оценка спектральной плотности вида

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

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

Также построены графики для центрированного случайного процесса.

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

1. Бриллинджер Д. Временные ряды. Обработка данных и теория. - М.: Мир, 1980. - 536 с.

2. Андерсон Т. Статистический анализ временных рядов. - М.: Мир, 1976. - 755 с.

3. Труш Н.Н. Асимптотические методы статистического анализа временных рядов. - Мн.: БГУ, 1999. - 218 с.

4. Журбенко И.Г. Спектральный анализ временных рядов. - М.: Изд-во МГУ, 1982. - 168 с.

5. Труш Н.Н., Мирская Е.И. Случайные процессы. Преобразования Фурье наблюдений. - Мн.: БГУ, 2000.

ПРИЛОЖЕНИЕ

Для исследования оценки (2.1) был исследован ряд, состоящий из 100 наблюдений за солнечной активностью по Вольфу с 1749 г. по 1901 г.

Рис. 1 - График оценки спектральной плотности (2.1) для окна Дирихле

Рис. 2 - График оценки спектральной плотности (2.1) для окна Дирихле для центрированного случайного процесса

Рис. 3 - График оценки спектральной плотности (2.1) для окна Фейера

Рис. 4 - График оценки спектральной плотности (2.1) для окна Фейера для центрированного случайного процесса

Рис. 5 - График оценки спектральной плотности (2.1) для окна вида 3

Рис. 6 - График оценки спектральной плотности (2.1) для окна вида 3 для центрированного случайного процесса

Рис. 7 - График оценки спектральной плотности (2.1) для окна Хэннинга

Рис. 8 - График оценки спектральной плотности (2.1) для окна Хэннинга для центрированного случайного процесса

Рис. 9 - График оценки спектральной плотности (2.1) для окна Хэмминга вида 5

Рис. 10 - График оценки спектральной плотности (2.1) для окна Хэмминга вида 5 для центрированного случайного процесса

Рис. 11 - График оценки спектральной плотности (2.1) для окна Хэмминга вида 6

Рис. 12 - График оценки спектральной плотности (2.1) для окна Хэмминга вида 6 для центрированного случайного процесса

Рис.13 - График оценки спектральной плотности (2.1) для окна Хэмминга вида 7

Рис. 14 - График оценки спектральной плотности (1) для окна Хэмминга вида 7 для центрированного случайного процесса

Рис. 15 - График оценки спектральной плотности (1) для окна Рисса

Рис. 16 - График оценки спектральной плотности (1) для окна Рисса для центрированного случайного процесса


Подобные документы

  • Исследование первого момента состоятельной оценки взаимной спектральной плотности. Задачи спектрального анализа временных рядов. Графики оценки для временного ряда, представляющего собой последовательность наблюдений температуры воздуха в городе Бресте.

    курсовая работа [324,9 K], добавлен 16.08.2011

  • Математическое ожидание и дисперсия случайного процесса. Спектральная плотность случайного процесса. Сглаживание значений на концах случайного временного ряда. График оценки спектральной плотности для окна Рисса, при центрированном случайном процессе.

    курсовая работа [382,3 K], добавлен 17.09.2009

  • Первые два момента состоятельной оценки спектральной плотности, исследование асимптотического поведения математического ожидания и дисперсии построенной оценки. Сравнительный анализ оценки спектральной плотности в зависимости от окон просмотра данных.

    курсовая работа [558,0 K], добавлен 12.04.2012

  • Понятие об основной тенденции ряда динамики, ее сущность и визуальное представление, методы анализа. Аналитическая оценка уравнения тренда. Характеристика, использование различных методов для выделения тренда временных рядов, прогнозирование показателей.

    курсовая работа [207,2 K], добавлен 04.03.2013

  • Точечное оценивание основных числовых характеристик, функции и плотности распределения компонент многомерного случайного вектора. Статистическая проверка характера распределения. Особенности корреляционного анализа признаков этой математической категории.

    курсовая работа [1,1 M], добавлен 01.10.2013

  • Преимущества и недостатки параметрических методов оценки. Процедура Роббинса-Монро, алгоритмы Литвакова и Кестена. Исследование стохастических аппроксимаций непараметрического типа. Непараметрическая оценка плотности вероятности и кривой регрессии.

    реферат [470,6 K], добавлен 22.04.2014

  • Ознакомление с математическим аппаратом анализа временных рядов и моделями авторегрессии. Составление простейших моделей авторегрессии стационарных временных рядов. Оценка дисперсии и автоковариации, построение графика автокорреляционной функции.

    лабораторная работа [58,7 K], добавлен 14.03.2014

  • Исследование задачи Дирихле для вырождающегося уравнения смешанного типа в прямоугольной области методами спектрального анализа. Обоснование корректности постановки нелокальных начально-граничных задач различных вырождающихся дифференциальных уравнений.

    курсовая работа [135,1 K], добавлен 06.05.2011

  • Описание сущности функции, которая была введена немецким математиком П.В. Дирихле как пример функции, свободной от аналитического задания значения. Характеристика и описание ряда ее свойств и области определения методами математического анализа.

    курсовая работа [44,8 K], добавлен 23.11.2011

  • Изучение изменений анализируемых показателей во времени как важнейшая задача статистики. Понятие рядов динамики (временных рядов). Числовые значения того или иного статистического показателя, составляющего ряд динамики. Классификация рядов динамики.

    презентация [255,0 K], добавлен 28.11.2013

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.