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

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

Утверждение 3.Решением задачи оптимизации

f_1(Q) = \cfrac{\mu g}{Q} + \cfrac{sQ}{2} \to min \\
Q \in \left \{ \cfrac{\mu T}{n}  n = 1, 2,\dots  \right \}

является либо Q _{1}, либо Q _{2}.

Действительно, из всех Q \in \left \{ \cfrac{\mu T}{n}  n = 1, 2,\dots  \right \} часть лежит правее Q_0, из них наименьшим является Q _{2}, а часть лежит левее Q_0, из них наибольшим является Q _{1}. Для построения оптимального плана обратим внимание на то, что производная функции f _{1}(Q ) отрицательна левее Q_0 и положительна правее Q_0, следовательно, функция средних издержек f _{1}(Q ) убывает левее Q_0 и возрастает правее Q_0. Значит, минимум по Q \in \left \{ \cfrac{\mu T}{n}  n = 1, 2,\dots  \right \} \cap (Q : Q \ge Q_0) достигается при Q = Q_{2}, а минимум по Q \in \left \{ \cfrac{\mu T}{n}  n = 1, 2,\dots  \right \} \cap (Q : Q < Q_0) - при Q = Q_1

Последнее утверждение эквивалентно заключению утверждения 3.

Итак, алгоритм построения оптимального плана таков.

  1. Найти Q_{0} по формуле квадратного корня (4.2).
  2. Найти n из условия (4.3).
  3. Рассчитать f_{1}(Q) по формуле (4.1) для Q = Q_{1} и Q = Q_{2}, где Q_{1} и Q_{2} определены в (4.3).
  4. Наименьшее из двух чисел f _{1}(Q _{1}) и f _{1}(Q_{2}) - искомый минимум, а то из чисел Q_{1} и Q_{2}, на котором достигается минимум, - решением задачи оптимизации. Обозначим его Q_{opt}.

Итак, оптимальный план поставки - это напряженный план, в котором объемы всех поставок равны Q_{opt}.

Замечание.Если f _{1}( Q _{1}) = f1(Q_{2}), то решение задачи оптимизации состоит из двух точек Q_{1} и Q_{2}. В этом частном случае существует два оптимальных плана.

Пример 1.На складе хранится некоторая продукция, пользующаяся равномерным спросом. За 1 день со склада извлекается 5 т продукции. Плата за хранение 1 т продукции в день - 50 руб. Плата за доставку одной партии - 980 руб. Горизонт планирования - 10 дней. Найти оптимальный план поставок.

В рассматриваемом случае \mu =5 (т/день), s = 50 (руб./тдень), g = 980 (руб./партия), Т = 10 (дней). По формуле (4.2) рассчитываем

Q_0 = \sqrt{\cfrac{2\mu g}{s}} = 
\sqrt{\cfrac{2 \cdot 5 \cdot 980}{50}} = \sqrt{196} = 14

Множество допустимых значений для Q имеет вид

\left \{ \cfrac{\mu T}{n}  n = 1, 2,\dots  \right \} = 
\left \{50; \cfrac{50}{2}; \cfrac{50}{3}; \cfrac{50}{4}; \dots  \right \} 
= \{50; 25; 16,67; 12,5;\dots \}

Следовательно, Q_{1} = 12,5 и Q_{2}= 16,67. Первое значение определяет напряженный план с четырьмя одинаковыми зубцами, а второе - с тремя. Поскольку

f_1(Q) = \cfrac{5\cdot 980}{Q} + \cfrac{50Q}{2} = \cfrac{4900}{Q} + 25Q

то

f_{1} (Q _{1}) = f _{1}(12,5) = \cfrac{4900}{12,5}+ 25\cdot 12,5 = 392 + 312,5 = 704,5

и

f _{1}( Q _{2}) = 
f _{1}\left ( \cfrac{50}{3}\right ) = 
\cfrac{4900 \cdot 3}{50} + 25\cdot \cfrac{50}{3} = 294 + 416,67 = 710,67.

Поскольку f1(Q1) < f _{1}(Q_{2}), то Q_{opt} = Q_1 = 12,5. Итак, оптимальным является напряженный план с четырьмя зубцами.

