Динамические модели управления производственно-финансовой деятельностью предприятия Текст научной статьи по специальности «Экономика и бизнес»

26 (281) - 2012

Экономико-математическое

моделирование

УДК 330.322.5

ДИНАМИЧЕСКИЕ МОДЕЛИ УПРАВЛЕНИЯ ПРОИЗВОДСТВЕННО-ФИНАНСОВОЙ ДЕЯТЕЛЬНОСТЬЮ ПРЕДПРИЯТИЯ

А. В. МИЩЕНКО,

доктор экономических наук, профессор кафедры логистики E-mail: nesterovich@gnext. ru Национальный исследовательский университет «Высшая школа экономики»

О. А. АРТЕМЕНКО,

аспирант кафедры математических методов в экономике

E-mail: o. a. artemenko@gmail. com Российский экономический университет имени Г. В. Плеханова

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

Ключевые слова: оптимизационная модель, конвейерная система, технологическая операция, заемные средства.

Оптимизационные задачи управления материальными и производственными ресурсами в конвейерных системах

Рассмотрим так называемые конвейерные системы обработки заявок в производственных

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

необходимы производственные ресурсы, к которым относятся оборудование, транспортные средства, компьютеры, станки, приборы.

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

Детерминированная модель управления производственными ресурсами. В простейшем случае конвейерная система обработки может быть графически задана ориентированным циклическим графом, состоящим из нескольких параллельных цепочек, задающих перечень операций по каждому виду заявок и последовательность обработки на этих операциях. Такой ориентированный граф еще называют П-сетью. На входе П-сети подается поток материальных ресурсов производства (или финансовый поток, предназначенный для закупки материальных ресурсов), который должен быть обработан с учетом имеющегося производственного арсенала на всех технологических операциях, и на выходе таким образом будет получена готовая продукция. Графически такой процесс обработки для М видов заявок представлен на рис. 1.

Рассмотрим задачу оптимизации управления потоками материальных и производственных ресурсов данной модели по критерию максимизации суммарного маржинального дохода от произведенной продукции в периоде (0,7), который будет получен после ее реализации. С учетом сделанных ранее предположений необходимо в этом случае решить следующую задачу оптимального управления:

M T

ZP. J 4щ(*) dt ^ max, (1)

1=1 о

где Р.. - маржинальный доход от выпуска одной единицы продукции вида 1;

Ч N (t) - искомая интенсивность обработки заявок вида i на последней N. операции (1 = 1, 2,..., M). 1

