Информационная система
«Ёшкин Кот»

XXXecatmenu

ОТРАСЛЕВОЙ СТАНДАРТ

ОТРАСЛЕВАЯ
АВТОМАТИЗИРОВАННАЯ СИСТЕМА УПРАВЛЕНИЯ.
ПОДСИСТЕМА УПРАВЛЕНИЯ КАЧЕСТВОМ

Построение математических моделей временных рядов показателей

ОСТ 1.00321-78

На 25 страницах

Введен впервые

Проверено в 1982 г.

Распоряжением Министерства от 26 декабря 1978 г. № 087-16
срок введения установлен с 1 июля 1979 г.

Настоящий стандарт распространяется на теоретические методы моделирования временных рядов показателей, закладываемых в отраслевой автоматизированной системе управления (ОАСУ).

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

1. ОБЩИЕ ПОЛОЖЕНИЯ

1.1. Временные ряды показателей строятся по результатам контроля изделий на этапах производства и эксплуатации. При этом вероятностные оценки значений показателей должны быть состоятельными, а временные ряды представительными.

1.2. Методы, используемые при построении математических моделей, инвариантны к видам показателей и этапам - «жизненного цикла» изделий.

1.3. Стандарт позволяет осуществлять построение моделей как стационарных, так и нестационарных со стационарными n-ми приращениями временных рядов.

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

1.5. Допускается самостоятельное использование алгоритмов и программ, соответствующих основным этапам построения математической модели. Основные этапы построения математической модели:

- проверка ряда на стационарность;

- идентификация пробной модели;

- оценка параметров модели;

- проверка адекватности модели.

2. МЕТОД МОДЕЛИРОВАНИЯ ВРЕМЕННЫХ РЯДОВ

2.1. Общая блок-схема процесса построения математической модели временного ряда представлена на черт. 1.

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


Черт. 1


Гипотеза о постоянстве значений первых двух моментов подтверждается, если их изменение во времени не превышает стандартной ошибки  (где N - длина временного ряда). Если гипотеза не подтверждается, то осуществляется взятие разностей

                                                             (1)

при t = 1, 2, ..., N,

где  - значения показателя, представленные в виде временного ряда;

Рt - значения, полученные в результате взятия разностей.

Число процедур взятия разностей до приведения ряда к стационарному виду определяет порядок нестационарности ряда d (d £ 4). При этом длина стационарного временного ряда становится равной

N = N0 - d.                                                                (2)

2.3. При построении математических моделей временных рядов используются вспомогательные операторы сдвига назад B и сдвига вперед F, определяемые как

                                                       (3)

В общем виде

                                                    (4)

где К - число шагов сдвига (К = 0, 1, ..., n).

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

2.4.1. Модель авторегрессии АР

                                                            (5)

где  - отклонение значений ряда от его среднего значения m;

at - импульс «белого шума»;

Ф (В) = 1 - Ф1 В - Ф2 В2 - ... - Фr Вr - оператор авторегрессии;

Ф1, ..., Фr - параметры авторегрессии;

r - порядок процесса авторегрессии.

2.4.2. Модель скользящего среднего СС

                                                              (6)

где q (В) = 1 - q1 В - q2 В2 - ... - qq Вq - оператор скользящего среднего;

q1, ..., qq - параметры скользящего среднего;

q - порядок процесса скользящего среднего.

2.4.3. Смешанная модель АРСС

                                                         (7)

Порядок смешанной модели определяется как сумма (r + q).

2.4.4. Модель авторегрессии проинтегрированного скользящего среднего АРПСС

                                                          (8)

где j (B) = Ф (B) (1 - B)d - обобщенный оператор авторегрессии;

d - порядок нестационарности процесса.

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

Гипотеза о равенстве нулю установившихся значений частной автокорреляционной функции подтверждается, если эти значения не превышают по модулю стандартной ошибки выборочной частной автокорреляции.

При проверке гипотезы о равенстве нулю установившихся значений автокорреляционной функции используется стандартная ошибка Бартлетта.