Как уже отмечалось, часто рекомендуют применять план поставок с Q=Q_{0}. Каков при этом проигрыш по сравнению с оптимальным планом?

Для плана с Q = Q_{0} интервал между поставками составляет \cfrac{Q_0}{\mu} = \cfrac{14}{5} = 2,8 \text{ дня}. Следовательно, партии придут в моменты t_{0} = 0 ; t _{1} = 2,8 ; t_{2} = 5,6 ; t_{3} = 8,4. Следующая партия должна была бы прийти уже за пределами горизонта планирования Т =10, в момент t_{4} = 11,2. Таким образом, график уровня запаса на складе в пределах горизонта планирования состоит из трех полных зубцов и одного неполного. К моменту Т =10 пройдет 10 - 8,4 = 1,6 дня с момента последней поставки, значит, со склада будет извлечено 5 \cdot 1,6 = 8 т продукции и останется 14 - 8 = 6 т. План с Q = Q_{0} не является напряженным, а потому не является оптимальным для горизонта планирования Т =10.

Подсчитаем общие издержки в плане с Q = Q_{0}. Площадь под графиком уровня запаса на складе равна сумме площадей трех треугольников и трапеции. Площадь треугольника равна \cfrac{14\cdot 2,8}{2}= 19,6 трех треугольников - 58,8. Основания трапеции параллельны оси ординат и равны значениям уровня запаса в моменты времени t_{3} = 8,4 и Т = 10, т. е. величинам 14 и 6 соответственно. Высота трапеции лежит на оси абсцисс и равна 10 - 8,4 = 1,6, а потому площадь трапеции есть \cfrac{(14+ 6)\cdot 1,6}{2}= 16. Следовательно, площадь под графиком равна 58,8 + 16 = 74,8, а плата за хранение составляет 50 \cdot 74,8 = 3740 \text{ руб}.

За 10 дней доставлены 4 партии товара (в моменты t_{0} = 0 ; t _{1} = 2,8 ; t_{2} = 5,6 ; t_{3} = 8,4 ), следовательно, затраты на доставку равны 4 \cdot 980 = 3920 руб. Общие издержки за 10 дней составляют 3740+3920 = 7660 руб., а средние издержки - 766 руб. Они больше средних издержек в оптимальном плане в 766 / 704,5 = 1,087 раза, т. е. на 8,7 \%.

Отметим, что

f_1(Q_0) = \cfrac{4900}{Q_0} + 25Q_0 = \cfrac{4900}{14} + 25 \cdot 14 = 350 +350 = 700

т. е. меньше, чем в оптимальном плане. Таким образом, из-за дискретности множества допустимых значений средние издержки возросли на 4,5 руб., т. e. на 0,64 \%. При этом оптимальный размер партии ( 12,5 т) отличается от Q_{0} = 14 т на 1,5 т, т. е. Q_{opt} / Q_{0} = 0,89 - различие на 11 \%. Достаточно большое различие объемов поставок привело к пренебрежимо малому изменению функции f _{1}(Q). Это объясняется тем, что в точке Q_{0} функция f _{1}(Q) достигает минимума, а потому ее производная в этой точке равна 0.

Оба слагаемых в f _{1}(Q_{0}) равны между собой. Случайно ли это? Покажем, что нет. Действительно,

\cfrac{\mu g}{Q_0} = 
\cfrac{\mu g}{\sqrt{\cfrac{2\mu g}{s}}} = 
\sqrt{\cfrac{\mu g s}{2}}; \, 
\cfrac{sQ_0}{2} = \cfrac{s\sqrt{\cfrac{2\mu g}{s}}}{2} = \sqrt{\cfrac{\mu g s}{2}}

Таким образом, составляющие средних издержек, порожденные различными причинами, уравниваются между собой.

Средние издержки в плане с Q = Q_{0} равны \sqrt{2\mu gs}. Интервал между поставками при этом равен

\cfrac{Q_0}{\mu} = \cfrac{\sqrt{\cfrac{2\mu g}{s}}}{\mu} = \sqrt{\cfrac{2g}{\mu s}}.

Издержки в течение одного интервала между поставками таковы:

\sqrt{2\mu g s}\cdot \sqrt{\cfrac{2g}{\mu s}} = 2g

при этом половина (т. е. g ) приходится на оплату доставки партии, а половина - на хранение товара.

Асимптотически оптимальный план.Из проведенных рассуждений ясно, что напряженный план с Q=Q_{0} является оптимальным тогда и только тогда, когда горизонт планирования Т приходится на начало очередного зубца, т. е. для

T = n \cfrac{Q_0}{\mu} = n \sqrt{\cfrac{2g}{\mu s}}, n = 1, 2, \dots ( 4.4)

Для всех остальных возможных горизонтов планирования Т этот план не является оптимальным. Оптимальным будет напряженный план с другим размером поставки. Для дальнейшего весьма существенно, что при изменении горизонта планирования Т от 0 до Т_{0} оптимальный план меняется на всем интервале [0; T_{0}].

Как происходит это изменение? При малых горизонтах планирования Т делается лишь одна поставка (в момент времени t = 0 ), график уровня запаса на складе состоит из одного зубца. При увеличении Т размер зубца плавно увеличивается. В некоторый момент Т(1) происходит переход от одного зубца к двум. В этот момент оптимальны сразу два плана поставки - с одним зубцом и с двумя. При переходе к планам с двумя зубцами размер зубца скачком уменьшается. При дальнейшем увеличении горизонта планирования оптимальный план описывается графиком с двумя одинаковыми зубцами, размер которых плавно растет. Далее в момент Т(2) становится оптимальным план с тремя зубцами, размер которых в этот момент скачком уменьшается (в компенсацию за увеличение числа скачков). И т. д.

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

Ответ можно указать, если горизонт планирования является достаточно большим. Оказывается можно использовать план, в котором все размеры поставок равны Q_{0}. Для него уровень запаса на складе описывается функцией y_{0}(t), 0 \le  t \le  +\infty , состоящей из зубцов высоты Q_{0}. Предлагается пользоваться планом, являющимся сужением этого плана на интервал [0; T). Другими словами, предлагается на интервале [0; T) использовать начальный отрезок этого плана. Он состоит из некоторого количества треугольных зубцов, а последний участок графика, описываемый трапецией, соответствует тому, что последняя поставка для почти всех горизонтов планирования не будет израсходована до конца. Такой план иногда называют планом Вильсона [11].

Ясно, что этот план не будет оптимальным (для всех Т, кроме заданных формулой (4.4)). Действительно, план Вильсона можно улучшить, уменьшив объем последней поставки. Однако у него есть то полезное качество, что при изменении горизонта планирования его начальный отрезок не меняется. Действительно, планы поставок для горизонтов планирования Т_{1} и Т_{2}, определенные с помощью функции y_{0}(t), 0 \le  t \le  +\infty , задающей уровень запасов на складе, совпадают на интервале [0; min \{Т _{1}, Т_{2}\}).

Определение.Асимптотически оптимальным планом называется план поставок - функция y:[0; +\infty) \to [0; +\infty) такая, что

\lim_{T \to \infty}{\cfrac{f(T; y_{opt}(T))}{f(T;y)}} = 1,

где y_{opt}(T) - оптимальный план на интервале [0; T).

В соответствии с определениями и обозначениями, введенными в начале раздела, f{T; y_{opt}(T)) - средние издержки за время Т для плана y_{opt}(T), определенного на интервале [0; T), а f(T;y) - средние издержки за время Т для плана y:[0; +\infty)  \to [0; +\infty).

Теорема 1.План y = y_{0} асимптотически оптимальный.

Таким образом, для достаточно больших горизонтов планирования Т планы y_{0}(t), 0 \le _t \le _T, все зубцы у которых имеют высоту Q_{0}, имеют издержки, приближающиеся к минимальным. Следовательно, эти планы Вильсона, являющиеся сужениями одной и той же функции y:[0; +\infty) \to [0; +\infty) на интервалы [0; T) при различных Т, можно использовать одновременно при всех достаточно больших Т.

Замечание.Согласно [11] решение проблемы горизонта планирования состоит в использовании асимптотически оптимальных планов, которые близки (по издержкам) к оптимальным планам сразу при всех достаточно больших Т.

