
- •Введение
- •1.Главное и текущие окна анализа данных
- •2.Электронная таблица
- •Пример 4
- •Пример 6
- •3.Графическое представление данных
- •4.Пользовательский интерфейс текущего окна анализа
- •Общее знакомство с пакетом Statgraphics. Генерация временных рядов
- •1.1 Генерация тренда временного ряда
- •Г) Рассчитайте размах тренда range и масштабный коэффициент (переменные range и scale соответственно).
- •Е) Сохраните график в отчетном документе (StatGallery).
- •1.2 Генерация реализаций абсолютно случайного временного ряда
- •1.3 Генерация временных рядов pk , sk , k
- •1.4 Моделирование грубых сбоев измерений
- •1.5 Генерация обобщенной реализации временного ряда
- •Статистический анализ стационарных временных рядов
- •2.1 Подготовка данных для статистических исследований
- •А) Напишите расчетные выражения для переменных x1 – x5 с учетом синтаксиса языка пакета statgraphics Plus 5.1 и числовых значений параметров алгебраических выражений своего индивидуального задания.
- •2.2 Анализ одномерного закона распределения вероятностей
- •2.3 Сравнение экспериментального и теоретического
- •2.4 Изучение описательных статистик стационарного
- •Для временных рядов x1 – x5
- •2.5 Исследование автокорреляционной функции стационарного временного ряда
- •Анализ независимости значений временного ряда по одной реализации
- •3.1 Подготовка данных для статистических исследований
- •Результат визуального анализа временных рядов
- •3.2 Анализ независимости временного ряда по критерию
- •Результаты анализа независимости временных рядов
- •3.3 Анализ независимости временного ряда по критерию
- •Результаты анализа независимости временных рядов
- •3.4 Анализ независимости временного ряда
- •Результаты анализа независимости временных рядов
- •3.5 Анализ интервала корреляции для принятия решения
- •Результаты анализа независимости временных рядов
- •3.6 Заключение о точности и области практического применения исследованных методов
- •Сглаживание временных рядов линейными и нелинейными фильтрами
- •4.1 Исследование возможностей медианного фильтра
- •4.2 Анализ чувствительности линейного фильтра к выбросам входных данных
- •4.3 Изучение импульсной реакции линейных сглаживающих фильтров
- •4.4 Исследование динамической ошибки
- •4.5 Исследование коэффициента подавления помехи линейными сглаживающими фильтрами
- •Оптимальная обработка данных линейным фильтром с конечной памятью
- •5.1 Расчет оптимального линейного фильтра с конечной памятью
- •5.2 Исследование свойств оптимального фильтра
- •Спектральный анализ временных рядов конечной длины
- •Исследование периодограммы импульсного сигнала
- •6.2 Спектральный анализ гармонической функции
- •Расчет амплитудно-частотных характеристик
- •Исследование периодограммы реализации белого шума
- •Спектральный анализ низкочастотного временного ряда, измеренного на фоне помех
3.5 Анализ интервала корреляции для принятия решения
о независимости временного ряда
Установите режим анализа автокорреляций в окне Tabular Options. По значениям, представленным в таблице, или по графику автокорреляционной функции (опция Autocorrelation Function в меню Graphical Options) определите экспериментальную оценку интервала корреляции анализируемого временного ряда. Заполните таблицу по результатам исследования автокорреляций всех временных рядов.
Результаты анализа независимости временных рядов
по значению интервала корреляции
№ п/п |
Имя переменной |
Теор. свойство (H0 – белый шум, __ H0 - не белый шум) |
Интервал корреляции |
Результат анализа (H0 – белый шум, __ H0- не белый шум) |
Ошибка принятия решения (ошибка отсутствует, ошибка первого рода, ошибка второго рода) |
1 |
noise1 |
|
|
|
|
2 |
x3 |
|
|
|
|
3 |
Slow |
|
|
|
|
4 |
x4 |
|
|
|
|
5 |
l_noise |
|
|
|
|
6 |
g_noise |
|
|
|
|
3.6 Заключение о точности и области практического применения исследованных методов
Из проведенных исследований следует, что
-
для временного ряда с линейным трендом рекомендуется применять методы _________________________________
__________________________________________________ ,
-
для временного ряда с гармоническим трендом рекомендуется применять методы ____________________
__________________________________________________ ,
-
универсальными с точки зрения фактических свойств исследуемого временного ряда являются методы ________
__________________________________________________ ,
-
наибольшее число ошибок наблюдалось при использовании методов _____________________________
_________________________________________________ .
Дата выполнения работы ______________________
Подпись преподавателя _______________________
Сглаживание временных рядов линейными и нелинейными фильтрами
Вычисление низкочастотной составляющей зашумленного временного ряда, которая, как правило, содержит информацию о его внутренних закономерных динамических свойствах, является одной из наиболее распространенных в приложениях задач обработки данных. Сложность этой задачи обусловлена тем, что подавление высокочастотных помех не должно сопровождаться существенными искажениями низкочастотного «полезного» сигнала. Дополнительную сложность в обработку данных могут вносить грубые сбои, которые должны автоматически исключаться при выполнении сглаживания временного ряда. В противном случае возникают недопустимые искажения его истинной динамики.
В лабораторной работе изучаются нелинейные фильтры, ориентированные на удаление выбросов в реализации временного ряда, и линейные сглаживающие фильтры. На модельных данных исследуются динамические и случайные ошибки обработки данных различными фильтрами, реализованными в пакете STATGRAPHICS Plus 5.1. Теоретические основы лабораторной работы изложены в разделах 5 и 8 учебного пособия.
В пакете STATGRAPHICS Plus 5.1 предоставляется возможность применить последовательно два фильтра для обработки временного ряда. Каждый из этих последовательно подключенных фильтров выбирается из следующего списка:
- Simple Moving Average (SMA) – простой фильтр скользящего
среднего,
- Spencer’s 15-Term MA – 15-точечный фильтр Спенсера,
- Spencer’s 21-Term MA – 21-точечный фильтр Спенсера,
- Henderson’s Weighted MA – фильтр Гендерсона взвешенного
скользящего среднего,
- EWMA – фильтр экспоненциально взвешенного скользящего
среднего (Exponential Weighted Moving Average, RC- фильтр
или фильтр Брауна),
- 3RSS, 3RSSH, 5RSS, 5RSSH, 3RSR – робастные медианные сглаживающие фильтры (Resistant Nonlinear Smoothing).
В обозначениях вариантов реализации робастных медианных сглаживающих фильтров первая цифра означает ширину окна просмотра входных данных, последующие символы имеют следующее содержание: R – сглаживание, S – процедура расщепления для устранения эффекта образования плоских участков на выходе медианного фильтра, SS – двукратное применение S-процедуры, H – сглаживание фильтром «взвешенного» среднего в окне ширины 3 с весовыми коэффициентами 0.25; 0.5; 0.25 (фильтр Хэннинга).
Фильтр EWMA, который в литературе часто называется фильтром Брауна (Brown), реализует следующий алгоритм сглаживания:
,
что соответствует решению разностного уравнения
.
Этот фильтр характеризуется коэффициентом усиления, равным 1.
В пакете STATSGRAPHICS Plus 5.1 параметром EWMA фильтра, определяющим его сглаживающие свойства, является = 1-.