Next Article in Journal
MIP-Based Sensors: Promising New Tools for Cancer Biomarker Determination
Next Article in Special Issue
Learning Traffic as Images: A Deep Convolutional Neural Network for Large-Scale Transportation Network Speed Prediction
Previous Article in Journal
Imprinting of Microorganisms for Biosensor Applications
Previous Article in Special Issue
Health Monitoring and Evaluation of Long-Span Bridges Based on Sensing and Data Analysis: A Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Density Estimation of an Urban Freeway Network Based on Dynamic Graph Hybrid Automata

College of Metropolitan Transportation, Beijing University of Technology, Beijing 100124, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(4), 716; https://doi.org/10.3390/s17040716
Submission received: 24 January 2017 / Revised: 20 March 2017 / Accepted: 23 March 2017 / Published: 29 March 2017
(This article belongs to the Special Issue Sensors for Transportation)

Abstract

:
In this paper, in order to describe complex network systems, we firstly propose a general modeling framework by combining a dynamic graph with hybrid automata and thus name it Dynamic Graph Hybrid Automata (DGHA). Then we apply this framework to model traffic flow over an urban freeway network by embedding the Cell Transmission Model (CTM) into the DGHA. With a modeling procedure, we adopt a dual digraph of road network structure to describe the road topology, use linear hybrid automata to describe multi-modes of dynamic densities in road segments and transform the nonlinear expressions of the transmitted traffic flow between two road segments into piecewise linear functions in terms of multi-mode switchings. This modeling procedure is modularized and rule-based, and thus is easily-extensible with the help of a combination algorithm for the dynamics of traffic flow. It can describe the dynamics of traffic flow over an urban freeway network with arbitrary topology structures and sizes. Next we analyze mode types and number in the model of the whole freeway network, and deduce a Piecewise Affine Linear System (PWALS) model. Furthermore, based on the PWALS model, a multi-mode switched state observer is designed to estimate the traffic densities of the freeway network, where a set of observer gain matrices are computed by using the Lyapunov function approach. As an example, we utilize the PWALS model and the corresponding switched state observer to traffic flow over Beijing third ring road. In order to clearly interpret the principle of the proposed method and avoid computational complexity, we adopt a simplified version of Beijing third ring road. Practical application for a large-scale road network will be implemented by decentralized modeling approach and distributed observer designing in the future research.

1. Introduction

A traffic flow network is a typical example of a complex network system [1,2]. Dynamic modeling of the traffic flow network is quite a challenging task in traffic flow analysis and control. Various traffic flow models have been proposed to study complex traffic phenomena on a road network, e.g., as the mathematically most consistent macroscopic traffic flow, the well-known LWR model proposed by Lighthill, Whitham [3] and Richards [4], which formulates the relationship among the key traffic flow parameters such as density, flow etc. using continuous time partial differential equations. The second-order model is presented by Payne [5]. From the view of application, Daganzo [6,7] first proposed the Cell Transmission Model (CTM) by spatially and temporally discretizing the LWR model. Later, many researchers further improved the CTM. For example, Lebcque [8] considered the merge law and the diverge law of traffic flow according to the proportional distributions rule. Fl o ¨ tter o ¨ d and Nagel [9] developed the model by considering the merge and diverge law of multiple cells. Gomes and Horwitz [10,11,12] proposed the Asymmetric Cell Transmission Model (ACTM) and applied it to the freeway on-ramp metering control. Lo and Lin [13,14] amended the CTM to describe road networks with intersections. Sumalee et al. [15] developed the Stochastic Cell Transmission Model (SCTM). Canudas-de-Wit [16] proposed a graph constrained CTM.
These modified CTMs mainly focus on rules of transmitted traffic flow but give too little care to the description of road network structure, which leads to a lack of modularization and expandability in modeling a large-scale traffic network. Moreover, the nonlinear relationship between traffic flow and density results in study difficulty. Although Mu n ˜ oz and Sun et al. [17,18] presented the Switching Mode Model (SMM) which is a piecewise linearized version of the CTM, this model may be improved further by introducing some modeling tools as solid theoretical basis. Hence, in the previous work [19,20,21] we proposed a model framework by introducing dynamic digraph and hybrid automata. Different from the traditional digraph method describing the road topology [22], we adopted dual digraph idea and introduced dynamic digraph [23] into the description of traffic flow over a road network, where road segments are marked vertex, and traffic flow directions are considered as the directed edges of digraph. In our opinion, it has at least two advantages. Firstly, dynamic weights of describing traffic state in a road segment can be added in vertices as customary as in the weight digraph theory. Secondly, traffic flow directions can be clearly indicated by the directed edges and transmitted volumes can be expediently described by dynamic weights and switchings on edges, where the edge switchings between vertices can describe the traffic light signals.
On the other hand, we point out that hybrid automata, as a kind of hybrid dynamic system models [24,25], are an appropriate model to describe a multi-mode dynamic process. For instance, references [26,27] applied the hybrid automaton theory to automated highway systems (AHS). Lei and Ozguner [28] presented a hybrid automaton model for a single intersection for the first time. Chen et al. [29,30,31] studied the optimal control of single intersection and the cooperative control between intersections by using hybrid automata. Mu n ˜ oz and Sun et al. [17,18] presented the Switching Mode Model (SMM) for the California Freeway. In the model framework proposed in [19,20,21] we introduced hybrid automata to describe traffic densities in road segments.
Consequently, by combining dynamic digraph with hybrid automata and applying the parallel composition law between hybrid automata [32,33], we can obtain a new networked model and name it Dynamic Graph Hybrid Automata (DGHA). As an independent model framework, the DGHA can be applied in various dynamic network systems with multi-mode switchings. When applying it to model traffic flow of a road network, one can adopt different types of traffic variables as the state of the hybrid automata, such as two-dimensional variables of traffic density and vehicle average velocity. One may also apply, in the place of the hybrid automata, more complex uncertain models, and more modes taken when piecewise linearizing the nonlinear fundamental diagram. However, in this paper, as well as in [19,20,21], we embedded CTM into DGHA and thus only the traffic density variable was used as the state.
In this paper, we first further improve the DGHA framework proposed in [19,20,21] and the modeling procedure of urban freeway networks. More in detail, we first restate the definition of DGHA model more clearly. Then we analyze types and number of combination modes of cells (the partitioned road segments) when embedding CTM into the DGHA. Compared with the existing literature [34], we not only consider the multi-mode switchings of vertex state, but also take into account the multi-mode switchings of edge weights. In other words, the nonlinear expression of the density in a cell is described by a linear hybrid automaton and the nonlinear expression of the transmitted traffic flow between two cells is transformed into piecewise linear functions in terms of multi-mode switchings. We distinguish two modes (free flow and congested flow) in each cell and two modes (upward traffic wave and downward traffic wave) for each of the three types of connections between cells (linear connection, merge connection and diverge connection). Although the number of the combination modes is larger than that in the existing literature [34], e.g., about 5400 possible combination modes for a road partitioned into 10 cells (see the proof in Section 3.5), the combination modes not only describe the traffic state in a road segment but also the relative flow feature between two road segments. Next we give a convex polyhedral description of the combination modes by introducing a partition of state subspace. Finally we deduce an expression of Piecewise Affine Linear System (PWALS) for the whole freeway network. We point out emphatically that the modeling procedure is modularized and rules-based. Thus we develop an algorithm to implement modeling procedure of traffic networks by accomplishing parallel composition of cells with the help of a computer program. So compared with the traditional manual calculation [17,18,34], the modeling procedure is easily-extensible with the help of the algorithm and can be applied to traffic networks with arbitrary topology structures and sizes, including the road networks with signal intersections.
The proposed PWALS model can be utilized in many aspects, such as traffic state estimation and traffic flow control. In this paper, we focus on the estimation problem of traffic densities. Over the past decades, abundant work on this problem has been done by using different approaches, including the Kalman filter [22,35,36,37], particle filtering [38], the set-valued estimation [39,40,41]. These approaches are not based on networked hybrid automata and thus cannot be directly applied to our model.
On the other hand, the simpler state observer has been well-known as a powerful estimation tool. However, for the DGHA, we need to consider a switching observer. In the existing literature, the design approaches of switching observer can be divided into two classes, depending on whether the active modes are known or not. When the active modes are known, one can utilize the switched Luenberger-type observer with the same switching rules as in the model. The corresponding observer gains can be found by using the Lyapunov function approach for the switched dynamics of estimation error system and solving Linear Matrix Inequalities (LMI) [42,43,44,45,46]. Literature [47] proposed a switching observer for the systems with continuous state jumps and based on a known and fixed switching signal sequence. In the case when the active modes are unknown, meanwhile one is required to identify the modes on the basis of the observation of the output of the system over a certain interval. Pettersson [48,49] proposed multiple quadratic Lyapunov functions to design observers for a switched linear system with continuous state jumps. Literature [50] also deals with the observer design problem for continuous-time switched linear systems with unknown switchings.
The authors in [16] proposed a graph constrained CTM observer. In [51] they further introduced a robust mode selector for the uncertain graph-constrained switching mode model and applied it to highway traffic density estimation via a switched state observer. Alvarez-Icaza et al. [52] designed adaptive observer to estimate vehicle density.
As to the sufficient condition of observer design, the system must be observable, this problem is closely related to the traffic sensor placement in the traffic network. References [53,54] presented an algebraic approach to understand the problem of identifying which subsets of OD-pair and link flows can be calculated based on a given subset of observed OD-pair and link flows, and further applied the algorithm on the Nguyen-Dupuis network problem. A new framework which investigates observability in terms of flow and routing information on network arcs was proposed in [55].
In this paper, we focus on the design of a switched state observer based on the deduced PWALS model to estimate traffic densities under the assumption that the active modes are arbitrary but known. We will discuss more practical issues when the active modes are unknown in the future research.
The main contributions of the paper can be summarized in the following two aspects: (1) The DGHA modeling framework is proposed and a PWALS model is deduced when embedding the well-known CTM into the DGHA. Particularly, this modeling procedure is modularized and easily-extensible and can be used to model traffic networks with arbitrary topology structures and sizes; (2) A switched state observer is designed to estimate traffic densities of an urban freeway network. The corresponding observer gains are found by using the Lyapunov function approach for the switched dynamics of estimation error system and solving Linear Matrix Inequalities (LMI).
The rest of the paper is organized as follows. In Section 2, we provide some relevant background material about dynamic graph and hybrid automata and propose the modeling framework of Dynamic Graph Hybrid Automata. In Section 3, the DGHA framework is used to model traffic flow network by embedding the CTM into it. In Section 4, we deal with the design of the switched state observer. In Section 5, we illustrate the obtained results by applying them to the simplified version of Beijing third ring road. The paper is concluded in Section 6.

