CC BY
70
УДК 519.246.8
ОЦЕНКА СОСТОЯНИЯ ОБЪЕКТА ПО РЕГРЕССИОННЫМ ЗАВИСИМОСТЯМ ПРИ ПРОГНОЗИРОВАНИИ ВХОДНЫХ
ПАРАМЕТРОВ
© 2014 Ю.Е. Кувайскова, Д.С. Бубырь, В.Н. Клячкин
Ульяновский государственный технический университет
Поступила в редакцию 28.05.2014
Предложена методика оценки выходных параметров, характеризующих состояние технического объекта, на основе построения кусочно-линейных регрессионных зависимостей от входных характеристик объекта, прогнозируемых с использованием адаптивного динамического регрессионного моделирования на основе смешанной модели авторегрессии-скользящего среднего, комплекса авторегрессионных моделей с условной гетероскедастичностью, или последовательного применения этих подходов.
Ключевые слова: моделирование, технический объект, прогнозирование, кусочно-линейная регрессия, локальность модели
Моделирование по поступающим результатам измерений и прогнозирование состояния технического объекта по полученным моделям с формированием сигнала предупреждения о возможной аномальной ситуации, при которой характеристики процесса выходят за допустимые пределы, является важной задачей [1-3]. Решение задачи прогнозирования состояния объекта по регрессионным моделям предполагает, что будущие значения входных параметров (регрес-соров) известны. Однако на практике при исследовании технических объектов значения входных и выходных параметров определяются по результатам измерений в режиме реального времени. Таким образом, при прогнозировании состояния объекта значения входных параметров неизвестны. Для решения этой задачи предлагается алгоритм прогнозирования состояния технического объекта по кусочно-линейным регрессионным зависимостям выходных характеристик объекта с использованием спрогнозированных значений входных параметров, получаемых по моделям временных рядов.
Рассмотрим технический объект У, состояние которого зависит от некоторых входных параметров, значения которых регистрируются через определенные промежутки времени и образуют временные ряды Х1,Х2,...,Хт, где
Кувайскова Юлия Евгеньевна, кандидат технических наук, доцент кафедры прикладной математики и информатики. E-mail: u.kuvaiskova@mail.ru Бубырь Дмитрий Сергеевич, аспирант Клячкин Владимир Николаевич, доктор технических наук, профессор кафедры прикладной математики и информатики. E-mail: v_kl@mail.ru
т - количество параметров. Алгоритм моделирования и прогнозирования состояния технического объекта состоит из следующих этапов:
- построение моделей входных параметров объекта, образующих систему временных рядов Х1,Х2,...,Хт, на основе методологии адаптивного динамического регрессионного моделирования [4-6];
- прогнозирование по полученным моделям временных рядов значений входных параметров исследуемого объекта;
- построение регрессионной зависимости ((Х1, Х2,..., Хт ) , характеризующей связь между состоянием объекта У и его входных параметров Х] = 1, т;
- использование полученной регрессии и спрогнозированных значений входных параметров для прогнозирования состояния технического объекта.
Согласно методологии динамического регрессионного моделирования временной ряд представляется в виде суммы составляющих:
У® = /(0 + ) + ) + е($) , (1)
где у(() - значения временного ряда, фиксируемые в моменты времени I = 1, 2,..., п; ((^)
- неслучайная (долговременная) функция тренда; е (^) - неслучайная периодическая функция; у/(() - случайная с элементами регулярности функция; е(1) - нерегулярная компонента (случайная величина, ошибка).
При выявлении значимой регулярности процесса методами фрактального и мультифрак-тального анализов [7] выделяется функция тренда f (t). Далее остатки исследуются на существование периодических компонент при помощи гармонического анализа и других методов. Остаточные колебания после выделения неслучайных функций сглаживаются случайной с элементами регулярности функцией iy(t) . В качестве данной функции могут выступать либо авторегрессионная модель подходящего порядка, либо смешанная модель авторегрессии-скользящего среднего (АРСС), либо одна из комплекса авторегрессионных моделей с условной гетероскедастичностью (ARCH) [8], либо последовательное применение этих подходов.
После структурно-параметрической идентификации модели временного ряда проверяется соблюдение условий применения регрессионного анализа: постоянство дисперсии, независимость регрессоров, нормальность распределения ошибок, нулевое значение математического ожидания ошибок, независимость ошибок. При нарушении этих условий используются специ-
альные методы адаптации [4]. Если основные условия соблюдаются, построенная комплексная модель временного ряда может быть использована для прогнозирования. После получения предсказанных значений входных параметров по моделям временных рядов, строятся регрессионные зависимости выходных характеристик объекта в зависимости от входных (регрессоров).
При неоднородности физических свойств объекта на области значений регрессоров, «глобальные» модели часто обладают недостаточно высокой точностью. В этом случае предлагается применять принцип «кусочности» или локальности модели, то есть вариации её параметров по области значений регрессоров. В качестве прогнозирующей зависимости можно использовать кусочно-линейную регрессию с разрывом по отклику и с добавлением элемента авторегрессии для учёта значения отклика в предыдущие дни. При этом свободный коэффициент в модель предполагается не включать. В этом случае регрессионная модель представляется в виде:
Y(t) = (РX, + AiX2 +... + ртХXm + p^Y(t - 1))(Y(t) < c) + (Д2X1 + P22X2 + ... + Pm2Xm + Pm+1,2Y(t - 1))(Y(t) > c)
(2)
где т - количество регрессоров модели; к - порядок авторегрессии; с - точка разрыва; (У^) < с), (У (¿) > с) - логические выражения, принимающие значения: 1 - если истинно, 0 -если ложно; У(1-1). У((-к) - значение отклика в предыдущие моменты времени. Фактически данная модель состоит из двух регрессий, каждая из которых применяется для определённого уровня отклика.
Одним из ключевых моментов при использовании кусочных регрессий является способ отнесения текущих данных к нужному «куску». Предлагается следующий алгоритм этой процедуры. Вначале вычисляется прогноз по обеим регрессиям. Затем полученные результаты сравниваются со значением в точке разрыва. Если прогноз по первому куску меньше либо равен значению в точке разрыва, а прогноз по второму куску - больше, то за итоговый результат принимается среднее арифметическое данных чисел. Если хотя бы для одного куска наблюдается «нарушение», то итоговым значением считается прогноз по тому куску, для которого выполняется неравенство. Если же условия нарушаются для обоих случаев, то также вычисляется среднее арифметическое.
Например, имеется кусочно-линейная регрессия:
У = 11 • (У < с) + 72 • (У > с), (3)
где У1 = (6Ш + Ьи • X + ■•• + Ьт1 • Хт) ;
У2 = (Ь02 + Ь12 • Х1 + ■■■ + Ьт2 • Хт ) ; С - точка
разрыва.
Пусть на основе регрессий У1 и У2 получены прогнозы У и У2, причём У1 < с и У2 < с. В данном случае наблюдается нарушение для второй регрессии, поэтому за итоговый прогноз берётся У1. По полученной модели кусочно-линейной регрессии с добавлением элемента авторегрессии и спрогнозированным значениям входных параметров прогнозируется будущее состояние технического объекта с целью формирования сигнала о возможной аварийной ситуации. Эффективность разработанной методики прогнозирования состояния технического объекта иллюстрируется на примере моделирования выходной характеристики объекта У и множества входных параметров, образующих временные ряды X, Х2,..., Х7 по результатам 35 наблюдений. Для системы временных рядов (входных параметров объекта) согласно методологии адаптивного динамического регрессионного моделирования построены следующие модели:
Х (0 = 0,47715 • Х (* -1) + 0,60453 • Х (* ~ 2),
Х2 (г) = 0,89898 • Х2 (г -1) + 0,09956 • Х2 (г - 2), X (г) = 0,79444 • Х2 (г -1) + 0,17374 • Х3 (г - 2), X (г) = 0,57316 • Х4 (г -1) + 0,42619 • Х4 (г - 2), X (г) = 0,87087 • Х5 (г -1) + 0,12289 • Х5 (г - 2) , Х6 (г) = 0,57535 • Х6 (г -1) + 0,41826 • Х6 (г - 2),
С использованием полученной регрессии и спрогнозированных значений входных параметров X, Х2,..., Х7 выполнен прогноз выходной характеристики объекта У на четыре наблюдения. На рис. 1 показан смоделированный временной ряд, а на рис. 2 - соответствующий прогноз на 4 измерения. Сплошной линией показаны опытные данные, штриховой - расчетные.
1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35
Рис. 1. Моделирование на основе КЛР-АР(1)
В табл. 1 представлены результаты сравнения кусочно-линейной модели с ЛР и ЛР-АР(1) по точности прогнозирования. В первом столбце таблицы представлены исходные значения отклика У; в следующих столбцах для каждой модели представлены значения спрогнозированных значений отклика У (прогноз), а также ошибка прогноза в процентах, вычисляемая по формуле (7).
I У - УI
е = --1 -100% , (7)
Х7 (г) = 0,90417 • Х7 (г -1) + 0,09523 • Х7 (г - 2) .
По полученным моделям временных рядов выполнен прогноз на четыре наблюдения. Для выходного параметра объекта У построена кусочно-линейная регрессия с авторегрессией 1 -го порядка (КЛР-АР(1)):
(4)
Рис. 2. Прогнозирование состояния объекта на основе КЛР-АР(1)
Для сравнительного анализа эффективности предлагаемой методики прогнозирования выходных характеристик объекта по той же выборке были построены модели множественной линейной регрессии (ЛР) и множественной линейной регрессии с добавлением элемента авторегрессии 1-го порядка (ЛР-АР(1)):
(5)
(6)
где У - исходное значение отклика; У- прогнозное значение отклика.
По данным таблицы следует, что при использовании модели кусочно-линейной регрессии с авторегрессией 1 -го порядка точность прогнозирования существенно повышается по сравнению с моделями обычной множественной линейной регрессии и множественной линейной регрессии с добавлением элемента авторегрессии 1 -го порядка.
у(г) = (-0,002228 • Х, - 0,034335 • Х2 + 0,028101 • Х3 - 0,295028 • Х4 + 0,505128 • Х5 + + 0,031676 • Х6 + 0,433677 • Х7 + 0,787707 • У (г -1)) • (У (г) < 1,0109) + (-0,028096 • Х --0,049394 • Х2 + 0,046701 • Х3 - 0,208433 • Х4 +1,182003 • Х5 + 0,010884• Х6 + + 0,520216 • Х7 - 0,110927 • У (г -1)) • (У (г) > 1,0109)
У (г) = -13,604 + 0,0112 • Х1 - 0,4476 • Х2 + 0,0422 • Х3 - 0,3325 • Х4 +1,1932 • Х5
+ 0,0216 • Х6 + 4,5121 • Х7
У (г) = -19,6682 - 0,0011 • Х - 0,6002 • Х2 + 0,0367 • Х3 - 0,2679 • Х4 + Х5 + 0,0229 • Х6 + 6,0737 • Х7 + 0,3024 • У (г -1)
Таблица 1. Сравнение моделей по точности прогнозирования
Исходные значения отклика Y Модель регрессии
КЛР-АР(1) ЛР ЛР-АР(1)
прогноз ошибка, % прогноз ошибка, % прогноз ошибка, %
0,99 1,02 2,87 0,97 1,75 1,27 28,77
1,06 1,07 0,84 0,60 43,82 0,81 24,05
1,03 1,06 2,76 0,58 43,48 0,64 37,44
1,03 1,05 2,14 0,61 40,75 0,63 39,24
Выводы: использование кусочно-линейных моделей регрессии позволяет повысить точность прогнозирования состояния технического объекта. Безусловным преимуществом кусочно-линейных регрессий относительно обычных линейных регрессий является учёт неоднородности свойств объекта за счёт своей «кусочности» и вариации уровней отклика. Это позволяет более точно моделировать и прогнозировать процесс, в особенности, когда поступающие данные отличаются значительной нестабильностью значений.
Работа выполнена в рамках задания Минобрнау-ки России №2014/232.
СПИСОК ЛИТЕРАТУРЫ:
1. Клячкин, В.Н. Информационно-математическая система раннего предупреждения об аварийной ситуации / В.Н. Клячкин, Ю.Е. Кувайскова, А.А. Алешина, ЮА. Кравцов // Известия Самарского научного центра РАН. 2013. №4(4). С. 919-923.
2. Клячкин, В.Н. Моделирование вибраций гидроагрегата на основе адаптивных динамических регрессий / В.Н. Клячкин, Ю.Е. Кувайскова, АА. Алешина // Автоматизация и современные технологии. 2014. №1. С. 30-34.
3. Клячкин, В.Н. Диагностика состояния объекта по наличию неслучайных структур на контрольных картах / В.Н. Клячкин, Ю.А. Кравцов // Приборы и системы. Управление, контроль, диагностика. 2013. №5. С.44-50.
4. Валеев С.Г. Регрессионное моделирование при обработке наблюдений. М. : Наука, 1991. 272 с.
5. Валеев, С.Г. Особенности построения регрессионных моделей при многомерном контроле технологического процесса / С.Г. Валеев, В.Н. Клячкин //
Радиоэлектроника. Информатика. Управление. 2002. №1. С. 48-52.
6. Валеев, С.Г. Программное обеспечение обработки временных рядов техногенных характеристик / С.Г. Валеев, Ю.Е. Кувайскова // Обозрение прикладной и промышленной математики. 2009. Т. 16, выпуск 6. С. 1037-1038.
7. Валеев, С.Г. Применение мультифрактального анализа при описании временных рядов в технике и экономике / С.Г. Валеев, Ю.Е. Кувайскова, С.А. Гу-байдуллина // Вестник Ульяновского государственного технического университета. 2008. №2. C. 2327.
8. Валеев, С.Г. Использование ARCH-структур и фильтра Калмана для моделирования динамики технико-экономических показателей / С.Г. Валеев, Ю.Е. Кувайскова // Вестник УлГТУ. 2007. №2. C. 29-33.
9. Кувайскова, Ю.Е. Прогнозирование состояния технического объекта на основе мониторинга его параметров / Ю.Е. Кувайскова, В.Н. Клячкин, Д.С. Бу-бырь // XII Всероссийское совещание по проблемам управления. Институт проблем управления им. Трапезникова РАН [Электронный ресурс] URL: http://vspu2014.ipu.ru/node/2940 (дата обращения: 16.05.2014).
10. Крашенинников, В.Р. Кусочно-квадратичное моделирование регрессионных зависимостей при оценке качества / В.Р. Крашенинников, Д.С. Бубырь // Мат-лы 3-й науч.-практ. internet-конференции «Междисциплинарные исследования в области математического моделирования и информатики» 2021 февраля 2014 г. / отв. ред. Ю.С. Нагорнов. -Ульяновск: SIMJET, 2014. С. 233-236.
11. Васильев, К.К. Статистический анализ многомерных изображений / К.К. Васильев, В.Р. Крашенинников. - Ульяновск: УлГТУ, 2007. 170 с.
FORECASTING OF THE TECHNICAL OBJECT STATE BASED ON THE PIECEWISE LINEAR REGRESSION MODELS
© 2014 Yu.E. Kuvayskova, D.S. Bubyr, V.N. Klyachkin
Ulyanovsk State Technical University
In this article the method of forecasting output parameters characterizing the state of the technical object, based on the construction of piecewise linear regression dependencies on input characteristics of the object.
Key words: modeling, technical object, forecasting, piecewise linear regression, locality of the model
Kuvayskova, Candidate of Technical Sciences, Associate Professor at the Department "AppliedMathematics and Computing Science". E-mail: u.kuvaiskova@mail.ru; Dmitriy Bubyr, Post-Graduate Student; Vladimir Klyachkin, Doctor of Technical Sciences, Professor at the Department "AppliedMathematics and Computing Science". E-mail: v_kl@mail.ru.
