Всё сдал! - помощь студентам онлайн Всё сдал! - помощь студентам онлайн

Реальная база готовых
студенческих работ

Узнайте стоимость индивидуальной работы!

Вы нашли то, что искали?

Вы нашли то, что искали?

Да, спасибо!

0%

Нет, пока не нашел

0%

Узнайте стоимость индивидуальной работы

это быстро и бесплатно

Получите скидку

Оформите заказ сейчас и получите скидку 100 руб.!


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

Тип Реферат
Предмет Математика
Просмотров
320
Размер файла
427 б
Поделиться

Ознакомительный фрагмент работы:

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

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

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

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

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

Брест 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,10.6631382.13239
0,20.4479861.48005
0,30.3081541.04694
0,40.2160920.7554
0,50.1547680.556644
0,60.1134830.41954
0,70.0854220.323925
0,80.066190.256576
0,90.05292130.208718
10.04372830.348932
α
0,10.6631381.63184
0,20.4479861.10052
0,30.3081540.755087
0,40.2160920.527538
0,50.1547680.375825
0,60.1134830.273535
0,70.0854220.203842
0,80.066190.155894
0,90.05292130.122613
10.04372830.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) для окна Рисса для центрированного случайного процесса


Нет нужной работы в каталоге?

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

Цены ниже, чем в агентствах и у конкурентов

Вы работаете с экспертами напрямую. Поэтому стоимость работ приятно вас удивит

Бесплатные доработки и консультации

Исполнитель внесет нужные правки в работу по вашему требованию без доплат. Корректировки в максимально короткие сроки

Гарантируем возврат

Если работа вас не устроит – мы вернем 100% суммы заказа

Техподдержка 7 дней в неделю

Наши менеджеры всегда на связи и оперативно решат любую проблему

Строгий отбор экспертов

К работе допускаются только проверенные специалисты с высшим образованием. Проверяем диплом на оценки «хорошо» и «отлично»

1 000 +
Новых работ ежедневно
computer

Требуются доработки?
Они включены в стоимость работы

Работы выполняют эксперты в своём деле. Они ценят свою репутацию, поэтому результат выполненной работы гарантирован

avatar
Математика
История
Экономика
icon
142645
рейтинг
icon
3069
работ сдано
icon
1329
отзывов
avatar
Математика
Физика
История
icon
140338
рейтинг
icon
5849
работ сдано
icon
2647
отзывов
avatar
Химия
Экономика
Биология
icon
94648
рейтинг
icon
2019
работ сдано
icon
1266
отзывов
avatar
Высшая математика
Информатика
Геодезия
icon
62710
рейтинг
icon
1046
работ сдано
icon
598
отзывов
Отзывы студентов о нашей работе
52 911 оценок star star star star star
среднее 4.9 из 5
игу
работа выполнена качественно, а самое главное быстро. всё требования были выполнены
star star star star star
ФЭК
Прекрасная Татьяна, выполнила все быстро, профессионально и раньше срока!!! Очень рекомендую 👌
star star star star star
Санкт-Петербургский государственный архитектурно-строительный университет
Выполнено быстро, качественно, преподаватель оценил по высшему баллу. Рекомендую.
star star star star star

Последние размещённые задания

Ежедневно эксперты готовы работать над 1000 заданиями. Контролируйте процесс написания работы в режиме онлайн

Сделать 15 вариант из методички

Решение задач, Математический анализ

Срок сдачи к 5 июня

только что

Моделирование бизнес-процессов операторов связи

Другое, Проектирование и эксплуатация систем связи

Срок сдачи к 17 июня

только что

Написать реферат

Реферат, Физика

Срок сдачи к 29 мая

только что
только что

Написать реферат с 70 процентами антиплагиата. Требования указала в файле.

Реферат, Юридическая аргументация

Срок сдачи к 28 мая

только что

Анализ главного героя одного из сериалов на выбор

Контрольная, Психоанализ

Срок сдачи к 31 мая

1 минуту назад

Современные концепции управления персоналом

Курсовая, Основы управления персоналом

Срок сдачи к 10 июня

1 минуту назад

решить задачу

Решение задач, Бухгалтерский учет в страховых организациях

Срок сдачи к 29 мая

1 минуту назад
1 минуту назад

Расчетно-аналитическая работа

Другое, Диджитализация, финтех и инновации в финансовых институтах

Срок сдачи к 9 июня

1 минуту назад

Решить задачи

Решение задач, Статистика

Срок сдачи к 29 мая

2 минуты назад

Тема Патриотическое воспитание школьников в учебно-воспитательном...

Курсовая, Педагогика и психология

Срок сдачи к 2 июня

2 минуты назад

написать реферат

Реферат, «Реконструкция зданий, сооружений и застройки»

Срок сдачи к 31 мая

2 минуты назад

Стилизовать текст

Другое, История русского литературного языка

Срок сдачи к 27 мая

2 минуты назад

4 задачи

Курсовая, Внутризаводские и электротехнические комплексы и системы

Срок сдачи к 7 июня

3 минуты назад

Часть работы уже написана, есть содержание, 1 и 2 главы

Диплом, государственное и муниципальное управление

Срок сдачи к 6 июня

3 минуты назад

Написание курсовой работы строго по методтички

Курсовая, анализ финансово-хозяйственной деятельности

Срок сдачи к 13 июня

3 минуты назад
3 минуты назад
planes planes
Закажи индивидуальную работу за 1 минуту!

Размещенные на сайт контрольные, курсовые и иные категории работ (далее — Работы) и их содержимое предназначены исключительно для ознакомления, без целей коммерческого использования. Все права в отношении Работ и их содержимого принадлежат их законным правообладателям. Любое их использование возможно лишь с согласия законных правообладателей. Администрация сайта не несет ответственности за возможный вред и/или убытки, возникшие в связи с использованием Работ и их содержимого.

«Всё сдал!» — безопасный онлайн-сервис с проверенными экспертами

Используя «Свежую базу РГСР», вы принимаете пользовательское соглашение
и политику обработки персональных данных
Сайт работает по московскому времени:

Вход
Регистрация или
Не нашли, что искали?

Заполните форму и узнайте цену на индивидуальную работу!

Файлы (при наличии)

    это быстро и бесплатно