CC BY
60
УДК 519.216.2
И.Ю. Кузнецова
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЭНЕРГОПОТРЕБЛЕНИЯ ПРИМЕНИТЕЛЬНО К ВУЗУ
Описана модель потребления электроэнергии на основе стохастического дифференциального уравнения. Приведены результаты численных расчетов по прогнозированию электропотребления на 2010 год на основании данных по потреблению электроэнергии общежитием №1 студенческого городка ТТИ ЮФУ за 2009 год. Применимость полученной модели к прогнозированию потребления электроэнергии оценена путем сравнения расчетных и реальных значений по двум показателям: средний процент ошибки и коэффициент несоответствия Теша.
Стохастические дифференциальные уравнения; численные методы; потребление.
I.U. Kuznetsova
MATHEMATICAL MODEL OF ENERGY CONSUMPTION FOR UNIVERSITY
The model of electricity consumption based on the stochastic differential equation is described in the article. Results of numerical calculations of energy consumption prognostication in 2010 based on data of energy consumption of hostel №1 of TTI SFU campus in 2009 are given in the article. The applicability of _ forecasting model of electricity consumption is estimated by comparing the calculated and actual values for the two indicators: the mean percentage error and Theil's variance coefficient.
Stochastic differential equations; numerical methods; consumption.
Моделирование спроса на электроэнергию является сравнительно новым направлением исследования в современной российской экономике, но представляется особенно актуальным в связи с реформированием электроэнергетики.
При моделировании процесса потребления электроэнергии необходимо учитывать три типа факторов: периодические (циклические), функциональные и слу-.
Периодические зависимости наиболее прогнозируемы и существенно влияют на потребление электроэнергии, например, продолжительность светового дня. При прогнозировании потребления электроэнергии в вузе огромное значение имеет , .
Функциональные факторы оказывают меньшее влияние на прогнозирование, чем циклические зависимости. К группе функциональных факторов можно отнести температуру воздуха, заселенность комнат и др.
Случайные тенденции составляют третью группу прогнозных факторов. Отличительной особенностью этой группы является то, что хотя их доля в процессе невелика, но амплитуда отклонений может быть довольно значительна.
Учитывая перечисленные выше факторы, была построена следующая модель :
jdP = (A(t) - aP(t ))dt + cr(t )dz,
I P (0) = po, (1)
где P(t) - функция потребления электроэнергии, p0 - значение потребления в начальный момент времени, параметр а представляет собой скорость возврата к среднему значению, функция A(t) учитывает сезонные факторы, a <j(t) - неопределенность, связанную с изменением потребления; z (t) — винеровский процесс.
Слагаемое c(t)dz позволяет учитывать зависимость изменения потребления : ( ), -ловеческий фактор (студент любит или нет яркое освещение) и т.д. Эти колебания
, , , -
.
Скорость возврата к среднему а можно статистически оценить с помощью регрессионного анализа, используя исторические данные [1]. Регрессионное уравнение для оценки параметра а имеет вид
Тогда a = -ln(1 + b).
Функцию неопределенности c(t) можно понимать как стандартное отклонение потребляемой энергии за год.
Для стохастического дифференциального уравнения (1) невозможно аналитически получить решение в терминах простых процессов. Метод Эйлера-Маруяма и метод Милштейна являются наиболее распространенными методами численного решения стохастических дифференциальных уравнений. Для уравне-(1) - , -онная часть c(t)dz не зависит от P(t).
Получим численное решение стохастического дифференциального уравнения (1) - .
Разобьем этот период времени на n равных интервалов
Переходя к конечным разностям из (1) и определения стохастического интеграла Ито, имеем
Отличие от случая детерминированного обыкновенного дифференциального уравнения заключается в том, что каждая последовательность {P0,...,PnJ, полученная методом Эйлера-Мару ямы, является приближенной реализацией стохастического процесса Pit), зависящего от выбора нормально распределенных случай-
i2)
0 = to < t1 < t2 < ... < tn = T ,
T
длина которых равна h = —. Введем следующие обозначения:
n
А zi = zi+1 - zi, i = 1 n.
где
ных величин ^ . Так как ) - стохастический процесс, то каждые его реализации будут различаться, а следовательно, будет отличаться и решение.
ных величин
Определение сходимости сходно с понятием сходимости решений обыкновенных дифференциальных уравнений, различия вызваны тем, что решение СДУ является случайным процессом и каждая вычисленная траектория является лишь одной возможной реализацией этого процесса. Каждая вычисленная траектория
решения P i = 1, n, полученная, например, методом Эйлера-Маруямы, дает случайное значение в точке T, таким образом, Pn также является случайной величиной. Разница между значениями в момент времени T, e(T) = P(T) — Pn, также является случайной величиной.
Говорят, что решение СДУ сходится сильно с порядком т, если математическое ожидание ошибки имеет m-й порядок от шага, т.е. для любого момента вре-T,
E {| P(T) — P, |}= O (hm)
h. -
терии сходимости для обыкновенных дифференциальных уравнений, и сводится к , . Определим порядок сходимости метода Эйлера-Маруямы, для этого разложим P(tM) в окрестности точки t = t;, используя формулу Тейлора-Ито, получим
P(ti J - P(ti) + —aP)h + CДZ, — дZi + O(h2),
тогда, учитывая д z- = S{,
h32
s P(t„)—P., =—C+O (h),
, ,
E ( P, = P(t,) ) = O(h3).
Отсюда можно заключить, что метод Эйлера-Маруямы для решения стохас-
(1) . Используя описанное численное решение стохастического дифференциаль-(1), 1 2009 год получено следующее прогнозное значение потребления на 2010 год.
Р, кВчас
V* / "
V'-./
-Q tr -Q _Q _Q 1— _Q _Q Q_ О ^ Й ПЗ X e; У Q- Q. □. CL m г! Й ™ £ g Q q £T\DvovOio f S & 5 £ = = * | e § § " -G- ra <u о x Jfc 1олученные прогнозныезначения ■■ Реау^ныеданные
Рис. 1. График потребления электроэнергии за 2010 год
Надежность моделей прогнозирования оценивается путем сравнения реальных и прогнозных значений. Эта разница позволяет проверить адекватность модели. Основными оценочными характеристиками качества прогнозной модели являются следующие показатели:
1.Средний процент ошибки MPE (the mean percentage error)
1 n nr — P
MPE =- ------li- -100%,
nY Pr
где Pr - фактическое значение потребления в i -й момент времени;
Pi - расчетное значение потребления в i -й момент времени; n - число наблюдений.
MPE .
2. :
n 0
I (—Pi )2
i=1_____________
n 1 n
I (Pr )2 + IP2
i=1 i=1
Pi r Pi ;
он ближе к нулю, тем ближе сравниваемые ряды.
На основании проведенных расчетов для рассматриваемой модели получены следующие значения оценочных характеристик: MPE = 11,8%, V = 0,09. По
результатам проведенных расчетов можно сделать вывод о применимости полученной модели к прогнозированию потребления.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Лукашов А. В. Риск-менеджмент // Управление корпоративными финансами. - 2005. - № 5. - С. 58-64.
2. . , : . англ. - 6-е изд. - М.: ООО «И.Д. Вильямс», 2008. - 1056 с.
3. Turner Wayne C., Doty S., 1942 - Energy management handbook / Library of Congress Cata-loging-in-Publication Data. - 6th ed., 2007. - 924 p.
Статью рекомендовал к опубликованию к.ф.-м.н. МЛ. Бородицкий.
Кузнецова Инна Юрьевна
Технологический институт федерального государственного автономного образовательного учреждения высшего профессионального образования «Южный федеральный университет» в г. Таганроге.
E-mail: inchikn55@rambler.ru.
347928, г. Таганрог, пер. Некрасовский, 44.
Тел.: +79185772911.
.
Kuznetsova Inna Urevna
Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”.
E-mail: inchikn55@rambler.ru.
44, Nekrasovskiy, Taganrog, 347928, Russia.
Phone: +79185772911.
Postgraduate Student.
