
- •1.1. Методы искусственного интеллекта в прикладных системах и системах принятии решений
- •1.2. Интеллектуальные информационные технологии в прикладных системах и системах принятия решений
- •1.3. Типология задач интеллектуализации систем
- •Лекция 2. Представление знаний в интеллектуальных системах
- •2.1. Модели представления знаний
- •2.2. Системы, основанные на правилах
- •2.3. Системы, основанные на автоматическом доказательстве теорем
- •2.4. Системы, основанные на автоматическом порождении (выдвижении) гипотез
- •Лекция 3. Структура и основные компоненты прикладных интеллектуальных систем
- •3.1. Прикладные системы, основанные на знаниях
- •3.2. Структура системы управления, основанной на знаниях
- •3.3. Структура интеллектуальных систем поддержки принятия решения
- •3.4. Обобщенная структура экспертной системы
- •Лекция 4. Классификация прикладных интеллектуальных систем
- •4.1. Классификация экспертных систем
- •4.2. Примеры прикладных интеллектуальных систем
- •Лекция 5. Основные понятия и определения теории принятия решений
- •5.1. Роли людей в процессе принятия решений
- •5.2. Альтернативы
- •5.3. Критерии
- •5.4. Основные этапы процесса принятия решений
- •5.5. Математические методы теории принятия решений
- •Лекция 6. Принятие решений с помощью статистической проверки гипотез
- •6.1. Статистические решения
- •6.2. Основные задачи статистических решений
- •6.3. Статистическая проверка гипотез
- •6.4. Ошибки решения
- •6.5. Решающее правило при проверке гипотез
- •Лекция 7. Байесовская и последовательная процедуры принятия решения.
- •7.1. Байесовские процедуры принятия решения
- •7.1.1. Байесовская процедура при проверке простой гипотезы
- •7.1.2. Байесовские процедуры в задаче классификации
- •7.2. Принятие решения с помощью последовательной процедуры Вальда
- •Лекция 8. Принятие решения методом дискриминантнного анализа
- •8.1. Классификация в случае, когда распределения классов определены полностью
- •8.1.1. Модель двух нормальных распределений с общей ковариационной матрицей (модель Фишера)
- •8.1.2. Модель двух нормальных распределений с разными ковариационными матрицами
- •8.1.3. Модель нескольких нормальных распределений с общей ковариационной матрицей
- •8.2. Классификация при наличии обучающих выборок
- •8.2.1. Подстановочный алгоритм в модели Фишера
- •8.2.3. Правила классификации
- •8.3. Ошибка решающего правила
- •Лекция 9. Древообразные классификаторы
- •9.1. Назначение древообразных классификаторов
- •9.1. Структура дерева классификации
- •9.3. Вычислительные задачи древообразных классификаторов
- •9.3.1. Определение качества предсказания
- •9.3.2. Выбор разбиений
- •9.3.3. Определение правила прекращения разбиения
- •Лекция 10. Деревья решений
- •9.1. Характеристики дерева решений
- •9.2. Построение дерева решений
- •Лекция 11. Методы прогнозирования
- •11.1. Анализ временных рядов
- •11.1.1. Модель временного ряда
- •11.1.2. Тренд, сезонная и циклическая компоненты
- •11.1.3. Декомпозиция временного ряда
- •11.1.4. Экспоненциальное сглаживание
- •11.2. Каузальные методы прогнозирования
- •11.3. Качественные методы прогнозирования
- •Лекция 12. Основная задача линейного программирования
- •12.1. Математическая модель основной задачи линейного программирования
- •12.2. Задача линейного программирования с ограничениями-неравенствами
- •12.3. Примеры задач линейного программирования
- •12.3.1. Транспортная задача
- •12.3.2. Задача о назначениях
- •Лекция 13. Симплекс-метода решения задачи линейного программирования
- •13.1. Характеристика симплекс–метода
- •13.2. Табличный алгоритм замены базисных переменных
- •13.3. Отыскание опорного решения основной задачи линейного программирования
- •13.4. Отыскание оптимального решения основной задачи линейного программирования
- •Лекция 14. Многокритериальные методы принятия решений при объективных моделях
- •14.1. Объединение критериев
- •14.2. Метод главного критерия
- •14.3. Метод последовательных уступок
- •14.4. Метод целевого программирования
- •14.5. Метод, использующий принцип гарантированного результата
- •14.6. Метод равных наименьших относительных отклонений
- •14.7. Процедура STEM поиска удовлетворительных значений критериев
- •Лекция 15. Выбор Парето–оптимальных решений
- •15.1. Основные определения
- •15.2. Графическая интерпретация
- •15.3. Постановка задачи
- •Лекция 16. Оценка многокритериальных альтернатив с помощью теории полезности
- •16.1. Теория полезности
- •16.2. Принятие решения на основе значения ожидаемой полезности
- •16.3. Многокритериальная теория полезности (MAUT)
- •Лекция 17. Сравнение альтернатив методом аналитической иерархии
- •17.1. Основные этапы метода аналитической иерархии
- •17.2. Декомпозиция задачи
- •17.3. Попарное сравнение критериев и альтернатив
- •17.4. Свойства идеальной матрицы сравнений
- •Лекция 18. Приоритеты для критериев и альтернатив и выбор наилучшей альтернативы в методе анализа иерархий
- •18.1. Вычисление собственных характеристик обратно симметричной матрицы
- •18.2. Вычисление величины приоритетов
- •18.3. Определение наилучшей альтернативы
- •18.4. Проверка согласованности
- •18.5. Пример применения метода анализа иерархий
- •Лекция 19. Оценка многокритериальных альтернатив методами ELECTRE
- •19.1. Этапы подхода, направленного на разработку индексов попарного сравнения альтернатив
- •19.2. Свойства бинарных отношений
- •19.3. Метод ELECTRE I
- •19.4. Метод ELECTRE II
- •19.5. Метод ELECTRE III
- •Лекция 20. Основные понятия и математическая модель игровых методов обоснования решений
- •20.1. Основные понятия теории игр
- •20.2. Математическая модель игры
- •20.3. Нижняя и верхняя цена игры. Принцип минимакса
- •Лекция 21. Методы решения игр
- •21.1. Решение игры в чистых стратегиях
- •21.2. Решение игры в смешанных стратегиях
- •21.3. Упрощение игр
- •21.4. Решение игры 2х2
- •21.5. Графический метод решения (2х2)-игр
- •Лекция 22. Игры 2 х п
- •Лекция 23. Решение игр т х 2 и т х п
- •23.1. Решение игр т х 2
- •23.2. Решение игр т х п
- •Лекция 24. Критерии принятия решений в условиях риска и неопределенности
- •24.1. Основные понятия. Математическая модель
- •24.3. Максиминный критерий Вальда
- •24.4. Критерий минимаксного риска Сэвиджа
- •24.5. Критерий пессимизма-оптимизма Гурвица
- •Литература
Выбор метода прогноза зависит от множества факторов. Отметим некоторые из них: наличие данных, планируемый момент исполнения и желаемая точность прогноза, а также временны'е и стоимостные затраты на его составление.
Многие методы прогноза требуют наличия значительного количества исходных данных. Другие разрабатывают прогноз при отсутствии достоверной количественной информации. Существующие методы прогнозирования можно условно разбить на две группы: 1) количественные и 2) качественные.
Качественные, или экспертные методы строятся на использовании мнений специалистов в соответствующей области (экспертов).
Количественные методы основаны на обработке числовых данных и, в свою очередь, разделяются на каузальные, или причинно-следственные методы и методы анализа временных рядов (ВР). Каузальные методы применяются в тех случаях, когда искомое состояние зависит не только от времени, но и от нескольких переменных. Отыскание математических связей (уравнений и/или неравенств между объясняющими переменными Х и результирующими показателями Y) является целью каузальных методов.
Остановимся на описании особенностей каждого из перечисленных выше методов прогноза.
11.1. Анализ временных рядов
11.1.1. Модель временного ряда
Последовательность случайных наблюдений y1, y2, …, yТ, полученная в фиксированные моменты времени 1, 2,…, Т, называется временным рядом. Его простейшей вероятностной моделью является модель
yt = dt + εt ,
где dt – детерминированная составляющая, εt – случайная составляющая.
Под детерминированной (закономерной) составляющей временного ряда y1, y2, …, yп
мы будем понимать числовую последовательность d1, d2, …, dп, элементы которой dt вычисляются по определенному правилу как функция времени t.
Детерминированная составляющая часто отражает действия каких-либо определенных факторов или причин. Однако для многих ВР в различных областях науки и техники причины, порождающие их закономерные составляющие, могут быть не столь ясными. Тем не менее, их совокупное влияние может быть устойчивым в течение довольно длительных промежутков времени. Это обеспечивает возможность прогноза для подобных временных рядов. Если мы полностью выявим закономерную составляющую в поведении ВР, то оставшаяся часть выглядит хаотично и непредсказуемо. Ее обычно именуют иррегулярной,
или случайной компонентой временного ряда. Обозначим эту компоненту через ε1 ,…, εt ,…,
εn .
Аддитивная и мультипликативная модели. Формы разложения (декомпозиции)
временного ряда на детерминированную и случайную компоненты могут различаться. Укажем наиболее простые из них.
Аддитивной моделью временного ряда называется представление ряда в виде суммы детерминированной и случайной компонент, а именно:
yt = dt + εt при t = 1,...., n или X = D + E.
Мультипликативной моделью ВР называется представление ряда в виде произведения детерминированной и случайной компонент, а именно:
yt = dt × εt при t = 1,...., n или X = D + E.
57
11.1.2. Тренд, сезонная и циклическая компоненты
В детерминированной компоненте ВР dt обычно выделяют три составляющих части: тренд trt , сезонную компоненту st и циклическую компоненту ct . Рассмотрим аддитивную модель ВР. Мы можем записать:
+ st + ct при t = 1,...., n.
Тренд. Анализ временного ряда обычно начинается с выделения именно этой компоненты. Выделение тренда позволяет перейти к дальнейшей идентификации других компонент ВР.
Трендом временного ряда trt при t = 1,...., n называют плавно изменяющуюся, не
циклическую компоненту, описывающую чистое влияние долговременных факторов, эффект которых сказывается постепенно.
Сезонная компонента. Сезонная компонента отражает повторяемость процессов во времени. Она часто присутствует в экономических, метеорологических и других ВР. Сезонная компонента чаще всего служит источником краткосрочных колебаний временного ряда, так что ее выделение заметно снижает вариацию остаточных компонент.
Сезонная компонента st временного ряда при t = 1,...., n описывает поведение,
изменяющееся регулярно в течение заданного периода (года, месяца, недели, дня и т.д.). Она состоит из последовательности повторяющихся циклов.
Главная идея подхода к анализу сезонных компонент заключается в переходе от сравнения всех значений временного ряда между собой к сравнению значений через определенные промежутки времени. Это позволяет снизить оценку вариации ВР около своего среднего значения.
Циклическая компонента занимает промежуточное положение между закономерной и случайной составляющими ВР. Тренд – это плавные изменения, проявляющиеся на больших временных промежутках, а сезонная компонента – это периодическая функция времени, ясно видимая, когда ее период много меньше общего времени наблюдения. Под циклической компонентой обычно подразумевают изменения временного ряда, достаточно плавные и заметные для того, чтобы не включать их в случайную составляющую, но такие, которые нельзя отнести ни к тренду, ни к периодической компоненте.
Циклическая компонента ct временного ряда описывает длительные периоды
относительного подъема и спада. Они состоят из циклов, которые меняются по амплитуде и протяженности.
Изучение циклической компоненты полезно для прогнозирования (особенно краткосрочного).
Модели тренда. Приведем модели трендов, наиболее часто используемые при анализе
временных рядов. Это простая линейная модель |
|
||||
trt = b0 + b1t. |
|
|
|
|
(11.1) |
Эта модель оказывается полезной во многих реальных задачах. Если нелинейный |
|||||
характер тренда очевиден, то может подойти одна из следующих моделей: |
|||||
полиномиальная: tr = b |
+ b t + b |
t2 +…+b tn |
, где значение степени полинома п в |
||
t |
0 |
1 |
2 |
n |
|
практических задачах редко превышает 5;
логарифмическая: trt = exp(b0 + b1t). Эта модель чаще всего применяется для данных,
имеющих тенденцию сохранять постоянные темпы прироста; логистическая: trt = а / (1 + b e−ct );
Гомперца: log(trt ) = a – b·rt , 0 < r < 1,
58
косинусоидальная модель: |
xt = a cos(ωt + θ) + εt . Здесь детерминированная |
компонента представляет собой косинусоидальную функцию с амплитудой a, частотой ω, периодом 2π/ω и фазой θ. Величины a, ω, θ являются константами.
11.1.3. Декомпозиция временного ряда
Декомпозиция временного ряда заключается в выделении во временном ряде его компонент: тренда, сезонной и циклической составляющих. После их исключения временной ряд должен стать стационарным.
Выделение тренда. Для оценки и удаления тренда из ВР чаще всего используется метод наименьших квадратов [1]. Значения временного ряда yt рассматривают как зависимую переменную (отклик), а время t – как независимую переменную (фактор, влияющий на отклик):
yt = f (ti , θ) + εi , i = 1,…, n,
где f – функция тренда (она обычно предполагается гладкой), θ – неизвестные нам параметры (параметры модели временного ряда), εi – независимые и одинаково
распределенные случайные величины, распределение которых мы предполагаем нормальным. Метод наименьших квадратов состоит в том, что мы выбираем функцию тренда так, чтобы
n |
− f (ti ,θ))2 → min |
∑( yti |
|
i=1 |
θ |
Для простого линейного тренда (11.1) с помощью метода наименьших квадратов по
|
|
|
|
ˆ |
ˆ |
: |
формулам находим оценки коэффициентов прямой b0 |
и b1 |
|||||
ˆ |
= |
y , где y – среднее значение yt , |
|
|
||
b0 |
|
|
||||
|
|
n |
|
|
|
|
ˆ |
|
∑( yt |
− y)(t −t(t +1) / 2) |
|
|
|
|
t=1 |
|
|
|
|
|
b1 |
= |
|
|
. |
|
|
n |
|
|
|
|||
|
|
∑(t −t(t +1) / 2)2 |
|
|
t=1
Выделение сезонных эффектов. Для выделения сезонной компоненты есть несколько способов. Их выбор обычно определяется моделью подбираемого временного ряда.
Рассмотрим две распространенных модели ВР. Первая из них включает в себя тренд (trt ),
сезонную ( st ) и случайную ( εt ) компоненты: |
|
yt = trt + st + εt . |
(11.2) |
Вторая модель, кроме перечисленных выше компонент, включает еще и циклическую |
|
компоненту ( ct ): |
|
yt = trt + st + ct + εt . |
(11.3) |
Сезонные эффекты на фоне тренда. Предположим, что рассматриваемый временной ряд y1, y2, …, yп может быть описан аддитивной моделью (11.2). Пусть р – период
последовательности st , так что st = st+ p для всякого t. Необходимо оценить значения st по наблюдениям yt при том, что величина р известна.
59

|
22 |
|
|
|
|
|
|
|
|
|
22 |
|
20 |
|
|
|
|
|
|
|
|
|
20 |
|
18 |
|
|
|
|
|
|
|
|
|
18 |
|
16 |
|
|
|
|
|
|
|
|
|
16 |
Урожайность |
14 |
|
|
|
|
|
|
|
|
|
14 |
12 |
|
|
|
|
|
|
|
|
|
12 |
|
10 |
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
8 |
|
6 |
|
|
|
|
|
|
|
|
|
6 |
|
4 |
|
|
|
|
|
|
|
|
|
4 |
|
2 |
5 |
10 |
15 |
20 |
25 |
30 |
35 |
40 |
45 |
2 |
|
0 |
50 |
|||||||||
|
|
|
|
|
|
Время |
|
|
|
|
|
|
|
|
|
|
Урожайность |
Тренд |
|
|
|
|
Рис.11.1. Исходный временной ряд и линейный тренд
Для этого сначала мы должны оценить тренд trt . Это можно сделать с помощью
метода наименьших квадратов. Обозначим через trˆt полученную оценку тренда. |
Затем для |
каждого сезона i, 1 ≤ i ≤ p, рассмотрим все относящиеся к нему разности |
|
yi – trˆi , yi+ p –trˆi+ p , …, yi+mp –trˆi+mp . |
(11.4) |
Каждое из этих отклонений yi от trˆi можно рассматривать как результат влияния сезонных изменений. Усреднение этих разностей дает нам оценку сезонной компоненты si . В качестве простейшей оценки можно взять среднее арифметическое, т.е. положить
m |
|
sˆi =1 / (т + 1) ∑( yi+lp −trˆi+lp ) для i = 1,…, p. |
(11.5) |
l=0
Вкачестве других оценок sˆi можно взять взвешенное среднее, медиану и т.д.
Перечисленные характеристики уменьшают влияние резко выделяющихся наблюдений. Метод скользящих средних. При наличии в ряде циклической компоненты расчет
сезонных эффектов несколько отличается от описанного выше. В этом случае для выяснения сезонных вкладов необходимо оценить не только тренд, но и циклическую компоненту. Проще всего одновременно оценить тренд и циклическую компоненту можно с помощью скользящего среднего (СС).
Метод скользящего среднего основан на переходе от начальных значений ряда к их средним значениям на интервале времени, длина которого выбрана заранее. При этом сам выбранный интервал скользит вдоль ряда. В методе СС можно использовать различные виды усреднения значений ряда: среднее арифметическое (простое или с весами), медианы и др. К медианному сглаживанию прибегают тогда, когда среди наблюдений есть выбросы.
60