Главная Контакты Добавить в избранное Авторы Вопросы и ответы
,

УДК 681.5

ОЦЕНКА АДЕКВАТНОСТИ АНАЛИТИЧЕСКОЙ МОДЕЛИ УПРАВЛЕНИЯ ЗАПАСАМИ МАТЕРИАЛЬНЫХ РЕСУРСОВ В МНОГОУРОВНЕВОЙ ЛОГИСТИЧЕСКОЙ СИСТЕМЕ

Становский А.Л., Кухаренко С.В., Колчин Р.В.

Введение

Одной из важнейших задач, возникающих в процессе проектирования АСУ запасами материальных ресурсов (МР) является разработка ее математического обеспечения. Данные, обеспечивающие подсистему АСУ, реализуется, как правило, в виде комплекса моделей управления запасами МР. Пример такой математической (аналитической) модели управления запасами МР в двухуровневой логистической системе (ЛС), был рассмотрен в [9]. В тоже время, необходимым условием применимости аналитической модели управления запасами МР является требование адекватности получаемых с ее помощью теоретических результатов данных практического опыта.

 

Анализ предшествующих публикаций

Известно, что адекватность предполагает воспроизведение моделью, с необходимой полнотой, всех свойств объекта, существенных для целей данного исследования. Так как всякая модель представляет собой упрощенное описание объекта моделирования, никогда нельзя говорить об абсолютной адекватности, при которой модель по всем параметрам соответствует оригиналу. Следовательно, оценка степени сходства может опираться только на оценку отличия модели от оригинала. Оценивание отличия наталкивается, естественным образом, на большие трудности, так как обычно невозможно использовать для сравнения объект во всей его действительной целостности [1].

Однако, зная структуру ЛС и ее параметры, некоторые характеристики внешней, по отношению к системе, среды, а также, используя опыт экспертов, можно ввести вполне определенные предположения относительно состояния системы в различных условиях функционирования. Таким образом, используя и анализируя существующую информационную базу и синтезируя, путем собственных рассуждений, новые знания о несуществующем объекте моделирования, дополняющие отсутствующую о нем информацию, образуется некоторый мысленный его образ или другими словами, концептуальная модель. В данном случае это концептуальная модель многоуровневой ЛС. Именно концептуальную модель ЛС приходится использовать в дальнейшем в качестве объекта моделирования. Следовательно, задача оценки адекватности аналитической модели управления запасами МР в ЛС реальной ЛС, сводится к задаче оценки ее адекватности концептуальной модели ЛС.

Как было отмечено, концептуальная модель ЛС  представляет собой совокупность представлений о свойствах объекта и его поведении в различных условиях под влиянием факторов самой различной природы. Однако, не всегда представляется возможным выяснить численные значения интересующих параметров системы, необходимых для оценки адекватности, в частности, показателей эффективности и затрат в тех или иных условиях функционирования ЛС, чтобы затем сравнить со значениями тех же параметров, только предоставляемых аналитической моделью в тех же условиях.

Для решения данной проблемы, предлагается заменить объект моделирования (концептуальную модель ЛС) его имитационной моделью. Известно, что в настоящее время по единодушному мнению многих специалистов в области моделирования сложных систем, к примеру [1, 2, 3, 4, 5, 6], метод имитационного моделирования предоставляет исследователю теоретически неограниченные возможности по учету всевозможных факторов, т.е. по реализации алгоритма функционирования моделируемой системы, в отличие от метода аналитического моделирования, в котором можно определить только явные зависимости, реализуемые в виде функциональных зависимостей. В свою очередь, недостатком метода имитационного моделирования, помимо неточности получаемых результатов, зависящей, в принципе, от количества прогонов модели, является также трудность, а иногда и невозможность ее исследования с целью параметрического синтеза системы (в данном случае оптимизации ее параметров) с требуемой точностью, распространяющегося на всю область факторного пространства. Однако в рассматриваемом случае данная задача решается при помощи аналитической модели, от имитационной же модели требуется лишь учет неучтенных в аналитической модели обстоятельств, а также получение необходимых первичных результатов, т.е. значений частных показателей эффективности, используемых впоследствии для оценки адекватности аналитической модели. К неучтенным обстоятельствам можно, в частности, отнести особенности закона распределения коэффициентов потерь, принятых в аналитической модели в виде детерминированных параметров.

 