Анализ спектров автокорреляционной и частной автокорреляционной функций состоит в отнесении их к одному из видов:

слабозатухающий спектр - спектр, имеющий на рассматриваемом интервале Кзад более половины ненулевых значений, причем первые 6 из них находятся в начале спектра;

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

2.6. Основные результаты анализа сводятся к следующему. Если спектры обеих функций содержат только первые составляющие, а остальные равны нулю, рассматриваемый процесс аппроксимируется «белым шумом».

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

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

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

При выборе вида математической модели должно учитываться то, что в случае высокого порядка (2 - 3) моделей авторегрессии или скользящего среднего предпочтительнее строить смешанную модель.

Примеры выбора вида математической модели приведены в справочном приложении 1.

2.7. Оценка параметров полученной модели производится в соответствии с критерием наименьших квадратов.

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

                                               (9)

где  - соответственно векторы параметров авторегрессии, скользящего среднего и значений временного ряда. Безусловная сумма квадратов находится суммированием квадратов всех значений последовательности импульсов at, вычисленных из системы уравнений

                        (10)

a-j = 0,

при j = T,

где et - последовательность независимо распределенных случайных импульсов, имеющих нулевое среднее значение и дисперсию

 - приведенные значения стационарного временного ряда;

T - момент времени, в который оценки , генерированные возвратной моделью  (F - оператор сдвига вперед) практически равны нулю;

Ф1, ...,Фr - параметры авторегрессии;

r - порядок процесса авторегрессии;

q1, ...,qq - параметры скользящего среднего;

q - порядок процесса скользящего среднего.

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

                                          (11)

где g1, ...,gк - параметры модели;

к = r + g - порядок модели;

                                                       (12)

2.8. Определение доверительных интервалов оценок наименьших квадратов осуществляется исходя из того, что приближенная (1 - e) % доверительная область оценок параметров ограничена изолинией суммы квадратов, для которой

                                                      (13)

где c2 - квантиль уровня (1 - e) распределения c2 с к = r + q, степенями свободы.

2.9. Гипотеза об адекватности полученной модели реальному временному ряду подтверждается, если величина Q распределена приблизительно как c2 (кy - r - q),

                                                        (14)

при n = N - d,

где Ra - автокорреляция импульсов аt.

Величина кy такая, при которой веса y модели, представленной в виде

                                               (15)

будут пренебрежимо малы.

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

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

3. АЛГОРИТМ ПРОВЕРКИ РЯДА НА СТАЦИОНАРНОСТЬ

3.1. Входные данные:

- число наблюдений временного ряда N;

- значения временного ряда {Pt}, t = 1, ..., N.

3.2. Вычисление статистических характеристик проводится следующим образом.

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

- математическое ожидание:

                                                       (16)

- автоковариационная функция:

                                                     (17)

где  - значения приведенного временного ряда;

l - сдвиг, l = 0, 1, ..., кзад,

- автокорреляционная функция:

Rl = Cl / C0, l = 0, 1, ..., кзад,                                           (18)

где С0 - дисперсия ряда.

3.3. Проверка временного ряда на стационарность по М0, C0, Rl для l = 1, ..., кзад производится по условиям

                                      (19)

где i = 1, ..., Ny - номера участков;

Ny - число участков.

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

3.4. Получение разностного ряда проводится в следующем порядке.

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

                                                  (20)

где j = 1, ..., d - номер процедуры взятия разностей.

Если d > 4, считается, что процесс не может быть приведен к стационарному виду.

3.5. Выходные данные:

- число наблюдений временного ряда N;

- значения временного ряда {Pt}, t = 1, ..., N;

- порядок разности временного ряда d.

4. АЛГОРИТМ ИДЕНТИФИКАЦИИ

4.1. Входные данные:

- число наблюдений временного ряда N;

- значения временного ряда {Pt}, t = 1, ..., N;

- порядок разности стационарного временного ряда d.

