Power Prediction-Based Model Predictive Control for Energy Management in Land and Air Vehicle with Turboshaft Engine

In recent years, the green aviation technology draws more attention, and more hybrid power units have been applied to the aerial vehicles. To achieve the high performance and long lifetime of components during varied working conditions, the eﬀective regulation of the energy management is necessary for the vehicles with hybrid power unit (HPU). In this paper, power prediction-based model predictive control (P 2 MPC) for energy management strategy (EMS) is proposed for the vehicle equipped with HPU based on turboshaft engine in order to maintain proper battery’s state of charge (SOC) and decrease turboshaft engine’s exhaust gas temperature (EGT). First, a modeling approach based on data-driven method is adopted to obtain the mathematical model of turboshaft engine considering time delay and inertial of states. An integrated power predictor consisting of the classiﬁcation of input status and the subpredictors are developed based on the deep learning method to improve the accuracy of the prediction model of the model predictive control (MPC). Subsequently, an EMS based on MPC using the proposed power predictor is introduced to regulate the SOC of battery and the EGTof turboshaft engine. The comparison with experimental results shows the high accuracy of mathematical model of turboshaft engine. The simulation results show the eﬀectiveness of the proposed EMS for the vehicle, and the eﬀects of diﬀerent weight coeﬃcients of objective function on the proposed EMS are discussed.