Пройдя обработку на операции 0т , на выходе будет получена конечная продукция вида i (/' = 1, 2,..., М). Учитывая, что ни на одной операции производственного цикла не может быть обработано большее количество заявок, чем поступило, необходимо учитывать следующее ограничение:

J q j (t ')dt' < V (0) + J q ^ _i(t ')dt', (2)

где = 1,

j = 1,

2,., М; 2,., N t е (0,7);

V.. (0) - объем незавершенного производства на операции 0.. в момент t = 0. Далее будем предполагать, что обработка заявок на каждой операции происходит с использованием производственных ресурсов, к которым относятся станки, оборудование, технологическая оснастка, обслуживающий персонал и др. Объем этого вида ресурсов задан с использованием вектора С = (С1,...,Ст). Чтобы на операции 0.. обеспечить единичную производительность, необходимы производственные ресурсы в количестве, заданном векторами а.. = (а1., а2,..., ат),

V V V V

.= 1, 2,., М;. = 1, 2,., N..

Если же необходимо обеспечить производительность qj. на операции 0.., то требуемый объем производственных ресурсов задается как qj. х а.. = ^ ха1,qij ха2,...,qij хат). При задании произво-дительностей qij (0 на каждой операции qij (. = 1, 2,., М;. = 1, 2,., N объем используемых ресурсов не должен превышать соответствующих компонент вектора С = (С1,...,Ст).Иными словами, должно выполняться следующее неравенство:

м N

ZZ4j(t)aj < C, l = 1, 2,., M; Vt e (0,T). (3)

=1 j =1

Рис. 1. Технологическая последовательность операций П-сети

Количество выпущенной на интервале (0,7) продукции каждого вида не должно быть выше объема спроса РХ , что определяется следующим

неравенством:

Т

| Ящ(Х Э < РХ,, г = 1, 2,..., М.

(4)

Рассмотрим алгоритм решения задачи (1) - (3) без ограничения на спрос (4) при условии, что V.. (0) > 0, 1 = 1, 2,..., М;. = 1, 2,..., N.

В момент Х =0 в этом случае необходимо назначить производительности дN (1 = 1, 2,., М), обеспечивающие максимальный маржинальный

доход за единицу времени, т. е.

Ер гЧ щ

^ тах

при ограничениях

М

ЕЯ^аЩ < С, I = 1, 2,., М,

(5)

(6)

дщ> 0, г = 1, 2,., М. (7)

Задача (5) - (7) является задачей линейного программирования, решение которой не представляет принципиальных трудностей. Получив решение задачи (5) - (7), которое обозначим = (#1ла, Яш2,. .,я\Щи), можно вычислить минимальный момент времени т1, в который закончатся заявки на одной из операций 01щ1,...,0МЩ .

Очевидно, что

(0)

т 1 = тт 1

г =1, 2,., М.

(7.1)

Ч

дкщп < дк1щК1_1, (8)

Чщ > ^ Чк1щК1_1 > 0, г = 1 = 1, 2,.,М. Ограничения (7.2) и (8) свидетельствуют о том, что производственные ресурсы для выпуска продукции вида Кхдолжны выделяться не только на последнюю операцию КК , но и на продукцию Кк -1. 1

Далее находится минимум следующих соотношений:

тт

г=1,М

(ТттК1 )

2

Ч гК

2

Ч Кк

22

Здесь и ЧКк11 - решение задачи (5), (7.2), (8). Таким образом, если V.. > 0, то решение задачи (1) - (3) (без ограничения на спрос) сводится к решению серии задач линейного программирования.

Если есть ограничения на спрос, то необходимо при решении задачи (1) - (4) отслеживать момент

тг

времени тг, для которого |Ящ (Х= Ьг, г = 1, 2,., М. ' 0 '

Далее с момента времени т выпуск продукции вида г останавливается.

Рассмотрим ситуацию, когда есть еще ограни-

, т. е. | Ящ^Ш, > Zl, г = 1, 2,., М.

Если т1 > Т, то задача (1) - (3) решена. При этом оно будет оптимальным не только для целевой функции (1), но и для целевой функции вида

М Х

Ев г]Як(Х') Л', УХ е (0,Т).

г=1 0

Таким образом, предложенный подход дает решение задачи (1) - (3) для любого временного интервала (0,Х) с (0,Т).

Пусть выполняется тх < Т и минимум в выражении (7.1) достигается на каком-либо 1 < К < М . В этом случае для организации выпуска продукции вида К необходимо выделять производственные ресурсы не только на операцию , но и на

предыдущую операцию ^. Следовательно, производительности на конечных операциях, как и ранее, назначаются исходя из того, что максимизируется целевая функция (5), а вместо ограничений (6) - (7) используются следующие ограничения:

М

ЕЯща'щ + Як1Кк1-1 аЩк1_1 < С, 1= 1, 2,.,М, (7.2)

чение на заказы

Здесь X . - объем продукции вида г, который необходимо поставить заказчику на интервале времени (0,Т). В этом случае необходимо оценить объем незавершенного производства по всем видам продукции, т. е. выполнить неравенство

N

ЕV. > Zl , г = 1, 2,., М. (9)

.=1

Если ограничение (9) не выполняется, то при распределении финансового потока F(Х), направляемого на закупку материальных ресурсов производства для всех видов продукции ^(Х) = F1(t) + F2(Х) +... + FM (Х), необходимо обеспечить закупку этих ресурсов в необходимых объемах для выполнения ограничения (9). Обозначим через у стоимость материальных ресурсов всех видов для выпуска одной единицы продукции вида и введем еще один показатель иы (Х) = Fi (Х) / у г, г = 1, 2,., М. Таким образом иы (Х) - это интенсивность потока материальных ресурсов, подаваемых на начало технологической цепи выпуска продукции вида (г = 1, 2,., М). В этом случае объем материальных ресурсов, поставляемых к моменту Х' е (0,Т) для вы-

К1-1

г=1

1=1

пуска продукции вида i рассчитывается по формуле

К,«') = ]Ч 0« №.

0

В данной ситуации необходимым условием решения задачи (1) - (4), (9) является существование

м

таких ит (г), где ит (í) = Fl (Г)/ у, и ДО = X Fl (г),

1=1

что

N т

XV, +1 и, оС Э > 2, V/ = 1, 2,..., N. (10)

./=1 о

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

В то же время выполнение условия (10) еще не гарантирует выполнения заказа в силу того, что, несмотря на то, что все материальные ресурсы поставлены в необходимых объемах, производственные мощности предприятия не могут обеспечить выпуск продукции каждого вида в объемах X. за период (0,7).

Рассмотрим ситуацию, когда

2 * V (П)

Выполнение ограничения (11), в частности, означает, что объем незавершенного производства на последней операции по каждому виду продукции достаточен для выполнения заказа. В этом случае решение задачи оптимального управления (1) - (4), (9) сведется к решению следующей задачи линейной оптимизации: м

ЕР^ж ^таХ (12)

=1

м

X С, 1 = 1, 2,., м, (13)

1

т х дт> X,, 1 = 1, 2,., м, (14)

Т х Р^, 1 = 1, 2,.,М, (15)

Чж> 0, 1 = 1, 2,.,М. (16)

Здесь Т - продолжительность директивного периода (0,7). Если задача (12) - (16) имеет решение, то производительности являются постоянными в течение всего интервала времени.

В том случае, если задача (12) - (16) не имеет решения, необходимо инвестировать средства в увеличение производственных мощностей, т. е. в приобретение дополнительного оборудования, станков, привлечение дополнительного количества обслуживающего персонала. Таким образом, необходимо перейти от производительных ресур-

сов, заданных вектором С = (С1,...,Ст), к увеличенным производственным ресурсам С + у = (С1 + у1,..., Ст + ут), где у, - дополнительное количество производственных ресурсов вида/(/=1, 2,., м). Естественным требованием при приобретении дополнительного производственного оборудования является минимизация затрат на его приобретение. Это требование можно выполнить, решив следующую оптимизационную задачу:

X у,у/ ^ т1П,

м

Xщ <С + у, 1 = 1, 2,.,М,

1

Т х дщ > Zi, 1 = 1, 2,., М, Т х Яж, < , 1 = 1, 2,.,М,

> 0, 1 = 1, 2,.,М, у > 0, у е 2+.

Решив эту задачу и получив в качестве решения объемы дополнительно закупаемого оборудования у,, переходим к решению задачи (12) - (16) с учетом того, что в правой части производственные ресурсы в количестве С1 +у1(1 = 1, 2,.,М).

Если ограничение (11) для какого-либо вида продукции Р не выполняется, то существует такое /, что

V// > , 1 < р < М, (17)

тогда задача оптимального управления (1) - (4), (9) сводится к решению следующей задачи линейного программирования:

М

Xе гЧж ^ тах (18)

=1

М Ыр

X Ц жаж + XЦрКаРк < С, 1 =1, 2,.,M, (19)

1=1 ¿=1

Ч р/ = Ч р/+1 =... = Ч ржр, (20)

ТхЧж >X, 1 = 1, 2,.,М, (21)

Т х Цж < Р(,, 1 = 1, 2,.,М, (22)

Чж > 0, Яр/ > 0, 1 = 1, 2,.,М,/ = 1, 2,., Ыр. (23)

Формулируя задачу оптимизации управления производственными ресурсами (17) - (23), необходимо учитывать, что по виду продукции Р следует осуществить обработку незавершенного производства не только на операциях ОрМр, но и на предшествующих операциях Оржр1,..., Ор1, что требует дополнительного привлечения производственных ресурсов - ограничение (19), а также согласования производительностей по операциям продукции вида Р - ограничение (20).

Оптимизация управления производственн ыми ресурсами с учетом риска. Формируя производственный план выпуска конечной продукции и загрузки производственного оборудования (задача (1) - (4), (9)), лицо, принимающее решение, не всегда может оценить будущий маржинальный доход Рг(г' = 1, 2,., М) по каждому виду выпускаемой продукции.

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

р; - Pi

-> Pi

^ РЧ - Рд

В этом случае при решении задачи (1) - (4), (9) естественно использовать в целевой функции (1) вместо вг математическое ожидание маржинального дохода вг-

в = Е в. Р.. (24)

з=1

Обозначим объем выпуска по каждому виду выпускаемой продукции в решении задачи (1) - (4), (9) через Ql.

Очевидно, что

Q, = JЯт(t)dt, i = 1, 2,...,M.

(25)

Соответственно объем переменных издержек, связанный с выпуском одной единицы продукции вида г через Ь,, (г = 1, 2,., п). Тогда объем переменных издержек продукции вида г при производственной программе Q = (0^..., Qm) определяется как = Ь,0,, г = 1, 2,.,М.

Обозначим долю затрат на выпуск продукции вида г при производственной программе 0 = (01,.., 0т) через Аг. Очевидно, что

Л,- =-

b,Q,

J b q т (t )dt

m t

ЕЬз0з Е\Ь,Ят (Х^Х

3 =1 г=1 0

В этом случае риск производственной программы 0 задачи (1) - (4), (9) может быть оценен следующим образом:

R = Е^2 Л2 + 2^1 covj Л ;Лj, (27)

i=1 i=1 i> j

2

где a, - дисперсия маржинального дохода по продукции вида i;

coVj - ковариация маржинального дохода продукции i и продукции j.

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

M _T

ЕР. J Ям, (t) dt ^ max, (28)

i=1 о

¡Я, (t ')dt' < Vj (0) + J q.j-i (t ')dt', (29) 00 Vt e (0,T) i = 1, 2,..., M; j = 1, 2,..., N,

M M

Ес2Л2 + 2ЕЕcoVj Л,Лj <Rq. (30) i=1 i=1 j >i Здесь Rq - допустимый уровень риска доходности производственной программы Q = (Q1,..., Qm).

M Np

Е Е ЯР (t)a'j < C, (31)

;=1 k=1

Vt e (0,T), l = 1, 2,..., m,

i

Jqm(t)dt < Ptt, i = 1, 2,...,M,

J qm(t )dt, > Zt, i = 1, 2,., М.

(32)

(33)

В модели (28) - (33) оптимизируется маржинальный доход от реализации продукции, выпущенной на директивном периоде (0,7), заданной соотношением (28). Напомним, что вг - это маржинальный доход от реализации одной единицы продукции вида г' .

Формула (29) задает балансовые ограничения на объем обработки незавершенного производства на каждой операции О.., когда технологическая схема выпуска продукции задается П-сетью (г = 1, 2,., М; 3 = 1, 2,., Кг).

Неравенство (30) задает ограничение на допустимый уровень риска доходности производственной программы, где а, определяются по формуле (26), (26) 52 - дисперсия маржинального дохода по продук-

ции вида г, ооу. - ковариация доходности продукции вида г и продукции вида. Неравенство (31) ограничивает производительность на каждой операции О.., что связано с существующими производственными мощностями предприятия. Соотношения (32) и (33) -

это соответственно ограничения на объем спроса и объем заказа на выпускаемую продукцию.

Анализ устойчивости в модели оптимального управления производственными ресурсами предприятия. Вернемся к модели (1) - (4), (9) и рассмотрим ситуацию линейного роста маржинальной доходности р, от инфляции, т. е. будем предполагать, что при уровне накопленной инфляции маржинальный доход Р; (Е) = Р; + п Р; Е, где Е - накопленная инфляция в долях. Будем также считать, что для Vдопустимого решения задачи (1) - (4), (9) Q = (0,...,Qm), компоненты вектора Q - целочисленны. Это соответствует реалиям практически любого промышленного предприятия, выпускающего конечную продукцию для потребителя. Обозначим множество всех допустимых производственных программ модели (1) - (4), (9) через 0 = (01,..., QL).

Учитывая целочисленность векторов 0/(/ = 1, 2,., L), число допустимых производственных программ будет конечным. Обозначим Р (Е) значение целой функции (1) на производственной программе 0 при уровне накопленной инфляции Е. Очевидно, что в этом случае Р (Е) вычисляется по следующей формуле:

М т

р/ © = £ (Р1 + п1р1 Е){ чж (г ж / = 1, 2,., М.

1=1 0 Функция Р (Е) является линейной относитель-

Е dF/ (Е) М но Е и-= > п

Ё n,ß, i qN (t) dt, j = 1, 2,..., M.

=1 0

Предположим, что при уровне накопленной инфляции Е = 0 оптимальной является производственная программа 01 = (011,...,). Будем считать также, что программы множества 0 упорядочены

dFi (Е)

по возрастанию производственных -— / = 1,

2,., Ь. ^

Таким образом, еслиp > q, то

dFp (Е) > dFj (Е)

для любых 1 < p < L и 1 < q < L и p Ф q.

В этом случае, если l > L, то рост маржинального дохода на производственной программе Qk(k > L) по мере накопления инфляции. Следовательно, каждое из уравнений вида F1 (Е) = Fk (Е), k = 1, 2,..., L будет иметь одно положительное решение. Обозначим эти решения как Ек и выберем min Ек = Ер (p > l).

утрачивает оптимальность и оптимальной становится производственная программа 0р (р > 1). Повторяя предыдущую цепочку рассуждений для множества производственных программ 0р, 0р+1,..., 0L, найдем следующую производственную программу 0 и соответствующее значение накопленной инфляции Е2 (/ > р), при котором оптимальной становится производственная программа . Очевидно, что через конечное число итерации и при некотором уровне накопленной инфляции оптимальной станет производственная программа 0 ^ Обозначим соответствующий уровень накопленной инфляции . Дальнейший переход на новые оптимальные производственные программы не произойдет в силу того,

что

dFL (Е) > dFj (Е)

для всех j= 1, 1 + 1,..., L - 1 и

для всех Е £ (ЕL,да).

Таким образом, доказано следующее утверждение: пусть множество 0 = (01,...,0) задает множество всех допустимых решений (1) - (4), (9); в том случае, если маржинальный доход по каждому виду выпущенной продукции растет по линейному закону Р1- (Е) = Р; + щ Р; Е, 1 = 1, 2,.,Мотносительно накопленной инфляции Е, то интервал изменения инфляции Е К е (0,да)] может быть таким образом разбит на конечное число отрезков (0, Ех), (Ер Е2)- • (Е^, да), что при изменении инфляции в рамках одного отрезка (Е р-1, Е р ) оптимальной остается одна и та же производственная программа е 0.

Интервальное задание маржинального дохода. Рассмотрим задачу (1) - (4), (9) в условиях, когда маржинальный доход по каждому виду выпускаемой продукции может принимать любое значение из заданного интервала [в1,Р2], т. е. р^ е[р1,р2]. В этом случае возникает проблема: какие решения из множества допустимых производственных программ 0 = {01,...,0L} могут быть оптимальными и на каком множестве возможных значений Р - они остаются оптимальными. Иными словами, можно ли разбить многомерный параллелепипед значений маржинальной доходности выпускаемой продукции

п

Р = П[в,в2] на области 01,., Оутак, чтобы вы-

1 =1

полнялись следующие два условия:

- 1^0/ = Р;

/=1

Очевидно, что при уровне накопленной инфля- - для каждой 0/ существует некоторая произ-

ции Е более чем Е1 производственная программа Q1

водственная программа Qj, которая остает-

ся оптимальной для всех значений Pe QJ, Р=Р1, Р2,..., Рм .

Решение этой задачи сводится к организации следующей процедуры.

1. Формирование множества потенциально оптимальных производственных программ. Для этого строится множество интервалов [ F/, F2j ] j = 1, 2,., L, где Fx3 - значение целевой функции (1) на производственной программе Q1 при условии, что маржинальный доход р,. = р'. (, = 1, 2,.,MM); F2 - значение целевой функции (1) на производственной программе Qj при условии, что маржинальный доход р,. = р2.

2. Формирование множества потенциально оптимальных производственных программ Quo с Q. Для этого необходимо сделать ряд шагов.

2.1. Определить maxFj = F24 (1 <l1 <L).

j=1,L

Включить в Q производственную программу Q4.

2.2. Определить maxF1 = F2 (1 <l2 <L).

j=1, L

Включить в Q производственную программу Q.

2.3. Включить в Qno все производственные программы QJ с Q, для которых F/ > F^.

3. После того как сформировано множество потенциально оптимальных производственных программ Q , определим, для каких значений рj e F (j = 1, 2,., M) производственная программа Q' является оптимальной. Легко видеть, что область изменения маржинального дохода Р^, на которой оптимальной будет производственная программа Q', определяется следующей системой неравенств:

'Р1, <Р, <Р2, j = 1,2,...,М

J м м (34)

IIP ,Qj jQ, 1 *'vQ1 ^ QL.'

.,=1 ,=1

Таким образом, система неравенств (34) задает многоугольник, при применении на котором значений маржинального дохода Pj/ = 1, 2,., M) оптимальной остается производственная программа

Q'.

Нелинейные модели управления производственными инвестициями

Рассмотрим подход, позволяющий разработать оптимальную стратегию использования заемных

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

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

Динамическая детерминированная модель учитывает фактор времени. Речь идет об интенсивности выпуска продукции в определенный момент времени, что может быть использовано при сезонных колебаниях спроса на продукцию. Также в модель заложено влияние внешних факторов на маржинальный доход предприятия и интенсивность спроса.

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

Динамическая детерминированная модель управления кредитными ресурсами. Данный вид модели содержит фактор времени для учета изменения основных параметров в различные периоды. В этом случае выражение целевой функции преобразуется к следующему виду: у, [Х, х, у(Х)] -маржинальный доход от единицы реализуемой продукции , в момент Х при выбранной производственной программе X = (Х1,..., Хп) и внешних условиях в момент Х, заданных вектором у(Х). К внешним условиям можно отнести такие факторы, как уровень инфляции, уровень безработицы, цены на энергоносители. С помощью данной задачи определяется х,(Х) - интенсивность выпуска (реализации)

продукции вида 1 в момент времени t. Тогда объем

выпуска (реализации) продукции 1 на отрезке врет

мени (0,Т) составляет х1 = |х1 (-)dt.

Здесь использовались следующие обозначения: L. - объем имеющихся материально-сырьевых ресурсов вида /; z. - объем дополнительных материально-сырьевых ресурсов вида /; У/ - стоимостная оценка материально-сырьевых ресурсов; в/[ z, х, у(-Е)] - стоимость материально-сырьевых ресурсов, зависящая от объема выпуска, объема закупок и момента в который эти закупки производились; т - эффективное время использования единицы оборудования вида на период планирования; ki - число единиц оборудования вида /; Ipti [-, х, у(t), (-, х, у(-)] - интенсивность спроса на продукцию 1, где у [-, х, у(-)] - функция, выражающая интенсивность продаж.

Учитывая описанные параметры, модель можно формализовать следующим образом: - целевая функция максимизации маржинального дохода

X {v,[t,х,y(t)] xt(t)dt ^ max; (35)

• T

j xt (t )dt;

1=1 о

x. =

X ljx, < Lj + Zj, j =1m;

(36)

условие на равномерную загрузку мощностей п г t

Xт(х) |х1 (t№' <Т1,Vt е (0,Т); (37)

1=1 0 т

ограничение на интенсивность использования производственных ресурсов

XT■1 (х)х < т,к,, 1 = 1,1, V- е (0,Т);

/=1

ограничение, учитывающее спрос на продукцию

х1 (-) < р(-,х,у(-),у;.[-,х,у(-)], V- е (0,Т); (38) ограничение на закупки материально-сырьевых ресурсов в пределах объема кредита

XZjPj [z, х,y(t%)] < V, Vt е (0,T),

(39)

j=i

где -е - момент закупки;

- ограничение на целочисленность переменных х е I. (40)

Решение задачи (35) - (40) можно осуществить путем разбиения периода (0,Т) на интервалы

(0, -:),...,(-ж-1,Т), и на каждом интервале решать нелинейную статическую задачу управления кредитными ресурсами. Далее решение динамической задачи получается путем слияния решений на всех временных интервалах вида (-/, +1). Обоснованность данного подхода базируется на теореме о том, что любую гладкую непрерывно дифференцируемую функцию можно со сколь угодно большой точностью аппроксимировать кусочно-постоянными функциями при выполнении условия

i 1 jv - (t )dt -jv K.n.(t)

<s,

(41)

- ограничение на объем материально-сырьевых ресурсов

где VK.n.(t) - кусочно-постоянная функция.

Рассмотрим механизм разбиения подробней. Детализируем детерминированную динамическую модель (35) - (40), сформулированную для любого интервала. Будем считать, что на этих интервалах маржинальный доход Vq[t,х,y(t)], q = 0,1,...,N- 1 -не зависит от времени, а интенсивность производства xf (t) постоянна для любого периода времени t из интервала [tq, tq+1]. Интенсивность спроса и отпускные цены на продукцию в этом промежутке времени также считаются постоянными, т. е. Iptf {t, х, y (t), ai [t, x, y (t)]} = const. Тогда указанную динамическую модель (35) - (40) можно представить в виде статической оптимизационной задачи. Производственная программа выпуска продукции в определенный период времени задается вектором xq = (x1q,...,xf),q = 1, 2,...,N. Задача формализуется следующим образом:

- максимизация целевой функции

X vf [xfy(t)]x,q ^ max; (42)

1 =1

- ограничение на объем материально-сырьевых ресурсов с учетом запасов прошлого периода

(43)

X xq lj < j + Zq;

ограничение на интенсивность использования производственных ресурсов

(44)

XTa (xq) xq <T? k;

ограничение, учитывающее спрос на продукцию

Xxq < ptq[xq,yk(t),a,(xq,y(t)];

(45)

ограничение на закупки материально-сырьевых ресурсов в пределах объема кредита м у

X zqJ в) < ж . (46)

1=1

Сделаем несколько замечаний к данной модели:

1) решая задачу (42) - (46), для любого я = 0, 1,.., N

ХЧ РХ Ч

примем хч (Х) = —!-, а 1рХ9г (Х) = —^;

г АтЧ АтЧ

2) если на интервале (ХЧ, ХЧ+1) кредит использован

V

в размере меньшем, чем — то остаток пере-

К'

носится на следующий период;

М

3) х, (Х)Ев- интенсивность денежного потока,

з=1

при интенсивности выпуска х. (Х), тогда если кредит выдается не единожды, а траншами, суммарный объем которых составляет величину V, ограничение (46) перепишется в виде

N М

Е X, (Х)ЕвА < V,(Х), для любого Х е (0,Т), (47)

,=1 .=1

(реализации) продукции i на отрезке времени (0, Т)

1

J vt (t )dt = V.

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

Здесь у, [Х, х, у(Х)] - случайная функция маржинального дохода при продаже одной единицы продукции вида ,, , = 1, 2,..., п. Как и в предыдущей модели, у (Х) - вектор внешних условий. Тогда целевая функция модели имеет вид

n

I J у, [t, х, y(t)] xt (t)dt ^ max,

(48)

i=1 0

где у, [Х, х, у(Х)] - математическое ожидание случайной функции у, [Х, х, у(Х)], рассматриваемой как случайный процесс;

у, [Х, х, у(Х)] х1 (Х) - интенсивность денежного потока в момент Х

Как и в предыдущих моделях, объем выпуска

х=

: | xi (х^х.

0

Ограничение на объем материально-сырьевых ресурсов

Е^. < ^ + ,. = (49)

,=1

Ограничение на производственные мощности Ет(x)х, <тА,i = 1Д (50)

Ограничение, учитывающее спрос на продук-

цию

X,(Х) < 1рХ,[Х,f(Х,X,у(Х)], УХ е (0,Т). (51) Ограничение на приобретение дополнительных материально-сырьевых ресурсов в пределах объема кредита

I j j [Z x, ) ]< V.

(52)

j=i

Это ограничение можно представить в не-

п

сколько ином виде. Пусть ZХ, =ЕI .в. [2,х,У(Х^)] -

=1

переменные затраты на производство единицы

п

продукции вида ,. Тогда ZХ = Е ZХ i х, < V - огра-

,=1

ничение на покупку ресурсов в пределах кредита, xlZХl

а и { = - доля кредита, потраченная на приобретение материально-сырьевых ресурсов для производства продукции вида , за весь период.

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

а) в пределах допустимого значения

п

Е°2 [Х, X, У(Х )]и2 +

+2IIcovj[t,x,y(t)] и ,uj < D, (53Л)

i=1 i>j

где о - дисперсия случайной величины у - маржинального дохода i-го актива,

I и, < 1, и, > 0, Vt е (0,T);

i=i

covj - ковариация маржинальной доходности

i-го и j-го продукта;

б) согласно минмаксному критерию

n

min max Va2[t, x, y(t)] u2 +

x t

,=1

n

+2IIcovj[t,x,y(t)] u,и^ < D. (53.2)

,=1 i >j

i=1

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

Если в модели (48) - (53) случайными являются еще и параметры в/, Т;1 и р-;, то дополнительно вводятся следующие ограничения:

Б(в /) < А, (54)

ЩТа (х)] < А (55)

D{Ipt[t,у (-,х,у(-)]} < Dз(-), - е (0,Т), (56)

где D{IptI.[t, у (-, х, у(-)]} - дисперсия случайной функции 1р-[-, , х, у(-)]. Данная модель - многокритериальная задача максимизации функционала (48) при ограничениях на количественные оценки риска доходности производственной программы, риска недостаточности производственной мощности при выпуске продукции, по выбранной производственной программе, риска недопроизводства, риска недостаточности кредитных ресурсов для обеспечения производства необходимыми материальными ресурсами.

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

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

производственной программы X = (Х1,..., Хп),

т

интенсивностей реализации х1 (-)/ х1 =| х1 (-^ и

0

вектора закупок Z, которые максимизировали бы функционал (48) при ограничениях (49) - (56).

Пример расчета оптимального использования производственных ресурсов в динамической модели управления производственными процессами

Рассмотрим пример применения детерминированной модели управления производственными ресурсами Красноярского металлургического завода «Старт». Схема производства и виды выпускаемой продукции завода представлены на рис. 2.

Пи волока для пор^го ц^нкоеажя

Г>в волока щишов^ная

ГТэоволока ДЛЯ гвоздей

Гвозди

Г)эовола<а для ятдаюй еьсэдки

Сймиарезакшуе ЕЖТЫ

Пэсеолока для СЕТКИ

(Ьтка плетеная

Г^ювслока оц|+коеэньвя

л 9

Г(в волока арлтатурная класса В 500

ГТэо волока от ожжен пая светгая

ГТэоволска обшрпо нззнанешя, СЕароная

Рис. 2. Схема производства и виды продукции КМЗ «Старт»

Размотка Чистка Отжиг Охлаждение Намотка Готовая

w W W продукция

Рис. 3. Укрупненная схема производства проволоки

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

Приведем основные параметры производства и реализации продукции, необходимые для применения модели. Маржинальный доход от реализации единицы продукции представлен в табл. 1.

Технологический процесс изготовления проволоки включает ряд классических операций (рис. 3).

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

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

В распоряжении завода находятся 24 намоточных станка, производительность которых представлена в табл. 2.

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

8,3 х 12 х1 + 6,8 х 14 х2 +

+5,6 х 19 х3 + 6,5 х 15,5 х4 ^ max,

Таблица 1

Маржинальный доход по видам продукции, тыс. руб./т

Таблица 2

Производительность намоточных станков, т/день

Продукция Производительность одного станка

Оцинкованная проволока 12

Отожженная проволока 14

Сварочная проволока 19

Арматурная проволока 15,5

xi + х2 + хз+ x4 ^ 24, X е Z +, 3 = 1,.., 4.

Решить эту оптимизационную задачу можно симплекс-методом с помощью Excel. При заданных ограничениях значение маржинального дохода будет равно 2 553,6 тыс. руб. При этом решение задачи следующее: x1 = 0, x2 = 0, x3 = 24, x4 = 0. Таким образом, все 24 станка используются для намотки сварочной проволоки. Очевидно, что общий объем готовой к продаже сварочной проволоки равен 456 т (рассчитан как произведение количества станков на их производительность).

Незавершенное производство на конечной стадии производства (намотка проволоки) задано значениями, представленными в табл. 3.

Тогда получаем, что через т1 = 820,8/456 = 1,8 дня полностью закончится незавершенное производство на стадии намотки сварочной проволоки.

Далее переходим ко второму этапу итерационной процедуры. Чтобы продолжить выпуск сварочной проволоки, необходимо выделить ресурсы на стадию охлаждения проволоки. Допустим, данная процедура длится в течение 0,5 дня. Охлаждение проводится в специальных ваннах, вместимость которых ограничена (200 т/день), поэтому учтем, что производительность сварочной проволоки не должна превышать 400 т/день.

Таблица 3

Объем незавершенного производства на конечной стадии VN(0), т

Продукция Маржинальный доход от реализации единицы продукции

Оцинкованная проволока 8,3

Отожженная проволока 6,8

Сварочная проволока 5,6

Арматурная проволока 6,5

•I

Продукция Объем незавершенного производства

Оцинкованная проволока 511,8

Отожженная проволока 213

Сварочная проволока 820,8

Арматурная проволока 709,4

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

8,3 х 12 х1 + 6,8 х 14 х2 +

+5,6х 0,5 х 19 х3 +6.5 х 15,5 х4 ^max,

x1 + x2 + x3+ x4 < 24, 0,5 х 19 х 3< 400,

xt е Z+, 1 = 1,.., 4.

Получим оптимальное решение задачи x1 = 0, x2 = 0, x3 = 0, x4 = 24, т. е. на данном этапе все станки используются для намотки арматурной проволоки. Маржинальный доход равен 2 418 тыс. руб. Данный этап продлится т2 = 709,4 / (24 х 15,5) = = 1,907 дня. То есть через это время закончится все незавершенное производство арматурной проволоки на этапе намотки.

Таким образом, получено оптимальное распределение производственных программ для первых 3,707 дня. Аналогичным образом можно продолжить расчеты до момента времени, когда будет обработан весь объем незавершенного производства на оставшихся операциях производственного процесса.

Список литературы

1. ВанХорн Д. Основы управления финансами / пер. с англ. М.: Финансы и статистика, 2004.

2 . Мищенко А. В. Методы управления ограниченными ресурсами в логистических системах. М.: Инфра-М, 2011.

3 . Шапиро Д. Моделирование цепей поставок / пер. с англ. СПб: Питер, 2007.