2. Preliminaries

First of all, we briefly review dynamic graph theory [23] and hybrid automata theory [25] and then combine dynamic digraph with hybrid automata to propose a new model framework which we call Dynamic Graph Hybrid Automata (in short, DGHA) [19,20].

2.1. Dynamic Graph

A graph [23] is a pair ( V , E ) , where V is a finite set of vertices or nodes, and E is a set of undirected edges, each being an unordered pair i , j and expressing a link between two vertices i , j V . A digraph is a pair ( V , E ) , where V is a finite set of vertices or nodes, and E is a set of directed edges e i j V × V , i j , expressing a link from vertex i V to vertex j V .
A vertex-weighted graph or digraph is a triple V , E , F , where ( V , E ) is a graph or digraph, F : V X is a function to assign a weight x i X for every vertex i, and X is a pre-specified set. An edge-weighted graph or digraph is a triple V , E , G , where ( V , E ) is a graph or digraph, G : E Y is a function to assign a weight y i j for every edge i , j or e i j , and Y is a pre-specified set. A fully-weighted graph or digraph V , E , F , G has weights assigned to both vertices and edges. Depending on the context, the weight sets X , Y may be real numbers, complex numbers, integers, or even elements of a group or a field, etc.
The set V of vertices of a graph or digraph ( V , E ) is said to be dynamic if its number | V | is time-varying. The set E of edges of a graph or digraph ( V , E ) is said to be dynamic if the edges are time-varying on-off switches. The vertex-weight function (respectively, edge-weight function) is called dynamic if the weights x i (respectively, y i j ) are time-varying. A fully-weighted graph or a digraph is called to be dynamic when any one of the four entities V , E , F , G is dynamic.
Remark 1.
Dynamic set V of vertices means that some vertices may be added or removed. When a vertex is added or removed, the corresponding edge linking is also appended or eliminated. Dynamic set E of edges implies that edges may be added or deleted over time. We also say that the graph or digraph has a switched set of edges.
Remark 2.
In what follows, for a digraph ( V , E ) , we denote by Pre ( i ) = { l V : e l i E } the set of all upstream neighbouring vertices and Post i = { j V : e i j E } the set of all downstream neighbour ones of vertex i, respectively. Thus N i = Pre ( i ) Post ( i ) is the set of all neighbour vertices of vertex i. We also introduce notation N ¯ i = i N i .

2.2. Hybrid Automaton

A hybrid automaton [25] with continuous input and output H is a collection H = S , Init , U , f , δ , Inv , G , R , Y , g , where
  • S = Q × X is a state space of hybrid state variables q , x , where Q is a finite set of discrete states q and X R n is a n-dimensional state space of continuous state vector x;
  • Init S is a set of initial hybrid states;
  • U R m is a set of continuous input variable u;
  • f : S X is a vector field describing the continuous state dynamics defined by
    x t + 1 = f ( q t , x t , u t ) ;
  • δ : S Q is a transition map describing the discrete state dynamics. One can also use a binary relation set Δ = Q × δ ( S ) Q × Q to express the discrete state transition;
  • Inv : Q 2 X defines the domain of continuous state vector under each discrete state. The domain is also called invariant set of continuous state vector;
  • G : Δ 2 X defines a guard condition for each discrete state transition;
  • R : Δ X is a reset map to assign a new initial state to the continuous state variable after the transition of the discrete state.
  • Y R p is a set of continuous output variable y;
  • g : S Y is an output map defined as
    y t = g ( q t , x t ) ;
A hybrid automaton can also be expressed by a digraph as shown in Figure 1 [25]. We notice that in this paper we focus on discrete time system defined by a difference equation for the continuous part of the hybrid automaton.

2.3. Dynamic Graph Hybrid Automata

Now we combine a dynamic graph with hybrid automata to propose a new model which we call Dynamic Graph Hybrid Automata (DGHA, for short). An early version of the model was given in references [19,20,21]. In what follows, we restate it in the way of parallel composition of hybrid automata.
Definition 1.
DGHA model is a hybrid automaton H over a digraph ( V , E ) which described by the following four parts:
(1)
Fully weighted dynamic digraph. In G = V , E , H , Φ , V = 1 , , N is a finite set of N fixed vertices, and E = e i j : i , j V is a set of directed edges, H and Φ are sets of the automata and the functions that assign the vertex weights and edge weights, respectively, which are described in detail as follows;
(2)
Vertex dynamics. The elements in H = H i : i V are hybrid automata H i = S i , Ini t i , U i , f i , δ i , Inv i , G i , R i , Y i , g i to describe the dynamics of vertex i V , whose components are similarly defined as in Subsection B.
(3)
Edge weight functions. The set Φ = ϕ i j : e i j E consists of the functions ϕ i j : Y i × Y j R m that assign a weight ϕ i j ( y i , y j ) for each edge e i j E ;
(4)
Composition of hybrid automata. In DGHA = V , E , H , H = S , Init , f , δ , Inv , G , R is the closed hybrid automaton over the digraph ( V , E ) obtained by parallel composition of hybrid automata in H , i.e., S = Q × X with Q = Q 1 × × Q N and X = X 1 × × X N ; Init = Ini t 1 × × Init N ; f : S X is defined as f = f ¯ 1 × × f ¯ N and f ¯ i : S i X i describes the dynamics of vertex i by substituting u i ( t ) in the the following equation
x i t + 1 = f i ( q i t , x i t , u i t ) ,
with the following relation:
u i ( t ) = l Pre ( i ) ϕ l i ( y l ( t ) , y i ( t ) ) j Post ( i ) ϕ i j ( y i ( t ) , y j ( t ) ) ,
where S i denotes the projected subspace of S on the state components with indices in N ¯ i , and for any i V the output y i is given by
y i t = g i ( q i t , x i t ) ;
transition map δ : S Q is defined by the following rule:
δ s 1 , , s N = q 1 , , q N δ i s i = q i , i V ;
Inv = In v 1 × × In v N ; G = G 1 × × G N ; R = R 1 × × R N .
A DGHA model can also be expressed by a digraph as shown in Figure 2.
Remark 3.
In Definition 1, the vertex set V and edge set E are assumed to be fixed. However, instead of the fixed edge set E, one can use a dynamically switching edge set E t = δ i j ( t ) e i j : i , j V by introducing a Boolean function δ i j : R { 0 , 1 } to indicate appearing or disappearing of edges, which can be used to express the traffic lights, i.e., δ i j ( t ) = 0 if the light is red, and δ i j ( t ) = 1 if the light is green. In this case, the edge weight function will be also affected by the switching function and thus become δ i j t ϕ i j y i t , y j t . In this paper, we focus on modeling urban freeway, so we do not further discuss this kind of models more in detail.
Remark 4.
In general, the functions ϕ i j y i , y j are nonlinear. When the hybrid automata of vertices are all linear or affine linear, one can further express these functions by linear ones in a way of multiple mode switchings. In this case, the vertex hybrid automata may be changed accordingly, depending on the edge modes and the modes of all neighbour vertices (see the application below in traffic flow model).

3. DGHA Model For Urban Freeway Network

As a generalized model framework, the DGHA can be applied to various complex network systems in different fields. When applying it to a traffic flow network, one can adopt different types of traffic variables to describe traffic state in road segments, e.g., two-dimensional traffic variables of traffic densities and vehicle average velocity. But in this paper we apply the DGHA to model traffic flow of an urban freeway network by embedding the well-known CTM into the DGHA. That means we only use the traffic densities to describe traffic state in road segments. As for urban street networks with intersections, we need DGHA with switching edges to describe traffic lights, which goes beyond our research in this paper.
We improve our previous work in [19,20,21] in the following aspects. Firstly, we introduce a partition of state subspace such that the combination modes of cells can be described by convex polytopes. Meanwhile, by adopting parallel composition algorithm of automata to obtain a hybrid automaton with high dimension over the whole road network. We finally deduce a Piecewise Affine Linear System (PWALS) of the overall road network. Based on the obtained PWALS, we further analyze the types and number of combination modes of cells and give a corresponding computational formula of the number of modes. An algorithm is developed on the basis of DGHA framework, and one can apply this algorithm to implement modelling of traffic networks with arbitrary topology structures and sizes.

3.1. The Dual Digraph Description of Road Network