Introduction
With the rapid urbanization and the continuous population growth, more and more vehicles are used to satisfy the requirements of transportation such as freight transportation and public transportation. However, the more vehicles bring the more traffic pressure to the transportation. Conventional vehicles only can run on the land, which cannot make full use of the public space. Land and air vehicle is a promising solution orienting the future transportation, which is capable of running on the land by wheels like conventional vehicle and flying in the air by rotary-wing or fixed-wing. Consequently, the integration of multiple functions of the land and air vehicle can contribute to relieving the traffic pressure and improving traffic efficiency.
However, more driving power is required by the land and air vehicle to meet the demands from different usage situations compared to the conventional vehicles. Generally, it is believed that the greenhouse gas emission, global warming, and air pollution are the main environmental problems as a result of the ongoing use of fossil fuel resources, and more driving power of vehicles means the more use of fossil fuel resources. If the full driving power of the land and air vehicle is still provided by the conventional combustion engine which uses fossil fuel without any improvement, it makes no contribution to energy conservation and pollution reduction. According to the research reports, the aviation sector is responsible for 2% of the total fuel consumption in the transportation sector. If this trend of using fossil fuel resources in aviation continues without any alternative methods, this number is expected to increase up to 11% within the next two decades [1]. To minimize the negative impacts such as air pollution and the greenhouse effect on environment, the green aviation technology which reduces the fossil fuel consumption and emissions in transportation industry is a current main concern. erefore, it is crucial to develop the green aviation technology on the land and air vehicles, and finding alternatives to fossil fuels in aviation is more important now than ever. To achieve low emissions, low fuel consumption, and high efficiency, it is necessary to improve the vehicle power systems greatly.
Hybrid and pure electric driving solutions have been applied in the automotive industry for many years and great success has been gained in the ground vehicle, recently these environmentally friendly solutions are being exploited for unmanned aerial vehicles and general aviation aircraft, and some research results about electric and hybrid power unit power demand estimation, reduction of required power, and battery sizes have been obtained [2][3][4][5][6][7][8]. Moreover, before the energy density of storage system can be increased greatly, the hybrid power unit (HPU) is the most excellent solution and it has more advantages in the endurance mileage because of the combined utilization of engine and energy storage system compared with pure electric power systems. erefore, the applications of HPU on the aerial vehicle are paid more attention. e configuration and architecture of HPU can be mainly divided into two types [5]: series HPU and parallel HPU. e series HPU decouples the engine from the power demand, which provides an opportunity to allow the engine continuously operating in high-efficient or less-emission area. erefore, the series HPU has a promising development on aerial vehicles.
Because of the success of hybrid-electric solution on conventional vehicles, piston engine is usually adopted in HPUs. In fact, turboshaft engine has higher power to weight ratio in comparison to piston engine, especially when it concerns high power requirement, and these advantages make the turboshaft engine have preferable performance in HPU [9]. ere are lots of obvious differences including but not limited to structure and aerospace features between the turboshaft engine and the piston engine, therefore the investigation on HPU with turboshaft engine is worth exploring for better performance of aerial vehicles [10]. To analyze the impacts of turboshaft engine on the HPU, an accurate engine model which can reflect the characteristics is necessary. Some methods have been developed to realize the modeling of turboshaft engine [11][12][13][14][15]. e modeling by mechanism based on aerothermodynamics is a common method, which is adopted widely to obtain the model of gas turbine. e modeling by mechanism can offer the internal parameters of engine such as pressure ratio of compressor and inlet temperature of turbine. However, it heavily relies on the characteristics of components and the knowledge of aerothermodynamics of the researcher. Besides, neutral network modeling method is also applied for modeling of engine. However, a lot of experimental data and time are necessary to train the neutral network model, and the convergent speed and accuracy of model are not guaranteed. Compared with the modeling methods mentioned above, data-driven modeling method uses the experimental data of objects to acquire the steady-state and dynamics characteristics, and fast calculation speed and high accuracy of model by data-driven method are the advantages when the internal parameters of objects are not concerned. erefore, the data-driven method is effective for the modeling of the turboshaft engine. Considering the application of turboshaft engine, the dynamics characteristics of turboshaft engine model are focused in this paper. In reference [16], related gain coefficients among rotate speed, fuel, and EGT are used to model the two-spool turboshaft engine based on the experimental data, and the selections of gain coefficients in simulation rely on the rotate speed. In references [14,17], the dynamics model of one-spool turboshaft engine is built for the startup process, but the dynamics characteristics in loading process at rated rotate speed are not discussed. In fact, the turboshaft engine with one-spool in HPU is used as the main energy source to output power, and the dynamics characteristics of the engine can have an impact on the performance of HPU. erefore, the dynamics characteristics of the engine in loading process at rated rotate speed are worth studying. e allocation of energy among energy sources has profound effect on improving efficiency of the components in the aerial vehicles. e appropriate energy management strategies can greatly develop the potential of HPU. Generally, the EMS of HPU can be divided into two types, namely, rule-based and optimization-based [18][19][20][21][22][23][24][25][26]. e rule-based control depends on the past and current system states to determine the modes of operation. erefore, it is easily implemented as the real-time control strategy. e rule-based energy management strategies (RBEMS) mainly consist of the deterministic rule-based control and fuzzy logic control. e deterministic rule-based control employs the classic set theory and elements to define different modes of operation, and the transitions between the different modes are determined by the designed rules. Different from the deterministic rule-based control, the fuzzy logic control deals with reasoning using an approximate way rather than the precise way. In general, the rule-based control strategies are developed based on the practical experience of researchers, and vehicles using the specific RBEMS can obtain the good performance when the vehicle runs on the specific driving condition. However, it is difficult that the specific RBEMS meets the requirements from different and varying working conditions. On the other hand, the optimizationbased EMS regulates the control variables based on calculation by minimizing a predefined cost function within feasible constraints. e optimization-based approaches can be divided into global optimization and real-time optimization. Dynamic programming (DP) is frequently applied in global optimization control strategy and the vehicle using DP can obtain the best performance of targets while dealing with the hybrid energy management problems. However, the overall working conditions of the future are required to be known in advance for global optimization control, which is extremely difficult in real applications. Among real-time optimization methods, model predictive control (MPC) is a promising method [27][28][29]. According to the MPC theory, by predicting the required parameters in the finite time horizon of the future, the control variables can be gained by minimizing the multiobjective optimization function. Compared to DP, the multiobjective optimization function 2 Complexity of MPC method can be solved online because of the latest prediction parameters during different working conditions, and the good control effect of the MPC method can be obtained if the prediction parameters are accurate enough. For hybrid power system unit, the power demand is the main key factor for accurate prediction model in MPC [30][31][32]. e conventional MPC method is also called as the frozen-time MPC, and the power demand is assumed to be constant equal to that in current time step over the prediction horizon for the prediction model simplification. However, the problem is that this assumption ignores possible changes to the power demands during the prediction horizon, and it may lead to inaccurate prediction model, which can bring poor control performance of MPC method, especially in sharp driving condition. e short sampling time interval in the prediction horizon is the potential method to increase the reliability of prediction model for conventional MPC. However, during actual aerial vehicle operation with HPU, the short sampling time generally means high sampling frequency of sensors and it can lead to high cost on sensors and large electrical energy consumption. erefore, how to obtain the accurate change information of future power demand in prediction horizon is necessary to improve the MPC-based EMS for vehicle with HPU. e power demands closely relate to the velocity of the vehicles. Hence, the velocity prediction method is focused in this paper.
ere are several speed prediction methods [21][22][23]31] including exponential function prediction, Markov chain prediction, and deep learning model prediction. In reference [21], a comparative study shows that an EMS with an ANN-based velocity predictor can achieve better performance in minimizing fuel consumption compared with Markov chain models. In reference [32], combined with the velocity predicted by a radius basis function neural network, an EMS based on the MPC can improve fuel economy better. erefore, the deep learning prediction method is a potential candidate due to its high accuracy. In fact, different driver behaviours and the trend of vehicle states can have an impact on different future speed distributions. However, conventional velocity prediction models usually ignore the characteristics of inputs mentioned above. e internal parameters of conventional velocity prediction model are fixed and cannot change based on the input driving states [22,31,32], which cannot guarantee prediction accuracy. Hence, an integrated velocity prediction model that can change the internal parameters to adapt to various input driving states should be further studied. erefore, in this paper, an energy management strategy based on power prediction MPC for vehicles with turboshaft engine is proposed. e following are the main contributions of this paper: (1) For simplifying model and guaranteeing the accuracy of model, based on experimental data, a modeling approach based on data-driven method is adopted to build mathematical models of turboshaft engine during startup process and loading process considering time delay and inertial of states.
(2) An integrated power predictor developed based on the deep learning method is proposed to forecast the future power demands. e proposed integrated power predictor consists of various subpredictors, and different predictive subpredictors are used for prediction based on the classification of input states from historical velocity and acceleration/brake pedal signal.
(3) To achieve coordinated deviation of SOC of battery and EGT of turboshaft engine, a power predictionbased MPC framework for EMS is designed with the proposed integrated power predictor, and the effect of different weight coefficients of objective function in the proposed EMS is discussed to determine the proper weights combination for desirable control performance.
is paper is organized as follows. In Section 2, the distributed structure of a land and air vehicle is introduced and the models of components such as turboshaft engine are established for EMS research. In Section 3, energy management problem is introduced; the power prediction method based on deep learning for MPC is proposed. In Section 4, the accuracy of the turboshaft engine model based on data-driven method is shown compared with experimental data. e performance comparison between conventional MPC and deep learning prediction-based MPC is performed. Besides, the effects of different weight coefficients in multiobjective function of the proposed EMS on control performance are discussed. e rule-based EMS and the proposed EMS are compared and discussed. Some conclusions are drawn in Section 5.