4.2. Приведение ряда к стационарному виду осуществляется при условии d ¹ 0. В этом случае производится процедура взятия разностей по формуле (20).

4.3. Вычисление статистических характеристик проводится следующим образом.

По формулам (16) - (18) осуществляется вычисление соответственно математического ожидания, дисперсии и автокорреляционной функции. Частная автокорреляционная функция вычисляется по формулам

                                                  (21)

при l > 1;

                                          (22)

при m = 1, 2, ..., l - 1.

l = 2, 3, ..., кзад.

4.4. Вычисление числа x, характеризующего отклонение от нуля установившегося значения частной автокорреляционной функции, осуществляется по формуле

                                                            (23)

4.5. Вычисление стандартной ошибки Бартлетта осуществляется по формуле

                                                  (24)

где i - предполагаемый порядок авторегрессии.

4.6. Анализ спектров автокорреляционной и частной автокорреляционной функций

4.6.1. Спектр считается быстрозатухающим, если выполняются следующие условия:

i £ 3 и к £ кзад / 2,

где i - номера начальных ненулевых составляющих спектра;

к - число ненулевых составляющих спектра.

Спектр считается слабозатухающим, если выполняются следующие условия:

5 £ i и к £ кзад / 2,

где i - номера начальных нулевых составляющих спектра;

К - число нулевых составляющих спектра.

4.7. Выбор математической модели (порядок модели не более двух)

4.7.1. Если оба спектра быстрозатухающие содержат только первую составляющую, то процесс описывается моделью «белого шума» с mа, sа2.

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

r = i - 1,                                                              (25)

где i - число первых ненулевых составляющих частной автокорреляционной функции.

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

q = i - 1,                                                             (26)

где i - число первых ненулевых составляющих автокорреляционной функции.

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

4.8. Оценка параметров модели авторегрессии

4.8.1. Для модели авторегрессии 1-го порядка при выполнении условия

/R1/ < 1,                                                              (27)

параметр авторегрессии равен Ф = R1. В противном случае рассматривается смешанная модель 2-го порядка.

4.8.2. Для модели авторегрессии 2-го порядка при выполнении условия

                                                    (28)

параметры авторегрессии находятся из уравнений

                                              (29)

В противном случае рассматривается смешанная модель 2-го порядка.

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

4.9.1. Для модели скользящего среднего 1-го порядка параметр скользящего среднего определяется по формуле

                                                   (30)

Из двух значений q выбирается то, которое удовлетворяет условиям

                                                     (31)

Если ни один из корней не удовлетворяет первому условию, рассматривается смешанная модель 2-го порядка.

4.9.2. Для модели скользящего среднего 2-го порядка при выполнении условий

                                                     (32)

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

                                                     (33)

В противном случае рассматривается смешанная модель 2-го порядка.

4.10. Оценка параметров смешанной модели 2-го порядка

4.10.1. Если выполняются условия

                                               (34)

параметры модели определяются из системы уравнений

                                                  (35)

В противном случае рассматривается смешанная модель выше 2-го порядка.

4.11. Оценка параметров смешанной модели выше 2-го порядка

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

                                                        (36)

при Вlj = C/q + i - j/,

Xi = Cq + i,

i = j = 1, ..., r.

4.11.2. Оценки параметров скользящего среднего определяются в соответствии с алгоритмом Ньютона-Рафсона, представляющего собой следующую итеративную процедуру:

- вычисление модифицированных автоковариаций временного ряда

                                               (37)

для j = 1, ..., q,

Ф0 = - 1, λ1 = λ2 = ... = λq = 0;

- вычисление fj

                                                   (38)

при

если

/fj/ < /β/,                                                                (39)

для j = 1, ..., q,

при β заданном, вычисляются оценки параметров скользящего среднего

                                                  (40)

В противном случае вычисляются новые значения λj, j = 1, ..., q

                                                     (41)

где , а матрицы Т составляются следующим образом

                                                (42)

                                              (43)

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

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