Постановка задачи

Таким образом, в ходе оценки адекватности аналитической модели необходимо учесть степень влияния неучтенных случайных факторов на частные показатели эффективности функционирования ЛС, путем сравнительно анализа с результатами имитационного моделирования. В данной задаче имитационная модель заменяет реальный объект моделирования, т.е. реальную ЛС.

 

Основная часть

Оценку адекватности разработанной методики исследования предлагается производить на примере ЛС, состоящей из одного склада первого уровня и трех складов второго уровня с соответствующим количеством потребителей. Такое упрощение структуры моделируемого объекта предположительно не окажет существенного влияния на результаты оценки адекватности, поскольку структура системы определяет характер влияния потребности потребителей, обеспечиваемых в первую очередь, на размер дефицита и затрат у потребителей, обеспечиваемых во вторую и последующую очереди, а для отражения данного свойства системы трех потребителей вполне достаточно.

Имитационная модель (ИМ) выполнена в виде самостоятельного приложения выполняющегося под управлением операционной системы семейства Windows. Исходный код программы описан на языке Object Pascal в среде Delphi 7. Центральным элементом программы является модуль, в котором описан класс объекта реализующего в себе имитационную модель.  В методы данного класса входят:

­                    процедура, реализующая моделирующий алгоритм, т.е. имитирующий процесс функционирования (ЛС). В данной процедуре используются значения случайных параметров модели, и по их значениям производится расчет значений дефицита и суммарных затрат;

­                    функции, генерирующие случайные значения потребности в МР и их потери, распределенные соответственно по нормальному и по показательному закону с заданными параметрами. Для формирования нормально распределенной случайной величины используется метод сложения 20-ти равномерно распределенных случайных величин с одинаковыми параметрами, а для формирования случайной величины распределенной по показательному закону используется  метод обратной функции [3, 7, 4].

­                    процедура, осуществляющая циклический прогон модели с накоплением статистики и подсчетом средних значений дефицита и затрат. Количество прогонов не ограничивается.

Состояние процесса имитации функционирования ЛС, а также значения выходных параметров модели отображаются в главном окне программы. Таким образом, в результате обработки программой исходных данных (входных параметров модели), пользователь определяет средние значения дефицита и затрат  возникающих в ЛС за планируемый период.

 Аналитическая модель, адекватность которой оценивалась, была разработана в соответствии с рассмотренной в [9] методикой для случая двухуровневой ЛС обеспечивающей трех потребителей. В частности, ЛС предлагается логически разделить на логистические цепи (линейная цепь складов с запасами МР между которыми осуществляется продвижение материального потока к потребителям) по количеству потребителей. При таком подходе, дефицит  и затраты  во всей системе складываются из дефицита и затрат возникающих в каждой логистической цепи т.е.

 

                                                                                                                      (1)

 

 

где  – математическое ожидание дефицита, возникающее при функционировании n-ой логистической цепи,

 – математическое ожидание затрат на подвоз МР (движение материального потока) между складами n-ой логистической цепи,

N – общее количество потребителей.

 

Дефицит у n-ого потребителя, при :

 

                        ,                                (2)

 

и затраты (затраты в n-ой логистической цепи)

 

       (3)

 

где  – потребность в течение периода T;

 – плотность распределения потребности;

 – запасы МР на складах второго уровня;

запас МР на складе первого уровня для удовлетворения потребности у n-ого потребителя;

 плотность распределения остатка МР на складе первого уровня оставшегося для обеспечения n-ого потребителя;

 – коэффициенты потерь, возникающие в процессе хранения запасов;

 – коэффициенты потерь, возникающие в процессе подвоза МР от склада второго уровня к потребителю;

 – коэффициенты потерь, возникающие в процессе подвоза МР от склада первого к складу второго уровня.

Здесь плотность распределения остатка  определяется с учетом следующего выражения:

 

     (4)

 

 

где   – остаток на складе системы для предыдущего потребителя по отношению к данному. Если n = 2 то . Под  подразумевается первоначальный запас на складе системы. Индекс n для каждого потребителя определяется исходя из очередности обеспечения установленной ЛПР применительно к сложившейся обстановки т.е. под n-ым потребителем подразумевается потребитель обеспечиваемый в n-ую очередь.

