Research Article | Open Access

# Distributed Model Predictive Control over Multiple Groups of Vehicles in Highway Intelligent Space for Large Scale System

**Academic Editor:**Wuhong Wang

#### Abstract

The paper presents the three time warning distances for solving the large scale system of multiple groups of vehicles safety driving characteristics towards highway tunnel environment based on distributed model prediction control approach. Generally speaking, the system includes two parts. First, multiple vehicles are divided into multiple groups. Meanwhile, the distributed model predictive control approach is proposed to calculate the information framework of each group. Each group of optimization performance considers the local optimization and the neighboring subgroup of optimization characteristics, which could ensure the global optimization performance. Second, the three time warning distances are studied based on the basic principles used for highway intelligent space (HIS) and the information framework concept is proposed according to the multiple groups of vehicles. The math model is built to avoid the chain avoidance of vehicles. The results demonstrate that the proposed highway intelligent space method could effectively ensure driving safety of multiple groups of vehicles under the environment of fog, rain, or snow.

#### 1. Introduction

A new research concept about highway intelligent space (HIS) is proposed in the literature [1]. The main objective of the HIS system is to create an intelligent driving space where sensors are arranged in some necessary places such as highway tunnel, road sections that are prone to accidents or poor visibility under adverse weather conditions. Meanwhile, advanced vehicle-to-vehicle (V2V) and vehicle-to-server (V2S) are adopted to inform drivers of road section situation ahead to avoid collision under the environment of space communication mode. The HIS system is introduced to provide vehicles with useful driving state to ensure that vehicles could have better access to road information ahead, especially multiple groups of vehicles of driving environment. Hence, studying vehicle driving state is important for the HIS system to build an intelligent space according to the above mentioned situations. This paper mainly focuses on multivehicle of driving state when passing through the highway tunnel.

The highway tunnel has its own complex driving features, for example, the channel tunnel fire resulting in multiple fatalities and injuries [2], several traffic accidents appearing in some tunnels [3], and vehicle crash caused by visual dysfunction of drivers [4]. When accidents happened in some highway tunnel, the effective vehicle state should be transmitted to vehicles towards highway tunnel in real time by the communication between vehicles and HIS system, especially multivehicle driving towards highway tunnel. Multivehicle has different driving modes such as single vehicle form, groups of vehicle forms, and vehicle organizing forms, among which groups of vehicle driving forms are the common. While HIS server could store large amounts of data information for vehicles, multivehicle can be decomposed to different groups of vehicles driving towards highway tunnels and different groups of vehicles need different driving states according to the specific driving environment. Each group of vehicles driving states is defined as information framework. Information framework mainly shows that the effective driving state obtained from HIS server can ensure different groups of vehicles of safety characteristics. Therefore, studying multiple groups of vehicles driving characteristics as well as accurate evaluations of the risk rear-end collisions is the key technologies under the environment of highway tunnel.

For the large scale of multi-intelligent vehicle systems, there are three MPC schemes for designing such large scale systems, namely, distributed MPC, centralized MPC, and decentralized MPC. The distributed MPC strategy is arguably the most promising one because it beats the centralized MPC one in terms of computational load and outperforms the decentralized MPC one in control performance [5]. In a DMPC design, different model predictive controllers communicate through a communication network to cooperate their actions in order to achieve optimal performance [6], for example, collision avoidance constraints [7]. The distributed framework of MPC is also gradually developing for the control of large scale system. There are two methods of distributed MPCs that appeared in the literature for the large scale system. One method is that each local controller exchanges estimation states with its neighbors and therefore improves the performance of closed-loop subsystem. However the performance of other subsystems is not considered in this optimization. The second method could achieve a good performance close to the centralized MPC. However, this strategy requires much more communication resources and the structure of controller is relatively complex [8]. In this paper, the neighbor optimization method is used for the large scale system in which each subsystem interacts in sequence by state.

