Model Predictive Control for Connected Hybrid Electric Vehicles

This paper presents a new model predictive control system for connected hybrid electric vehicles to improve fuel economy. The new features of this study are as follows. First, the battery charge and discharge profile and the driving velocity profile are simultaneously optimized. One is energy management for HEV for ; the other is for the energy consumption minimizing problem of acc control of two vehicles. Second, a system for connected hybrid electric vehicles has been developed considering varying drag coefficients and the road gradients. Third, the fuel model of a typical hybrid electric vehicle is developed using the maps of the engine efficiency characteristics. Fourth, simulations and analysis (under different parameters, i.e., road conditions, vehicle state of charge, etc.) are conducted to verify the effectiveness of the method to achieve higher fuel efficiency. The model predictive control problem is solved using numerical computation method: continuation and generalized minimum residual method. Computer simulation results reveal improvements in fuel economy using the proposed control method.


Introduction
In recent years, the energy and environmental problems are emphasized.In particular, energy consumption of vehicles accounts for a substantial amount in the transportation sector.There are various approaches to reduce the fuel consumption of vehicles [1][2][3][4][5].High efficient vehicles are being developed to increase fuel economy using lightweight automobiles, efficient power train systems, electric vehicles, and hybrid vehicles [1].On the other hand, the socalled ecodriving can also reduce the fuel consumption [5][6][7][8][9].Ecodriving can be characterized as avoiding aggressive acceleration or braking at any road-traffic situations, cruising at steady speed, decelerating smoothly at stops with little or no braking, and maintaining an optimal distance from the preceding vehicle.An ecological control of a single vehicle on a road with up-down shapes [2] and efficient spacing control of multiple vehicles [10] were presented.
A lot of works have been published on the energy management problem of hybrid electric vehicle (HEV) and plugin hybrid electric vehicle (PHEV) systems.These approaches are typical in a family of optimal control techniques.They can be subdivided into four categories: numerical optimization, analytical optimal control theories, instantaneous optimization, and heuristic control techniques [11].The most representative of numerical optimization is dynamic programming (DP) [11,12].However DP is based on fixed speed patterns which are impossible to get in reality.A kind of analytical optimal control techniques is Pontryagin's minimum principle [13].It gives necessary conditions that the optimal solution must satisfy.It also needs to know the entire driving cycle in advance.The convex optimization method [14] is also a kind of analytical optimal control techniques.The global optimality is guaranteed and the optimal solution can be rapidly and efficiently attained by solvers available.The instantaneous optimization includes the equivalent consumption minimization strategy (ECMS) [2,15].It is based on instantaneous optimization and is easy to implement in real-time.However it cannot guarantee the optimality over the whole driving cycle.Heuristic control techniques like rule-based control strategies [2] are robust, but they are impossible to guarantee the optimality.
To obtain even more fuel economy improvements, connected hybrid electric vehicles can be considered to reduce the air resistance.The air resistance of a vehicle is proportional to the square of the vehicle speed.When a vehicle runs at 100 km/h, its aerodynamic drag reaches more than sixty percent of total motion resistance forces [16].It is obvious that its air resistance causes high fuel consumption.However, the air resistance can be reduced by maintaining a short spacing between two vehicles.Connected vehicles in an automated highway system can lead to increased driver safety, decreased road congestion, and improved fuel economy [17].Connected vehicles can improve fuel economy through reduced wind resistance [18][19][20].
A low pressure area distributes in the rear of the lead vehicle.The size of the area can be different by changing the spacing between the vehicles.If the following vehicle runs at the back of the lead vehicle with a short spacing, its air resistance is decreased owing to improved airflow profile between the vehicles.Furthermore, the air resistance of the lead vehicle is also decreased by the smooth airflow [20].Therefore, fuel consumption of both vehicles can be reduced.However, it is difficult to follow the lead vehicle with a short spacing at high speed by a human driver.Automated cruise control of the vehicle should be introduced to achieve this.Various conventional longitudinal control systems have been proposed such as vehicle following method using information of other vehicles [10] and point following method using a certain decided phase point [21].A control law for internal combustion engine vehicles is proposed which uses relative speed and spacing information from the preceding and following vehicles in order to choose the proper control action for smooth vehicle following and for maintaining a desired intervehicle spacing specified by the driver [22].Connected automatic guided electric vehicles to solve problems of traffic saturation, relying on GPS sensors and intervehicle communication, are addressed in [23].However, these conventional methods consider string stability only.The quantitative effect of road shape and air resistance on fuel consumption for hybrid electric vehicles (HEVs) has not been researched.
For connected hybrid electric vehicles, it is necessary to compute the optimum control inputs of the vehicles by anticipating the future situations including road shape, vehicles' states, and road loads.Therefore, model predictive control (MPC) method can be used.
This paper extends HEV energy management research by adding two novel contributions.First, the battery charge and discharge profile and the driving velocity profile are simultaneously optimized.We make the two connected problems together: one is energy management for HEV for the battery; the other is for the energy consumption minimizing problem of speed control of two vehicles.In reality, the two connected problems are coupled together and affect each other always.The speed of the vehicle affects the charge and discharge profile of the battery.The charge and discharge profile of the battery affects the speed of the vehicle.Second, a new policy between the global optimization method and the instantaneous optimization method is developed.The global optimization method like dynamic programming needs all the information in the future to compute the global optimal control input.The instantaneous optimization method needs no information in the future to compute the control input.The easiest way to deal with the complicated control system is to divide the longitudinal vehicle control system into an upper and lower level controller.The upper level controller determines the desired acceleration of the vehicle on the basis of the position and velocity relative to the other vehicles in the string.The lower level controller determines the input commands to the engine and the braking system, to accomplish the desired acceleration.Also, there is possible to consider road slope, wind, and so forth as a disturbance for the problem.However, in this work we intended to optimize the fuel economy and the speed profile for high fuel efficiency and safety simultaneously.In the HEV operation it is desirable to charge or discharge the battery properly according to the road loads.There is a problem between the fast dynamics components like the engine and the slow dynamics components like the battery.The prediction horizon of the battery state is limited.We developed a new policy to predict the battery state in a longer future for better performance.The desired battery state of charge is designed according to the road slopes for better recuperation of free braking energy.The battery state of charge profile is scheduled systematically to improve fuel economy inside the HEV considering the effect of different parameters, that is, road conditions, battery state of charge, and real-time implementation ability.The quantitative analysis of the vehicle spacing influence and the battery state of charge profile influence for the fuel economy is presented.Performance of the proposed system has been evaluated by computer simulation.The proposed system is found to be more fuel efficient and safer for running over several typical roads with up-down slopes.
The rest of this paper is organized as follows.In Section 2, the nonlinear model of two connected power-split HEVs is derived.Section 3 formulates the nonlinear model predictive control algorithm.Section 4 presents comparative simulation results.Section 5 provides conclusions.

