Распределение Q средств между N предприятиями.

Пусть х _n – средства, выделенные n-му предприятию; они приносят в конце года прибыль с_n(х _n).

Будем считать, что х _n принимает только целые значения, прибыль с_n(х _n) не зависит от вложения средств в другие предприятия и суммарная прибыль равна сумме прибылей, полученных от каждого предприятия. Тогда модель имеет вид:

Найти целочисленные неотрицательные переменные х _n (n=1,2,…,N), удовлетворяющие ограничению:

∑_n х _n = Q, (2.8.2)

и обращающие в максимум функцию

С=∑_n с_n(х _n). (2.8.3)

Здесь процесс распределения средств можно рассматривать как многошаговый, номер шага совпадает с номером предприятия; состояние будет определяться величиной s_n – количество средств, подлежащих распределению на n-м шаге (с конца).

Обозначим f_n(s_n) – условную оптимальную прибыль, полученную от последних n предприятий при распределении между ними s_n средств и вычисляемую в соответствие с динамическим рекуррентным соотношением:

f_n(s_n)=mах_"х(с_n(х _n) + f_n-1(s_n-1)), n=1,2,…,N. (2.8.4)

Пример 2.8.2. Пусть N = 4, Q =5, значения с_n(х _n) заданы в табл. 2.8.1.

Таблица 2.8.1.

Как и в предыдущем примере начинаем анализ с последнего предприятия. Индекс «1» соответствует последнему предприятию, а индекс «4» –первому. Для n=1 прибыль проставлена в последней колонке.

Для n=2

f₂(0)=mах[с₂(0)+f₁(0)]=0 при x ₂(0)=0,

f₂(1)=mах[с₂(1)+f₁(0),с₂(0)+f₁(1)]=mах[3+0,0+4]=4 при x ₂(1)=0,

f₂(2)=mах[с₂(2)+f₁(0),c₂(1)+f₁(1),с₂(0)+f₁(2)]=

=mах[4+0,3+4,0+6]=7 при x ₂(2)= 1,

f₂(3)=mах[с₂(3)+f₁(0),с₂(2)+f₁(1),с₂(1)+f₁(2),с₂(0)+f₁(3)]=

=mах[7+0,4+4,3+6,0+8]=9 при x ₂(3)=1,

f₂(4)=mах[с₂(4)+f₁(0),с₂(3)+f₁(1),с₂(2)+f₁(2),с₂(1)+f₁(3),с₂(0)+f₁(4)]=

=mах[11+0,7+4,4+6,3+8,0+13]=13 при х ₂(4)=0,

f₂(5)=mах[с₂(5)+f₁(0),с₂(4)+f₁(1),с₂(3)+f₁(2),с₂(2)+f₁(3),с₂(1)+f₁(4),с₂(0)+f₁(5)]

=mах[18+0,11+4,7+6,4+8,3+13,0+16]=18 при x₂(5)=5.

Для n=3

f₃(0)=mах[с₃(0)+f₂(0)]=0 при x₃(0)=0,

f₃(1)=mах[с₃(1)+f₂(0),с₃(0)+f₂(1)]=mах[6+0,0+4,]=6 при x ₃(1)=1,

f₃(2)=mах[с₃(2)+f₂(0),c₃(1)+f₂(1),с₃(0)+f₂(2)]=

=mах[9+0,6+4,0+7]=10 при x ₃(2)=1,

f₃(3)=mах[с₃(3)+f₂(0),с₃(2)+f₂(1),с₃(1)+f₂(2),с₃(0)+f₂(3)]=

=mах[11+0,9+4,6+7,0+9]=13 при x ₃(3)=1 или 2,

f₃(4)=mах[с₃(4)+f₂(0),с₃(3)+f₂(1),с₃(2)+f₂(2),с₃(1)+f₂(3),с₃(0)+f₂(4)]=

=mах[13+0,11+4,9+7,6+9,0+13]=16 при х ₃(4)=2,

