УДК 021.8 + 025.1 ББК 78.34
МОДЕЛЬ ЧИСЛЕННОЙ ОПТИМИЗАЦИИ СТРУКТУРЫ ОПЕРАЦИОННОГО ЯДРА ОРГАНИЗАЦИИ
1 2 Воронин А. А. , Харитонов М. А.
(Волгоградский государственный университет, Волгоград)
Построена модель оптимизации операционного ядра организации, которое состоит из базового технологического модуля и сети модулей вспомогательных производств. Производственная функция операционного ядра представлена в виде суперпозиции производственных функций Леонтьева, отвечающих каждому из модулей. Задача оптимизации сведена к задаче линейного программирования с параметром, описывающим структуру операционного ядра. Разработан алгоритм автоматического построения исходных уравнений для каждого значения параметра. Представлены результаты численного решения оптимизационной задачи в широком диапазоне переменных и параметров модели. Данная модель может служить основой синтеза обобщенных механизмов управления организационной системой на большом интервале времени.
Ключевые слова: организационная система, операционное ядро, оптимизация структуры, производственная функция, линейное программирование.
1. Введение
Построение моделей управления организационной системой (ОС) на большом интервале времени на базе представлен-
1 Александр Александрович Воронин, доктор физико-математических наук, профессор (a.voronin@volsu.ru).
2 Михаил Алексеевич Харитонов, студент (kharitonov. mihail@gmail. com).
ной в [10] системы механизмов управления невозможно без учета ее структурно-функциональной изменчивости. Изменчивость внутренних и внешних связей, краткосрочных, среднесрочных и долгосрочных целей и критериев развития ОС мотивирует обобщение не только моделей базовых механизмов [10], но и обобщение термина «механизм управления» путем введения в его состав в дополнение к собственно методам принятия решений также и комплекса организационных структур, норм, параметров ОС, изменяющихся в течение долгосрочного периода управления согласно актуальной и проектируемой динамике ОС с учетом реализации временной последовательности простых и комплексных механизмов управления [10]. Такие обобщенные механизмы управления представляют собой сложноструктурированные системы, подсистемами и элементами которых являются обобщенные, комплексные или простые механизмы, оптимальность каждого из которых, таким образом, определяется уже в контексте задачи оптимизации всей системы - обобщенного механизма управления.
Эта концепция фактически совпадает с концепцией организационной структуры Г. Минцберга [8], поэтому термины «организационная структура» (в смысле Г. Минцберга) и «обобщенный механизм управления» являются синонимичными, а структурные конфигурации Г. Минцберга можно назвать оптимальными (в рамках формалистики менеджмента) обобщенными механизмами управления. Структурные и организационные изменения ОС можно таким образом трактовать как переходы между условно оптимальными структурными конфигурациями или обобщенными механизмами управления, выделяя в них (в соответствии с введенными признаками) траектории адаптации, развития и др. Развитие ОС при этом можно связывать с взаимодействием и развитием ее обобщенных механизмов управления.
Этапами синтеза обобщенных механизмов, т.е. этапами решения задач моделирования структурной динамики и динамической структурной оптимизации ОС являются построение семейств соответствующих уравнений динамики и критерия эффективности в структурнозависимой форме, описание пространства структурных состояний, корректное введение «прин-
ципов» динамики, разработка эффективных алгоритмов оптимизации. В качестве первых успешных шагов в этом направлении можно указать на работы [1-3, 5, 6, 9], в которых исследуются задачи статической и динамической оптимизации иерархий управления ОС с критерием эффективности в виде структурно зависимой функции затрат. Для синтеза обобщенных механизмов управления необходимо сопряжение этих моделей со структурно зависимой производственной функцией (ПФ) ОС с теоретико-игровыми моделями иерархического управления. Сказанное мотивирует актуальность синтеза системы многопараметрических динамических моделей оптимизации и управления, в совокупности образующих иерархическую структуру обобщенных механизмов управления ОС.
В настоящей работе представлена математическая модель оптимизации операционного ядра организации, под которым согласно терминологии Г. Минцберга [8] будем понимать совокупность структурных единиц ОС, непосредственно осуществляющих основное и вспомогательные производства (включая персонал и оборудование). Структура операционного ядра представлена в виде ориентированного графа специального вида с переменным числом слоев. Вершины графа обозначают технологические модули базового и вспомогательных производств, ребра - входы и выходы технологических модулей -«потоки» факторов производства. Целевая функция - производственная функция операционного ядра - представлена в структурнозависимой форме в виде суперпозиции производственных функций Леонтьева, отвечающих каждому из технологических модулей.
При заданных входах (факторах производства) и структуре операционного ядра значение производственной функции является решением задачи линейного программирования.
Увеличение числа вспомогательных производств, с одной стороны, обеспечивает устойчивость выпуска при нарушении оптимальных пропорций между факторами производства, с другой - увеличивает себестоимость. Альтернативой усложнения структуры операционного ядра является покупка недостающих факторов на рынке (при наличии такой возможности). В экономиках с высоко развитой иституциональной средой
последний способ оказывается наиболее эффективным, что приводит к упрощению структуры ОС и появлению специализированных аутсорсинговых организаций, реализующих не только вспомогательные производства, но и ранее традиционные функции ОС («универсальная бухгалтерия», «универсальный офис-менеджер» и т.п.). А в экономиках с низким уров-немразвития институциональной среды преобладает противоположная тенденция - вертикальная интеграция путем поглощения смежных и вспомогательных производств. Однако в предложенной в настоящей работе модели не исследуется задача поиска оптимального баланса между указанными тенденциями, и оптимизация структуры операционного ядра заключается в поиске минимального числа слоев ее структуры, обеспечивающего максимальное для заданных входов (факторов производства) значение производственной функции (то есть не увеличивающееся при дальнейшем наращивании числа слоев).
Для решения задачи разработан алгоритм автоматического построения исходных уравнений для каждого числа слоев. Представлены результаты численного решения задачи в широком диапазоне переменных и параметров модели.
2. Структура и производственная функция операционного ядра ОС
Построенная в работе модель структурно зависимой ПФ основана на моделях [4, 6]. В качестве элементарной ПФ каждого элемента структуры операционного ядра с неизменной технологией (простого преобразователя - ПП) используется ПФ Леонтьева F = kтт^/а\, ...,УП/а„) где fi - величины аргументов - трансформационных факторов производства; аь ..., ап -технологические коэффициенты; k - нормирующий множитель. В настоящей работе будем использовать агрегированную макроэкономическую трехфакторную модель трансформационных факторов (природный - N, технический - Т, человеческий - Н). Учет векторной природы факторов производства влечет лишь увеличение числа переменных и уравнений модели.
В условиях устойчивых межфакторных диспропорций или флуктуаций трансформационных факторов среднее значение ПФ IIII значительно ниже ее максимума. В рамках ОС эластичность ПФ по аргументам обеспечивается построением вертикальной структуры вспомогательных производств, производящих недостающие для максимальной эффективности базового II части трансформационных факторов. Таким образом, структуру факторных производственных потоков операционного ядра ОС можно представить в виде сложного преобразователя (СП) -многоуровневой структуры II, в котором выходы одних II являются входами других. Структура СП имеет вид ориентированного графа, вершинам которого отвечают ПП, а ребрам -факторные потоки. Последняя вершина отвечает базовому, а остальные - вспомогательным производствам. Исходные факторы производства являются входами всех ПП, выходом последнего (базового) ПП является продукт ОС.
Задача оптимизации факторных пропорций ПФ ПП имеет
вид
(1) F = Smin\ N,T, — max, N + T + H = R, S = A + B + C
\ ABC J nt ,—
и очевидное решение:
n=RA, t=RB, — = RC, f=r.
S S S
На рис. 1 приведена структура ПФ, на рис. 2 - графики ПФ
ПП для оптимальных и неоптимальных факторных пропорций. На рис. 2 приведены некоторые графики функции (1).
N Т Н
Рис. 1. Простой преобразователь
а) Ь)
Рис. 2. Графи^(1, Т, И): a)A = B = C = 1; Ь)А = 1; В = 2; С = 3
3. Задача оптимизации структуры операционного ядра ОС.
Рассмотрим СП с факторами производства N НЩТ^
HN соответственно части факторов Т и Н, участвующие в производстве недостающей части фактора ^и возможностью производства части фактора N на п последовательных 1111 (рис. 3).
Рис. 3. Сложный преобразователь с производством части фактора N
ПФ этого СП задается следующей системой уравнений:
£№ + N0= N, £Т + тр = тм, Xн + Ир = ИМ,
I=1 I=1 I=1
/■
г-1
т я.
Л
(2)
N = I N1 = (« + Ь + с)Шп -£N1 ,-і
ь С у
к=г+1
V «-г=0
^ = («п + Ьп + Сп )Ш1П
^ 1 п-1 т я ^
1 I N п
п ’ 7 ’
«п к=0
Ьп Сп
у "У 5 ’ С
Здесь первые три уравнения уравнения описывают распределение исходных факторов производства по структуре СП. Четвертое уравнение дает балансовые соотношения для /-го вспомогательного ПП; N - общая величина фактора, произведенного на /-м ПП; - величина фактора, произведенного на /-м ииспользуемого на к-м ПП, / = 1, ..., п - 1. Пятое уравнение системы отвечает последнему вспомогательному ПП (/ = п).
Из переменных системы (2) составим вектор факторных потоков СП
Ф
N N0,..., N0,Ii, Hi, Nl, Nl,
..., Nn -1, I„i, Hni, NF1, NF, If , Hf ,pp
., <, Io, H 2,.
Задача оптимизации ПФ СП (2) с переменной структурой при фиксированных N, TN, HN имеет вид:
( м" + N0 т
(3) ^ = 5тп + М т
Б ’ C
^ max.
Ф ,n
Заметим, что задача (1) эквивалентна следующей задаче линейного программирования (ЛП):
S S S
---N + F < 0, — T + F < 0, — H + F < 0, F ^ max.
ABC n,t ,н
Таким же образом задача (2)-(3) сводится к следующей задаче параметрического ЛП:
p
(4) фр(n) ^ max, £Мг}Ф} < А,i = 1,-,u(n),
«,Ф ;=!
где р = p(n) - число переменных; u(n)- число неравенств, определяемых следующими формулами:
p(n) =
n + 5n + S 2
, u(n) = Зп + 6.
В частности, для n
1,a,
b, = с, = 1 имеем:
M =
г i G G G 1 G G G "
G 1 G G G 1 G G
G G 1 G G G 1 G
G G G -З -З G G 1
G G G G G -З G 1
G G G G G G -З 1
-З G G 1 G G G G
G -З G 1 G G G G
G G -З 1 G G G G
Ф = ( N0,71, H1, NF, NF,If , HF ,pP ), в = ( N ,m, HN ,0,0,0,0,0,0).
Обозначим фр*(п) - решение задачи (4). Оптимальной структурой СП, описываемого системой (2), будем считать величину
n = min I Arg max p* I.
n \ n J
Обозначим X = (N + TN + HN)/F. Внешними параметрами задачи (4) являются исходные межфакторные пропорции N: TN: HNи пропорции технологических коэффициентов ai : bi : c. Для каждого СП величины ai,bi, ^предполагаются заданными. В численных экспериментах эти величины задавались случайным образом, с сохранением их значений при варьировании числа слоев.
Результаты решения задачи (4) для некоторых значений переменных N, TN, HN и параметров a, bi, ci, принимающих случайные значения из отрезка K с R+, представлены на рис. 4.
Рис. 4: График X(n): a) N = 1, TN = HN = 10, K = {1}, n* = 2, b)N = 1, TN = HN = 10, K = [1, 2], n* = 3
Структура полного СП с возможностью вспомогательного производства каждого из трансформационных факторов представлена на рис. 5.
ПФ СП на рис. 5 состоит из трех СП (см. рис. 3), каждый из которых соответственно производит недостающие части факторов Ы, Т, Н. На рис. 5 обозначено: ТЫ,НЫ- соответственно части факторов Т и Н, участвующие в производстве недостающей части фактора Ы; ТТ,НТ - соответственно части факторов Т и Н,
участвующие в производстве недостающей части фактора Т; ТН,НН - соответственно части факторов Т и Н, участвующие в производстве недостающей части фактора Н.
(5)
Рис. 5. Полный сложный преобразователь ПФ этого СП имеет вид:
* = Х шт ( N, Т, Н
У ЛВС N = Ы* + N0,
Т = Тп + Т*. + ТЫ* + ТТ* + ТН*,
Н = Н*3 + Н* + НЫ* + НТ* + НН*.
+ №р = ы, + тыр = ты, ^т1 + ныр = ны,
і=1
і=1
і=1
АП '”1 АП / 1 г,1 г.Л - Г 1 ^ ТЫ{ НК1
Ы = Х К = (а + Ь) + Ь))тіП — ж,~гг, г
ч а і=о
ь2 ’ с,2
ЫР2=(< + ь;+ <)тт
( 1 ”1 -1 , НЫ„ л 1 ■у' ы1_______________”2_________”і
1 ^ і ’ г.2 ’ 2
а” к=0 Ь С
ч ”1 ”1 ” у
I:
Ет0 + tF = ть, Етт + ttf = тт, EhT + htf = ht ,
l=1 l=1 l=1
n2 / ч Г 1 І-1
тг - Е T'k = (a2 + b2 + b2 )mln -2 £?
TT2 HT2Л
к=i+1
v a l=o
b2 с2
^ =(an2 + bn2 + c2)mln
1 n2-1 TTn2 HTn к n2 n2
г У
2
T- Е N
2 n2 2
an к=o bn
v n2 n2
n у
III:
EH + н* = Hb, E?Hl + thf = тн, Енн/ + нн; = нн,
і=1 і=1 l=1
h: - Е нк=(a; + Ы + Ы) mm Г4; eh, H, SLл
к=і+1
v a l=o
HFl =(al + bi + сі )mln
F v n n П '
^ 1 n -1 THn HHn Л
1 ттк n n
Е П^'
і Lu n; ,i
an к=o bn
V Пі Пі
"і у
В обозначении технологических коэффициентов нижний индекс n, j = 1, 2, 3, обозначает номер слоя в j-м СП, что соответствует нумерации I, II, Шна рис. 5. Верхний индекс (так же как и j) обозначает номер СП.
При заданных входах задача оптимизации структуры ОЯ с ПФ (5) имеет вид
^ тг с ■ Г NTH Л
(6) F = S mrn I—,—,— 1^ max,
I ABC) Ф ,n’
Ф = (Ф1, Ф 2, Фі, Фр), где
Ф1 =
(No,N”,...,N°,TN1,HN1,N1,N1,...,N\,TN2,HN2,
Ф2 =
v...,Nn-1,TNni,HNn,NF1,NF,TNf,HNF
^o, to , .,Tn02, TT1, h?1, t'?1 ,., t12 ,tt2, HT2,. ...,t”2-1,tt ,нт ,тП2,t,o,tt.,HT.
2
2
Фі =
г Hi H 20,
,S TH1,HH1,H2,H1,..
,H1 ,TH2,HH2,...,Hn~X,TH ,HH ,HF3,Hl,THF,HHF
’Пз’ 2 5 2 5 «3 ’ «3’ П3 5 F’F’ F 5 F ^
n = (n1, n2, n3) где n1,n2,n3- число слоев каждого СП (I, II, II).
Сводя задачу построения полного СП (6) к задаче линейного параметрического программирования с векторным параметром n, получим:
p
(7) Фр(n) ^ max, «Ф/ ^ ^, г' = 1,...,U(n),
n ,ф j=i
где P = P(n) - число переменных; U(n)- число неравенств, определяемых следующими формулами:
P(n) = p(n1) + p(n2) + p(n3) + 1, U(n) = u(n1) + u(n2) + u(n3);
N ,TN, HN ,0,...,T ,TT, HT ,0,., H ,TH, HH ,0,
Г M1
u ( П1)
oo
0 M2 0 f2
0 0 M F,
(матриц^іМ1, M2, Жістроятся аналогично задаче (4)),
F1= (0,0,0,A,B,C,0,..,0І F2= fo,0,0,B,B,C,0,...,0
1 I S S S у 2 V S S’S
;і= I 0,0,0,C,B,C ,0,...,0|, dlmF - u(n,), i-1,2,і, S - A+В+C
і I SSS у г
Для исследования зависимости решения задачи (7) от векторного параметра был разработан алгоритм автоматического вывода неравенств (7) для различных n.
Обозначим Фр* - решение задачи (7), Q - Arg max Фр (n).
n
Множеством n оптимальных структур полного СП для задачи (7) будем считать множество Парето пп многокритериальной задачи оптимизации:
n., ^ mln, i - 1,2,і.
neg*
T
Y
4. Результаты численных экспериментов
Внешними параметрами задачи (7) являются пропорции входных переменных N : Т : Ни технологических коэффициентов вспомогательных ПП. Далее представлены результаты решения двух частных оптимизационных задач (7):
7(а): задача (7) при фиксированных входах.
7(б): задача (7) при п = п1 = п2 = п3.
Некоторые результаты решения параметрической задачи 7(а) симплекс-методом [11] представлены на рис. 6.
а) Ь)
Рис. 6. График зависимости А(пь п2, п3) при щ = 10,
Ыь = Ть = Нь = 1,ты = ны = тт = нт = тн = нн = 10 а) К = {1}, пп = (10, 1, 1); Ь) К = [1, 2], пп = (10, 8, 4); с) а, = 1, Ь = 2, сг = 3, Пп = (10, 2, 1), ф) К = [1, 5], пп = (10, 10, 4)
На рис. 7, 8 представлены результаты решения задач 7(а, б) для некоторых значений параметров.
На рис. 10-12 приведены некоторые решения задачи 7(б): графики изменения X при ТЬ є [0, 1], НЬ є [0, 1] с шагом А = 0,1 и значениях параметра п = щ = п2 = п3от 1 до 8.
Рис. 7. График зависимости Х(п) а)Ыь = Ть = Нь = 1, ТЫ = НЫ = ТТ = НТ = ТН = НН = 10, К = {1}, п = 2; Ь)Ыь = Ть = Нь = 1,
_ _____ ___ *
ТЫ = НЫ = ТТ = НТ = ТН = НН = 100, К = {1}; п = 4
Рис. 8. График зависимости Х(п). п = 5; приЫь = 1,0 иТЫ = НЫ = ТТ = НТ = ТН = НН = 10, К = [1, 5]. а)Ть = 0,3; Нь = 0,3; Ъ)Ть = 0,3; Нь = 0,6; с)Ть = 0,3;Нь = 0,9; Ф)Ть = 0,6; Нь = 0,3
0 2 4 п 6 8 10 с)
Г
0 2 4 П 6 8 9 10 е)
Г , і—■
4 Ц 6 8 10
5'
0 2 4 П 6 8 10 Ф
5- Г г-1-
“0 2 4 11 6 8 10
П 6 8 10
Рис. 9: ГрафикзависимостиЛ(п). п = 5; приЫЬ = 1,0 иТЫ = НЫ = ТТ = НТ = ТН = НН = 10, К = [1, 5]. а) Ть = 0,6;Нь = 0,6; Ь) Ть = 0,6;Нь = 0,9; с) Ть = 0,9;Нь = 0,3; ф) Ть = 0,9;Нь = 0,6; е) Ть = 0,9;Нь = 0,9; /) Ть = 1,0;Нь = 1,0
Рис. 10. Графи^(1, Т, Н)при ТЫ = НЫ = ТТ = НТ = ТН = НН = 5,0 иприК = {1}
п=1 п=2
Рис. 11. График F(1, Т, Н) при ТЫ = НЫ = ТТ = НТ = ТН = НН = 10,0 иприК = [1, 5].
п=5 п=б
Рис. 12. Графи^(1, Т, Н)при ТЫ = НЫ = ТТ = НТ = ТН = НН = 10,0 иприК=[1, 5]
5. Заключение
Результаты численного исследования модели оптимизации структуры операционного ядра ОС показали, что
- значительные запасы факторов во вспомогательных производствах и достаточно большое число 1111 во всех линиях полного СП позволяет достичь устойчивого максимума ПФ при значительных отклонениях пропорций факторов базового производства от оптимальных и разбросах технологических коэффициентах вспомогательных производств.
- увеличение разброса технологических коэффициентов вспомогательных производств при больших запасах факторов компенсируется ростом числа слоев ПП.
- уменьшение запасов факторов и числа ПП вспомогательных производств резко снижает устойчивость выпуска при отклонении пропорций факторов базового производства от идеальных.
Введение в данную модель параметров, характеризующих внутреннюю и внешнюю среду ОС, позволяет строить на ее основе модели оптимизации комплексной структуры ОС и механизмов управления. Второе направление непосредственного развития модели - модели структурной динамики и оптимального управления. Синтезом двух указанных направлений являются модели управления адаптацией и развитием ОС.
Литература
1. ВОРОНИН А.А., МИШИН С П. Алгоритмы поиска оптимальной структуры организационной системы // Автоматика и телемеханика. - 2002. - №5. -С. 120-132.
2. ВОРОНИН А.А, МИШИН С П. Оптимальные иерархические структуры. - М. ИПУ РАН, 2003. - 214 с.
3. ВОРОНИН А.А., МИШИН СП. Модель оптимального управления структурными изменениями организационной системы // Автоматика и телемеханика. - 2002. - №8. -С.136-150.
4. ВОРОНИНА И.Д. Задача управления организационной структурой в условиях глобального инновационного процес-
са // Управление большими системами. -2006. - №12-13. -С.51-59.
5. ГУБКО М.В. Математические модели оптимизации иерархических структур. - М.: ЛЕНАНД, 2006. - 264 с.
6. ГУБКО М.В. Структура оптимальной организации континуума исполнителей // Автоматика и телемеханика. - 2002. -№12. -С. 116-130.
7. КЛЕИНЕР Г.Б. Производственные функции. Теория, методы, применение. - М.: Финансы и статистика, 1986. - 239 с.
8. МИНЦБЕРГ Г. Структура в кулаке. Создание эффективной организации / Пер. с англ. под ред. Ю.Н. Каптуревско-го. - СПб.: Питер, 2002. - 512 с.
9. МИШИН С.П. Оптимальные иерархии управления в экономических системах. - М.:ПМС0ФТ,2004. - 205 с.
10. НОВИКОВ Д.А. Теория управления организационными системами. - М.: Московский психолого-социальный институт, 2005. - 584 с.
11. ПАПАДИМИТРИУ Х., СТАЙГЛИЦ К. Комбинаторная оптимизация. Алгоритмы и сложность. Пер. сангл. - М.: Мир, 1985. - 512 с.
CONSTRAINED OPTIMIZATION MODEL OF ORGANIZATION'S OPERATING CORE
Alexander Voronin, Volgograd State University, Volgograd, Doctor of Science, professor (a.voronin@volsu.ru).
Mikhail Kharitonov, Volgograd State University, Volgograd, student (kharitonov.mihail@gmail. com).
Abstract: We suggest a constrained optimization model of an organization’s operating core, whose structure consists of a basic technological module and modules of support facilities. The production function of the operating core is represented as a superposition of Leontief production functions corresponding to each of the modules. We reduce the optimization problem to a linear program with a parameter describing operational structure of the core. We develop an algorithm of automatic construction of basic equations for each value of the parameter. We solve the problem numerically for a wide range of model variables and parameters. The proposed model may serve as the basis for generalized mechanisms synthesis of organizational system control on the long time interval.
Keywords: organizational system, operating core, structure optimization, production function, linear program.
Статья представлена к публикации членом редакционной коллегии В. В. Клочковым