Доказательство.По определению оптимального плана

\cfrac{f(T;y_{opt}(T))}{f(T;y)}\le 1

Найдем нижнюю границу для рассматриваемого отношения. При фиксированном Т можно указать неотрицательное целое число n такое, что

\cfrac{nQ_0}{\mu} \le T \le \cfrac{(n+1)Q_0}{\mu}

Так как Tf(T; y_{opt} (T)) и \cfrac{nQ_0}{\mu}f \left ( \cfrac{nQ_0}{\mu}; y_{opt}(T )\right ) - общие издержки на интервалах (0; Т) и \left (0; \cfrac{nQ_0}{\mu}\right ) соответственно при использовании оптимального на интервале времени (0; Т) плана, то, очевидно, поскольку второй интервал - часть первого (или совпадает с ним), первые издержки больше вторых, т. е.

Tf(T; y_{opt} (T)) \ge  \cfrac{nQ_0}{\mu}f \left ( \cfrac{nQ_0}{\mu}; y_{opt}(T )\right )

Далее, т. к. на интервале (0; nQ_{0}/ \mu ), включающем целое число периодов плана у _{0}, оптимальным является начальный отрезок этого плана у_{0}(nQ_{0} / \mu ), то

\cfrac{nQ_0}{\mu}f \left ( \cfrac{nQ_0}{\mu}; y_{opt}(T )\right ) \ge 
\cfrac{nQ_0}{\mu}f \left ( \cfrac{nQ_0}{\mu}; y_{0}(T )\right )

В правой части последнего неравенства стоит \cfrac{nQ _{0}}{\mu}\sqrt{2\mu gs} (здесь использована формула для минимального значения средних издержек f(T; y) при Т, кратном nQ_{0}/ \mu ). Из проведенных рассуждений вытекает, что

Tf(T; y_{opt}(T)) \ge \cfrac{nQ _{0}}{\mu}\sqrt{2\mu gs} ( 4.6)

Для общих издержек на интервалах (0; Т) и (0; (n + 1)Q_{0} / \mu  ) при использовании плана у _{0}, очевидно, справедливо следующее неравенство

Tf(T; y_{0}(T)) \le \cfrac{(n+1)Q_0}{\mu} f \left (\cfrac{(n+1)Q_0}{\mu}; y_0(T) \right )

Следовательно,

Tf(T; y_{0}(T)) \le \cfrac{(n+1)Q_0}{\mu} \sqrt{2\mu gs} ( 4.7)

Из неравенств (4.6) и (4.7) вытекает, что

\cfrac{ f(T;y_{opt}(T))}{f(T;y_0)} \ge 
\cfrac{n}{n+1} = 1 - \cfrac{1}{n+1} \ge 1 - \cfrac{Q_0}{\mu T}

Так как \cfrac{Q_0}{\mu T}\to 0 при Т  \to \infty, то, учитывая неравенство (4.5), из последнего неравенства выводим справедливость заключения теоремы 1. Таким образом, асимптотическая оптимальность плана у_{0} доказана.

При небольшом Т средние издержки в плане Вильсона могут существенно превышать средние издержки в оптимальном плане. Превышение вызвано скачками функции f (T; y_{0}(T)), связанными с переходами через моменты прихода очередных поставок (и увеличением общих издержек скачком на величину платы за доставку партии). Величину превышения средних издержек в плане Вильсона по сравнению с оптимальными планами можно рассчитать.

Пусть горизонт планирования T = t_{k} + \varepsilon, где t_{k} - момент прихода (k+1) -й поставки в плане Вильсона, \varepsilon \ge  0. Тогда, как можно доказать,

\lim_{\varepsilon \to \infty}{\cfrac{f(T;y_0(T))}{f(T; y_{opt}(T))}} = 
\lim_{\varepsilon \to \infty}{\cfrac{f(t_k + \varepsilon;y_0(t_k + \varepsilon))}{f(t_k + \varepsilon; y_{opt}( t_k + \varepsilon))}} = 1 + \cfrac{1}{2k}