Configuration of Vehicle.
e distributed structure of land and air vehicle with HPU is designed and produced by Vehicle Research Center in Beijing Institute of Technology. e topological structure is shown in Figure 1. ere are four wheels and distributed driving motors mounted on the chassis to propel the vehicle on land, and eight main propellers with motors provide the lift for air mode. Electric power for moving on the land and flying in the air is from the HPU. e HPU adopts the series hybrid-electric configuration to simplify the vehicle structure, thus there is no mechanical connection between the engine and the drive motor; the configuration of HPU is shown in Figure 2. e turboshaft engine and the generator are connected mechanically as the main energy source of HPU, and the power battery is used as an auxiliary energy source for voltage stabilization and power compensation. e electric power from HPU is transmitted to the driving motors for meeting the power requirement. e application of CAN bus in the vehicle realizes communication among controllers of systems and supports the controllers monitored by the master vehicle control unit, which contributes to the simplification and lightweight design for the vehicle.

Modeling of Vehicle.
Because of the feature of application environments, the vehicle usually runs in land mode during most of the working time. To achieve the high efficiency and long life of components, the research about EMS of the vehicle with HPU under the complex road working condition is important, and the models of components of vehicle such as turboshaft engine and dynamic battery should be established firstly. e rotor wings and the motors which drive rotor wings are not used in land mode, therefore these components are treated as mass loads in the model of vehicle.