Наиболее важным моментом при расчете выходных параметров аналитической модели является определение плотности распределения остатка на складе верхнего уровня. Так как распределение потребности подчинено нормальному закону, то и остаток на складе соединения является также нормально распределенной величиной. Для определения параметров закона распределения остатка, использовались формулы, приведенные в работе [7]. Согласно этих формул, параметры закона распределения остатка для второго

 

                                                                                                    (5)

 

и для третьего потребителей

 

 

                                                                                                   (6)

 

 

где   – математические ожидания остатка на складах ЛС соответственно для второго и третьего потребителей;

– средние квадратические отклонения этих остатков;

 – математические ожидания расходов МР соответственно у первого и второго потребителей;

– средние квадратические отклонения этих расходов.

Плотности распределения остатка на складах будут иметь вид

 

                                           ,                                                   (7)

 

                                           .                                                  (8)

 

Все расчеты, относительно определения численных значений выходных параметров модели, производились с использованием специализированного математического пакета прикладных программ – Maple 9 компании Waterloo Maple.

В качестве меры адекватности (идентичности) аналитической и имитационной модели был выбран коэффициент корреляции между значениями выходных параметров данных моделей [1].

Для оценки коэффициента корреляции необходимо наличие результатов ряда испытаний обеих моделей. Практика показывает, что, как правило, рациональное количество должно составлять порядка десяти испытаний [8]. При проведении испытания, значения входных параметров обеих моделей в каждом опыте естественно должны совпадать. В каждом испытании выбирался уникальный (не повторяющийся) набор входных параметров.

Коэффициент корреляции случайных величин X и Y определяется по формуле [7]:

 

                                                     ,                                                             (9)

 

где  – математическое ожидание произведения случайных величин X и Y;

 – математические ожидания случайных величин X и Y соответственно;

 – средние квадратические отклонения случайных величин X и Y соответственно.

В данном случае в качестве X и Y выступают две пары случайных величин – дефицита ,  и затрат , . Таким образом оценку адекватности моделей необходимо производить по значениям двух коэффициентов корреляции – дефицита и затрат.

Для проведения практических расчетов использовались инструменты статистического анализа пакета Microsoft Excel.

В ходе расчетов были получены следующие значения коэффициентов корреляции для дефицитов и затрат:

 

,

.

 

Выводы

Таким образом, в результате оценки адекватности аналитической модели управления запасами в многоуровневой ЛС объекту моделирования с использованием рассмотренного в статье подхода, можно сделать вывод, что аналитическая модель адекватна имитационной, а значит, учитывая ранее принятые допущения, и самому объекту моделирования.

Используемый подход к оценке адекватности аналитических моделей заключается в сопоставлении данных полученных в результате аналитического моделирования с результатами имитационного моделирования, считая при этом имитационную модель рассматриваемой ЛС абсолютно адекватной самой моделируемой системе. Данный подход обусловлен физическим отсутствием, по объективным причинам, объекта моделирования, по сравнению с которым возможно оценить адекватность аналитической модели. В тоже время неограниченные возможности имитационного моделирования относительно учета всевозможных свойств моделируемого объекта, чего нельзя сделать в аналитических моделях, позволили заменить рассматриваемый объект исследования его имитационной моделью.

Рассмотренный подход применим к оценке адекватности аналитических моделей не только применительно к ЛС но и к другим сложным объектам моделирования  непосредственное проведение экспериментов с которыми, не возможно по тем или иным причинам.

 

In the article one of possible approaches is presented to estimation of adequacy of analytical model to control by the supplies of material resources in the multilevel logistic system.

 

1.        Молчанов А.А. Моделирование и проектирование сложных систем. Учебное пособие. – К.: Вища школа. 1988. – 56 с.

2.        Бусленко Н.П. Моделирование сложных систем.  – М.: Наука, 1968. – 352 с.

3.        Максимей В.М. Имитационное моделирование  на ЭВМ. – М.: Радио и связь, 1988. – 223 с.

4.        Советов Б.Я., Яковлев С.А. Моделирование систем. Учебник. – М.: Высшая школа. 1985. – 271 с.

