- •Лукьяненко Михаил Васильевич чурляева Наталья Петровна моделирование технических систем и процессов
- •Оглавление
- •6. Использование теории Марковских процессов и временных рядов
- •7. Использование теории очередей при моделировании работы атс
- •8. Использование метода сетевого планирование при моделировании
- •Предисловие
- •1. Основные этапы моделирования систем
- •1.1. Построение концептуальной модели системы и её формализация
- •1.2. Алгоритмизация модели и ее компьютерная реализация
- •1.3. Получение и интерпретация результатов моделирования
- •Контрольные вопросы
- •2. Моделирование систем массового обслуживания.
- •2.1. Системный анализ смо
- •2.2. Статистический анализ смо.
- •2.3. Операционный анализ смо.
- •Контрольные вопросы
- •3. Имитационное моделирование.
- •3.1. Моделирование работы сборочного цеха с программированием на языке высокого уровня.
- •3.2. Моделирование работы ремонтного цеха с использованием языка имитационного моделирования систем.
- •Контрольные вопросы
- •4. Моделирование процессов во времени.
- •4.1 Моделирование эволюции систем на основе теории Марковских процессов
- •4.2. Анализ процессов с помощью временных рядов
- •4.3. Оценка точности регрессионных моделей.
- •Контрольные вопросы
- •Глава 5. Моделирование сетевых структур.
- •5.1. Сетевое моделирование
- •5.2. Сетевое планирование.
- •5.3. Динамическое программирование при моделировании в сетях.
- •Контрольные вопросы
- •6. Использование теории Марковских процессов и временных рядов при моделировании работы блоков шб3Бт и шбт4Бт.
- •6.1. Паспортные данные, схемы исследуемых блоков и анализ возможных неисправностей.
- •6.2 Анализ и прогноз для блока шб3Бт
- •Выводы по блоку шбт3Бт.
- •6.3. Анализ и прогноз работоспособности для блока шб4Бт
- •6.4. Подведение итогов моделирования и выдача рекомендаций. Общие выводы по обоим блокам.
- •Подведение итогов и выдача практических рекомендаций
- •7. Использование теории очередей при моделировании работы атс Нicоm 353
- •7.1 Описание объекта моделирования.
- •7.2 Цель моделирования.
- •7.3. Концептуальная модель системы и методы исследования.
- •7.4. Получение результатов моделирования для группы №1.
- •7.5. Получение результатов моделирования для группы № 2.
- •7.6. Получение результатов моделирования для группы № 5.
- •7.7. Интерпретация результатов моделирования, практические выводы и рекомендации.
- •8. Использование метода сетевого планирования при моделировании регламентных работ перед техобслуживанием.
- •8.1 Введение.
- •8.2 Основные регламентные работы перед проведением техобслуживания.
- •8.3 Краткое описание последовательности основных регламентных работ
- •8.4. Построение сетевого графика без учёта аккумуляторных и карбюраторных работ.
- •8.5. Расчёт сетевого графика.
- •8.6. Сетевой график с включением аккумуляторных и карбюраторных работ.
- •8.7. Анализ полученных результатов с учётом мнения руководителя автохозяйства.
- •8.8. Общий вывод по проведенному исследованию
- •Варианты заданий для моделирования условных объектов.
- •Заключение
4.2. Анализ процессов с помощью временных рядов
Временной ряд представляет собой ряд наблюдений за каким-то определенным параметром изучаемой системы в дискретные, равноотстоящие моменты времени. Особый интерес при этом представляет проблема прогнозирования будущих состояний системы исходя из текущих и прошлых значений временного ряда. В зависимости от исходных предположений относительно прогнозирующей функции и вероятностных интервалов прогнозирования результаты прогнозирования будущих состояний системы могут сильно различаться.
Методы прогнозирования можно разделить на статистические и причинно-следственные методы. Причинно-следственные метода весьма сложны и включают помимо модели самой системы построение модели окружающей её среды. Прогноз, сделанный на основе такой модели, позволяет объяснять будущее системы, исходя из её внутренних свойств. Статистические методы прогноза исходят исключительно из внешних проявлений поведения системы.
Обычно статистические методы включают в рассмотрение: 1) исходные данные; 2) модели прогноза; 3) методы сглаживания; 4) ошибки прогноза.
1). При рассмотрении исходных данных следует исключить выбросы, т.е. те наблюдения, которые не характерны для прогнозируемого процесса. После этого удобно принять предположение о том, что ошибки в исходных данных подчиняются закону нормального распределения.
2). Если в момент Т задана последовательность наблюдении Хt для t < Т, то прогнозирующая модель задает множество выходных переменных в векторной форме XT+t = АTF(t), где вектор АT представляет собой вектор коэффициентов модели, получаемых по результатам наблюдений за поведением системы до наступления момента Т; матрица F – некоторый набор аппроксимирующих функций. Многие модели могут быть представлены полиномами Fi:
F1(t)=1, F2(t)=1/2t(t-1) и т.д.
Довольно часто для анализа локальных изменений достаточно использовать полином 1-й степени (прямую).
3). Система сглаживающих множителей включает весовые множители Wj при каждом j-м наблюдении. Простейший ряд равных весовых множителей Wj = 1 для всех 0 < j < Т придает одинаковую значимость каждому члену ряда. Скользящее среднее представляет собой оценку по методу наименьших квадратов единственной константы модели прогноза для представления исходных данных с одинаковыми весовыми множителями. Этому случаю соответствует аппроксимирующая функция F(t) = 1 для всех t.
4). Ошибки прогноза оцениваются аналогично погрешностям физических изменений.
Как форму прогноза, статистические методы позволяют выявить тенденции временных рядов. В качестве примера рассмотрим выявление сложившихся тенденций изменения отдельного параметра производственной системы. Тенденция может быть явной либо неявной. Явная тенденция не требует дополнительных методов исследования в отличие от неявной тенденции. Для выявления неявной основной тенденции используются такие приемы, как укрупнение интервалов наблюдений, расчет средних уровней, построение средней, аналогическое выравнивание рядов динамики.
Рассмотрим, как выявляется тенденции на примере изменения некоторого производственного показателя (выпуска продукции) с помощью метода скользящей средней. Сущность метода скользящей средней заключается в том, что вычисляются средние значения показателя сначала, например, из первых четырех значений ряда динамики, т.е. с первого по четвертый месяц, затем также из четырех значений, но начиная со второго по пятый месяц, далее берется среднее с третьего по шестой месяцы и т. д.
месяц |
Выпуск (т) продукции |
Ср. за месяцы |
Выпуск продукции |
Центрированная скользящая средняя |
Центральное звено | |
за 4 месяца |
в среднем за месяц | |||||
1 |
1,4 |
|
|
|
|
|
2 |
2,2 |
1-4 |
10,4 |
2,6 |
|
|
3 |
4,6 |
2-5 |
10,9 |
2,72 |
(2,6+2,72)/2=2,66 |
3 |
4 |
2,2 |
3-6 |
11,0 |
2,75 |
2,73 |
4 |
5 |
2,0 |
4-7 |
14,1 |
3,53 |
3,14 |
5 |
6 |
2,2 |
5-8 |
14,6 |
3,65 |
3,6 |
6 |
7 |
7,7 |
6-9 |
21,2 |
5,3 |
4,48 |
7 |
8 |
2,7 |
7-10 |
26,7 |
6,7 |
(5,3+6,7)/2=6,0 |
8 |
9 |
2,6 |
и т.д. |
|
|
|
|
Каждое звено скользящей средней характеризует средний уровень показателя за соответствующий период. При графическом изображении каждое звено относят к тому периоду, где проводилось усреднение.
Другим распространенным методом выявления тенденций является выравнивание данных по прямой. Прямая xt=a0+a1t, аппроксимирующая исходные данные, определяется методом наименьших квадратов.
Для иллюстрации рассмотрим следующий пример аппроксимации:
Месяц |
Выпуск продукции Xt (н.час) |
Отсчет времени, t |
Xt*t |
t2 |
Выровненные значения | ||
1 |
36.9 |
-9 |
-242.1 |
81 |
23.83 | ||
2 |
25.1 |
-7 |
-175.7 |
49 |
25.17 | ||
3 |
30.9 |
-5 |
-154.5 |
25 |
26.51 | ||
4 |
29.3 |
-3 |
-77.9 |
9 |
27.85 | ||
5 |
22.7 |
-1 |
-22.7 |
|
29.79 | ||
6 |
26.9 |
1 |
26.9 |
1 |
30.53 | ||
7 |
32.3 |
3 |
96.9 |
9 |
31.87 | ||
8 |
24.0 |
5 |
220.0 |
25 |
33.21 | ||
9 |
27.1 |
7 |
259.7 |
49 |
34.55 | ||
10 |
33.4 |
9 |
300.6 |
81 |
35.85 | ||
Итого |
298.6 |
|
221.2 |
330.0 |
298.6 |
Методом наименьших квадратов определяется свободный член аппроксимирующей линейной зависимости а0 и её наклон а1.
;
В результате имеем следующий временной ряд, выровненный по прямой:
В некоторых случаях целесообразно использовать временной ряд с экспоненциальным затуханием весовых множителей wj=aj, а<1. Эти множители придают большой вес более поздней информации. Величина а определяется для различных уровней и тенденций путем систематического исследования точности прогнозов, получаемых с разными весовыми функциями. Уточнение прогноза производится по принципу обратной связи, т.е. новые прогнозы корректируют с учетом ошибок в предыдущих. Существуют теоретические обоснования величины a=0.9.
Возможен подход к прогнозированию, ставящий каждому наблюдению до его проведения в соответствие ряд первичных вероятностных коэффициентов модели. После того, как получены результаты наблюдений, по теореме Байеса определяют апостериорные вероятности. Теорема Байеса имеет вид:
где р(Е) – априорная вероятность; S - информационная выработка (данные); р'(Е) - апостериорная вероятность события Е; P(S|E) - вероятность S при данном событии Е.
Основываясь на апостериорных вероятностях можно вычислить распределение вероятностей прогнозируемой величины. Такой подход позволяет игнорировать кратковременные изменения прогнозируемой величины и в то же время выявлять изменения уровня и тенденции временного ряда.