Modeling of Turboshaft Engine Based on Data-Driven
Method.
e turboshaft engine is mainly composed of the inlet, compressor, combustion chamber, turbine, gearbox, and other subsystems. e airflow from inlet is compressed by compressor, and the fuel is injected into the airflow to produce the mixture of fuel and air. e mixture of fuel and air is ignited in combustion chamber. e fuel gas from combustion offers the power for turbine to rotate. e turbine can drive the compressor and provide the power to balance loads for stable rotate speed.
Generally, there are two mainly modeling methods for turboshaft engine: the analytical modeling method and the experimental modeling method. e analytical modeling method for turboshaft engine needs accurate components performance data, such as compressor and turbine performance data, and the researchers who use analytical modeling method are required to know lots of knowledge on thermodynamics. However, the turboshaft engine model by the analytical modeling method usually costs lots of time to calculate the parameter value in the simulation test because of a lot of iterations, which is not helpful for the simulation process.
In fact, modeling method by the data-driven method focuses on the input-output parameters based on the experimental data, and it is easier to model the engine for researchers without detailed information about turboshaft engine construction. In this paper, to establish the turboshaft engine model reflecting dynamic response characteristics, a modeling method for dynamics model of turboshaft engine combined with dynamic coefficient based on experimental results data is adopted. e engine starting experiment and load experiment have been conducted on the turboshaft engine to collect the performance data, and the experimental results including speed, fuel flow, torque, and EGT are shown as Figures 3 and 4. In the engine starting experiment, the engine rotation speed increases from 0 r/min to 62000 r/min gradually and the fuel flow also increases from 0 L/h to 47.3 L/h. In the engine load experiment, the output power increases from 0 kW to 115 kW gradually, and similarly the fuel flow   e rotation speed stabilizes at around 62000 r/min during loading experiment because of the engine controller. When the power loads vary, the engine can adjust the output power for balancing the power loads to achieve the purpose of stable speed.
According to turboshaft engine theory, at certain atmospheric temperature and pressure, the output performances of turboshaft engine, such as rotation speed and torque, are approximately determined by the fuel flow when the vane angles of compressor and turbine do not change. As shown in Figures 3 and 4, when the engine rotation speed is stable at around 31000 r/min in the engine starting experiment, the fuel flow also stabilizes at around 20.0 L/h. Similarly, when the output torque is stable at around 7.85 Nm in engine loading experiment, the fuel flow also stabilizes at around 63.5 L/h. e steady-state data relationship is shown in Figure 5.
Based on the data in Figure 5, the steady-state fuel flow is described by where n e is the engine spool rotate speed in engine starting process and T e is the engine output torque in engine loading process, which is different from the work in reference [16] where the n e is used rather than T e in engine loading process. e reason for choosing the T e is that the rotate speed of turboshaft engine in this paper is constant in normal working state. W fs−start and W fs−load are the steadystate fuel flows calculated by interpolation algorithm. When the actual fuel flow changes from the "steady-state point," the actual output performances such as speed or torque also change from the "steady-state point" into dynamic state, which means the change of fuel flow has an impact on dynamic state of engine. e difference between actual fuel flow and steady-state fuel flow usually is known as the actual residual fuel flow. e actual residual fuel flow can be expressed as where W f is the actual fuel flow and W f−x indicates the steady-state fuel flow W fs−start in the engine starting process or W fs−load in the engine loading process, which can be calculated by equation (1); the corresponding change of output torque of turboshaft engine can be expressed as where T e is the actual output torque and T es-x is the steadystate output torque T es-start in the engine starting process or T es-load in the engine loading process, which can be obtained from experimental results in Figures 3 and 4. e dynamics coefficient K t-x can be described respectively as follows: where K t-start is determined by engine rotate speed in engine starting process and K t-load is determined by engine output torque in engine loading process. In addition, the time delay of output torque should be concerned because of the combustion process and power delivery. Assuming time delay parameter is τ, based on the experiment result data, the τ is fitted to polynomial by least squares method, as shown in Figure 6. erefore, the rotor dynamics of output shaft of turboshaft engine are denoted as where T e0 is the initial torque, K t-x is the dynamics coefficient such as K t-start or K t-load , ΔW is the actual residual fuel flow, T load is the actual load torque, J is the rotor inertia of output shaft, and τ is the time delay parameter. Similarly, the steady-state relationship between EGT and fuel flow is shown in Figure 7. e dynamics model of EGT is treated as the first order inertia plant with time delay; the transfer function is described by where K is the steady-state gain coefficient, τ egt is the time delay parameter, and T egt is the time constant, which are fitted to polynomial by least squares method, as shown in Figure 8.

Modeling of Electric Generator and
Motor. e dynamics models of generator and drive motor are described as the first order inertia plants, which are expressed as where Tq g and Tq m are the output torques of the generator and drive motor, respectively, Tq g-cmd and Tq m-cmd are the reference torques of the generator and drive motor, respectively, and Tg and Tm are the time constants.

Modeling of Power Battery.
e internal resistance model is applied to simulate the characteristics of the battery; the equations are as follows: where U ocv is the open-circuit voltage of the battery, U out is the output voltage, I is the battery current, R is the internal resistance of the battery, P motor is the power of the motor and P e-g is the power of engine-generator set, and SOC represents the state of charge of battery.

Modeling of Driving Wheel
Dynamics. e vehicle obtains the driving force from the land by the wheels. e driving wheel dynamics equation is where Tq motor is the output torque of drive motor, J wheel is the equivalent rotor inertia of the driving wheel, i is the drive motor to drive wheel ratio, and Tq f is the driving resistance torque.

