Презентация Случайные процессы (лекция 15). Параметрические модели временных рядов. Сглаживание и фильтрация онлайн

Содержание слайда: Лекция 15 Случайные процессы Параметрические модели временных рядов. Сглаживание и фильтрация. Моделирование искусственных гидрологических рядов (Ахметов С.К.)

Содержание слайда: Параметрические модели временных рядов Непараметрические методы - это методы описания СП с помощью корреляционных функций (это то, что делалось до сих пор) Параметрические методы – описание СП с помощью моделей авторегрессии и их комбинаций Модель авторегрессии В этой модели текущие значения СП выражаются в виде линейной комбинации предыдущих его значений и белого шума:   X’(t) = φ1X’(t-1) + φ2X’(t-2) + …… φpX’(t-p) + a(t)   X’(t) – центрированный СП; X’(t) = X(t) – mx a(t) – белый шум с нулевым математическим ожиданием и среднеквадратическим отклонение σа φ1, φ2…… φp – коэффициенты модели (константы) mx – математическое ожидание

Содержание слайда: Параметрические модели временных рядов (2) Модель содержит р + 2 параметров: φ1, φ2…… φp, mx, σа Эта модель авторегрессии р-го порядка. Она обозначается как АР(р) При р=1 X’(t) = φ1X(t-1) + a(t) Это модель называется моделью авторегрессии первого порядка АР(1) или модель Марковского процесса. Для этой модели коэффициент φ1 и ординаты автокорреляционной функции связаны соотношением   rk = φ1rk-1 при k > 0 Так как r0 = 1, то rk = φ1k   Таким образом, для АР(1) автокорреляционная функция полностью определяется своей первой ординатой. При этом φ1 = r1

Содержание слайда: Модель скользящего среднего Модель скользящего среднего м.б. получена из общей линейной модели (ОЛМ) при предположении, что ОЛМ содержит конечное число членов. При этом текущие значения СП X(t) выражаются через предыдущие значения белого шума at-1, at-2 …at-q   X’(t) = a(t) – η1 a(t-1) – η2 a(t-2) -….- ηq a(t-q) Модель содержит q + 2 параметров: η1, η2 ….- ηq, mx, σа X’(t) – центрированный случайный процесс, X’(t)= X(t) - mx a(t) - белый шум с нулевым математическим ожиданием и среднеквадратическим отклонение σа η1, η2 ….- ηq - коэффициенты модели (константы) Это выражение называется моделью скользящего среднего q - го порядка и обозначается CC(q)

Содержание слайда: Смешанные модели Иногда целесообразно объединять модели АР и СС. В этом случае получается смешанная модель ССАР (p,q), где р – порядок авторегрессии, q - порядок скользящего среднего. Выражение ССАР имеет вид   X’(t) = φ1X’(t-1) + φ2X’(t-2) + …… φpX’(t-p) + a(t) + a(t) – η1 a(t-1) – η2 a(t-2) -….- ηq a(t-q)   Такая модель может быть полезна, например, когда наблюдаемый временной ряд является суммой двух или более независимых составляющих, каждая из которых описывается либо моделью АР, либо моделью СС, но которые непосредственно не измеряются. При p = 1 и q=1, модель имеет вид   X’(t) = φ1X’(t-1) + a(t) – η1 a(t-1)

Содержание слайда: Модель авторегрессии и проинтегрированного скользящего среднего (АРПСС) Модели АР, СС и АРСС относятся к классу стационарных моделей, которые описывают процессы с постоянными МО и дисперсиями Если ряды нестационарны, то от нестационарности можно избавиться, заменяя исходный ряд на ряд разностей: Y(t) = X(t) – X(t-1)   Если от нестационарности не удается избавиться, то можно взять разность повторно После проведенных преобразований к исходному ряду можно применить модель АРСС Такую модель называют моделью АРПСС (p,d,q). В этой модели: p - порядок авторегрессии d - порядок разности q - порядок скользящего среднего

Содержание слайда: Сглаживание и фильтрация Если временной ряд содержит некоторые частоты и периоды, которые в данный момент не представляют интереса, то амплитуда этих волн может быть уменьшена с помощью статистической фильтрации. При этом изменяется спектр исходного ряда. Одной из форм статистической фильтрации может быть сглаживание, в которой спектральные компоненты с высокой частотой уменьшены. Такой фильтр называется низкочастотным. Простейшим статистическим фильтром является скользящая средняя с равными весами. Скользящее среднее рассчитывается путем суммирования n последовательных величин временного ряда и делением полученной суммы на n.

Содержание слайда: Оценка характеристик СП по эмпирическим данным Определение характеристик СП по множеству реализаций  

Содержание слайда: Оценка характеристик СП по эмпирическим данным Для стационарного СП математическое ожидание, дисперсия и СКО являются константами и их можно оценить по любому сечению. Корреляционная функция стационарного процесса не зависит от моментов времени рассматриваемого процесса и зависит только от расстояния между сечениями τ. Поэтому ее оценивают по формуле

Содержание слайда: Определение характеристик стационарного эргодического процесса по одной реализации Если для стационарного СП принимается гипотеза об эгродичности, то МО, дисперсию и СКО можно оценить по одной (достаточно продолжительной) реализации. Тогда получим, что

Содержание слайда: Моделирование искусственных гидрологических рядов В основе моделирования искусственных рядов лежит метод Монте – Карло. Это метод решения математических задач при помощи моделирования случайных чисел Процесс моделирования включает следующие шаги: Нужно получить последовательность случайных чисел, равномерно распределенных в интервале от 0 до 1 Каждое значение случайного числа рассматривается как вероятность не превышения и по нему рассчитывается соответствующий квантиль заданного закона распределения. Можно это сделать по графику, можно по таблицам или по готовым компьютерным программа.

Содержание слайда: СПАСИБО ЗА ВНИМАНИЕ!