f₃(5)=mах[с₃(5)+f₂(0),с₃(4)+f₂(1),с₃(3)+f₂(2),с₃(2)+f₂(3),с₃(1)+f₂(4),с₃(0)+f₂(5)]

=mах[15+0,13+4,11+7,9+9,6+13,0+18]=19 при x ₃(5)=1.

И, наконец, для n=4

f₄(0)=mах[с₄(0)+f₃(0)]=0 при x ₄(0)=0,

f₄(1)=mах[с₄(1)+f₃(0),с₄(0)+f₃(1)]=mах[8+0,0+6,]=8 при x ₄(1)=1,

f₄(2)=mах[с₄(2)+f₃(0),c₄(1)+f₃(1),с₄(0)+f₃(2)]=

=mах10+0,8+6,0+10]=14 при x ₄(2)=1,

f₄(3)=mах[с₄(3)+f₃(0),с₄(2)+f₃(1),с₄(1)+f₃(2),с₄(0)+f₃(3)]=

=mах[11+0,10+6,8+10,0+13]=18 при x ₄(3)=1,

f₄(4)=mах[с₄(4)+f₃(0),с₄(3)+f₃(1),с₄(2)+f₃(2),с₄(1)+f₃(3),с₄(0)+f₃(4)]=

=mах[12+0,11+6,10+10,8+13,0+16]=21 при х ₄(4)=1,

f₄(5)=mах[с₄(5)+f₃(0),с₄(4)+f₃(1),с₄(3)+f₃(2),с₄(2)+f₃(3),с₄(1)+f₃(4),с₄(0)+f₃(5)]

=mах[18+0,12+6,11+10,10+13,8+16,0+19]=24 при x ₄(5)=1.

Теперь соберем оптимальное решение (при последовательном рассмотрении всех состояний оптимальные переходы подчеркивались):

Для первого предприятия, когда s₄=5, видим, что x ₄(5)=1, значит, первое предприятие получает 1 и остается s₃=s₄– x ₄(5)=5–1=4. Находим лучшее размещение средств для второго предприятия (на третьем с конца шаге) при s₃=4. Это х ₃(4)=2, остается s₂=s₃– x ₃(4)=4–2=2. На втором (с конца) шаге x ₂(2)=1 и на последнее предприятие (первый с конца шаг) остается s₁= s₂ – x ₂(2)=2–1=1 и x₁(1)=1.

Максимум суммарной прибыли равен 24 у.е.

Динамическая задача управления запасами.

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

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

d_n – спрос на отрезке n от конца;

c_n(x,s) – затраты на отрезке n, связанные с выпуском х единиц изделия и с содержанием запасов, объем которых на конец отрезка равен s единиц. В этой системе обозначений подстрочный индекс «1» соответствует конечному, а «Т» – начальному состоянию.

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

Пусть f_n(s) – стоимость, отвечающая стратегии минимальных затрат на n оставшихся отрезках при уровне запасов s на начало n-го от конца отрезка;

x _n(s) – объем выпуска, обеспечивающий достижение f_n(s).

Пусть уровень запасов на конец планового периода равен нулю, тогда при уровне запасов s на начало последнего (1-го от конца) отрезка выпуск x ₁(s)=d₁–s и

f₁(s)= c₁(x,0)= c₁(d₁ – s,0), s=0,1,…,d₁.

Заметим, что если начальный уровень запасов отрезка n равен s, а объем выпуска – х, то величина (s+ х –d_n) – есть уровень запасов на конец данного отрезка, отсюда получаем общее рекуррентное соотношение в виде:

f_n(s) = min_x[c_n(x, s+ х –d_n)+ f_n-1(s+ х –d_n)], n=1,…,Т, s=0,1,…,d₁+…+ d_n.

