Линейный аналог системы
Система станет линейной, если в аналогичной нелинейной системе пренебречь физическими ограничениями на склады и мощности производства, то есть считать их неограниченными.
Программная реализация линейного аналога системы:
У СКЛ1.Н=СКЛ1.П+Т0.ПН-Т1.ПН
У ПР1.Н=ПР1.П+Т1.ПН-Т2.ПН
У СКЛ2.Н=СКЛ2.П+Т2.ПН-Т3.ПН
У ПР2.Н=ПР2.П+Т3.ПН-Т4.ПН
У СКЛ3.Н=СКЛ3.П+Т4.ПН-Т5.ПН
У ОБЕМ1.Н=ОБЕМ1.П+Т4.ПН
У ОБЕМ2.Н=ОБЕМ2.П+100+20*ВРЕМЯ
Д П1.Н=100+20*ВРЕМЯ
Д ПОСТ.Н=П1+RAND(-15,15)
Т Т0.НБ=ПОСТ.Н
Д Т11.Н=DELAY(П1.Н,1)+М*(DELAY(СКЛ1.Н,1)-DELAY(СКЛ1.Н,2))
Д Т21.Н=CLIP(Т11.Н,СКЛ1.Н,Т11.Н,СКЛ1.Н)
Т Т1.НБ=Т21.Н
Д Т12.Н=CLIP(ПР1.Н,Т11.Н,ПР1.Н,Т11.Н)
Т Т2.НБ=Т12.Н
Д Т13.Н=DELAY(П1.Н,1)+К*(DELAY(ОБЕМ2.Н,2)-DELAY(ОБЕМ1.Н,2))
Д Т23.Н=CLIP(СКЛ2.Н,Т13.Н,СКЛ2.Н,Т13.Н)
Т Т3.НБ=Т23.Н
Д Т14.Н=CLIP(Т13.Н,ПР2.Н,Т13.Н,ПР2.Н)
Т Т4.НБ=Т14.Н
Д Т15.Н=П1.Н+RAND(-15,15)
Т Т5.НБ=CLIP(Т15.Н,СКЛ3.Н,Т15.Н,СКЛ3.Н)
E
И М=0.4
И К=0.01
* ВРЕМЕННЫЕ ПАРАМЕТРЫ
И DT=1
И ВРЕМЯ=0
И ДЛИНА=30
E
Г ПОСТ,СКЛ1,Т11,ПР1,СКЛ2,Т13,ПР2,СКЛ3,Т5,Т15
E
Результат работы программы - см. приложение 2.
Нелинейная модель системы
Программная реализация:
У СКЛ1.Н=СКЛ1.П+Т0.ПН-Т1.ПН
У ПР1.Н=ПР1.П+Т1.ПН-Т2.ПН
У СКЛ2.Н=СКЛ2.П+Т2.ПН-Т3.ПН
У ПР2.Н=ПР2.П+Т3.ПН-Т4.ПН
У СКЛ3.Н=СКЛ3.П+Т4.ПН-Т5.ПН
У ОБЕМ1.Н=ОБЕМ1.П+Т4.ПН
У ОБЕМ2.Н=ОБЕМ2.П+100+20*ВРЕМЯ
Д П1.Н=100+20*ВРЕМЯ
Д ПОСТ.Н=П1+RAND(-15,15)
Д П2.Н=А-СКЛ1.Н
Т Т0.НБ=CLIP(ПОСТ.Н,П2.Н,ПОСТ.Н,П2.Н)
Д Т11.Н=DELAY(П1.Н,1)+М*(DELAY(СКЛ1.Н,1)-DELAY(СКЛ1.Н,2))
Д Т21.Н=CLIP(Т11.Н,СКЛ1.Н,Т11.Н,СКЛ1.Н,1)
Д П3.Н=А-ПР1.Н
Д Т31.Н=CLIP(П3.Н,Т21.Н,П3.Н,Т21.Н)
Т Т1.НБ=Т31.Н
Д Т12.Н=CLIP(ПР1.Н,Т11.Н,ПР1.Н,Т11.Н)
Д П4.Н=А-СКЛ2.Н
Д Т22.Н=CLIP(Т12.Н,П4.Н,Т12.Н,П4.Н)
Т Т2.НБ=Т22.Н
Д Т13.Н=DELAY(П1.Н,1)+К*(DELAY(ОБЕМ2.Н,2)-DELAY(ОБЕМ1.Н,2))
Д Т23.Н=CLIP(СКЛ2.Н,Т13.Н,СКЛ2.Н,Т13.Н)
Д П5.Н=А-ПР2.Н
Д Т33.Н=CLIP(Т23.Н,П5.Н,Т23.Н,П5.Н)
Т Т3.НБ=Т33.Н
Д Т14.Н=CLIP(Т13.Н,ПР2.Н,Т13.Н,ПР2.Н)
Д П6.Н=А-СКЛ3.Н
Д Т24.Н=CLIP(Т14.Н,П6.Н,Т14.Н,П6.Н)
Т Т4.НБ=Т24.Н
Д Т15.Н=П1.Н+RAND(-15,15)
Т Т5.НБ=CLIP(Т15.Н,СКЛ3.Н,Т15.Н,СКЛ3.Н)
E
И А=600
И М=0.4
И К=0.01
* ВРЕМЕННЫЕ ПАРАМЕТРЫ
И DT=1
И ВРЕМЯ=0
И ДЛИНА=30
E
Г ПОСТ,П2,СКЛ1,Т11,П3,ПР1,П4,СКЛ2,Т13,П5,ПР2,П6,СКЛ3,Т5,Т15
E
Результат работы программы - см. приложение 2
Анализ полученных результатов
На основании проведенных экспериментов с моделью были получены некоторые результаты (см. приложение 2), которые помогут нам ответить на ряд вопросов:
во-первых, как влияют запаздывания по выработке и реализации решения на скорость сходимости фактической интенсивности к директивной?
Для ответа на этот вопрос построим графики.
Для первой фазы производства
(с запаздыванием)
(без запаздываний)
Для второй фазы производства.
(с запаздыванием)
(без запаздываний)
Из графиков видно, что в системе без запаздываний фактическая интенсивность сходится к директивной. А в системе, где присутствуют запаздывания фактическая интенсивность отклоняется от директивной.
во-вторых, какова роль ограничений на складские помещения и мощности производства?
Роль значительна, т.к. ограничения не позволяют нам загружать склады так, как мы этого хотим, производить столько, сколько мы запланировали (см приложение 2 - результат работы программы для нелинейной системы) . Гораздо проще дело обстоит с линейным аналогом системы. Здесь нет таких ограничений, а следовательно, не надо беспокоиться об объеме свободного места на складе или о свободных мощностях на фазах производства.