Based on the above state, a new vehicle safety distance among groups based on the HIS system is proposed to develop a warning system. Three time warning distances will be studied considering the complex tunnel space. In nature, human beings suffer from perception limitations with a typical reactive time of 0.75 seconds to 1.5 seconds on emergency events. Therefore, it is highly important that vehicles should work in the three warning distances to guarantee its immediate and effective stopping to avoid collisions. The first warning distance is calculated by the time, the driver reaction time, and the message propagation delay. As to the second warning distance, the compensated safety distance is proposed as the buffer distance to reduce its acceleration to the range of the vehicle near the tunnel. The third warning distance is the tracking safety distance. The proposed HIS system is adapted for each group and can effectively achieve vehicle’s safety. Figure 1 shows the work principle among the groups. Distributed model prediction control approach is used to solve groups of vehicles collision avoidance. Meanwhile, the method based on neighbor optimization is used to solve the large scale system. Therefore, by predicting the states based on the prediction model over the finite horizon, the third warning distance which is collision-free at discrete time steps is planned. In other words, since the collision avoidance is considered only at the discrete prediction time steps, a collision may occur in the intervals between the prediction time steps. In particular, since it is difficult to keep the sampling steps small enough to ensure the convergence of the trajectory, the collisions between prediction time steps become quite significant [9]. Therefore, as to the multiple groups of vehicles collision avoidance problem, adopting the above mentioned could solve vehicle environment perception.

The organization of the paper is as follows. System description and objective are presented in Section 2. In Section 3, distributed model predictive control approach is studied. In Section 4, the stability of vehicles in each group is studied. In Section 5, simulation and analysis are presented. Finally, conclusions are given.

#### 2. System Description and Objective

In this section, we define the system dynamics and pose an integrated cost function for every group to ensure system stabilization and vehicles safety driving. Meanwhile, the three time warning distances are studied addressing the system dynamics.

##### 2.1. System Description

Figure 2 shows the basic work principle of overall system. The set is used to denote the multivehicle driving towards highway tunnel. Autonomous vehicles organize different groups according to the specific driving environment or driving state. Therefore, the symbol is used to denote group-1, group-2,, group-. For each group-, the leading vehicle should track the three time warning distances. The last vehicle acts on the intervehicle dynamics model between it and the next group. For each following vehicle in each group, the reference trajectory should track the leading vehicles. Therefore, distributed model predictive control approaches used for the leading vehicles of each group can achieve the vehicles safety driving characteristics. This study is divided into several groups. This class of serially connected subgroup is composed of many similar subgroups placed after one another, in such a way that each subgroup is connected with dynamic state between its neighbors. Some algorithms and assumptions are made as follows.

*Assumption 1. *The vehicles in each group are self-organization form; the following vehicles in each group can favorably track the leading vehicle. The time delay communication of the vehicle driving state transmission between vehicles is not considered.

*Assumption 2. *Based on highway tunnel characteristics, the overtaking situation is not considered in the research range and the velocity change of vehicles is small.

*Algorithm 3. *Distributed model predictive controller for the leading vehicle : initialization: the driving state of the vehicle is sent to the first group by HIS server, starting at time . In addition, other groups will adjust to the driving states by the communication between the HIS system and the groups.

##### 2.2. Analysis of the Three Time Warning Distances

The safety of the driving behavior for each group is typically related to vehicle driving environment. Hence, regarding each group of safety aspect, the information framework includes intervehicle distance, relative velocity, and acceleration. The basic principle of the three time warning distances can be seen from Figure 2. When an accident has happened in a highway tunnel, the vehicle near the tunnel has obtained the message and lowers its velocity subconsciously. The leading vehicle received the road information ahead by the HIS server and passed through the three time warning distances. The latter groups will adjust to the driving states based on the vehicle and the former groups. The overall groups of distance could be shown as

The first time warning distance for each leading vehicle of safety distance is calculated before arriving at its second warning and mainly includes two parts of time, the reaction time of the driver and message propagation delay . Message propagation delay is composed of two parts, and . is the time when the message about an accident in tunnel is sent to the HIS server and is the time when a message is sent to related leading vehicles of the group. With this notation, we can define the first warning distance:

Furthermore, we can solve the first time warning distance for each group denoted by

The second time warning distance is calculated before arriving at its third time warning distance. The goal for the second time warning distance is to lower the leading vehicle of velocity; after arriving at the third time warning distance, the leading vehicle of driving states could satisfy the tracking mode.

*Definition 4. *
Compensated safety distance in our work is defined as the buffer distance among groups to avoid collisions.

Based on the above assumptions and definition, the model for compensated safety distance for each group is denoted as follows: where is the compensated safety distance and is the desired deceleration for the leading vehicle.

*Definition 5. *The vehicle is in the safety driving when the desired acceleration is defined in the range . is a special state that denotes the leading vehicles entering the tracking mode after arriving at the third warning distance.

The second time warning distance is a buffer range, so we introduce the virtual vehicle concept to calculate the distance. Additional advantage of the virtual vehicle scheme is that the motion of the vehicle can be smoothly controlled when a new leading vehicle cuts in or the current vehicle cuts out [10].

*Definition 6. *The virtual vehicle distance is defined as a constant safety distance. For every group , the compensated safety distance should be set in the range as follows:

The constant safety distance of virtual vehicle distance can be shown:

The second time warning distance is built to reduce vehicle’s acceleration during some adverse weather or road sections that are prone to accidents.

The third warning time is the tracking mode before arriving at its permitted minimum distance. The vehicle is uncontrollable, so the leading vehicle in each group should be controlled to ensure a suitable relative distance. As can be seen from Figure 2, the intervehicle dynamics model is designed among groups as follows:

The longitudinal dynamics of the leading vehicle are nonlinear. According to the vehicle dynamics in [11], the longitudinal dynamics is transferred as follows:

If the parameters in (7) are exactly known, the following feedback linearizing control law could be adopted:
where is intervehicle distance, is the error among groups, and is the error between actual intervehicle distance and desired vehicle safety distance. , are vehicle3’s speed and vehicle1’s speed. is the input signal that makes the closed-loop system satisfy certain performance criteria. In controller (8), we achieve the following objectives: the feedback linearization results in a linear system as discussed above: however uncertainties in parameters can potentially make the linearization process inexact; the study of such a case would be an interesting topic to be considered in future research; the simplification of the system model by excluding some characteristic parameters (e.g., the mechanical drag, mass and air resistance) from the vehicle dynamics. Manipulating (6) through (8), the equation becomes
where is the engine time constant and its value is 0.25, a single parameter that describes the dynamics of the propulsion system and internal disturbances. can be viewed as the throttle/brake input causing acceleration/deceleration in the controlled vehicle. The system thus takes the form which can be described by the following standard equations:

Writing the above equation as standard state-space equations, we have
where state variables are and is the acceleration of the last vehicle .

By (12), the overall system of vehicle dynamics can be expressed as

Simplifying the above equation,
where , , , and = .

*Criterion 1. *The dynamics model (12) is stable, controllable, and observability.

The stability system is analyzed by eigenvalue criteria. The controllability is verified by using rank and the observability is studied by using rank .

As we have defined it, the three time warning distances are designed to imply that, combined with HIS space, vehicles safety driving could be guaranteed under the environment of some uncertain weather conditions, such as fog, rain, or snow.

#### 3. Distributed Model Predictive Control Approach

##### 3.1. Induction Process

In this section, we introduce notation and define the optimal control problem and the distributed model predictive control approach (DMPC) for the leading vehicles in each group. Combined with the above schematics derivation, the DMPC only is developed in the leading vehicle ; the following vehicles in each group can maintain the safety driving distance in the form of the platoon. The DMPC algorithm captures an important class of practical problems, including, for example, maneuvering a group of vehicles from one point to another while maintaining relative formation and/or avoiding collisions [12]. Figure 3 shows the principle of the optimal control application for each leading vehicle.

As stated above, each group exchanges information with the HIS server system. Each group sequentially achieves performance development in distributed MPC algorithms under the environment of the highway tunnel.

Because the latter groups irregularly drive on the highway tunnel, the interconnections between different subsystems are assumed to be weak and are considered as disturbances which can send the relative driving information to the latter group via HIS server system. In addition, the driving states of the last vehicle in each group are uncontrollable but can be measured by sensors, so the acceleration of the last vehicle can be regarded as the measured disturbance signal when calculating the latter groups. The overall of longitudinal dynamics system that is composed of interconnected subsystems can be descried by (14). The global optimization problem can be decomposed into a number of local optimization subproblems and the whole control performance can be efficiently improved [13]. The cooperation between subsystems is achieved by exchanging information between each subsystem and its neighbors in a distributed structure via network between HIS server and vehicles.

By (12), according to the vehicle driving characteristics in highway intelligent space, subsystem interacts with , and the output state acceleration of subsystem is affected by subgroup . In this case, is called input neighboring subgroup of . is called the output neighboring subsystem. In (12), the state acceleration of the is regarded as disturbance.

