CC BY
94
Математическое моделирование. Оптимальное управление Вестник Нижегородского университета им. Н.И. Лобачевского, 2012, № 5 (2), с. 28-32
УДК 517.587+519.651.3+534.4
ПРОГНОЗ ВРЕМЕННЫХ ПОСЛЕДОВАТЕЛЬНОСТЕЙ С ИСПОЛЬЗОВАНИЕМ ОБОБЩЁННОГО СПЕКТРАЛЬНО-АНАЛИТИЧЕСКОГО МЕТОДА
© 2012 г. А.К. Бритенков1, Ф.Ф. Дедус2
1 Нижегородский госуниверситет им. Н.И. Лобачевского 2Институт математических проблем биологии РАН, Пущино
jkd@ya.ru
Поступила в редакцию 10.09.2012
На основе обзора регрессионных моделей рассмотрено применение обобщённого спектрально-аналитического метода для прогноза временных рядов, эффективного в задачах описания и распознавания сигналов. Предлагается использование ортогональных базисов на основе классических полиномов непрерывного аргумента в задачах прогноза как частного случая экстраполяции. Проанализировано применение дополнительных данных для улучшения достоверности прогноза.
Ключевые слова: прогноз временных рядов, регрессия, обобщённый спектрально-аналитический метод, классические ортогональные полиномы, метод наименьших квадратов, экстраполяция, системный анализ.
Краткий обзор методов прогнозирования
Методы прогнозирования временных рядов
[1] условно можно разделить на несколько направлений [2]. К первому направлению относятся аппроксимативные методы прогноза [3], например динамический анализ, модели экспоненциального сглаживания (ЛШМЛ), метод декомпозиции временных рядов (тренд, сезонность, случайность, выделение циклической компоненты), выделение скрытых периодичностей (спектральный анализ, анализ периодограммы) [4]. Другое направление, применяя комплексный подход к решению задач прогноза, предполагает использование нейросетевых методов, метода парадигм, регрессионного анализа (ЛЯ(р), ЛКМЛ(р, д) и т.д.) и статистических оценок [3]. В описательном анализе [5] используются выборочные характеристики временных рядов: среднее, дисперсия, автокорреляционная функция, частная автокорреляционная функция, распределения выборочных характеристик, стационарность временного ряда [6]. Прогнозирование на основе пропорциональных зависимостей и комбинированные методы [7] могут включать в себя также экспертные оценки и применение методов множественной регрессии
[2], анализ распределенных лагов (сетевых задержек, когда на прогноз оказывают влияние более ранние события, чем происходящие непосредственно накануне прогноза) и оценки адекватности модели прогноза. В последнем случае результаты, полученные с помощью модели,
сравниваются с фактическими значениями через какой-то промежуток времени, когда те появляются [5]. Недостаток такого подхода состоит в том, что проверка «беспристрастности» модели может занять много времени, так как по-настоящему проверить модель можно только на продолжительном временном отрезке [8].
Параметры прогноза временных рядов
Регрессионные модели прогноза используют описание независимой переменной у с помощью зависимой переменной х [2]. Прогноз в дальнейшем - результат подстановки новых значений х в найденное функциональное уравнение у = /(х). В качестве такой функции часто используются аппроксимирующие функции. Аппроксимация с помощью классических ортогональных функций [9], где расчёт коэффициентов предполагается по методу наименьших квадратов (МНК), наиболее удовлетворяет перечисленным критериям [10, 11]. В регрессионных моделях, где функциональная зависимость также ищется в виде функции, наилучшим образом приближающей искомую кривую с помощью МНК, вводится понятие функции невязки
1 м
^)=тЕ(^-4)2, а)
1 к=1
где М - объём выборки. Поиск минимума невязки отклонений реально наблюдаемых У от их оценок 7 (имеются в виду оценки с помощью
кривой, претендующей на то, чтобы представлять искомую регрессионную зависимость) основан на том, что согласно МНК сумма в выражении (1) принимает минимальное значение именно тогда, когда У = у(хі,х2,...хл). Временная последовательность s(t), описываемая ортогональным рядом вида
N
(2)
где ф, - ортогональный базис (полученный модификацией классических ортогональных полиномов, которые являются частным случаем решения гипергеометрического уравнения [12, 13]), как было отмечено выше, так же подразумевает невязку вида (1) для вычисления А,:
»ц—2 <•
Ф , | |фг ^^)p(t)dt .
(3)
Таким образом, под параметрами прогноза подразумеваются длина выборки М, глубина прогноза (количество шагов), порядок модели (для АЯ - р), ошибка прогноза [14].
550
500
450
400
350
300
250
А
Г V* ,/
/ 1
\
/ і \
V, у
2115
2125 2135
214-5 2155
21(55
Рис. 1. Прогноз энергопотребления авторегрессионной моделью (АЩр), р = 480), длина выборки N = 2100 точек (48 отсчетов в сутки), прогноз составлен на 6 шагов вперед (3 часа)
Прогноз временных рядов ортогональными базисами непрерывного аргумента
Применение обобщенного спектрально-аналитического метода (ОСАМ) для прогнозирования потоков данных, определение задачи прогноза как частого случая экстраполяции и
фильтрации помех исходит из того, что, как показано выше, МНК, применяемый в регрессионных методах прогноза, идентичен МНК при вычислении коэффициентов ряда Фурье [11, 15]. Предполагается, что разработка методов с использованием свойств ортогональных базисов на основе классических полиномов непрерывного аргумента (в задачах, где тригонометрические ряды Фурье менее эффективны, например [16]) позволяет расширить возможности математического аппарата ОСАМ, разработанного во второй половине XX века [17].
Проводя аналогию регрессии с физическими моделями, рассмотрим систему, содержащую N отражательных поверхностей для акустических волн [18]. Математически такая система напоминает регрессионную модель. Действительно, в такой системе исходный звуковой сигнал S(?) преобразуется в S (?) так, что в общем случае S (?) можно представить суммой дискретных задержанных импульсов и ослабленного после прохождения через среду исходного сигнала
N N
S * (?) = (V (?) + X (? -т,) = £ (? - т,), (4)
,=1 , = 0
где а, - коэффициенты отражения; т, - задержка во времени одного отражения, N - число пришедших в приемник отражений (первичных, вторичных и т.д.). Для сигнала, достигшего приемника без отражений, т0 = 0. Решение с помощью обобщенного спектрально-аналитического метода обратной задачи о распространении звука в условиях дискретных отражений [19] показывает возможность использования подобной технологии как компонента регрессионного моделирования.
Определение параметров прогноза квазипериодических сигналов с аддитивным шумом
В каждом реальном сигнале присутствует шум. Особенно это относится к задачам обработки временных рядов, когда шум присутствует естественным образом в самой системе или механизме измерений. В финансовых приложениях, статистических параметрах социальных задач данные зашумлены особенно сильно. Например, данные о совершении сделок или котировки ценных бумаг и другие показатели могут регистрироваться с запозданием, причем в различных случаях - с разным. Пропуск значений или неполная информация также рассматривается как шум. Использование в таких случаях на коротком отрезке среднего значения вместо отсутствующего отсчета приводит к зашумле-
, 0
нию данных. Таким образом, для успешного прогнозирования необходима эффективная обработка входных данных, в частности минимизация случайных флуктуаций и шума.
Фильтрации шума посвящено большое количество работ. Условно можно разделить методы шумоподавления на спектральные [18, 20], и статистические [21], не пренебрегая, разумеется, и их комбинацией [20, 22]. В ряде работ спектральные методы используют вейвлет-анализ, как, например, в работах [20, 23]. В этих случаях также возможно применение ортогональных базисов на основе классических полиномов непрерывного аргумента как базисных вейвлетов. Примером такого вейвлета является функция Эрмита (рис. 2):
н„(>) = —— X (-1)" гТ (5)
2Г(-п) £ I 2 ) к
где Г - гамма-функция, п - порядок функции Эрмита.
Рис. 3. Прогноз-экстраполяция потребления электроэнергии предприятия (кВт/ч) функциями Лежандра (Рп(?), п = 100, число шагов прогноза более 1000, число «обучающих точек» р = 675)
получение статистических соотношений на основе коэффициентов разложения (2) оптимизирует вычислительные процедуры [10, 17].
Алгоритмы оптимизации прогноза и дополнительные данные
Прогнозирование одномерных сигналов и временных данных можно существенно улучшить, используя дополнительные данные, например корреляцию прогнозируемых сигналов с другими параметрами, испытывающими колебания на той же временной оси. В случае с энергопотреблением предприятия такими данными могут быть план производства, температура окружающей среды и собственно время (время суток, день недели, время года). В этих графиках спектральный анализ позволяет выделить характерные периоды, что дает дополнительные данные для уточнения прогноза. Синергетический подход показывает, что колебания прогнозируемой величины и дополнительных данных представляют с разных сторон исследуемую единую систему [24], поэтому корреляция между ними очевидна (рис. 4, 5).
Рис. 2. Функция Эрмита (Ип(?), п = 4, п = 9)
Идея использования классических ортогональных полиномов для прогноза (экстраполяции) основана на том, что фактически такая экстраполяция является частным случаем выделения сигнала на фоне помех (рис. 3).
Однако простое продление аппроксимирующей кривой за пределы доступных данных не всегда приводит к успешному предсказанию. С другой стороны, использование ортогональных базисов на основе полиномов непрерывного аргумента из числа классических в задаче со статистическими методами также является логичным ввиду математического аппарата ОСАМ:
Рис. 4. Пример потребления электроэнергии предприятия (кВт/ч) в течение 3 дней
Рис. 5. Циклическое представление времени на оси отсчетов данных электропотребления предприятия в течение 3 дней
З1
Обсуждение результатов
Актуальность классических ортогональных функций и обобщенных рядов Фурье подтверждается использованием их в описании случайных процессов [20], решении дифференциальных уравнений [25], моделировании оптических негармонических импульсов [26], эмпирических базисах в гидрофизике [16], развитии вейвлет-анализа [27] и обобщенного спектрально-аналитического метода [5, 10, 14, 15].
Исходя из природы используемых данных (котировок ценных бумаг, показателей электропотребления предприятия), анализ дополнительных условий на входные данные, а также конкретизация «физики» колебаний показателей временных последовательностей для составления прогноза является важнейшим условием повышения точности прогноза. Математически это определение таких критериев решения, как ошибка аппроксимации; точность предсказания и выбор параметров прогнозирования; оценка статистических характеристик исходных данных, фильтрация шума, корреляция с параллельными потоками данных. На следующем этапе реализуется оптимизация способа аппроксимации дискретных данных разложением по классическим ортогональным базисам. Однако построение прогноза временного ряда, использующее экстраполяцию аппроксимирующей функцией, требует комбинации с другими методами (например, регрессионным анализом), поскольку в соответствии со скалярным произведением (3) нужны ненулевые данные /(?) (или модели) для вычисления коэффициентов разложения на участке прогноза.
Важность системного подхода продиктована требованием целостного решения задач с элементами выбора, принятия решений и поиска наилучшей стратегии [7]. Математические методы описания и анализа сложных экономических, социальных, экологических систем показали, что кроме объективных законов в функционировании смешанных систем существенную роль играют субъективные представления [8]. Поскольку в таком случае не представляется возможным формализовать (или количественно описать) реальную систему и, таким образом, построить адекватную математическую модель, только системный подход и интеграция различных методов (в т.ч. эвристических приемов) позволяет по-настоящему эффективно решать подобные задачи [28].
Авторы выражают глубокую признательность д.ф.-м.н. Н.Е. Никитиной и к.ф.-м.н.
А.Н. Панкратову за внимание к работе и высказанные замечания.
Работа выполнена при поддержке гранта РФФИ 1001-00609.
Список литературы
1. Вучков И., Бояджиева Л., Солаков Е. Прикладной линейный регрессионный анализ. М.: Финансы и статистика, 1987. 239 с.
2. Robert K. O., Loren E. Applied Time Series Analysis. N. Y. - Chichester - Brisbane - Toronto: A Wiley-Interscience Publication John Wiley and Sons, 1978. 428 p.
3. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. Множественная регрессия. М.: Диалектика, 2007. 912 с.
4. Дженкинс Г., Ватте Д. Спектральный анализ и его приложения. М.: Мир. Ч. 1. Вып. 1. 1971. 316 с.; Вып. 2. 1972. 288 с.
5. Britenkov A. Forecasting of Time Series by Classical Ortohonal Functions in Regression Models // Proceеdings of 8th Open German-Russian workshop «Pattern Recognition and Image Understanding OGRW-8-2011». N.Novgorod: Nizhny Novgorod Lobachevsky State University, 2011. P. 28-30.
6. Пугачев В.С. Теория вероятностей и математическая статистика. М.: Наука, 1979. 496 с.
7. Бритенков А.К., Артюхин И.В. Описание, анализ и прогнозирование котировок ценных бумаг с помощью обобщенного спектрально-аналитического метода для динамического управления портфелем ценных бумаг// Докл. X Всероссийской конференции ММР0-10 «Математические методы распознавания образов». М.: Изд-во ВЦ РАН, 2001. C. 165.
8. Бритенков А.К., Воинов Б.С., Сорокин Ю.М. Системный анализ (программа курса и компьютерный практикум). Экосоциум. Проблемно-ориентированная программа подготовки для специальности «Информационные системы (радиофизика, телекоммуникации)»/ Под общ. ред. Ю. М. Сорокина. Н. Новгород: Изд. ННГУ, 1997. 63 с.
9. Никифоров А.Ф., Уваров В.Б. Специальные функции математической физики М.: Наука, 1984. 344 с.
10. Britenkov A.K., Pankratov A.N. Stable Algorithms of Adaptive Approximation for Acoustic Signals Description by Orthohonal Polynomials // Physics of Wave Phenomena. 2004. V. 12 (3). P. 168-174.
11. Линник Ю.В. Метод наименьших квадратов и основы теории обработки наблюдений. М.-Л.: Физ-матгиз, 1962. 352 с.
12. Справочник по специальным функциям: Пер. с англ. / Под ред. М. Абрамовица, И. Стиган. М.: Наука, 1979. 832 c.
13. Kazuhiko Aomoto, Michitake Kita. Theory of Hypergeometric Functions / Transl. by Kenji Iohara. Springer, 2011. Vol. 305. 317 p.
14. Бритенков А. К., Артюхин И. В. Динамическое управление портфелем ценных бумаг на основе оптимального прогноза вектора наблюдений // Докл.
X Всерос. конф. ММР0-10 «Математические методы распознавания образов». М.: ВЦ РАН, 2001. C. 167.
15. Pankratov A. Long-Range Monitoring of Fuel Consumption of Car Based on Generalized Spectral-Analytical Method // Book of Proceedings «Driver Car Interaction & Interface 2010». Prague, 2010. P. 49-50.
16. Вировлянский А. Л., Казарова А. Ю., Любавин Л.Я., Стромков А.А. Эмпирические ортогональные функции в гидроакустических экспериментах. Методы и средства измерения физических величин // Тез. докл. 3-й Всероссийской науч.-техн. конф. Н.Новго-род: Изд. ИПФ РАН, 1998. С. 33.
17. Дедус Ф.Ф., Махортых С.А., Устинин М.Н., Дедус А.Ф. Обобщенный спектрально-аналитический метод обработки информационных массивов. Задачи анализа изображений и распознавания образов / Под общ. ред. Ф. Ф. Дедуса. М.: Машиностроение, 1999. 357 с.
18. Зверев В.А., Стромков А.А. Выделение сигналов из помехи численными методами. Н. Новгород: ИПФ РАН, 2001. 188 с.
19. Бритенков А.К., Панкратов А.Н. Подавление мультипликативных помех с помощью обобщенного спектрально-аналитического метода в условиях дискретных отражений // Вестник ННГУ им. Н. И. Лобачевского. Серия Радиофизика. 2006. Вып. 1(4).
20. Рытов С. М. Введение в статистическую радиофизику. Часть 1. Случайные процессы. М.: Наука, 1976. 491 с.
21. Grafov B.P., Grafova I.B. Theory of the Wavelet Analysis for Electrochemical Noise by Use of Laguerre Functions // Electrochemistry Communications. 2000.
(2). P. 386-389.
22. Уидроу Б., Стирнз С. Адаптивная обработка сигналов: Пер. с англ. М.: Радио и связь, 1989. 440 с.
23. Шварцбург А.Б. Оптика нестационарных сред // Успехи физических наук. 2005. Т. 175. №8. С. 833-861.
24. Shen J. Stable and Efficient Spectral Methods in Unbounded Domains Using Laguerre Functions // SIAM Journal on Numerical Analysis. 2000. V. 38 (4). P. 1113-1133.
25. Лоскутов А.Ю., Михайлов А.С. Введение в синергетику: Учеб. руководство. М.: Наука, 1990. 272 с.
26. Холод В.П. Новая вейвлет-образующая функ-
ция // [Электронный ресурс «Scientific World» (дата обращения: 20.02.2012)]. - Режим доступа:
http://www.sworld.com.ua/konfer24/802.htm
27. Астафьева Н. М. Вейвлет-анализ: основы теории и примеры применения // Успехи физических наук. 1998. Т. 166. №11. С. 1145-1170.
28. Бритенков А. К., Бугров В. Н. Экспертная и вычислительно-поисковая системы синтеза новых решений. Общие положения // В кн. «Информационные технологии и системы» (в 2-х кн.) / Под общ. ред. Б.С. Воинова. Н.Новгород: Изд-во ННГУ, 2001. Кн. 1. С. 231-241.
TIME SERIES FORECAST USING GENERALIZED SPECTRAL-ANALYTICAL METHOD
A.K. Britenkov, F.F. Dedus
Based on a review of regression models, the application of the generalized spectral-analytical method for time series forecasting is considered. This prediction is effective in the tasks of signal description and recognition. Classical orthogonal polynomials of a continuous variable are proposed to be used as a particular case of extrapolation. The use of additional data has been analyzed to improve forecast accuracy.
Keywords: time series forecasting, regression, generalized spectral-analytical method, classical orthogonal polynomials, method of least squares, extrapolation, systems analysis.