In most literature, a road network is modeled as a digraph according to its original structure [22], i.e., each intersection or bottleneck point is represented as a vertex and the road section between adjacent intersections or bottleneck points is represented as the directed edge. However, in contrast to this traditional methodology, in our DGHA model, we adopt the dual digraph of the original structure of the road network. More precisely, each road section connecting two adjacent intersections or bottleneck points is partitioned into one or several segments and these segments, which are called cells, are modeled as vertices of the digraph. Two vertices are connected by an directed edge if there is traffic flow transition between them along the arrowed direction.
In our opinion, the dual digraph description has at least three advantages. Firstly, the network is easily extended by adding vertices and edges of the digraph. Secondly, weights describing dynamic procedures can be added in vertices as customary as in the weight digraph theory. Finally, traffic flow directions can be clearly indicated by the directed edges of the digraph and transited flows can be expediently described by dynamic weights and switchings, and especially, the edge switchings between vertices can describe the traffic light signals.
In this paper, we use fully weighted dynamic digraph G = V , E , H , Φ to describe urban freeway network, where V = 1 , 2 , , N denotes the set of all partitioned road segments of a given road network, E = e i j : i , j V denotes the set of all the directed edges indicating the transition of traffic flow, and the set H of automata and the set Φ of edge weighted functions will be illustrated later. Compared with the traditional modeling method, directed edges clearly indicate the transition of traffic flows and their weights express the traffic volumes transited during a pre-specified time interval. For the description of urban street network with intersections, we can express the switchings of traffic lights, instead of the fixed set E, by adopting a switching edge set E ( t ) .

3.2. Dynamics of Traffic Flow in Road Segment

To the best of our knowledge, the CTM [6,7] can well describe traffic flow, where the triangular fundamental diagram or the trapezoid fundamental diagram is adopted to approximatively express flow-density relationship [17,18,22]. In this paper, we use traffic density ρ i to express the traffic state of road segment i, and also adopt the triangular fundamental diagram shown in Figure 3 for road segment i, where C i is the traffic capacity (vph), V i is the free flow speed (mph), W i is the traffic wave speed (mph), ρ i 0 is the critical density (vpm), and ρ i m is the maximum/jam density (vpm).
To express the dynamics of vertex i, we use the hybrid automaton H i = S i , Init i , U i , f i , δ i , In v i , G i , Y i , g i , whose entries are illustrated as follows:
  • S i = Q i × X i , where Q i = F , C , “ F = Free flow , C = Congested flow ”, and the continuous state space is X i = R since only density ρ i is adopted as continuous state variable;
  • Init i S i can be given arbitrarily;
  • The input variable u i represents the changed flow of road segment i during the time period t T , ( t + 1 ) T , where T is the sample time period. Thus we have the input space U i = R ;
  • The vector field f i : S i × U i X i is described by the following equation:
    ρ i + = ρ i + a i u i .
    where a i = T / L i and L i is the length of the road segment (cell) i. Here and henceforth for simplicity we omit the time t in various variables and use ρ i + to express ρ i ( t + 1 ) ;
  • The transition map δ i is defined according to the guard conditions, i.e., a transition happens if and only if the corresponding invariant set is damaged and the guard condition is satisfied; The reset map R i is identical and thus is omitted;
  • The invariant sets are Inv i F = ρ i : 0 ρ i ρ i 0 and Inv i C = ρ i : ρ i 0 < ρ i ρ i m ;
  • The guard conditions are G i ( F , C ) = ρ i : ρ i > ρ i 0 and G i ( C , F ) = ρ i : ρ i ρ i 0 ;
  • The output space is Y i = R 2 ;
  • The output function g i : S i Y i is defined as y i = g i ( q i , x i ) = [ r i , s i ] T , where r i = min C i , W i ρ i m ρ i is the flow that can be received by cell i over the interval [ t T , t + 1 T ] ), s i = min V i ρ i , C i is the flow that can be supplied by cell i over the interval [ t T , t + 1 T ] ), i.e.,
    y i = g i F , ρ i = C i V i ρ i ,
    y i = g i C , ρ i = W i ρ i m ρ i C i .

3.3. Multi-Mode Description of Transited Traffic Flows

The edge weight function ϕ i j : Y i × Y j R in the set Φ = ϕ i j : e i j E expresses traffic flow practically transited from upstream cell i to downstream cell j. According to the rule of the CTM, the transited traffic flow is defined by a nonlinear function
ϕ i j y i , y j = min s i , r j ,
where r j is the flow that can be received by cell j and s i is the flow that can be supplied by cell i, both during the time interval t T , t + 1 T .
In what follows, we give a two-mode switching description of linear functions via two modes “D = Downward wave, U = Upward wave”:
D : ϕ i j y i , y j = s i , if s i r j , U : ϕ i j y i , y j = r j , if s i > r j .
Remark 5.
Two modes “D” and “U” describe the relation of the traffic flow between the road segments i and j. “D” mode means that the road segment j can release the traffic flow from i, while “U” mode indicates that the road segment j blocks the traffic flow from i.

3.4. Composition for DGHA Model of Traffic Flow Network

In literature [56], used a combination of LWR model with Godunov based numerical solutions [53] to solve the Riemann problem where the initial condition is a piecewise constant function with two values ρ and ρ r for the upstream (left) and downstream (right) densities. Either a shockwave or a rarefaction wave originates from the junction of the two densities. A shockwave develops if f ρ > f ρ r . Similarly, a rarefaction develops if f ρ < f ρ r . However, in the paper the hybrid switching automata approach proposed is to solve the traffic characteristics and system evolution by considering a Riemann problem with a piecewise constant initial condition.
Now we discuss the composition rules of the vertex automata affected by the edge weight functions described in Section 3.2 and Section 3.3 above. That is, we combine the dynamics of traffic flows in road segments with the transited traffic flows described by the two-mode switching linear functions to deduce the DGHA model of traffic flow network.
We assume that the merging interrelationship of flows in segments l Pre ( i ) and dividing interrelationship of flows in segments j Post ( i ) have been separated according to the rules via merging or dividing ratio [6,8], so that we can independently deal with the transited flows between any two segments ( l , i ) or ( i , j ) . In other words, assuming that ζ l i and η i j are the merging and dividing ratios respectively, the Equation (7) is modified by the merging and dividing rules
ϕ l i y l , y i = min s l , ζ l i r i ,
ϕ i j y i , y j = min η i j s i , r j .
Thus, according to the composition rules of hybrid automata in Equation (2), the function u i in Equation (4) is given as follows:
u i = l Pre ( i ) min s l , ζ l i r i j Post ( i ) min η i j s i , r j .
Especially, when only linear connections exist, e.g., in a urban ring freeway, we have Pre i = l and Post i = j in Equation (11), and thus ζ l i = 1 and η i j = 1 , which is the case we assume for simplicity in the following statement.
Now we use the way of modes in Equation (8) to express the Equation (11) into a switching procedure between a series of linear equations. In other words, we analyze the modes of vertex i in connected network by combining discrete states Q i and modes of all the edges connected with vertex i. It is obvious that these modes depend on the states of the vertices in N ¯ i , and each mode can be expressed by a group of inequalities of state variables ρ k , k N ¯ i . Thus, we use a partition of state subspace Θ i = k N ¯ i X k to express these modes so that the Equation (4) with Equation (11) is expressed by a piecewise affine linear system (PWALS).
For the special case of Pre i = l and Post i = j , Equation (11) into Equation (4) gives the following equation:
ρ i + = ρ i + a i ( min s l , r i min s i , r j )
In order to Equation (12) into a PWALS, we partition the subspace Θ i = { ρ ¯ i = ( ρ l , ρ i , ρ j ) T R 3 } into a group of subsets defined by convex polytopes Θ i k = { ρ ¯ i : H i k ρ ¯ i h i k } , and thus in each of the subsets the flows q l i = min { s l , r i } and q i j = min { s i , r j } can be expressed by linear relations.
First of all, we show that the cell i generally possesses seven modes with relation to upstream l shown in Table 1, where α l i = min { C i / V l , ρ l 0 } and β l i = max { ρ i 0 , ρ i m C l / W i } .
As seen in Table 1, the practical modes depend on the traffic capacities C l , C i . In fact, we can list the following three cases: C l > C i , C l < C i and C l = C i . In the case of C l > C i , there exist six combined modes “FDF, FUF, CUF, CUC, FUC, FDC”; In the case of C l < C i , there exist six combined modes “FDF, CDF, CDC, CUC, FUC, FDC”; In the case of C l = C i , mode “FUF” (resp. “CDC”) will be merged into mode “FDF” (resp. “CUC”) and thus there exist only five combined modes “FDF, CDF, CUC, FUC, FDC”, where “FDF” and “FDC” indicate free flow, “CUC” and “FUC” are congested flow, and “CDF” is saturated flow. Similarly, one can analyze the modes of i with respect to downstream j.
Next, we combine the modes of i with relation to l and j to obtain the modes of cell i in the connected network by using the corresponding subsets Θ i k . In this case, Equation (12) is expressed into a piecewise affine linear system with the following twelve linear subsystems:
ρ i + = ( 1 a i V i ) ρ i + ( a i V l ) ρ l , if FDFDF : H i 1 ρ ¯ i h i 1 ; or if FDFDC : H i 2 ρ ¯ i h i 2 ; ρ i + = ρ i + ( a i V l ) ρ l + ( a i W j ) ρ j a i W j ρ j m , if FDFUC : H i 3 ρ ¯ i h i 3 ; or if FDCUC : H i 4 ρ ¯ i h i 4 ; ρ i + = ρ i + ( a i V l ) ρ l a i C , if FDCDF : H i 5 ρ ¯ i h i 5 ; ρ i + = ( 1 a i V i ) ρ i + a i C , if CDFDF : H i 6 ρ ¯ i h i 6 ; or if CDFDC : H i 7 ρ ¯ i h i 7 ; ρ i + = ρ i + ( a i W j ) ρ j + a i ( C W j ρ j m ) , if CDFUC : H i 8 ρ ¯ i h i 8 ; ρ i + = ( 1 a i W i ) ρ i + a i ( W i ρ i m C ) , if CUCDF : H i 9 ρ ¯ i h i 9 ; or if FUCDF : H i ( 10 ) ρ ¯ i h i ( 10 ) ; ρ i + = ( 1 a i W i ) ρ i + ( a i W j ) ρ j + a i ( W i ρ i m W j ρ j m ) , if FUCUC : H i ( 11 ) ρ ¯ i h i ( 11 ) ; or if CUCUC : H i ( 12 ) ρ ¯ i h i ( 12 ) .
where the matrices H i k , h i k , k = 1 , , 12 are listed in Appendix A to express Θ i k = { ρ ¯ i : H i k ρ ¯ i h i k } . We add Θ i 0 = Θ i \ k = 1 12 Θ i k , and thus have Θ i = k = 0 12 Θ i k .
Similarly, we can deal with all the cells in an urban ring freeway which is partitioned into N cells. Finally, we can obtain all modes of the freeway network by finding intersections of all the convex polytopes Θ i k , i = 1 , , N , k = 1 , , 12 , in the state space R N .
Let S = S ( N ) be the number of all the combined modes of the freeway network with N cells, and D s R N , s = 1 , , S , be the convex polytopes describing the combined modes. We use x = [ ρ 1 , , ρ N ] T R N to denote the traffic density vector of all the cells, use the vector u R M to represent the traffic demand from on-ramps of the freeway network, and assume that the dividing ratio of the cell i with an off-ramp to the downstream cell j is η i j . Thus traffic flow of the freeway network can be modeled by a PWALS with S subsystems and a switching function σ : [ 0 , + ) { 1 , 2 , , S } among these subsystems:
x ( t + 1 ) = A σ ( t ) x ( t ) + B σ ( t ) u ( t ) + F σ ( t ) ,
where the switching function σ ( t ) is determined by the convex polytopes D s , i.e., σ ( t ) = s if and only if x ( t ) D s ; for s = 1 , , S , vector F s , matrices A s and B s consist of the parameters in the fundamental diagrams and the dividing ratios of all the cells.