Problem Statement.
According to vehicle dynamics, the power provided by HPU should be in balance with the environment resistance power and acceleration resistance power. e power demand P de of the vehicle can be expressed as where P resist_env is the environment resistance power and P acc is the acceleration resistance power. In order to meet the power demand, the total output power of HPU should be equal to power demand P de . e equation is described by where P batt is the power of battery. To achieve the low fuel consumption and good performance, how to effectively distribute the power between P e-g and P batt is focused. Compared to other energy management strategies, MPC is a model-based closed-loop optimal control method whose control action is obtained by solving the open-loop optimal control problem for a prediction horizon. e controlled state variables are monitored and delivered to the MPC objective function for the optimal control sequence calculation based on the predictive model, and only the first element of the control sequence is applied. e other elements are ignored considering the error of predictive model. According to the MPC theory, the accuracy of the predictive model plays a vital role in improving the control performance. For the vehicle with HPU, the power demand from environment is the key factor in the predictive model 8 Complexity [22,33]. In this paper, a power prediction-based MPC (P 2 MPC) method is introduced into the EMS for vehicles with HPU. e operation process of the proposed EMS is displayed in Figure 9. As shown in Figure 9, P 2 MPC method adopts the power predictor with integrated deep neural network models to obtain the future change of power demands. Subsequently, the predictive power demands and the current vehicle states including the SOC and EGT of turboshaft engine are delivered into MPC controller to achieve the optimal control command. Different from the conventional MPC which assumes the power demands to be constant over the prediction horizon, the accurate power demands from the power predictor can improve the accuracy of the prediction model in MPC controller, which can contribute to the high performance of the EMS on the vehicle with HPU.

Deep Learning-Based Power Prediction.
According to the theory of MPC, the accurate future change of state variables of controlled objects in a prediction horizon should be estimated in order to achieve the accurate prediction model. e future change of power demand P de plays a vital role on MPC and is usually determined by vehicle velocity.
ere are several popular vehicle velocity prediction methods: exponential function prediction, Markov chain prediction, deep learning model, and so on. In this paper, the deep neural networks learning is adopted to predict velocity due to its high accuracy. e predictive power demands in upcoming period of time can be obtained by vehicle dynamics calculation.
A deep neural network determines the relationship between its inputs and outputs through training method. Its essence is to learn and form the rules between input and output data. A deep neural network model uses the error value between network output and expected output to train the network structure; the parameters in the neural network can be adjusted to make the mean square error reach the minimum possible value, which makes the network output reach the expected value as closely as possible. Finally, the training process for deep neural network is over and the parameters such as weight and bias in deep neural network can be determined. Different from reference [23], in this paper, a deep neural network predictor (DNNP) uses the historical velocity and the acceleration/brake pedal position as inputs to predict the future velocity in upcoming period of time. Besides, different from reference [21], in order to improve the accuracy of predictor, before the input data are delivered to neural network predictor, a state identifier is applied to identify the state of the input data. e state identifier uses the first order polynomial to fit the historical velocity, and the state of the input data is determined by both of the range of slope of the fitted polynomial and the signal of acceleration/brake. e classification rules for state of input data are shown in Table 1; the α1 and α2 are reference slope values which are determined by predictive effect. e "+" means the acceleration pedal signal is effective, and the "−" means the brake pedal signal is effective. Some different single basic predictors based on deep neural deep network can be used for training according to the state of the input data, and all of the trained single basic predictors are applied to construct an integrated predictor. e single basic predictor scheme is shown in Figure 10, and the integrated predictor scheme for P 2 MPC is shown in Figure 11.
For accurate prediction output, deep neural networks predictor should be trained through a lot of training data. Training sample data consist of the 5 driving cycles of CSHVR, NEDC, 1015_6PRIUS, HWFET, and UNIF01, and these driving cycles are arranged randomly in each training sample to increase training information, as shown in Figure 12.
e training process uses every 6 historical velocities and the current acceleration/brake pedal position as the inputs and the next 5 velocities as the output label, where each time interval is 1 s. Subsequently, the predictive velocity can be used to calculate the power demands in prediction horizon. In this paper, DNNP adopts 2 hidden layers and there are 12 nodes in each layer. In order to compare the accuracy between the single basic predictor and proposed integrated predictor, the root-mean-square error (RMSE) is applied. e calculation equation of RMSE in each second over prediction horizon is expressed as follows: where V p denotes the predictive velocity, V a denotes the actual velocity, i denotes the i-th second in prediction horizon, and P is the total seconds of prediction horizon. e overall RMSE over driving cycle is described as where N denotes the total amount of seconds in driving cycle. e compared results are shown in Table 2. e results indicate that the proposed integrated predictor can obtain better prediction performance. In order to show the accuracy of the predictor in each second over prediction horizon, the calculation equation of RMSE for each second is expressed as follows: where V p x denotes the predictive velocity in the x-th second in prediction horizon, V a x denotes the actual velocity in the x-th second in prediction horizon, and i denotes the i-th second in driving cycle, N denotes the total amount of seconds in driving cycle. e RMSEs between the predictive velocity and actual velocity in 1st, 2nd, 3rd, 4th, and 5th seconds in prediction horizon are compared in Table 3. Based on the trained DNNP, the predictive velocity and actual velocity in the 1st to 5th seconds over prediction horizon are compared in Figure 13. As shown in Figure 13 and Table 3, the prediction velocity in the 1st second of prediction horizon matches actual    Figure 9: Power prediction-based MPC for energy management scheme.   10 Complexity result more closely than that in the 5th second. e compared results show that the prediction accuracy decreases as the sampling moment in prediction horizon increases, and the longer prediction horizon is not necessary. e overall trend of prediction results is consistent with the variation trend of the actual results, and it also indicates the DNNP is suitable for the future power demand prediction of the vehicle in proper prediction horizon.

Prediction Model and Objective Function for P 2 MPC.
e control objectives of P 2 MPC are mainly to regulate SOC around reference value and to decrease the EGT of turboshaft engine, which contribute to improving the power quality and long engine life. erefore, according to equation (8), the dynamic characteristics of the battery SOC are described by e variation of P de over the prediction horizon can be obtained by the prediction of DNNP. According to equation (6), the EGT of turboshaft engine is expressed as the following equation: e η denotes the energy conversion efficiency between turboshaft engine and generator. According to equation (7), the P e-g satisfies the first order inertial process and can be expressed as where P e-g-cmd denotes the reference control command and P e-g is the actual control variable which has been applied to the system. e nonlinear prediction model for P 2 MPC for the vehicle is expressed as follows: where x is the state variable, u is the control variable, v is the measured input variable, and y is the output variable. e small variation of SOC from the reference value can contribute to the good power supply quality and long life of power battery. erefore, the objective optimization function for SOC is expressed as follows: In addition, the low EGT can reduce the damage of components of turboshaft engine and extend the life of HPU. erefore, the objective optimization function for EGT is expressed as follows: In addition, in order to reduce disturbance among the control sequences, the change of the control sequences is also concerned. e objective optimization function for the control sequences is expressed as follows: Different from the works in references [21][22][23], linear weighting method is applied to establish the multiobjective optimization function in this paper, and the multiobjective optimization function for P 2 MPC is described by In equations (20)- (22), where P e-g0 is the initial output power of HPU, SOC r is the SOC reference value, W soc , W EGT , and W u are the weight coefficients, t 0 is the starting time and t f is the ending time in the prediction horizon, and J soc , J EGT , J u are the reference values for normalization. e control sequences can be obtained by minimizing the value of the objective cost function J cost .
To improve the real-time performance of online optimization of P 2 MPC, the Taylor expansion technology is applied for equation (18) to obtain the linear model of nonlinear prediction model and the linear model is expressed as where Here, u 0 , x 0 , and v 0 are the current state variables in current step, u is P e-g . According to equation (17) and the response characteristics of the first order inertial process, assume the control reference command P e-g-cmd does not change between current time t k and next time t k+1 in the prediction horizon, and the average value of P e-g-cmd is regarded as the actual applied control variable. erefore, the correction factor δ is determined by where the ts denotes time interval in the prediction horizon of MPC. To improve the calculation speed, the discrete model of linear prediction model is obtained and expressed as where k represents the current step and k + i represents the i-th step in prediction horizon. e discrete objective cost function is expressed as e interior-point method can be applied to solve equation (28) to obtain the control sequences, and the first element in the control sequences is regarded as the control variable to be employed for the HPU.

Simulation Validation of Turboshaft Engine Dynamics
Model. Turboshaft engine is different from the diesel engine which is widely adopted by hybrid vehicle, and it plays a vital role in the research on EMS. For verifying the accuracy of the turboshaft engine model based on the proposed method, according to the experimental data and modeling process in Section 2.2.1, the simulation model of turboshaft engine has been established and the simulation results are compared with experimental data in Figures 14-16.
As shown in Figure 14, the simulation input is experimental fuel flow data in engine starting progress; the output speed results of model show that the maximum rotate speed relative error between the model output and the experimental data is 4.2% in the engine starting process. As shown in Figures 15 and 16, the simulation inputs are experimental fuel flow data and outputs are engine torque and EGT. e simulation torque results show that the maximum relative error between the output steady-state torque of model and the experimental data is 5%, and the dynamics error of torque is larger than steady-state error, but the dynamics tendency agrees well with the experimental data. Furthermore, the simulation EGT results show that the maximum relative error between the model output EGT and the experimental data is 1.7% in the engine loading process. e simulation results show that the turboshaft engine model by the proposed dynamics coefficient modeling method in this paper has good accuracy, and the proposed data-driven method is very efficient and practical for the modeling of turboshaft engine dynamics.

Simulation of EMS Based on P 2 MPC.
In order to verify the performance achieved by the proposed EMS based on P 2 MPC, the simulation model for the vehicle model with HPU in land mode is established in MATLAB/SIMULINK based on the subsystem models such as turboshaft engine and battery. Besides, the UDDS driving cycle is used and the Complexity 13 driver model is also modelled by PI control method to simulate the power demand. e main parameters of model are shown in Table 4, and the simulation results of the proposed EMS are shown as follows. As shown in Figure 17(a), the performance of the velocity tracking between the velocity from the vehicle simulation model and the reference velocity from the UDDS driving cycle is good. e positive value of the output power P motor of distributed drive motors in Figure 17(b) means the power from HPU is delivered to increase or maintain the velocity of the vehicle, and the negative value of power indicates the power is recovered to reduce the velocity.

Comparison between Conventional MPC and P 2 MPC.
e simulation for comparing the control performance between the conventional MPC and P 2 MPC has been carried out. In this paper, the prediction horizon is set as the future 5 seconds. P dep_x represents the predictive power demand in x-th second in the prediction horizon. Besides, in order to reveal the effect of the predictive power demands on the P 2 MPC, the different configurations of predictive power demand over prediction horizon in MPC approach are set, and the simulations of mean relative deviation of SOC and EGT are compared. e simulation results of conventional MPC method are used as the    Table 5. Besides, in order to better demonstrate the advantage of the proposed predictor in this paper, the simulations for the predictor designed in reference [23] where only historical velocity is used as input without the pedal position are carried out, and the results are shown in Table 6.
As shown in Table 5, the No. 1 configuration shows the results of the conventional MPC where the power demands over the prediction horizon are assumed to be constant P de0, which is the power demand in the current calculation step. No. 2-No. 6 configurations show the simulation results of P 2 MPC. As shown in Table 5, compared to conventional MPC, the increased mean relative errors of SOC of different power demand configurations for P 2 MPC are all negative, which suggest the mean relative error of SOC by P 2 MPC is smaller than that by conventional MPC. Besides, the mean relative errors of EGT by P 2 MPC increase a little, which are tolerable considering the large decrease in the mean relative errors of SOC. e results show that compared to conventional MPC, the higher control performance can be obtained by the proposed P 2 MPC.
As shown in Table 6, due to low accuracy of the predictor without pedal position as input, only No. 2 power demand configuration can achieve better control performance than conventional MPC; the worse control performance of the predictor without pedal position as input is obvious compared with the results in Table 5 where the predictor of P 2 MPC adopts historical velocity and pedal position as inputs, which also indicate the advantage of the proposed power predictor.
Besides, as shown in Table 5, the control performances of No. 3 and No. 4 configurations are worse than that of No. 2 configuration, and the control performances of No. 5 and No. 6 configurations are little better. In addition, as shown in Table 6, the increased mean relative error of SOC of the No. 2 configuration, where the power demands are assumed to be constant P dep_1 over the prediction horizon, is the smallest compared with No. 3-No. 6. e obvious better control performance cannot be furthermore obtained although more predictive power demand results are used over prediction horizon. e reason for the obvious better control performance cannot be obtained is that the prediction error increases with the lapse of time over the prediction horizon. According to Section 3.2, the predicted error of the predictive velocity in the first second over the prediction horizon is the smallest, therefore predictive power demand P dep_1 is closer to the actual power demand in the future, thus the No. 2 configuration of power demands in Tables 5 and 6 can achieve the good performance most possibly. It indicates that the predictive results in the time steps which are closer from the current time moment can improve control performance better. Based on the above analysis, the power predictor of P 2 MPC-based EMS only predicts the power demand in the next second to reduce the calculation time cost, and the No. 2 configuration of power demands in Table 5 is adopted in the calculation of optimal control.

Weight Coefficient of P 2 MPC.
In order to analyze the effect of weight coefficient value on the performance of the proposed EMS, different weight coefficient values W soc , W EGT , and W u are set in multiobjective function according to equation (28). According to Section 3.3, the linear weighting method is used to set the values of W soc , W EGT , and W u . e sum of weight coefficient satisfies the following equation: In this paper, maintaining SOC around the reference value, low EGT and small disturbance between control sequences are concerned. For comparing the control performance, the following 5 indexes are introduced: the maximum relative error (MARE) of SOC from the reference value, the mean relative error (MERE) of SOC from the reference value, the MARE of EGT from the reference value, the MERE of EGT from the reference value, and the fluctuation coefficient (FC) of turboshaft engine power. e fluctuation coefficient of turboshaft engine power h fc is expressed as follows: where N represents the total samples of simulation results, i represents the i-th step, and P e denotes the engine power. In order to compare the effect of the different Table 5: Simulation results of different configuration of predictive power demands for the proposed predictor which adopts velocity and pedal position as inputs. No.

Power demand configuration in prediction horizon
Increased mean relative error compared to conventional MPC method Table 6: Simulation results of different configuration of predictive power demands for the predictor which only adopts velocity as inputs. No.

Power demand configuration in prediction horizon
Increased mean relative error compared to conventional P dep_1 P dep_2 P dep_3 P dep_3 P dep_3 0.2978 −0.0132 5 P dep_1 P dep_2 P dep_3 P dep_4 P dep_4 0.2985 −0.0130 6 P dep_1 P dep_2 P dep_3 P dep_4 P dep_5 0.2958 −0.0129         weight coefficients in the multiobjective function in equation (28), W soc , W u , and W EGT are set according to Table 7. e simulation tests are carried out and the results are shown in Figure 18 and Table 8. As shown in Figure 18 and Table 8, the results suggest that the error of the SOC from reference value 0.6 increases with the weight coefficient W soc decreasing, and large error of the SOC means bad efficiency and short life of the battery. On the other hand, the small weight coefficient W soc means the restriction on SOC is reduced, and more power from power battery is used to satisfy the requirements of driving power from driver, which decreases the power from turboshaft engine and contributes to the low EGT. On the contrary, the larger W EGT has more benefits for low EGT rather than stable SOC.

Complexity
In addition, when W EGT , the W soc is set as the same value in No. 5-6 in Table 7; the control performance is poor because of the large error of SOC from the reference value 0.6. e results indicate that the W soc should be set larger than W EGT to maintain the SOC around the reference value. Although the larger W soc can lead to higher EGT, the increase in relative error of EGT is lower than the decrease in relative error of SOC. erefore, large W soc and small W EGT should be selected for P 2 MPC-based EMS. erefore, in order to determine the final weight coefficient, the W soc , W EGT , and W u are set shown in Table 9, and the simulation results are shown in Figure 19 and Table 10. As shown in Figure 19 and Table 10, when the W soc � 0.9, W EGT � 0.05, and W u � 0.05, the control performance is acceptable. Compared to the results of other configurations, there erefore, based on the above comparison for the effect of different weight coefficients on all control performances, the weight coefficients of multiobject cost function of P 2 MPC are set as W soc � 0.9, W EGT � 0.05, and W u � 0.05.

Comparison between Rule-Based MPC and P 2 MPC.
To better demonstrate the performance achieved by the proposed EMS, the power following rule-based EMS is introduced for comparisons. Besides UDDS, driving cycles of HWFET and LA92 are also used as the simulation references, as shown in Figure 20. In the rule-based EMS, the power of the engine-generator unit will follow the power demand from the signal of accelerator/brake pedal when the SOC belongs to the per-designed range. e engine-generator unit will offer the maximum output power or the minimum output power when the SOC is lower or higher than the limiting value and will follow the power demand from driver when SOC reaches the reference value. e simulation results are shown in Figures 21-23 and Table 11. e maximum EGTs of turboshaft engine by the rulebased strategy are approximate 20 degrees Kelvin higher than that by P 2 MPC over the three driving cycles, and the FCs of engine power by rule-based EMS are larger than that by P 2 MPC-based EMS. e results shown in Figures 21(c), 22(c), and 23(c) also show the turboshaft engine power by rule-based EMS change more sharply than that by P 2 MPC-based EMS, which can lead to the unstable working state of turboshaft engine. e larger deviation of SOC in P 2 MPC-based EMS is tolerable considering maintaining the turboshaft engine in stable working state. e results show P 2 MPC-based EMS can offer the better control performance for vehicle with HPU.

Conclusions
In this paper, a new land and air vehicle powered by HPU with turboshaft engine is introduced, and in order to improve the efficiency and performance of the vehicle, the full model of the vehicle is established and the P 2 MPC-based EMS is proposed. First, a modeling method of turboshaft engine based on data-driven method is proposed, and the great accuracy of turboshaft engine model compared with experimental data is demonstrated by simulation results. To improve the accuracy of the predictive model of proposed P 2 MPC, an integrated power demand predictor with state  where the first order inertial response of control variables is considered. e low EGT, small deviation of SOC from the reference, and small disturbance of control variables are concerned in the multiobjective optimization function of P 2 MPC. e effects of the different weight coefficients on P 2 MPC-based EMS are discussed, and a series of simulation results are demonstrated to determine the appropriate weight coefficients of optimization function. e simulation results between rule-based EMS and P 2 MPCbased EMS demonstrate the advantages of the proposed EMS in this paper.

Data Availability
Readers can obtain the datasets used in this work by contacting the corresponding author via e-mail: yuema@ bit.edu.cn.