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

УДК 519.612.4: 519.865: 331.554

01.00.00 Физико-математические науки

ЧИСЛЕННЫЙ АЛГОРИТМ В ЗАДАЧЕ САМООРГАНИЗАЦИИ ТРУДОВЫХ РЕСУРСОВ

Невечеря Артём Павлович аспирант

Кубанский государственный университет, Краснодар, Россия

Предложен численный алгоритм решения задачи самоорганизации трудовых ресурсов. Под этой задачей подразумевается задача отыскания вероятностей устройства и увольнения специалистов в определённых отраслях внутри исследуемого рынка труда. При этом используется математическая модель динамики трудовых ресурсов. Полученная задача является некорректной, так как уравнений в системе, описывающей данную задачу, меньше, чем неизвестных. Разработанный с учётом предметной области алгоритм, позволяет гарантированно найти нормальное решение за конечное число итераций. Сам алгоритм поделён на два ключевых этапа. Изначально модифицированным методом Гаусса для недоопределённых систем находится безусловное нормальное решение поставленной задачи. В дальнейшем найденное решение проецируется в подпространство допустимых значений, после чего с помощью метода проекции градиентов находится нормальное решение задачи с учётом условий неотрицательности искомых значений. Предложенный алгоритм был успешно использован при разработке в среде программирования С++ приложения, ориентированного на решение задачи самоорганизации трудовых ресурсов. Сравнительный анализ быстродействия данного приложения и надстройки «Поиск решений» табличного процессора MS Excel показал, что одна и та же задача в реализованном автором приложении решается значительно быстрее, чем в табличном процессоре MS Excel при использовании указанной надстройки, что свидетельствует о высокой эффективности предложенного в данной работе алгоритма

Ключевые слова: РЫНОК ТРУДА, САМООРГАНИЗАЦИЯ, ДИНАМИКА ТРУДОВЫХ РЕСУРСОВ, ЧИСЛЕННЫЙ АЛГОРИТМ, МЕТОД ПРОЕКЦИИ ГРАДИЕНТА

UDC 519.612.4: 519.865: 331.554

Physical and mathematical sciences

A NUMERICAL ALGORITHM IN THE PROBLEM OF SELF-ORGANIZATION OF LABOR RESOURCES

Nevecherja Artjom Pavlovich postgraduate student

Kuban State University, Krasnodar, Russia

In this article, there is a numerical method of solving the problem of self-organization of the labor resources. The problem deals with finding probabilities of hiring and the layoffs of specialists from the sectors of the labor market. A mathematical model of labor resources dynamics is used to solve this problem. The initial problem is incorrect, because number of equations of the descriptive system is less than number of unknown variables. A special algorithm is designed for guaranteed finding the normal solution in finite number of iterations. The algorithm is separated into two key stages. Initially, unconditional normal solution of the problem is found by applying the modified method of Gauss for underdetermined systems. Later, this solution is projected in the subspace of permissible values. After that, the normal solution of the problem with consideration of non-negativity of the desired values is being found by using the gradient projection method. The proposed algorithm has been successfully used to develop application in programming environment C++. This application is focused on solving of the problem of self-organization of the labor resources. Comparative analysis of speed of the application and add-ins MS Excel " Solver" showed that the same problem is solved much faster in the application designed by the author than in a table processor MS Excel when using the add-in " Solver". This demonstrates the high efficiency of the proposed method

Keywords: LABOR MARKET, SELF-ORGANIZATION, DYNAMICS OF LABOR RESOURCES, NUMERICAL ALGORITHM, GRADIENT PROJECTION METHOD

Введение

Задача самоорганизации трудовых ресурсов сводится к отысканию вероятностей устройства и увольнения работников на исследуемом рынке труда. На данный момент уже существует несколько подходов к решению этой задачи [1-3]. Ранее [4] автором была предложена математическая модель самоорганизации трудовых ресурсов, с помощью которой на основе статистических данных о количестве занятых и безработных, предоставляемых Федеральной службой государственной статистики [5], можно решать одну из актуальных задач современной межотраслевой экономики - задачу прогнозирования динамики трудовых ресурсов. Для решения этой задачи необходимы значения вероятностей динамики трудовых ресурсов за исследуемый промежуток времени. Для подсчёта вероятностей использовалась надстройка «Поиск решений» табличного процессора MS Excel. Получаемые результаты вычислений имели минимальную невязку, тем не менее, допустимая размерность задачи являлась ограниченной (не более 10-ти отраслей), время вычислений превосходило 30 минут.

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

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

1. Математическая модель самоорганизации трудовых ресурсов

Балансовая математическая модель самоорганизации трудовых ресурсов выглядит следующим образом [4]:

N(l) (t +1) = N"' (t) + X N2j) (t) • pj,г) (t) + AN20) (t) + N20) (t)

j=1

• Pl{0'l)( t )■

- Nг) (г) • Г Р2(г) (г) + Р^,"+1) (г) ], г = М, (1)

N2') (г +1) = n2г) (г) + ^) (г) • Р^ (г) - ^) (г) • £ р^3) (г), г = 1, п, (2)

з=1

, п+1

N20)(г +1) = N20)(г) + Д^0)(г)- Д^0)(г) + N20)(г) • £Р(0,3)(г), (3)

3=1

0 < Р2(г)(г)< 1, г = 1, п, (4)

0 < Р3(г,п+1)(г)< 1, г = 1, п, (5)

р1('3)(г )> 0, г = 0~П, з = 1, п +1. (6)

п+1 _

£р;(г,3)(г)< 1, г = 0,п. (7)

3 =1

Здесь г - номер года, п - количество отраслей на исследуемом рынке труда. Количественные величины: ^(г)(г) - количество занятых в г-ой

отрасли; ^г)(г) - количество безработных, последнее место работы

которых было в г-ой отрасли; ^0)( г) - количество безработных, которые

ранее не имели занятости на исследуемом рынке труда, г = 1, п; дN20)( г) -

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

оценить [2]. Вероятностные величины: Р1(3,г)(г) - вероятность того, что

безработный, последнее место работы которого было в 3-ой отрасли,

найдёт работу в г-ой отрасли; Р2()(г) - вероятность того, что специалист,

работающий в г-ой отрасли, будет уволен; Р(0,г)(г) - вероятность того, что

безработный, не имевший занятости на исследуемом рынке труда с момента последнего появления на данном рынке, найдёт работу в г-ой

отрасли;

Р(г,п+1)( г)

- вероятность того, что безработный, последнее место

работы которого было в г-ой отрасли, покинет рынок труда; ^(0,й+1)( г) -

вероятность того, что безработный, ранее нигде не занятый на

исследуемом рынке труда, покинет данный рынок; Р3(г ,и+1)( г) - вероятность

того, что специалист, работающий в момент времени г в г-ой отрасли, покинет рынок труда.

Представим математическую модель (1) - (7) в векторном виде. Для

этого введём дополнительные понятия: р(г ,г+2)( г) - вероятность того, что

безработный, последнее место работы которого было в г-ой отрасли, за

период времени (г,г +1) останется безработным; Р3(г,г)(г) - вероятность

того, что специалист, занятый в г-ой отрасли в момент времени г, в момент времени г +1 также будет занят в г-ой отрасли. Тогда справедливо следующее:

N = А • Р, (9)

Р >0. (10)

Здесь N следующий вектор:

N =

г \

N

V 2г+1

Т

(11)

где

N = (д1, д2,..., дг, д2,..., д 2, Д2 )Т

д; = N1')(г +1)-^(0(г), г = 1,г, д2 = N2')(г +1)-№(г), г = й

д02 = N 20 )(г +1)- N 20 )(г)- ДN 20 )(г).

Вектор Р представляется следующим образом:

Р

: ( Р, Р3(1,"+1) ( X ) ,..., Р^ (X ) , Р3(1Д) ( X ) ,..., Р3М

(12)

где

Р - (Р(0Д) (X),..., Р(0п+2) (X),..., Р(йД) (X),..., Р(пп+2) (X), Р2(1) (X),..., Р2

Матрица А определяется видом векторов N и Р:

Г А А1,1 А1,2 А1,3 А Л А1,4

А - А2,1 А2,2 А2,3 А2,4

V А3,1 А3,2 А3,3 А3,4 У

где

-(А

-4,1 А1,

(п+1)х(п+2)

А

(0) 1,1

А

(1)

1,1

А1),

ГN20)(X) + Л^0) (X)

А(0) . А1,1 (2п+1)х(п+1)

(0),

0

N20)(X ) + AN20)( X)

(0)