3.5. Analysis of the Number of Combined Modes

As the basis of other work, such as density estimation and on-ramp metering, it is essential to analyze the complexity of combined modes. Now we take an urban freeway with N partitioned cells as an example to investigate the combined mode number S ( N ) .
We prove that S ( N ) satisfies the following relation:
S ( N ) = 2 S ( N 1 ) + S ( N 2 ) , N 3 ,
with initial conditions S ( 1 ) = 2 and S ( 2 ) = 5 .
The proof is as follows. As seen in Figure 4, if the current cell is in mode “F”, there are only three possible modes combined with the downstream cells, i.e., “FDF”, “FUC”, “FDC”. Otherwise, there exist only two possible combination modes when the current cell is in mode “C”, i.e., “CDF”, “CUC”. Hence, we have the following equation:
S ( N ) = 3 S F N 1 + 2 S C N 1
where S F N 1 and S C N 1 are the mode numbers of the first N 1 cells when the ( N 1 )th cell is in mode “F” and “C”, respectively.
On the other hand, we have S F N 1 = S ( N 2 ) since “F” can appear after any modes. Moreover, because of S F N 1 + S C N 1 = S ( N 1 ) , we get S C N 1 = S ( N 1 ) S F N 1 . Substituting these two expressions into Equation (16), we obtain the following equation:
S ( N ) = 3 S F ( N 1 ) + 2 S C ( N 1 ) = 3 S ( N 2 ) + 2 S ( N 1 ) S ( N 2 ) = 2 S ( N 1 ) + S ( N 2 ) ,
and thus verify the Equation (15).
Although we focus on the analysis of urban freeway networks, we can use similar method to investigate other types of road networks. However, since these combination modes are affected by various factors, such as traffic lights, the types and the number of merge connections and diverge connections, the situation may become more complex. On the other hand, we should point out that although we list all the combination modes for a given road network theoretically, further studies are needed to check whether all the modes exist practically and how these modes change during traffic flow evolution.
In the literature [34], the number of modes is 2 N for a road segment with N cells. They only considered the modes in the vertices and ignored the multi-mode switchings of the edges. In our DGHA model, the modes of vertices and edges are both considered. Consequently, the modes are more than the those in [34].

4. State Observer Design

State observers have been considered as effective tools of state estimation in the control theory, and various types of state observers have been developed and widely applied in practical fields. In this section, we choose the Luenberger-type switching observer in [43] to estimate traffic densities based on the PWALS Equation (14) for a freeway network divided into N cells.
We suppose that M fixed or floating sensors are installed on some road segments and aim to estimate the densities on all the segments of the whole freeway network via the PWALS Equation (14) and the traffic data obtained from the sensors.
Let y ( t ) R M be the measured output vector of the sensors and assume that it has the following relation to the densities:
y ( t ) = C x ( t )
where C is the output matrix depending on the type, number and location of the sensors. We also assume that A s , C , s = 1 , 2 , S , are observable or detectable.
In this paper, we assume the switching rule of the observer is the same with that of the PWALS Equation (14). Thus, based on the Equation (14) and the measured output Equation (23), we construct the following piecewise affine linear observer (Figure 5):
x ^ ( t + 1 ) = A σ ( t ) x ^ ( t ) + B σ ( t ) u ( t ) + F σ ( t ) + K σ ( t ) y ( t ) C x ^ ( t ) ,
where x ^ t R N is the estimation of the state x ( t ) , K s , s = 1 , , S , are the observer gain matrices which will be designed below, the switching function σ : [ 0 , + ) { 1 , 2 , , S } is the same with that of the PWALS Equation (14) and thus is determined by the convex polytopes D s , s = 1 , , S , i.e., σ ( t ) = s if and only if x ( t ) D s .
Therefore, the state estimation error e ( t ) = x ( t ) x ^ ( t ) can be described by the following switched linear system:
e ( t + 1 ) = A σ ( t ) K σ ( t ) C e ( t ) .
In order to ensure lim t e ( t ) = 0 for any switching signal σ ( t ) , the switched linear Equation (19) must be asymptotically stable for any switching signal σ ( t ) . We apply the common quadratic Lyapunov function approach for the following subsystems to compute the observer gain matrices K s :
e ( t + 1 ) = ( A s K s C ) e ( t ) , s = 1 , , S .
First of all, we state two well-known lemmas as follows.
Lemma 1.
Assume that A s , C , s = 1 , 2 , S , are observable or detectable. If there exists a positive definite symmetric matrix P satisfying the following Lyapunov inequalities:
( A s K s C ) T P ( A s K s C ) P s < 0 , s = 1 , , S
the gain matrices K s , s = 1 , , S , ensure that, for any switching signal σ ( t ) , the estimation error of the observer Equation (18) exponentially converge to zero.
Lemma 2.
Let X s = P K s . Then the inequalities Equation (21) are equivalent to the following LMIs with respect to P and X s :
P ( P A s X s C ) T P A s X s C P > 0 , s = 1 , , S .
Thus, we can find the gain matrices K s = P 1 X s , s = 1 , , S by solving the LMIs Equation (22).
In fact, the condition in Equation (21) or (22) provides a sufficient and necessary one for the existence of the observer Equation (18). In other words, if the condition Equation (21) or (22) is not satisfied, one has to change the structure of the observer or even to find another type of state estimation approach. The same case happens in other estimation methods such as Kalman filter, where one cannot always find the Kalman gain matrix.

5. Application to Beijing Third Ring Freeway

5.1. Test Data and Parameter Settings

Beijing third ring freeway (Figure 6) is approximately 48 km long and includes 62 on-ramps and 62 off-ramps. We apply the proposed model and state observer to the outer ring freeway. According to consistency rule of cells [21], the outer ring freeway can be partitioned into 128 cells. As shown in Section 3, if a centralized Equation (15) is used for the freeway with 128 cells, there exist 8.44342 × 10 48 combination modes, which results in high computational complexity. Therefore, in this paper, in order to clearly interpret the principle of the proposed method and avoid computational complexity, we take into account a simplified version of the outer ring freeway. For the practical application we will consider a decentralized model and a distributed observer in the future research.
The simplified outer ring freeway (Figure 7) contains 20 cells, 4 on-ramps and 4 off-ramps. Each cell has three lanes and the length of the cells is listed in Table 2. The parameters of all cells are given in Table 3.
Moreover, since Modelica language [20,57] has the advantage of modeling hybrid systems, we choose the open source software OpenModelica to develop a program library to implement the traffic flow model and state estimation.
For the experimental data, we collected 6 h (from 6:00 a.m. to 12:00 a.m.) traffic flow data of 4 on-ramps, 4 off-ramps and traffic densities of 4 road segments of the outer ring freeway on 19 May 2015. Then we used a microscopic model of the outer ring freeway built by Paramics [58] to simulate the traffic flow from 6:00 a.m. to 12:00 a.m. and collect the densities (we call them the simulated densities) by virtual sensors in all the segments, where the real traffic flows of 4 on-ramps and 4 off-ramps were used as the traffic demand and to determine the dividing ratios. The real initial densities of the 4 road segments were used as the initial traffic state of all the cells. In simulation procedure, by setting the dividing ratio, and adjusting the whole microscopic simulator, the modes of traffic flow change in the following sequence order: 20 F 16 F 4 C 12 F 2 C 2 F 4 C 6 F 3 C 3 F 2 C 2 F 4 C 6 F 3 C 3 F 8 C 16 C 7 F 1 C 1 F 3 C 1 F 3 C 3 F 1 C , where k F or k C denotes continuous k cells in mode F or C .
The corresponding matrices A s = [ a i , j ] and vectors F s = [ f i ] , s = 1 , 2 20 , are listed in Appendix B. For the outer ring freeway with 4 on-ramps the matrix B s = [ b i , j ] is a constant matrix with the entries b 3 , 1 = 0 . 0149 , b 8 , 2 = 0.0116 , b 13 , 3 = 0.0109 , b 18 , 4 = 0.0238 , and b i , j = 0 for the other entries in B s . Here the sample time period T = 5 s such that the condition V i T L i , i = 1 , 2 20 , is satisfied.
Nowadays various sensor technologies can be applied to collect traffic data, including floating vehicle, loop detector, microwave and video. Different source data can provide different measurement outputs and thus one gets different types of matrix C in the measurement output Equation (23). One can even utilize the fused data or multi-dimensional data of several types of sensors at the same time.
However, in this paper we used the loop detector to collect traffic densities. In the example, the virtual loop detectors in Paramics installed in the cells 1, 3, 5, 7, 9, 11, 13, 15, 17, 19 were used to collect traffic densities (see Figure 7). Hence, corresponding to these detectors, C = [ c i , j ] is a 10 × 20 matrix with entries c 1 , 1 = c 2 , 3 = c 3 , 5 = c 4 , 7 = c 5 , 9 = c 6 , 11 = c 7 , 13 = c 8 , 15 = c 9 , 17 = c 10 , 19 = 1 and c i , j = 0 for the other i , j . Moreover, one can verify that the pair ( A s , C ) is observable for all s = 1 , , 20 .

5.2. Analysis of Simulation Results

The simulation results for the whole simplified outer ring freeway are shown in Figure 8 in the way of time-space diagram, where Figure 8a shows the real traffic densities from 6:00 a.m. to 12:00 a.m., and Figure 8b presents the estimated densities by the observer. It is clear that some segments with red color represent a real congestion state. Especially, heavy traffic jams happen in the morning rush hour from about 7:00 a.m. to 10:00 a.m. and the congestion time periods are different for different road segments. Error curve(see Figure 9) shows that the error converges to zero at about 200 s, the speed is acceptable and the gain matrices are appropriate.
One can further compare estimated densities with the simulated ones in all road segments. For example, the cells 14 and 20 are presented to exhibit more in detail the the density comparison (see Figure 10). From Figure 8 and Figure 10, it is clear that the proposed modeling approach and the observer are feasible, and the estimated densities well approximate to the simulated ones.
Table 4 shows the performance indicator of the state observer by using the mean square error (MSE) which is defined by the Equation (23)
MSE = i = 1 n ρ ^ i ρ i 2 n
From the MSE of each cell, we can further compute the mean value of 20 cells is approximately 11.34 % . The results indicate that the designed model-based switched state observer has a good estimated performance.

6. Conclusions and Future Work

In the paper, a modeling framework named Dynamic Graph Hybrid Automata (DGHA) has been proposed by combining a dynamic graph with hybrid automata. Then this framework was applied to model traffic flow over an urban freeway network by embedding the Cell Transmission Model (CTM) into the DGHA.
In the modeling procedure, the dual digraph of road network structure was adopted to describe the road topology, hybrid automata were used to describe multi-modes of dynamic densities in road segments, and the nonlinear expressions of the transmitted traffic flow between two road segments were transformed into piecewise linear functions in terms of multi-mode switchings. By using a combination algorithm for the dynamics of traffic flow, the modeling procedure is modularized and rule-based, and thus is easily-extensible. Thus it can be used to model an urban freeway network with arbitrary topology structures and sizes. Furthermore, mode types and number in the model of the whole freeway network were analyzed, and, consequently, a Piecewise Affine Linear System (PWALS) model of urban freeway was deduced. Finally, based on the PWALS model, a switched state observer was designed to estimate the traffic densities of the urban freeway, where a set of observer gain matrices were computed by using the Lyapunov function approach. The proposed modeling approach and the observer were validated to be feasible by the simplified example of Beijing third ring freeway.
However, there are some drawbacks which require further improvements in the future research. First, we adopt centralized model and observer, which, as shown above, results in huge number of combined modes and computational difficulty. So from the view of practical application, decentralized model and decentralized observer are required to reduce computational complexity.
Second, the assumption that the switching function of the observer is the same with that of the model means that the switching of the observer is driven by the model states. However, in practice it may be difficult to implement and one needs to identify mode switchings by designing a discrete event observer, which will result in unsynchronized mode switchings between the practical object and the observer and thus bring about larger estimation error. So, in this case, another type of switched observer is required to be design for the state estimation.
Third, in the application example, we took the simplified version of Beijing third ring road, and used the incompletely collected traffic data and the simulation data obtained by Paramics. That means we need real experiments for the practical application to further verify our proposed approach.

Acknowledgments

The work was supported by the National Natural Science Foundation of China (Grant No. 61273006, 61573030, 61511130044).

Author Contributions

Yangzhou Chen and Yuqi Guo modeled the urban freeway network; Yangzhou Chen and Yuqi Guo conceived and designed the experiments; Yuqi Guo and Ying Wang performed the experiments and analyzed the data; Yangzhou Chen and Yuqi Guo wrote the paper.

Conflicts of Interest

The authors declare no conflict if interest.

Appendix A

H i 1 = 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1 , h i 1 = 0 ρ l 0 0 ρ i 0 0 ρ j 0
H i 2 = 1 0 0 1 0 0 0 1 0 0 0 1 0 V i W j , h i 2 = 0 ρ l 0 0 ρ j 0 W j ρ j m
H i 3 = 1 0 0 1 0 0 0 1 0 0 0 1 0 V i W j , h i 3 = 0 ρ l 0 ρ i 0 ρ j m W j ρ j m
H i 4 = 1 0 0 0 1 0 V l W i 0 0 0 1 0 0 1 , h i 4 = 0 ρ i 0 W i ρ i m ρ j 0 ρ j m
H i 5 = 1 0 0 0 1 0 V l W i 0 0 0 1 0 0 1 , h i 5 = 0 ρ i 0 W i ρ i m 0 ρ j 0
H i 6 = 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1 , h i 6 = ρ l 0 ρ l m 0 ρ i 0 0 ρ j 0
H i 7 = 1 0 0 1 0 0 0 1 0 0 0 1 0 V i W j , h i 7 = ρ l 0 ρ l m 0 ρ j 0 W j ρ j m
H i 8 = 1 0 0 1 0 0 0 1 0 0 0 1 0 V i W j , h i 8 = ρ l 0 ρ l m ρ i 0 ρ j m W j ρ j m
H i 9 = 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1 , h i 9 = ρ l 0 ρ l m ρ i 0 ρ i m 0 ρ j 0
H i 10 = 1 0 0 0 1 0 V l W i 0 0 0 1 0 0 1 , h i 10 = ρ l 0 ρ i m W i ρ i m 0 ρ j 0
H i 11 = 1 0 0 0 1 0 V l W i 0 0 0 1 0 0 1 , h i 11 = ρ l 0 ρ i m W i ρ i m ρ j 0 ρ j m
H i 12 = 1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1 , h i 12 = ρ l 0 ρ l m ρ i 0 ρ i m ρ j 0 ρ j m

Appendix B

Mode 1: FFFFFFFFFFFFFFFFFFFF (20F)
a 1 , 1 = 0.2491 , a 1 , 20 = 0.7509 , a 2 , 1 = 0.1558 , a 2 , 2 = 0.8442 , a 3 , 2 = 0.2481 , a 3 , 3 = 0.7519 , a 4 , 3 = 0.6040 , a 4 , 4 = 0.3960 , a 5 , 4 = 0.1422 , a 5 , 5 = 0.8578 , a 6 , 5 = 0.5954 , a 6 , 6 = 0.4046 , a 7 , 6 = 0.2966 , a 7 , 7 = 0.7034 , a 8 , 7 = 0.1938 , a 8 , 8 = 0.8062 , a 9 , 8 = 0.1667 , a 9 , 9 = 0.8333 , a 10 , 9 = 0.2348 , a 10 , 10 = 0.7652 , a 11 , 10 = 0.2466 , a 11 , 11 = 0.7534 , a 12 , 11 = 0.2348 , a 12 , 12 = 0.7652 , a 13 , 12 = 0.1820 , a 13 , 13 = 0.8180 , a 14 , 13 = 0.3256 , a 14 , 14 = 0.6744 , a 15 , 14 = 0.2778 , a 15 , 15 = 0.7222 , a 16 , 15 = 0.2381 , a 16 , 16 = 0.7619 , a 17 , 16 = 0.2511 , a 17 , 17 = 0.7489 , a 18 , 17 = 0.3969 , a 18 , 18 = 0.6031 , a 19 , 18 = 0.2271 , a 19 , 19 = 0.7729 , a 20 , 19 = 0.1820 , a 20 , 20 = 0.8180 , and the other entries a i , j are zero. F = 0 .
Mode 2: FFFFFFFFFFFFFFFFCCCC (16F4C)
a 1 , 1 = 0.2491 , a 2 , 1 = 0.1558 , a 2 , 2 = 0.8442 , a 3 , 2 = 0.2481 , a 3 , 3 = 0.7519 , a 4 , 3 = 0.6040 , a 4 , 4 = 0.3960 , a 5 , 4 = 0.1422 , a 5 , 5 = 0.8578 , a 6 , 5 = 0.5954 , a 6 , 6 = 0.4046 , a 7 , 6 = 0.2966 , a 7 , 7 = 0.7034 , a 8 , 7 = 0.1938 , a 8 , 8 = 0.8062 , a 9 , 8 = 0.1667 , a 9 , 9 = 0.8333 , a 10 , 9 = 0.2348 , a 10 , 10 = 0.7652 , a 11 , 10 = 0.2466 , a 11 , 11 = 0.7534 , a 12 , 11 = 0.2348 , a 12 , 12 = 0.7652 , a 13 , 12 = 0.1820 , a 13 , 13 = 0.8180 , a 14 , 13 = 0.3256 , a 14 , 14 = 0.6744 , a 15 , 14 = 0.2778 , a 15 , 15 = 0.7222 , a 16 , 15 = 0.2381 , a 16 , 16 = 0.7619 , a 17 , 16 = 0.2511 , a 17 , 17 = 1 , a 17 , 18 = 0.0837 , a 18 , 18 = 0.8676 , a 18 , 19 = 0.1324 , a 19 , 19 = 0.9243 , a 19 , 20 = 0.0757 , a 20 , 20 = 0.9363 , and the other entries a i , j are zero.
f 1 = 0.0338 , f 17 = 0.0151 , f 20 = 0.0027 , and the other entries f i are zero.
Mode 3: FFFFFFFFFFFFCCFFCCCC (12F2C2F4C)
a 1 , 1 = 0.2491 , a 2 , 1 = 0.1558 , a 2 , 2 = 0.8442 , a 3 , 2 = 0.2481 , a 3 , 3 = 0.7519 , a 4 , 3 = 0.6040 , a 4 , 4 = 0.3960 , a 5 , 4 = 0.1422 , a 5 , 5 = 0.8578 , a 6 , 5 = 0.5954 , a 6 , 6 = 0.4046 , a 7 , 6 = 0.2966 , a 7 , 7 = 0.7034 , a 8 , 7 = 0.1938 , a 8 , 8 = 0.8062 , a 9 , 8 = 0.1667 , a 9 , 9 = 0.8333 , a 10 , 9 = 0.2348 , a 10 , 10 = 0.7652 , a 11 , 10 = 0.2466 , a 11 , 11 = 0.7534 , a 12 , 11 = 0.2348 , a 12 , 12 = 0.7652 , a 13 , 12 = 0.1820 , a 13 , 13 = 1 , a 13 , 14 = 0.0607 , a 14 , 14 = 0.8914 , a 15 , 15 = 0.7222 , a 16 , 15 = 0.2381 , a 16 , 16 = 0.7619 , a 17 , 16 = 0.2511 , a 17 , 17 = 1 , a 17 , 18 = 0.0837 , a 18 , 18 = 0.8676 , a 18 , 19 = 0.1324 , a 19 , 19 = 0.9243 , a 19 , 20 = 0.0757 , a 20 , 20 = 0.9363 , and the other entries a i , j are zero.
f 1 = 0.0338 , f 13 = 0.011 , f 14 = 0.0049 , f 15 = 0.0125 , f 17 = 0.0151 , f 20 = 0.0027 , and the other entries f i are zero.
Mode 4: FFFFFFCCCFFFCCFFCCCC (6F3C3F2C2F4C)
a 1 , 1 = 0.2491 , a 2 , 1 = 0.1558 , a 2 , 2 = 0.8442 , a 3 , 2 = 0.2481 , a 3 , 3 = 0.7519 , a 4 , 3 = 0.6040 , a 4 , 4 = 0.3960 , a 5 , 4 = 0.1422 , a 5 , 5 = 0.8578 , a 6 , 5 = 0.5954 , a 6 , 6 = 0.4046 , a 7 , 6 = 0.2966 , a 7 , 7 = 1 , a 7 , 8 = 0.0989 , a 8 , 8 = 0.9353 , a 8 , 9 = 0.0647 , a 9 , 9 = 0.9444 , a 10 , 10 = 0.7652 , a 11 , 10 = 0.2466 , a 11 , 11 = 0.7534 , a 12 , 11 = 0.2348 , a 12 , 12 = 0.7652 , a 13 , 12 = 0.1820 , a 13 , 13 = 1 , a 13 , 14 = 0.0607 , a 14 , 14 = 0.8914 , a 15 , 15 = 0.7222 , a 16 , 15 = 0.2381 , a 16 , 16 = 0.7619 , a 17 , 16 = 0.2511 , a 17 , 17 = 1 , a 17 , 18 = 0.0837 , a 18 , 18 = 0.8676 , a 18 , 19 = 0.1324 , a 19 , 19 = 0.9243 , a 19 , 20 = 0.0757 , a 20 , 20 = 0.9393 , and the other entries a i , j are zero.
f 1 = 0.0338 , f 7 = 0.018 , f 9 = 0.0025 , f 10 = 0.0106 , f 13 = 0.011 , f 14 = 0.0049 , f 15 = 0.0125 , f 17 = 0.0151 , f 20 = 0.0027 , and the other entries f i are zero.
Mode 5: FFFFFFCCCFFFCCCCCCCC (6F3C3F8C)
a 1 , 1 = 0.2491 , a 2 , 1 = 0.1558 , a 2 , 2 = 0.8442 , a 3 , 2 = 0.2481 , a 3 , 3 = 0.7519 , a 4 , 3 = 0.6040 , a 4 , 4 = 0.3960 , a 5 , 4 = 0.1422 , a 5 , 5 = 0.8578 , a 6 , 5 = 0.5954 , a 6 , 6 = 0.4046 , a 7 , 6 = 0.2966 , a 7 , 7 = 1 , a 7 , 8 = 0.0989 , a 8 , 8 = 0.9353 , a 8 , 9 = 0.0647 , a 9 , 9 = 0.9444 , a 10 , 10 = 0.7652 , a 11 , 10 = 0.2466 , a 11 , 11 = 0.7534 , a 12 , 11 = 0.2348 , a 12 , 12 = 0.7652 , a 13 , 12 = 0.1820 , a 13 , 13 = 1 , a 13 , 14 = 0.0607 , a 14 , 14 = 0.8914 , a 14 , 15 = 0.1086 , a 15 , 15 = 0.9073 , a 15 , 16 = 0.0927 , a 16 , 16 = 0.9206 , a 16 , 16 = 0.0794 , a 17 , 17 = 0.9163 , a 17 , 18 = 0.0837 , a 18 , 18 = 0.8676 , a 18 , 19 = 0.1324 , a 19 , 19 = 0.9243 , a 19 , 20 = 0.0757 , a 20 , 20 = 0.9393 , and the other entries a i , j are zero.
f 1 = 0.0338 , f 7 = 0.018 , f 9 = 0.0025 , f 10 = 0.0106 , f 13 = 0.011 , f 20 = 0.0027 , and the other entries f i are zero.
Mode 6: CCCCCCCCCCCCCCCCCCCC (16C)
a 1 , 1 = 0.7495 , a 1 , 2 = 0.2505 , a 2 , 2 = 0.9480 , a 2 , 3 = 0.0520 , a 3 , 3 = 0.9173 , a 3 , 4 = 0.0827 , a 4 , 4 = 0.7986 , a 4 , 5 = 0.2014 , a 5 , 5 = 0.9526 , a 5 , 6 = 0.0474 , a 6 , 6 = 0.8014 , a 6 , 7 = 0.1986 , a 7 , 7 = 0.9011 , a 7 , 8 = 0.0989 , a 8 , 8 = 0.9353 , a 8 , 9 = 0.0647 , a 9 , 9 = 0.9444 , a 9 , 10 = 0.0556 , a 10 , 10 = 0.9217 , a 10 , 11 = 0.0783 , a 11 , 11 = 0.9178 , a 11 , 12 = 0.0822 , a 12 , 12 = 0.9217 , a 12 , 13 = 0.0783 , a 13 , 13 = 0.9393 , a 13 , 14 = 0.0607 , a 14 , 14 = 0.8914 , a 14 , 15 = 0.1086 , a 15 , 15 = 0.9073 , a 15 , 16 = 0.0927 , a 16 , 16 = 0.9206 , a 16 , 17 = 0.0794 , a 17 , 17 = 0.9163 , a 17 , 18 = 0.0837 , a 18 , 18 = 0.8676 , a 18 , 19 = 0.1324 , a 19 , 19 = 0.9243 , a 19 , 20 = 0.0757 , a 20 , 20 = 0.9393 , and the other entries a i , j are zero.
f 4 = 0.7347 , f 6 = 0.8036 , f 9 = 0.2167 , f 11 = 0.2630 , f 14 = 0.3771 , f 16 = 0.1885 , f 19 = 0.0719 , and the other entries f i are zero.
Mode 7: FFFFFFFCFCCCFCCCFFFC
(7F1C1F3C1F3C3F1C)
a 1 , 1 = 0.2491 , a 2 , 1 = 0.1558 , a 2 , 2 = 0.8442 , a 3 , 2 = 0.2481 , a 3 , 3 = 0.7519 , a 4 , 3 = 0.6040 , a 4 , 4 = 0.3960 , a 5 , 4 = 0.1422 , a 5 , 5 = 0.8578 , a 6 , 5 = 0.5954 , a 6 , 6 = 0.4046 , a 7 , 6 = 0.2966 , a 7 , 7 = 1 , a 7 , 8 = 0.0989 , a 8 , 8 = 0.9353 , a 8 , 9 = 0.0647 , a 9 , 9 = 1 , a 9 , 10 = 0.0556 , a 10 , 10 = 0.9217 , a 10 , 11 = 0.0783 , a 11 , 11 = 0.9178 , a 11 , 12 = 0.0822 , a 12 , 12 = 0.9217 , a 13 , 13 = 1 , a 13 , 14 = 0.0607 , a 14 , 14 = 0.8914 , a 14 , 15 = 0.1086 , a 15 , 15 = 0.9073 , a 15 , 16 = 0.0927 , a 16 , 16 = 0.9206 , a 17 , 16 = 0.9206 , a 17 , 17 = 0.7489 , a 18 , 17 = 0.3969 , a 18 , 18 = 0.6031 , a 19 , 18 = 0.2271 , a 19 , 19 = 1 , a 19 , 20 = 0.0757 , a 20 , 20 = 0.9393 , and the other entries a i , j are zero.
f 1 = 0.0338 , f 7 = 0.018 , f 8 = 0.0029 , f 9 = 0.0025 , f 12 = 0.0035 , f 12 = 0.0027 , f 16 = 0.0036 , f 17 = 0.0113 , f 13 = 0.0035 , f 19 = 0.0136 , f 20 = 0.0027 , and the other entries f i are zero.