In every group of distributed optimal control strategy, we assume that the same constant prediction horizon and constant update period are used. In practice, the update period is the sample interval. The common update times are denoted by , where and . The leading vehicle in each group sequentially solves an optimal control problem at the update period and applies the optimal control trajectory until its next update time. We have that and are the actual error state and control input, respectively. For each leading vehicle at any time , over any prediction interval , , associated with updated time , we denote two trajectories: : the optimal control trajectory, : the predicted state trajectory,where .

The predicted value is transmitted to all other followers as soon as the optimal control problem at is solved to take account of collision avoidance. The overall system running process is as follows.

*Step 1. *At the time = , the leading vehicle building the model predictive controller is defined.

Input State.;

Disturbance Variance . . Because of the uncontrollable vehicle , the can be regarded as a disturbance variance.

*Step 2. *The DMPC used for the leading vehicle obeys the following implementation strategy.(1)At the time , all the controllers receive the full state measurement from the sensors.(2)Each controller evaluates its own future input trajectory based on and . Based on the input trajectories, each controller calculates the current decided set of inputs trajectories.(3)Each controller updates its entire future input trajectory and sends the future state to the following vehicles in each group.(4)When a new measurement is received, go to Step 1 ().

At each iteration, for , each controller solves the following optimization problem:
with
where and are the weighting matrices that tune the relative importance of the output vector’s elements as well as the magnitude of the control effort. formulation is used to solve the model predictive control problem.

##### 3.2. Control Objectives and Constraints

Typically, the primary control objective of the HIS server system is to keep vehicles at a safety distance. So there exit some constraints on the information framework for each group and the information unit for each vehicle.

*Step 1. *The first parameter is the intervehicle distance. The intervehicle distance includes two parts. One part is the intervehicle distance in each group, and the other is the intervehicle distance between the last vehicle and the leading vehicle . About the first part of the intervehicle distance, the so-called safety distance can be defined as (6). As to constraints of relative distance error, when the vehicle runs uniformly, larger or smaller intervehicle distance may occur in real time. The inequality for subobjective is as follows:
where m is the lower boundary, which is again obtained from the driver experimental data, and m is the higher boundary in the literature [14].

The intervehicle distance between the last vehicle and the leading vehicle can be studied on the three time warning distances. The second time warning distance constraint can be shown in Definition 5.

The third time warning distance is the tracking mode and the constraint is the same as the first part of intervehicle distance.

*Step 2. *The second parameter is the intervehicle velocity. The intervehicle velocity also includes two parts. One is the intervehicle velocity in each group, and the other is the intervehicle distance between the last vehicle and the leading vehicle . About the first part of the intervehicle velocity, we define the initial vehicle velocity as . The relative velocity in each group should be minimized in the range . The second intervehicle velocity is researched based on the third time warning distance, because, during the first and the second time warning distances, the vehicle is in acceleration or deceleration without any rule. The constraint is in the range .

*Step 3. *The third parameter is the acceleration of the leading vehicle in each group. First the vehicle of acceleration constraint is defined in the range , . The absolute value of being bigger than that of can accommodate larger braking degree to prevent rear-end collisions [15].

*Step 4. *The fourth parameter is the weight and the control input. The weight which is the weight of the error between the desired and the actual distance has to be considered. The larger is, the smaller the time reaches a steady-state situation. Although the focus is on safety, it has to be remarked that, for increasing , the acceleration will increase as well. The control input constraint includes throttle input or brake input. So, we could define the control formulation as .

#### 4. The Dynamics Model in Each Group

We consider groups of vehicles that travel in a straight line towards a highway tunnel. In this section, we focus on the stability analysis for vehicles in each group . Meanwhile, we suppose the intervehicle distance following a constant time headway policy to be (6). In this way, the actual distance between vehicle and vehicle should be defined as follows: where is vehicle length and , stand for the position for vehicle and vehicle . So, respectively, the continuous vehicle of vehicle distance error is calculated by the actual distance and desired distance as follows:

From the view of vehicle platoon stability, the objective for the vehicle-to-HIS framework is that vehicle interdistance should tend to 0 which holds the desired interdistance between the leading vehicle and the following vehicle. Equation (19) can be defined as follows:

According to slid mode control method principle, when the variable satisfies the following expression:

If the equation is established, means that the th vehicle is said to provide individual vehicle stability. The spacing error derivation process for vehicle platoon is shown in [16], which is used in the paper as follows: where and . By using the above transfer function, string stability of vehicle platoon can be analyzed. It is shown that condition ensures system stability.

#### 5. Simulation and Analysis

In this section, we mainly research the simulation process of the groups. The three time warning distances have been used to ensure multivehicle safety. The first and the second time warning distances mainly lower the deceleration of the leading vehicle in each group. The third time warning distance mainly provides the safety driving distance with groups, for example, deceleration state of safety distance or track mode of safety distance. When the leading vehicle in each group is in the deceleration state during the third warning distance, other vehicles in each group are in a normal safety driving. Therefore, we mainly study the tracking mode for each group during the third warning distance using DMPC. It is assumed that groups , , , , and . The prediction horizon and control horizon for the DPMC are 300 sampling times. At time , the vehicle acceleration is defined:

Meanwhile, we assume that, after the first and the second time warning distance, vehicle of velocity is 10.5 m/s. When the leading vehicle in each group lowers its acceleration in the third time warning distance, the relative distances between the following vehicles in each group and the distance between the last vehicle and the leading vehicle are 3 m, 5 m, 10 m, and 15 m. The following vehicles in each group are tracking the reference trajectories. Hence, the driving state of the leading vehicles can represent each group of information framework. The leading vehicles can only be studied. Four groups of simulation results are shown as shown in Figures 4, 5, 6, and 7.

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

From Figure 4, it can be seen that the intervehicle5-1 distance between the leading vehicle and the vehicle can maintain stability until the time 15 s. When the vehicle decelerated irregularly, the leading vehicle adjusts to its velocity subsequently. Meanwhile, the following vehicles have received the road information ahead from the HIS server system; they still retain the constant safety distance. As to the acceleration curve, the reference acceleration m/s^{2}; by using Definition 5, we can know that only if , the system is controllable. Figure 5 shows the intervehicle distance between vehicles can retain stability until the time 10 s. The change range of the intervelocity is larger than the former group. Generally speaking, the first group and the second group are stability. The third and the fourth group signify that the intervehicle distance between vehicles can retain stability until the time 10 s to 15 s. These groups show that the intervelocity between vehicles changes in real time. The acceleration also adapts to the environment tracking characteristics. The simulated information framework demonstrated that by the first and the second time warning distances to groups of vehicles deceleration system can realize effectively vehicle safety characteristics. Figure 8 shows the intervehicle distances among groups. Through comparative analysis of the intervehicle5-1 distance in the first group, the intervehicle3-1 distance in the second group, the intervehicle4-1 distance in the third group, and the intervehicle5-1 distance in the fourth group, we can conclude that the distributed optimization based DMPC controller and the three time warning distances are carried out to verify the performance of the vehicles passing the highway tunnel environment.

#### 6. Conclusion

The paper presents the distributed model predictive control approach for groups of vehicles of information framework analysis on the highway tunnel environment. The three time warning distances are proposed to ensure vehicles of safety driving. The framework can be adapted for a highway tunnel and some adverse weather such as fog, rain, and snow when drivers cannot distinguish the road ahead because of weather and road reasons. The system is different from traditional vehicle system of active electric functions; its goal is to provide a warning system based on the HIS system for drivers to obtain the road situation in advance, especially multiple groups of vehicles towards some road section prone to accidents or some adverse places. First, the three time warning distances are introduced to guide the mechanism. Second, the distributed model predictive control approach is studied to develop the optimal control sequence. Last, in order to verify the feasibility for the system, four different groups are simulated. The simulation results demonstrate that the proposed three time warning distances could achieve vehicles of safety and stability. Furthermore, it is shown that the local optimization of distributed model predictive control approach could ensure that the whole control performance is effective.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The present work is supported by the National Foundation of Beijing (3133040) and China Postdoctoral Science Foundation (124609).

#### References