Для упрощения вычислений предположим, что производственные мощности и складские площади ограничены, пусть х =0,1,…,5 и s=0,1,…,4. Допустим также, что спрос и затраты постоянны во времени, и пусть d_n=3, а c_n(x, s)= c(x)+hs, где первое слагаемое относится к производству, а второе определяется стоимостью содержания запасов (арендная плата за складские помещения, проценты за кредит для создания запасов, страховые взносы и собственно расходы по содержанию запасов). Пусть с(0)=0, с(1)=15, с(2)=17, с(3)=19,с(4)=21, с(5)=23; h=1.

Для n=1 f₁(0)=с(3)=19 при x ₁(0)=3,

f₁(1)=с(2)=17 при x ₁(1)=2,

f₁(2)=с(1)=15 при x ₁(2)=1,

f₁(3)=с(0)=0 при x ₁(3)=0.

Для n=2 f₂(0)=min[с(3)+0+f₁(0),c(4)+1+f₁(1),c(5)+2+f₁(2)]=

=min[19+19,21+1+17,23+2+15]=38 при x ₂(0)=3,

f₂(1)=min[с(2)+0+f₁(0),c(3)+1+f₁(1),c(4)+2+f₁(2),c(5)+3+f₁(3)]=

=min[17+19,19+1+17,21+2+15,23+3+0]=26 при x ₂(1)=5,

f₂(2)=min[с(1)+0+f₁(0),c(2)+1+f₁(1),c(3)+2+f₁(2),c(4)+3+f₁(3)]=

=min[15+19,17+1+17,19+2+15,21+3+0]=24 при x ₂(2)=4,

f₂(3)=min[с(0)+0+f₁(0),c(1)+1+f₁(1),c(2)+2+f₁(2),c(3)+3+f₁(3)]=

=min[0+19,15+1+17,17+2+15,19+3+0]=19 при x ₂(3)=0,

f₂(4)=min[с(0)+1+f₁(1),c(1)+2+f₁(2),c(2)+3+f₁(3)]=

=min[0+1+17,15+2+15,17+3+0]=18 при x ₂(4)=0.

Для n=3 f₃(0)=min[с(3)+0+f₂(0),c(4)+1+f₂(1),c(5)+2+f₂(2)]=

=min[19+38,21+1+26,23+2+24]=48 при x ₃(0)=4,

f₃(1)=min[с(2)+0+f₂(0),c(3)+1+f₂(1),c(4)+2+f₂(2),c(5)+3+f₂(3)]=

=min[17+38,19+1+26,21+2+24,23+3+19]=45 при x ₃(1)=5,

f₃(2)=min[с(1)+f₂(0),c(2)+1+f₂(1),c(3)+2+f₂(2),c(4)+3+f₂(3),c(5)+4+f₂(4)]=

=min[15+38,18+26,21+24,24+19,23+4+18]=43 при x ₃(2)=4,

f₃(3)=min[с(0)+0+f₂(0),c(1)+1+f₂(1),c(2)+2+f₂(2),c(3)+3+f₂(3),c(4)+4+f₂(4)]=

=min[0+38,16+26,19+24,22+19,25+18]=38 при x ₃(3)=0,

f₃(4)=min[с(0)+1+f₂(1),c(1)+2+f₂(2),c(2)+3+f₂(3),c(3)+4+ f₂(4)]=

=min[1+26,17+24,20+19,23+18]=27 при x ₃(4)=0.

И, наконец, для n=4

f₄(0)=min[с(3)+0+f₃(0),c(4)+1+f₃(1),c(5)+2+f₃(2)]=

=min[19+48,21+1+45,23+2+43]=67 при x ₄(0)=3 или 4,

f₄(1)=min[с(2)+0+f₃(0),c(3)+1+f₃(1),c(4)+2+f₃(2),c(5)+3+f₃(3)]=

=min[17+48,19+1+46,21+2+43,23+3+38]=64 при x ₄(1)=5,

f₄(2)=min[с(1)+f₃(0),c(2)+1+f₃(1),c(3)+2+f₃(2),c(4)+3+f₃(3),c(5)+4+f₃(4)]=

=min[15+48,18+45,21+43,24+38,23+4+27]=54 при x ₄(2)=5,