References

  1. Newman, M.E.J. The structure and function of complex networks. SIAM Rev. 2003, 45, 167–256. [Google Scholar] [CrossRef]
  2. Barrat, A.; Barthelemy, M.; Vespignani, A. Dynamical Processes on Complex Networks; Cambridge University Press: New York, NY, USA, 2008; pp. 242–264. [Google Scholar]
  3. Lighthill, M.J.; Whitham, G.B. On kinematic waves. II. A theory of traffic flow on long crowded roads. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci. 1955, 229, 317–345. [Google Scholar] [CrossRef]
  4. Richards, P.I. Shock waves on the highway. Oper. Res. 1956, 4, 42–51. [Google Scholar] [CrossRef]
  5. Payne, H.J. Models of Freeway Traffic and Control; SIMULATION COUNCILS, INC.: San Diego, CA, USA, 1971; pp. 51–61. [Google Scholar]
  6. Daganzo, C.F. The cell transmission model: A dynamic representation of highway traffic consistent with the hydrodynamic theory. Transp. Res. Part B Methodol. 1994, 28, 269–287. [Google Scholar]
  7. Daganzo, C.F. The cell transmission model, part II: Network traffic. Transp. Res. Part B Methodol. 1995, 29, 79–93. [Google Scholar] [CrossRef]
  8. Lebacque, J.P. The Godunov scheme and what it means for first order traffic flow models. In Proceedings of the 13th Internaional Symposium on Transportation and Traffic Theory, Lyon, France, 24–26 July 1996; pp. 647–677. [Google Scholar]
  9. Flötteröd, G.; Nagel, K. Some practical extensions to the Cell Transmission Model. In Proceedings of the 2005 8th IEEE international Conference on Intelligent Transportation Systems(ITSC), Vienna, Austria, 13–16 September 2005; pp. 172–177. [Google Scholar]
  10. Gomes, G.; Horowitz, R. Globally optimal solutions to the on-ramp metering problem-Part I. In Proceedings of the 7th IEEE international Conference on Intelligent Transportation Systems(ITSC), Washington, DC, USA, 3–6 October 2004; pp. 509–514. [Google Scholar]
  11. Gomes, G.; Horowitz, R. Optimal freeway ramp metering using the asymmetric cell transmission model. Transp. Res. Part C Emerg. Technol. 2006, 14, 244–262. [Google Scholar] [CrossRef]
  12. Gomes, G.; Horowitz, R. Globally optimal solutions to the on-ramp metering problem-Part II. In Proceedings of the 7th IEEE international Conference on Intelligent Transportation Systems(ITSC), Washington, DC, USA, 3–6 October 2004; pp. 515–520. [Google Scholar]
  13. Lo, H.K. A cell-based traffic control formulation: strategies and benefits of dynamic timing plans. Transp. Sci. 2001, 35, 148–164. [Google Scholar] [CrossRef]
  14. Lin, W.H.; Wang, C. An enhanced 0-1 mixed-integer LP formulation for traffic signal control. IEEE Trans. Intell. Transp. Syst. (ITS) 2004, 5, 238–245. [Google Scholar] [CrossRef]
  15. Sumalee, A.; Zhong, R.; Pan, T.; Lam, W.H.K. Stochastic cell transmission model (SCTM): A stochastic dynamic traffic model for traffic state surveillance and assignment. Transp. Res. Part B Methodol. 2011, 45, 507–533. [Google Scholar]
  16. Canudas-de-Wit, C.; Ojeda, L.L.; Kibangou, A.Y. Graph constrained-CTM observer design for the Grenoble south ring. In Proceedings of the 13th IFAC Symposium on Control in Transportation Systems, Sofia, Bulgaria, 12–14 September 2012; Volume 45, pp. 197–202. [Google Scholar]
  17. Muñoz, L.; Sun, X.; Horowitz, R.; Luis, A. Piecewise-linearized cell transmission model and parameter calibration methodology. Transp. Res. Rec. 2006, 1965, 183–191. [Google Scholar]
  18. Muñoz, L.; Sun, X.; Horowitz, R.; Luis, A. Traffic density estimation with the cell transmission model. In Proceedings of the 2003 American Control Conference (ACC), Denver, CO, USA, 4–6 June 2003; Volume 5, pp. 3750–3755. [Google Scholar]
  19. Chen, Y.; He, Z.; Shi, J.; Han, X. Dynamic graph hybrid system: A modeling method for complex networks with application to urban traffic. In Proceedings of the 10th Intelligent Control and Automation (WCICA), Beijing, China, 6–8 July 2012; pp. 1864–1869. [Google Scholar]
  20. Chen, Y.; Li, W.; Guo, Y.; Wu, Y. Dynamic graph hybrid automata: A modeling method for traffic network. In Proceedings of the 18th IEEE international Conference on Intelligent Transportation Systems(ITSC), Canary Islands, Spain, 15–18 September 2015; pp. 1396–1401. [Google Scholar]
  21. Chen, Y.; Guo, Y.; Wang, Y.; Li, W. Modeling freeway network by using dynamic graph hybrid automata and estimating its states by designing state observer. In Proceedings of the Chinese Automation Congress (CAC), Wuhan, China, 27–29 November 2015; pp. 237–242. [Google Scholar]
  22. Sun, X.; Muñoz, L.; Horowitz, R. Mixture Kalman filter based highway congestion mode and vehicle density estimator and its application. In Proceedings of the American Control Conference (ACC), Boston, MA, USA, 30 June–2 July 2004; Volume 3, pp. 2098–2103. [Google Scholar]
  23. Harary, F.; Gupta, G. Dynamic graph models. Math. Comput. Model. 1997, 25, 79–87. [Google Scholar]
  24. Lunze, J.; Lamnabhi-Lagarrigue, F. Handbook of Hybrid Systems Control: Theory, Tools, Applications; Cambridge University Press: New York, NY, USA, 2009. [Google Scholar]
  25. Henzinger, T.A. The theory of hybrid automata. Verification of Digital and Hybrid Systems; Springer: Berlin/Heidelberg, Germany, 2000; pp. 278–292. [Google Scholar]
  26. Varaiya, P. Smart cars on smart roads: problems of control. IEEE Trans. Autom. Control 1993, 38, 195–207. [Google Scholar] [CrossRef]
  27. Lygeros, J.; Godbole, D.N.; Sastry, S. Verified hybrid controllers for automated vehicles. IEEE Trans. Autom. Control 1998, 43, 522–539. [Google Scholar]
  28. Lei, J.; Ozguner, U. Decentralized hybrid intersection control. In Proceedings of the 40th IEEE Conference on Decision and Control, Orlando, FL, USA, 4–7 December 2001; Volume 2, pp. 1237–1242. [Google Scholar]
  29. Zhao, X.; Chen, Y. Traffic light control method for a single intersection based on hybrid systems. In Proceedings of the 2003 IEEE Intelligent Transportation Systems, Shanghai, China, 12–15 October 2003; Volume 2, pp. 1105–1109. [Google Scholar]
  30. Zhao, X.; Chen, Y.; Li, Z.; Chen, Y. An optimal control method for hybrid systems based on Q-learning for an intersection traffic signal control. Chin. High Technol. Lett. 2007, 17, 498–502. [Google Scholar]
  31. Jiang, G.; Chen, Y.; Zhang, L.; Li, H.F. Modeling and Reachability Analysis for Single Intersection Based on Rectangular Hybrid Automata. J. Transp. Syst. Eng. Inf. Technol. 2009, 9, 120–126. [Google Scholar]
  32. Grahlmann, B. Combining finite automata, parallel programs and SDL using petri nets. In Proceedings of the International Conference on Tools and Algorithms for the Construction and Analysis of Systems; Springer: Berlin/Heidelberg, Germany, 1998; pp. 102–117. [Google Scholar]
  33. Uygur, G.; Sattler, S.M. Parallel Composition-A practical solution. In Proceedings of the 1st International Electric Drives Production Conference (EDPC), Nuremberg, Germany, 28–29 September 2011; pp. 235–242. [Google Scholar]
  34. Rinaldi, M.; Capisani, L.; Ferrara, A.; Nunez, A.; Hajiahmadi, M.; De Schutter, B. Distributed identification of the cell transmission traffic model: A case study. In Proceedings of the American Control Conference (ACC), Montréal, QC, Canada, 27–29 June 2012; pp. 6545–6550. [Google Scholar]
  35. Vigos, G.; Papageorgiou, M.; Wang, Y. Real-time estimation of vehicle-count within signalized links. Transp. Res. Part C Emerg. Technol. 2008, 16, 18–35. [Google Scholar] [CrossRef]
  36. Wang, Y.; Papageorgiou, M. Real-time freeway traffic state estimation based on extended Kalman filter: A general approach. Transp. Res. Part B Methodol. 2005, 39, 141–167. [Google Scholar] [CrossRef]
  37. Anand, A.; Ramadurai, G.; Vanajakshi, L. Data fusion-based traffic density estimation and prediction. J. Intell. Transp. Syst. 2014, 18, 367–378. [Google Scholar] [CrossRef]
  38. Mihaylova, L.; Boel, R.; Hegyi, A. Freeway traffic estimation within particle filtering framework. Automatica 2007, 43, 290–300. [Google Scholar] [CrossRef]
  39. Kurzhanskiy, A.A.; Varaiya, P. Guaranteed prediction and estimation of the state of a road network. Transp. Res. Part C Emerg. Technol. 2012, 21, 163–180. [Google Scholar] [CrossRef]
  40. Kurzhanskiy, A.A. Set-valued estimation of freeway traffic density. In Proceedings of the IFAC Symposium on Transportation Systems, Redondo Beach, CA, USA, 2–4 September 2009; Volume 42, pp. 271–277. [Google Scholar]
  41. Kurzhanskiy, A.A.; Varaiya, P. Active traffic management on road networks: A macroscopic approach. Philos. Trans. R. Soc. Lond. A Math. Phys. Eng. Sci. 2010, 368, 4607–4626. [Google Scholar] [CrossRef] [PubMed]
  42. Alessandri, A.; Coletta, P. Design of Luenberger observers for a class of hybrid linear systems. In Proceedings of the International Workshop on Hybrid Systems: Computation and Control, Rome, Italy, 28–30 March 2001; Volume 2034, pp. 7–18. [Google Scholar]
  43. Alessandri, A.; Coletta, P. Switching observers for continuous-time and discrete-time linear systems. In Proceedings of the American Control Conference (ACC), Arlington, VA, USA, 25–27 June 2001; Volume 3, pp. 2516–2521. [Google Scholar]
  44. Bara, G.I.; Daafouz, J.; Kratz, F.; Iung, C. State estimation for a class of hybrid systems. In Proceedings of the 4th International Conference on Automation of Mixed Processes, Dortmund, Germany, 18–19 September 2000; pp. 313–316. [Google Scholar]
  45. Juloski, A.L.; Heemels, W.; Weiland, S. Observer design for a class of piece-wise affine systems. In Proceedings of the 41st IEEE Conference on Decision and Control, Las Vega, NV, USA, 10–13 December 2002; Volume 3, pp. 2606–2611. [Google Scholar]
  46. Juloski, A.L.; Heemels, W.; Boers, Y.; Verschure, F. Two approaches to state estimation for a class of piecewise affine systems. In Proceedings of the 42nd IEEE Conference on Decision and Control, Maui, HI, USA, 9–12 December 2003; Volume 1, pp. 143–148. [Google Scholar]
  47. Tanwani, A.; Shim, H.; Liberzon, D. Observability for switched linear systems: Characterization and observer design. IEEE Trans. Autom. Control 2013, 58, 891–904. [Google Scholar] [CrossRef]
  48. Pettersson, S. Observer design for switched systems using multiple quadratic Lyapunov functions. In Proceedings of the IEEE International Symposium on Intelligent Control, Limassol, Cyprus, 27–19 June 2005; pp. 262–267. [Google Scholar]
  49. Pettersson, S. Designing switched observers for switched systems using multiple Lyapunov functions and dwell-time switching. In Proceedings of the 2nd IFAC Conference on Analysis and Design of Hybrid Systems, Alghero, Italy, 7–9 June 2006; Volume 39, pp. 18–23. [Google Scholar]
  50. Gómez-GutiÉrrez, D.; Čelikovský, A.; Ramĺrez-Trevinños, A.; Castillo-Toledo, B. On the observer design problem for continuous-time switched linear systems with unknown switchings. J. Frankl. Inst. 2015, 352, 1595–1612. [Google Scholar] [CrossRef]
  51. Morbidi, F.; Ojeda, L.L.; Canudas-de-Wit, C.; Bellicot, I. A new robust approach for highway traffic density estimation. In Proceedings of the European Control Conference (ECC), Strasbourg, France, 24–27 June 2014; pp. 2575–2580. [Google Scholar]
  52. Alvarez-Icaza, L.; Muñoz, L.; Sun, X.; Horowitz, R. Adaptive observer for traffic density estimation. In Proceedings of the American Control Conference, Boston, MA, USA, 30 June–2 July 2004; pp. 2705–2710. [Google Scholar]
  53. Castillo, E.; Conejo, A.J.; Menéndez, J.M.; Jiménez, P. The observability problem in traffic network models. Comput. Aided Civ. Infrastruct. Eng. 2008, 23, 208–222. [Google Scholar] [CrossRef]
  54. Castillo, E.; Jiménez, P.; Menéndez, J.M.; Conejo, A.J. The observability problem in traffic models: Algebraic and topological methods. IEEE Trans. Intell. Transp. Syst. 2008, 9, 275–287. [Google Scholar] [CrossRef]
  55. Agarwal, S.; Kachroo, P.; Contreras, S. A dynamic network modeling-based approach for traffic observability problem. IEEE Trans. Intell. Transp. Syst. 2016, 17, 1168–1178. [Google Scholar]
  56. Agarwal, S.; Kachroo, P.; Contreras, S.; Sastry, S. Feedback-coordinated ramp control of consecutive on-ramps using distributed modeling and Godunov-based satisfiable allocation. IEEE Trans. Intell. Transp. Syst. 2005, 16, 2384–2392. [Google Scholar]
  57. Navarro-Lopez, E.M.; Carter, R. Hybrid automata: An insight into the discrete abstraction of discontinuous systems. Int. J. Syst. Sci. 2011, 42, 1883–1898. [Google Scholar] [CrossRef]
  58. Boriboonsomsin, K.; Barth, M. Impacts of freeway high-occupancy vehicle lane configuration on vehicle emissions. Transp. Res. Part D Transp. Environ. 2008, 13, 112–125. [Google Scholar] [CrossRef]
