УДК 519.237.5
И. М. Якимов, А. П. Кирпичников, О. В. Маркина
МОДЕЛИРОВАНИЕ ДВИЖЕНИЯ АВТОТРАНСПОРТА
ПО ЗАГОРОДНЫМ ПРОТЯЖЁННЫМ ТРАССАМ
Ключевые слова: автотранспорт, легковой автомобиль, грузовик, грузоподъёмность, габариты, имитационное моделирование, планирование экспериментов, корреляционный анализ, регрессионный анализ, система имитационного моделирования
GPSSW.
Приведено описание результатов имитационного моделирования (ИМ) движения автотранспорта по загородным трассам большой протяжённости. По результатам ИМ построена математическая модель, состоящая из совокупности уравнений регрессии, связывающих результативные показатели движения автотранспорта с влияющими на них факторами. Оценена степень влияния факторов на результативные показатели движения по математической модели и коэффициентам линейной корреляции.
Keywords: motor transport, motor-car, auto truck, capacity, external dimensions, simulation modeling, experiment design, correlation
analysis, regression analysis, system simulation GPSS W.
Describes the results of simulation modeling (SM) of traffic on suburban routes a large extent. In accordance with any findings of the SM mathematical model was constructed. The model consists of a set of regression equations, linking indicators of traffic with the factors, influencing them. Describes the assessment of the degree of the factor's influence on indicators of effective movement on the mathematical model and the coefficients of linear correlation.
Введение
Намечаемое на ближайшее время значительное расширение экономических связей РФ с Китаем и создание высокоскоростной автотрассы Западная Европа - Западный Китай делает актуальным исследование движения автотранспорта по протяжённым магистралям. Исследование эффективности движения по трассам целесообразно проводить с помощью имитационного моделирования (ИМ). За объект исследования в данной работе было взято движение автотранспорта по участку строящейся международной трассы Западная Европа -Западный Китай [1], проходящей по Татарстану, а точнее участок трассы от моста в устье Камы до города Альметьевска. Длина участка трассы составляет примерно 160 км. В дальнейшем результаты исследования планируется использовать и для моделирования других участков строящейся трассы. Строительство трассы началось в 2008 году, окончание планируется к 2019 году [1]. Реализация проекта позволит значительно сократить время доставки грузов с Запада на Восток по сравнению с существующими альтернативными возможностями железнодорожному по Транссибу и морскому по Суэцкому каналу.
Основные характеристики международной
трассы:
• Общая протяженность трассы- 8445 км.
• Техническая категория - автомагистраль, скоростная дорога, шоссе, трасса.
• Количество полос движения - 4 или 6.
• Расчетная скорость движения - 80 - 150
км/ч.
• Климатические условия -континентальные, резко континентальные.
Дополнительные характеристики участка трассы, взятогов качестве объекта исследования:
• Протяженность - от 20 до 160 км.
• Количество полос движения - 4, по две в каждом направлении.
• Количество съездов/въездов на трассе - 4.
• Виды транспортных средств - легковые, одногабаритные, двухгабаритные и трехгабаритные грузовые.
Методика исследования
Целью данного исследования является определение основных показателей движения грузового автотранспорта при перевозке грузов тремя видами грузовиков различной грузоподъёмности: одногабаритного,
двухгабаритного и трёхгабаритного. Методика, предлагаемая для исследования движения автотранспорта, основанная на ИМ описана в [2]. Примеры применения методики приведены в [3], [4]. В методике, выделены следующие 9 этапов.
1. Выбор результативных показателей эффективности движения автотранспорта по трассам и влияющих на них факторов,
2. Постановка задач.
3. Разработка блок-схемы и программы имитационной модели.
4. Разработка стратегического плана проведения имитационного эксперимента.
5. Вычисление коэффициентов линейной корреляции и их анализ.
6. Проведение имитационного моделирования по стратегическому плану.
7. Построение математической модели, состоящей из совокупности уравнений регрессии, представляющих собой зависимости результативных показателей эффективности от влияющих на них факторов.
8. Вычисление коэффициентов эластичности влияния факторов на изменение результативных показателей эффективности движения автотранспорта.
9. Выводы по результатам исследования. Далее приведём описание всех выделенных
этапов.
Выбор результативных показателей эффективности движения автотранспорта и влияющих на них факторов
Для исследования было отобрано 15 результативных показателей эффективности (табл. 1).
Таблица 1 - Перечень отобранных переменных
№ Код перемен-ной Наименование переменной
Изменяемые факторы
1 Х1 Среднее время между поступлением на трассу одногабаритных грузовиков в мин.
2 Х2 Среднее время между поступлением на трассу двухгабаритных грузовиков в мин.
3 хз Среднее время между поступлением на трассу трехгабаритных грузовиков в мин.
4 х4 Среднее время между поступлением на трассу легковых автомобилей в мин.
5 Х5 Протяженность трассы в км.
Результативные показатели
6 У1 Среднее время проезда по полной длине трассы легковых автомобилей в мин.
7 У2 Стандартное отклонение времени проезда по полной длине трассы легковых автомобилей в мин.
8 Уз Средняя скорость движения по трассе легковых автомобилей в км/ч.
9 У4 Среднее время проезда по полной трассе одногабаритных грузовиков в мин.
10 У5 Стандартное отклонение времени проезда по полной длине трассы одногабаритных грузовиков в мин.
11 Уб Средняя скорость движения по трассе одногабаритных грузовиков в км/ч.
12 у7 Среднее время проезда по полной трассе двухгабаритных грузовиков в мин.
13 у8 Стандартное отклонение времени проезда по полной длине трассы двухгабаритных грузовиков в мин.
14 у9 Средняя скорость движения по трассе двухгабаритных грузовиков в км/ч.
15 у10 Среднее время проезда по полной трассе трехгабаритных грузовиков в мин.
16 у11 Стандартное отклонение времени проезда по полной длине трассы трехгабаритных грузовиков в мин.
17 у12 Средняя скорость движения по трассе трехгабаритных грузовиков в км/ч.
18 у13 Коэффициент заполнения трассы в долях от единицы
19 У14 Общий вес всех перевезенных грузов
20 У15 Общая стоимость перевозки всех грузов в тыс. руб.
21 У1б Стоимость перевозки всех грузов в тыс. руб., отнесённая к 1 тонне на 100 км.
Постановка задач
Главная задача исследования, основанного на ИМ - построение математической модели, представляющей собой совокупность уравнений регрессии, связывающих результативные показатели эффективности движения автотранспорта с влияющими на них факторами.
у} _ f] х2,--хт); 7 = 1к, (1)
где У7 - }-й показатель эффективности движения автотранспорта по трассе; х, - 1-й фактор, влияющий на движение автотранспорта по трассе; к— количество результативных показателей эффективности движения автотранспорта по трассе; т— общее количество факторов.
По зависимостям (1) можно произвести оценку степени влияния факторов на результативные показатели эффективности движения автотранспорта по трассе по коэффициентам эластичности.
Перед проведением регрессионного анализа целесообразно вычислить коэффициенты линейной корреляции между переменными, отобранными для исследования, по которым также можно произвести оценку степени влияния факторов на результативные показатели эффективности движения автотранспорта по трассе.
Коэффициенты линейной корреляции вычисляются по п экспериментальным данным по следующей формуле:
_ Щху - т1*х -Чу (2)
стХ-стУ
В формуле (2) оценки математических ожиданий переменных х,у и их произведения вычисляются по формулам:
1 N
'1*х _ — V
1 х N ^
х'; т,* _ — V у,
У N '
1
"1хУ _ N £
Оценки вторых начальных моментов требуются для вычисления средних квадратических отклонений. Для этого используются следующие формулы:
N
V 1 _1 х,2; т 2* у
»2
'2 у - И1у-
4 V у-2
'2х "Ч х =
По статистическим таблицам находится критическое значение коэффициента линейной корреляции.
1 2 (3)
г _ +1 крит ■ крит + п - 2
Имитационная модель движения автотранспорта по протяженным трассам
При разработке имитационной модели автотранспортной трассы за основу берется клеточная структура представления транспортной сети. При таком подходе вся трасса разбивается на ячейки, где размер каждой ячейки соответствует размеру стандартного легкового автомобиля плюс допустимая величина интервала между автомобилями 5+3=8 метрам. Размер крупных
х, - у,
, _ 1
т
2 х
1 _ 1
у
грузовых транспортных средств будет соответствовать одной, двум или трем ячейкам трассы в зависимости от своих габаритов. Исходя из клеточной структуры представления
автотранспортной сети, движение автомобиля подразумевает его перемещение с одной клетки на другую. Необходимо также учитывать, что автомобили перемещаются с различной скоростью. Трасса двухполосная. Поэтому обгон может совершаться по аналогичному принципу «с клетки на клетку». Для разработки имитационной модели использована система GPSSW.
Стратегическое планирование
имитационных экспериментов
Стратегическое планирование - определяет количество проводимых экспериментов, порядок их выполнения и значения изменяемых факторов в каждом эксперименте. Для стратегического планирования используется концепция «черного ящика», суть которой - абстрагирование от физической сущности процессов, происходящих в моделируемой системе и выдаче заключений о ее функционировании только на основании входных и выходных переменных. Входные, независимые переменные называются факторами. Выходные - откликами, их величина зависит от значений факторов и параметров объекта исследования (ОИ).
В стратегическом плане значения факторов задаются в кодированном или натуральном виде. В кодированном виде для каждого фактора указывается интервал от -1 до +1, в пределах которого проводится исследование. Выскажем допущение, что уравнения регрессии (1) являются нелинейными, нелинейность которых не превышает второй степени. Для их построения используем стратегический план, состоящий из центральной точки, вариантов составленных по плану дробного факторного эксперимента (ДФЭ) плюс «звёздные» точки. В плане ДФЭ первых три фактора берём за основные и для них составляем план ПФЭ, а два дополнительных фактора меняем по закону изменения произведений основных факторов: х4 —х1*х2,'х5 —х1*хз. Тогдаколиче-ство вариантов стратегического плана составит:
N = 1 + 2г + 2т = 1 + 23 + 2 • 5 = 19, (4)
где г- количество основных факторов; т - общее количество факторов.
Имитационное моделирование движения автотранспорта по загородным трассам
По стратегическому плану проведено имитационное моделирование процесса по 19 вариан-там.По результатам проведенного имитационного моделирования можно сделать следующий вывод: с увеличением длины трассы (х5) скорости движения любого автотранспорта: у3, у6, у9,у12увеличиваются, а коэффициент заполнения трассы (у13) уменьшается. При уменьшении среднего времени между поступлением автотранспорта (другими словами, автотранспорт поступает чаще), скорость автомобилей падает, и загруженность трассы растет. Полученные результаты логичны и не противоречат здравому смыслу.
Вычисление коэффициентов линейной корреляции
Перед проведением регрессионного анализа целесообразно вычислить коэффициенты линейной корреляции между переменными, отобранными для исследования, по которым можно оценить степень влияния факторов на результативные показатели эффективности движения автотранспорта по трассе.На рис.1 приведены диаграммы коэффициентов линейной корреляции между факторами х' = 1,5 и средней скоростью движения
трёхгабаритных грузовиков - у12 и коэффициентом заполнения автотранспортом трассы - у 13.
-0,05
-од
-0Д5 -0,2 -0.25 -0,3 -0,35 -0,4 -0,45 -0,5
У
б
Рис. 1 - Коэффициенты линейной корреляции междуфакторами х' = 1,5 и результативными показателями эффективности движения автотранспорта: а) средней скоростью движения трёхгабаритных грузовиков - у12; б) коэффициентов заполнения автотранспортом трассы - у13
Рассмотрим построенные диаграммы.На рис.1 а для результативного показателя - скорости движения автотранспорта у12наблюдается положительное влияния всех факторов. Наибольшее сильное влияние имеет фактор х5 - длина трассы. На рис.1б для результативного показателя у13отметим отрицательное влияние всех факторов. Вычисленные значения коэффициентов линейной корреляции позволяют сделать вывод, чтоввиду того, что сила связи между результативными показателями и влияющими на них факторами по абсолютной величине не превышает 0,5, то связь между ними нелинейная. Поэтому в уравнения регрессии введём все факторы в первой и второй степени.
Построение математической модели
Так как все переменные, используемые в исследовании, являются случайными, независимыми и непрерывными, то для построения математиче-
а
ской модели будем использовать регрессионный анализ, основанный на методе наименьших квадратов (МНК) [2]. Приведём математическую модель движения автотранспорта по загородной трассе, состоящую из 16 уравнений регрессии. Для построения уравнений регрессии использован пакет прикладных программ 81аИ8Йса 6.0 [5].
При построении математической модели поставлены следующие требования. Отношения стандартных ошибок к средним значениям не должно превышать 0.15, уровни значимости по коэффициентам множественной детерминации и критериям Фишера не должны превышать 0.05.
Приведём построенную математическую модель движения автотранспорта по загородной трассе, состоящую из 16 уравнений регрессии.
у1 =77,60085 -5,27110*х1 - 5,05549* х2 -5,71830* х3 - 9,19535* х4 + 0,68850* х5 + 0,33237*х1*х1 + 0,33935* х2* х2 + 0,37017* х3* х3 + 0,68360* х4* х4- 0,00002* х5* х5;
у2 =35,08224 -2,38796* х1 - 2,17986* х2 -2,39806* х3 - 4,39816* х4 + 0,02044* х5 + 0,14736* х1* х1 + 0,14434* х2* х2 + 0,14463* х3* х3 + 0,34242* х4*х4- 0,00001*х5*х5;
у3 = -12,0206+5,0354* х1 +6,5281* х2 + 7,1698* х3 +10,6807* х4 + 0,2038* х5 - 0,4949* х1* х1 -0,5010* х2* х2 -0,5365* х3* х3 -0,8467* х4* х4 -0,0001* х5* х5;
у4 = 79,24260 -5,25694* х1 - 5,27947* х2 -5,74293* х3 - 9,31465* х4 + 0,67476* х5 + 0,33509* х1* х1 + 0,35950* х2* х2 + 0,36999* х3* х3 + 0,69627* х4* х4 + 0,00002* х5* х5;
у5 = 34,95741 -2,21942* х1 - 2,25450* х2 -2,43878* х3 - 4,42267* х4 + 0,02361* х5 + 0,13054* х1* х1 + 0,15115* х2* х2 + 0,14952* х3* х3 + 0,34209* х4*х4- 0,00002*х5*х5;
у6 =-13,2983+5,0141* х1+ 6,7269* х2 + 7,1625* х3 +10,7385* х4 + 0,2164* х5 - 0,5009* х1* х1 -0,5174* х2* х2 - 0,5286* х3* х3 - 0,8575* х4* х4 -0,0002* х5* х5;
у7 = 76,45790 -5,03301* х1 - 5,19193* х2 -5,51027* х3 - 9,10047* х4 + 0,68377* х5 + 0,31673* х1* х1 + 0,35955* х2* х2 + 0,35506* х3* х3 + 0,68877* х4*х4- 0,00003*х5*х5;
у8 = 34,73600 -2,17974* х1 - 2,18539* х2 -2,42276* х3 - 4,41426* х4 + 0,01800* х5 + 0,13422* х1* х1 + 0,14745* х2* х2 + 0,14859* х3* х3 + 0,34382* х4*х4- 0,00002*х5*х5;
у9 = -12,3458 +4,9605* х1 +6,7887* х2 + 7,0253* х3 +10,6662* х4 + 0,2169* х5 - 0,4867* х1* х1 -0,5351* х2* х2 -0,5277* х3* х3 -0,8627* х4* х4 -0,0001* х5* х5;
у10 = 79,30235 -5,30177* х1 - 5,32621* х2 -5,56427* х3 - 9,46047* х4 + 0,68202* х5 + 0,33054*
х1* х1 + 0,36332* х2* х2 + 0,35646* х3* х3 + 0,70989* х4* х4 - 0,00001* х5* х5;
у11 =35,55981 -2,229861* х1 - 2,31788* х2 -2,56178* х3 - 4,32307* х4 + 0,01755* х5 + 0,14197* х1* х1 + 0,15671* х2* х2 + 0,15997* х3* х3 + 0,33516* х4* х4 - 0,00000* х5* х5;
у12 = -13,5124 +5,1076* х1 +6,7839* х2 + 7,0217* х3 +10,8625* х4 + 0,2153* х5 -0,4960* х1* х1 -0,5308* х2* х2 -0,5232* х3* х3 -0,8767* х4* х4 -0,0001* х5* х5;
у13 =1,557847 -0,057116* х1 - 0,109675* х2 -0,152679* х3 - 0,096475* х4 -0,002712* х5 + 0,005403* х1* х1 + 0,008383* х2* х2 + 0,011485* х3* х3 + 0,007240*х4*х4 +0,000005*х5* х5;
у14 =31658,13 +2681,7* х1 + 138,07* х2 -3037,22* х3 - 152,34* х4+5,16* х5 - 187,08* х1* х1 -1,7* х2* х2 +213,41* х3* х3 -9,78* х4* х4 -0,03* х5* х5;
у15 =-170671 +378301 * х1 - 207569* х2 -182641* х3+ 13507* х4+59130* х5 - 30203 * х1* х1 +15477* х2* х2 +20601* х3* х3 -1160* х4* х4 +3* х5* х5;
у16=-9,15135 +1,02140* х1 -7,36574* х2 + 6,42650* х3+ 0,77142* х4+1,94263* х5 + 0,06476* х1* х1 +0,50137* х2* х2 -0,38480* х3* х3 -0,01504* х4* х4 + 0,00002 * х5* х5. (5)
Показатели качества полученных уравнений регрессии приведены в таблице 2.
Таблица 2 - Показатели качества полученных уравнений регрессии
Код Средн. УГ Станд. ошибка От- эш. средн. Коэф. множ. детер. Крит. Фишера Уров. знач. по крит. Фишера Р^
у1 71,864 3,8891 0,05412 0,997 181,4 ,00000
у2 6,473 1,8091 0,27948 0,985 27,828 ,00004
у3 73,948 10,297 0,13924 0,943 6,481 ,00708
у4 72,068 3,9759 0,05516 0,997 171,6 ,00000
у5 6,463 1,7837 0,27595 0,986 28,631 ,00004
у6 73,738 10,556 0,14315 0,942 6,3618 ,00751
у7 71,791 3,4324 0,04781 0,998 227,08 ,00000
у8 6,436 1,7287 0,26857 0,986 29,135 ,00003
у9 73,76 10,275 0,13929 0,944 6,5854 ,00672
у10 72,108 3,6861 0,05111 0,998 200,1 ,00000
уц 6,5303 1,7948 0,27484 0,986 28,503 ,00004
у12 73,596 10,368 0,14087 0,944 6,5693 ,00678
у13 0,368 0,089 0,14323 0,976 16,188 ,00030
у14 30582 1056,4 0,034542 0,985 27,590 ,00004
у15 5313,8 33,64 6,33Е-06 0,997 156,10 ,00000
у16 174,17 9,2388 0,053044 0,998 220,43 ,00000
Результаты нелинейной множественной регрессии позволяют сделать следующие выводы.
Во-первых, получены высокие значения коэффициентов множественной детерминации в пределах от 0,942 до 0,998.Критическое значение коэффициента множественной детерминации для рекомендуемого уровня значимости а=0.05 при двух степенях свободы к1—10 и к2—19-1—18 равно 0,573 [5].
Во-вторых, критерий Фишера лежит в пределах от 6,3618 до 227,08 с уровнем значимости от 0 до 0,00751.
В-третьих, отношение стандартных ошибок к средним значениям для наиболее существенных результативных показателей движения автотранспорта по загородной трассе - скоростей движения: у3, у6, у9, у12, у 13 не превышает 0,15, что следует считать вполне приемлемым.Отношение стандартных ошибок к средним значениям для стандартных от-клоненийвремени движения автотранспорта по трассе для результативных показателей: у2, у5, у8, у11 превышает значение 0,15, но с этим можно смириться ввиду меньшей существенности этих показателей по сравнению с остальными.
Таким образом, применение нелинейной множественной регрессии позволило достичь вполне приемлемых результатов. По всем ограничениям, наложенным на уравнения регрессии, получены удовлетворительные значения для основных результативных показателей движения автотранспорта. Следовательно, построенная регрессионная модель (5) адекватно отображает зависимость результативных показателей эффективности движения автотранспортных средств по трассе от факторов. Построенная математическая модель позволяют оценить влияние изменяемых факторов = 1,5 на результативные показатели движения автотранспорта, например, на скорость движения трёхгабаритного автотранспорта.
Вычисление коэффициентов эластичности влияния факторов на изменение результативных показателей эффективности
Коэффициент эластичности г - го фактора в ] - ом результативном показателе эффективности вычисляется по формуле:
А у
£а(
; г = 1, т; у = 1, к,
(6)
где А у - изменение у-го результативного показателя эффективности, вызываемое изменением г-го фактора; т - общее количество изменяемых факторов.
Вычисленные значения коэффициентов эластичности факторов в средней скорости движения трёхгабаритных грузовиков - у 12 и в коэффициенте заполнения автотранспортом трассы - у13 приведены на рис. 2.
Гистограмма, приведенная на рис.2а, позволяет сделать вывод, что на повышениесредней скорости движения по трассе трехгабаритного грузового автотранспорта - у12 положительное влияние оказывают все факторы. Наиболее сильное влияние оказывает фактор х5- длина трассы. Гистограмма,
приведенная на рис.2б позволяет сделать вывод, что на повышение значения показателя заполнения трассы - у13 отрицательное влияние оказывают все факторы. Так как влияние изменяемых факторов на результативные показатели эффективности по коэффициентам линейной корреляции (рис.1) и коэффициентам эластичности (рис.2) идентичны в качественном плане делаем заключение об их достоверности.
0,35
0,3
г 0,25
ч 0,2
¡5
ОД 5
0,1
-е-
1 0,05
я
Факторы
а
-0,05 -ОД -0,15 -0,2 -0,25 -0,3 -0,35 -0,4 -0,45
Факторы б
Рис. 2 - Коэффициенты эластичности факторов х^г = 1,5: а) в средней скорости движения трёхгабаритных грузовиков - у12; б) в коэффициенте заполнения автотранспортом трассы - у13
Введём интегрированный показатель - интенсивность поступления на трассу всех видов автотранспорта, выражаемая в количестве автотранспорта поступающих на трассу за одну минуту.
х6—1/х1 + 1/х2+1/х3+1/х4. Примем х1—х2—х 3—х4 тогда х6—4/х1.
Будем менять значения х6—4/х1 с одинаковым шагом 0,5 начиная с 0,5 до 5 и по математической модели вычислим у12и у13. Затем построим два графика: у12—/(х6); у13—/(х6), приведённые на рис.3.
По рис.3.б можно определить нежелательную интенсивность поступления автотранспорта на трассу выше которой скорость автотранспорта падает ниже нежелательного значения 60 км/час для трасс трёх различных протяжённостей: 20, 90, 160 км.
По рис.3 отметим имеющуюся разницу в показателях на трассах различной протяжённости.
На рис.4 приведены графики зависимостейу12—((у13) при общей интенсивности движения автотранспорта.
По рис.4 можно определить коэффициент заполнения трассы, при котором скорость движения автотранспорта падает ниже нежелательного значения 60 км/час.для трасс трёх различных протяжённостей: 20, 90, 160 км. Эти значения можно пропускной способностью трассы. Аналогично можно найти допустимые значения интенсивности поступ-
ления автотранспорта на трассу для других допустимых значений скоростей.
а
б
Рис. 3 - Графики зависимости результативных показателей эффективности движения автотранспорта по трассе от общей интенсивности движения автотранспорта: а) скорости движения трёхгабаритных грузовиков, б) коэффициента заполнения трассы
трёхгабаритных грузовиков от коэффициента заполнения трассы
Заключение
В работе получены следующие основные результаты:
1. Проведен предварительный анализ работы движения автотранспорта по трассе, определен перечень результативных показателей и влияющих на них факторов, произведена постановка задач исследования и предложена методика ее решения.
2. Разработан стратегический план имитационного эксперимента. Проведено имитационное моделирование движения автотранспорта по загородной трассе по стратегическому плану.
3. Проведен корреляционный анализ.
4. Разработана статистическая модель, представляющая собой совокупность уравнений регрессии, которые отображают зависимость результативных показателей от влияющих на них факторов.
5. Вычислены коэффициенты эластичности влияния факторов на изменение результативных показателей эффективности.
6. Проведенное исследование показало высокую эффективность применения статистических методов для оценки влияния факторов на движение автотранспорта по загородной трассе. Использование прикладных программ показало их пригодность для обеспечения выполнения всех этапов предложенной методики.
7. Предложенная и опробованная в работе методика вполне пригодна для проведения аналогичных исследований различных участков трасс, и трасс в целом.
8. Результаты работы докладывались на двух международных конференциях [7], [8].
Литература
1. www.intermost.ru
2. Якимов И.М. Компьютерное моделирование: Учебное пособие - Казань: Изд-во Казан. гос. техн. ун-та, 2008. 220 с.
3. ЯкимовИ.М., Бетретдинов Р.З., Моделирование и оптимизация процессов разработки систем программного обеспечения. Вестник Казан. технол. ун-та, 15, 18, С. 261-265 (2012).
4. Якимов И.М., Кирпичников А.П. Сытник А.С., Математическая модель работы грузоподъёмных кранов и оптимизация их обслуживания. Вестник Казан. технол. унта. 16,24. С.182-187 (2013).
5.Боровиков В.П., Боровиков И.П. Statistica: Статистический анализ и обработка данных в среде Windows. - М.: Информационно-издательский дом "Филин", 1997. - 608с.
6.Шмойлова Р.А., Минашкин В.Г., Садовников Н.А., Шувалова Е.Б. Теория статистики. 4-е издание.: М.: Финансы и статистика, 2004. - 656с.
7.Маркина О.В., «Информационные технологии на службе общества», материалы международной науч.-практ. конф., Нижнекамск, 18 апреля 2014 г.
8.Маркина О.В., Якимов И.М., «Наука в современном информационном обществе», материалы III международной научно-практической конференции, том 2,NorthCharleston, USA, 10-11 апреля 2014 года.
© И. М. Якимов - канд. техн. наук, проф. каф. автоматизированных систем обработки информации и управления КНИГУ им А.Н.Туполева; А. П. Кирпичников - д-р физ.-мат. наук, зав. кафедрой интеллектуальных систем и управления информационными ресурсами КНИТУ, [email protected]; О. В. Маркина - магистр каф. автоматизированных систем обработки информации и управления КНИТУ им А.Н.Туполева, [email protected].
© 1 М. Yakimov- PhD, Professor of the Department of Automated Information Processing Systems & Control, KNRTU named after A.N. Tupolev; А. P. Kirpichnikov - Dr. Sci, Head of the Department of Intelligent Systems & Information Systems Control, KNRTU, [email protected]; O. V. Markina-- Master of the Department of Automated Information Processing Systems & Control, KNRTU named after A.N. Tupolev, [email protected].