Sensors in Connected Vehicle Technology: How Sensors Play a Critical RoleView this Special Issue
Research Article | Open Access
Continuous and Discrete-Time Optimal Controls for an Isolated Signalized Intersection
A classical control problem for an isolated oversaturated intersection is revisited with a focus on the optimal control policy to minimize total delay. The difference and connection between existing continuous-time planning models and recently proposed discrete-time planning models are studied. A gradient descent algorithm is proposed to convert the optimal control plan of the continuous-time model to the plan of the discrete-time model in many cases. Analytic proof and numerical tests for the algorithm are also presented. The findings shed light on the links between two kinds of models.
Along with the fast increase of auto population, urban streets are becoming more crowded nowadays. To relieve congestions and reduce accidents, various traffic control methods have been proposed since the late 1950s .
As a typical traffic scenario, oversaturated intersections attracted consistent interest during the last six decades [2, 3]. The term “oversaturated” means the following: the vehicles that remained since the last cycle plus the vehicles that newly arrived exceed the capacity of the intersection. This leads to the carryover of vehicle queues (at least in one leg of the intersection) to the next cycle.
Discussions on the road networks that consist of many oversaturated intersections can be found in researches done by Chang and Sun , Di Febbraro and Sacco , Dotoli and Fanti , Ma [7, 8], and Sun et al. [9, 10]. In the research done by Varaiya  and Le et al. , the study of pressure-based signal control developed stability properties of a decentralized signal timing policy for networks with stochastic arrivals. But for a real-time signal timing optimization problem, the data that could be used is the arriving information of vehicles in the recent several signal cycles (data could be gotten by connecting vehicle technology, etc.). The optimization objects and scenarios are different between the model in this paper and pressure-based policies. This paper will focus on the isolated oversaturated intersection.
Usually, researchers aim to find an optimal signal timing plan that minimizes the total delay of vehicles passing this intersection. The total delay is often defined as the time integral of the sum of all queue lengths for all legs of the intersection over a given time horizon. However, the total delay is a nonlinear and nonconvex function of control variables (e.g., green phases), which makes it difficult to optimize. One promising approach is to apply heuristic algorithms to solve the formulated optimization problem. For example, the genetic algorithm was applied by Park et al.  to optimize the total delay. An alternative approach is to first approximate the nonconvex total delay with some convex functions and then solve the newly formulated optimization problem. The rest of this paper will focus on the second approach within a typical traffic scenario: an isolated intersection with only two movements.
There are mainly two kinds of convexified models for this scenario. The first kind of models originated from Gazis  who used continuous-time differential equations to describe the traffic dynamics. The cycle length, departing flow rates, and arriving flow rates in Gazis  were all assumed to be constant. Michalopoulos and Stephanopoulos [14, 15] extended the continuous-time model by including the maximum queue lengths constraints and time-variant arrival flow rates. Such formulations led to a classical control problem that can be solved via the Pontryagin Maximum Principle (PMP) . However, the obtained continuous-time signal timing plan should be discretized into the corresponding discrete-time signal timing plan that can be executed in practice.
The second kind of models uses discrete-time difference equations to describe the traffic dynamics [7, 8, 17–19]. The corresponding design problem can then be formulated as a linear programming (LP) problem [20–25]. One interesting question that naturally arises is how to depict the difference and connection between the discrete-time model and the continuous-time model.
Recently, Ioslovich et al.  studied the formulated LP problem by considering the corresponding continuous-time approximation model and gave an elegant approximate solution in continuous-time forms. However, it was not verified how this approximate solution differs from the accurate solution (the solution obtained by the continuous-time model). Whether a discretized version of this continuous-time approximate solution is still optimal to the LP problem also needs further discussions.
Zou et al.  have given a preliminary result on the relationship between the continuous-time model and the discrete-time model. Zou et al. applied a graphical method to adjust the continuous-time approximate solution to a discrete-time accurate solution. They assumed that both streams are dispatched simultaneously and the adjustment method does not influence the clearance cycle. It is shown that the approximate solution can be adjusted to the optimal solution by changing the green ratio in the switching cycle. However, the discrete-time LP problem in Ioslovich et al.  is more complex than the one considered in , since the two streams are allowed to be cleared at different times. Merely changing the green ratios at the switching cycle may not obtain the accurate solution.
In this paper, the LP problems proposed by Ioslovich et al.  are directly attacked using the strong duality theorem . It is first shown that the approximate solution and the accurate solution do not coincide in many situations. Then, the relationship between the approximate solution and the accurate solution will be discussed. The errors introduced by discretization are carefully studied. It is shown that, in many cases, the discretized approximate solution can be converted to the accurate solution within a few gradient descent adjustments. Finally, an algorithm is proposed to implement this conversion. These findings shed light on the connection between the continuous-time and discrete-time signal timing models.
To give a detailed analysis, the rest of this paper is arranged as follows. Section 2 introduces the LP problem proposed by Ioslovich et al.  and lists the nomenclature used in this paper. Section 3 proposes two counterexamples to show that the optimality of the discretized approximate solution may not hold. Section 4 discusses the difference and connection between the discrete-time model and the continuous-time model. Finally, Section 5 concludes the paper.
2. Problem Presentation
2.1. Nomenclature and Assumptions
For presentation simplicity, the nomenclature in this paper follows Ioslovich et al.  as shown in Nomenclature List. Similar to Ioslovich et al. , consider an isolated intersection with two one-way streams and governed by signal lights, as shown in Figure 1. Note that the two-stream isolated intersection is chosen as an initial building block to understand the connection between the continuous-time model and the discrete-time model. This is a widely applied treatment even in recent works (e.g., ). Nevertheless, this assumption can be relaxed to consider general intersections. Both the continuous model and the discrete model are able to handle cases with more generalized intersections. Although it is hard to directly compare the two types of models in generalized intersections, we expect that the findings in this paper can shed light on more generalized intersections.
Likewise, the following assumptions are imposed as in Ioslovich et al. .
Assumption 1. (a) Each signal cycle starts with a green light for .
(b) There is no lost time for each cycle. Suppose that the green ratio for stream in cycle is denoted as ; the green ratio for stream is then equal to .
(c) The arrival flow rates and the saturation departure flow rates of both streams are constant over time. Denote as the arrival flow rates for and as the departure flow rates.
(d) Without loss of generality, is assumed to be the main stream (i.e., ).
(e) For some reasons (e.g., physical constraints), both streams are constrained by minimum green ratios. Define and as the minimum green ratios for and , respectively; then, the constraint on is represented as .
(f) There exists a signal timing strategy such that queues in both streams can be cleared finally but may not be in the same cycle. Define and . This assumption is written as , , and .
(g) Cycle lengths are known to be constant .
In Assumption 1, (a)–(f) are inherited from Ioslovich et al.  directly and (g) is added to simplify the discussion. In fact, if cycle lengths are not fixed, a similar conclusion can still be drawn but the detailed timing plan is much more complicated.
2.2. The Discrete-Time LP Model and Its PMP Solution
Similar to Ioslovich et al. , the discrete-time model is formulated as the following problem to minimize the total delay of vehicles:for .
In (1)–(5), the objective function is a linear approximation for total delay. The decision variables are , . For presentation convenience, it is assumed that the number of total cycles is large enough such that both streams are cleared before cycle . Equations (2) and (3) represent the evolution of both queues in time. Equation (4) gives the upper and lower bound of the green ratio in each cycle. Equation (5) is the initial queue lengths. It is easy to verify that this LP problem is equivalent to the Relaxed Discrete-Event Max-Plus Problem proposed by Haddad et al. .
Particularly, the following four cases were discussed with respect to the order of , and : (I) ; (II) ; (III) ; (IV) .
Table 1 lists the solution. The optimal continuous-time solution is denoted as. is the switching time where the optimal control switches. and are represented as
For example, Figure 2 illustrates the evolution of queue lengths under the optimal solution for Case I(a). It is observed that the optimal solution is a two-stage strategy, namely, bang-bang control. The green ratio remains at first, causing the length of stream to decrease sharply and stream to increase. After the switching time , the green ratio changes into and both streams start to decrease. The optimal switch-over point where the solution switches from maximum to minimum green split is given in Table 1.
The discretized solution of is obtained as follows. Suppose the cycle length is ; the discretization is done as
For presentation simplicity, the solution , , is called the discretized approximate solution in the rest of this paper.
2.3. The Steady State and the Steady-State Solution
According to Haddad et al. , the steady state is defined as the state where the queue lengths are the same in the beginning and end of each cycle. This means, given the green ratio for a certain stream, the queue lengths of this stream keep the same in a certain period. Denote the queue lengths in the steady state as , and the green ratio in the steady state as . So, considering a single-cycle version of discrete-time model (see (1)–(5)), the steady state can be formulated into the following model:
The optimal green ratio is called the steady-state solution. Haddad et al.  gave a neat analytical form of as follows:
According to Haddad et al. , the two streams that follow the discrete-time model (see (1)–(5)) will eventually enter the steady state under Assumptions (a)–(g). Therefore, in the accurate solution, the green ratios in the steady state are the steady-state solution .
3. The Nonoptimality of the Discretized Approximate Solution
In this section, it is shown that the discretized approximate solution does not always coincide with the accurate solution. The term “coincide” here means that the discretized approximate solution is not always optimal for the discrete-time model. In particular, two counterexamples are constructed. The first counterexample is used to deal with Case I and Case II where the minimum green ratio plays an important role in the approximate solution. In contrast, the second counterexample is used to handle Case III and Case IV where the minimum green ratio does not influence the approximate solution. The two counterexamples are as follows.
Counterexample 1 (for Case I and Case II). The parameters are set as = , = , = , = , = , = , = , = , = , and = . Then, the accurate solution and the discretized approximate solution are shown in Figure 3(a).
Counterexample 2 (for Case III and Case IV). The parameters are set as = , = , = , = , = , = , = , = , = , and = . Then, the accurate solution and the discretized approximate solution are shown in Figure 3(b).
As is shown in Figure 3, the accurate solution deviates from the discretized approximate solution in two ways. In Counterexample 1, the accurate solution roughly remains as a bang-bang form before the system enters steady state. The difference concentrates around the switching cycle and steady state. However, in Counterexample 2, the accurate solution exhibits oscillation and differs from the discretized approximation solution in almost every cycle.
Such deviations are caused by discretization and can be classified into two categories.
The first category of deviations is caused by the discretization of the model constraints. For the continuous-time model, the queue length for stream is lower-bounded by the constraint , while for the discrete-time model, the corresponding constraint is . When cycle length is large, the discretization of model constraints may introduce huge errors. In Counterexample 1, the switching time locates in the middle of a cycle. As shown in Figure 4, the discretization leads to waste of green ratio around the switching cycle due to the restriction of the minimum green ratio. When is small enough, the solution of the continuous-time model and the solution of the discrete-time model are coherent. A proper should be selected to ensure consistency.
The second category deviation is caused by the discretization of the approximate solution. In Counterexample 2, the clearance time of stream locates in the middle of a cycle. So, the discretization causes oscillation in the accurate solution.
4. The Relationship between the Discretized Approximate Solution and the Accurate Solution
In this section, the connection between the discrete approximate solution and the accurate solution is further discussed.
For Case I and Case II, the errors are relatively concentrated. The discretized approximate solution does not deviate largely from the accurate solution in many cases. In fact, the discretized approximate solution can be modified to the accurate solution by some minor adjustments; see Algorithm 1. The details of adjustments will be discussed in Sections 4.1–4.3.
For Case III and Case IV, however, the dispersed errors make it hard to adjust the solution in a few steps. However, for some carefully chosen cycle length, the accurate solution coincides with the discretized approximate solution.
For example, in Case III, if the cycle length is chosen as , where is any positive integer, the accurate solution can be written as
It is easy to show that this accurate solution coincides with the corresponding discretized approximate solution.
In the rest of this section, the discussion will focus on the adjustments for Case I, and the adjustments for Case II are similar.
For presentation simplicity, define several denotations before starting the discussion. Define and , respectively, as the clearance cycle indices in which streams and are cleared exactly:
Apparently, in Case I(a), it holds that . So, both streams enter steady state after cycle . Note that there is also a switch in the discretized approximate solution. Define as the switching cycle index where switch occurs; that is,
The adjustments contain two parts. Section 4.1 deals with the adjustment in the steady state. Section 4.2 includes the adjustment around the switching cycle. Both parts are summarized as Algorithm 1. The optimality of Algorithm 1 will be proven in the Appendix. Finally, Counterexample 1 in Section 3 is revisited to better illustrate the adjustments.
4.1. The Steady-State Adjustment
The adjustment for the steady state is shown in lines (34)–(38) in Algorithm 1. Since the steady-state solution is optimal for the discrete-time model after the two streams enter steady state, can be set as
However, in the clearance cycle , the optimal green ratio is not due to discretization. So, can be set as follows:
4.2. The Adjustment around the Switch
The adjustment around the switching cycle is shown in lines ()–(33) in Algorithm 1. There are two cases of adjustment: increasing the green ratios after the switch and decreasing the green ratios before the switch. For presentation simplicity, this subsection only handles the case of increasing the green ratios after the switch. The adjustment to the other case is similar except for some details; see Algorithm 1, lines ()–().
This adjustment can be viewed as a gradient descent approach to gradually reduce the total delay. It operates in three steps: () determine the direction of adjustment; () update the green ratio; () check feasibility of the new green ratio. It will be proven in the Appendix that such adjustment operations can be finished in a few steps.
Due to constraint equations (3) and (4), the discrete-time model is a nonsmooth model. Thus, a subdifferential will be used instead of a gradient. Without loss of generality, assume that both streams are cleared after the switching cycle (i.e., ). The case of can be transformed easily to the case of (Algorithm 1, lines ()–()).
In the first step, the direction of adjustment is determined by calculating the subdifferential. Denote the subdifferential of the total delay with respect to as . Define the function asso as to represent the subdifferential in a more convenient way.
Clearly, the case of increasing the green ratios after the switch corresponds to the condition of . There are three situations based on the queue evolution in clearance cycles and .
() Stream is cleared exactly in cycle ; that is, . The subdifferential in this situation is calculated aswhere “” means the closed interval between real numbers and .
The physical meaning of the subdifferential is shown as follows. The total delay changes if a slight disturbance is imposed on . Denote as the increment in and as the corresponding change in total delay. If , then can be calculated as
Similarly, if , the corresponding changes in total delay are written as .
Therefore, if , the total delay can be reduced by increasing . If and , neither increasing nor decreasing will render a better total delay and thus the optimal is achieved.
() Stream is cleared exactly in cycle ; that is, . The subdifferential can be calculated similarly as .
Then, if , the total delay can be reduced by increasing . If and , the optimal choice of is obtained.
() Both streams are not cleared exactly in their clearance cycles. Then, the subdifferential is . If , the total delay can be reduced by increasing .
In the second step, the green ratio is updated by reducing and increasing . The green ratio is calculated such that either stream or stream is cleared exactly in its clearance cycle. The detailed calculation is shown in Algorithm 1 (lines ()–()).
The third step is the feasibility check of whether the new green ratio satisfies the constraint equation (4). In the case of increasing , the only possible violation is . Therefore, if the feasibility check fails, is set as and a new bang-bang solution is obtained. Then, the algorithm goes back to the first step and determines the direction of adjustment.
4.3. The Numerical Examples
In this subsection, Counterexample 1 in Section 3 is revisited to illustrate the steps of Algorithm 1. Recall that the parameters are = 0.55, = 0.3, = 0.15, = 0.10, = 0.4, = 0.8, = 60, = 20, , and . Algorithm 1 will be applied to adjust the discretized approximate solution. The detailed operation is shown in Table 2.
The adjustment of the discretized approximate solution is also shown in Figure 5.
This paper studies the continuous-time and discrete-time signal timing models for the classical isolated signalized intersection with only two one-way vehicle flows. Against intuitions, the nonequivalence between the solutions of the two models is first shown by two counterexamples. Then, the differences between the solutions of the two models are explained. Finally, an algorithm is proposed to transform the discretized continuous-time solution to the discrete-time solution for many cases. Via this gradient descent based algorithm, a more concrete link is set up between the continuous-time and discrete-time signal timing models.
Proof of the Optimality of Algorithm 1
The optimality is stated in the following theorem.
Theorem A.1. The solution derived by Algorithm 1 is optimal for Case I(a).
Lemma A.2 (strong duality theorem for linear programming). A primal feasible solution and a dual feasible solution to a linear programming problem are optimal if and only if the corresponding primal and dual objectives are the same.
Proof of Theorem A.1. For presentation simplicity, only the case of increasing the green ratios after the switch will be considered in this proof. The other case can be similarly proven. Since the adjustment for the switching cycle (lines ()–()) is a gradient descent approach and the discrete LP model is apparently feasible, the algorithm will eventually terminate.
The dual problem to the discrete LP model (see (1)–(5)) is written asThere are three situations when the iteration of lines ()–() terminates.(1)The iteration terminates at line (). Then, is given such that stream is cleared exactly. The corresponding , , and satisfy and .(2)The iteration terminates at line () and is set as . The corresponding , , and satisfy .(3)The iteration terminates at line () and is given such that stream is cleared exactly. The corresponding , , and satisfy and .The proofs of the three situations is similar. Only the first situation is handled here. Moreover, assume . The proof of is the opposite.
Denote and , respectively, as and . Then the dual solution is rewritten asSubstituting (A.5) and (A.6) into (A.4), can be represented asWith simple calculation, the sign of is determined asFinally, it is obtained by comparing the primal and dual objectives thatThis proves that the primal and dual objective values are equal. By Lemma A.2, both primal and dual solutions are optimal. Therefore, Algorithm 1 gives the optimal solution for Case I(a).
Nomenclature ListKnown Constants
|:||The total number of cycles in the observation period|
|:||The cycle length of the isolated intersection|
|:||The two streams,|
|:||The arrival rate of stream ,|
|:||The saturation departure flow rate of stream ,|
|:||The initial queue length of stream ,|
|:||The green ratio that keeps stream unchanged (|
|:||The green ratio that keeps stream unchanged (|
|:||The minimum and maximum bounds of green ratio|
|:||The steady-state solution of the discrete-time model.|
|:||The switching time of signal timing plan in the continuous-time optimal solution|
|:||The final time in the continuous-time model|
|:||The objective of the continuous-time model, representing the approximate total delay|
|:||The green time ratio of stream at time , the decision variable for continuous-time model|
|:||The approximate solution, that is, the optimal solution to the continuous-time model.|
|:||The green time ratio of stream in cycle , the decision variable for the discrete-time model|
|:||The queue lengths of stream in the beginning of cycle|
|:||The index of the cycle in which the stream is cleared|
|:||The index of the cycle in which the stream is cleared|
|:||The index of the cycle in which switches from to|
|:||The objective of the discrete-time model, representing approximate total delay|
|:||The accurate solution, that is, the optimal solution to the discrete-time model|
|:||The discretized solution of .|
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Support by Beijing Natural Science Foundation (4164083), National Natural Science Foundation of China (61603005, 71303269), Fundamental Research Funds for the Central Universities, and Central University of Finance and Economics Major Research Task of Fostering Project (14ZZD006) is gratefully acknowledged.
- F. Webster, Traffic signal settings. Road Research Technical Paper 39. Great Britain Road Research Laboratory, London, UK, 1958.
- D. Green, “Control of oversaturated intersections,” Operational Research Quarterly, vol. 18, no. 2, pp. 161–173, 1967.
- D. C. Gazis, “Optimum control of a system of oversaturated intersections,” Operations Research, vol. 12, no. 6, pp. 815–831, 1964.
- T.-H. Chang and G.-Y. Sun, “Modeling and optimization of an oversaturated signalized network,” Transportation Research Part B: Methodological, vol. 38, no. 8, pp. 687–707, 2004.
- A. Di Febbraro and N. Sacco, “On modelling urban transportation networks via hybrid Petri nets,” Control Engineering Practice, vol. 12, no. 10, pp. 1225–1239, 2004.
- M. Dotoli and M. P. Fanti, “An urban traffic network model via coloured timed Petri nets,” Control Engineering Practice, vol. 14, no. 10, pp. 1213–1229, 2006.
- X. Ma, H. Yu, Y. Wang, Y. Wang, and J. Gomez-Gardenes, “Large-scale transportation network congestion evolution prediction using deep learning theory,” PLoS ONE, vol. 10, no. 3, Article ID e0119044, 2015.
- X. Ma, Z. Tao, Y. Wang, H. Yu, and Y. Wang, “Long short-term memory neural network for traffic speed prediction using remote microwave sensor data,” Transportation Research Part C: Emerging Technologies, vol. 54, pp. 187–197, 2015.
- W. Sun, Y. Wang, G. Yu, and H. X. Liu, “Quasi-optimal feedback control for a system of oversaturated intersections,” Transportation Research Part C: Emerging Technologies, vol. 57, pp. 224–240, 2015.
- W. Sun, Y. Wang, G. Yu, and H. X. Liu, “Quasi-optimal feedback control for an isolated intersection under oversaturation,” Transportation Research Part C: Emerging Technologies, vol. 67, pp. 109–130, 2016.
- P. Varaiya, “Max pressure control of a network of signalized intersections,” Transportation Research Part C: Emerging Technologies, vol. 36, pp. 177–195, 2013.
- T. Le, P. Kovács, N. Walton, H. L. Vu, L. L. H. Andrew, and S. S. P. Hoogendoorn, “Decentralized signal control for urban road networks,” Transportation Research Part C: Emerging Technologies, vol. 58, pp. 431–450, 2015.
- B. Park, C. J. Messer, and T. Urbanik II, “Traffic signal optimization program for oversaturated conditions: genetic algorithm approach,” Transportation Research Record, no. 1683, pp. 133–142, 1999.
- P. G. Michalopoulos and G. Stephanopoulos, “Oversaturated signal systems with queue length constraints—I: single intersection,” Transportation Research, vol. 11, no. 6, pp. 413–421, 1977.
- P. G. Michalopoulos and G. Stephanopoulos, “Oversaturated signal systems with queue length constraints-II. Systems of intersections,” Transportation Research, vol. 11, no. 6, pp. 423–428, 1977.
- L. Pontryagin, V. Boltyanskii, R. Gamkrelidze, and E. Mishchenko, The mathematical theory of Optimal processes, vol. 43, Wiley-Interscience, 1962.
- R. Allsop, “SIGSET: A computer program for calculating traffic signal settings,” Transportation Engineering Control, vol. 13, no. 2, pp. 58–60, 1971.
- R. E. Allsop, “Estimating the traffic capacity of a signalized road junction,” Transportation Research, vol. 6, no. 3, pp. 245–255, 1972.
- R. E. Allsop, “SIGCAP: A computer program for assessing the traffic capacity of signal-controlled road junctions,” Transportation Engineering Control, vol. 17, no. 819, pp. 338–341, 1976.
- J. Haddad, B. De Schutter, D. Mahalel, and P.-O. Gutman, “Steady-state and N-stages control for isolated controlled intersections,” in Proceedings of the 2009 American Control Conference, ACC 2009, pp. 2843–2848, June 2009.
- J. Haddad, B. De Schutter, D. Mahalel, I. Ioslovich, and P.-O. Gutman, “Optimal steady-state traffic control for isolated intersections,” in Proceedings of the 6th IFAC symposium on robust control design, 2009, pp. 96–101.
- J. Haddad, B. De Schutter, D. Mahalel, I. Ioslovich, and P.-O. Gutman, “Optimal steady-state control for isolated traffic intersections,” Institute of Electrical and Electronics Engineers. Transactions on Automatic Control, vol. 55, no. 11, pp. 2612–2617, 2010.
- L. Li, K. Yang, Z. Li, and Z. Zhang, “The optimality condition of the multiple-cycle smoothed curve signal timing model,” Transportation Research Part C: Emerging Technologies, vol. 27, pp. 46–57, 2013.
- K. Han, V. V. Gayah, B. Piccoli, T. L. Friesz, and T. Yao, “On the continuum approximation of the on-and-off signal control on dynamic traffic networks,” Transportation Research Part B: Methodological, vol. 61, pp. 73–97, 2014.
- F. Zhu and S. V. Ukkusuri, “A linear programming formulation for autonomous intersection control within a dynamic traffic assignment and connected vehicle environment,” Transportation Research Part C: Emerging Technologies, vol. 55, pp. 363–378, 2015.
- I. Ioslovich, J. Haddad, P.-O. Gutman, and D. Mahalel, “Optimal traffic control synthesis for an isolated intersection,” Control Engineering Practice, vol. 19, no. 8, pp. 900–911, 2011.
- B. Zou, J. Hu, and Y. Zhang, “The optimal discretized timing plan for individual oversaturated intersections,” in Proceedings of the 2012 15th International IEEE Conference on Intelligent Transportation Systems, ITSC 2012, pp. 1656–1660, September 2012.
- D. G. Luenberger and Y. Ye, Linear and Nonlinear Programming, vol. 116 of International Series in Operations Research & Management Science, Springer, New York, NY, USA, 3rd edition, 2008.
Copyright © 2017 Jiyuan Tan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.