Figure 1. Digraph expression of hybrid automaton.
Figure 1. Digraph expression of hybrid automaton.
Sensors 17 00716 g001
Figure 2. A general DGHA model.
Figure 2. A general DGHA model.
Sensors 17 00716 g002
Figure 3. Triangular fundamental diagram.
Figure 3. Triangular fundamental diagram.
Sensors 17 00716 g003
Figure 4. Combination modes.
Figure 4. Combination modes.
Sensors 17 00716 g004
Figure 5. The schematic diagram of state observer.
Figure 5. The schematic diagram of state observer.
Sensors 17 00716 g005
Figure 6. Beijing third ring freeway (from Google Map). The points A and B are marked as the first on-ramp and the first off-ramp, respectively. The segment between A and B is labeled as cell 1. (Note: The Chinese word in this map is just the name of some buildings and will not affect the meaning of this image.)
Figure 6. Beijing third ring freeway (from Google Map). The points A and B are marked as the first on-ramp and the first off-ramp, respectively. The segment between A and B is labeled as cell 1. (Note: The Chinese word in this map is just the name of some buildings and will not affect the meaning of this image.)
Sensors 17 00716 g006
Figure 7. The simplified version of Beijing third ring freeway.
Figure 7. The simplified version of Beijing third ring freeway.
Sensors 17 00716 g007
Figure 8. Time-space diagram of traffic density.
Figure 8. Time-space diagram of traffic density.
Sensors 17 00716 g008
Figure 9. Error curve.
Figure 9. Error curve.
Sensors 17 00716 g009
Figure 10. Simulated and estimated densities of cells 14 and 20.
Figure 10. Simulated and estimated densities of cells 14 and 20.
Sensors 17 00716 g010
Table 1. Modes of i with relation to upstream l.
Table 1. Modes of i with relation to upstream l.
ModeCondition of ( ρ l , ρ i ) Flow q li
FDF 0 ρ l α l i , 0 ρ i ρ i 0 V l ρ l
FUF C i / V l ρ l ρ l 0 , 0 ρ i ρ i 0 C i
CDF/CUF ρ l 0 ρ l ρ l m , 0 ρ i ρ i 0 C l or C i
CDC ρ l 0 ρ l ρ l m , ρ i 0 ρ i ρ i m C l / W i C l
CUC ρ l 0 ρ l ρ l m , β l i ρ i ρ i m W i ( ρ i m ρ i )
FUC 0 ρ l ρ l 0 , ρ i 0 ρ i ρ i m , V l ρ l + W i ρ i W i ρ i m W i ( ρ i m ρ i )
FDC 0 ρ l ρ l 0 , ρ i 0 ρ l ρ i m , V l ρ l + W i ρ i W i ρ i m V l ρ l
Table 2. Cell length of simplified third ring freeway.
Table 2. Cell length of simplified third ring freeway.
NumberLengthNumberLengthNumberLength
1111 m8430 m15300 m
2535 m9500 m16350 m
3336 m10355 m17332 m
4138 m11338 m18210 m
5586 m12355 m19367 m
6140 m13458 m20458 m
7281 m14256 m
Table 3. Road segments parameters.
Table 3. Road segments parameters.
NumberVWC ρ 0 ρ m
2,3,46520280046185
7,8,96321285045180
12,13,146520276044182
17,18,196219265043180
1,206019245040170
5,65818235041165
10,115520220040150
15,165019210041155
Table 4. Mean square error of estimated densities.
Table 4. Mean square error of estimated densities.
Number of cell12345678910
MSE0.10950.11260.11320.12230.12050.11760.13010.10050.11080.1201
Number of cell11121314151617181920
MSE0.11320.10870.11360.10090.11430.12090.11200.10130.11470.1106

Share and Cite

MDPI and ACS Style

Chen, Y.; Guo, Y.; Wang, Y. Modeling and Density Estimation of an Urban Freeway Network Based on Dynamic Graph Hybrid Automata. Sensors 2017, 17, 716. https://doi.org/10.3390/s17040716

AMA Style

Chen Y, Guo Y, Wang Y. Modeling and Density Estimation of an Urban Freeway Network Based on Dynamic Graph Hybrid Automata. Sensors. 2017; 17(4):716. https://doi.org/10.3390/s17040716

Chicago/Turabian Style

Chen, Yangzhou, Yuqi Guo, and Ying Wang. 2017. "Modeling and Density Estimation of an Urban Freeway Network Based on Dynamic Graph Hybrid Automata" Sensors 17, no. 4: 716. https://doi.org/10.3390/s17040716

APA Style

Chen, Y., Guo, Y., & Wang, Y. (2017). Modeling and Density Estimation of an Urban Freeway Network Based on Dynamic Graph Hybrid Automata. Sensors, 17(4), 716. https://doi.org/10.3390/s17040716

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop