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

УДК 69.05 DOI: 10.22227/1997-0935.2018.1.80-86

ПРИМЕНЕНИЕ МЕтОДА лИНЕЙНОГО

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

БРИГАД

Е.В. Тарарушкин

Российский университет транспорта (МИИТ), 127994, г. Москва, ул. Образцова, д. 9, стр. 9

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

Цели: назначение составов строительных бригад из условия максимально возможного качества составов, которое определяется квалификацией, опытом рабочего и необходимостью назначения рабочего для определенного вида работ. Задача поиска оптимальных составов бригад поставлена в виде задачи линейного программирования. Методы: для решения поставленной задачи использовался симплекс-метод совместно с методом ветвей и границ для получения целочисленных решений. Для реализации методов решения использовался язык программирования математического пакета МаИаЬ.

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

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

КЛЮчЕВыЕ СЛОВА: качество строительных работ, рабочий, бригада, захватка, задача линейного программирования, симплекс-метод, метод ветвей и границ

Благодарности: выражаю благодарность Фетисову Д.А. (к.ф.-м.н., доцент кафедры математического моделирования МГТУ им. Н.Э. Баумана) за помощь при написании выпускной квалификационной работы на тему «Решение задачи о назначении экипажей воздушных судов» по специальности «Прикладная математика». Также выражаю благодарность независимым рецензентам за проверку текста статьи.

ДЛЯ ЦИТИРОВАНИЯ: Тарарушкин Е.В. Применение метода линейного программирования для назначения оптимального состава строительных бригад // Вестник МГСУ. 2018. Т. 13. Вып. 1 (112). С. 80-86.

<N

О >

APPLICATION OF THE LINEAR PROGRAMMING METHOD FOR CHOOSING OPTIMAL STRUCTURE OF CONSTRUCTION TEAM

E.V. Tararushkin

Russian University of Transport (MUT), 9, bld. 9, Obraztsov st., Moscow, 127994, Russian Federation

Subject: creation of a mathematical model based on the theory of operations research for appointing the members of construction teams when executing construction, installation and finishing works in civil engineering.

Research objectives: optimization of the construction team structure from the condition of the highest possible quality of teams, which is determined by qualification, experience of the worker and the need to appoint a worker for a certain type of IQ works. The problem of finding the optimal composition of the team is formulated as a linear programming problem.

ft Materials and methods: to solve the given problem, a simplex algorithm was used in combination with the branch and bound

T" method to obtain an integer value solution. The programming language of the mathematical package Matlab was used to

«g implement the indicated methods of solution.

O Results: a mathematical model was created for appointing workers of construction teams. An example of realization of the

H obtained mathematical model is given.

^ Conclusions: this mathematical model can be used for practical purposes for planning and organization of the construction

O process. The use of the model by designers will help to reduce their time expenditures, reduce the labor intensity of work and

^ better perform the assignment of workers to construction teams, especially when designing large-scale construction projects.

^ KEY WORDS: quality of construction works, worker, construction team, construction staging, linear programming problem,

5 simplex algorithm, branch and bound method

H FOR CITATION: Tararushkin E.V. Primenenie metoda lineynogo programmirovaniya dlya naznacheniya optimal'nogo sos-

tava stroitel'nykh brigad [Application of the linear programming method for choosing optimal structure of construction team].

^ Vestnik MGSU [Proceedings of the Moscow State University of Civil Engineering]. 2018, vol. 13, issue 1 (112), pp. 80-86.

80 © Е.В. Тарарушкин

Применение метода линейного программирования для назначения

- С. 80—86

оптимального состава строительных бригад

ВВЕДЕНИЕ

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

ОБЗОР ЛИТЕРАТУРЫ

Задача о назначениях является частным случаем задачи линейного программирования теории исследований операций [1, 2]. Некоторые виды задач линейного программирования применяются при планировании и организации строительного производства [3, 4]: транспортная задача, задача оптимального размещения строительной базы, распределительная задача и др.

Задачи о назначении оптимальных составов рабочих групп является модифицированной задачей о назначениях. Известны задачи о назначении составов рабочих групп в таких отраслях, как железнодорожные и авиаперевозки [5-8].

МЕТОДЫ

Составим математическую задачу поиска оптимального состава строительных бригад для работ на захватках строительного объекта с учетом временных интервалов выполнения работ. Пусть необходимо составить бригады на строительно-монтажные работы на захватках с определенным количеством рабочих с учетом временных интервалов выполнения работ, так чтобы качество работ на захватках объекта было максимально возможным. Работы бригадами на захватках выполняются последовательно. Автор рассмотрел задачу назначение состава рабочих звеньев методом линейного программирования, но без учета временного интервала исполнения работ [9].

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

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

2. Имеется п. рабочих г-го типа, г = 1, N, так что каждый конкретный рабочий определяется парой (г, /), где г — тип рабочего;/ — порядковый номер рабочего среди рабочих г-го типа, / = 1, п1.

3. Квалификация рабочего (г,/) д (сюда входят опыт работы, разряд рабочего и др.) может принимать значение из отрезка [0,1].

4. Учет квалификации должен зависеть от степени необходимости рабочего для определенного вида работ, поэтому вводится вес а, с которым учитывается степень необходимости -го типа рабочего (а может принимать только натуральные значения), так что квалификация рабочего 0,/), которой оперирует проектировщик, вычисляется по формуле ад., г = 1, N, / = 1, п.

5. В соответствии с количественным составом назначаемой бригады и временем выполнения работ строительный объект разбит на М захваток; каждая захватка однозначно определяется парой чисел (к, I), где k — порядковый номер захватки; I — порядковый номер временного интервала для выполнения строительно-монтажных работ (месяц, неделя и др.); k = 1, М, ! = 1, тк, тк — количество временных интервалов необходимых для завершения работ на захватке с номером к.

6. На захватку под номером к необходимо на-

. к

значить о. сотрудников г -го типа; двоичная пере-

к!

менная X/ принимает значение единицы, если сотрудник (г , /) назначен на захватку (к, I), и значение нуля в противном случае; таким образом, в постановке задачи присутствует ограничение

V хи = Ь, i = 1, N, к = 1,М, ! = 1,

т,г

/=1

7. Некоторые рабочие могут быть закреплены за каким-либо видом работ на определенный временной интервал работ на захватке по желанию проектировщика; некоторым рабочим может быть запрещено работать на каком-либо виде работ на определенном временном интервале по желанию проектировщика, т.е. в постановке задачи возможны ограничения вида хкк = 1 для некоторых наборов г , /, к, I и ограничение вида хкк = 0 для некоторых наборов г,/, к, I.

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

= 2

х/ + х/

для некоторых значений к, I и некоторых пар (г,

Оу Л).

9. Должна быть предусмотрена возможность учета пожеланий рабочих при работе на захватках,

00

Ф

0 т

1

*

О У

Т

0

г

1

В

г

3

у

о *

К)

N

О >

с во

2 О

н *

О

X 5 I н о ф ю

т.е. для некоторых рабочих (г, /) допускается наличие ограничений вида

м т

X м < < ^

хх<2 х« < а,

»,2 '

м т ^ ^ м

ХХС х» < А

4,3'

»,1'