0 0

) + AN20)( X) N20)(X ) + А^0)( X)

(0)

(0)

(0)

0 0

N20)( X ) + А^0)( X) 0

N20)( X ) + А^0)( X) N20)( X ) + А^0)( X) 0

(0)

(0)

0 0

0 0

(0)

Г у(') Л А1,1

А

(■)

1,1

А

*('■)

1,1

лМ1) VА1 У

■ -1, п, а",1 -

пх( п+2)

N 2[)( X), N2 2)( X) ,..., N2п)( X)) ),

а"(1) А1,1

пх(п+2)

0

- ^)(X) - N21)(X)

0

0

0

0

0

0

0

-N21)(x) -N21)(x) 0 1-ая строка

0

0

А"?) - А

1,1

1х( п+2)

1,(п+2)

А -( А(0) А(1) , А2,1 - А2,1 А2,1

(п+1)х( п+2)

?),

т

т

0

0

0

0

0

0

0 0

42 = 1 1

(п+1)х( п+2)

ч 0 0

А3,1 = вп,(п+2)(п+1), А1,2 пх( п+2)( п+1)

000

111 (I +1) -ая строка

000

I = 0, п

г А \

А1,2

А"

Ат

V А1,2 )

4,2 = dlag (- Г), - *<2)( Г) ,..., - *<п)( г))

А12 = diag (г), м(2\ г) ,..., м(п\г)), а;;; = вм, А2,2 =^(п+1),п, АЗ,2 = Еп,

(п+1)хп пхп

1хп

А1,3 =

г А \

А1,3

А

V а1,З ]

, А1 ,3 А1 ,2 , А1 , 3 = в(п+1),п , А2,3 = в(п+1),п , А3,3 = Еп >

пхп

1,3 1,2 : пхп (п+1)хп

( п +1 )хп

А1,4 = в(2п+1),п, А2,4 = в(п+1),п , А3,4 Еп ■ (2п+1)хп (п+1)хп пхп

Здесь в1к - нулевая матрица с размерностью / х к, Е1 - квадратная

единичная матрица с размерностью /х/, diag(а1,а2,...,а1) - квадратная

диагональная матрица с размерностью / х / и элементами а1, а2,..., а1 на

главной диагонали.

Заметим, что матрица А является матрицей полного ранга. Другими словами, все её строки линейно независимы.

2. Задача нахождения вероятностных параметров модели

Математическая модель (9), (10) позволяет по статистическим данным

в моменты времени г и г +1 получить значения вероятностей JP1(l,])(г),

Р2(к)(г), Р3(/,п+1)(г), I = 0, п, ] = 1, п +1, к, / = 1, п, что, в свою очередь,

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

Система, соответствующая задаче (9), (10) содержит 4п + 2 строки. Количество вероятностных параметров в (9), (10) равняется

(п +1)(п + 2) + 3п - п2 + 6п + 2 - т . Таким образом, количество строк матрицы А при любом п будем меньше количества столбцов данной матрицы. Следовательно, задача отыскания Р^1'])(X), Р2(к^ (X), Р3(/,п+1)(X),

1 - 0, п, ] -1, п +1, к, / -1, п, соответствующих (9), (10) является некорректной. Её нормальным решением будем называть вектор [6]:

Р0 -нш Ра. (14)

н®0

Здесь РН доставляет на множестве М минимум функции ^Н(Р), где

М-{Ре Ят,Р >0,1 - 1т},

Ра(Р)-\\АР - N112 +н-||Р||2 н> 0. (15)

Здесь |. 1 и | . ||2 - евклидовы нормы в пространстве Я4п+2 и Ят

соответственно.

Для решения данной задачи воспользуемся следующим подходом:

1) Найдём безусловное нормальное решение, то есть минимальный по норме вектор Р, удовлетворяющий равенству (9).

2) Проекцию данного безусловного нормального решения на множество М, будем считать начальным приближением при решении задачи (9), (10).

4) Для отыскания нормального решения, удовлетворяющего (9), (10), будем использовать метод проекции градиента.

2.1. Нахождение начального приближения

В задаче (9), (10) содержатся нижние ограничения (10) - области допустимых значений вероятностей построенной математической модели. Найдем начальное приближение задачи (9), (10), которое будет являться ее безусловным нормальным решением. Другими словами, рассмотрим СЛАУ с прямоугольной матрицей полного ранга