Максимальный порядок смешанной модели равен 4.

4.12. Выходные данные:

- вид и порядок модели временного ряда r, q, d;

- значения предварительных оценок параметров модели

Фi, i = 1, ..., r; qj, j = 1, ..., q.

5. АЛГОРИТМ ОЦЕНКИ ПАРАМЕТРОВ МОДЕЛИ

5.1. Входные данные:

- число значений временного ряда N;

- значения временного ряда {Pt}, t = 1, ..., N;

- вид и порядок модели временного ряда r, q, d;

- значения предварительных оценок параметров модели

Фi, i = 1, ..., r; qj, j = 1, ..., q;

- таблица процентных точек распределения c2;

- критерий сходимости γ;

- приращение d;

- уровень значимости c2 - распределения e.

5.2. Процедура определения случайной последовательности импульсов at в процессе вычислений производится несколько раз, поэтому она должна быть выделена отдельным блоком (блок А).

5.2.1. Вычисление случайной последовательности импульсов et производится по формуле

et = Pt - Ф1 Рt + 1 - ... - Фr Pt + r + q1 et + 1 + ... + qq et + q.                           (44)

5.2.2. Вычисление значений временного ряда Рt для t £ 0 производится по формуле

Pt = et + Ф1 Рt + 1 + ... + Фr Pt + r - q1 et + 1 - ... - qq et + q,                           (45)

при t = 0, - 1, ..., Т,

где Т = t, при котором /Pt/ £ 0,01; e-t = 0, t = 0, 1, ..., n.

5.2.3. Вычисление последовательности at производится по формуле

at = Pt - Ф1 Рt - 1 - ... - Фr Pt - r + q1 at - 1 + ... + qq at - q,                              (46)

при t = T, ..., 0, 1, ..., N - r.

a-t = 0 при t < 1 - T.

5.3. Вычисление последовательности at0 (блок А) производится для предварительных оценок параметров Фi, i = 1, ..., r; qj, j = 1, ..., q, полученных при идентификации модели.

5.4. Вычисление последовательности atб (блок А) с одним возмущенным параметром производится по формуле

gl¢ = gl + d,                                                           (47)

при l = 1, ..., (r + q),

где gl - параметр модели.

5.5. Вычисление производной yl, t производится по формуле

                                                       (48)

при t = Т, ..., N,

l = 1, ..., (r + q).

5.6. Определение приращения производится по формуле

                                                     (49)

при l = 1, ..., (r + q).

5.7. Вычисление нового значения параметра производится по формуле

                                                       (50)

5.8. Сходимость итеративной процедуры определяется из условия /gпре/ > /γ/ для l = 1, ..., (r + q) Все параметры заменяются на вновь полученные и процедура вычисления последовательности ats повторяется.

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

5.9. Определение доверительных интервалов по критерию наименьшей суммы квадратов

5.9.1. Вычисление последовательности at для полученных оценок параметров проводится в блоке А.

5.9.2. Вычисление минимальной суммы квадратов проводится по формуле

                                                        (51)

5.9.3. Определение доверительного интервала оценок параметров проводится по формуле

                                                   (52)

5.10. Выходные данные:

- количественные оценки параметров математической модели временного ряда

Фi, i = 1, ..., r; qj, j = 1, ..., q.

- доверительный интервал оценок наименьших квадратов Sдов.

6. АЛГОРИТМ ПРОВЕРКИ АДЕКВАТНОСТИ МОДЕЛИ

6.1. Входные данные:

- число наблюдений временного ряда N;

- значения временного ряда {Pt}, t = 1, ..., N;

- вид и порядок модели временного ряда r, q, d;

- значения параметров модели Фi, i = 1, ..., r; qj, j = 1, ..., q;

- величина, определяющая допустимое отклонение Q от соответствующего квантиля c2 - распределения a.

6.2. Вычисление автокорреляционной функции остаточных ошибок

6.2.1. Вычисление коэффициентов j производится по формуле

                                      (53)

