Математическое моделирование транспортных потоков на основе теории клеточных автоматов
ОГЛАВЛЕНИЕ
Введение
Актуальность
Современное состояние исследований
Цель и задачи диссертационной работы
Научная новизна
Теоретическая и практическая значимость работы
Методика исследования
Достоверность полученных результатов
Апробация работы, внедрение
Основные положения, выносимые на защиту
Объем и структура диссертации
Глава 1. Обзор основных методов и подходов к транспортному
моделированию
1.1 Различные подходы к моделированию трафика.
1.2 Макроскопические модели.
1.3 Мезоскопические (кинетические) модели.
1.4 Микроскопические модели.
1.5 История развития клеточных автоматов.
1.6 Теория клеточных автоматов применительно к транспортному
моделированию.
1.7 Обзор моделей на основе теории клеточных автоматов.
1.7.1 Правило 184 Вольфрама
1.7.2 Модель Нагеля-Шрекенберга
1.7.3 Модель Кернера-Кленова
1.7.4 Подход Буслаева – моделирование на разветвленных структурах
1.7.5 Другие модели
1.8 Современные исследования транспортных потоков
Выводы по Главе 1
Глава 2. Математическая модель транспортных потоков на основе теории
клеточных автоматов для моделирования трафика на городской улично-
дорожной сети. Основные алгоритмы движения. Различные элементы УДС,
примеры тестовых расчетов
2.1 Построение двумерной модели
2.2 Общий алгоритм движения
2.2.1 Алгоритм а1- общий алгоритм смены полосы
2.2.2 Алгоритм а2 – общий алгоритм движения вперед
2.3 Достижение цели
2.4 Движение на элементах улично-дорожной сети
2.4.1 Алгоритм а3 перестроения с достижением цели (движение по
участку дороги с выездом)
2.4.2 Алгоритм а4 прохождения четырехстороннего перекрестка со
светофором
2.4.3 Алгоритм а5 – обеспечение безопасности на перекрестке
2.4.4 Алгоритм а6 прохождения Т-образного перекрестка со светофором
2.4.5 Алгоритм а7 прохождения участка с сужением/расширением
2.4.6 Алгоритм а8 въезда с второстепенной дороги на главную
2.4.7 Алгоритм а9 проезда участка с разворотом
2.4.8 Алгоритм а10 проезда участка с пешеходным светофором
2.4.9 Алгоритм а11 объезда широкого неподвижного препятствия
2.5 Различные стратегии поведения водителей
2.5.1 Алгоритм а12 «вежливый водитель» для смены полосы
2.5.2 Алгоритм а13 «вежливый водитель» для въезда с второстепенной
дороги на главную
2.5.3 Алгоритм а14 снижения скорости при наличии затора
2.6 Модель транспортных потоков с «медленным стартом»
2.6.1 Алгоритм а15 «медленный старт»
2.7 Граничные условия в модели
2.7.1 Алгоритм а16 задания граничных условий с возможностью
образования очереди на въезд в расчетную область
Выводы по Главе 2
Глава 3. Численная реализация. Программный комплекс CAM-2D
Пользовательский интерфейс и визуализация. Параллельная реализация
Расчеты на сети
3.1 Структура последовательной программы
3.2 Основные элементы и структуры
3.3 Параллельная реализация
3.3.1 Алгоритмы а17 и а18 выезда и въезда автомобиля на расчетный
элемент
3.4 Интеграция с модулями Визуализации и Пользовательского интерфейса
3.4.1 Модуль пользовательского интерфейса и визуализации данных
3.4.2 Интеграция расчетного модуля CAM-2D с модулем
пользовательского интерфейса и визуализации
3.5 Расчеты на сети
Выводы по Главе 3
Глава 4. Верификация модели: Сравнение с Aimsun TSS, воспроизведение
экспериментально наблюдаемых пространственно-временных структур
4.1 Сравнение с Aimsun TSS
4.1.1 Тестовая задача 1 о пропускной способности перекрестка
4.2 Сравнение с экспериментальными данными
4.2.1 Различные способы осреднения данных
4.2.2 Тестовая задача 2 о въезде на автомагистраль с постоянными
входящими потоками. Сравнение методов осреднения.
4.2.3 Тестовая задача 3 о въезде на магистраль с переменным потоком Q1
на главной дороге. Сравнение результатов с данными, полученными при
помощи модели Кернера-Конхойзера
4.2.4 Тестовая задача 4 о движении на участке магистрали с выездом и
въездом
Выводы по Главе 4
Заключение
Список литературы
Приложение 1 Лицензионный договор и акт приема-передачи по проекту
«Разработка архитектуры, принципов построения и программного
обеспечения суперкомпьютерной информационно-аналитической платформы
прикладного моделирования, прогнозирования и экспертиз транспортных
процессов и систем в виде сетевой компьютерной лаборатории»
Приложение 2 Свидетельство о государственной регистрации программы
CAM-2D
Во введении приведены краткие сведения об актуальности выбранной
темы исследования, сделан обзор современного состояния исследований по
теме, указаны цели и задачи диссертационной работы, научная новизна,
методика исследований, теоретическая и практическая значимость работы,
достоверность полученных результатов, приведен список конференций, на
которых представлялась работа, и список грантов РФФИ, включающих
полученные в рамках работы над диссертацией результаты. Приведены
основные положения, выносимые на защиту, а также объем и структура
диссертационной работы.
В первой главе содержится обзор основных методов и подходов к
транспортному моделированию – макроскопического, мезоскопического и
микроскопического.
Первый подход – макроскопический, основанный на аналогии движения
потока автомобилей и потока сжимаемой жидкости или газа. В основе моделей
– адаптированные уравнения гидрогазодинамики, чаще всего одномерные. В
уравнения добавляют дополнительные члены, которые отвечают за
«человеческий фактор». При этом работа идет с усредненными величинами,
характеризующими поток, такими как плотность и скорость потока. Отдельные
автомобили в данном подходе не рассматриваются.
Главный плюс макромоделей – в экономии вычислительных ресурсов,
поскольку время расчета не растет с ростом числа автомобилей в системе (то
есть, плотности). Однако, приближение сплошной среды, к которому относятся
макромодели, накладывает свои ограничения: для того, чтобы считать поток
сплошной средой, расстояние между автомобилями должно быть порядка
размеров самих автомобилей. Это соответствует синхронизированному, то есть
достаточно плотному, потоку, где автомобили должны подстраивать свою
скорость под скорость потока и скорость движения меньше свободных
скоростей для данной дороги. Другое ограничение, вытекающее из
приближения сплошной среды – это то, что возможно рассмотрение только
достаточно протяженных участков дорог, ведь в масштабах десятков и сотен
метров физическими размерами автомобилей нельзя пренебречь.
Помимо перечисленных ограничений, есть еще ограничения, связанные со
сложностью описания реакции водителей на светофоры и разметку, что важно
при моделировании небольших участков с целью решений локальных задач
управления дорожным движением.
Наконец, стратегии водителей в макромоделях должны быть одинаковыми.
Второй возможный подход – мезоскопический, он основан на уравнениях,
выводящихся из кинетических соображений, аналогично уравнениям
физической кинетики. Основная переменная в этом подходе – фазовая
плотность потока, представляющая собой распределение автомобилей по
координате и скорости. Основные недостатки модели – те же, что и для
макроскопических моделей.
И, наконец, третий существующий подход – микроскопический. В нем мы
имеем дело с каждым автомобилем в отдельности. Существенная часть
микроскопических моделей относятся к моделям следования за лидером, в
центре внимания в них взаимодействие ведомого и головного автомобиля. Как
правило, модели этого класса используют уравнения движения автомобилей, но
есть и исключение – модели на основе теории клеточных автоматов.
Работа с каждым автомобилем отдельно позволяет сделать модель более
гибкой для описания процессов на «микроуровне» – то есть в масштабах
отдельных перекрестков. Микроскопический подход не имеет ограничений на
плотность потока, можно описывать как очень разреженное, так и плотное
движение. Запись уравнений движения в явном виде позволяет точнее
описывать паттерны торможения и ускорения автомобилей в потоке в
зависимости от скоростей и расстояний между ними, что важно для
воспроизведения экспериментальных свойств трафика. С развитием
компьютерных технологий проблема нехватки вычислительных ресурсов для
расчетов с большим числом автомобилей на масштабных сетях уже не
представляет проблемы.
Приводятся основные модели, относящиеся к каждому из подходов, а
также рассматривается история развития клеточных автоматов как
самостоятельной математической теории. Отдельный подраздел главы
посвящен применению теории клеточных автоматов к транспортному
моделированию.
Впервые теорию клеточных автоматов к транспортному моделированию
применили Нагель и Шрекенберг в 1992 году. Трасса в модели Нагеля-
Шрекенберга представляется в виде одномерной решетки, каждая ячейка
которой может быть либо пустой, либо содержать частицу, обозначающую
транспортное средство. Частицы перемещаются из одной ячейки в другую
(свободную) в одном направлении. В случае однополосного движения они не
могут обгонять друг друга. Вся система – пространство, время, скорость –
дискретна. Предметом рассмотрения является состояние всей системы в
отдельные последовательные моменты времени. Скорость показывает, на
сколько ячеек автомобиль перемещается за один шаг по времени. Ускорение
происходит мгновенно между шагами. На каждом слое по времени происходит
обновление состояния системы по определённым правилам с определёнными
вероятностями. Обновление в модели Нагеля-Шрекенберга состоит из четырех
шагов:
1) Ускорение. Скорость автомобиля i увеличивается на единицу, если
максимальная разрешенная скорость не достигнута:
Vi → min(Vi + 1, Vmax )
2) Торможение. Скорость автомобиля уменьшается на единицу, если есть
угроза столкновения с впереди идущим автомобилем:
→ min( , − 1),где – расстояние до впереди идущего
автомобиля.
3) Случайные возмущения. Если скорость автомобиля положительна, то она
может быть уменьшена на единицу с некоторой вероятностью:
→ max( − 1, 0) с вероятностью p.
4) Движение. Каждый автомобиль продвигается вперед на количество ячеек,
соответствующее его новой скорости после выполнения предыдущих
шагов:
→ + .
Для упрощения записи считаем, что скорость и расстояние измеряются в
ячейках, а время безразмерно. По этой причине величины можно складывать,
вычитать и сравнивать друг с другом.
Каждый шаг по времени состояние ячеек обновляется, автомобили
передвигаются вдоль решетки на целое число ячеек.
Длина ячейки равна 7.5 м, что соответствует средней длине участка
дороги, занимаемым автомобилем, стоящим в пробке (исходя из
ячеек
экспериментальных значений). Максимальная скорость = 5.
шагповремени
Исходя из проведенных вычислительных экспериментов, средняя скорость,
ячейки
отвечающая свободному движению, оценивается как 4.5; этому
шагповремени
значению ставится в соответствие экспериментальное значение скорости
свободного движения – 120 км/ч (33,3 м/с). Простой расчет позволяет, исходя
из этих данных, оценить шаг по времени, который получается примерно
метровячеек
7.5×4.5
ячейкушагповремени
равным 1 секунде: ≈ 1с.
33.3м/с
∆
Отсюда следует, что «шаг по скорости» ∆ = = 27км/ч.
∆
Эти значения используются в большинстве созданных впоследствии
моделей.
Модель Нагеля-Шрекенберга до сих пор остается одной из важнейших
моделей транспортных потоков на основе теории клеточных автоматов,
множество существующих современных моделей основано на четырех
правилах Нагеля-Шрекенберга. В первой главе диссертационной работы далее
приводится обзор современных моделей на основе теории клеточных
автоматов.
Во второй главе представлена созданная автором диссертационной
работы двумерная математическая модель транспортных потоков на основе
теории клеточных автоматов, а также алгоритмы, входящие в нее.
В основу многополосной модели положена одномерная модель Нагеля-
Шрекенберга, которая обобщается на двумерный случай. Трасса в этом случае
представляет собой двумерную решетку, где количество ячеек в поперечном
направлении соответствует числу полос трассы, как показано на рисунке 1:
Рисунок 1 – Расчетная область в обобщенной двумерной микроскопической
модели
В многополосной модели разрешены перестроения машин из полосы в
полосу и обгоны. Процесс обновления состояний ячеек делится на два шага:
– смена полосы при необходимости и наличии возможности;
– движение вперед по правилам однополосного движения.
Размер ячейки расчетной области равен 7,5 метрам в длину и 3,5 метрам в
ширину. Эти параметры соответствуют средней площади, занимаемой
легковым автомобилем, стоящим на дороге, с учетом минимального
необходимого расстояния до соседних автотранспортных средств. Ширина
ячейки соответствует стандартной ширине полосы автомобильной дороги.
В двумерную микроскопическую модель в качестве параметров входят две
вероятностные характеристики: вероятность спонтанного снижения скорости из
модели Нагеля-Шрекенберга и вероятность остаться в своей полосе в случае,
когда перестроение возможно и выгодно водителю. Первый параметр
соответствует ситуациям, когда водитель, например, посчитал состояние
дорожного полотна потенциально опасным или вынужден отвлечься на что-то
во время движения и по этой причине принимает решение временно снизить
скорость. Второй параметр характеризует наличие водителей, не склонных
лишний раз менять полосу движения, даже если это позволит ехать быстрее.
Другой причиной отказа от перестроения может быть нерешительность
водителя или отсутствие опыта.
Как и в случае одномерного движения, АТС характеризуется текущей и
максимальной скоростью. Скорость принимает целые значения от 0 до 4, а при
высоких скоростных режимах, и более. Максимальная скорость равна =
ячейки
4. Если шаг по времени равен одной секунде, то дискретному
шагповремени
значению скорости 4 соответствует скорость 108 километров в час, то есть
скорость свободного движения автомобиля по автомагистрали. Так как
правилами дорожного движения в России скорость ограничена 110 км/ч, нам
ячеек
нет смысла рассматривать = 5= 135км/ч, как в модели
шагповремени
Нагеля-Шрекенберга.
Ускорение в модели – это число единиц, на которое автомобиль может
поменять скорость за один временной шаг. Важно сопоставлять это значение с
реальным возможным ускорением автомобиля при движении. С этой точки
зрения для более адекватного описания физического движения автомобиля при
заданных параметрах ячеек и шага по времени не имеет смысла менять
скорость автомобиля более чем на одну единицу за один шаг по времени.
Модель обеспечивает возможность достижения автомобилями цели.
Целями автомобилей при многополосном движении могут быть, например,
съезд с дороги или поворот в определенную сторону на светофоре. В обоих
случаях машины, начиная с определенного момента времени, стремятся
перестроиться в целевую полосу, игнорируя значения плотности и скорости на
ней. Таким образом, в описываемой модели для каждой машины необходимо
хранить параметр цели. Некоторые автомобили не достигают цели, если не
смогли вовремя перестроиться в полосу, из которой разрешен поворот или
съезд. Это происходит из-за того, что нужная полоса занята другими
транспортными средствами. Модель предусматривает гибкий подход к
достижению автомобилем его цели: число автомобилей, которые не достигли
цели, может уменьшаться за счет разрешения поворотов из большего числа
полос, усложнения алгоритма перестроения или за счет реализации алгоритма
«вежливый водитель». Алгоритм заставляет водителей в целевой полосе
пропускать машины, которые хотят перестроиться.
Каждый автомобиль в модели имеет набор параметров, таких как
уникальный номер (ID), текущая скорость, максимальная скорость, конечная
цель. Имеются также вспомогательные параметры, упрощающие выполнение
логических операций. В модели, однако, мы оперируем именно с ячейками, а не
с автомобилями. Наличие автомобиля и набор его характеристик – это свойство
данной ячейки на данном шаге по времени. Все параметры, кроме скорости, в
пустых ячейках равны значениям по умолчанию, мы их не используем. Если
скорость в ячейке с автомобилем принимает значение от 0 до 4 и соответствует
скорости автомобиля, то параметр скорости пустой ячейки полагается равным –
1, ячейки с неподвижным препятствием -10. Эти значения используются для
построения алгоритмов движения.
Рассмотрим сначала самый простой и самый общий случай – движение по
участку дороги без перекрестков, светофоров, сужений/расширений и т.д. Итак,
на каждом шаге по времени обновление состояния ячеек в модели происходит в
два этапа:
А) Для каждого АТС выясняется возможность и необходимость смены
полосы. Производится смена полосы.
Б) Производится движение вперед по выбранной полосе по правилам
однополосного движения.
В самом общем виде алгоритм смены полосы А) можно записать так:
1) Проверка четности номера шага по времени. На четных шагах
происходит смена полосы направо, на нечетных – налево: если step%2=0
wmr=true, wml=false; иначе wmr=false, wml=true.
Здесь wml(r)=true/false – переменная, обозначающая, состоится или не
состоится смена полосы налево/направо соответственно. Изначально ее
значение – true (перестроение состоится). Дальнейшая последовательная
проверка условий переопределяет значение переменной (или может оставлять
без изменений).
2) Выгодно ли перестроение: если dm,k±1≤ dm,k или vm,k±1< vm,k то wml(r)=false
(здесь dm,k – расстояние до а/м в данной полосе, dm,k±1 – расстояние до а/м в
целевой полосе, vm,k и vm,k±1 – скорости этих автомобилей, m – номер
ячейки, k – номер полосы).
3) Есть ли препятствие на этой полосе: если существует m: i < m ≤ i+10,
vm,k=-10 то wml(r)=true (m – номер ячейки впереди, vm,k=-10 соответствует
неподвижному препятствию).
4) Свободна ли целевая ячейка: если vi,k±1≠ -1 то wml(r)=false (vi,k±1= -1
соответствует пустой ячейке)
5) Условие безопасности: если dprevm,k±1 < V то wml(r)=false (dprevm,k±1 -
расстояние до предыдущего автомобиля в целевой полосе, vprevi,k±1 - его
скорость, V – величина, выбираемая исходя из стиля вождения данного
водителя. V=vmax или V=vprevi,k±1 для «осторожных» или «агрессивных»
водителей соответственно, более подробно об этом будет рассказано
ниже).
6) Стохастический элемент: если p
vi,k,j=min(vi,k,j -1, 0), при di,k,j=0 то vi,k,j= 0 (reg – режим светофора, R
соответствует красному сигналу для данного направления)
4) Поворот/смена дороги j: если vi,k,j= 0 и reg ≠ R и k = kgoal и vi*,k*,j* =-1
(целевая ячейка свободна) то vi,k,j→ vi*,k*,j* (здесь kgoal – целевая полоса, из
которой разрешен поворот)
5) Движение прямо, если текущая полоса не предполагает поворот и целевая
ячейка свободна: если vi,k,j= 0 и reg ≠ R и k ≠ kgoal и vi+1,k,j =-1 то vi,k,j→ vi+1,k,j
Если цель машины поворот на перекрестке, и она находится в пределах
D2 клеток от поворота, то ее скорость уменьшается, пока не достигнет 1. При
достижении точки поворота машина останавливается.
Если горит красный сигнал светофора, автомобиль, подъезжающий к
перекрестку, замедляется, в ячейке перед перекрестком он останавливается.
Когда загорается зеленый, автомобиль снова начинает движение.
Если цель автомобиля – поворот на перекрестке (номер дороги,
соответствующий цели, предполагает поворот) и он находится в полосе, из
которой данный поворот разрешен, он поворачивает. Это означает, что он
попадает в соответствующую ячейку другой дороги, проверив, что она не
занята.
Алгоритм обеспечения безопасности на перекрестке
Рассмотрим реализацию поворотов более подробно. Расчетная область
всех элементов УДС состоит из четырех независимых участков. Для
четырехстороннего перекрестка эти участки – направления движения
(обозначены цифрами 0-3 на рисунке 3). Обозначения направлений следующие:
движение снизу вверх – j=0, сверху вниз – j=1, справа налево – j=2, слева
направо – j=3.
Каждое направление, или дорога, может иметь любое количество полос.
Рисунок 3 – схема движения на четырехстороннем перекрестке
Если моделируется, скажем, прямой участок УДС, то мы «выключаем»
лишние направления (просто взяв входящие потоки по данным направлениям
равными нулю) и взаимодействие этих направлений с актуальными.
При описании пересечений дорог различных типов необходимо,
разумеется, описать взаимодействие пересекающихся направлений. В данном
случае при расчете движения на четырехстороннем перекрестке у направлений
0-3, 0-2, 1-3, 1-2 имеются общие области. Они вместе с «разделительными
полосами» в одну клетку, имеющимися между направлениями 0-1 и 2-3,
образуют область перекрестка.
Для учета нахождения автомобилей на перекрестке создается матрица
cross[x][y]=vj1 + vj2, где vj1 и vj2 – значения скоростей в соответствующих
ячейках по направлениям j1 и j2, x и y – переменные, нумерующие элементы
матрицы перекрестка. Если ячейка расчетной области пустая, значение
скорости в ней равно -1. Таким образом, элемент матрицы равен -2 тогда и
только тогда, когда ячейка перекрестка пустая.
При подъезде автомобиля к границе перекрестка происходит обращение к
соответствующему элементу матрицы cross[x][y] для проверки наличия
автомобилей. Если целевая ячейка занята, автомобиль останавливается и ждет,
когда она освободится. На следующем шаге он может попробовать
перестроиться в другой ряд и повторить попытку проехать.
Ряд алгоритмов отвечает за описание особенностей поведения водителей.
Так, водители могут быть «осторожными», «агрессивными» и «вежливыми».
Характеристика «осторожный/агрессивный» отвечает за минимальное
расстояние в целевой полосе, необходимое для того, чтобы данный водитель
счел перестроение безопасным. Характеристика «вежливый» предназначена
для описания ситуаций, в которых водитель пропускает другие автомобили.
Рассмотрим данный алгоритм подробнее.
Алгоритм «вежливый водитель» для смены полосы
Алгоритм «вежливый водитель» (“cooperative driver”) состоит из двух
частей. Первая часть – модифицированный алгоритм смены полосы при
наличии в системе вежливых водителей (рисунок 4). Вторая часть – присвоение
водителю статуса «вежливый».
Вежливые водители в системе отмечены при помощи дополнительного
параметра, описывающего состояние каждой ячейки – cd, который принимает
значения да/нет (1 или 0). Напомним, что каждая ячейка имеет набор
параметров, позволяющий сохранить всю необходимую информацию об
автомобиле: его скорость, максимальная скорость, идентификационный номер
ID, точка назначения (номер выезда), время въезда в расчетную область, а
также ряд других вспомогательных параметров. Среди этих параметров есть
индикаторы желания перестроиться направо/налево (lchr(l), lane change
right/left), которые не учитывают условия безопасности, то есть, они
определяют не окончательное решение о перестроении, а лишь намерение
поменять полосу. Эти параметры в модели выполняют ту же роль, что и
«поворотники» в настоящем автомобиле. Другие водители могут их видеть и,
при определенных обстоятельствах, принимать во внимание и замедляться или
пропускать автомобиль, который хочет перестроиться, а могут и игнорировать.
Рисунок 4 – Алгоритм перестроения при наличии «вежливых» водителей
Первая часть алгоритма является алгоритмом смены полосы при наличии
препятствия. Отличие начинается после проверки, занимает ли препятствие
меньше полос с той стороны, куда предполагается перестроение. В данном
месте параметр wmr/wml сохраняется в переменную, обозначающую желание
перестроиться. Если в результате указанных действий параметр имеет значение
true, автомобиль «включает поворотник».
Далее происходит проверка безопасности перестроения. Она выглядит по-
разному для расчета с «вежливыми» водителями (cdonoff=1) и без них
(cdonoff=0). Если в системе нет «вежливых» водителей, то проверяется условие
dprev>V*, где dprev – расстояние до предыдущего автомобиля в целевой полосе, а
V* – скорость для сравнения, зависящая от особенностей стратегий водителей
(«агрессивных» и «осторожных»). При наличии в системе «вежливых»
водителей (cdonoff=1) условие безопасности dprev>V* уже может быть
избыточно жестким, так как вежливые водители специально тормозят и
останавливаются, чтобы пропустить других, и водитель с включенным
«поворотником» может безопасно перестроиться вплотную перед стоящим
автомобилем с «вежливым» водителем. Таким образом, совместно с условием
безопасности проверяется, является ли водитель предыдущего автомобиля в
целевой полосе «вежливым» (cd=1) на данный момент. О том, как водители в
системе становятся «вежливыми», будет изложено далее.
Соответственно, если условие безопасности не выполнено, смена полосы
отменяется.
Следующим шагом в алгоритме является рандомизация. Водитель,
желающий поменять полосу и имеющий такую возможность, тем не менее,
может передумать перестраиваться с некоторой заданной вероятностью p.
Затем, если водитель не передумал, происходит перестроение.
Далее рассмотрим, в каких условиях водителю присваивается статус
«вежливый». Это переменная величина, которая переопределяется каждый шаг
по времени, однако можно задать определенное правило, по которому только
часть водителей может стать вежливыми в процессе расчета. В нашем случае в
качестве определяющего параметра служит идентификационный номер (ID),
присваиваемый при добавлении автомобиля в систему. Половина водителей (с
четными ID) могут стать «вежливыми», а другая половина – нет.
Присвоение параметра cd=1 происходит непосредственно перед сменой
полосы, после выполнения всех шагов рис. 1. Таким образом, использоваться
этот параметр в алгоритме будет только на следующем шаге по времени. Для
случая перестроения условия для того, чтобы водитель стал «вежливым»,
следующие:
• step%2=0 для wmr (step%2=1 для wml)
• «Поворотник» автомобиля в полосе с препятствием включен: lchr(l)=true.
• Скорость автомобиля в полосе с препятствием vprevi,k±1=0. В полосе с
препятствием уже образовался затор.
• Скорость данного автомобиля -1< vi,k≤ 1. Он либо уже остановился, либо
его скорость достаточно мала для полной остановки.
• Номер данного автомобиля в системе ID%2=0 (или другое условие,
определяющее соотношение «вежливых» водителей и всех остальных),
• cd≠-1 (cd=-1 после того, как водитель пропустил автомобиль из соседнего
ряда и тот перестроился).
Последнее условие служит для того, чтобы один и тот же водитель, став
вежливым, не пропускал сразу много машин из соседнего ряда, а сам начинал
движение, пропустив только одну. Таким образом моделируется так
называемое перестроение «елочкой», когда водители из свободного ряда
пропускают по одному автомобилю из ряда с затрудненным движением.
Практика показывает, что такой подход проезда участка с неподвижным
препятствием наиболее эффективен.
Все созданные алгоритмы верифицированы с помощью тестовых задач,
показывающих адекватность их работы.
Третья глава посвящена численной реализации созданной модели и
алгоритмов. Программный комплекс, получивший название CAM-2D (Cellular
Automata 2-Dimensional), написан на языке С/С++ и использует библиотеку Glut
для визуализации расчетов. Визуализация при помощи Glut происходит
одновременно с расчетом и позволяет отлаживать программу и анализировать
результаты. Пример такой визуализации представлен, например, на рисунке 3.
Программный комплекс имеет модульную структуру, каждый модуль
отвечает за расчет одного элемента УДС, реализуя соответствующие
алгоритмы, представленные в Главе 2: прямолинейный участок, перекресток
(включая различные типы и настройки светофоров), участок с
расширением/сужением, въезд/выезд с автомагистрали, участок с разворотом,
участок с пешеходным светофором и т. д.
Для расчетов на дорожных сетях была создана параллельная версия
программы. Так как логические операции плохо поддаются распараллеливанию
на графических ускорителях (GPU), был выбран вариант реализации
параллельной версии на CPU. Параллельная версия программного комплекса
CAM-2D написана на языке С/С++. Отладка и тестовые расчеты выполнены на
суперкомпьютере К-100 ЦКП ИПМ им. М.В. Келдыша РАН.
Созданной модели транспортных потоков присущ как геометрический
параллелизм, так и параллелизм данных. Для реализации был выбран первый
вариант. Расчетная область – часть городской улично-дорожной сети –
разделяется на отдельные элементы: перекрестки, прямые участки дороги,
участки с въездами и съездами и т.д. Границы проходят по участкам без
особенностей, обеспечивая простоту передачи данных.
Из модулей, как из конструктора, можно собирать любые сети,
соответствующие реальной карте города.
Для того чтобы программный комплекс был универсальным, модули,
отвечающие за моделирование на отдельных участках УДС, стандартизированы
по размеру. Каждый модуль рассчитывается отдельно на своем процессоре, на
границах происходит обмен данными при помощи библиотеки MPI.
Автомобиль, выезжающий с одного элемента, появляется на другом, в месте,
соответствующем его скорости на момент перехода.
В отличие от многих других задач, при моделировании с помощью
клеточных автоматов недостаточно просто передать данные с одного
процессора на другой. Перед передачей нужно убедиться, что автомобили,
которые покидают один элемент, действительно могут попасть на следующий:
их целевая ячейка, соответствующая текущей скорости, также как и все
предшествующие ячейки, должна быть свободна.
Для удобства использования программного комплекса совместно с
системными программистами был разработан модуль пользовательского
интерфейса и визуализации данных. Веб-интерфейс позволяет задавать
дорожную сеть, вводить начальные данные и параметры моделирования, а
также получать результаты расчета, как в виде численных характеристик, так и
в виде визуализации (дорог со светофорами и движущимися по ним
автомобилями). Пример визуализации расчета для сети, состоящей из двух
перекрестков (Т-образного и четырехстороннего) представлен на рисунке 5.
Данный расчет выполнен при помощи параллельной версии программного
комплекса.
Рисунок 5 – результаты расчета на сети. Т-образный + четырехсторонний
перекрестки
Четвертая глава посвящена верификации созданной модели, алгоритмов
и программного комплекса. Для этого решается ряд тестовых задач и
проводится сравнение полученных результатов с экспериментальными
данными, результатами других исследователей и результатами, полученными
при помощи коммерческого пакета для транспортного моделирования Aimsun
TSS.
В качестве первой из тестовых задач была выбрана задача о пропускной
способности регулируемого четырехстороннего перекрёстка. Целью этой
задачи является поиск оптимального светофорного режима, а именно
определениедлительностисигналовсветофоров,обеспечивающей
минимальное время прохождения перекрёстка всеми участниками движения.
Перекрёсток состоит из двух пересекающихся дорог, по которым
автомобили движутся в двух направлениях и каждое направление имеет по две
полосы. Светофор на перекрестке имеет 4 фазы работы, которые включаются
последовательно. Меняя длительности фаз, при заданных потоках можно
получить разные значения пропускной способности. Результаты расчетов
представлены на рисунке 6.
По оси абсцисс отложены номера режимов светофоров, по оси ординат –
пропускная способность перекрестка в автомобилях в час. Верхняя кривая
(красные маркеры в виде квадратов) – расчет при помощи Aimsun, нижняя
(синие маркеры в виде ромбов) – расчет при помощи CAM-2D.
Из рисунка можно сделать вывод о том, какие режимы обеспечивают
наибольшую и наименьшую пропускную способность перекрестка.
Рисунок 6 – Пропускная способность перекрестка при различных режимах
работы светофора
Полученные результаты для микроскопической модели CAM-2D и для
Aimsun TSS качественно совпадают. Количественного совпадения мы не
ожидаем, этот вопрос решается калибровкой моделей.
Для сравнения с экспериментальными данными на рисунке 7 (слева)
представлен фрагмент пространственно-временной диаграммы реального
транспортного потока по автомагистрали, взятый из книги Б. Кернера4. По оси
абсцисс отложено время (ч), по оси ординат – расстояние вдоль трассы (км),
цветом представлена скорость автомобилей. Красные участки соответствуют
затору (~0-10 км/ч), желтые – синхронизированный поток (~60-70 км/ч), белые
участки – свободный поток (~110 км/ч).
Для сравнения выбрана задача о въезде на автомагистраль с переменным
потоком, она взята из книги М. Трайбера5, стр. 141. Результат моделирования,
взятый из данной книги, выполнен при помощи модели Кернера-Конхойзера и
Kerner B. The physics of Traffic – Berlin: Springer, 2004. – 682 p.
Treiber M., Kesting A. Traffic Flow Dynamics. Data, Models and Simulation –Berlin-Heidelberg: Springer,
2013. – 503 p.
представлен на рисунке 7 посередине. Результат моделирования при помощи
пакета CAM-2D представлен на рисунке 7 справа.
Сравнивая рисунки, можно сделать вывод, что, хотя качественно
результаты моделирования совпадают, однако, результат, полученный при
помощи СА-модели гораздо ближе к характерным экспериментально
наблюдаемым паттернам скоростей на пространственно-временной диаграмме.
Так как данная задача является модельной, и в реальных условиях потоки
меняются по-другому, мы не ожидаем полного аналога полученных при
помощи CAM-2D результатов эксперименту, где широкие движущиеся
кластеры появляются и исчезают с течением времени. Еще одна причина
отличия – участок реальной трассы содержит не один въезд, как в задаче, а
несколько въездов и выездов. В остальном, полученная при моделировании
картина отвечает всем ожидаемым критериям соответствия эксперименту.
Рисунок 7 – Экспериментальная пространственно-временная диаграмма
транспортного потока (слева), результат моделирования при помощи модели
Кернера-Конхойзера (посередине), результат моделирования при помощи
CAM-2D (справа).
В заключении приводятся основные результаты диссертационной работы.
В рамках диссертационной работы было выполнено исследование истории
развития транспортного моделирования и теории клеточных автоматов,
проведен анализ методов и подходов к моделированию транспортных потоков,
выполнено обоснование выбора направления исследования.
Были разработаны:
Двумерная многополосная микроскопическая модель транспортных потоков
на основе теории клеточных автоматов CAM-2D,
Алгоритмы движения на различных элементах улично-дорожной сети для
созданной модели
Различные поведенческие алгоритмы водителей для созданной модели
Комплекс программ для моделирования движения автотранспорта на
городскойулично-дорожнойсети,адаптированныйдля
высокопроизводительных вычислительных систем, интегрированный с
модулем пользовательского интерфейса и визуализации результатов
расчетов.
Были проведены:
Вычислительные эксперименты, верифицирующие созданные модель,
алгоритмы и программный комплекс.
Результаты, представленные в диссертации, были опубликованы в работах
[1-17].
Обсуждаются дальнейшие перспективы развития созданных модели,
алгоритмов и программного комплекса. В будущем планируется дальнейшее
усложнение созданной модели, например, адаптация ее для описания движения
грузового и общественного автотранспорта. Планируется при помощи
пользовательского интерфейса реализовать интерактивное управление
непосредственно в процессе расчета: смену светофорных режимов, добавление
неподвижных препятствий и так далее. Планируется провести исследование
оптимизации загрузки процессоров при параллельных расчетах.
Актуальность
Высокая загруженность автомобильных дорог транспортом – проблема, с
которой сталкивается каждый житель большого города и его пригородов.
Маятниковая миграция населения, вызванная сосредоточением рабочих мест
в крупных городах, вместе с высокой стоимостью жилья, переносит
проблемы заторов далеко за пределы мегаполиса. Растущее число
автомобилей в расчете на душу населения только усугубляет эту проблему.
Для улучшения транспортной ситуации требуется инструмент, позволяющий
как управлять существующей транспортной инфраструктурой, так и
планировать строительство новой, а также прогнозировать эффективность
ограничительных мер, которые обычно бывают непопулярными (например,
платная парковка на улицах города, ограничение въезда и др.)
Во многих мегаполисах мира успешно работает так называемая
интеллектуальная транспортная система (ИТС), позволяющая в реальном
времени отслеживать и вносить изменения в управление дорожным
движением, выбирая самые эффективные решения.
В связи с появлением более мощных вычислительных ресурсов появилась
возможность строить более сложные и приближенные к жизни
математические модели транспортных потоков, учитывающие
взаимодействие водителей с дорожными знаками, разметкой, светофорами и
так далее, а также между собой. Моделирование на высокопроизводительных
системах может позволить охватывать дорожные сети в масштабах целого
города.
Появление в последние годы множества новых данных о потоках
автотранспорта, собираемых датчиками и камерами видеонаблюдения, во-
первых, открывает новые возможности для изучения закономерностей
динамики потоков, а во-вторых, позволяет проводить точную калибровку
моделей для использования в конкретном месте, что может улучшить
качество прогнозов, полученных в рамках вычислительных экспериментов.
Это особенно важно для управления потоками на «микроуровне» – то есть в
масштабе нескольких соседних перекрестков или кварталов – например,
вблизи мест скопления людей, таких как торговые центры, концертные залы
и пр.
Все вышеперечисленное позволяет сделать вывод об актуальности и
перспективности выбранного автором диссертации направления
исследований.
Современное состояние исследований
Несмотря на более чем 70-летнюю историю развития моделирования
движения автотранспорта, данная тема по-прежнему остается
малоизученной. Первые модели, созданные в рамках макроскопического
подхода, не обеспечивают необходимую точность и не воспроизводят многие
свойства реальных транспортных потоков, которые наблюдаются
экспериментально. Связано это, в первую очередь, с вынужденной
экономией на вычислениях. По мере развития этой области моделирования,
развивались и усложнялись применяемые методы и подходы одновременно с
ростом производительности компьютеров. С этим же связан и переход от
макроскопических моделей, описывающих транспортный поток как
сплошную среду, к микроскопическим моделям, которые отображают более
тонкие взаимодействия, описывая движение каждого отдельного автомобиля,
и дальнейшее усложнение этих моделей. Этот процесс идет до сих пор.
В целом, из-за большого разнообразия методов и подходов, применяемых в
транспортном моделировании, а также большого разнообразия задач,
которые необходимо решить для реального моделирования трафика на
городской дорожной сети, помимо непосредственно описания движения
автомобилей, область по-прежнему остается недостаточно исследованной.
Выбранный для диссертационной работы подход на основе теории
клеточных автоматов берет свое начало со статьи Нагеля и Шрекенберга [1]
1992 года. За прошедшие годы выделилось два основных направления
развития данного подхода: моделирование на городских сетях при помощи
сложных структур на основе ячеек клеточных автоматов [2, 3] и
исследования в рамках теории трех фаз Б. Кернера [4], дополняющие правила
перемещения вдоль дороги из модели Нагеля-Шрекенберга [5-10].
Подход на основе теории клеточных автоматов очень удобен для описания
различных поведенческих стратегий водителей, включающих ситуации
дискретного выбора, поскольку основан не на уравнениях в частных
производных, а на логических алгоритмах и условиях «если – то».
Моделирование водительских стратегий и принятия решений во время
движения в потоке – область в настоящее время мало исследованная и
потому представляющая интерес.
Также до сих пор при решении задач транспортного моделирования редко
встречается выполнение расчетов на высокопроизводительных системах, что,
безусловно, является важным для масштабных расчетов на городских сетях.
Цель и задачи диссертационной работы
Целью работы является создание новой математической модели
транспортных потоков, учитывающей реальные особенности поведения
водителей и их взаимодействия с дорожно-транспортной инфраструктурой, и
комплекса программ, адаптированного для высокопроизводительных
вычислительных систем, реализующего эту модель.
Для этого автором диссертации решены следующие задачи:
Построена двумерная (многополосная) модель, включающая алгоритмы для
движения автотранспорта на различных элементах улично-дорожной сети:
регулируемых и нерегулируемых перекрестках, участках дороги с
расширением и сужением, прямых участках дороги и т.д. В модели учтены
разные стили вождения и стратегии смены полос участников движения,
рассмотрен объезд препятствий (например, автомобилей, попавших в ДТП).
Разработан модульный программный комплекс, позволяющий делать
расчеты на указанных элементах УДС по отдельности (в последовательном
режиме), а также собирать из данных элементов сети различной
конфигурации и делать расчеты на них (в параллельном режиме).
Проведено тестирование алгоритмов и верификация модели и комплекса на
тестовых задачах. Проведено сравнение полученных результатов с
экспериментальными данными, взятыми из литературы, а также с расчетами,
выполненными при помощи коммерческого пакета для моделирования
движения автотранспорта Aimsun TSS [11].
Научная новизна
Работа сочетает в себе два свойства, редко встречающихся в работах других
исследователей (даже по отдельности): во-первых, учет при моделировании
стратегий вождения, а во-вторых – адаптацию модели для расчетов на
суперкомпьютерах. В этом основное отличие представленных модели,
алгоритмов и программного комплекса от других.
Созданная модель на основе теории клеточных автоматов является новой и
оригинальной. Входящие в модель алгоритмы проезда различных элементов
УДС, а также алгоритмы взаимодействия водителей с дорожной
инфраструктурой и с другими водителями, алгоритмы, реализующие
различные стратегии поведения на дороге, и программный комплекс для
высокопроизводительных вычислительных систем, реализующий модель,
также являются новыми и оригинальными.
Теоретическая и практическая значимость работы
В теоретической части, созданные новые алгоритмы, учитывающие разные
стратегии поведения в транспортном потоке, а также реализация этих
алгоритмов в виде комплекса программ для высокопроизводительных
систем, представляет интерес как основа для создания новых подходов к
моделированию автотранспорта.
В практической части, созданная модель и программный комплекс могут
служить базой для интеллектуальных транспортных систем городов,
позволяя управлять транспортными потоками в режиме реального времени, а
также решать прогнозные задачи при планировании новой и модификации
существующей дорожно-транспортной инфраструктуры.
Методика исследования
Подход, основанный на теории клеточных автоматов, максимально хорошо
подходит для описания сложных разнородных систем, какими являются
транспортные потоки. Одно из преимуществ подхода – его дискретная
логическая основа, которая позволяет легко описывать принятие решений
водителями при движении и взаимодействии с другими участниками
дорожного движения.
Созданная модель основана на одномерной модели Нагеля-Шрекенберга.
Модель обобщена на многополосный случай путем добавления набора
правил для смены полосы в различных условиях. Вместо одномерной
расчетной области имеем двумерную; автомобили могут перестраиваться и
обгонять друг друга. Правила обновления состояния ячеек из модели Нагеля-
Шрекенберга, состоящие из четырех шагов, описывающих движение вдоль
дороги, дополнены разветвленными логическими алгоритмами для
перестроений. Для моделирования движения в городе в модель добавлены
элементы управляющего воздействия, такие, как дорожные знаки, сигналы
светофора и разметка.
Поведение водителей на дороге может существенно отличаться в
зависимости от их личностных качеств, стажа вождения и других
параметров. Учет этих стратегий облегчает калибровку модели при
практическом моделировании на участках дорожных сетей реальных городов
(например, в больших городах стиль вождения, как правило, носит более
агрессивный характер, а в городах с неинтенсивным движением – более
спокойный). Также, как известно, частая смена полосы отдельными
водителями может существенно влиять на весь транспортный поток,
порождая волны уплотнений при стесненном движении, даже при отсутствии
физических препятствий на дороге.
Для учета разной манеры вождения в модель включены «агрессивные» и
«осторожные» водители, а для реалистичного моделирования проезда
участка с сужением, при котором обычно водители пропускают автомобили
из ряда с препятствием по одному (“zipper merge” или «елочка») – алгоритм
«вежливый водитель».
Для верификации созданной модели и алгоритмов решен ряд тестовых задач,
выполнено сравнение с экспериментальными данными, данными,
полученными с использованием других моделей (взятыми из литературы) и
результатами моделирования при помощи коммерческого пакета для
транспортного моделирования. Выполнены расчеты на небольших дорожных
сетях в параллельном режиме на суперкомпьютере К-100 (ИПМ им.
Келдыша РАН) [12].
Достоверность полученных результатов
Тестовые расчеты согласуются с экспериментальными данными. Сравнение с
расчетами, полученными другими исследователями и с помощью
коммерческого пакета Aimsun TSS, демонстрируют сходство основных
закономерностей, однако, обеспечивают более детализированные данные и
большее сходство с экспериментом.
Апробация работы, внедрение
Работа была апробирована на ряде научных конференций и семинаров:
1. “Different Approaches to the Multilane Traffic Flows Simulation”. 17th
International Conference Mathematical Modelling and Analysis, June 6 – 9,
2012, Tallinn, Estonia.
2. «Макро- и микроскопические модели транспортных потоков на
многополосных магистралях». Всероссийская конференция
«Теоретические основы и конструирование численных алгоритмов
решения задач математической физики», посвященной памяти К.И.
Бабенко, 10-16 сентября 2012 г., Россия, п. Абрау-Дюрсо.
3. “Different Approaches to the Multilane Traffic Simulation”. 15th International
Conference “Traffic and Granular Flow -13” (TGF-13), 25-27 September,
2013, Juelich Germany.
4. “Mathematical simulation of multilane traffic flows using macro- and
microscopic approaches” 6th Int. Conf. on Computational Methods in Applied
Mathematics (CMAM-6), Sept. 28 – Oct. 4, 2014, St. Wolfgang, Austria.
5. “Two approaches to the vehicular traffic flows simulation using high-
performance computer systems”. International Interdisciplinary Seminar
“Mathematical Models and Modeling in Laser-Plasma Processes and Advanced
Science Technologies”, May 30 – June 6, 2015, Petrovac, Montenegro.
6. “Simulation of multilane vehicular traffic on the basis of cellular automata
theory”, International Interdisciplinary Seminar “Mathematical Models and
Modeling in Laser-Plasma Processes and Advanced Science Technologies”
July 4-9, 2016, Moscow, Russia.
7. “Multilane Traffic Flow Modeling Based on Cellular Automata Theory Using
High-performance Computer Systems”. 5th European Seminar on Computing
(ESCO 2016), June 5-10, 2016, Pilsen, Czech Republic.
8. “Multilane traffic flow modeling using cellular automata theory”. International
Conference on Mathematical Modeling and Computational Physics
(MMCP’2017), July 3-7, 2017, Dubna, Russia.
9. “Traffic flow modelling on road networks using cellular automata theory”,
International Conference on Communication, Management and Information
Technology (ICCMIT2018), April 02-04, 2018, Madrid, Spain.
10. “Cellular automata in application to traffic flow simulation”, Актуальные
проблемы прикладной математики и информационных технологий – Аль-
Хорезми 2018, 13–15 сентября 2018 г., Ташкент, Узбекистан.
11. «Моделирование и визуализация транспортных потоков на элементах
дорожной сети с использованием теории клеточных автоматов», ХII
мультиконференция по проблемам управления (МКПУ-2019), 23-28
сентября 2019, Дивноморское, Геленджик, Россия.
12. “Simulation and visualization of vehicular traffic on road networks using high
performance computing systems”. 19th International Conference of
Computational and Mathematical Methods in Science and Engineering, June
30-July 6, 2019, Rota, Cadiz, Spain.
13. “Reproduction of experimental spatio-temporal structures in traffic flows using
mathematical model based on cellular automata theory”, International
Conference on Communication, Management and Information Technology
(ICCMIT2019), March 26-28, 2019, Vienna, Austria.
14. “Supercomputer Technology for Traffic Simulation in a Metropolis”. 2020
International Conference on Engineering Management of Communication and
Technology (EMCTECH), Oct. 20-22, 2020, Vienna, Austria.
15. «Модель транспортных потоков на основе теории клеточных автоматов
для решения задач управления движением на городских дорожных сетях»
13 Конференция по проблемам управления МКПУ-2020, Конференция
«Математическая теория управления и ее приложения» (МТУиП-2020). 6-
8 октября 2020г., Санкт-Петербург, Россия.
16. «Алгоритмы поведения водителей в математической модели
транспортных потоков на основе теории клеточных автоматов» V
Международная конференция Моделирование нелинейных процессов и
систем (MNPS-2020), 16-20 ноября 2020 г., Москва, Россия.
17. «Алгоритмы поведения водителей в математической модели
транспортных потоков на основе теории клеточных автоматов», 79-я
международная научно-методическая и научно-исследовательская
конференция МАДИ, 25-30 января 2021г., Москва, Россия.
Гранты РФФИ, включающие результаты работы:
13-01-00781-а (2013-2015)
Математическое моделирование потоков автотранспорта с использованием
высокопроизводительных вычислительных систем на основе методов
математической физики и нейросетевой технологии, руководитель –
Чурбанова Н. Г.
13-01-12008-офи_м (2013-2015)
Различные подходы к математическому моделированию многополосных
транспортных потоков на высокопроизводительных многопроцессорных
вычислительных системах, руководитель – Чурбанова Н. Г.
13-01-12046-офи_м (2013-2015)
Разработка фундаментальных основ, методов и средств прикладного
моделирования транспортных процессов и систем на вычислительных
комплексах гибридной архитектуры, руководитель – Колесниченко А. В.
14-01-90009-Бел_а (2014-2015)
Разработка и исследование сеточных алгоритмов для эффективных расчётов
некоторых эволюционных задач математической физики, гидродинамики и
динамики транспортных потоков, руководитель – Четверушкин Б. Н.
15-01-03445-а (2015-2017)
Разработка эффективных вычислительных технологий для решения
крупномасштабных задач математической физики на многопроцессорных
системах сверхвысокой производительности, руководитель – Кулешов А. А.
16-01-00347-а (2016-2018)
Моделирование потоков автотранспорта на больших транспортных сетях с
использованием микроскопического и мультиагентного подходов,
руководитель – Чурбанова Н. Г.
16-31-00087-мол-а (2016-2017)
Математическое моделирование многополосных транспортных потоков на
основе теории клеточных автоматов с использованием
высокопроизводительных вычислительных систем, руководитель – Чечина А.
А.
18-01-00405-а (2018-2020)
Развитие современных вычислительных технологий для моделирования
крупномасштабных физических процессов с использованием
высокопроизводительных вычислительных систем, руководитель – Кулешов
А. А.
18-51-41001-Узб_т (2018-2019)
Высокоточные вычислительные алгоритмы для решения задач газовой
динамики, фильтрации и динамики транспортных потоков с использованием
современных суперкомпьютеров, руководитель – Якобовский М. В.
20-07-00528-а (2020-2022)
Разработка вычислительных технологий для исследования и визуальной
интерпретации динамики транспортных потоков на улично-дорожной сети
мегаполиса на базе систем сверхвысокой производительности, руководитель
– Чурбанова Н. Г.
Созданный программный комплекс CAM-2D был внедрен при создании
платформы поддержки принятия решений для проекта «Разработка
архитектуры, принципов построения и программного обеспечения
суперкомпьютерной информационно-аналитической платформы
прикладного моделирования, прогнозирования и экспертиз транспортных
процессов и систем в виде сетевой компьютерной лаборатории» в рамках
ФЦП «Исследования и разработки по приоритетным направлениям развития
научно-технологического комплекса России на 2014-2020 годы»,
Соглашение о предоставлении субсидии от 30.06.2014 г. № 14.604.21.0052.
(см. Приложение 1).
На программный комплекс CAM-2D (упрощенная версия) было получено
свидетельство на программу для электронных вычислительных машин №
2016662572 от 15.11.2016 (см. Приложение 2).
Основные положения, выносимые на защиту
1. Двумерная многополосная микроскопическая модель транспортных
потоков на основе теории клеточных автоматов
2. Алгоритмы движения на различных элементах улично-дорожной сети
для созданной модели
3. Различные поведенческие алгоритмы водителей для созданной модели
4. Комплекс программ CAM-2D для моделирования движения
автотранспорта на городской улично-дорожной сети, адаптированный
для высокопроизводительных вычислительных систем,
интегрированный с модулем пользовательского интерфейса и
визуализации результатов
5. Результаты вычислительных экспериментов, полученные с помощью
комплекса программ.
Объем и структура диссертации
Объем диссертационной работы – 136 страниц машинописного текста.
Диссертация состоит из введения, четырех глав, заключения,
библиографического списка из 90 наименований и двух приложений. Работа
содержит 50 рисунков. Глава 1 содержит обзор основных существующих
методов и подходов к транспортному моделированию, главы 2-4 содержат
собственные исследования автора диссертации.
В рамках диссертационной работы было выполнено исследование истории
развития транспортного моделирования и теории клеточных автоматов,
проведен анализ методов и подходов к моделированию транспортных
потоков, выполнено обоснование выбора направления исследования.
Были разработаны:
1. Двумерная многополосная микроскопическая модель транспортных
потоков на основе теории клеточных автоматов CAM-2D,
2. Алгоритмы движения на различных элементах улично-дорожной сети для
созданной модели
3. Различные поведенческие алгоритмы водителей для созданной модели
4. Комплекс программ для моделирования движения автотранспорта на
городской улично-дорожной сети, адаптированный для
высокопроизводительных вычислительных систем, интегрированный с
модулем пользовательского интерфейса и визуализации результатов
расчетов.
Были проведены:
5. Вычислительные эксперименты, верифицирующие созданные модель,
алгоритмы и программный комплекс.
Результаты, представленные в Главах 2-4, были опубликованы в работах [73-
90].
В будущем планируется дальнейшее усложнение созданной модели,
например, адаптация ее для описания движения грузового и общественного
автотранспорта. Планируется при помощи пользовательского интерфейса
реализовать интерактивное управление непосредственно в процессе расчета:
смену светофорных режимов, добавление неподвижных препятствий и так
далее. Планируется провести исследование оптимизации загрузки
процессоров при параллельных расчетах, поскольку на данном этапе загрузка
каждого процессора недостаточна для высокоэффективных вычислений.
Сложность данной задачи в том, что расчетные модули должны быть
универсальными, а чем больше и сложнее каждый отдельный модуль, тем
больше вариантов модулей нужно создать, чтобы моделировать движение на
городской дорожной сети. Для решения этой проблемы необходимо найти
«золотую середину» между размером модуля (и, соответственно, более
оптимальной загрузкой вычислительного узла) и его универсальностью.
Один из возможных путей решения – анализ конфигураций связанных
элементов УДС, наиболее часто встречающихся на карте города.
Созданная модель, алгоритмы и комплекс программ могут служить основой
для создания интеллектуальной транспортной системы города, позволяя
решать как краткосрочные, так и долгосрочные задачи транспортного
управления и планирования.
1. Nagel K., Schreckenberg M. A Cellular automaton model for freeway traffic //
J. Phys. I France. – 1992. – V. 2. – P. 2221-2229.
2. Buslaev A.P., Yashina M.V., Strusinskiy P.M., Sokolov P.A. On Qualitative
Properties of Cluster Model for Flows on Regular Networks // Composites,
New Challenges for the Greening of Transport, Springer-ECCOMAS, CM3
ECCOMAS Thematic Conference. – 2015. – P. 24-27.
3. Buslaev A.P., Tatashev A.G., Yashina M.V. On cellular automata, traffic and
dynamical systems in graphs // Int. J. of Eng. & Techn. – 2018. – V.7 (2.28) –
P. 351–356.
4. Kerner B. The physics of Traffic – Berlin: Springer, 2004. – 682 p.
5. Kerner B., Klenov S., Schreckenberg M. Simple cellular automaton model for
traffic breakdown, highway capacity, and synchronized flow // Phys. Rev. E 84
– 2011. – P. 046110.
6. Kerner B., Klenov S., Hermanns G., Schreckenberg M. Effect of driver over-
acceleration on traffic breakdown in three-phase cellular automaton traffic flow
models // Physica A: Statistical Mechanics and its Applications. – 2013. –
V.392 (18) – P. 4083-4105.
7. Vranken T.P.E., M. Schreckenberg M. Cellular Automata Intersection Model //
Collective Dynamics. – 2020. – V.5 A80 – P. 1-25.
8. Zhao H.T., Liu X.R., Chen X.X., Lu J.C. Cellular automata model for traffic
flow at intersections in internet of vehicles // Physica A: Statistical Mechanics
and its Applications. – 2018. – V.494 – P. 40-51.
9. Zhao H.T., Yang S., Chen X.X. Cellular automata model for urban road traffic
flow considering pedestrian crossing street // Physica A: Statistical Mechanics
and its Applications. – 2016. – V.462 – P. 1301-1313.
10. Gao K., Jiang R., Wang B.-H., Wu Q.-S. Discontinuous transition from free
flow to synchronized flow induced by short-range interaction between vehicles
in a three-phase traffic flow model // Physica A: Statistical Mechanics and its
Applications. – 2009. – V.388 (15–16) – P. 3233-3243.
11. https://www.aimsun.com
12. https://www.kiam.ru/MVS/resourses/k100.html
13. Lighthill M.H., Witham G.B. On kinematic waves: A theory of traffic flow on
long crowded roads // Proc. Royal Soc. Ser. A. – 1955. – V.229 – P. 317–345.
14. Kerner B., Konhäuser P. Cluster effect in initially homogeneus traffic flow //
Physical review. E, Statistical physics, plasmas, fluids, and related
interdisciplinary topics. – 1993. – V.48 – P. R2335-R2338.
15. Kerner, B., Konhäuser, P., Schilke, M. Deterministic spontaneous appearance
of traffic jams in slightly inhomogeneous traffic flow // Physical review. E,
Statistical physics, plasmas, fluids, and related interdisciplinary topics. – 1995.
– V. 51 – P. 6243-6246.
16. Сухинова А. Б., Трапезникова М. А., Четверушкин Б. Н., Чурбанова Н. Г.
Двумерная макроскопическая модель транспортных потоков // Матем.
моделирование. – 2009. – т. 21 №2 – С.118–126.
17. Prigogine I., Andrews F. C. A boltzmann-like approach for traffic flow // Opns
Res. – 1960. – Vol. 8 – P. 789–797.
18. Paveri-Fontana S.L. On Boltzmann like treatments for traffic flow. //
Transportation Research. – 1975. – V.9 No 225.
19. Helbing D. Gas-kinetic derivation of Navier-Stokes-like traffic equations. //
Phys. Rev. E. – 1996. – V.53 No 2366.
20. Helbing D. Structure and instability of high-density equations for traffic flow. //
Phys. Rev. E. – 1998. – V.57 No 6176.
21. Bando, M., Hasebe, K., Nakayama, A., Shibata, A. and Sugiyama, Y. Structure
Stability of Congestion in Traffic Dynamics // Japanese Journal of Industrial
and Applied Mathematics. – 1994. – V. 11 – P. 203-223.
22. Bando, M., Hasebe, K., Nakanishi, K., Nakayama, A., Shibata, A., Sugiyama,
Y. Phenomenological Study of Dynamical Model of Traffic Flow // Journal of
Physics I France. – 1995. – V. 5 – P. 1389-1399.
23. Treiber M., Kesting A. Traffic Flow Dynamics. Data, Models and Simulation –
Berlin-Heidelberg: Springer, 2013. – 503 p.
24. Treiber M., Hennecke A., Helbing D. Congested traffic states in empirical
observations and microscopic simulations // Physical Review E. – 2000. –
V. 62 (2) – P. 1805–1824.
25. Von Neumann J. The general and logical theory of automata // Cerebral
Mechanisms in Behavior. – 1948. – P. 1–41.
26. Фон Нейман Дж. Теория самовоспроизводящихся автоматов – М.: УРСС,
2010. – 384 c.
27. Gardner M. Mathematical Games – The Fantastic Combinations of John
Conway’s New Solitaire Game ‘Life’ // Scientific American. – 1970. – V. 223 –
P. 120–123.
28. S. Wolfram. A new kind of science. https://www.wolframscience.com/nks/
29. http://atlas.wolfram.com/
30. Бугаев А. С., Буслаев А. П., Козлов В. В., Таташев А. Г., Яшина М. В.
Обобщенная транспортно-логистическая модель как класс динамических
систем // Матем. моделирование. – 2015. – т. 27 №12. – С. 65–87.
31. Maerivoet S., De Moor B. Cellular automata models of road traffic // Physics
Reports. – 2016. – V. 419 – P. 1-64.
32. Li X., Sun J.-Q. Effects of turning and through lane sharing on traffic
performance at intersections // Physica A: Statistical Mechanics and its
Applications. – 2016. – V. 444 – P. 622-640.
33. Vasic J., Heather J., Ruskin H.J. Cellular automata simulation of traffic
including cars and bicycles // Physica A: Statistical Mechanics and its
Applications. – 2012. – V. 391 (8) – P. 2720-2729.
34. Lárraga M.E., Alvarez-Icaza L. Cellular automaton model for traffic flow
based on safe driving policies and human reactions // Physica A: Statistical
Mechanics and its Applications. – 2010. – V. 389 (23) – P. 5425-5438.
35. Combinido J.S.L., Lim M.T. Modeling U-turn traffic flow // Physica A:
Statistical Mechanics and its Applications. – 2010. – V. 389 (17) – P. 3640-
3647.
36. Zhu H.B. Numerical study of urban traffic flow with dedicated bus lane and
intermittent bus lane // Physica A: Statistical Mechanics and its Applications. –
2010. – V. 389 (16) – P. 3134-3139.
37. Qian W.L., Wang B., Lin K., Machado R.F., Hama Y. A mesoscopic approach
on atability and phase transition between different traffic flow states //
International Journal Of Non-Linear Mechanics. – 2017. – V. 89 – P. 59-68.
38. Майоров Н.Н., Романек В.А. Вопросы выбора математических моделей
для исследования пассажирских потоков в транспортных системах //
Системный Анализ И Логистика. – 2017. – т. 1 №14 – C. 39-45.
39. Мазурин Д.С. Комплексная методика моделирования транспортных
систем на основе моделей динамического распределения потоков // Труды
Института Системного Анализа Российской Академии Наук. – 2017. – т.
67 №3 – C. 3-12.
40. Кармадонов В.Ю., Косов А.А. Моделирование транспортного потока на
основевиртуальнойвычислительнойсреды//Информационные
Технологии И Проблемы Математического Моделирования Сложных
Систем. – 2017. – № 18 – C. 80-87.
41. Сорокин А.А., Яковлев С.В., Орлова А.Ю. Имитационная модель
транспортного узла с режимом оптимизации параметров // Вестник
Севкавгти. – 2017. – т. 4 №31 – C. 246-250.
42. Разумов Д.С., Катаев М.Ю., Шелестов А.А. Алгоритмы управления
потоком автотранспорта в городских условиях // Вестник Современных
Исследований. – 2018. – т. 6.1 №21 – C. 471-473.
43. Нургалиев Е.Р. Имитационное моделирование улично-дорожной сети
города // Актуальные Направления Научных Исследований XXI Века:
Теория и Практика. – 2017. – т.5 №7-1 (33-1) – C. 172-176.
44. Zhou J., Zhang H.L., Wang C.P., Shi Z.K. A new lattice model for single-lane
traffic flow with the consideration of driver’s memory during a period of time.
// International Journal Of Modern Physics C. – 2017. – V. 28 No 7. – P.
1750086
45. Jin D., Zhou J., Zhang H.L., Wang C.P., Shi Z.K. Lattice hydrodynamic model
for traffic flow on curved road with passing. // Nonlinear Dynamics. – 2017. –
V. 89 No 1 – P. 107-124.
46. Kaur R. Sharma S. Analysis of driver’s characteristics on a curved road in
lattice model // Physica A: Statistical Mechanics And Its Applications. – 2017.
– V.471 – P. 59-67.
47. Qian Y., Zeng J., Wang N., Zhang J., Wang B. A traffic flow model
considering influence of car-following and its echo characteristics // Nonlinear
Dynamics. – 2017. – V.89 No 2 – P. 1099-1109.
48. КурцВ.В.,АнуфриевИ.Е.Модельавтомобильноготрафикас
запаздывающим аргументом – исследование устойчивости на кольце //
Математическое Моделирование. – 2017. – т. 29 № 4 – C. 88-100.
49. Kuang H., Xu Z.P., Li X.L., Lo S.M. An extended car following model
accounting for the honk effect and numerical tests // Nonlinear Dynamics. –
2017. – V. 87 No 1 – P. 149-157.
50. Wang P., Yu G., Wu X., Wang Y. Linear and nonlinear stability analysis of an
extended car following model considering pedestrians on adjacent lane //
Nonlinear Dynamics. – 2017. – V. 88 No 1 – P. 777-789
51. Qi L., Zheng Z., Gang L. A cellular automation model for ship traffic flow in
waterways // Physica A: Statistical Mechanics And Its Applications. – 2017. –
V. 471 – P. 705-717.
52. Chen J., Jiang R., Lin L. Assigning on ramp flows to maximize capacity of
highway with two on-ramps and one off-ramp in between // Physica A:
Statistical Mechanics And Its Applications. – 2017. – V.465 – P. 347-357
53. Qian Y.S., Feng X., Zeng J.W. A cellular automata traffic flow model for three
phase theory // Physica A: Statistical Mechanics And Its Applications. – 2017.
– V. 479 – P. 509-526.
54. Guzman H.A., Larraga M.E., Alvarez-Icaza L., Carvajal J. A cellular automata
model for traffic flow based on kinetics theory, vehicles capabilities and driver
reactions. Physica A: Statistical Mechanics And Its Applications. –2018. –
V.491 – P. 528-548.
55. Dong P., Wang X., Yun L., Fan H. Research on the characteristics of mixed
traffic flow based on an improved bicycle model // Simulation, Sage
Publications. – 2018. – V.94 N 5 – P. 451-462.
56. Zeng J.W., Qian Y.S., Wang H., Wei X.T., modeling and simulation of traffic
flow under different combination setting of taxi stop and bus stop // Modern
Physics Letters B. – 2018. – V. 32 N 25 – P. 1850301
57. Kaur R.1, Sharma S. Modeling and simulation of driver’s anticipation effect in
a two lane system on curved road with slope // Physica A: Statistical Mechanics
And Its Applications. –2018. – V. 499 – P. 110-120.
58. Su Z., Liu S., Deng W., Li W., Cai X., transportation dynamics on networks of
heterogeneous mobile agents // Physics Letters A. – 2019. – V. 523 – P. 1379-
1386.
59. Yao W., Jia N., Zhong S., Li L., best response game of traffic on road network
of non-signalized intersections // Physica A: Statistical Mechanics And Its
Applications. – 2018. – V. 490 – P. 386-401.
60. Herty M., Fazekas A., Visconti G., a two-dimensional data-driven model for
traffic flow on highways // Networks And Heterogeneous Media. – 2018. – V.
13 N 2 – P. 217-240.
61. СтепановЕ.П.Анализэффективностидемультиплексирования
транспортных потоков // Моделирование и анализ информационных
систем. –2019. – т.26 №1(79) – C. 170-190.
62. Морозов И. И., Гасников А. В., Тарасов В. Н., Холодов Я. А., Холодов
А. С. Численноеисследованиетранспортныхпотоковнаоснове
гидродинамическихмоделей//Компьютерныеисследованияи
моделирование. – 2011. – т.3:4 – C. 389–412.
63. ХолодовЯ.А.,АлексеенкоА. Е.,ВасильевМ. О.,Холодов
А. С. Построение математической модели дорожного перекрестка на
основе гидродинамического подхода // Компьютерные исследования и
моделирование. – 2014. – т.6:4 – C. 503–522.
64. Прокопцев Н. Г.,Алексеенко А. Е., Холодов Я. А., Использование
сверточныхнейронныхсетейдляпрогнозированияскоростей
транспортного потока на дорожном графе // Компьютерные исследования
и моделирование. – 2018. – т. 10:3 – C. 359–367.
65. Веренцов С. И., Магеррамов Э. А., Виноградов В. А., Гизатуллин Р. И.,
Алексеенко А. Е., Холодов Я. А. Байесовская вероятностная локализация
автономного транспортного средства путем ассимиляции сенсорных
данных и информации о дорожных знаках // Компьютерные исследования
и моделирование. – 2018. – т. 10:3 – C. 295–303.
66. Введение в математическое моделирование транспортных потоков. /
Гасников А.В., Кленов С.Л., Нурминский Е.А., Холодов Я.А., Шамрай
Н.Б.; Приложения: Бланк М.Л., Гасникова Е.В., Замятин А.А. и Малышев
В.А., Колесников А.В., Райгородский А.М; Под ред. А.В. Гасникова. —
М.:МФТИ, 2010. — 362 с.
67. Smirnova M., Bogdanova A., Smirnov N., Kiselev A., Nikitin, V., Manenkova
A. Multi-Lane Unsteady-State Traffic Flow Models // Journal of Mechatronics.
– 2014. – V.2.
68. Smirnov N., Kiselev A., Nikitin V., Silnikov M., Manenkova A. Hydrodynamic
traffic flow models and its application to studying traffic control effectiveness
// WSEAS Transactions on Fluid Mechanics. – 2014. – V. 9. – P. 178-186.
69. Barlovic R. e. a. Metastable states in cellular automata for traffic flow // The
European Physical Journal B-Condensed Matter and Complex Systems. – 1998.
– V. 5 No. 3 – P. 793–800.
70. Герман М.С., Ермаков А.В. Архитектура и функциональные возможности
системы визуализации потоков автотранспорта на элементах улично-
дорожной сети // Препринты ИПМ им. М.В.Келдыша. – 2017. – № 145 –
21 с.
71. Измайлова М. О., Рахманкулов И. Ш. Категория «средняя величина» и ее
методологическое значение в научном исследовании. – Казань: Изд.
Казанского Университета, 1982. – 143 с.
72. Helbing D., Hennecke A., Shvetsov V., Treiber M. MASTER: Macroscopic
traffic simulation based on a gas-kinetic, non-local traffic model //
Transportation Research Part B: Methodological. – 2001. – V. 35 – P. 183-211.
73. Churbanova N.G., Chechina A.A., Furmanov I. R., Trapeznikova M.A.
Microscopic Model for Simulation of Traffic Flows on Multilane Highways
and Crossroads // ECCOMAS 2012 – European Congress on Computational
Methods in Applied Sciences and Engineering, e-Book Full Papers. – 2012. –
P. 11 – 17.
74. Трапезникова М.А., Чечина А.А., Чурбанова Н.Г., Поляков Д.Б.
Математическое моделирование потоков автотранспорта на основе макро-
и микроскопических подходов // Вестник АГТУ Сер.: Управление,
вычислительная техника и информатика. – 2014. – №1 – C. 130-139.
75. Chechina A., Churbanova N., Trapeznikova M. Different Approaches to the
Multilane Traffic Simulation // Traffic and Granular Flow ‘13. – 2015. – No. 19
– P. 361-368.
76. Трапезникова М.А., Чечина А.А., Чурбанова Н.Г. Описание динамики
транспортныхпотоковнаэлементахулично-дорожнойсетис
использованием двумерных математических моделей // Препринты ИПМ
им. М.В. Келдыша. – 2016. – №93 – 20 с.
77. Чечина А.А., Герман М.С., Ермаков А.В., Трапезникова М.А., Чурбанова
Н.Г.Моделированиеивизуализацияпотоковавтотранспортана
элементах улично-дорожной сети с использованием комплекса программ
САМ-2D // Препринты ИПМ им. М.В.Келдыша. – 2016. – № 124 – 17 с.
78. Трапезникова М.А., Чечина А.А., Чурбанова Н.Г. Двумерная модель
клеточных автоматов для описания динамики транспортных потоков на
элементах улично-дорожной сети // Математическое моделирование. –
2017. – т. 19 № 9 – С. 110-120.
79. ЧечинаА.А.Новыеалгоритмыперестроенияавтомобилейдля
микроскопической модели транспортных потоков на основе теории
клеточных автоматов // Препринты ИПМ им. М.В. Келдыша. – 2017. –
№136 – 14с.
80. Chechina A., Churbanova N., Trapeznikova M. Multilane Traffic Flow
Modeling Using Cellular Automata Theory // EPJ Web of Conferences. – 2018.
– V. 173 – P. 06003.
81. Chechina A., Churbanova N., Trapeznikova M., Ermakov A., German M.
Traffic flow modelling on road networks using cellular automata theory // Int.
J. of Eng. & Techn. – 2018. – V. 7 N. 2.28 – P. 225-227.
82. Чечина А.А. Воспроизведение экспериментальных пространственно-
временныхструктурвтранспортныхпотокахприпомощи
математической модели на основе теории клеточных автоматов //
Препринты ИПМ им. М.В.Келдыша. – 2018. – № 236 – 16 с.
83. Chechina A., Churbanova N., Trapeznikova M. Reproduction of experimental
spatio-temporal structures in traffic flows using mathematical model based on
cellular automata theory // Periodicals of Engineering and Natural Sciences. –
2019. – V.7 No. 1 – P. 76-81.
84. Churbanova N., Chechina A., Trapeznikova M., Sokolov P. Simulation of
traffic flows on road segments using cellular automata theory and
quasigasdynamic approach // Mathematica Montisnigri. – 2019. – V. XLVI – P.
72-90.
85. Churbanova N., Chechina A., Trapeznikova M., Ermakov A., German M.,
Sokolov P., Bozorov O. CMMSE-2019: Simulation and visualization of
vehicular traffic on road networks using high performance computing systems
// Computational and Mathematical Methods. – 2020 – V. 2 No. 3 – e1082.
86. Chechina A., Churbanova N., Trapeznikova M. Modelling traffic on road
junctions on parallel computing systems using cellular automata approach //
International Journal of Interactive Mobile Technologies (iJIM). – 2020. – V.
14 No. 10 – P. 178-185.
87. Trapeznikova M.A., Churbanova N. G., Chechina A.A., Ermakov A.V.,
German M.S. Supercomputer Technology for Traffic Simulation in a
Metropolis // 2020 International Conference on Engineering Management of
Communication and Technology (EMCTECH, Vienna, Austria, 2020) IEEE. –
2020. – P. 1-4.
88. Chechina A., Churbanova N., Trapeznikova M. Driver Behaviour Algorithms
for the Cellular Automata-Based Mathematical Model of Traffic Flows // EPJ
Web Conf. – 2021. – V. 248 – P. 02002.
89. ЧечинаА.А.,ЧурбановаН.Г.,ТрапезниковаМ.А..Сравнение
воспроизведения пространственно-временных структур транспортных
потоков при использовании различных способов осреднения данных //
Математическое моделирование. – 2021. – т. 33(1) – С. 25-31.
90. Чечина А.А. Алгоритмы поведения водителей на нерегулируемых
перекресткахсприоритетомиприобъездепрепятствий//
Математическое моделирование. – 2021. – т. 33(9) – С. 47-59.
Публикации автора в научных журналах
- 2 Инженерное дело, технологии и технические науки
- 09.00.00 Информатика и вычислительная техника
- 09.04.02 Информационные системы и технологии
- 02.00.00 Компьютерные и информационные науки
- 01.00.00 Математика и механика
- 1 Математические и естественные науки
- 02.04.02 Фундаментальная информатика и информационные технологии
- 2 Инженерное дело, технологии и технические науки
- 09.00.00 Информатика и вычислительная техника
- 09.04.02 Информационные системы и технологии
- 02.00.00 Компьютерные и информационные науки
- 01.00.00 Математика и механика
- 1 Математические и естественные науки
- 02.04.02 Фундаментальная информатика и информационные технологии
Помогаем с подготовкой сопроводительных документов
Хочешь уникальную работу?
Больше 3 000 экспертов уже готовы начать работу над твоим проектом!