«Нижний предел»= 0,0943 t 2 |
3,9529 t |
72,886 . |
«Верхний предел»= 0,2289 t 2 |
4,7651 t |
75,337. |
Будем исследовать процесс производства одного вида продукции с целью максимизации прибыли от его реализации. Перепроизводство приводит к очевидным убыткам, так как резко начинают увеличиваться переменные затраты.
Рассмотрим процесс производства продукции в дискретные моменты времени t = 0,1,...,T, где Т – плановый период. Спрос на продукцию в эти моменты времени определяется, по предположению, заданной функцией
Qd t
5,2338
t 72,07 .
При эффективном планировании производства предложение должно стремиться полностью удовлетворить спрос: Qs Qd , поэтому необходимо, чтобы на каждом этапе t спрос был полностью удовлетворён, причём это может достигаться необязательно через выпуск продукции, равный спросу на каждом этапе t, но и через накапливаемый остаток нереализованной продукции, выпущенной в предыдущие периоды. Поэтому важным моментом в решении задачи будет не равенство выпуска спросу в каждом месяце, а равенство суммарного годового спроса и суммарного годового выпуска продукции. Функция переменных затрат VC(Qs) на единицу продукции зависит от объёма выпуска и имеет зависимость второй степени VC(Qs)
0,107Q 26,846. Оптовая цена единицы продукции Р составляет 55 руб. и остаётся постоянной в течение всего планируемого периода. Постоянные затраты на выпуск этого вида продукции FC составляют 50 тысяч руб. в месяц. Известно, что налог на прибыль τ составляет 24%. С учётом всего этого можно определить функцию прибыли после налоговых отчислений так:
f2
f1 ( (Qs))
(1
) , f1 ( (Qs)) FC
P VC Qs
Qs t
FC .
Для производителя продукции наиболее предпочтительным является
уровень |
постоянной интенсивности |
производства, |
то есть когда |
|
Qs t const , или, другими словами, u t |
Qs t |
1 Qs t |
0 . И в случае уве- |
|
личения |
выпуска продукции (u(t)>0), и |
в |
случае уменьшения выпуска |
|
(u(t)<0) функция прибыли изменяется на определённую величину, которую можно определить следующим образом:
(P Vc(Qs u))
(Qs u) (P Vc(Qs))
Qs .
21
В этом случае величина изменения прибыли будет иметь следующий функциональный вид:
f3 (u) |
f |
2 |
( |
(Qs u)) |
f |
2 |
( |
(Qs)) |
28,154 u 0,0048 Qs 2 u 0,0048 Qs u 2 |
|
|
|
|
|
|
|
|
||
0,0016 |
|
u3 |
0,214 Qs |
|
u |
|
0,107 |
u 2 . |
|
Задача планирования выпуска продукции в общем виде формулируется в данном случае следующим образом: найти функцию объёма производства продукции Qs(t), t=1,2,...,T и динамику необходимого изменения этого объёма, выражаемую функцией u(t), t=1,2,...,T-1, чтобы свести к максимуму суммарную прибыль от реализации данной продукции в течение планового периода Т. Математическая модель задачи в нашем случае примет вид:
|
T 1 |
|
|
|
|
J |
f1 (Qs(t)) |
f3 (u(t)) (1 |
) f1 (Qs(T )) (1 |
) min. |
(2.1) |
|
t 1 |
|
|
|
|
Исходное количество продукции задаёт начальное условие |
|
||||
|
Qs(1) |
Qs1 , или |
76,7446 Qs(1) |
80,331. |
(2.2) |
Уравнение процесса записывается в канонической форме:
Qs(t+1)=Qs(t)+u(t). |
(2.3) |
||
0,0943 t 2 3,9529 t |
72,886 |
Qs(t) 0,2289 t 2 4,7651 t 75,337 . |
(2.4) |
T |
|
T |
|
|
Qst |
Qd t . |
(2.5) |
t |
1 |
t 1 |
|
В соответствии с общей постановкой задачи для многошаговых процессов соотношение (2.4) определяет в качестве состояния системы Qs(t), в качестве управления используется u(t). Кроме того, Qs(t)≥0 – объём выпуска продукции не может быть отрицательным. Однако, для того чтобы решать задачу (2.1) – (2.5) методом Лагранжа, необходимо освободиться от ограничения на состояние Qs(t).
Имея в виду такую возможность ликвидации отдельных ограничений с включением в функционал «штрафа» за их нарушение, тем не менее функционал (2.1) предпочтительнее оставить без изменения. Наложение таких строгих ограничений на объём выпуска может привести к невозможности решения данной задачи методом Лагранжа, поэтому вместо явного учёта ограничений Qs(t)≥0 проще воспользоваться численным методом прямой прогонки при решении краевой задачи. При этом при каждой итерации бу-
22
дет учитываться только обобщённое ограничение о неотрицательности объёма выпуска, если окажется, что Qs(t)<0 при первом же t, будут прекращаться вычисления и переходить к новой итерации. Таким образом, можно ограничиться более простой постановкой задачи в отличие от использования метода «штрафных функций», а неучтённые ограничения на состояние Qs(t)≥0 будут учитываться непосредственно в процессе вычислений. Также для удобства решения поставленной задачи не будут учитываться налоговые отчисления и постоянные затраты, так как величина налоговой ставки и постоянных затрат постоянны, исключение их из функционала не изменит найденный в процессе решения задачи оптимальный план управления. Просто на конечном этапе при представлении найденного решения прибыль будет указываться с учётом постоянных затрат и налоговых отчислений.
Решение задачи
Итак, известно, что планируемый расчётный период исчисляется 12 месяцами Т=12, спрос на каждом расчётном периоде вычисляется по формуле Qd(t)=72,07+5,2338*t. Исходные данные отражены в таблице 2.2.
Таблица 2.2 – Данные по спросу продукции на каждом этапе планирования
Этап t |
1 |
2 |
3 |
4 |
5 |
6 |
(мес) |
|
|
|
|
|
|
Спрос |
77,3038 |
82,5376 |
87,7714 |
93,0052 |
99,239 |
103,4728 |
(тыс. бут.) |
|
|
|
|
|
|
7 |
8 |
9 |
10 |
11 |
12 |
|
108,7066 |
113,9404 |
119,1792 |
124,408 |
129,6418 |
134,8756 |
1273,076 |
Итерация 1. T=12 |
|
|
|
|
|
||
Чтобы задать необходимым образом Qs*(T ) |
или Qs*(12) , можно решить на этапе |
||||||
t=12 |
задачу |
максимизации |
прибыли |
при |
ограничениях: |
||
0,0943 |
12 2 |
3,9529 |
12 72,886 |
Qs (12) |
0,2289 |
12 2 |
4,7651 12 75,337 . |
Таким образом, необходимо решить задачу:
|
f1 |
Qs 12 |
FC |
P VC Qs 12 Qs 12 FC |
max |
при 106,74 ≤Qs |
12 ≤ 165,48. Или же |
|
|||
f |
1 |
Qs 12 |
FC |
28,154 Qs 12 0,0016 Qs 12 3 |
0,107 Qs 12 2 50 max , |
|
|
|
|
|
|
|
|
|
|
23 |
|
при 106,74 ≤Qs
12 ≤ 165,48.
В результате решения поставленной задачи Qs*(12) =106,74 (тыс.бут.). В
этом случае по формуле (3.2.18) можно получить y(12)= 3,69.
t=11. Теперь, подставляя значение y(t+1)= 3,69, Qs(12)=106,74 и представляя u как разность Qs(t+1)-Qs(t), неизвестным в этом неравенстве остаётся только Qs(t). Таким об-
разом, на этапе t=11 решается система неравенств:
3,69 |
28,154 |
0,0048 |
Qs(11)2 |
0,0096 Qs(11) (106,74 |
Qs(11)) |
||||||
0,0048 |
(106,74 |
Qs(11))2 |
0,214 |
Qs(11) |
0,214 |
(106,74 |
Qs(11)); |
||||
_ при _ Qs(11) |
104,96 |
|
|
|
|
|
|
|
|||
3,69 |
28,154 |
0,0048 |
Qs(11)2 |
0,0096 Qs(11) (106,74 |
Qs(11)) |
||||||
0,0048 |
(106,74 |
Qs(11))2 |
0,214 |
Qs(11) |
0,214 |
(106,74 |
Qs(11)); |
||||
_ при _ 104.96 |
Qs(11) |
155.45 |
|
|
|
|
|
||||
3,69 |
28,154 |
0,0048 |
Qs(11)2 |
0,0096 Qs(11) (106,74 |
Qs(11)) |
||||||
0,0048 |
(106,74 |
Qs(11))2 |
0,214 |
Qs(11) |
0,214 |
(106,74 |
Qs(11)); |
||||
_ при _ Qs(11) |
155.45. |
|
|
|
|
|
|
|
|||
Решение данного неравенства даёт результат Qs(11)= 106,74 и u =0, при том, что правая часть данного неравенства равна левой. По имеющимся ре-
зультатам можно определить y(11):
y(11) |
3,69 |
28,154 |
0,0048 106,742 |
0,214 |
106,74 |
||||
0,0096 106,74 |
0 |
0,0048 |
02 |
0,214 |
0 |
0. |
|
||
На следующем шаге в формулу вместо значения y(t+1), подставляя зна-
чение y(11), решается новая система неравенств
t=10. Проделав аналогичную шагу 11 операцию, получим объём выпус-
ка (Qs(10)) в данной период равный 102,985 тыс.бут., а u=3,755, и можно получить y(10), равный:
y(10) |
0 |
28,154 |
0,0048 102,9852 |
0,214 102,985 |
|
||||
0,0096 102,985 |
3,755 |
0,0048 |
3,7552 |
0,214 |
3,755 |
3,69. |
|||
|
|
|
|
24 |
|
|
|
|
|
Теперь имеются необходимые данные для решения системы неравенств
на следующем этапе t=9.
t=9
Qs(9)=100,824 тыс.бут. u = 2,161
y(8) = -3,472
Проделывая аналогичные операции на последующих этапах можно получить следующие результаты:
t=8 |
t=7 |
Qs(8)= 98,474 тыс.бут. |
Qs(7)= 95,936 тыс.бут. |
u = 2,35 |
u = 2,538 |
y(7) = -0,79 |
y(6) = 3,716 |
t=6 |
t=5 |
Qs(6)= 112,168 тыс.бут. |
Qs(5)= 90,293 тыс.бут. |
u = -16,232 |
u = 21,8750 |
y(5) = -4,518 |
y(4) = 3,825 |
t=4 |
t=3 |
Qs(4)= 98,06 тыс.бут. |
Qs(3)= 91,692 тыс.бут. |
u = -7,7670 |
u = 6,368 |
y(3) = 6,808 |
y(2) = 14,229 |
t=2 |
t=1 |
Qs(2)= 85,783 тыс.бут. |
Qs(1)= 80,331 тыс.бут. |
u = 5,909 |
u = 5,452 |
y(1) = 25,418 |
|
Суммарный годовой объём выпуска по приближенно оптимальному плану равен 1 170,026 тыс. бут., однако он не полностью удовлетворяет годовой спрос. Необходимо провести ещё пару итераций, чтобы попытаться повысить годовой объём выпуска.
Итерация 2. На данной итерации можно попробовать принять значение выпуска в период t = 12 равному верхнему значению доверительного интервала для данного периода, то есть 165,48 тыс. бут., тогда результат решения примет другой вид. Суммарный годовой объём выпуска равен 1
25