Опубликован: 01.09.2009 | Уровень: для всех | Доступ: платный
Лекция 4:

Организационно-экономическая система управления материальными запасами промышленных корпоративных систем

Аннотация: Данная лекция посвещена организационно-экономической системе управления материальными запасами промышленных корпоративных систем. Рассмотрены такие разделы, как Базовые модели управления запасами; Теоретические положения создания организационно-экономической системы управления запасами; даны определения, приведены примеры.
Ключевые слова: ПО, оптимизация, принятия решений, товар, интервал, функция, плата, дифференциал, затраты, оптимальный план, расходы, запись, область определения, стоимость, предел, график, графика, минимизация, натуральное число, Гипербола, минимум, целое число, производная функции, алгоритм, значение, площадь, высота, множества, отрезок, неравенство, оптимальность, отношение, тождество, равенство, Приращение, устойчивость, точность, погрешность, определение, индикатор, обобщение, обобщенная функция, параметр, множитель, Абсолютной погрешностью, опыт, математическим ожиданием, случайная величина, вывод, анализ, моделирование, модель оптимизации, принципиальная схема, деятельность, финансовые инструменты, информация, эффективное управление, работ, стабильность, ликвидность, краткосрочные обязательства, заемные источники, коэффициенты, поток, пространство, место, очередь, торговля, длина, критерий эффективности, стохастическая модель, вероятность, ресурс, вес, класс, оптимальная стратегия, нижняя граница, сравнительные оценки, приемка товара, метод решения, объект, состояние системы, пошаговая оптимизация, функция выигрыша, максимум, оптимальное решение задачи, управление системой, размерность, множество состояний, блок-схема, Типовая, e-procurement, inventory management, operations management, решение игры, операционный, менеджмент, цикла, группа, service support, автоматизация, регулирование, контроль, математическая модель, подразделения, IDEF0, диаграмма, декомпозиция, выход, инициация, постановка задачи

4.1. Базовые модели управления запасами

Управление запасами (материальнотехническое снабжение) - неотъемлемая часть работы фирм и организаций. Речь идет о запасах сырья, топлива, материалов, инструментов, комплектующих изделий, полуфабрикатов, готовой продукции на промышленном (или сельскохозяйственном) предприятии, о запасах товаров на оптовых базах, складах магазинов, на рабочих местах продавцов, наконец, у потребителей. Запасы постоянно расходуются и пополняются по тем или иным правилам, принятым на предприятии. Оптимизация этих правил, т. е. оптимальное управление запасами, дает большой экономический эффект.

Математическая теория управления запасами - крупная область экономико-математических исследований, получившая свое развитие начиная с пятидесятых годов ХХ века. Предложенная еще в 1915 г. Ф. Харрисом классическая модель теории управления запасами, называемая также моделью Вильсона (поскольку получила известность после публикации работы Р.Г. Вильсона в 1934 г.) - один из наиболее простых и наглядных примеров применения математического аппарата для принятия решений в экономической области. Формула оптимального размера заказа, полученная в модели Вильсона, широко применяется на различных этапах производства и распределения продукции. Она практически полезна для принятия решений при управлении запасами, в частности, приносит заметный экономический эффект [11]. Рассмотрим эту модель подробнее.

Классическая модель управления запасами.Пусть y(t) - величина запаса некоторого товара на складе в момент времени t, t \ge 0. Дефицит не допускается, т. е. y(t) \ge 0 при всех t. Товар пользуется равномерным спросом с интенсивностью \mu, т. е. за интервал времени \Delta t со склада извлекается и поступает потребителям часть запаса величиной \mu \Delta t . В моменты времени t_{0} = 0, t_{1}, t_{2},\dots пополняется запас на складе - приходят поставки величиной Q_{0}, Q_{1}, Q_{2},\dots соответственно. Таким образом, изменение во времени величины запаса y(t) товара на складе изображается зубчатой ломаной линией (рис. 4.1), состоящей из наклонных и вертикальных звеньев, причем наклонные отрезки параллельны.

График изменения величины y(t) запаса на складе в момент времени t

Рис. 4.1. График изменения величины y(t) запаса на складе в момент времени t

Таким образом, в момент t_{i} величина запаса на складе y(t) скачком увеличивается на Q_{i}. Следовательно, функция y(t) имеет разрывы в точках t_{1}, t_{2},\dots Для определенности будем считать, что эта функция непрерывна справа.

Пусть s - плата за хранение единицы товара в течение единицы времени. Поскольку можно считать, что величина запаса y(t) не меняется в течение интервала времени (t; t+dt), где dt - дифференциал, т. е. бесконечно малая, то плата за хранение всего запаса в течение этого интервала времени равна sy(t)dt. Следовательно, затраты за хранение в течение интервала времени [0;T), где T - интервал планирования, пропорциональны (с коэффициентом пропорциональности s ) площади под графиком уровня запаса на складе y(t) и равны

s \int_{0}^{T}{y(t)dt}

Пусть g - плата за доставку одной партии товара. Примем для простоты, что она не зависит от размера поставки. Позже покажем, что если эта плата равна g+g_{1}Q, где Q - размер поставки, то оптимальный план поставки - тот же, что и при отсутствии линейного члена. Будет проанализирована и более сложная модель, в которой предусмотрена скидка с ростом поставки, приводящая к выражению g+g_{1}Q+ g_{2}Q^{2} для платы за доставку одной партии товара размером Q.

Пусть n(T) - количество поставок, пришедших в интервале [0;T). При этом включаем поставку в момент t = 0 и не включаем поставку в момент t = T (если такая поставка происходит). Тогда суммарные издержки на доставку товара равны gn(T). Следовательно, общие издержки (затраты, расходы) за время T равны

F(T; y) = F(y(t), 0 \le t \le T) = gn(T) + s \int_{0}^{T}{y(t)dt}.

Запись F(T; y) = F(y(t), 0 \le t \le T) означает, что общие издержки зависят от значений функции y = y(t) при всех . Символ у обозначает функцию как целое. То есть область определения F(T;y) при фиксированном T - не множество чисел, а множество функций.

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

f ( T; y) = f ( y ( t ), 0\le t \le T) = \cfrac{1}{T}F(T;y) = \cfrac{1}{T} \left \{gn(T) + s \int_0^{T}{ y(t)dt}  \right \}.

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

Если задать моменты прихода поставок и величины партий, то будет полностью определена функция y = y(t) при всех 0 \le  t \le  T. Верно и обратное - фиксация функции y = y(t), 0 \le  t \le  T рассматриваемого вида (рис. 4.1) полностью определяет моменты прихода поставок и величины партий. И то, и другое будем называть планом поставок или планом работы системы управления запасами. Для ее оптимизации необходимо выбрать моменты времени t_{0} = 0, t_{1}, t_{2},\dots  пополнения запаса на складе и размеры поставляемых партий товара Q_{0}, Q_{1}, Q_{2},\dots  так, чтобы минимизировать средние издержки f_{T}(y) при фиксированном Т. Модель производственной ситуации (т. е. работы склада) описывается четырьмя параметрами - \mu (интенсивность спроса), s (стоимость хранения единицы продукции в течение единицы времени), g (стоимость доставки партии товара), Т (горизонт планирования).

Решение задачи оптимизации.Поставленная задача оптимизации работы склада интересна тем, что неизвестно число 2n(T)-1 параметров, определяющих план поставок. Поэтому ее решение не может быть проведено с помощью стандартных методов теории оптимизации.

Решим эту задачу в три этапа. На первом установим, что оптимальный план следует искать среди тех планов, у которых все зубцы доходят до оси абсцисс, т. е. запас равен 0 в момент доставки очередной партии. Цель второго этапа - доказать, что все зубцы должны быть одной и той же высоты. Наконец, на третьем находим оптимальный размер поставки.

Оптимальный план.Найдем наилучший план поставок. План, для которого запас равен 0 (т. е. y(t) = 0 ) в моменты доставок очередных партий, назовем напряженным.

Утверждение 1.Для любого плана поставок, не являющегося напряженным, можно указать напряженный план, для которого средние издержки меньше.

Покажем, как можно от произвольного плана перейти к напряженному плану, уменьшив при этом издержки. Пусть с течением времени при приближении к моменту t _{1} прихода поставки Q_1 уровень запаса не стремится к 0, а лишь уменьшается до y(t _{1}-) \ne 0 (где знак "минус" означает предел слева функции y(t) в точке t _{1} ). Тогда рассмотрим новый план поставок с теми же моментами поставок и их величинами, за исключением величин поставок в моменты t = 0 и t = t_{1}. А именно, заменим Q_{0} на Q_{01} = Q_{0} y(t_{1}-), а Q_{1} на Q_{11} = Q_{0} + y(t_{1}-). Тогда график уровня запаса на складе параллельно сдвинется вниз на интервале (0; t _{1}), достигнув 0 в t_{1}, и не изменится правее точки t_{1}. Следовательно, издержки по доставке партий не изменятся, а издержки по хранению уменьшатся на величину, пропорциональную (с коэффициентом пропорциональности s ) площади параллелограмма, образованного прежним и новым положениями графика уровня запаса на интервале (0; t _{1}) (см. рис. 4.2).

Первый шаг перехода к напряженному плану

увеличить изображение
Рис. 4.2. Первый шаг перехода к напряженному плану

В результате первого шага перехода получен план, в котором крайний слева зубец достигает оси абсцисс. Следующий шаг проводится аналогично, только момент времени t = 0 заменяется на t = t_{1}. Если есть такая возможность, второе наклонное звено графика уровня запаса на складе параллельно сдвигается вниз, достигая в крайней правой точке t_{2 } оси абсцисс.

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

Утверждение 1.Оптимальный план следует искать только среди напряженных планов, т. е. план, не являющийся напряженным, не может быть оптимальным.

Утверждение 2.Среди напряженных планов с фиксированным числом поставок минимальные издержки имеет тот, в котором все интервалы между поставками равны.

При фиксированном числе поставок затраты на доставку партий не меняются. Следовательно, достаточно минимизировать затраты на хранение.

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

Q_{i-1} = \mu (t_{i} -t_{i-1})  i = 1, 2,\dots , n(T)-1, Q_{n(T)-1} = \mu (T- t_{n(T)-1}).

Действительно, очередная поставка величиной Q_{i-1} совпадает с размером запаса на складе в момент t_{i-1}, расходуется с интенсивностью \mu единиц товара в одну единицу времени и полностью исчерпывается к моменту t_{i} прихода следующей поставки.

Для напряженного плана издержки по хранению равны

s \int_0^T{y(t)dt} = s \sum_{i=1}^{n(T)}{\cfrac{Q_{i-1}(t_i - t_{i-1})}{2}} = 
s \sum_{i=1}^{n(T)}{\cfrac{\mu (t_i - t_{i-1})^2}{2}} = 
s \sum_{i=1}^{n(T)}{\cfrac{\mu \Delta_i^2}{2}} = 
\cfrac{\mu s}{2} \sum_{i=1}^{n(T)}{ \Delta_i^2 }

где \Delta _{i} = t_{i} - t_{i-1}, i=1,2,\dots , n(T), 
t_{n(T)} = T.

Ясно, что \Delta_i, \, i=1, 2, \dots ,n(T) - произвольные неотрицательные числа, в сумме составляющие Т. Следовательно, для минимизации издержек среди напряженных планов с фиксированным числом поставок достаточно решить задачу оптимизации

\left \{
\begin{array}{l}
\Delta_{1}^{2} +\Delta_2^{2} +\dots  + \Delta^{2}_{n} \to min, \\
\Delta_{1}+\Delta _{2}+\dots  + \Delta_{n}=T, \\
\Delta_i \ge 0, i = 1, 2,\dots ,n,
\end{array}

где n = n(T).

Полученная задача оптимизации формально никак не связана с логистикой, она чисто математическая. Для ее решения целесообразно ввести новые переменные \alpha_i = \Delta_i - \cfrac{T}{n}, i = 1,2,\dots , n.

Тогда

\sum_{i=1}^{n}{\alpha_i} = 
\sum_{i=1}^{n}{\Delta_i - \cfrac{T}{n}}
= \sum_{i=1}^{n}{\Delta_i} -n\cfrac{T}{n} = T-T = 0

Поскольку \Delta_i = \cfrac{T}{n} + \alpha_i то \Delta_i^2 = \cfrac{T^2}{n^2} + 2 \cfrac{T}{n} \alpha_i + \alpha_i^2, следовательно, с учетом предыдущего

\sum_{i=1}^{n}{ \Delta_i^2} = n\cfrac{T^2}{n^2} + 2 \cfrac{T}{n} \sum_{i=1}^{n}{\alpha_i} + \sum_{i=1}^{n}{\alpha_i^2} = \cfrac{T^2}{n} + \sum_{i=1}^{n}{\alpha_i^2}

Сумма квадратов всегда неотрицательна. Она достигает минимума, равного 0, когда все переменные равны 0, т. е. при \alpha_1 = \alpha_2 = \dots = \alpha_n = 0. Тогда

\Delta_i = \cfrac{T}{n}, \, i = 1, 2, \dots, n

При этих значениях \Delta_{i} выполнены все ограничения оптимизационной задачи. Итак, утверждение 2 доказано.

Для плана с равными интервалами между поставками все партии товара имеют одинаковый объем. Для такого плана издержки по хранению равны

s \int_{0}^{T}{y(t)dt}=\cfrac{\mu s}{2} \sum_{i=1}^{n(T)}{ \Delta_i^2 }=
\cfrac{\mu s T^2}{n(T)} ( T)})

Средние издержки (на единицу времени) таковы:

f(T;y) = \cfrac{1}{T}\left \{ gn(T) + \cfrac{\mu s T^2}{n(T)}\right \}
= g \cfrac{n(T)}{T} + \mu s \cfrac{T}{2n(T)}. ( T)})

Итак, минимизация средних издержек - это задача дискретной оптимизации. На третьем этапе построения оптимального плана необходимо найти натуральное число n(T) - самое выгодное число поставок.

Поскольку к моменту Т запас товара должен быть израсходован, то общий объем поставок за время T должен совпадать с общим объемом спроса, следовательно, равняться \mu T. Справедливо балансовое соотношение (аналог закона Ломоносова-Лавуазье сохранения массы при химических реакциях):

Qn(T) = \mu T.

Из балансового соотношения следует, что

\cfrac{n(T)}{T} = \cfrac{\mu}{Q}

Средние издержки (на единицу времени) можно выразить как функцию размера партии

f(T;y) = g \cfrac{n(T)}{T} + \mu s \cfrac{T}{2n(T)} = 
f_1(Q) = \cfrac{\mu g}{Q} + \cfrac{sQ}{2} ( 4.1)

Задача состоит в минимизации f _{1}(Q) по Q. При этом возможная величина поставки принимает дискретные значения, Q \in \left \{ \cfrac{\mu T}{n}  n = 1, 2,\dots  \right \}.

Изучим функцию f _{1}(Q ), определенную при Q\ge 0. При приближении к 0 она ведет себя как гипербола, при росте аргумента - как линейная функция. Производная имеет вид

\cfrac{df_1(Q)}{dQ}= - \cfrac{\mu g}{Q^2} + \cfrac{s}{2}

Производная монотонно возрастает, поэтому рассматриваемая функция имеет единственный минимум в точке, в которой производная равна 0, т. е. при

Q_0 = \sqrt{\cfrac{2\mu g}{s}} ( 4.2)

Получена знаменитая "формула квадратного корня".

В литературе иногда без всяких комментариев рекомендуют использовать напряженный план, в котором размеры всех поставляемых партий равны Q_0. К сожалению, получаемый таким путем план почти всегда не является оптимальным,т. е. популярная рекомендация неверна или не вполне корректна. Дело в том, что почти всегда

Q_0 \notin \left \{ \cfrac{\mu T}{n}  n = 1, 2,\dots  \right \}

Всегда можно указать неотрицательное целое число n такое, что

Q_1 = \cfrac{\mu T}{n +1} < Q_0 \le \cfrac{\mu T}{n} = Q_2 ( 4.3)
Михаил Агапитов
Михаил Агапитов

Не могу найти  требования по оформлению выпускной контрольной работы по курсу профессиональной переподготовки "Менеджмент предприятия"

Подобед Александр
Подобед Александр

Я нажал кнопку "начать курс" и почти его уже закончил, но для получения диплома на бумаге, нужно его же оплатить? Как оплатить? 

Александр Кулаков
Александр Кулаков
Россия, Челябинск
Максим Овцин
Максим Овцин
Россия, Санкт-Петербург, НИУ ИТМО