Имитационное моделирование
Принципы и подходы к формированию входных данных для имитационных моделей. Построение имитационной модели производственного участка и оценка стабильности его работы, рекомендации по совершенствованию. Принципы прогнозирования объема выпуска изделий.
Рубрика | Экономико-математическое моделирование |
Вид | контрольная работа |
Язык | русский |
Дата добавления | 09.11.2016 |
Размер файла | 157,6 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Размещено на http://www.allbest.ru/
Задача №1
Формирование входных данных для имитационных моделей одна из важнейших задач, так как упрощение, пренебрежение или сведение входных данных к каким-либо аналитическим зависимостям делает модель, как правило, неадекватной моделируемому объекту. К входным данным по терминологии систем обычно относят входные сигналы, управляющие сигналы, параметры системы и выходные сигналы от одних блоков, поступающие на вход каких-либо других блоков данной модели.
В общем случае входные сигналы из внешней среды можно представить в виде динамических рядов, фиксирующих значение какого-то показателя в определенные моменты времени, или какого-либо потока событий, появляющихся в заранее неизвестные моменты времени. Характеристика или значение того или иного события могут быть как одинаковыми, так и разными.
Разберем четыре возможных типа входных сигналов: Y1 поток случайных событий во времени, где событие представлено лишь фактом его появления в какие-то неизвестные заранее моменты времени; Y2 поток случайных событий во времени, у которого событие характеризуется не только фактом его появления, но и конкретным неодинаковым во времени числовым значением данного события; Y3 дискретный ряд, характеризующий значение показателя в определенные регулярные моменты времени t, t,…, t; Y4 значение показателя, которое может быть получено в любой момент времени (непрерывная функция).
Все эти типы сигналов, различные по форме, как регулярные, так и не регулярные во времени можно записать в форме единого динамического ряда с событиями, происходящими в регулярные моменты времени t', t',…, t',…, t'. Для входных сигналов вида Y3 и Y4 это очевидно, так как или они сами являются регулярными или их значение может быть измерено в регулярные промежутки времени, а для сигналов Y1 и Y2 всегда можно на оси времени подобрать такие регулярные интервалы [t, t], при которых будет меньше допустимой наперед заданной ошибки появления случайного события, если оно произошло не в регулярные моменты времени t или t. Следовательно, данные потоки событий можно всегда сделать регулярными, подобрав соответствующие интервалы.
Необходимо построить имитационную модель производственного участка. Оценить стабильность работы участка. Выработать рекомендации для улучшения стабильности работы участка. Входными данными является число деталей, прошедших обработку.
Входные данные
Месяц |
Годы |
||||||
2011 |
2012 |
2013 |
2014 |
2015 |
2016 |
||
Январь |
29 |
29 |
27 |
29 |
38 |
2 |
|
Февраль |
12 |
35 |
43 |
47 |
1 |
15 |
|
Март |
5 |
3 |
14 |
13 |
18 |
20 |
|
Апрель |
46 |
27 |
49 |
10 |
8 |
||
Май |
16 |
39 |
1 |
36 |
42 |
||
Июнь |
12 |
36 |
9 |
40 |
30 |
||
Июль |
39 |
30 |
20 |
9 |
38 |
||
Август |
10 |
9 |
40 |
0 |
43 |
||
Сентябрь |
18 |
14 |
43 |
9 |
34 |
||
Октябрь |
27 |
44 |
21 |
9 |
1 |
||
Ноябрь |
30 |
34 |
31 |
3 |
6 |
||
Декабрь |
18 |
6 |
10 |
30 |
3 |
Рассмотрим график изменения значения количества деталей, на нем отчетливо просматриваются значительные циклические колебания.
График изменения значения количества деталей
Они отражают, прежде всего, циклический характер загрузки участка. Существует способ фильтрации таких циклических явлений. В таблице 2 приведены предварительные расчеты сглаживания модельного ряда данных по трем точкам и найдены относительные колебания.
Расчеты сглаживания модельного ряда
Месяц |
Год |
y |
y1 (сглаж) |
Относ. колеб |
Колеб. по годам |
Невыр. волна |
Выров. Волна |
Отклон. от 100% |
Отклон. в днях |
Реком. |
|
Январь |
2011 |
29 |
509,06 |
101,81 |
104,93 |
4 |
|||||
Февраль |
12 |
15,33 |
78,26 |
673,86 |
112,31 |
115,75 |
15 |
3 |
Добавить |
||
Март |
5 |
21,00 |
23,81 |
332,99 |
66,60 |
68,64 |
-31 |
-6 |
Отдыхать |
||
Апрель |
46 |
22,33 |
205,97 |
639,19 |
127,84 |
131,75 |
31 |
6 |
Добавить |
||
Май |
16 |
24,67 |
64,86 |
467,74 |
93,55 |
96,41 |
-3 |
||||
Июнь |
12 |
22,33 |
53,73 |
469,58 |
93,92 |
96,79 |
-3 |
||||
Июль |
39 |
20,33 |
191,80 |
556,56 |
111,31 |
114,72 |
14 |
2 |
Добавить |
||
Август |
10 |
22,33 |
44,78 |
324,40 |
64,88 |
66,87 |
-33 |
-6 |
Отдыхать |
||
Сентябрь |
18 |
18,33 |
98,18 |
565,68 |
113,14 |
116,60 |
16 |
3 |
Добавить |
||
Октябрь |
27 |
25,00 |
108,00 |
453,68 |
90,74 |
93,52 |
-6 |
-1 |
Отдыхать |
||
Ноябрь |
30 |
25,00 |
120,00 |
592,86 |
118,57 |
122,20 |
22 |
4 |
Добавить |
||
Декабрь |
18 |
25,67 |
70,13 |
348,43 |
69,69 |
71,82 |
-28 |
-5 |
Отдыхать |
||
Январь |
2012 |
29 |
27,33 |
106,10 |
среднее=> |
97,03 |
100,00 |
||||
Февраль |
35 |
22,33 |
156,72 |
||||||||
Март |
3 |
21,67 |
13,85 |
||||||||
Апрель |
27 |
23,00 |
117,39 |
||||||||
Май |
39 |
34,00 |
114,71 |
||||||||
Июнь |
36 |
35,00 |
102,86 |
||||||||
Июль |
30 |
25,00 |
120,00 |
||||||||
Август |
9 |
17,67 |
50,94 |
||||||||
Сентябрь |
14 |
22,33 |
62,69 |
||||||||
Октябрь |
44 |
30,67 |
143,48 |
||||||||
Ноябрь |
34 |
28,00 |
121,43 |
||||||||
Декабрь |
6 |
22,33 |
26,87 |
||||||||
Январь |
2013 |
27 |
25,33 |
106,58 |
|||||||
Февраль |
43 |
28,00 |
153,57 |
||||||||
Март |
14 |
35,33 |
39,62 |
||||||||
Апрель |
49 |
21,33 |
229,69 |
||||||||
Май |
1 |
19,67 |
5,08 |
||||||||
Июнь |
9 |
10,00 |
90,00 |
||||||||
Июль |
20 |
23,00 |
86,96 |
||||||||
Август |
40 |
34,33 |
116,50 |
||||||||
Сентябрь |
43 |
34,67 |
124,04 |
||||||||
Октябрь |
21 |
31,67 |
66,32 |
||||||||
Ноябрь |
31 |
20,67 |
150,00 |
||||||||
Декабрь |
10 |
23,33 |
42,86 |
||||||||
Январь |
2014 |
29 |
28,67 |
101,16 |
|||||||
Февраль |
47 |
29,67 |
158,43 |
||||||||
Март |
13 |
23,33 |
55,71 |
||||||||
Апрель |
10 |
19,67 |
50,85 |
||||||||
Май |
36 |
28,67 |
125,58 |
||||||||
Июнь |
40 |
28,33 |
141,18 |
||||||||
Июль |
9 |
16,33 |
55,10 |
||||||||
Август |
0 |
6,00 |
0,00 |
||||||||
Сентябрь |
9 |
6,00 |
150,00 |
||||||||
Октябрь |
9 |
7,00 |
128,57 |
||||||||
Ноябрь |
3 |
14,00 |
21,43 |
||||||||
Декабрь |
30 |
23,67 |
126,76 |
||||||||
Январь |
2015 |
38 |
23,00 |
165,22 |
|||||||
Февраль |
1 |
19,00 |
5,26 |
||||||||
Март |
18 |
9,00 |
200,00 |
||||||||
Апрель |
8 |
22,67 |
35,29 |
||||||||
Май |
42 |
26,67 |
157,50 |
||||||||
Июнь |
30 |
36,67 |
81,82 |
||||||||
Июль |
38 |
37,00 |
102,70 |
||||||||
Август |
43 |
38,33 |
112,17 |
||||||||
Сентябрь |
34 |
26,00 |
130,77 |
||||||||
Октябрь |
1 |
13,67 |
7,32 |
||||||||
Ноябрь |
6 |
3,33 |
180,00 |
||||||||
Декабрь |
3 |
3,67 |
81,82 |
||||||||
Январь |
2016 |
2 |
6,67 |
30,00 |
|||||||
Февраль |
15 |
12,33 |
121,62 |
||||||||
Март |
20 |
Далее рассчитываем циклическую волну работы производственного участка по средней арифметической. Относительные колебания группируются по месяцам (по модельному времени t), и их сумма занесена в графу таблицы 2 «Колеб. по годам».
Затем, были вычислены не выровненные колебания нахождением средней арифметической и занесены в столбец «Невыр. волна».
Значения величин в столбце «Выровненная волна» находились через пропорцию.
В итоге построен график выровненной волны.
Выровненная волна
Формирование управляющей компоненты для имитационной модели это подбор факторов управления, влияющих на траекторию существования объекта моделирования.
Из графика изменения выровненной волны на рисунке 2, заметно, что колебания уровня загрузки участка составляет около 60%. В этом случае необходимо разработать мероприятия, связанные с выравниванием уровня загрузки. Для этого необходимо в имитационную модель ввести факторы, управляющие загруженностью участка, например, изменение коэффициента сменности.
Фактор изменения коэффициента сменности предполагает, что в периоды перегрузок и спадов следует соответственно увеличивать или уменьшать коэффициент сменности: в период перегрузок запретить отпуск рабочих и служащих или нанять дополнительное число работников, а в периоды спадов, наоборот, предоставлять максимально возможное число отпусков, осуществить профилактику оборудования и т.п.
В таблице 2 приведено расчетное число рабочих дней (смен), на которое нужно изменить величину рабочего времени помесячно, что бы сгладить пики уровня загрузки. Выполнено это следующим образом.
Число (фактическое) рабочих дней взято без суббот, воскресений и официальных праздников. Полагаем, что расчетное увеличение нагрузки не может превысить 4-х дней в месяц (4 субботы), расчетное уменьшение не может превышать 2 дней. Для упрощения расчетов можно принять среднее число рабочих дней в месяце - 20 дней, но для повышения точности имитационной модели следует брать фактическое число рабочих дней в каждом месяце.
Тогда 20 дней соответствует 100% загрузке, а 1 день соответствует х% нагрузке, т.е. х=5%, что означает, что увеличение числа рабочих дней в месяце на 1 день изменит уровень нагрузки предприятия на 5%. По выровненной волне рассчитали отклонение уровня загрузки от 100%.
Из таблицы 2 видно, что в феврале, апреле, июле, сентябре и ноябре, выровнять загрузку производственного участка можно одним фактором - увеличением коэффициента сменности. А в марте и августе спад настолько велик, что недостает 6 дней для выравнивания до 100%. Эти 6 дней соответствуют 33% от уровня загрузки.
Эти расчеты проделаны без учета структуры заказов производственного участка, например, того что изготавливают различные по своей сложности и конфигурации детали.
Задача №2
При прогнозировании различных показателей экономических или производственных систем часто оказывается, что тенденция изменения прогнозируемого показателя меняется. Такое явление может происходить из-за изменения интенсивности факторов, действующих на показатели, или в связи с влиянием новых факторов, которые не оказывали действие на начальные значения исследуемого динамического ряда. Например, моменты времени, в которые были измерены значения некоторого показателя. Допустим, что для значений времени на наш показатель действовала определенная группа факторов , а для моментов времени другая группа факторов . Если для группы факторов можно кривой роста вида
, (1)
то для второй части ряда, начиная с момента времени и до момента времени можно построить модель кривой роста вида
. (2)
Если бы группы факторов и оказывали одинаковое действие на показатель то приросты и были бы одинаковы и можно было бы построить единую модель кривой роста:
. (3)
Если указанные приросты не близки друг к другу, т.е. приросты для начальных членов ряда сильно отличаются от приростов для конечных его членов, то целесообразно учесть изменение тенденций при помощи специальной адаптивной процедуры. Имитационная модель, имеющая механизм изменения параметров в зависимости от оценки своих текущих значений, носит название адаптивной.
Механизмы изменения параметров могут быть различными. Выбор того или иного из них и его настройка зависят от целей исследования и вида объекта моделирования. Наиболее часто для адаптивных моделей краткосрочного прогнозирования используются две схемы: скользящего среднего (СС - модели) и авторегрессии (АР - модели).
Схема СС - моделей обеспечивается путем вычисления взвешенных средних всех предыдущих значений, причем весовые коэффициенты, как правило, убывают по мере удаления от последнего значения ряда, т.е. последним наблюдениям придается большее значение или больший вес.
Такой подход соответствует гипотезе о том, что последние значения ряда более информативны для прогнозируемых значений, чем начальные. Вместе с тем, это всего лишь гипотеза, которая тем более вероятна, чем меньше случайная или циклическая составляющие у данных членов ряда.
Для АР - моделей оценка текущего значения ряда определяется взвешиванием не всех, а некоторых предшествующих уровней, и весовые коэффициенты не ранжируются. В данном случае весомость значений ряда или их информационная ценность, или близость к прогнозируемым значениям определяется исходя из тесноты связи между ними.
Рассмотрим пример СС - моделей.
В практике имитационного прогнозирования наиболее часто применяются СС - модели Брауна и Хольта. Данные модели прогнозируют, предполагая, что тенденция исследуемого процесса моделируется линейной зависимостью с переменными коэффициентами.
Процесс прогнозирования для такого типа моделей описывается формой
, (4)
имитационный производственный прогнозирование
где оценка текущего (t-го) значения ряда; оценка текущего прироста; прогнозная оценка значения ряда в момент времени t на k шагов вперед.
Требуется построить имитационную адаптивную модель для прогнозирования объема выпуска изделий на предприятии на 2 месяца вперед, имея данные выпуска за 9 мес.
Данные выпуска за 9 месяцев
Y(t) |
31 |
30 |
33 |
36 |
35 |
36 |
37 |
38 |
39 |
Этап 1. Примем значение l=5. Применим метод наименьших квадратов (МНК) для вычисления и .
Этап 2. k=1; как эксперты, принимаем значения =0,6; Вычисляем Yp(t).
Этап 3. Находим ошибку E(t).
Этап 4. Откорректируем коэффициенты с учетом полученной ошибки.
Этап 5. Осуществляем прогноз для второго шага.
Этап 6. Повторяем этапы 3…5 до тех пор, пока t<n, т.е. t<9. Процедура заканчивается при t=9.
Этап 7. Исследуем качество модели при =0,6 с помощью средней относительной ошибки аппроксимации.
,
где n - количество наблюдений / экспериментов.
Выполненные расчеты
Y(t) |
А0 (t) |
А1 (t) |
Yр(t) |
E(t) |
в |
E(относ) |
|
31 |
28,80 |
1,40 |
30,20 |
0,80 |
0,6 |
2,58 |
|
30 |
30,71 |
1,53 |
32,24 |
-2,24 |
7,47 |
||
33 |
30,81 |
1,17 |
31,98 |
1,02 |
3,10 |
||
36 |
32,63 |
1,33 |
33,96 |
2,04 |
5,65 |
||
35 |
35,27 |
1,66 |
36,93 |
-1,93 |
5,50 |
||
36 |
35,69 |
1,35 |
37,04 |
-1,04 |
2,90 |
||
37 |
36,38 |
1,18 |
37,56 |
-0,56 |
1,51 |
||
38 |
37,20 |
1,09 |
38,30 |
-0,30 |
0,78 |
||
39 |
38,11 |
1,05 |
39,15 |
-0,15 |
0,39 |
||
39,06 |
1,02 |
40,08 |
3,32 |
Расчетное значение Еотн. при изменении значения бетта
в |
E(относ) |
|
0 |
4,71 |
|
0,1 |
4,48 |
|
0,2 |
4,14 |
|
0,3 |
3,76 |
|
0,4 |
3,54 |
|
0,5 |
3,37 |
|
0,6 |
3,32 |
|
0,7 |
3,32 |
|
0,8 |
3,37 |
|
0,9 |
3,51 |
|
1 |
3,75 |
Промежуточные расчеты
t |
Y(t) |
A+B*t |
|
1 |
31 |
30,2 |
|
2 |
30 |
31,6 |
|
3 |
33 |
33 |
|
4 |
36 |
34,4 |
|
5 |
35 |
35,8 |
|
A |
28,8 |
||
B |
1,4 |
График отклонения
Вывод: Из таблицы видно, что наиболее точный прогноз даст имитационная адаптивная модель при в=0,6, так как она имеет наименьшую погрешность Еотн=3,32%. При этом значении в будем осуществлять прогноз.
В следующем месте система предприятие сделает предположительно 39 деталей, в дальнейшем 40.
Список используемой литературы
1. Кобелев Н.Б. Основы имитационного моделирования сложных систем: Учебное пособие. - М.: Дело, 2003. - 336 с.
2. Орлова И.В., Половников В.А. Экономико-математические методы и модели: компьютерное моделирование. Учебное пособие. Гриф УМО вузов России. - М: Издательство: Вузовский учебник, 2014. - 389 с.
3. Бережная Е.В., Бережной В.И. Математические методы моделирования экономических систем: Учеб. пособие. - 2-е изд., перераб. и доп. - М.: Финансы и статистика, 2006. - 432 с: ил.
4. Лычкина Н.Н. Имитационное моделирование экономических процессов: Учеб. пособие. - М.: Академия АйТи, 2005.
Размещено на Allbest.ru
Подобные документы
Метод имитационного моделирования в разработке экономико-математических моделей для учета неопределенности статистики предприятий. Функционирование имитационной модели изготовления малогабаритного стула: время работы и коэффициенты загрузки оборудования.
курсовая работа [2,0 M], добавлен 16.11.2010Процедура проведения имитационных экспериментов с моделью исследуемой системы. Этапы имитационного моделирования. Построение концептуальной модели объекта. Верификация и адаптация имитационной модели. Метод Монте-Карло. Моделирование работы отдела банка.
курсовая работа [549,5 K], добавлен 25.09.2011Сущность понятия термина "имитация". Сущность этапов имитационного эксперимента. Основные принципы и методы построения имитационных моделей. Типы систем массового обслуживания. Логико-математическое описание, выбор средств и анализ работы модели.
реферат [7,5 M], добавлен 25.11.2008Методы построения имитационных моделей экономических объектов. Проведение анализа по результатам численных экспериментов на имитационной модели оптового магазина. Выявление закономерностей, которые помогут в проведении кадровой политики предприятия.
курсовая работа [389,0 K], добавлен 28.11.2010Построение имитационной модели бизнес-процесса "Управление инцидентами" компании "МегаФон" с целью прогнозирования совокупной стоимость ИТ-сервиса по обслуживанию инцидентов. Разработка моделирующих алгоритмов для реализации компьютерных программ модели.
курсовая работа [2,6 M], добавлен 09.04.2012Построение имитационной модели технологического процесса методом Монте-Карло, ее исследование на адекватность. Оценка и прогнозирование выходных характеристик технологического процесса с помощью регрессионных моделей. Разработка карт контроля качества.
курсовая работа [1,2 M], добавлен 28.12.2012Изучение и отработка навыков математического моделирования стохастических процессов; исследование реальных моделей и систем с помощью двух типов моделей: аналитических и имитационных. Основные методы анализа: дисперсионный, корреляционный, регрессионный.
курсовая работа [701,2 K], добавлен 19.01.2016Разработка программной имитационной модели работы билетной кассы железнодорожного вокзала на языке GPSS World. Описание пошаговой работы программы и плоскости отклика модели. Исследование функционирования модели на чувствительность изменения факторов.
курсовая работа [1,3 M], добавлен 22.06.2015Теоретические основы математического прогнозирования продвижения инвестиционных инструментов. Понятие системы имитационного моделирования. Этапы построения моделей экономических процессов. Характеристика ООО "Брянск-Капитал". Оценка адекватности модели.
курсовая работа [2,1 M], добавлен 20.11.2013Базовые принципы и приемы, используемые при имитационном моделировании доходности финансового актива. Построение модели, способной прогнозировать доходность акции компании "РосНефть" через индекс MICEX и нефть марки Brent. Проверка модели на адекватность.
контрольная работа [415,5 K], добавлен 11.12.2014