5.        Шарашканэ А.С., Железнов Н.Г., Иваницкий В.А. Сложные системы. – М.: Высшая школа, 1977. – 229 с.

6.        Нейлор Т. Машинный имитационный эксперимент с моделями. – М.: Мир. 1975.

7.        Вентцель Е.С. Овчаров Л.А. Теория вероятности и ее инженерное приложение. – М.: Наука, 1988. – 358 с.

8.        Вентцель Е.С. Теория вероятности. – М.: Наука, 1969.

9.        Козак Ю.А., Колчин Р.В. Математическое обеспечение АСУ запасами материальных ресурсов в двухуровневой логистической системе со случайным спросом. // Збірник наукових праць №28.Одеса: ОНПУ, 2004.

 





Ответы на вопросы [_Задать вопроос_]

Моделирование объектов и систем управления

Соколов А.Е., Махова Е.О. Моделирование процесса принятия педагогического решения при компьютеризированном обучении

Славко О.Г. Порівняльний аналіз керування регулятором на основі локальної моделі керованого процесу та П-регулятором

Войтенко В.В., Дикусар Е.В, Ситников В.С. Определение частоты среза устройства сглаживания данных на основе метода скользящего среднего

Передерій В.І. Алгоритм визначення та оцінки характеристик ефективності комп’ютерних систем на початковій стадії проектування в умовах невизначенності

Ляшенко С.А, Ляшенко А.С. Оценка модели псевдолинейной регрессии

Ладієва Л.Р. Математична модель процесу газової мембранної дистиляції

Носов П.С., Косенко Ю.І. Нечіткі моделі і методи ідентифікації та прогнозу стану інформаційної моделі студента

Китаев А.В., Глухова В.И. Анализ работы синхронного двигателя с неявнополюсным ротором по данным каталога

Дорошкевич В.К., Пироженко А.В., Хитько А.В., Хорольский П.Г. К определению требований к системам увода космических объектов

Голінко І.М., Ковриго Ю.М., Кубрак А.І. Настройка системи керування за імпульсною характеристикою об’єкта

Яшина К.В., Садовой А.В. Комплексная математическая модель тепловых процессов, происходящих в дуговых электросталеплавильных печах

Шейник С.П., Рудакова А.В. Использование функций принадлежности для моделирования параметров распределенных объектов

Хомченко А.Н., Литвиненко Е.И. Метод барицентрического усреднения граничных потенциалов электростатического поля

Селяков Е. Б. Моделирование требований к техническим системам методами математической логики

Тодорцев Ю.К., Ларіонова О.С., Бундюк А.М. Математична модель контура теплопостачання когенераційної енергетичної установки

Кириллов О.Л. , Якимчук Г.С. Моделирование процесса управления системой перегрузки углеводородных жидких топлив

Шеховцов А.Н., Козел В.Н. Построение математической модели формирования распределенных систем

Китаев А.В., Глухова В.И. Анализ поведения генератора постоянного тока по данным каталога

Хомченко А.Н., Козуб Н.О. Задачі наближення функцій: від лагранжевих до серендипових поліномів

Хобин В.А., Титлова О.А. Определение температуры парожидкостной смеси в дефлегматоре АДХМ по результатам измерений температуры его поверхности

Григорова Т.М., Усов А.В. Вероятностно-статистическое моделирование маршрутизированных пассажиропотоков в крупных городах

Горач О.О., Тернова Т.І. Моделювання технологічного процесу одержання трести при використані штучного зволоження з урахуванням складу мікрофлори

Дубік Р.М., Ладієва Л.Р. Математична модель розділення неоднорідних рідких систем

Казак В.М, Лейва Каналес Родриго, Яковицкая Е.Ю. Моделирование динамики полета магистрального самолета на исследовательском стенде

Завальнюк И.П. Исследование процесса торможения автомобиля как критического режима динамической системы

Дмитриев С.А., Попов А.В. Построение портрета неисправностей проточной части газотурбинного двигателя на примере АИ-25

Русанов С.А., Луняка К.В., Клюєв О.І., Глухов Г.М. Математичне моделювання робочого процесу в апаратах з віброкиплячим шаром та розробка систем автоматизованого моделювання гідродинаміки віброкиплячих шарів

Боярчук В.П., Сыс В.Б. Экспериментальные исследования влияния технологии шлихтования на изменение жесткости текстильных нитей

Селін Ю.М. Використовування контекстних марківських моделей для аналізу дії промислових вибухів на будівельні конструкції

Рудакова А.В. Проблемы интеграции сложных систем

Передерій В.І., Касап А.М. Математична модель та алгоритм автоматизації розрахунку параметрів комп’ютеризованих систем працюючих у реальному часі

Передерий В.И., Еременко А.П. Математические модели и алгоритмы принятия релевантных решений пользователями автоматизированных систем с учетом личностных и внешних факторов на базе генетических алгоритмов

Михайловская Т.В., Михалев А.И., Гуда А.И. Исследование правил клеточных автоматов для моделирования процессов затвердевания квазиравновесных бинарных сплавов

Хомченко А.Н., Колесникова Н.В. Явление «сверхсходимости» в задаче Прандтля для уравнения Пуассона

Китаев А.В., Глухова В.И. Анализ работы трансформатора по данным каталога

Квасницкий В.В., Ермолаев Г.В., Матвиенко М. В., Бугаенко Б.В., Квасницкий В.Ф. Оценка применимости метода компьютерного моделирования к исследованию напряженно-деформиррованного состояния цилиндрических узлов

Китаев А.И., Глухова В.И. Анализ работы асинхронного двигателя по данным каталога

Шелестов А.Ю Имитационная модель взаимодействия GRID-узлов с очередью доступа к общей памяти

Chizhenkova R.A. Mathematical Aspects of Bibliometrical Analysis of Neurophysiological Investigations of Action of Non-ionized Radiation (Medline-Internet)

Хомченко А.Н., Козуб Н.А. Геометрическое моделирование дискретных элементов с криволинейными границами

Славич В.П. Модель автоматизованої системи управління потоками транспортних засобів

Маркута О.В., Мысак В.Ф. Программная реализация и исследование особенностей метода группового учета аргументов

Степанкова Г.А., Баклан І.В. Побудова гібридних моделей на основі прихованих марківських моделей та нейронних мереж

Бакшанська Т.Д., Рижиков Ю.Г., Тодорцев Ю.К. Математична модель процесу горіння природного газу з рециркуляцією продуктів згорання для цілей управління

Хомченко А.Н. Новые решения обобщенной задачи Бюффона

Передерий В.И., Еременко А.П. Математические модели и алгоритмы определения релевантности принимаемых решений с учетом психофункциональных характеристик пользователей при управлении автоматизированными динамическими системами

Ложечников В.Ф., Михайленко В.С., Максименко И.Н. Аналитическая много режимная математическая модель динамики газовоздушного тракта барабанного котла средней мощности

Ковриго Ю.М., Фоменко Б.В., Полищук И.А. Математическое моделирование систем автоматического регулирования с учетом ограничений на управление в пакете Matlab

Исаев Е.А., Наговский Д.А. Математическое описание влияния кривизны контактирующих тел на угол смачивания жидкости в межчастичном пространстве

Бідюк П.І., Литвиненко В.І., Кроптя А.В. Аналіз ефективності функціонування мережі Байєса

Тищенко И.А., Лубяный В.З. Математическое моделирование вокодера для определения оптимальной формы импульса сигнала возбуждения.

Николаенко Ю.И., Моисеенко С.В. Моделирование гармонического полиномиального базиса гексагона.

Козуб Н.А., Манойленко Е.С., Хомченко А.Н. Температурный тест для модифицированных базисов бикубической интерполяции.

Клименко А.К. Об упрощенном численном конструировании обратной модели динамического объекта.

Китаев А.В., Сушич Е.Ф. Расчет погрешностей измерительных трансформаторов.

Передерій В.І.,Касап А.М. Математична модель та алгоритм автоматизації розрахунку параметрів комп’ютеризованих систем працюючих у реальному часі

Шпильовий Л.В. Математична модель та алгоритм екстремального управління процесом осадження дисперсної фази суспензії.

Тулученко Г.Я. Інформаційний модуль експрес-пошуку точок еквівалентності процесу нейтралізації.

Тернова Т.І. Урахування морфогенетичного рівняння в математичній моделі тканини.

Попруга А.Г. Теоретические и экспериментальные исследования электрических нагревателей по критерию экономии энергии.