- F. Gao, J. Liu, G. Xu, K. Guo, and Y. Cui, “Vehicle state detection in highway intelligent space based on information fusion,” in
*Proceedings of the Twelfth COTA International Conference of Transportation Professionals (CICTP '12)*, pp. 2307–2318, Beijing, China, August 2012. View at: Google Scholar - J. R. Mawhinney, “Fixed fire protection systems in tunnels: issues and directions,”
*Fire Technology*, vol. 49, no. 2, pp. 477–508, 2013. View at: Publisher Site | Google Scholar - T.-T. Chen and Y. T. Hsu, “The research on the operating performance assessment for the automatic incident detection system of the highway tunnel-as an example of ba-guah mountain tunnel in Taiwan,” in
*Proceedings of the GeoHunan International Conference—Tunnel Management, Emerging Technologies, and Innovation*, pp. 40–47, Hunan, China, June 2011. View at: Publisher Site | Google Scholar - B. Yan, H. Chen, and L. Wang, “Visual characteristics of the driver to tunnel group traffic safety,” in
*Proceedings of the 3rd International Conference on Transportation Engineering (ICTE '11)*, pp. 3009–3014, July 2011. View at: Publisher Site | Google Scholar - H. Li and Y. Shi, “Distributed model predictive control of constrained nonlinear systems with communication delays,”
*Systems and Control Letters*, vol. 62, no. 10, pp. 819–826, 2013. View at: Publisher Site | Google Scholar | Zentralblatt MATH - J. Zhang and J. Liu, “Distributed moving horizon state estimation for nonlinear systems with bounded uncertainties,”
*Journal of Process Control*, vol. 23, pp. 1281–1295, 2013. View at: Publisher Site | Google Scholar - W. B. Dunbar and R. M. Murray, “Distributed receding horizon control for multi-vehicle formation stabilization,”
*Automatica*, vol. 42, no. 4, pp. 549–558, 2006. View at: Publisher Site | Google Scholar | MathSciNet - Y. Zheng, S. Li, and N. Li, “Distributed model predictive control over network information exchange for large-scale systems,”
*Control Engineering Practice*, vol. 19, no. 7, pp. 757–769, 2011. View at: Publisher Site | Google Scholar - K. Kon, S. Habasaki, H. Fukushima, and F. Matsuno, “Model predictive based multi-vehicle formation control with collision avoidance and localization uncertainty,” in
*Proceedings of the IEEE/SICE International Symposium on System Integration (SII '12)*, pp. 212–217, Kyushu University, Fukuoka, Japan, December 2012. View at: Publisher Site | Google Scholar - S. G. Kim, M. Tomizuka, and K. H. Cheng, “Smooth motion control of the adaptive cruise control system by a virtual lead vehicle,”
*International Journal of Automotive Technology*, vol. 13, no. 1, pp. 77–85, 2012. View at: Publisher Site | Google Scholar - Y. F. Peng, “Adaptive intelligent backstepping longitudinal control of vehicleplatoons using output recurrent cerebellar model articulation controller,”
*Expert Systems with Applications*, vol. 37, no. 3, pp. 2016–2027, 2010. View at: Publisher Site | Google Scholar - S. Li, K. Li, R. Rajamani, and J. Wang, “Model predictive multi-objective vehicular adaptive cruise control,”
*IEEE Transactions on Control Systems Technology*, vol. 19, no. 3, pp. 556–566, 2011. View at: Publisher Site | Google Scholar - Y. Zhang and S. Li, “Networked model predictive control based on neighbourhood optimization for serially connected large-scale processes,”
*Journal of Process Control*, vol. 17, no. 1, pp. 37–50, 2007. View at: Publisher Site | Google Scholar - G. J. L. Naus, J. Ploeg, M. J. G. van de Molengraft, P. M. H. Heemels, and M. Steinbuch, “Design and implementation of parameterized adaptive cruise control: an explicit model predictive control approach,”
*Control Engineering Practice*, vol. 18, no. 8, pp. 882–892, 2010. View at: Publisher Site | Google Scholar - T. Petrinic and I. Petrovic, “Longitudinal spacing control of vehicles in a platoon for stable and increased traffic flow,” in
*Proceedings of the IEEE International Conference on Control Applications (CCA '12)*, pp. 178–183, Dubrovnik, Croatia, October 2012. View at: Publisher Site | Google Scholar - P. D. Christofides, R. Scattolini, D. Muñoz de la Peña, and J. Liu, “Distributed model predictive control: a tutorial review and future research directions,”
*Computers and Chemical Engineering*, vol. 51, pp. 21–41, 2013. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2014 Tang Xiaofeng 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.