(!".,, 4', е{-1, 0, 1} д, „, А,.,.

Абсолютное качество работ на захватке (к, I) определяется выражением

_ХХаАи

1=1 4=1

Максимальное качество работ на захватке к-го номера, которое возможно достичь, задается формулой

Ок _ХаЬк.

_ щах / , г г

,=1

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

0ге1 о» .

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

о:

г (х) _ X 0

(к,') 0

к'

щах

щах.

Поставленная выше задача оптимального назначения рабочих в бригады в описанной постановке является задачей линейного программирования. Для решения данной задачи используется симплекс-метод в совокупности с методом ветвей и границ для получения целочисленных решений [10-12].

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

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

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

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

Пусть необходимо выполнить строительно-монтажные работы на строительном объекте. Объект состоит из двух захваток (М = 2). Часть работ одной бригадой на одной захватке выполняется за один месяц. Всего на окончание работ на одной захватке отводится 3 мес. (тк = 3). Для выполнения определенных работ на захватках строительного объекта у строительной организации имеется три типа профессий рабочих, например монтажник, бетонщик и арматурщик (Ы = 3). Для первого типа пять рабочих (п = 5), для второго — восемь рабочих (п = 8), для третьего — шесть рабочих (п3 = 6). Необходимо назначить составы рабочих бригад для работ на захватках. Назначение состава бригад необходимо выполнить с учетом максимизации качества работ. Все работы выполняются последовательно.

Векторы квалификации для рабочих следующие:

^ = [0,8 0,7 0,6 0,5 0,4] — для рабочих первого типа;

д2 . = [0,8 0,3 0,4 0,3 0,7 0,9 0,8 0,4] — для рабочих второго типа;

д3 . = [0,6 0,8 0,4 0,6 0,7 0,3] — для рабочих третьего типа.

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

4,

0,8 0,7 0,6 0,5 0,4 0,0 0,0 0,8 0,8 0,3 0,4 0,3 0,7 0,9 0,8 0,4 0,6 0,8 0,4 0,6 0,7 0,3 0,0 0,0

Весовой вектор квалификаций для рабочих разных типов: а = [2 1 1]

Вектор ограничений Ь выглядит следующим образом:

"2 3 4" 2 4 3

Ь _

Целевая функция

Применение метода линейного программирования для назначения

- С. 80—86

оптимального состава строительных бригад

I (х)=—х;; +— х1121 +— х;з+... +— х? +

55

55

..23

55

55

-х" +--х36 ^ тах.

110 35 110 36

Система ограничений из п. 6. постановки задачи:

х11 =1

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

Система ограничений из п. 8 постановки задачи имеет вид

< + х12 + < + х14 + х15 = 2

х11+х21 + х21 + х14 + х15 + х11 + х27 + х11 = 3 х» + хи + хи + х» + хи + хи = 4

31 32 33 34 35 36

х13 + х12 + Х3 + х14 + = 2 х21 + х13 + х23 + + х25 + х16 + х23 + х23 = 3

х13 + х13 + х13 + х13 + х13 + х13 = 4

31 32 33 34 35 36

х21 + х2 + х21 + X,2; + х21 = 2

х 2' + х + х 2' + х + х + х + х 21 + х 21 = 4 21 22 23 24 25 26 27 А28 ^ х 21 + х 21 + х 21 + х21 + х 21 + х 21 = 3

31 32 33 34 35 36

х" 21 +х23 = 2

х12 23 +х25 = 2

х 21 22 +х241 = 2

х 23 35 +х263 = 2

х23 + х2223 + х23 + х243 + х23 = 2

х23 + х23 + х23 + х23 + х 23 + х 23 + х 23 + х 23 = 4

21 22 23 24 25 26 27 28 х 23 + х 23 + х23 + х23 + х23 + х23 = 3

31 32 33 34 35 36

Пояснение ограничений в системе на примере следующего ограничения:

х» + х12 + хЦ + хЦ + хЦ = 2

Пояснение ограничений в системе на примере следующего ограничения:

х21 + х23 = 2

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

Система ограничений из п. 8 постановки задачи имеет вид

х'2 + х'2 + х^ < 1 х12 - х122 < 0

х141 + х142 + х14 < 2

Х4 - х^3 < 0

х1131 + Х32 + Х33 + х21 + х22 + х23 < 2

х? - х22 < 0.

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

Система ограничений из п. 7 постановки задачи имеет вид

х11 11 =0

х11 13 =0

х11 32 =1

х12 31 =0

х21 12 =0

х21 35 =0

х22 21 =1

х22 35 =1

х23 12 =0

х23 28 = 1

х„ = 0

Пояснение ограничений в системе на примере следующих ограничений:

Х2 + х2 + Х23 < 1

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

х12 - х122 < 0

(второй рабочий из первого типа не может быть одновременно назначен на первый месяц работы на первой захватке и на второй месяц на первой захватке);

х14 + х14 + х14 < 2

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

Требование целочисленности имеет вид

Пояснение ограничений в системе на примере следующих ограничений:

хк е:

(первый рабочий из первого типа не может быть назначен на первую захватку на первый месяц работ);

Е {0,1} , при

к = 12, I = 13, i = 1Д у = 1,5, у2 = 18, у3 = 1,6.

На языке программирования математического пакета Matlab [13, 14] была написана реализации симплекс-метода и метода ветвей и границ для решения задачи о назначении оптимального состава

00

Ф

0 т

1

*

О У

Т

0

1

В

г 3

у

о *

К)

Дерево решения примера

строительных бригад в постановке задачи линейного программирования. Дерево решения задачи представлено на рис.

На итерации под номером 1 все значения переменных оказались целочисленными (нуль или единица), кроме следующих переменных:

_ 0,5, Х1 _ 0,5, х;22 _ 0,5, х;32 _ 0,5,

х2 _ 0,25, х2 _ 0,75, х^2 _ 0,75, х^2 _ 0,25,

сч

Значение целевой функции составило / = = 881/220.

Переход на итерацию под номером 2 осущест-^ влялся при х^1 _ 0. На этой итерации все значения ^ переменных оказались целочисленными (нуль или 2 единица), кроме следующих переменных: 10

х23 _ 0,5, х2 _ 0,5, х22 _ 0,5, хЦ _ 0,5,

Значение целевой функции составило / = I- = 439/110.

^ Переход на итерацию под номером 3 осуществлялся при х1'2 _ 1. На этой итерации множество 2 всех решений оказалось пустым. X Переход на итерацию под номером 4 осуществлялся при х1231 _ 0. На этой итерации все значения ¡^ переменных оказались целочисленными (нуль или Ф единица). Значение целевой функции составило/= М = 219/55.

Переход на итерацию под номером 5 осуществлялся при х2 = 1. На этой итерации все значения переменных оказались целочисленными (нуль или единица). Значение целевой функции составило/= = 219/55.

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

выводы

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

Данная постановка задачи может найти практическое применение при планировании и организации строительного производства, при этом вышеприведенная постановка задачи о формиро-

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

С. 80-86

Составы строительных бригад на захватки

Номер захватки

Номер типа рабочего 1 2

Номер месяца Номер месяца

1 2 3 1 2 3

i ii iii iv v vi vii

1 4 1 1 3 1 1

1 5 2 4 4 3 4

2 1 3 1 1 1 1

2 3 5 6 2 5 6

2 6 6 7 4 6 7

2 — — — 6 7 8

3 1 2 1 1 1 2

3 2 3 2 2 2 5

3 4 4 4 4 5 6

3 5 5 5 — — —

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

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

захватке (к, 12) и рабочего (г,/2) на захватке (к, 12), т.е. возможны ограничения вида

х/ + х/2 < 1

для некоторых захваток (к, /1), (к, /2) и некоторых рабочих (г, .Д), (г, /).

ЛИТЕРАТУРА

1. Hiller F.S., Lieberman G.J. Introduction to operations research. 9th ed. McGraw Hill, 2013. 1102 p.

2. Волков И.К., Загоруйко Е.А. Исследование операций / под ред. B.C. Зарубина, А.П. Крищен-ко. М. : Изд-во МГТУ им. Н.Э. Баумана, 2000. 436 с.

3. Сырцова Е.Д. Математические методы в планировании и управлении строительным производством. М. : Высш. шк., 1972. 336 с.

4. Кадыров А.С., Бестембек Е.С., Сунгатол-лакызы А. и др. Решение транспортной задачи на примере строительства протяженных объектов // Международный журнал прикладных и фундаментальных исследований. 2015. № 6. С. 224-227.

5. Caprara A., Fischetti M., Toth P. et al. Algorithms for railway crew management // Mathematical Programming. October 1997. Vol. 79. Iss. 1. Pp. 125-141.

6. Suyabatmaz A., Sahin G. Railway crew capacity planning problem with connectivity of schedules // Transportation research. Part E. 2015. Vol. 84. Pp. 88-100.

7. Klabjan D., Johnson E., Nemhauser G. et al. Solving large airline crew scheduling problems: random pairing generation and strong branching // Compu-

tational Optimization and Applications. October 2001. Vol. 20. Iss. 1. Pp. 73-91.

8. Hoffman K., Padberg M. Solving airline crew scheduling problems by branch-and-cut // Management Science. June 1993. Vol. 39. Issue 6. Pp. 657-682.

9. Тарарушкин Е.В. Назначение состава рабочих звеньев методом линейного программирования // до Вестник БГТУ им. В.Г. Шухова. 2017. № 7. С. 71-74. с

10. Eiselt H.A., Sandblom C.-L. Linear н Programming and its Applications. Springer-Verlag s Berlin Heidelberg, 2007. 380 p. *

11. Вентцель Е.С. Исследование операций: за- Г дачи, принципы, методология. 2-е изд., стер. М. : q Наука, 1998. 208 с. У

12. Конюховский П.В. Математические методы О исследований операций в экономике. СПб : Питер, 2 2000. 208 с. (Краткий курс).

13. Поршнев С.В. Matlab 7. Основы работы и Я программирования. М. : ООО «Бином-Пресс», 2006. ы 320 с.

14. Кетков Ю.Л., Кетков А.Ю., Шульц М.М. С Matlab 6.x.: программирование численных методов. Я СПб. : БХВ-Петербург, 2004. 672 с.

Поступила в редакцию 2 октября 2017 г. Принята в доработанном виде 10 октября 2017 г. Одобрена для публикации 12 декабря 2017 г.

Об авторе: тарарушкин Евгений Викторович — ассистент кафедры строительных материалов и технологий, Российский университет транспорта (МИИт), 127994, г. Москва, ул. Образцова, д. 9, стр. 9; tarщstu@ yandex.ru.

REFERENCE

<N

О >

с

10

«

s о

H >

о

1. Hiller F.S., Lieberman G.J. Introduction to operations research. 9th ed. McGraw Hill, 2013. 1102 p.

2. Volkov I.K., Zagoruyko E.A. Issledovanie operatsiy [Operations Investigation]. Moscow, Moscow State Bauman Technical University, 2000. 436 p.

3. Syrtsova E.D. Matematicheskie metody v planirovanii i upravlenii stroitel 'nym proizvodstvom [Mathematical methods in planning and management of construction production]. Moscow, Vysshaya shkola Publ., 1972. 336 p.

4. Kadyrov A.S., Bestembek E.S., Sungatol-lakyzy A. et al. Reshenie transportnoy zadachi na pri-mere stroitel'stva protyazhennykh ob"ektov [Solution of the transport problem by the example of building long objects]. Mezhdunarodnyy zhurnalprikladnykh i fundamental'nykh issledovaniy [International Journal of Applied and Fundamental Research]. 2015, no. 6, pp. 224-227.

5. Caprara A., Fischetti M., Toth P. et al. Algorithms for railway crew management. Mathematical Programming. October 1997, vol. 79, iss. 1, pp. 125-141.

6. Suyabatmaz A., Sahin G. Railway crew capacity planning problem with connectivity of schedules. Transportation research. Part E. 2015, vol. 84, pp. 88-100.

7. Klabjan D., Johnson E., Nemhauser G. et al. Solving large airline crew scheduling problems: random pairing generation and strong branching. Computational Optimization and Applications. October 2001, vol. 20, iss. 1. pp. 73-91.

Received October 2, 2017.

Adopted in revised form October 10, 2017.

Approved for publication on December 12, 2017.

8. Hoffman K., Padberg M. Solving airline crew scheduling problems by branch-and-cut. Management Science. June 1993, vol. 39. issue 6, pp. 657-682.

9. Tararushkin E.V. Naznachenie sostava rabo-chikh zven'ev metodom lineynogo programmirovaniya [Assignment of the composition of working links by the method of linear programming]. Vestnik BGTU im. V.G. Shukhova [Bulletin of Belgorod State Technological University named after. V.G. Shukhov]. 2017, no. 7, pp. 71-74.

10. Eiselt H.A., Sandblom C.-L. Linear programming and its applications. Springer-Verlag Berlin Heidelberg, 2007. 380 p.

11. Venttsel' E.S. Issledovanie operatsiy: zadachi, printsipy, metodologiya [Research of operations: tasks, principles, methodology]. 2nd ed. Moscow Nauka Publ., 1998. 208 p.

12. Konyukhovskiy P.V. Matematicheskie metody issledovaniy operatsiy v ekonomike [Mathematical methods of research of operations in the economy]. Saint-Petersburg, Piter Publ., 2000. 208 p.

13. Porshnev S.V. Matlab 7. Osnovy raboty i programmirovaniya [Fundamentals of work and programming]. Moscow, Binom-Press Publ., 2006. 320 p.

14. Ketkov Yu.L., Ketkov A.Yu., Shul'ts M.M. Matlab 6.x.: programmirovanie chislennykh metodov [Matlab 6.x .: programming of numerical methods]. Saint-Petersburg, BKhV-Peterburg Publ., 2004. 672 p.

About the author: Tararushkin Evgeniy Victorovich — Assistant, Building Materials and Technologies Department, Russian University of Transport (MIIT), 9, bldg 9, Obraztsova str., Moscow, 127994, Russian Federation, tarmstu@yandex.ru.

S

I

h

О

Ф

10