Ах = /

где

А = (ау), I = 1, п, у = 1, т , т > п, х е Ят, / е Яп.

Рассмотрим множество 5 = {х е Ят : Ах = /}. Пусть множество £ не пусто, тогда задача

ш1п х хе£ 11 11

имеет единственное решение.

Воспользуемся методом Гаусса, модифицированным для недоопределённых систем.

Матрицу системы (16) приведём к ступенчатому виду, эти же преобразования нужно применять к столбцу свободных членов - вектору N (11). Получим:

1 0 0 1

00

0 а

1,п+1

а

1,т

0 а,

2,п+1

а,

2,т

1 а

п , п +1

а„

/1Л /;

/

п

Обозначим полученную систему - А, правую часть - /', £' = {х е Ят : Ах = /'}. Так как £ = £', то

ш1п|| х| |2 = ш1п|| х||2

хе£ хе£'

(17)

2

/

||хц = Е х;=Е х; + Е х;=Е /- Е

1=1 I=1 I=п+1 I=1 V у=п+1

\

Найдем решение задачи (17). Пусть хе Ы', тогда

т

+ Ех2 = Ф(

г=1 V у=п+1 у

ш1п||х2 = шпФ(Xn+1,K,хт ) .

ауху

х

п+1

,х ),

т

1=п+1

хе£'

хе£'

2

Запишем достаточные условия минимума линейного функционала Ф( хт):

ЭФ

Эх,,

0, к = п +1, т.

Вычислим

ЭФ

Эх,

II

-2Е

к ¿=1

/ - Е

V ;=п+1

а;х;

а

¿к

+ 2 хк = 0.

Для отыскания хк, к = п +1, т, получаем следующую СЛАУ:

т

Е ЕЕ'

Обозначив:

ам^х. + хк = £ /о*, к = п +1, т.

И ¿к ] к ;=п+1 ¿=1 ¿=1

Ък]=Е оо , ьк=Е -^к, к=п+1, m,

¿=1

¿=1

получим

Е V; + хк = Ьк, к = п +1, т.

}=п+1

или

(5 + Е) у = Ь ,

'Ь — Ь л

(п+1),(п+1) (п+1),т

здесь 5

V Ьт,(п+1)

Ь

(Ь л ип+1 (х л п+1

, Ь= ... , У = —

V Ьт у V хт )

Заметим, что

Ь; = Е 00 = (О, ак), К = Е /О* = (/, ак),

¿=1

¿=1

где

/ л а1,к 1 (лл

Ок = VОт,к , , / = V /т )

(18)

п

п

Таким образом, решив СЛАУ (18), найдем минимальное по норме решение задачи (9), (10).

2.2. Градиентный метод решения задачи (9), (10)

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

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

1. Задаём а> 0, 0 <Ь< 1, £> 0; полагаем х(0) равным спроецированному на множество М начальному приближению, найденному при решении (18); вычисляем по (15) /(х(0)), /'(х(0)); полагаем к =1.

2. Проецируем антиградиент (- / (х(к))) в область допустимых

/'( х^)) . Полагаем 1к = а.

значений и вычисляем

г'(х(к-1))

3. Вычисляем Ах(к) = -1 ■^-, х(к) = х(к-1) +Ах(к), / (х(к)).

кГ (х(к-1)) , , /1 )

4. Проверяем условие выбора Як: / (х^к)) < / (х(к 1)).

Если оно выполняется, то осуществляем переход к п. 5 данного алгоритма. Иначе полагаем 1 = (51к и осуществляем переход к п. 3.

5. Полагаем I = 0, х(к) = х{к).

6. Полагаем I = I +1, вычисляем х(к) = х^ +1.

7. Если I > 2, тогда переходим к п. 8, в противном случае переходим

п. 6.

8. Проверяем следующие условия:

/(*й)</(х-2), (19)

/(*£>)</(хк)), (20)

ь(*«)< 0. (21)

Здесь Ь (х) принимает отрицательное значение тогда и только тогда,

когда вектор х содержит хотя бы одну отрицательную компоненту. Если условия (19) - (21) выполняются, то, очевидно, оптимальное решение (для

к-й итерации) будет находиться на интервале (х(-2, х(к)). Опишем алгоритм

нахождения решения на данном интервале:

8.1. Для дальнейшего приближения обозначим: х(-2 как хр^; х^ —

х(к) • х(к) - х{к)

8.2. Вычисляем

Х(кХ(к) Х(к)_ х(к)

х(к) _ now pred х(к) _ next лnow midl 2 ' mid 2 2 '

Если L (xp)d )< 0, тогда переходим к нахождению неотрицательного

(к)

решения методом дихотомии с краевыми точками c — xKpJed, a _ х% -(хШ_ )_ 2xfred _ хЩ. Считая полученный вектор вектор-

решением к-й итерации х^к), переходим к п. 9. Если L (хр}еЛ) > 0 переходим к п. 8.3.

8.3. В полученной последовательности приближений хр}ес1, х^х^,

х^тш2, х^ могут найтись три соседних точки х(_1, х(к), х^, для которых будут выполняться условия:

f (У))< f (хП), f (х,(к))< f (х«).

В данном случае принимаем http://ej .kubagro.ru/2016/04/pdf/86.pdf

Иначе:

х(к) _ х(к) х(к) _ „(к) х(к) _ „(к) ртей Л2-1 ' Апо» _ Аг > Лпех1 ~ г+1 '

х(к) _ х(к) и х(к) _ х(к) ртей лпоп п Лпок тгй 2 •

8.4. Если выполняется условие

f ( Хп1х1 ) f ( Хр}еа )

(24)

тогда вектор-решением к-й итерации х( к) будем считать вектор х^, если

е (хпо! )< 0, и хП01, в противном случае. Если (24) не выполняется, переходим к шагу 8.2.

9. Вычисляем х(к)).

10. Проверяем условия:

1

п .

X«

(к)

п—1

1 п

п~1

XI яр(к)

<е,

<е,

где g(к) - компоненты вектора антиградиента (g(к)_-^ (х^к))), gp(k) -

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

Если они выполняются, то полагаем х* @ х^к), f * @ f (х(к)) и

завершаем вычисления. Иначе полагаем к _ к +1 и осуществляем переход к п. 2.

Решением поставленной задачи будет вектор х*.

По представленному алгоритму нахождения нормального решения задачи (9), (10) было написано приложение в среде программирования С++. Проведённые тесты показали, что при исследовании 10-и отраслевого рынка труда невязка решения при а_ 0,001 (9), полученная в

реализованном приложении на С++, практически нулевая - не превышает 0.0001 (отношение невязки к объёму исследуемого рынка труда), и незначительно превосходит невязку решения, полученного с помощью табличного процессора MS Excel. Однако следует отметить скорость нахождения решения - для реализованного приложения она занимает около 15 секунд, в то время как при использовании надстройки «Поиск решений» (метод обобщённого приведённого градиента) табличного процессора MS Excel решение вычисляется около 30 минут. Также следует отметить, что реализованное приложение способно решать задачи с размерностью, превосходящую максимально допустимую размерность задачи при её решении в табличном процессоре MS Excel.

3. Применение разработанного метода на практике

Воспользуемся данными, предоставляемыми федеральной службой государственной статистики Российской Федерации [5]. В соотношениях (1) - (3) межотраслевой балансовой модели в качестве количественных характеристик используем данные по рынку труда, состоящему из следующих 10-и отраслей экономки РФ: №1 - «Сельское и лесное хозяйство, охота, рыболовство и рыбоводство», №2 - «Добыча полезных ископаемых», №3 - «Обрабатывающие производства», №4 -«Производство и распределение электроэнергии, газа и воды», №5 -«Строительство», №6 - «Оптовая и розничная торговля, ремонт автотранспортных средств, мотоциклов, бытовых изделий и предметов личного пользования, гостиницы и рестораны», №7 - в «Транспорт и связь», №8 - «Финансовая деятельность, операции с недвижимым имуществом, аренда и предоставление услуг», №9 - «Государственное управление и обеспечение военной безопасности, социальное обеспечение», №10 - «Образование». Распределение трудовых ресурсов по

отраслям, участвовавших в производственном процессе в 2011 — 2013 гг., приведено в таблице 1. Таблица 1. Распределение работающих специалистов по отраслям (в тыс.

чел.) в 2011 — 2013 гг.

1 И(1) И(2) И(3) И(4) И(5) И(6) И(7) И(8) И(9) И<10)

2011 5456,0 1417,1 10628,5 2267,4 5101,7 12754,2 6660,5 6164,5 5456,0 6518,8

2012 5222,8 1430,9 10731,8 2361,0 5294,4 13021,3 6725,3 6224,5 5365,9 6596,5

2013 4997,4 1570,6 10565,9 2284,5 5425,8 13100,3 6746,5 6425,2 5247,3 6532,3

Численность безработных, последнее место работы которых было в ¿-ой отрасли, оценим через соотношения численности выбывших работников списочного состава в процентах к среднесписочной численности работников в Российской Федерации по видам экономической деятельности за 2011 — 2013 гг. Приведённые оценочные значения содержаться в таблице 2.

Таблица 2. Оценочное количество безработных, последнее место работы

которых было в ¿-ой отрасли (в тыс. чел.) в 2011 — 2013 гг.

1 И(1) И(2) И(3) И(4) И(5) И(6) И(7) И(8) И(9) И<10)

2011 917,78 238,46 274,54 232,3 449,65 1067,37 265,74 520,05 160,15 174,23

2012 788,02 191,3 227,36 181,74 367,15 886,61 228,09 443,67 121,4 143,48

2013 843,99 205,16 219,76 180,33 378,92 841,07 213,19 443,17 126,31 162,08

Количество безработных на момент времени 1, которые не были заняты каким-либо видом деятельности с момента последнего попадания

на исследуемый 10-и отраслевой рынок труда, (t), в 2011 году

равняется, примерно, 100 тысячам человек, в 2012 году — 113 тысячам (оценка данного фактора возможна за счёт анализа количества выпускников и мигрантов, поступивших к этому периоду и ориентированных на представленные отрасли). В течение 2011 — 2012 гг. на исследуемый рынок труда поступило (без учёта распределения этих безработных между отраслями за данный промежуток времени),

приблизительно, 300 тысяч человек (АИ20)( 2011) = 300).

По приведённым данным за 2011 - 2012 гг. оценим значения вероятностей математической модели (9), (10). Результат приведён в таблицах 3, 4.

Таблица 3. Значения вероятностей устройства специалистов, последнее место работы которых было в ¿-ой отрасли, ву-ую отрасль (округлено до

третьего знака после запятой)

Р(1. у) 11 1 2 3 4 5 6 7 8 9 10 11 12

0 0,089 0,086 0,066 0,084 0,097 0,087 0,068 0,087 0,027 0,056 0,004 0,248

1 0,116 0,109 0,063 0,103 0,134 0,111 0,067 0,110 0,004 0,040 0,000 0,142

2 0,105 0,103 0,092 0,102 0,110 0,104 0,093 0,104 0,068 0,086 0,029 0,005

3 0,089 0,087 0,073 0,085 0,095 0,088 0,075 0,087 0,046 0,067 0,004 0,204

4 0,112 0,110 0,098 0,108 0,116 0,111 0,099 0,110 0,075 0,093 0,038 0,000

5 0,099 0,096 0,073 0,093 0,108 0,097 0,075 0,096 0,029 0,062 0,004 0,169

6 0,045 0,037 0,017 0,030 0,067 0,039 0,013 0,039 0,004 0,012 0,000 0,697

7 0,100 0,098 0,085 0,096 0,105 0,099 0,086 0,098 0,059 0,078 0,015 0,080

8 0,099 0,095 0,070 0,092 0,110 0,097 0,072 0,096 0,018 0,057 0,004 0,190

9 0,112 0,111 0,103 0,110 0,115 0,111 0,103 0,111 0,087 0,099 0,061 0,000

10 0,100 0,098 0,090 0,097 0,103 0,099 0,090 0,099 0,072 0,085 0,044 0,022

В таблице 3 вероятность, находящаяся в ячейке, расположенной на пересечении строки, крайняя левая ячейка которой содержит число 1, и столбца, крайняя верхняя ячейка которого содержит число у, соответствует

вероятности р( ] 'г) (?) - вероятность того, что безработный, последнее место

работы которого было в у-ой отрасли, найдёт работу в ¿-ой отрасли. Уход в 11-ую отрасль соответствует уходу из исследуемого 10-и отраслевого рынка труда. В 12-ом столбце содержаться вероятности того, что безработный специалист, последнее место занятости которого было в ¿-ой отрасли, останется безработным.

В таблице 4 приведены вероятности того, что занятые на исследуемом

