CC BY
17
УДК 519.6, 656.11
К.К. Глухарев1'2, Н.М. Улюков1'2
1 Институт машиноведения им. А.А. Благонравова РАН 2 Московский физико-технический институт (Национальный исследовательский университет)
К теории автомобильных потоков
Рассматривается цепочка движущихся друг за другом автомобилей в однорядной полосе. Автомобили представляются точечными частицами, которые взаимодействуют с соседними в соответствии с законом безопасной дистанции. Формулировка закона безопасной дистанции основывается на известных данных наблюдения за автомобильными потоками [1]. Выводятся дифференциально-функциональные уравнения потока частиц. Строится решение уравнений потока частиц. При этом учитываются следующие обстоятельства. 1. Каждая внутренняя частица в цепочке одновременно является как лидером для следующей частицы, так и преследователем предыдущей частицы. 2. Гладкие траектории автомобилей приближённо могут представляться кусочно-линейными неубывающими функциями (запрещено движение «назад»). 3. Алгоритм смены скорости преследователя исходит из критерия поддержки безопасной дистанции. Приводятся результаты численного интегрирования уравнений, иллюстрирующие типовые потоковые явления. Формулируются направления сетевого развития теории дискретного потока.
Ключевые слова: транспортные потоки, дискретные потоки, закон безопасной дистанции, дифференциально-функциональные уравнения потока, типовые потоковые явления.
Дорожные сети мегаполисов все более насыщаются автомобилями. Следствием этого на магистралях наблюдаются протяженные очереди, а на перекрёстках спонтанно возникают скопления блокирующих друг друга автомобилей. Эти недостаточно изученные к настоящему времени явления все более тормозят потоки автомобилей и тем самым создают все более осязаемые проблемы — социальные, экологические и экономические.
Миграцию очередей в реальном времени на московской дорожной сети можно наблюдать на интернет-сайте «Карты дорожных пробок на Яндексе». Сайт экспонирует карту дорог г. Москва, на которую с интервалом 4-5 минут наносятся изменения в распределении дорожных «пробок».
Рассмотрим отрезок магистрали между перекрёстками. При красном сигнале светофора, установленном на выезде из магистрали в перекрёсток, поток автомобилей тормозится с образованием очереди из останавливающихся автомобилей. С включением зеленого сигнала головной автомобиль начинает движение, а следующий за ним становится головным. Если в некоторый момент времени тело очереди не пусто, то в этот момент к очереди может присоединиться следующий автомобиль, который становится хвостовым. Такую очередь из стоящих автомобилей назовем пробкой. Ее голова и хвост перемещаются по магистрали в направлении, противоположном движению.
Существуют два случая:
— пробка рассеивается за конечное время на достаточно длинном отрезке магистрали, когда интенсивность потока, входящего в тело пробки, меньше интенсивности отделения головных автомобилей;
— длина пробки не убывает с течением времени, когда интенсивность входящего потока больше или равна интенсивности выходящего потока из тела пробки.
Как в первом случае (при определённых отношениях длины отрезка магистрали и длительности красного сигнала), так и всегда во втором случае пробка за конечное время достигает соседнего перекрёстка и может тормозить на нём потоки с других магистралей.
Полное или частичное покрытие площади перекрёстка телами пробок, при котором отсутствует возможность отделения головных автомобилей от тел пробок, назовем затором.
Рассмотренный механизм образования заторов может приводить как к полному или частичному рассеиванию автомобилей на исходящих от перекрёстка магистралях, так и к росту длин пробок на входящих магистралях с захватом заторами соседних перекрёстков.
Отметим двойственную роль пробок на дорожной сети. С одной стороны, образование пробок является необходимым технологическим следствием светофорного регулирования. С другой стороны, первичная пробка индуцирует эффект затора, который в свою очередь приводит к неконтролируемому росту длин вторичных пробок, поражающих дорожную сеть.
Одной из целей построения модели и теории дискретного потока является установление условий, при которых первичная пробка успевает рассеяться в пределах магистрали, не достигнув соседнего перекрёстка. В этом контексте модель однорядного дискретного потока является первичным объектом для исследования.
В общем контексте модель дискретного потока частиц, которая сформулирована ранее в работах [2-4] и развитая в данной работе, ближе всего примыкает к модели В.А. Марченко [5], в которой цепочки частиц укладываются в систему каналов со структурой графа.
Принципиальное отличие моделей состоит в следующем. Модель В.А. Марченко предусматривает упругое взаимодействие соседних частиц. А модель данной работы предусматривает взаимодействие соседних частиц в соответствии с законом безопасной дистанции. Ещё раз подчеркнем, что этот закон основан на известных экспериментальных данных о зависимости дистанции между частицами от их скорости.
Другие представления о моделях потока автомобилей на магистральных отрезках дорог можно получить из работ [6-8]. Их более полный обзор содержится в [9].
Также заметим, что в исследованиях [1], по-видимому, впервые, проблема организации дорожного движения в крупных городах была сформулирована как системная, с попыткой создания сетевой концепции светофорного регулирования потоков автомобилей. Комбинаторный анализ различных схем расщепления потоков на регулируемом перекрёстке представлен в [10].
Данная работа, содержащая вывод ДФ-урав-нений и их интегрирование, является первой частью теории дискретного потока. Вторая часть этой теории — динамика очереди — будет представлена в следующих номерах журнала.
I. Модель дискретного потока
Схема дискретного потока. Рассмотрим однорядную полосу для движения автомобилей, с которой свяжем координатную ось S (рис. 1). Полоса наполнена движущимися друг за другом автомобилями, образующими то, что будем называть дискретным потоком. Относительно этого потока предположим следующее:
— автомобили представляются точечными частицами;
— запрещены въезды и выезды автомобилей через боковые ограничивающие линии;
— скорость каждого автомобиля в произвольный момент времени направлена по координатной оси или равна нулю;
— взаимодействие соседних частиц в потоке подчиняется закону безопасной дистанции.
ф ф © ф 1 I ф ф ф ф
S+ si+1 st s~ S
Рис. 1. Однорядный дискретный поток частиц. S+, Si, S- е S
Изложим содержание закона безопасной дистанции.
Предварительно заметим, что начало исследований механизмов взаимодействия автомобилей в потоке относится к 30-м годам прошлого века (см. [1]). В этот период было экспериментально установлено, что максимальная интенсивность потока автомобилей в однорядной полосе составляет величину 2000 авт/ч, которая достигается на скоростях 40-50 км/ч. Ясно, что величина 2000 авт/ч определяет пропускную способность однорядной полосы.
Закон безопасной дистанции. В рассматриваемом потоке частиц обозначим через si(t) координату г-й частицы (автомобиля) в момент времени t, si(t) G S, t G R, г G N. Дистанцию между соседними частицами с номерами г и г + 1 обозначим через Ai+i(t) = sj(t) — si+i(t), г G N и свяжем её со скоростями частиц. С этой целью рассмотрим стационарный однородный поток, в котором
— равны и постоянны скорости всех частиц: si(t) = s^+i (t) = v = const, t G R, г G N;
— равны и постоянны дистанции между соседними частицами: Ai(t) = Д = const, t G R, г ^ 2.
В плоскости R x S такому потоку соответствует семейство параллельных прямых с равными дистанциями (рис. 2). Очевидны формулы
Д
г' = д? q = v^
где fj, = q = здесь Д — дистанция между соседними частицами, Дt — интервал времени, в течение которого частица проходит дистанцию Д со скоростью v, ¡л — плотность частиц, q — интенсивность движения частиц.
t' teR
Рис. 2. Графики движения частиц, образующих однородный стационарный дискретный поток
Дистанция Д при прочих равных условиях является функцией скорости V. Эта зависимость в основном определяется двумя факторами — длиной тормозного пути и запаздыванием реакции водителя при изменении скорости впереди движущегося автомобиля.
На рис. 3 представлены графики описанных выше следующих функций: ц(у), ^¿(у) и
Ф).
Рис. 3. Графики функций: а — функция безопасной дистанции б — плотность потока ß(v), в — функ-
ция временного интервала At(v) иг — интенсивность потока q(v)
В каждой паре соседних частиц из однорядной полосы выделяем лидера и преследователя, относительно которых считаем, что
— преследователь в каждый момент времени имеет возможность сохранять свою скорость либо переключаться на скорость лидера, то есть: ii+iCO е Vi+1(i), Vi+1 (t) = {¿i+1(t_),Si(t)},i e N, где Si — скорость лидера с номером i, si+i — скорость преследователя;
— преследователь, имея выбор из двух скоростей, стремится поддерживать безопасную дистанцию, следуя правилам:
(а) если Ai+i(t) > p(Si(t)), то преследователь стремится сократить дистанцию Ai+i, выбирая максимальную скорость Si+i(t) = maxVi+i(t),
(б) если A^it < Si(t)), то преследователь стремится увеличить дистанцию Ai+i, выбирая минимальную скорость Si+i(t) = minVi+i(t),
(в) если выполняется равенство Ai+i (t) = ц>(Si(t)), то преследователь для поддержания безопасной дистанции должен двигаться со скоростью лидера, то есть Si(t) = Si+i(t);
— действия каждого лидера не зависят от действий преследователя;
— движение лидера с номером i = 1 считается заданным, например, указателями скорости, размещенными по длине полосы; в этом случае движению лидера с номером i = 1 может быть сопоставлено обыкновенное дифференциальное уравнение первого порядка c кусочно-постоянной правой частью. При этом скачки сопоставляются точкам полосы, где установлены указатели скорости.
Дифференциально-функциональные уравнения (ДФ-уравнения) однорядного
дискретного потока. Учитывая изложенное, уравнения дискретного потока с безопасной дистанцией представимы в следующем виде:
¿1 = f (si), si е S,
(1)
Sj+1 (t)
maxVi+i(t), Ai+i(t) > y(Si(t)),
Si(t), A+i(t) = p(Si (t)),
minVi+i(t), Ai+i(t) < p(Si(t)).
(2)
Здесь введены следующие обозначения: f(si) — кусочно-постоянная функция со скачками в точках полосы, где установлены указатели скорости;
Vi+i (t) = {Si+i(t-),Si(t)}; Ai+i(t) = si(t) - si+i(t), t е R, i е N; si(t)) — заданная неубывающая функция от скорости лидера — представление закона безопасной дистанции; si(t) — координата i-й частицы в однорядной полосе, si (t) — её скорость, si(t) е S, t е R, i е N.
Краевые условия. К уравнениям дискретного потока (1), (2) следует присоединить начальные и граничные условия. Запишем общее представление этих условий.
Начальные условия для уравнений (2) в начальный момент времени задаются распределением частиц в полосе, а также начальным распределением их скоростей; то есть в момент времени t = to задаются следующие два множества:
{si(to)}, si(to) > si+i (to) е [S+,S-), {si(to-)}, si(to-) > 0,
(3)
где г € 2,/о и 1о — число частиц в отрезке полосы [Б+,Б_) в момент времени ¿о, 1о € N, 1о < а Б +, 5_ — граничные точки этого отрезка.
Для интегрирования уравнения (1), как для обыкновенного дифференциального уравнения, в начальный момент времени £ = ¿о достаточно задать положение первой частицы в полосе
51 (¿о) € [б+,б_).
Граничное условие в сечении Б+ задается моментами времени проезда с ненулевой скоростью сечения Б + ¿-ми частицами, г = 1,1 +; то есть считаются заданными следующие множества:
{¿+ы <¿++1,
{¿¿г+_)}Мг+_) > 0.
Ясно, что при этом следует считать
(4)
5.;Д+) = Б+, г = 1,1+, /+<+ оо.
Здесь I + — общее число частиц, вошедших в отрезок полосы [Б+,Б_) за время [¿о, + го). Отметим, что множества из (4) считаются конечными.
Граничное условие в сечении Б_ задается аналогично граничному условию в сечении Б+ следующими множествами:
{t-},t- < ^v {si(t-)},si(t-) > 0
при условии
si(t- ) = S , i = 1,II <
Здесь I- — общее число частиц, вышедших из отрезка полосы [5 + ,5-) за время [¿о, + то). Также отмечаем, что множества (5) считаются конечными. Особенности других представлений краевых условий (5) излагаются далее.
II. Интегрирование уравнений потока частиц
Программное движение лидера. Пусть в соответствии с п. I движению лидера сопоставлено обыкновенное дифференциальное уравнение с кусочно-постоянной и неотрицательной правой частью следующего вида: в = I(в),в е [5+,5-) с Б,
I(в) = vn,s е [5(п-1),Б(п)) С [5 +,5-),
1,n+1,
где п — число скачков в интервале [5+,5-), п < (п) — координаты точек скачков ско-
рости, vn — значение скорости лидера (г = 1) на интервале [5(п-1) ,Б(п)). При этом 5(0) = 5+
и 5(п+1) = 5"Я+Я- <
s
S+ Sin-i) S(n) S- s(ES
Рис. 4. Скоростная программа движения лидера в однорядной полосе
Обозначим через s(t) решение уравнения на интервале [t',t''), t" < при начальных данных
si(t') = S + .
Считая функцию f (s) непрерывной справа, решение si(t) строим в классе кусочно-линейных неубывающих функций. Решение представимо следующим интегралом:
si(t) =
S + + vi(t — t'),t G [t',m),
S + + E m=i
vm ДПт + Vn(t)(t — Vn(t)), t G [r,i,t"),
nn = nn-i + Д'Пп, n+i
t" = nn = t' +Y, ДПт, (6)
m=i
п(£) = шах{п : г/п ^ где п — число скачков функции I(в) на интервале
[5+,5-), в е [5+,5-).
График функции I(в) представлен на рис. 4. Подчеркнем, что здесь все vn ^ 0.
Заметим, что если формально снять ограничение, связанное с неотрицательностью функции I(в), в е [5+ ,Б-), то возникают определённые проблемы с определением решения в точках типа в* (рис. 5), то есть в точках смены знака скорости частицы.
Рис. 5. Точка в* скачка скорости лидера со сменой знака, приводящая к «скользящему режиму»
Такая ситуация может возникать в случаях, когда частицам разрешен «задний ход».
Определение решения ДФ-уравнения.
Решение ДФ-уравнений (2) будем строить в классе кусочно-линейных неубывающих функций. Тогда под решением ДФ-уравнения (2) для каждого преследователя с номером г + 1 будем понимать функцию — интегральную ломаную, в£+1: И, ^ Б, г е N из указанного класса, которая для любого £ е [¿', ¿''), ¿'' < удовлетворяет уравнению (2), где [£', £'') — интервал определения решения для соответствующего лидера с номером г, а величина ' даётся представлением (6).
Напомним, что при г = 1 движение лидера задается обыкновенным дифференциальным уравнением с кусочно-постоянной правой частью, и поэтому может быть проинтегрировано независимо от всех других уравнений из (2) (г ^ 2).
Пусть для пары лидер-преследователь известен алгоритм построения решения для преследователя при известной интегральной ломаной лидера. Очевидно, что в этом случае решения всех уравнений из (2) могут быть построены индуктивно. Построим общий шаг индукции.
Вычисление точек скачков скорости преследователя. Обозначим через {1п}г+1 множество моментов времени на действительной оси К, где потенциально возможны скачки скорости преследователя с номером г + 1, {¿п}г+1 С [¿'¿''), г е N .С каждой точкой ¿п е {п}г+1 могут быть связаны следующие типы событий:
(а) событие типа «смена скорости лидера»;
(б) событие типа «достижение безопасной дистанции»;
(в) событие типа «нарушение безопасной дистанции».
В силу (1) и (2) легко устанавливаются следующие свойства:
n
если соседние точки скачков скорости преследователя обозначить через tn-i и tn (tn-i < tn), то на интервале [tn-i,tn) график интегральной ломаной преследователя есть отрезок прямой;
каждая точка скачка скорости преследователя имеет место при наступлении либо события типа (а), либо события типа (б), либо при одновременности событий типа (а) и (б);
событие типа (в) наступает только одновременно с событием типа (а).
Заметим, что при наступлении события типа «нарушение безопасной дистанции» скачок скорости преследователя отсутствует, так как безопасная дистанция при смене скорости лидера всегда терпит скачок, а траектория преследователя является непрерывной функцией.
Таким образом, алгоритм вычисления момента времени tn, соответствующего скачку скорости преследователя (при известном tn-i), основан на сопоставлении моментов времени события типа (а) и события типа (б).
Далее для вычисления точки tn e {tn}i+i при известной точке tn-i e {tn}i+i вводим в рассмотрение следующие множества:
Ti+i(tn-i) = [п : п> tn-i, п е {n}i} и [t"} = г, r;+i(tn-i) = [t : Ai+i(t) = v(Si(tn-i)}, t e [tn-i, minTi+i(tn-i))} и [t"} =
tn-i e [tn}i+i с [t',t"), i e N,
где [t"} — одноэлементное множество; [n}i — множество точек скачков скорости лидера с номером i,{n}i С [t',t''); а множество
[t : Ai+i(t) = y>(S i(tn-i)),...
t e [tn-i, minTi+i(tn-i))} выделяет интервал нулей уравнения
Ai+i(t) = Si(tn-i)),
t e [tn-i, min Ti+i(tn-i)).
Заметим, что введённые множества Tri и Tv не являются пустыми. Также заметим, что операция вычисления minTi}+i(tn-i) выделяет либо точку t'', либо следующую за точкой tn-i справа точку скачка скорости лидера с номером i. Операция minT^+i(tn-i) выделяет либо точку t'', либо точку, соответствующую событию типа «достижение безопасной дистанции», которая является ближайшей справа к точке tn-i или совпадает с tn-i.
Пусть вычислена точка tn-i e [tn}i+ii e N. Тогда, учитывая выделенные свойства для событий типа (а), (б) и (в), замечания о вычислении min Tn и min Tv, а также, учитывая, что на любом интервале [tn-i,tn) С [t',t'') интегральная ломаная есть отрезок прямой, для вычисления следующей точки tn потенциального скачка скорости
преследователя с номером i + 1 получаем рекуррентную формулу
tn =
min[min Ti+i(tn-i), min Tlv+i(tn-i)}, Ai+i(tn-i) = Si(tn-i)),
min Ti+i(tn-i)
П \"n-
Ai+i(tn-i) = ф(Si(tn-i)).
(7)
При n = 1 считается, что to = t'.
Алгоритм преследования лидера.
Считаем, что на отрезке [t',t''), t'' < +( построена интегральная ломаная для i-й частицы. Тогда алгоритм построения интегральной ломаной для преследователя с номером i + 1 на отрезке [t',t'') представим в следующем виде:
строится множество Vi+i(t') = [Si+i(t'_),Si(t')}; в силу (2) однозначно вычисляется скорость преследователя Si+i(t');
строится первый отрезок ломаной с конечной точкой Si+i(ti) = Si+i(t') + Si+i(t')(ti - t'), где
min [min Ti]+i(t'), min T^+i(t')}, t =, Ai+i(t') = ф&')), ti = S min Ti1+i(t'),
Ai+i(t') = <p(Si(t'));
далее считается, что построен отрезок ломаной с номером n - 1; тогда, учитывая, что множество Vi+i (tn-i) = [Si+i(t(n-i)-),Si(tn-i)}, а скорость Si+i(tn-i) вычисляется в силу (2), конечная точка n-го отрезка ломаной вычисляется по формуле Si+i(tn) = Si+i (tn-i ) + Si+i(tn-i )(tn - tn-i), а
tn определяется выражением (7). Процесс построения ломаной заканчивается, если при некотором n' < (
min [min Tn+i(tn'), min ^(tn)} = t''.
Ясно, что в силу непрерывности (Si(t'') = Si(t''_)) имеем
Si+i(t'') = Si+i (tn' ) + S i+i(tn' )(t'' - tn' ),
а точка с координатами (t'',Si+i(t'')) — крайняя справа точка построенной ломаной на плоскости R x S, может быть начальной для продолжения решения при t ^ t''.
В результате искомая ломаная определяется интегралом
Si+i(t') + Si+i(t')(t - t'),t e [t',ti),
n(t)
E r
Si+i(t) = { Si+i(t') Si+i(tm-i)(tm - tm-i)+
i
+ Si+i(tn(t))(t - tn(t)),t e [ti,t''),
е п,т = 1 ,М1+\ АР+1 = |{*„}*+1| < 00,
где п(Ь) = тах{п : ¿п ^ ¿}, ¿п € {¿п}г+1, а моменты времени ¿п потенциальных скачков скорости преследователя определяются формулой (7). То, что построенная функция Si+l(t) для каждого £ € [¿',Ь") удовлетворяет системе (2), проверяется её подстановкой в систему (2).
Некоторые свойства частиц в потоке.
1. Пусть траектории г-го лидера сопоставлена функция si(t) — кусочно-линейная неубывающая функция с конечным числом отрезков прямых, удовлетворяющая (2) для всех £ е [¿',£''), ¿'' < то. Тогда движение преследователя в.1+1(1) с номером г + 1 с начальным положением (£') < в^') и скоростью si+l(t'_) ^ 0 однозначно продолжаемо в силу (2) на отрезок £ е \Ъ'¿''),£'' < то в классе неубывающих ломаных с конечным числом отрезков прямых.
2. Переключение скорости преследователя может происходить в моменты времени переключения скорости лидера, а также в моменты времени достижения преследователем безопасной дистанции.
3. На каждом интервале постоянства скорости лидера может быть не более одной точки смены скорости преследователя.
4. Пусть при £ =
Дi+l(t') = в^') — Si+l(t') > 0,
тогда для £ е [¿'¿'')
Дi+l (*) > 0.
5. Пусть движение частиц в однорядной полосе подчиняется уравнениям (1), (2). Тогда для конечного отрезка однорядной полосы в цепочке двигающихся друг за другом частиц для каждой их пары сохраняется порядок их следования. Таким образом, множество частиц в отрезке однорядной полосы есть параметризованное множество с отношением строгого порядка (транзитивность и антисимметричность), где параметр — текущее время г е [£',£''),£'' < то.
Примеры численного моделирования. На рис. 6 (а) и (б) приведены примеры численного моделирования дискретного потока, представленного ДФ-уравнениями (1) и (2).
О 10 ь
Рис. 6. Результаты численного моделирования очереди частиц: а — случай роста и б — случай стабилизации длины очереди
В плоскости К х 5 на рис. 6 рассмотрим последовательность £-сечений. Такое рассмотрение легко обнаруживает ранее рассмотренные свойства пробки.
Поток с достаточно большим скоплением частиц становится прерывистым (стоп-старт движение). В таком потоке периодически образуется очередь частиц. На рис. 7 представлен поток в режиме стоп — старт.
III. Замечания о других моделях дискретного потока
О дискретном потоке с импульсным управлением. Исходим из теоремы об изменении количества движения:
mX(t) = mX(t — At) +
г-Аг
где т — масса частицы, х — её скорость, Г — ак-
г
тивная сила, действующая на частицу, J —
г-Аг
импульс силы на интервале [£ — Д£,£], Д£ > 0. Пусть равномерно по £ существует предел:
Рис. 7. Поток в режиме стоп-старт
lim (X(t — At) +
Fd£) =
г- Дг
t
max V (t),
= 4 s(t),
min V (t),
Д^) >^(s *(t)), Д(t) = Ф * (t)), Д(t) <^(s *(t)),
у (£) = {¿(^т,
Д(£) = в*(Ь) — в(£),
где в* — скорость лидера, в — скорость преследователя, в* и в — положения частиц в полосе, а значение предела является правой частью уравнения (2).
В таком представлении наша модель дискретного потока частиц является системой конечного числа материальных точек в канале без трения с импульсным управлением, которые поддерживают безопасную дистанцию между лидером и преследователем.
В случаях ограниченности активных сил, Г е [Г- ,Г], где Г- и Г конечны и характеризуются ресурсом двигателя автомобиля и сцеплением колес и дороги, траектории частиц как материальных точек должны иметь непрерывные первые и ограниченные вторые производные.
При таком представлении вопрос построения уравнений для «гладких» потоков может быть сведен к задаче о сплайновом сглаживании ранее построенных ломаных ([11], с. 168). Особенностью этой задачи (как математической) является следующее. Заданное поле ломанных, образующих дискретный поток, приближается (в некоторой метрике) на классе сплайнов, в котором:
— сглаживающие сплайны должны иметь неотрицательную и непрерывную первую производную;
— вторые производные — ускорения материальных точек в гладком потоке, должны удовлетворять неравенствам
Г < тв(£) < Р.
Заметим, что вычислительная сложность задачи связана как с выбором сетки узлов, так и с возможностью построения решения на сплайнах второго порядка.
Далее построенные сглаживающие сплайны рассматриваем в качестве заданных движений материальных точек в потоке. В результате вопрос построения соответствующих управляющих сил может быть сведен к решению так называемой первой (основной) задаче динамики — по заданному движению определить силы его вызывающие ([12], с. 321). Обзор таких задач, трактуемых как построение программных движений, содержится
в [14].
Другая трактовка вопроса о преследовании лидера на безопасной дистанции может быть связана с задачей о наложении сервосвязи, введённой А. Бегеным ([13], с. 344). Реализация такой связи (иначе — обратной связи) предусматривает наличие измерителей, силовых приводов и автоматических регуляторов. Обзор алгоритмов стабилиза-
ции обратных связей, в частности, использующие скользящие режимы, содержится в работе [15].
Таким образом, с одной стороны, обсуждаемая модель потока на классе кусочно-линейных траекторий является приближением (в указанных смыслах) к потокам механических частиц. С другой стороны, использование на автомобилях форсированных двигателей допускает режимы с быстрым изменением скорости. Кроме того, построение программируемых потоков даёт перспективу перехода на автоматическое вождение автомобилей.
О связи дискретного потока с моделью конечного автомата. Вводим следующие равномерные сетки:
— на координатной оси S строим сетку с шагом 5, 5 = const > 0;
— на оси времени R строим сетку с шагом Д^
— вводим конечное множество {vm}, vm = mAv, Av = const > 0, m = 0,M, M < oo.
Устанавливаем между шагами сеток следующую связь: ДtДv = 5. Кроме того, функцию заменяем на её дискретное приближение
В этом случае предположим, что начальные положения каждой i-й частицы si (to), i G N, соот-
Sj(t о)
о)
ветствуют узлам ¿-сетки, то есть г е N. Кроме того, будем считать, что скачки скорости лидера цепочки (г = 1) допускаются только в момент времени Ьк = кД£, к е N, = 0. Также предположим, что скорости лидера цепочки и всех других частиц (г ^ 2) в сечении входа в отрезок полосы считаются кратными шагу Дv; тогда при всех других £к, к е N, все частицы в силу уравнений (1) и (2) с дискретным временем £к будут располагаться в узлах ¿-сетки на координатной оси 5.
Пусть М — число узлов в ¿-сетке отрезка полосы, а т — число двоичных разрядов, требуемых для записи числа М. Пусть также имеет место равенство log2 М = т. Тогда число возможных представлений цепочки частиц в конечной полосе с числом узлом М равно числу 22 . Это утверждение следует из факта: множество булевских функций с т переменными равномощно множеству различных цепочек частиц на сетке с М узлами, если только т и М связаны равенством т = к^2 М.
IV. О развитии теории дискретного потока
Алгоритм построения сети. Эскизно наметим алгоритм построения дорожной сети из типовых элементов с их топографическим представлением на декартовой плоскости.
Первый тип элементов — отрезки ориентированных однорядных полос (отрезки полос или кратко ОП). ОП допускают искривления в плоскости и малые их деформации при выводе их из
плоскости с целью скрещивания двух и большего числа ОП. Связь отдельных ОП осуществляется через узлы коммутации (кратко — коммутаторы, например, в виде зоны регулируемого перекрёстка двух дорог). В коммутаторах входы одних ОП соединяются с выходами других ОП. Считается, что каждый ОП связан с двумя коммутаторами. Коммутаторы так же, как и ОП, имеют топографическую привязку.
Построенную односвязную область, состоящую из коммутаторов и магистральных отрезков (параллельных и одинаково ориентированных ОП), будем называть магистральной сетью. Эта сеть отделяет друг от друга области, называемые Р-зонами. Считается, что в Р-зонах допускаются остановки, стоянки и парковки транспортных средств. Также считается, что в Р-зонах может существовать собственное внутреннее дорожное устройство. Р-зоны формально следует считать открытыми множествами; тогда область, занятая магистральной сетью, будет замкнутой. Заметим, что к Р-зонам относится внешняя относительно сети часть плоскости. В магистральную сеть дополнительно врезаются коммутаторы, связывающие её с Р-зонами, и, таким образом, решается вопрос связи между собой всех Р-зон.
Построение сети может быть завершено прокладкой допустимых маршрутов между Р-зонами. Маршруты на карте сети выделяются маршрутными нитками, которые можно рассматривать как допустимые траектории движения транспортных единиц.
Отметим принципиальную особенность коммутатора. Его коммутационная зона покрыта «паутиной» ниток — траекторий, соединяющих входящие и выходящие ОП, примыкающие к коммутатору. Поэтому тактовое действие светофорного регулятора может быть представлено следующим образом: попеременное открытие и закрытие движения по отдельным группам траекторий, при котором в каждый момент времени исключаются точки пересечения ниток.
Маршрутизация как алгоритмическая проблема. Построенную сеть со светофорным регулированием наполним подвижной дискретной средой — потоками частиц, которые локально подчиняются закону безопасной дистанции. В результате получаем новый объект. Этот объект, с одной стороны, является топографическим представлением дорожной сети, которое напоминает ориентированный граф. С другой стороны, дискретная среда частиц, помещённая в сеть, создаёт совокупный поток, действие которого должно подчиняться следующим целевым установкам:
— каждая частица должна реализовывать свою, установленную ею самой, последовательность прохождения Р-зон;
— управление совокупным потоком этих частиц должно устойчиво блокировать образования заторов.
Дискретная среда может содержать порядка 107 частиц, поэтому вопрос выбора маршрутов для каждой частицы является весьма трудоемкой задачей. В этом контексте проблему маршрутизации следует рассматривать как алгоритмическую, например, как задачу о существовании некоторого единого правила прокладки маршрутных ниток между любыми парами Р-зон.
Тезисно выделим ряд базовых вопросов в теоретических исследованиях нового объекта.
— Аксиоматика сети, связанной с её топографическим и теоретико-графовым представлениями, а также алгоритмов сетевой сборки.
— Способ варьирования сети; в том числе деформации сети, например, преобразование поперечно-продольной схемы в радиально-кольцевую; преобразование сети при включениях и удалении сетевых элементов; топологические характеристики и критерии сравнения сетей с оценкой сложности их маршрутно-ниточного покрытия.
— Теория коммутатора, включающая модели прохождения и преобразования частиц в коммутаторах потоков; модель светофорного регулирования в коммутаторе; механизмы формирования заторов в площади коммутатора.
— Краевые условия на границах ОП и коммутаторов; общие уравнения дискретных потоков, условия существования сетевых напряжённых потоков.
— Полнота данных в задаче управления потоками частиц; принципы и оценки возможных ресурсов управления.
В заключение подчеркнем, что новый объект — сеть с вложенными в неё дискретными потоками частиц — может быть полезен при разработке других общесетевых тематик, в том числе при проектировании вычислительных наноструктур.
Литература
1. Иносе Х, Хамада Т. Управление дорожным движением. — М.: Транспорт, 1983. — 248 с.
2. Глухарев К.К. К теории управления потоками автомобилей на дорожной сети // Тезисы докладов VII Международного семинара «Устойчивость и колебания нелинейных систем управления». — М., 2002. — С. 53.
3. Улюков Н.М. Численное моделирование потока автомобилей в однорядной полосе // Труды Х1^Ш научной конференции МФТИ, Ч. III. — М.: МФТИ, 2005. — С. 56--58.
4. Глухарев К.К., Улюков Н.М. Светофорное регулирование однородного дискретного потока // Труды научной конференции МФТИ, факультет Аэрофизики и космических исследований. — М.: МФТИ, 2006. — С. 98--100.
5. Марченко В.А. Прямые и обратные задачи многоканального рассеивания: доклад на семинаре МИАНа. — 26 мая 2006. — www.mi.ras.ru.
6. Lee H.K., Lee H.-W., Kim D. Macroscopic traffic models from microscopic car-following models // Phys. Rev. E64. 2001.
7. Lubashevsy I.A., Mahnke R. Order parameter model to unstable multilane traffic flow. — arxiv. org, November 13, 2000. — V. 2. — http://arxiv.org/abs/cond-mat/9910268.
8. Ковшов Н.В. Математические модели многополосного дорожного движения / / Труды XLVII научной конференции МФТИ, Ч. III. — М.: МФТИ, 2004. — С. 197.
9. Семенов В.В. Математическое моделирование динамики транспортных потоков мегаполиса. http:// www.keldysh.ru/papers/2004/prep34/prep2004_34. html.
10. Яблонский С.В. [и др.\. Устройство для исследования алгоритмов регулирования уличного движения // ДАН CCCP. — 1960. — Т. 132, № 1.
11. Стечкин С.Б., Субботин Ю.Н. Сплайны в вычислительной математике. — М.: Наука, 1976, 248 с.
12. Кильчевский Н.А. Курс теоретической механики. Т. I. — М.: Наука, 1972. — 456 с.
13. Аппель П. Теоретическая механика. Т. II. — М.: Изд. Физико-математической литературы, 1960. — 487 с.
14. Галиуллин А.С. [и др.\. Построение систем программного движения. — М.: Наука, 1971. — 352 с.
15. Теория систем с переменной структурой / под ред. С.В. Емельянова. — М.: Наука, 1970. — 592 с.
16. Глухарев К.К., Улюков Н.М. Об одной модели однорядного потока автомобилей. Вывод уравнений и их интегрирование // Проблемы машиноведения и надежности машин (машиноведение). — М., 2008. — № 4.
17. Глухарев К.К., Улюков Н.М. Об одной модели однорядного потока автомобилей. Динамика очереди // Проблемы машиноведения и надежности машин (машиноведение). — М., 2009. — № 1.
18. Глухарев К.К., Улюков Н.М. Введение в теорию автомобильных потоков. Динамика однорядного движения: учеб.-метод. пособие. — М.: МФТИ, 2008. — 44 с.
Поступила в редакцию 24.03.2009.