Modeling of Two Connected HEVs
The configuration of the HEV system is shown in Figure 1.FD represents the final drive.The power-split device (PSD) is the key component of the power-split HEV system and has both functionality of speed coupler and continuously variable transmission (CVT).There are five dynamic components: the engine, the battery, two motor/generators (/), and the wheels in this power-split HEV system.The only dynamic state to be considered in the optimal control problem based on known driving cycle is the battery state of charge (SOC) which can simplify the MPC algorithm for implementation.This simplification is possible because this paper introduces four constraints: the road load, the torque and speed relationship of the speed coupler, the power flow relationship among the five components, and the engine optimal operating line (OOL) using CVT.In this work, we assume that the engine works along its OOL using CVT.For simplicity, we assume the two vehicle configurations are the same.It is assumed that the central controller set in the lead vehicle controls the two vehicles.The central controller computes the control inputs of the two vehicles.The control inputs of the two vehicles are fed into the two vehicles, respectively.The states of the two vehicles are measured and sent to the central controller.In this way a closed control loop is formed.Here, we call it central control system of connected vehicles.In a distributed control system of connected vehicles [18][19][20], the individual vehicles are controlled separately by its own controller.It cannot predict other vehicles precisely.In a central control system of connected vehicles, all the information of the vehicles is shared, and the global optimality of all the vehicles can be obtained.The distance between the two vehicles changes (which affect the air drag coefficient) and the slope changes; therefore the stability of this controller is very important.However, the control scheme proposed in this work is brand new; the stability of this controller is completely different from that of the distributed control system of connected vehicles.Hence, we would like to add the stability problem as our future directions because of its complexity.The control signals are transmitted to the vehicles through intervehicle communication.It is assumed that there is no delay of the communication.The proposed methodology will work independently of other kinds of vehicles on the roadway in the network if the vehicle has the functionality of CVT.This paper divided the optimal control problem into two levels.The high-level controller determines the optimal battery power and the low-level controller determines the optimal torque and speed of the engine and the motor/generators.This paper focuses on the high-level controller.
The torque and speed relationship of the speed coupler can be expressed as [24]  eng () = − ( where  and  are the number of sun gear and ring gear teeth, respectively,  /1 ,  /2 ,  req , and  eng are the torques of /1, /2, the road load, and the engine, respectively, and  /1 ,  /2 , and  eng are the angular speeds of /1, /2, and the engine, respectively.The power flow relationships among the five components at the inverter and the power-split device in Figure 1 are given as where  batt ,  /1 ,  /2 ,  eng , and  req are the power of the battery, /1, /2, the engine, and the road load.This paper assumes that the engine always works along its OOL using CVT which can also be considered as a constraint.When the engine power is known, by looking up the table of OOL, the engine speed and torque can be obtained.
This paper evaluates the fuel consumption using Willans line method to reduce the complexity of the engine fuel consumption model.It was found that good approximations are obtained using the Willans line method [25].The fuel consumption can be expressed as where   is a constant.The detailed explanation of this fuel consumption model is included in Appendix A.
The road loads which are the vehicle speed and the required power at the wheels are known when the driving cycle is known.From the configuration of the power-split HEV system, /2 speed is also known as where  /2 is the speed of /2,   is the final drive gear ratio,   is the wheel radius, and V req is the required vehicle speed by the driving cycle.
For simplicity, it is assumed that there are two cars in this central control system.When the driving cycle is unknown, the system dynamics includes the battery and the vehicle dynamics.Both the fuel economy and the driving profile are optimized.The system model is then represented by where , V, and  are the vehicle position, speed, and acceleration or deceleration converted from the traction force or brake force.The parameters ,   , , , , , and () are the air density, the air drag coefficient, the frontal area of the vehicle, the vehicle mass, the gravity acceleration, the rolling resistance coefficient, and the road grade.  ,  ℎ , and   are the vehicle acceleration or deceleration control inputs and the delay constant. OC ,  batt , and  batt are the open-circuit voltage, the internal resistance, and the capacity of the battery.
The suffixes  and ℎ denote the parameters of the preceding vehicle and the host vehicle.The slope information from GPS or the digital map is approximated by the sigmoid functions.This modeling method of road slope is firstly proposed in this work.A description of the modeling method is provided in Appendix B.

Model Predictive Control
The driving control inputs are derived using model predictive control algorithm.The optimal control problem is defined as where  is the prediction horizon and min and max denote the minimum and maximum bounds of the parameters.
The following objectives are considered in this optimal control problem.
The term   : acceleration or deceleration of vehicles is moderated.
The term   : the vehicle speed is kept near to its desired value.
The term   : the fuel consumption is minimized.The term   : the battery SOC is kept near to its desired value.This is one of the cores of the proposed approach.This paper adapts the battery energy to the vehicle future energy requirements by setting the desired battery SOC as a function of road slopes which represent the main part of the future road load.
The term   : the battery energy is made best use of.This is one of the cores of the proposed approach.The battery energy is firstly used to satisfy the required road load.If it is not enough, the engine energy should be used, and the engine can work along its OOL.
The term   : the battery SOC constraint is kept satisfied.The term   : the desired vehicle spacing is kept.This is one of the cores of the proposed approach.The following distance constraint is kept in a predictive controller structure.The following distance is varied above the minimum following distance, which improves the freedom of ecodriving car following control to optimize the driving profile for better fuel economy.
The cost function  is defined as follows: where   ,   ,   ,   ,   ,   , and   are the weights and V  is the desired vehicle speed.The parameters ,   , and   are the vehicle spacing, the desired vehicle spacing, and the length of the preceding vehicle.The parameter SOC  is the desired SOC value.The sigmoid function is chosen to evaluate the vehicle brake fuel consumption.The log barrier function is used as a penalizing term for violations of state constraints.
The structure of the nonlinear model predictive control system is shown in Figure 2. The system inputs contain the control inputs.The system outputs consist of the vehicle states.The predictive controller uses terrain information from a digital map to calculate SOC  () and ().
At each time , the optimal control input is computed by solving the above optimal control problems during the prediction horizon .Only the first element of the optimal control sequence is applied.At the next time step, the prediction horizon moves forward, and the process is repeated.

Comparison Controllers.
There are two simulations in this work.They are the MPC approaches with fixed desired battery SOC and variable desired battery SOC.The aim is to demonstrate how the desired battery SOC affects the fuel economy, the power-split profile, and the drag coefficients.
The desired battery SOC value is set according to the road elevation.The authors think it is reasonable to utilize the road elevation information since this future road load information is known already.The desired battery SOC is assumed to use the function as where  SOC and SOC  are constant parameters set as  SOC = −2 and SOC  = 0.7, respectively.

Simulation Conditions.
In these simulations, the parameters of both HEVs are used from ADVISOR 2002 Toyota Prius data (see Table 1).Seven tuning weights are used in this cost function; this makes the performance very subjective to choice of these weights.The goal is to minimize total fuel used, so the real cost function should be integral of  for minimizing the equivalent fuel cost.  does not need to be tuned ordinarily.The reason is that it is for the state constraint.In reality, it is a physical constraint which cannot be violated.The model predictive control problem is solved using the numerical computation method: the continuation and generalized minimum residual (C/GMRES) method [26].The C/GMRES method uses forward difference approach and discretizes the HEV plant with a sampling interval ℎ  to implement the nonlinear real-time optimal control algorithm.A brief description of the solution of the model predictive control problem using the C/GMRES method is included in Appendix C. The flowchart of the nonlinear model predictive control algorithm implementation is shown in Figure 3. Since the optimization problem is nonlinear and nonconvex, its solution can be local optimal; it is assumed that it can only be solved numerically.A detailed mathematical analysis of robust stability and performance for the proposed method can be found in [26].The fuel economy is calculated using the engine fuel consumption map which is obtained from ADVISOR 2002.
The MPC algorithm is realized by utilizing the C MEX S-function builder in MATLAB/Simulink.First, the optimal battery power is calculated by the high-level controller.Next, this optimal value is fed into the low-level controller where the optimal torque and speed of the engine and /s are determined.Finally, these actual control input signals are applied to the vehicle.The fuel economy is calculated using the quasi-static map of ADVISOR.The backward simulation approach has been employed in this work, as ADVISOR software based on quasi-static maps of power train elements is utilized.However, (3) presents a different formula for fuel consumption estimation which is not consistent with the way The parameters predicted are the road slope based on the GPS data and the traffic conditions.They are not embedded in ADVISOR, as in this software the slope is assumed to be zero.Quasi-static maps of power train elements in ADVISOR are used only for the output evaluation.The traffic conditions and their impact on the fuel consumption are predicted using the connected vehicle model to achieve optimal vehicle spacing for reducing air drag.
A set of data representing the relationship of the aerodynamic drag coefficient and the spacing obtained from a wind tunnel experiment [20] is adopted.By using the seventh-degree polynomial representation, the aerodynamic drag coefficient   can be represented as where parameters  1 - 8 are shown in Table 2.The approximation results of   are shown in Figure 4.

Simulation Results.
The driving profile of the HEV using the MPC algorithm with fixed desired battery SOC and unfixed desired battery SOC (see Figures 5 and 6) shows that the MPC algorithm can use the road slope information well to reduce the fuel consumption.The rows of Figures 5 and 6 from the top are the slope of the road, the speed of the preceding vehicle and the host vehicle, the battery SOC of the preceding vehicle and the host vehicle, and the vehicle distance between the preceding vehicle and the host vehicle.The MPC algorithm simultaneously controls both vehicles by predicting their states, and fast convergence of their spacing is achieved.The vehicle accelerates before the up slope to make use of the kinetic energy.The battery recuperates vehicle braking power during the vehicle down slope driving.The preceding vehicle SOC variation range is smaller than that of the host vehicle.The host vehicle SOC variation range with unfixed desired battery SOC is smaller than that with fixed desired battery SOC.Since the engine needs to charge the battery too often, the overuse of the battery leads to worse fuel economy.The energy profile of the HEV using the MPC algorithm with fixed desired battery SOC and unfixed desired battery SOC (see Figures 7 and 8) shows that the MPC algorithm can make the vehicle drag coefficients converge to minimum.The rows of Figures 7 and 8 from the top are the air drag power of the preceding vehicle and the host vehicle, the drag coefficients of the preceding vehicle and the host vehicle, the fuel consumption rate of the preceding vehicle and the host vehicle, and the total cumulative fuel consumption.The drag coefficients of both vehicles are reduced significantly as a result of the fast convergence of the vehicle distance.Therefore, the drag force is reduced, and fuel savings are achieved.The air drag power of the host vehicle is smaller than that of the preceding vehicle because of the vehicle platooning.The fuel consumption of the HEV with unfixed desired battery SOC is smaller than that with fixed desired battery SOC.The reason is that, during the down slope driving, the MPC algorithm with unfixed desired battery SOC makes better use of the motor regenerative braking than that with fixed desired battery SOC.
All the constraints are satisfied in the simulation.The overall fuel economy results are presented in Table 3.It is shown that the MPC approach with unfixed desired battery SOC can improve fuel economy and keep the final SOC near the initial SOC compared to that with fixed desired battery SOC.There are two reasons.The first reason is that the MPC approach with unfixed desired battery SOC makes the battery use less often than that with fixed desired battery SOC.
The second reason is that the MPC approach with unfixed desired battery SOC reduces the duration and magnitude of the braking during the down slope driving period compared with that with fixed desired battery SOC.To sum up, the reason why MPC approach improves the fuel economy while keeping the SOC near initial status is the variation of the battery SOC.This variation can facilitate charging when the vehicles decelerate.The proposed method can improve the fuel economy of both the preceding vehicle and the host vehicle.The fuel economies with different desired vehicle speed, prediction horizon, and control horizon using the proposed method were shown in Figures 10,11,and 12.It is shown that the best fuel economy occurs when the desired vehicle speed is equal to 50 km/h, the prediction horizon is equal to 2 s, and the control horizon is equal to 0.2 s.
MPC is sensitive to process-model mismatch.Processmodel mismatch is always present: wheels, weather and road conditions, sensor accuracy, and so on.Therefore performance and stability under nominal conditions do not guarantee the robust performance and stability of the real car.In order to analyze the robust stability and performance of the MPC method, errors in slope sensor are assumed to exist while traveling.The proposed MPC method is found to be very robust against slope-sensing error.Figure 9 shows the deviation of the driving and power-split profiles due to errors in slope sensor for the lead vehicle.An error of 20% means that the sensor provides 1.2 times the actual value, whereas −20% error means that the sensor provides 0.8 times the actual value of the slope.Due to a sensing error of 20%, the vehicle speed is a bit lower; the battery SOC is a bit higher; and compared with the vehicle with no slope-sensing error, fuel savings dropped 3%.Similarly, due to a sensing error of −20%, the vehicle speed is a bit higher; the battery SOC is a bit lower; and compared with the vehicle with no slope-sensing error, fuel savings increased 0.8%.Therefore, it can be concluded that, within a reasonable sensing error, the system is robust to maintain its ecological performance without significant deviation.Since the optimization problem is nonlinear and nonconvex, its solution can be local optimal; it is assumed that it can only be solved numerically.A detailed mathematical analysis of robust stability and performance for the proposed method can be found in [26].
The process-model mismatch can also give rise to important unfeasibility issues when computing the MPC control law.A policy to recover from unfeasibility issues is needed to be developed.The numerical computation method cannot manage unfeasibility by itself.Since we are dealing here with vehicles (i.e., people), a policy to recover from unfeasibility issues cannot be developed by reformulating the MPC problem with soft constraints.A switch policy is proposed to cope with the unfeasibility issues.It is summarized as follows.
Step 1. Switch the MPC algorithm to the rule-based algorithm [1] when unfeasibility issues occur.
Step 2. Use the rule-based algorithm for 5 s, and switch back to the MPC algorithm.
Step 3. Go back to Step 1 when unfeasibility issues occur.

Conclusions
A model predictive control system for two connected powersplit HEVs considering the fuel economy, the aerodynamic drag varied by vehicle spacing, and the road shape information has been presented.The performance of the proposed control system was confirmed by the computer simulations.The proposed control method has produced the fast convergence of the vehicle spacing.The excessive acceleration and deceleration have been avoided by predicting the road shapes.The results revealed improvements of the fuel economy considering the effect of different parameters, that is, road conditions, battery state of charge, and real-time implementation ability.Since experiments of vehicles are expensive, we would like to conduct experiments in the future.For the sake of simplicity, we have considered only two cars.In the future, we will add more vehicles to make the model more realistic.

A. Engine Fuel Consumption Model
The proposed engine fuel consumption modeling method is a special method using both Willans line method and the assumption of operating the engine along the engine optimal operating line and is introduced as follows.The HEV parameters are used from the ADVISOR 2002 Toyota Prius HEV data [27].
The Willans line model consists of an affine representation relating the available energy, that is, the energy that is theoretically available for conversion, to the useful energy that is actually present at the output of the energy converter [25].Formally   where the parameter  represents the peak intrinsic energy conversion efficiency of the converter and  loss represents external (parasitic) losses.In fact, this model of energy conversion efficiency is nonlinear, in that the parameters  and  loss are represented as explicit functions of the output flow variable (e.g., engine speed) and are also implicit functions of the effort variable.The modeling method given above is for general engines.However, in this work, the electric CVT can realize idle stop, so  loss becomes zero.When it is assumed that the engine operating points are maintained at the best efficiency, the parameters  can be approximated as a constant.In this case, the fuel consumption rate corresponding to the optimal operating line can be fitted using a linear function.
The engine optimal operating line can be plotted on the engine map as shown in Figure 13.The engine optimal operating points provide the highest efficiency for a given power level.The engine best efficiency related to the engine power according to the engine characteristics is shown in Figure 14.
The fuel consumption rate is estimated as (see Figure 15) where  is the calorific value of the gasoline, which is equal to 34.5 × 10 6 [J/l], and  is the engine efficiency.

B. Road Slope Modeling Method
A brief description of the road slope modeling method is provided as follows.In this research, the sigmoid function is used to model the road slope.The general sigmoid function to model the road slope is expressed as follows: where  1 ,   When it is the up slope case, the sigmoid function is expressed as follows: The figure of the sigmoid function is showed in Figure 16.
When it is the down slope case, the sigmoid function is expressed as follows: The figure of the sigmoid function is showed in Figure 17.When it is the up-down slope case, the sigmoid function is expressed as follows: The figure of the sigmoid function is showed in Figure 18.

C. Solution of the Model Predictive Control Problem
A brief description of the solution of the model predictive control problem is provided as follows.
To implement the model predictive control algorithm, the horizon  is divided into  steps, and the optimal control problem is discretized.The general discretized optimal control problem is formulated as where  is the control input,  is the state, and  is the cost function.(, ) is the state equation.(, ) is the inequality constraint.
The inequality constraint in the optimal control problem is converted to an equality constraint by introducing a dummy input   for computation simplicity as follows: where  max denotes the upper bound of the control input.
To solve this optimal control problem with the calculus of variation method, the Hamiltonian function is defined by  (, , , ) =  (, ) +    (, ) +    (, ) , (C. 3) where  denotes the costate and  denotes the Lagrange multiplier associated with the equality constraint.
The first-order necessary conditions for the optimal control input , the multiplier , and the costate  are obtained using the calculus of variation as where  0 is the initial state.
To solve this optimal control problem, the continuation and GMRES (C/GMRES) method is employed for computation cost reduction.The necessary conditions of optimality for the constrained control input can be expressed as the following equation: The above differential equation can be solved by the GMRES method.The presented approach is also a kind of continuation method.The solution curve () is traced by integrating the above differential equation.Because there is no need to calculate the Jacobians and the linear equation iteratively, C/GMRES method assures the real-time optimal control ability because of small computational cost.The detailed description of the solution for the model predictive control algorithm can be found [26].

Notations and Abbreviations
HEVs: Hybrid electric vehicles MPC: Model predictive control CVT: Continuously variable transmission OOL: Optimal operating line SOC: State of charge PSD: Power-split device /: Motor/generator FD: Final drive : F r o n t a l a r e a   : Aerodynamic drag coefficient : Gravitational acceleration   : Final drive gear ratio : V e h i c l e m a s s  batt : B a t t e r yp o w e r  eng : E n g i n ep o w e r  batt : Battery nominal capacity   : W h e e l r a d i u s  batt : Battery internal resistance  OC : B a t t e r yo p e n -c i r c u i tv o l t a g e : Road inclination : Rolling friction coefficient : D e n s i t y o f a i r : The number of sun gear teeth : Th e n u m b e r o f r i n g g e a r t e e t h  /1 : The torque of motor/generator 1  /2 : The torque of motor/generator 2  req : The torque of the road load  eng : Th et o r q u eo ft h ee n g i n e  /1 : The angular speed of motor/generator 1  /2 : The angular speed of motor/generator 2  eng : The angular speed of the engine  req : The power of the road load   : A constant for the fuel consumption : Th e v e h i c l e p o s i t i o n V: The vehicle speed : The acceleration or deceleration converted from the traction force or brake force   : The acceleration or deceleration control input of the preceding vehicle  ℎ : The acceleration or deceleration control input of the host vehicle   : Th e d e l a y c o n s t a n t  1 ,  2 ,  3 ,  4 ,  5 , and  6 : Slope shape parameters : The prediction horizon min and max : The minimum and maximum bounds of the parameters   ,   ,   ,   ,   ,   , and   : The weights The desired vehicle speed : The vehicle spacing   : The desired vehicle spacing   : The length of the preceding vehicle SOC  : The desired SOC value  SOC and SOC  : Constant parameters  1 - 8 : Constant parameters for the aerodynamic drag coefficients : The calorific value of the gasoline : Th e c o n t r o l i n p u t : Th e s t a t e : Th e c o s t f u n c t i o n (, ): Th es t a t ee q u a t i o n (, ): The inequality constraint   : Th e d u m m y i n p u t  max : The upper bound of the control input : Th e c o s t a t e : The Lagrange multiplier associated with the equality constraint  0 : Th e i n i t i a l s t a t e   : A s t a b l e m a t r i x .

Figure 2 :
Figure 2: Structure of the model predictive control system.

Figure 5 :
Figure 5: Driving profile of the HEV using the MPC algorithm with fixed desired battery SOC.

Figure 6 :
Figure 6: Driving profile of the HEV using the MPC algorithm with unfixed desired battery SOC.

Figure 7 :
Figure 7: Energy profile of the HEV using the MPC algorithm with fixed desired battery SOC.

Figure 8 :
Figure 8: Energy profile of the HEV using the MPC algorithm with unfixed desired battery SOC.
out =  in −  loss ,

Figure 9 :
Figure 9: Driving and power-split profile deviation due to the slope-sensing error.

Figure 10 :
Figure 10: Effects of the desired vehicle speed on fuel economy.

Figure 11 :
Figure 11: Effects of the prediction horizon on fuel economy.

Figure 12 :Figure 13 :
Figure 12: Effects of the control horizon on fuel economy.

Figure 14 :Figure 15 :
Figure 14: The engine efficiency curve to the best engine operating points.

Figure 16 :Figure 17 :
Figure 16: The sigmoid function of the up slope.

Figure 18 :
Figure 18: The sigmoid function of the up-down slope.
the fuel economy and driving performance.  is tuned for minimal fuel consumption by dynamic acceleration of the vehicle.  is tuned for minimal fuel consumption by static maps of the engine which cannot evaluate the dynamic fuel consumption of the engine.  and   are tuned for safety of driving to avoid real-ends collision.  and   are tuned

Table 2 :
Comparison of the energy used by motion resistance forces.
ADVISOR calculates the fuel consumption.The reason is that (3) is for the control input calculation of the model predictive control which needs to be continuous and have derivative, and quasi-static maps of power train elements are used for the output evaluation.

Table 3 :
Fuel economy comparison results.
Figure 4: Approximation of  1 and  2 with respect to the spacing .