при j = 1, ..., (r + q); j0 = 1;

6.2.2. Вычисление коэффициентов y производится по формуле

                                                    (54)

для j < 0 yj = 0;

для j > 0  = 0, y0 = 1.

Вычисления продолжаются до выполнения условия

yj £ 0,01.                                                                 (55)

6.2.3. Максимальная задержка автокорреляции равна индексу j при коэффициенте y, удовлетворяющему условию (55).

6.2.4. Вычисление последовательности остаточных ошибок at проводится в блоке А.

6.2.5. Вычисление автоковариаций остаточных ошибок производятся по формуле

                                                           (56)

6.2.6. Вычисление автокорреляций остаточных ошибок проводится по формуле

Ral = Cal / Ca0,                                                             (57)

при l = 0, 1, ..., кy.

6.3. Проверка модели временного ряда по статистике Q.

6.3.1. Вычисление статистики Q проводится по формуле

                                                              (58)

6.3.2. Определение числа степеней свободы проводится по формуле

U = кy - r - q.                                                              (59)

6.3.3. Проверка адекватности по величине Q производится по условию

Q < cl - a.

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

6.4. Проверка адекватности модели по верхней границе стандартной ошибки автокорреляции.

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

                                                             (60)

при l = 1, ..., кy.

Когда это условие выполняется, Na = Na + 1.

6.4.2. Сравнение полученного числа Na с допустимым числом выбросов Ng = кy / 3 проводится по условию Nа < Ng.

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

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

6.5. Выходные данные:

- подтверждение или неподтверждение гипотезы об адекватности модели по статистике Q;

- подтверждение или неподтверждение гипотезы об адекватности модели по верхней границе стандартной ошибки автокорреляции.

6.6. Блок-схемы алгоритмов построения математической модели временного ряда приведены в рекомендуемом приложении 2.

ПРИЛОЖЕНИЕ 1

Справочное

Примеры выбора вида математической модели

1. Модель «белый шум»

1.1. Спектр автокорреляционной функции приведен на черт. 1.

Черт. 1

1.2. Спектр частной автокорреляционной функции приведен на черт. 2.

Черт. 2

2. Модель авторегрессии 2-го порядка

2.1. Спектр автокорреляционной функции приведен на черт. 3.

Черт. 3

2.2. Спектр частной автокорреляционной функции приведен на черт. 4.

Черт. 4

3. Модель скользящего среднего 1-го порядка

3.1. Спектр автокорреляционной функции приведен на черт. 5.

Черт. 5

3.2. Спектр частной автокорреляционной функции приведен на черт. 6.

Черт. 6

4. Смешанная модель авторегрессии - скользящего среднего

4.1. Спектр автокорреляционной функции приведен на черт. 7

Черт. 7

4.2. Спектр частной автокорреляционной функции приведен на черт. 8.

Черт. 8

ПРИЛОЖЕНИЕ 2

Рекомендуемое

Блок-схемы алгоритмов построения математической модели временного ряда

1. Блок-схема алгоритма проверки ряда на стационарность приведена на черт. 1.

Черт. 1

2. Блок-схема алгоритма идентификации приведена на черт. 2.

Черт. 2

3. Блок-схема вычисления последовательности случайных импульсов at (блок А) приведена на черт. 3.

Черт. 3

4. Блок-схема алгоритма оценки параметров модели приведена на черт. 4.

Черт. 4

6. Блок-схема проверки адекватности модели приведена на черт. 5.

Черт. 5

 

 

ЛИСТ РЕГИСТРАЦИИ ИЗМЕНЕНИЙ

№ изм.

Номера страниц

Номер «Изв. об изм.»

Подпись

Дата

Срок введения изменения

Измененных

Замененных

Новых

Аннулированных

1

1

-

-

-

9071

 

22.05.84

01.07.84

2

1

-

-

-

9340

 

3.12.85

 

 

СОДЕРЖАНИЕ

 



© 2013 Ёшкин Кот :-)