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

УДК 519.6

Е.Н. Барышева, Л.А. Сараев*

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ СТРУКТУРНЫХ ИЗМЕНЕНИЙ ПЕРСОНАЛА ПРЕДПРИЯТИЙ

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

Ключевые слова: система бонус-малус, производственная функция, качество персонала, класс, тарифная сетка, стационарное распределение, предельное распределение.

1. Постановка задачи

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

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

В работе [2] рассмотрена система бонус-малус (СБМ), применяемая в автомобильном страховании. Суть СБМ заключается в распределении водителей по классам качества и пересмотре установленного класса по итогам каждого года. Основное достоинство СБМ заключается в минимуме информации: по каждому водителю нужно знать класс качества в начале года и количество страховых случаев (аварийных ситуаций) в течение года.

В настоящей работе предлагается рассмотреть возможность применения СБМ для оценки качества персонала в целях прогнозирования расходов на оплату труда, влияния на производительность труда и управления структурой персонала, понимая под структурой распределение персонала по классам.

Более подробно рассмотрим следующие вопросы.

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

Во-вторых, в производственную функцию предприятия, описывающую процесс производства, входят два фактора — капитал и труд. Труд учитывается через количество трудовых часов или, в размере их стоимости, последнее позволяет в некоторой степени отделять квалифицированный труд от неквалифицированного. Изменение распределения персонала по категориям должно иметь адекватное отражение в производственной функции. Можно ли считать постепенное

* © Барышева Е.Н., Сараев Л.А., 2011

Барышева Евгения Николаевна ([email protected]), Сараев Леонид Александрович ([email protected]), кафедра математики и бизнес-информатики Самарского государственного университета, 443011, Россия, г. Самара, ул. Акад. Павлова, 1.

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

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

В-четвертых, каким образом оценивать качество персонала (ненаблюдаемая величина) через наблюдаемые параметры — дисциплина, тестирование, производительность, выполнение плановых или контрольных заданий и прочие наблюдения.

2. Градация персонала по классам и премиальная шкала

Тарификация персонала является атрибутом любой бюджетной организации. На предприятиях для основной категории рабочих (служащих) также применяется тарифная или квалификационная сетка.

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

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

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

Предположим, установлено событие, на основании которого принимается решение об изменении класса, назовем это событие как «отклонение показателя от норматива».

Тогда в зависимости оттого, сколько раз было зафиксировано отклонение в течение квартала, принимается решение об изменении класса.

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

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

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

СБМ определяется тремя элементами:

а) премиальной шкалой уь ] = 1, 2,..., К, где К — количество классов;

б) начальным классом;

в) переходными правилами, которые определяют переход из одного класса в другой при условии, что число случаев за отчетный период известно;

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

Система классов эквивалентна в некотором роде тарифной сетке, в частности, размер заработной платы или премии зависит от номера класса, например, для каждого класса возможно применение скидки / надбавки в % от оклада, от основной части заработной платы и т. д.

В данной работе рассмотрим вариант СБМ, когда изначально (например, при приеме на работу или при введении СБМ-системы) служащий, рабочий попадает в 7-й класс и получает 65 % от оклада (величина оклада — это максимальная величина заработной платы, примем ее за 100 %).

При анализе СБМ-системы размер надбавок и скидок можно указывать в общем виде, например, для у-го класса скидка / надбавка установлена в размере / %, где в случае, когда / > 0, это надбавка, иначе скидка.

Другой вариант — установление заработной платы в определенном проценте от базового оклада, ставки для каждого класса.

Премиальная шкала вводимой СБМ может применяться не только для оценки расходов на оплату труда персонала, но и для оценки вклада в производственную функцию или в производительность труда (табл. 1).

Таблица 1

Премиальная шкала

№ п/п 1 2 3 4 5 6 7

Надбавки / скидки, % /1 /2 /з /4 /5 /б /7

Процент от оклада, % У1 ^2 Уз У4 У5 Уб У7

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

F = ^ К аЬв. (1)

Пусть весь персонал распределен по классам, максимальная производительность в том случае, когда весь персонал в 1-м классе, и минимальная, когда в последнем, тогда при одних и тех же численных значениях трудозатрат, выраженных в единицах времени или в стоимостном выражении, обозначим эту величину Ь0, следует учесть распределение по классам.

Наиболее просто это учитывается введением поправочного множителя:

1 к

уь = — V Н:У:. Следовательно, производственная функция (1) с учетом градации пр ; ;

по классам будет иметь вид:

F = КаЬв = ГьКаЬв, ¥ь = уьЕ0. (2)

3. Переходные правила СБМ

Очевидно, что количество возможных СБМ безгранично, и выбор оптимального из них зависит, прежде всего, от цели или целей вводимой системы.

Рассмотрим СБМ, достоинство которой заключаются прежде всего в ее простоте, что облегчает проведение анализа (табл. 2).

Таблица 2

Система бонус-малус

Класс Класс в конце периода после 0, 1 , 2 ... 6 и более случаев

Вначале Ь 0 1 2 3 4 5 >= 6

7 65 6 7 7 7 7 7 7

6 70 5 7 7 7 7 7 7

5 75 4 6 7 7 7 7 7

4 80 3 5 6 7 7 7 7

3 85 2 4 5 6 7 7 7

2 90 1 3 4 5 6 7 7

1 100 1 2 3 4 5 6 7

В данной системе 7 классов, входной класс — это 7-й класс (табл. 2). Наивысший класс — это 1-й класс. Если в течение отчетного периода не зафиксировано отклонений показателей от норматива, то сотрудник переводится в более высокий класс.

Предположим, что в результате исследований получены оценки вероятности наступления события «отклонение параметра от норматива» — #.

На основе данной оценки введем процесс Пуассона как простой поток наступления такого рода событий (в отношении одного лица) с параметром X таким, что вероятность наступления события к раз за период равна:

чк

Р(у = к) = — е Л к!

где V — случайная величина возможного числа наступления событий за период.

Для параметра X, характеризующего интенсивность процесса поступления событий, можно получить оценку вида: так как Р^ = 0) = 1 — Р^ > 0) = 1 — q =

= в—л, то отсюда

— = - 1п(1 - ц).

Переходные правила можно ввести в виде преобразований Тк таких, что Тк(1) = у, если сотрудник переходит из класса Сг в класс С при условии, что зарегистрировано к случаев отклонений от параметра [2].

Преобразование Тк может быть представлено в виде матрицы Тк = (гг(к)), где

^) = 1, если Тк(1) = у и г() = 0 в противном случае.

Вероятность того, что через год сотрудник переместится из класса С{ в класс

С, имеет вид: р1} (—) = £ Рк — ^).

к=0

Здесь Рк (—) — вероятность того, что сотрудник с частотой X будет иметь к случаев в течение года.

7

Очевидно, что р г (Л) > 0 и что Е р^ (Л) = 1 или переместится в другой класс, У 7=1

или останется на месте.

Матрица М (Л) = (р(Л)) = Е рк (Л)Тк , где Тк = Е17) есть переходная матрица

к=0 1,7

цепи Маркова (первого порядка). Это процесс без памяти, состояниями которого являются различные классы СБМ.

В качестве примера формирования матрицы рассмотрим ее элемент

Р12(Л) = Ро(Л)40) + МЛ)^ + Р2(Л)42) +...Р7(Л)47) +... = Ро(Л), так как $5 = 8к,о , где ёу — символ Кронекера.

Так же просто заполняются все остальные элементы матрицы. Переходная матрица будет иметь вид (табл. 3).

Таблица 3

Переходная матрица

М (Л) =

И 7 7 6 5 4 3 2 1

7 1- Ро Ро о о о о о

6 1- Ро о Ро о о о о

5 1 1- Е р 1 1=о Р1 о Ро о о о

4 2 1- Е р 1 г=о Р2 Р1 о Ро о о

3 3 1- Е р 1 г=о Рз Р2 Р1 о Ро о

2 4 1- Е р 1 г=о Р 4 Рз Р2 Р1 о Ро

1 5 1- Е р 1 г=о Р5 Р 4 Рз Р2 Р1 Ро

4. Динамика численности персонала по классам.

Стационарное распределение

Пусть задан вектор распределения персонала по классам, например, это может быть полное сосредоточение всего персонала в 7-м классе в начальный момент времени:

«г(/) = N • wi(t) = N5^ 7, здесь иг(?) — численность персонала в /-м классе, wi(t) — удельный вес /-го класса в общей численности персонала, в момент времени t, где t = к — это номер отчетного периода, так как решение о смене класса производится по итогам периода.

Таким образом, wг (к + 1) = м>т(к)М(А), wт(k), м>т(к + 1) — вектор-строка удельного распределения персонала по классам в начале и в конце периода. Совокупность м^(к), у = 1, 2,..., К есть динамика изменения распределения численности у-го класса.

На ее основе можно определить, например, фонд оплаты труда:

7

Ф(к) = ФоVУ(к)Г} = Ф„ ), ] =1

где Ф0 — базовая ставка заработной платы (точнее ее предельное или максимальное значение).

На ее основе можно определить также значение поправочного множителя,

1 К

входящего в производственную функцию (1): ^ = yLF , уь = — Vп¡/:.

Данная СБМ регулярна, все ее состояния эргодические, то есть можно переходить из одного класса в другой, и цепь не является циклической. В этом случае число 1 — простое собственное значение переходной матрицы М(1). Отметим, что матрица М(1) — стохастическая матрица, собственный вектор для данной матрицы есть в то же время предельное распределение персонала по классам.

Начиная с номера км вектор распределения практически не меняется: у(к +1) « у(к), к > км и, следовательно:

цт (тс) = цт (тс) м (Л) или Мт (Л)ц(тс) = ^(тс).

Отметим, что предельные распределения являются функциями от 1: у(тс) = у(тс;Л) .

Обозначим предельные вероятности распределения персонала по классам:

Ут (тс; Л) = ат (Л) = {а7(Л),а6 (Л),...а1(Л)}.

—т

Соответствующий левый собственный вектор-строка а = {а1(Л, а21(Л,...а7(Л)}

_ _ 7

определяется из уравнений а (Л) = ат (Л)М (Л) и V а1 (Л) = 1 и называется ста-

1=1

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

Из системы ат(Х) = ат(Х)М(Х) находим для заданной СБМ предельные вероятности распределения персонала по классам [1]: а2 (Л) = а1 (Л)(1 - Ро ) / Ро ; а3 (Л) = [а2 (Л) - а! (Л) р1 ]/ ро ;

а4(Л) = [аз(Л) - а2(Л)Р1 - а1 (Л)Ро] /Ро ;

а5 (Л) = [а4 (Л) - а3 (Л)р1 - а2(Л)р2 - ах(Л)Р3]/ Ро;

аб (Л) = [а5 (Л) - а4 (Л) Р1 - аз(Л) Р2 - а2(Л)Рз - а1(Л)Р4]/ Ро ;

а7 (Л) = [аб (Л) - а5 (Л) Р1 - а4 (Л) Р2 - аз (Л) Р3 - а2 (Л) Р4 - а1 (Л) Р5]/ Ро.

7

Так как V а (Л) должна быть равной 1, то в данной системе выбирается

1=1

7

произвольное значение а1(Х), затем производятся суммирование а+ = V а] (Л) и

1=1

нормировка: а] (Л) = а] (Л) / а+ (знак тильда в дальнейшем опускаем).

Предельное распределение не зависит от начального распределения, а зависит от переходной матрицы М(Х) и параметра Х (табл. 4).

Таблица 4

Предельное распределение персонала по классам

Щ/Л о,о1, % о,о5, % о,1, % о,3, % о,5, % 1, %

щ 98,99о 94,744 88,948 59,8о1 26,313 о,722

Щ о,995 4,858 9,355 2о,922 17,о7о 1,241

Щ о,о15 о,369 1,444 1о,3о1 14,987 2,652

Щ4 о,ооо о,о27 о,215 4,938 12,885 5,6о6

Щ5 о,ооо о,оо2 о,о32 2,364 11,о68 11,847

Щ6 о,ооо о,ооо о,оо5 1,132 9,5о9 25,о33

щ7 о,ооо о,ооо о,оо1 о,542 8,169 52,898

Как можно видеть, для большинства реальных значений Л < о,2 предельное распределение характеризуется сосредоточением в первом классе, что сводит на нет значение премиальной шкалы. Скорость сходимости различна, но в целом после 10 периодов близка к предельному.

5. Эластичность СБМ

Одной из целей введения СБМ является поощрение хороших сотрудников (увеличение оклада или премии) и уменьшение размера оклада, зарплаты, премии для сотрудников низших классов.

Очевидно, что размер зарплаты должен быть уменьшающейся функцией от X. В идеальном случае эта зависимость должна быть линейной. Относительное увеличение частоты случаев должно приводить к такому же изменению размера зарплаты. Например, пусть есть два сотрудника с частотой событий 0,10 и 0,11. На большом временном интервале второй должен получить сумму на 10 % меньше первого, СБМ, которая обладает этим свойством, называется совершенно эластичной.

п

Пусть Р(Л) = Е ai (Л)у, — средняя стационарная заработная плата, где у, —

¡=1

премиальная шкала. В идеальном случае увеличение йХ/Х для частоты случаев должно приводить к такому же изменению йР(Х)/Р(Х). Введем эластичность СБМ: П(Х) = [йР(Х)/Р(Х)]/[йХ/Х].

В идеальном случае эластичность должна быть близка к 1, для наиболее характерных значений X, например при X = 0,1, эластичность всего лишь 0,037742.

Для вычисления эластичности необходимо определить

йР(Л) = ^ йа1 Л) . (3)

йЛ ¡=1 йЛ

Уравнения, необходимые для определения йai (Л) из (3), могут быть получены

йЛ

дифференцированием системы, определяющей стационарное распределение

^=ёИЛ М (Л)+ат (Л) МЛ. йЛ йЛ йЛ

1 ( 2 )

При нахождении а1 ( ' из (3) одно из уравнений данной системы нужно за-

¿Л

менить на у ¿а1(Л) = о. Отметим, что р'0=-р0, р = р0 - р1 и так далее, то есть £ ¿Л

Рк = Рк-1 - Рк ■ Поэтому матрица (Л) будет иметь вид (табл. 5).

¿Л

Таблица 5

Расчетная матрица

1 / ] 7 6 5 4 з 2 1

7 Ро - Ро о о о о о

6 Ро 0 - Ро о о о о

5 Р1 Ро - Р1 о - Ро о о о

4 Р 2 Р1 - Р2 Ро - Р1 о - Ро о о

3 Рз Р2 - Рз Р1 - Р 2 Ро - Р1 о - Ро о

2 Р 4 Рз - Р4 Р2 - Рз Р1 - Р2 Ро - Р1 о - Ро

1 Р5 Р4 - Р5 Рз - Р4 Р2 - Рз Р1 - Р2 Ро - Р1 - Ро

Численный анализ эластичности в диапазоне 0,1 >Л> 0,2 показывает, что она не превышает значения 0,08, в то время как для идеальной системы эластичность должна быть близка к единице.

6. Управление структурой

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

Можно взять другое количество классов, другой входной первоначальный класс, другую премиальную шкалу, установить другие переходные правила или изменить условие критерия события, рассматриваемого в качестве отклонения от норматива, что эквивалентно изменению вероятности его наступления, сделать СБМ более «жесткой», сдвинуть в область больших значений q и, соответственно, Л.

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

больших значений q и, соответственно, Л.

Таким образом, возможности администрации ограничены, прежде всего, математическими свойствами СБМ. Рассмотрим еще один фактор, который оказывает существенное влияние на структуру, — это увольнение сотрудников по тем или иным причинам и прием в связи с этим новых при условии, что общее количество сотрудников должно остаться неименным.

Пусть £ — вероятность выхода из системы для сотрудника из /-го класса (увольнение). Поскольку в течение периода сотрудник может из /-го класса пек

рейти в у-й класс или уволиться, то Е Ру + / = 1. Если, как и раньше, wi(k) —

7 =1

удельный вес /-го класса в момент / = к, то доля уволившихся по всем классам:

к

Е Щ (к)/г .

;=1

Прием на работу производится в соответствии с определенными требованиями, пусть г — это доля от общего числа принимаемых на работу зарезервирован-

к

ных для /-го класса, при этом Е Г = 1, и, следовательно, полная доля прини-

.=1 к

маемых в /-й класс составит: ГуЕ /кЩк .

к =1

Обозначим новые переходные вероятности с учетом увольнений и приема на работу как рГ] (Л), нетрудно видеть, что р..(Л) = (1 — /)р^ (Л). Действительно,

7 7

так как Е Рц (Л) = 1, то Е р(Л) = 1 — /1.

7=1 7=1

Вероятности распределения персонала по классам с учетом увольнения и

к к

приема примут вид: щ (к +1) = Е Щ (к)Р^ + г) Е Щ (к)/ или в векторной фор-

¡=1 г =1

ме: м>(к + 1) = мг(к)[Ы(Х) + / Т г], где / г — соответствующие векторы.

Выберем в качестве примера следующие значения / г (табл. 6).

Таблица 6

Распределение персонала по классам

га) о,6 о, 15 о,о5 о,о5 о,о5 о,о5 о,о5

Ю) о,о5 о,о5 о,о5 о,о5 о,о5 о,о5 о,о5

Другими словами, предполагается равномерное увольнение по тем или иным причинам в размере 5 % от численности данного класса и прием в размере 60 % — в первый класс, в размере 15 % — во второй класс и в остальные классы по 5 % от общего количества принимаемых на работу.

Предельное распределение, получаемое как собственный вектор уравнения:

аТ = аТ[М(Л) + /тг] при Я = 0,1, приведено в табл. 7.

Таблица 7

Предельное распределение

3 7 6 5 4 3 2 1

Щ, % 3,76 4,38 4,44 4,60 5,41 10,69 66,72

Для сравнения в табл. 8 приведено предельное распределение без учета увольнений и приема.

Таблица 8

Предельное распределение (без учета увольнений и приема)

j 7 6 5 4 3 2 1

Wj , % 0,00 0,00 0,00 0,00 0,015 0,995 98,99

Таким образом, налицо изменение структуры. Библиографический список

1. Бартоломью Д. Дж. Управление структурой преподавательского состава в университете // Математическое моделирование. М.: Мир, 1975.

2. Лемер Ж. Системы бонус-малус в автомобильном страховании. М.: Янус-К, 2003.

E.N. Barysheva, L.A. Saraev*

MATHEMATICAL MODELLING OF STRUCTURAL CHANGES OF THE PERSONNEL OF THE ENTERPRISES

In the article the method of governing the structure and quality of personnel on the industrial enterprise in the conditions of market relations is described. Dynamics of structure of personnel and opportunity of management are considered on the basis of system of bonus-malus with a view of forecasting and increasing of efficiency.

Key words: system of bonus-malus, production function, quality of the personnel, class, scale of charges, stationary distribution, the limiting distribution.

* Barysheva Evgeniya Nikolaevna ([email protected]), Saraev Leonid Alexandrovich ([email protected]), the Dept. of Mathematics and Business-Informatics, Samara State University, Samara, 443011, Russia.