рынке труда в 2011 году специалисты в 2012 году либо уволятся (Р2(1)), либо покинут исследуемый рынок труда (Р3(11)).

Таблица 4. Значения вероятностей увольнения работающего специалиста из г-ой отрасли и вероятности его ухода из исследуемого рынка труда

(округлено до четвёртого знака после запятой)

р г

1 2 3 4 5 6 7 8 9 10

2 0,1204 0,1338 0,0160 0,0873 0,0570 0,0111 0,0309 0,0558 0,0258 0,0213

р(ш) Г3 0,0000 0,1400 0,0017 0,0405 0,0000 0,0000 0,0037 0,0005 0,0144 0,0040

р('.') 13 0,8796 0,7262 0,9823 0,8722 0,9430 0,9889 0,9653 0,9437 0,9598 0,9747

Заключение

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

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

Автор благодарит Анастасию Белашову за помощь при реализации алгоритма и проведении тестовых расчётов.

Литература

1. Семенчин Е.А., Невечеря А.П. Об обратной задаче в математической модели самоорганизации рынка труда // Фундаментальные исследования. - М.: Академия Естествознания, 2014. - № 6. - С. 1184 - 1190.

2. Семенчин Е.А., Зайцева И.В. Математическая модель самоорганизации рынка труда для двух отраслей экономики // Экономика и математические методы. - М.: Наука, 2004. - Т. 40. В. 4. - С. 137 - 139.

3. Семенчин Е.А., Зайцева И.В. Математическая модель самоорганизации рынка труда для нескольких отраслей экономики // Экономика и математические методы. -М.: Наука, 2007. - Т. 43. В. 1. - С. 133 - 136.

4. Невечеря А.П. Прогнозирование динамики трудовых ресурсов с помощью межотраслевой математической модели / А.П. Невечеря // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета (Научный журнал КубГАУ) [Электронный ресурс]. - Краснодар: КубГАУ, 2015. -

№05 (109). С. 560 - 572. - IDA [article ID]: 1091505033. - Режим доступа: http://ej .kubagro.ru/2015/05/pdf/33.pdf

5. Федеральная служба государственной статистики. [Электронный ресурс]. -Режим доступа: http://www.gks.ru

6. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. - М: Наука, 1979. - 285 с.

7. Васильев Ф.П. Численные методы решения экстремальных задач. - М.: Наука, 1988. - 522 с.

References

1. Semenchin E.A., Nevecherja A.P. Ob obratnoj zadache v matematicheskoj modeli samoorganizacii rynka truda // Fundamental'nye issledovanija. - M.: Akademija Estestvoznanija, 2014. - № 6. - S. 1184 - 1190.

2. Semenchin E.A., Zajceva I.V. Matematicheskaja model' samoorganizacii rynka truda dlja dvuh otraslej jekonomiki // Jekonomika i matematicheskie metody. - M.: Nauka, 2004. -T. 40. V. 4. - S. 137 - 139.

3. Semenchin E.A., Zajceva I.V. Matematicheskaja model' samoorganizacii rynka truda dlja neskol'kih otraslej jekonomiki // Jekonomika i matematicheskie metody. - M.: Nauka, 2007. - T. 43. V. 1. - S. 133 - 136.

4. Nevecherja A.P. Prognozirovanie dinamiki trudovyh resursov s pomoshh'ju mezhotraslevoj matematicheskoj modeli / A.P. Nevecherja // Politematicheskij setevoj jelektronnyj nauchnyj zhurnal Kubanskogo gosudarstvennogo agrarnogo universiteta (Nauchnyj zhurnal KubGAU) [Jelektronnyj resurs]. - Krasnodar: KubGAU, 2015. -№05 (109). S. 560 - 572. - IDA [article ID]: 1091505033. - Rezhim dostupa: http://ej .kubagro.ru/2015/05/pdf/33.pdf

5. Federal'naja sluzhba gosudarstvennoj statistiki. [Jelektronnyj resurs]. - Rezhim dostupa: http://www.gks.ru

6. Tihonov A.N., Arsenin V.Ja. Metody reshenija nekorrektnyh zadach. - M: Nauka, 1979. - 285 s.

7. Vasil'ev F.P. Chislennye metody reshenija jekstremal'nyh zadach. - M.: Nauka, 1988. - 522 s.