Таким образом, затраты в плане Вильсона минимальные (относительно оптимального плана) при T = t_{k}, k = 1, 2, \dots , где t_{k} - моменты прихода поставок. Напомним, что план Вильсона является оптимальным при указанных Т. Однако при Т, бесконечно близком к t_{k}, но превосходящем t_{k}, затраты увеличиваются по сравнению с затратами в оптимальном плане в \{1+1 / (2k)\} раз. При дальнейшем возрастании Т отношение издержек (средних или общих) в плане Вильсона к аналогичным издержкам в оптимальном плане постепенно уменьшается, приближаясь к 1 при приближении (снизу) к моменту t_{k+1} прихода следующей поставки. А там - новый скачок, но уже на меньшую величину \{1+1 / (2k+2)\}. И т. д.

Сразу после прихода первой поставки отношение затрат составляет 1,5 (превышение на 50%), после прихода второй - 1,25 (превышение на 25%), третьей - 1,167 (превышение на 16,7%), четвертой - 1,125 (превышение на 12,5%), пятой - 1,1 (превышение на 10%), и т. д. Таким образом, при небольших горизонтах планирования Т превышение затрат может быть значительным, план Вильсона отнюдь не оптимальный. Но чем больше горизонт планирования, тем отклонение меньше. Уже после сотой поставки оно не превышает 0,5%.

Влияние отклонений от оптимального объема партии.В реальных производственных и управленческих ситуациях часто приходится принимать решения об использовании объемов партии, отличных от оптимальной величины Q_{0}, рассчитанной по формуле квадратного корня (4.2). Например, при ограниченной емкости склада или для обеспечения полной загрузки транспортных средств большой вместимости. Это возможно также в ситуации, когда величина партии измеряется в целых числах (штучный товар) или даже в десятках, дюжинах, упаковках, ящиках, контейнерах и т. д., а величина Q_{0} не удовлетворяет этому требованию и, следовательно, не может быть непосредственно использована в качестве объема поставки.

Поэтому необходимо уметь вычислять возрастание средних издержек при использовании напряженного плана с одинаковыми поставками объема Q, отличного от Q_{0}, по сравнению со средними издержками в оптимальном плане. Будем сравнивать средние издержки за целое число периодов. Как показано выше, они имеют вид

f_1(Q) = \cfrac{\mu g}{Q} +  \cfrac{sQ}{2}

где Q - объем партии. Тогда

\cfrac{f_1(Q) - f_1(Q_0) }{f_1(Q_0) } = 
\cfrac{1 }{2 } \left (\cfrac{Q-Q_0 }{Q }
\right )
\left (\cfrac{Q-Q_0 }{Q_0}
\right ) ( 4.8)

Это тождество нетрудно проверить с помощью простых алгебраических преобразований.

Пример 2.Пусть используется план с Q = 0,9 Q_{0}. Тогда

\cfrac{f_1(Q) - f_1(Q_0) }{f_1(Q_0) } = 
\cfrac{1 }{2 } \left (\cfrac{-0,1Q_0 }{0,9Q_0 }
\right ) \left (\cfrac{-0,1Q_0 }{Q_0 }
\right ) = \cfrac{0,01}{1,8 } = 0,0056

Таким образом, изменение объема партии на 10% привело к увеличению средних издержек лишь на 0,56%.

Пример 3.Пусть используемое значение объема поставки Q отличается от оптимального не более чем на 30%. На сколько могут возрасти издержки?

Из формулы (4.8) вытекает, что максимальное возрастание издержек будет в случае Q = 0,7 \cdot Q_{0}. Тогда

\cfrac{f_1(Q) - f_1(Q_0) }{f_1(Q_0) } = 
\cfrac{1 }{2 } \left (\cfrac{-0,3Q_0 }{0,7Q_0 }
\right ) \left (\cfrac{-0,3Q_0 }{Q_0 }
\right ) = \cfrac{0,09}{1,4 } = 0,0643

Таким образом, издержки могут возрасти самое большее на 6,43%.

Михаил Агапитов
Михаил Агапитов

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

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

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

Владимир Левадный
Владимир Левадный
Россия, г. Санкт-Петербург
Екатерина Куцак
Екатерина Куцак
Москва, МИОО, 2007