f₄(3)=min[с(0)+0+f₃(0),c(1)+1+f₃(1),c(2)+2+f₃(2),c(3)+3+f₃(3),c(4)+4+f₃(4)]=

=min[0+48,16+45,19+43,22+38,25+27]=48 при x ₄(3)=0,

f₄(4)=min[с(0)+1+f₃(1),c(1)+2+f₃(2),c(2)+3+f₃(3),c(3)+4+ f₃(4)]=

=min[1+45,17+43,20+38,23+27]=46 при x ₄(4)=0.

Сведем результаты вычислений в таблицу 2.8.2.

Таблица 2.8.2.

Читателю предоставляем возможность проверить результаты вычислений для n = 5 ¸ 8. Теперь, задаваясь различным уровнем запасов на начало планового периода, можно определить оптимальные стратегии для любого Т от 1 до 8. Так, например, если исходный уровень запасов на начало планового периода равен нулю, то оптимальный календарный план при Т=4 будет:

x ₄(0)=3, x ₃(0)=4, x ₂(1)=5, x ₁(3)=0 или

x ₄(0)=4, x ₃(1)=5, x ₂(3)=0, x ₁(0)=3

и минимальная общая сумма затрат составит 67.

Пусть Т=8, тогда минимальная общая сумма затрат составит 127 и x ₈(0)=5, останется на следующий отрезок 5 - 3=2, имеем x ₇(2)=5, значит на следующий отрезок останется 2+5 - 3=4 и x ₆(4)=0, останется 4 - 3=1 и x ₅(1)=5, останется 1+5 - 3=3, x ₄(3)=0, останется 3 - 3=0, x ₃(0)=4, останется 4 - 3=1, x ₂(1)=5, останется 1+5 - 3=3 и x ₁(3)=0.

В таблицу 2.8.3 сведем оптимальные стратегии (планы выпуска) для плановых периодов длительностью Т и начальным уровнем запасов равным нулю.

Обратим внимание на Т=4. Здесь два оптимальных решения, дающих минимальные затраты 67. При Т=7 имеем три решения с одинаковым значением целевой функции 115.

Таблица. 2.8.3

Анализ оптимальных вариантов производственной программы, приведенных в табл. 2.8.3, свидетельствует о том, что январский выпуск поначалу возрастает с ростом Т, а затем колеблется около 5, аналогично, среднемесячные затраты испытывают колебания (при минимальном их значении 15.8). Проводя подобный анализ для следующих месяцев (при достаточно большом Т), получим оптимальное решение для бесконечного периода планирования. В условиях настоящего примера таким решением будет повторение производственного цикла (5,5,0,5,0,…) при среднемесячных затратах 15.8. (вычисления опускаем).

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

1) Управляемые переменные и соответствующие ограничения группируются по шагам, и многошаговый процесс принятия решений исследуется в определенной последовательности.

2) Для выбора оптимальных значений переменных на рассматриваемом шаге используется значение состояния только на данном шаге.

3) Решение, принимаемое при заданном текущем состоянии системы, оказывает прогнозируемое влияние на состояние системы на последующем шаге.

4) Оптимальность текущего решения оценивается в терминах прогнозируемого экономического эффекта для рассматриваемого шага и всех последующих шагов.

Что делать, если нет взаимности? А теперь спустимся с небес на землю. Приземлились? Продолжаем разговор...

Что вызывает тренды на фондовых и товарных рынках Объяснение теории грузового поезда Первые 17 лет моих рыночных исследований сводились к попыткам вычис­лить, когда этот...

Что будет с Землей, если ось ее сместится на 6666 км? Что будет с Землей? - задался я вопросом...

ЧТО И КАК ПИСАЛИ О МОДЕ В ЖУРНАЛАХ НАЧАЛА XX ВЕКА Первый номер журнала «Аполлон» за 1909 г. начинался, по сути, с программного заявления редакции журнала...

Не нашли то, что искали? Воспользуйтесь поиском гугл на сайте: